Table of contents for issue 2, volume 811, The Astrophysical Journal

Numerical simulations of active galactic nuclei (AGNs) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation, and stellar feedback, focusing on the interplay between cooling, AGN heating, and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the intracluster medium (ICM) and reduces its cooling rate. Within 1–2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and . The star formation rate (SFR) covers a wide range, from 0 to a few hundred , with an average of . The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars.

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by active galactic nuclei (AGNs). Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel/PACS observation of an extremely high-J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3σ upper limit of . We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO spectral line energy distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities of , located at 1.6–5 pc from the AGN, with column densities of at least . At such high column densities, however, dust absorbs most of the CO(40-39) line emission at . Therefore, even if NGC 1068 has a molecular torus that radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

The nearby dwarf starburst galaxy NGC 5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the “radio nebula”). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC 5253 with wavelength coverage from 1500 Å to 1.9 μm in 13 filters. These include Hα, Pβ, and Pα, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the nine optically brightest clusters (M_V < −8.8) and the two young radio nebula clusters. The clusters have ages ∼1–15 Myr and masses ∼1 × 10⁴–2.5 × 10⁵M_⊙. The clusters’ spatial location and ages indicate that star formation has become more concentrated toward the radio nebula over the last ∼15 Myr. The most massive cluster is in the radio nebula; with a mass ∼2.5 × 10⁵M_⊙ and an age ∼1 Myr, it is 2–4 times less massive and younger than previously estimated. It is within a dust cloud with A_V ∼ 50 mag, and shows a clear near-IR excess, likely from hot dust. The second radio nebula cluster is also ∼1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

We present a study of the recent star formation (SF) of 30 Doradus in the Large Magellanic Cloud (LMC) using the panchromatic imaging survey Hubble Tarantula Treasury Project. In this paper we focus on the stars within 20 pc of the center of 30 Doradus, the starburst region NGC 2070. We recovered the SF history by comparing deep optical and near-infrared color–magnitude diagrams (CMDs) with state-of-the-art synthetic CMDs generated with the latest PAdova and TRieste Stellar Evolution Code (PARSEC) models, which include all stellar phases from pre-main-sequence to post-main-sequence. For the first time in this region we are able to measure the SF using intermediate- and low-mass stars simultaneously. Our results suggest that NGC 2070 experienced prolonged activity. In particular, we find that the SF in the region (1) exceeded the average LMC rate ≈ 20 Myr ago, (2) accelerated dramatically ≈ 7 Myr ago, and (3) reached a peak value 1–3 Myr ago. We did not find significant deviations from a Kroupa initial mass function down to . The average internal reddening is found to be between 0.3 and 0.4 mag.

Using all-sky maps obtained with COBE/DIRBE, we reanalyzed the diffuse sky brightness at 1.25 and which consists of zodiacal light, diffuse Galactic light (DGL), integrated starlight (ISL), and isotropic emission including the extragalactic background light. Our new analysis including an improved estimate of the DGL and the ISL with the 2MASS data showed that deviations of the isotropic emission from isotropy were less than 10% in the entire sky at high Galactic latitude (). We derived the DGL to 100 μm brightness ratios of ∼4.79 and ∼1.49 n W m⁻² MJy⁻¹ at 1.25 and 2.2 μm, respectively. The result of our analysis revealed a significantly large isotropic component at 1.25 and with intensities of 60.15 ± 16.14 and respectively. This intensity is larger than the integrated galaxy light, upper limits from γ-ray observation, and potential contribution from exotic sources (i.e., Population III stars, intrahalo light, direct collapse black holes, and dark stars). We therefore conclude that the excess light may originate from the local universe: the Milky Way and/or the solar system.

We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC) using UV absorption spectroscopy. We measure P ii and Fe ii along 84 spatially distributed sightlines toward the MCs using archival Far Ultraviolet Spectroscopic Explorer observations. For 16 of those sightlines, we also measure Si ii, Cr ii, and Zn ii from new Hubble Space Telescope Cosmic Origins Spectrograph observations. We analyze these spectra using a new spectral line analysis technique based on a semi-parametric Voigt profile model. We have combined these measurements with H i and H₂ column densities and reference stellar abundances from the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). Of our 84 P and 16 Zn measurements, 80 and 13, respectively, are depleted by more than 0.1 dex, suggesting that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the SMC. Si, Cr, and Fe are systematically less depleted in the SMC than in the Milky Way (MW) or LMC. The minimum Si depletion in the SMC is consistent with zero. We find GDR ranges of 190–565 in the LMC and 480–2100 in the SMC, which is broadly consistent with GDRs from the literature. These ranges represent actual location to location variation and are evidence of dust destruction and/or growth in the diffuse neutral phase of the interstellar medium. Where they overlap in metallicity, the gas-phase abundances of the MW, LMC, and SMC and damped Lyα systems evolve similarly with metallicity.

We present an analysis of multi-wavelength observations from various data sets and Galactic plane surveys to study the star-formation process in the W42 complex. A bipolar appearance of the W42 complex is evident due to the ionizing feedback from the O5–O6 type star in a medium that is highly inhomogeneous. The Very Large Telescope/NACO adaptive-optics K and L^′ images (resolutions ∼02–01) resolved this ionizing source into multiple point-like sources below ∼5000 AU scale. The position angle ∼15° of the W42 molecular cloud is consistent with the H-band starlight mean polarization angle, which in turn is close to the Galactic magnetic field, suggesting the influence of the Galactic field on the evolution of the W42 molecular cloud. Herschel sub-millimeter data analysis reveals three clumps located along the waist axis of the bipolar nebula, with the peak column densities of ∼(3–5) × 10²² cm⁻² corresponding to visual extinctions of A_V ∼ 32–53.5 mag. The Herschel temperature map traces a temperature gradient in W42, revealing regions of 20 K, 25 K, and 30–36 K. Herschel maps reveal embedded filaments (length ∼1–3 pc) that appear to be radially pointed to the denser clump associated with the O5–O6 star, forming a hub-filament system. A total of 512 candidate young stellar objects (YSOs) are identified in the complex, ∼40% of which are present in clusters distributed mainly within the molecular cloud, including the Herschel filaments. Our data sets suggest that the YSO clusters, including the massive stars, are located at the junction of the filaments, similar to those seen in the Rosette Molecular Cloud.

The C ii lines are important observables for the NASA/SMEX mission Interface Region Imaging Spectrograph. To make three-dimensional (3D) non-LTE radiative transfer computationally feasible, it is crucial to have a model atom with as few levels as possible while retaining the main physical processes. We here develop such a model atom and we study the general formation properties of the C ii lines. We find that a nine-level model atom of C i–C iii with the transitions treated assuming complete frequency redistribution (CRD) suffices to describe the C ii lines. 3D scattering effects are important for the intensity in the core of the line. The lines are formed in the optically thick regime. The core intensity is formed in layers where the temperature is about 10 kK at the base of the transition region. The lines are 1.2–4 times wider than the atomic absorption profile due to the formation in the optically thick regime. The smaller opacity broadening happens for single peak intensity profiles where the chromospheric temperature is low with a steep source function increase into the transition region, the larger broadening happens when there is a temperature increase from the photosphere to the low chromosphere leading to a local source function maximum and a double peak intensity profile with a central reversal. Assuming optically thin formation with the standard coronal approximation leads to several errors: neglecting photoionization severly underestimates the amount of C ii at temperatures below 16 kK, erroneously shifts the formation from 10 kK to 25 kK, and leads to too low intensities.

We use 3D radiation magnetohydrodynamic models to investigate how the thermodynamic quantities in the simulation are encoded in observable quantities, thus exploring the diagnostic potential of the C ii 133.5 nm lines. We find that the line core intensity is correlated with the temperature at the formation height but the correlation is rather weak, especially when the lines are strong. The line core Doppler shift is a good measure of the line-of-sight velocity at the formation height. The line width is both dependent on the width of the absorption profile (thermal and non-thermal width) and an opacity broadening factor of 1.2–4 due to the optically thick line formation with a larger broadening for double peak profiles. The C ii 133.5 nm lines can be formed both higher and lower than the core of the Mg ii k line depending on the amount of plasma in the 14–50 kK temperature range. More plasma in this temperature range gives a higher C ii 133.5 nm formation height relative to the Mg ii k line core. The synthetic line profiles have been compared with Interface Region Imaging Spectrograph observations. The derived parameters from the simulated line profiles cover the parameter range seen in observations but, on average, the synthetic profiles are too narrow. We interpret this discrepancy as a combination of a lack of plasma at chromospheric temperatures in the simulation box and too small non-thermal velocities. The large differences in the distribution of properties between the synthetic profiles and the observed ones show that the C ii 133.5 nm lines are powerful diagnostics of the upper chromosphere and lower transition region.

Detailed observational characterization of transiting exoplanet systems has revealed that the spin-axes of massive () stars often exhibit substantial misalignments with respect to the orbits of the planets they host. Conversely, lower-mass stars tend to only have limited obliquities. A similar trend has recently emerged within the observational data set of young stars’ magnetic field strengths: massive T-Tauri stars tend to have dipole fields that are ∼10 times weaker than their less-massive counterparts. Here we show that the associated dependence of magnetic star–disk torques upon stellar mass naturally explains the observed spin–orbit misalignment trend, provided that misalignments are obtained within the disk-hosting phase. Magnetic torques act to realign the stellar spin-axes of lower-mass stars with the disk plane on a timescale significantly shorter than the typical disk lifetime, whereas the same effect operates on a much longer timescale for massive stars. Cumulatively, our results point to a primordial excitation of extrasolar spin–orbit misalignment, signalling consistency with disk-driven migration as the dominant transport mechanism for short-period planets. Furthermore, we predict that spin–orbit misalignments in systems where close-in planets show signatures of dynamical, post-nebular emplacement will not follow the observed correlation with stellar mass.

The afterglows of gamma-ray bursts (GRBs) are believed to originate from the synchrotron emission of shock-accelerated electrons produced by the interaction between the outflow and the external medium. The accelerated electrons are usually assumed to follow a power-law energy distribution with an index of p. Observationally, although most GRB afterglows have a p larger than 2, there are still a few GRBs suggestive of a hard () electron spectrum. The well-sampled broad-band afterglow data for GRB 091127 show evidence of a hard electron spectrum and strong spectral evolution, with a spectral break moving from high to lower energies. The spectral break evolves very fast and cannot be explained by the cooling break in the standard afterglow model, unless evolving microphysical parameters are assumed. In addition, the multi-band afterglow light curves show an achromatic break at around 33 ks. Based on the model of a hard electron spectrum with an injection break, we interpret the observed spectral break as the synchrotron frequency corresponding to the injection break, and the achromatic break as a jet break caused by the jet-edge effect. It is shown that the spectral evolution and the multi-band afterglow light curves of GRB 091127 can be well reproduced by this model.

We consider the stability of axially unbounded cylindrical flows that contain a toroidal magnetic background field with the same radial profile as their azimuthal velocity. For ideal fluids, Chandrasekhar had shown the stability of this configuration if the Alfvén velocity of the field equals the velocity of the background flow, i.e., if the magnetic Mach number . We demonstrate that magnetized Taylor–Couette flows with such profiles become unstable against non-axisymmetric perturbations if at least one of the diffusivities is finite. We also find that for small magnetic Prandtl numbers the lines of marginal instability scale with the Reynolds number and the Hartmann number. In the limit the lines of marginal instability completely lie below the line for and for they completely lie above this line. For any finite value of , however, the lines of marginal instability cross the line , which separates slow from fast rotation. The minimum values of the field strength and the rotation rate that are needed for the instability (slightly) grow if the rotation law becomes flat. In this case, the electric current of the background field becomes so strong that the current-driven Tayler instability (which also exists without rotation) appears in the bifurcation map at low Hartmann numbers.

The Cygnus OB2 Association is one of the nearest and largest collections of massive stars in the Galaxy. Situated at the heart of the “Cygnus X” complex of star-forming regions and molecular clouds, its distance has proven elusive owing to the ambiguous nature of kinematic distances along this ℓ ≃ 80° sightline and the heavy, patchy extinction. In an effort to refine the three-dimensional geometry of key Cygnus X constituents, we have measured distances to four eclipsing double-lined OB-type spectroscopic binaries that are probable members of Cyg OB2. We find distances of 1.33 ± 0.17, 1.32 ± 0.07, 1.44 ± 0.18, and 1.32 ± 0.13 kpc toward MT91 372, MT91 696, CPR2002 A36, and Schulte 3, respectively. We adopt a weighted average distance of 1.33 ± 0.06 kpc. This agrees well with spectrophotometric estimates for the Association as a whole and with parallax measurements of protostellar masers in the surrounding interstellar clouds, thereby linking the ongoing star formation in these clouds with Cyg OB2. We also identify Schulte 3C (O9.5V), a 4″ visual companion to the 4.75 day binary Schulte 3(A+B), as a previously unrecognized Association member.

The next time a core-collapse supernova (SN) explodes in our galaxy, various detectors will be ready and waiting to detect its emissions of gravitational waves (GWs) and neutrinos. Current numerical simulations have successfully introduced multi-dimensional effects to produce exploding SN models, but thus far the explosion mechanism is not well understood. In this paper, we focus on an investigation of progenitor core rotation via comparison of the start time of GW emission and that of the neutronization burst. The GW and neutrino detectors are assumed to be, respectively, the KAGRA detector and a co-located gadolinium-loaded water Cherenkov detector, either EGADS or GADZOOKS!. Our detection simulation studies show that for a nearby SN (0.2 kpc) we can confirm the lack of core rotation close to 100% of the time, and the presence of core rotation about 90% of the time. Using this approach there is also the potential to confirm rotation for considerably more distant Milky Way SN explosions.

Semiempirical atmospheric modeling attempts to match an observed spectrum by finding the temperature distribution and other physical parameters along the line of sight through the emitting region such that the calculated spectrum agrees with the observed one. In this paper we take the observed spectrum of a sunspot and the quiet Sun in the EUV wavelength range 668–1475 Å from the 2001 SUMER atlas of Curdt et al. to determine models of the two atmospheric regions, extending from the photosphere through the overlying chromosphere into the transition region. We solve the coupled statistical equilibrium and optically thick radiative transfer equations for a set of 32 atoms and ions. The atoms that are part of molecules are treated separately, and are excluded from the atomic abundances and atomic opacities. We compare the Mg ii k line profile observations from the Interface Region Imaging Spectrograph with the profiles calculated from the two models. The calculated profiles for the sunspot are substantially lower than the observed ones, based on the SUMER models. The only way we have found to raise the calculated Mg ii lines to agree with the observations is to introduce illumination of the sunspot from the surrounding active region.

Nonthermal velocities measured in the solar corona imply a strong damping of upward-propagating low-frequency Alfvén waves at heliocentric distances from 1.02 to 1.4 solar radii. We propose a vector Alfvén wave decay as a feasible mechanism for the observed Alfvén wave damping. Contrary to the extensively studied scalar decay, the vector decay does not depend on the wave frequency and can be efficient for low-frequency coronal Alfvén waves. We show that the vector decay is much stronger than the scalar decay and can provide the observed damping of 0.01 Hz coronal Alfvén waves with perpendicular wavelengths of or less. Fully three-dimensional (3D) numerical simulations are needed to capture this decay, whose growth rate is proportional to the vector product of interacting wave vectors.

We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l, b) = (0°, 0°) and (329°, 0°) at both Q (43 GHz) and W bands (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hr of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 μK deg, 2.3–6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 273 and 128 FWHM at Q and W bands, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ≈0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM = −4000 ± 200 rad m⁻²), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information.

We report small-scale clustering measurements from the PRIsm MUlti-object Survey (PRIMUS) spectroscopic redshift survey as a function of color and luminosity. We measure the real-space cross-correlations between 62,106 primary galaxies with PRIMUS redshifts and a tracer population of ∼545,000 photometric galaxies over redshifts from z = 0.2 to z = 1. We separately fit a power-law model in redshift and luminosity to each of three independent color-selected samples of galaxies. We report clustering amplitudes at fiducial values of z = 0.5 and The clustering of the red galaxies is times as strong as that of the blue galaxies and as strong as that of the green galaxies. We also find that the luminosity dependence of the clustering is strongly dependent on physical scale, with greater luminosity dependence being found between and , compared to the to range. Moreover, over a range of two orders of magnitude in luminosity, a single power-law fit to the luminosity dependence is not sufficient to explain the increase in clustering at both the bright and faint ends at the smaller scales. We argue that luminosity-dependent clustering at small scales is a necessary component of galaxy-halo occupation models for blue, star-forming galaxies as well as for red, quenched galaxies.

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at using high signal-to-noise ratio (S/N) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass , and velocity dispersion of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass for each object. The quasars are preferentially hosted by massive galaxies with characterized by stellar ages around 1 billion yr, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past 1 billion yr, which was subsequently quenched or suppressed. The derived and relations agree with our past measurements and are consistent with no evolution from the local universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-S/N fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

The finding of residual gas in the large central cavity of the HD 142527 disk motivates questions regarding the origin of its non-Keplerian kinematics and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6–5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuths. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low-mass companion, whose orbit may be contained in the plane of the inner disk. The companion’s high inclination with respect to the massive outer disk could drive Kozai oscillations over long timescales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD 142527.

We constrain the minimum variability timescales for 938 gamma-ray bursts (GRBs) observed by the Fermi/Gamma-ray Burst Monitor instrument prior to 2012 July 11. The tightest constraints on progenitor radii derived from these timescales are obtained from light curves in the hardest energy channel. In the softer bands—or from measurements of the same GRBs in the hard X-rays from Swift—we show that variability timescales tend to be a factor of two to three longer. Applying a survival analysis to account for detections and upper limits, we find median minimum timescale in the rest frame for long-duration and short-duration GRBs of 45 and 10 ms, respectively. Less than 10% of GRBs show evidence for variability on timescales below 2 ms. These shortest timescales require Lorentz factors and imply typical emission radii cm for long-duration GRBs and cm for short-duration GRBs. We discuss implications for the GRB fireball model and investigate whether or not GRB minimum timescales evolve with cosmic time.

The vortical motions of the baryonic gas residing in large-scale structures are investigated by cosmological hydrodynamic simulations. Proceeding in the formation of the cosmic web, the vortical motions of baryonic matter are pumped up by baroclinity in two stages, i.e., the formation of sheets and filaments. The mean curl velocities are about <1, 1–10, 10–150, and 5–50 km s⁻¹ in voids, sheets, filaments, and knots at z = 0, respectively. The scaling of the vortical velocity of gas can be well described by the She–Leveque hierarchical turbulence model in the range of l < 0.65(1.50) h⁻¹ Mpc in a simulation with a box of size 25(100) h⁻¹ Mpc. The fractal Hausdorff dimension of vortical motions, d, revealed by velocity structure functions, is ∼2.1–2.3(∼1.8–2.1). It is slightly larger than the fractal dimension of mass distribution in filaments, D^f ∼ 1.9–2.2, and smaller than the fractal dimension of sheets, D^s ∼ 2.4–2.7. The vortical kinetic energy of baryonic gas is mainly transported by filaments. Both scalings of mass distribution and vortical velocity increments show distinctive transitions at the turning scale of ∼0.65(1.50) h⁻¹ Mpc, which may be closely related to the characteristic radius of density filaments.

We conduct a pilot investigation to determine the optimal combination of color and variability information to identify quasars in current and future multi-epoch optical surveys. We use a Bayesian quasar selection algorithm to identify 35,820 type 1 quasar candidates in a 239 field of the Sloan Digital Sky Survey (SDSS) Stripe 82, using a combination of optical photometry and variability. Color analysis is performed on 5-band single- and multi-epoch SDSS optical photometry to a depth of From these data, variability parameters are calculated by fitting the structure function of each object in each band with a power-law model using 10 to observations over timescales from ∼1 day to ∼8 years. Selection was based on a training sample of 13,221 spectroscopically confirmed type-1 quasars, largely from the SDSS. Using variability alone, colors alone, and combining variability and colors we achieve 91%, 93%, and 97% quasar completeness and 98%, 98%, and 97% efficiency, respectively, with particular improvement in the selection of quasars at where quasars and stars have similar optical colors. The 22,867 quasar candidates that are not spectroscopically confirmed reach a depth of 21,876 (95.7%) are dimmer than coadded i-band magnitude of 19.9, the cutoff for spectroscopic follow-up for SDSS on Stripe 82. Brighter than 19.9, we find 5.7% more quasar candidates without confirming spectra in sky regions otherwise considered complete. The resulting quasar sample has sufficient purity (and statistically correctable incompleteness) to produce a luminosity function comparable to those determined by spectroscopic investigations. We discuss improvements that can be made to the process in preparation for performing similar photometric selection and science on data from post-SDSS sky surveys.

The nearby, middle-aged PSR B1055−52 has many properties in common with the Geminga pulsar. Motivated by the Geminga's enigmatic and prominent pulsar wind nebula (PWN), we searched for extended emission around PSR B1055−52 with Chandra ACIS. For an energy range 0.3–1 keV, we found a flux enhancement in a annulus around the pulsar. There is a slight asymmetry in the emission close, , to the pulsar. The excess emission has a luminosity of about 10²⁹ erg s⁻¹ in an energy range 0.3–8 keV for a distance of 350 pc. Overall, the faint extended emission around is consistent with a PWN of an aligned rotator moving away from us along the line of sight with supersonic velocity, but a contribution from a dust scattering halo cannot be excluded. Comparing the properties of other nearby, middle-aged pulsars, we suggest that the geometry—the orientations of rotation axis, magnetic field axis, and the sight-line—is the deciding factor for a pulsar to show a prominent PWN. We also report on an flux decrease of PSR B1055−52 between the 2000 XMM-Newton and our 2012 Chandra observation. We tentatively attribute this flux decrease to a cross-calibration problem, but further investigations of the pulsar are required to exclude actual intrinsic flux changes.

We demonstrate that ∼10 s after the core-collapse of a massive star, a thermonuclear explosion of the outer shells is possible for some (tuned) initial density and composition profiles, assuming that the neutrinos failed to explode the star. The explosion may lead to a successful supernova, as first suggested by Burbidge et al. We perform a series of one-dimensional (1D) calculations of collapsing massive stars with simplified initial density profiles (similar to the results of stellar evolution calculations) and various compositions (not similar to 1D stellar evolution calculations). We assume that the neutrinos escaped with a negligible effect on the outer layers, which inevitably collapse. As the shells collapse, they compress and heat up adiabatically, enhancing the rate of thermonuclear burning. In some cases, where significant shells of mixed helium and oxygen are present with pre-collapsed burning times of ≲100 s (≈10 times the free-fall time), a thermonuclear detonation wave is ignited, which unbinds the outer layers of the star, leading to a supernova. The energy released is small, ≲10⁵⁰ erg, and negligible amounts of synthesized material (including ⁵⁶Ni) are ejected, implying that these 1D simulations are unlikely to represent typical core-collapse supernovae. However, they do serve as a proof of concept that the core-collapse-induced thermonuclear explosions are possible, and more realistic two-dimensional and three-dimensional simulations are within current computational capabilities.

The collapse of the primordial gas in the density regime ∼10⁸−10¹⁰ cm⁻³ is controlled by the three-body H₂ formation process, in which the gas can cool faster than free-fall time—a condition proposed as the chemothermal instability. We investigate how the heating and cooling rates are affected during the rapid transformation of atomic to molecular hydrogen. With a detailed study of the heating and cooling balance in a 3D simulation of Pop III collapse, we follow the chemical and thermal evolution of the primordial gas in two dark matter minihalos. The inclusion of sink particles in modified Gadget-2 smoothed particle hydrodynamics code allows us to investigate the long-term evolution of the disk that fragments into several clumps. We find that the sum of all the cooling rates is less than the total heating rate after including the contribution from the compressional heating (pdV). The increasing cooling rate during the rapid increase of the molecular fraction is offset by the unavoidable heating due to gas contraction. We conclude that fragmentation occurs because H₂ cooling, the heating due to H₂ formation and compressional heating together set a density and temperature structure in the disk that favors fragmentation, not the chemothermal instability.

The complete understanding of the stellar abundances of lithium, beryllium, and boron represents one of the most interesting open problems in astrophysics. These elements are largely used to probe stellar structure and mixing phenomena in different astrophysical scenarios, such as pre-main-sequence or main-sequence stars. Their different fragility against (p,α) burning reactions allows one to investigate different depths of the stellar interior. Such fusion mechanisms are triggered at temperatures between T ≈ (2–5) × K, thus defining a corresponding Gamow energy between ≈ 3–10 keV, where S(E)-factor measurements need to be performed to get reliable reaction rate evaluations. The Trojan Horse Method is a well defined procedure to measure cross sections at Gamow energies overcoming the uncertainties due to low-energy S(E)-factor extrapolation as well as electron screening effects. Taking advantage of the measure of the ⁹Be(p,α)⁶Li and ¹⁰B(p,α)⁷Be cross sections, the corresponding reaction rates have been calculated and compared with the evaluations by the NACRE collaboration, widely used in the literature. The impact on surface abundances of the updated ⁹Be and ¹⁰B (p,α) burning rates is discussed for pre-MS stars.

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500, and 850 μm. The disk is resolved at 70, 160, and 450 μm. In addition to the planetesimal belt, we detect thermal emission from AU Mic’s halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is and its milimeter-grain dust mass is (±20%). We create a simple spatial model that reconciles the disk spectral energy distribution as a blackbody of 53 ± 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best-fit model is consistent with the “birth ring” model explored in earlier works, i.e., an edge-on dust belt extending from 8.8 to 40 AU, but with an additional halo component with an surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass-loss rates of 10–100 times solar, compact (zero porosity) grains can only be removed if they are very small; consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 μm can be removed via corpuscular forces (i.e., the stellar wind).

3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our graphics processing unit hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet’s Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the “transient” horseshoe flow. The flow continues as it rolls up into a pair of up–down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of three slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.

Discoveries from the prime Kepler mission demonstrated that small planets (<3 ) are common outcomes of planet formation. While Kepler detected many such planets, all but a handful orbit faint, distant stars and are not amenable to precise follow up measurements. Here, we report the discovery of two small planets transiting K2-21, a bright (K = 9.4) M0 dwarf located pc from Earth. We detected the transiting planets in photometry collected during Campaign 3 of NASA’s K2 mission. Analysis of transit light curves reveals that the planets have small radii compared to their host star, / = % and %, respectively. We obtained follow up NIR spectroscopy of K2-21 to constrain host star properties, which imply planet sizes of 1.59 ± 0.43 and 1.92 ± 0.53 , respectively, straddling the boundary between high-density, rocky planets and low-density planets with thick gaseous envelopes. The planets have orbital periods of 9.32414 days and 15.50120 days, respectively, and a period ratio = 1.6624, very near to the 5:3 mean motion resonance, which may be a record of the system’s formation history. Transit timing variations due to gravitational interactions between the planets may be detectable using ground-based telescopes. Finally, this system offers a convenient laboratory for studying the bulk composition and atmospheric properties of small planets with low equilibrium temperatures.

We report the discovery of a low-mass companion to HR 3549, an A0V star surrounded by a debris disk with a warm excess detected by WISE at 22 μm (10σ significance). We imaged HR 3549 B in the L band with NAOS-CONICA, the adaptive optics infrared camera of the Very Large Telescope, in January 2013 and confirmed its common proper motion in 2015 January. The companion is at a projected separation of ≃80 AU and position angle of ≃157°, so it is orbiting well beyond the warm disk inner edge of r > 10 AU. Our age estimate for this system corresponds to a companion mass in the range 15–80 M_J, spanning the brown dwarf regime, and so HR 3549 B is another recent addition to the growing list of brown dwarf desert objects with extreme mass ratios. The simultaneous presence of a warm disk and a brown dwarf around HR 3549 provides interesting empirical constraints on models of the formation of substellar companions.

The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement of the spatial distribution of density, however, is difficult, and indirect seismology inversion methods have been suggested as an alternative. We applied Bayesian inference, model comparison, and model-averaging techniques to the inference of the cross-field density structuring in solar magnetic waveguides using information on periods and damping times for resonantly damped magnetohydrodynamic transverse kink oscillations. Three commonly employed alternative profiles were used to model the variation of the mass density across the waveguide boundary. Parameter inference enabled us to obtain information on physical quantities such as the Alfvén travel time, the density contrast, and the transverse inhomogeneity length scale. The inference results from alternative density models were compared and their differences quantified. Then, the relative plausibility of the considered models was assessed by performing model comparison. Our results indicate that the evidence in favor of any of the three models is minimal, unless the oscillations are strongly damped. In such a circumstance, the application of model-averaging techniques enables the computation of an evidence-weighted inference that takes into account the plausibility of each model in the calculation of a combined inversion for the unknown physical parameters.

Spectra of the cellular photospheric flows are determined from full-disk Doppler velocity observations acquired by the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory spacecraft. Three different analysis methods are used to separately determine spectral coefficients representing the poloidal flows, the toroidal flows, and the radial flows. The amplitudes of these spectral coefficients are constrained by simulated data analyzed with the same procedures as the HMI data. We find that the total velocity spectrum rises smoothly to a peak at a wavenumber of about 120 (wavelength of about 35 Mm), which is typical of supergranules. The spectrum levels off out to wavenumbers of about 400, and then rises again to a peak at a wavenumber of about 3500 (wavelength of about 1200 km), which is typical of granules. The velocity spectrum is dominated by the poloidal flow component (horizontal flows with divergence but no curl) at wavenumbers above 30. The toroidal flow component (horizontal flows with curl but no divergence) dominates at wavenumbers less than 30. The radial flow velocity is only about 3% of the total flow velocity at the lowest wavenumbers, but increases in strength to become about 50% at wavenumbers near 4000. The spectrum compares well with the spectrum of giant cell flows at the lowest wavenumbers and with the spectrum of granulation from a 3D radiative-hydrodynamic simulation at the highest wavenumbers.

Advanced three-dimensional (3D) radiative MHD simulations now reproduce many properties of the outer solar atmosphere. When including a domain from the convection zone into the corona, a hot chromosphere and corona are self-consistently maintained. Here we study two realistic models, with different simulated areas, magnetic field strength and topology, and numerical resolution. These are compared in order to characterize the heating in the 3D-MHD simulations which self-consistently maintains the structure of the atmosphere. We analyze the heating at both large and small scales and find that heating is episodic and highly structured in space, but occurs along loop-shaped structures, and moves along with the magnetic field. On large scales we find that the heating per particle is maximal near the transition region and that widely distributed opposite-polarity field in the photosphere leads to a greater heating scale height in the corona. On smaller scales, heating is concentrated in current sheets, the thicknesses of which are set by the numerical resolution. Some current sheets fragment in time, this process occurring more readily in the higher-resolution model leading to spatially highly intermittent heating. The large-scale heating structures are found to fade in less than about five minutes, while the smaller, local, heating shows timescales of the order of two minutes in one model and one minutes in the other, higher-resolution, model.

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/Solar Optical Telescope Spectro-Polarimeter scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

Using high-resolution 3D and 2D (axisymmetric) hydrodynamic simulations in spherical geometry, we study the evolution of cool cluster cores heated by feedback-driven bipolar active galactic nuclei (AGNs) jets. Condensation of cold gas, and the consequent enhanced accretion, is required for AGN feedback to balance radiative cooling with reasonable efficiencies, and to match the observed cool core properties. A feedback efficiency (mechanical luminosity where is the mass accretion rate at 1 kpc) as small as 6 × 10⁻⁵ is sufficient to reduce the cooling/accretion rate by ∼10 compared to a pure cooling flow in clusters (with ). This value is much smaller compared to the ones considered earlier, and is consistent with the jet efficiency and the fact that only a small fraction of gas at 1 kpc is accreted onto the supermassive black hole (SMBH). The feedback efficiency in earlier works was so high that the cluster core reached equilibrium in a hot state without much precipitation, unlike what is observed in cool-core clusters. We find hysteresis cycles in all our simulations with cold mode feedback: condensation of cold gas when the ratio of the cooling-time to the free-fall time () is ≲10 leads to a sudden enhancement in the accretion rate; a large accretion rate causes strong jets and overheating of the hot intracluster medium such that further condensation of cold gas is suppressed and the accretion rate falls, leading to slow cooling of the core and condensation of cold gas, restarting the cycle. Therefore, there is a spread in core properties, such as the jet power, accretion rate, for the same value of core entropy or . A smaller number of cycles is observed for higher efficiencies and for lower mass halos because the core is overheated to a longer cooling time. The 3D simulations show the formation of a few-kpc scale, rotationally supported, massive () cold gas torus. Since the torus gas is not accreted onto the SMBH, it is largely decoupled from the feedback cycle. The radially dominant cold gas (T < 5 × 10⁴ K; ) consists of fast cold gas uplifted by AGN jets and freely infalling cold gas condensing out of the core. The radially dominant cold gas extends out to 25 kpc for the fiducial run (halo mass and feedback efficiency 6 × 10⁻⁵), with the average mass inflow rate dominating the outflow rate by a factor of ≈2. We compare our simulation results with recent observations.

Soft lags from the emission of the lower kilohertz quasi-periodic oscillations (kHz QPOs) of neutron star low-mass X-ray binaries have been reported from 4U1608-522 and 4U1636-536. Those lags hold prospects for constraining the origin of the QPO emission. In this paper, we investigate the spectral-timing properties of both the lower and upper kHz QPOs from the neutron star binary 4U1728-34, using the entire Rossi X-Ray Timing Explorer archive on this source. We show that the lag-energy spectra of the two QPOs are systematically different: while the lower kHz QPO shows soft lags, the upper kHz QPO shows either a flat lag-energy spectrum or hard variations lagging softer variations. This suggests two different QPO-generation mechanisms. We also performed the first spectral deconvolution of the covariance spectra of both kHz QPOs. The QPO spectra are consistent with Comptonized blackbody emission, similar to the one found in the time-averaged spectrum, but with a higher seed-photon temperature, suggesting that a more compact inner region of the Comptonization layer (boundary/spreading layer, corona) is responsible for the QPO emission. Considering our results together with other recent findings, this leads us to the hypothesis that the lower kHz QPO signal is generated by coherent oscillations of the compact boundary layer region itself. The upper kHz QPO signal may then be linked to less-coherent accretion-rate variations produced in the inner accretion disk, and is then detected when they reach the boundary layer.

The B-type hypergiant Cygnus OB2 no. 12 is a popular target for studies of interstellar phenomena at visible-infrared wavelengths because of its exceptional brightness for a star dimmed by some 10 mag of visual extinction. A lack of detectable ice absorption has led investigators to regard the line of sight as a standard for studies of the “diffuse” interstellar medium (ISM), an assumption challenged both by observations of molecular gas toward the star and by uncertainties concerning the degree to which such a luminous object may affect its local environment. This paper presents a reassessment of the nature of the material responsible for extinction toward Cyg OB2 no. 12. The excess relative to other cluster members appears to occur in translucent clumps within an extensive network of clouds in the region. Attenuation of the ambient radiation field is sufficient in the cores of the clumps to support the presence of gas-phase molecules, but not to sustain detectable ice formation. In general, the optical properties of dust in the clumps are closely similar to those observed in typical diffuse interstellar material, with the notable exception of an unusually low value for the wavelength of maximum polarization. The implied enhancement of polarization by small grains is attributed to increased alignment efficiency in an enhanced magnetic field. This caveat apart, the results of the current paper provide reassurance that Cyg OB2 no. 12 is, indeed, an appropriate choice for studies that target diffuse and translucent phases of the ISM.

We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344-4243). Deep optical imaging reveals previously undetected filaments of star formation, extending to radii of ∼50–100 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy reveals a massive (2 × 10⁹M_⊙), young (∼4.5 Myr) population of stars, consistent with a time-averaged star formation rate of 610 ± 50 M_⊙ yr⁻¹. We report a strong detection of O viλλ1032,1038, which appears to originate primarily in shock-heated gas, but may contain a substantial contribution (>1000 M_⊙ yr⁻¹) from the cooling intracluster medium (ICM). We confirm the presence of deep X-ray cavities in the inner ∼10 kpc, which are among the most extreme examples of radio-mode feedback detected to date, implying jet powers of 2–7 × 10⁴⁵ erg s⁻¹. We provide evidence that the active galactic nucleus inflating these cavities may have only recently transitioned from “quasar-mode” to “radio-mode,” and may currently be insufficient to completely offset cooling. A model-subtracted residual X-ray image reveals evidence for prior episodes of strong radio-mode feedback at radii of ∼100 kpc, with extended “ghost” cavities indicating a prior epoch of feedback roughly 100 Myr ago. This residual image also exhibits significant asymmetry in the inner ∼200 kpc (0.15R₅₀₀), reminiscent of infalling cool clouds, either due to minor mergers or fragmentation of the cooling ICM. Taken together, these data reveal a rapidly evolving cool core which is rich with structure (both spatially and in temperature), is subject to a variety of highly energetic processes, and yet is cooling rapidly and forming stars along thin, narrow filaments.

The study of kinetic effects in heliospheric plasmas requires representation of dynamics at sub-proton scales, but in most cases the system is driven by magnetohydrodynamic (MHD) activity at larger scales. The latter requirement challenges available computational resources, which raises the question of how large such a system must be to exhibit MHD traits at large scales while kinetic behavior is accurately represented at small scales. Here we study this implied transition from kinetic to MHD-like behavior using particle-in-cell (PIC) simulations, initialized using an Orszag–Tang Vortex. The PIC code treats protons, as well as electrons, kinetically, and we address the question of interest by examining several different indicators of MHD-like behavior.

We have analyzed a sample of 27,258 fundamental-mode RR Lyrae variable stars (type RRab) detected recently toward the Galactic bulge by the Optical Gravitational Lensing Experiment (OGLE) survey. The data support our earlier claim that these metal-poor stars trace closely the barred structure formed of intermediate-age red clump giants. The distance to the Galactic center (GC) inferred from the bulge RR Lyrae stars is kpc. We show that their spatial distribution has the shape of a triaxial ellipsoid with a major axis located in the Galactic plane and inclined at an angle of to the Sun–GC line of sight. The obtained scale-length ratio of the major axis to the minor axis in the Galactic plane and to the axis vertical to the plane is 1:0.49(2):0.39(2). We do not see evidence for the bulge RR Lyrae stars forming an X-shaped structure. Based on the light curve parameters, we derive metallicities of the RRab variables and show that there is a very mild but statistically significant radial metallicity gradient. About 60% of the bulge RRab stars form two very close sequences on the period–amplitude (or Bailey) diagram, which we interpret as two major old bulge populations: A and B. Their metallicities likely differ. Population A is about four times less abundant than the slightly more metal-poor population B. Most of the remaining stars seem to represent other, even more metal-poor populations of the bulge. The presence of multiple old populations indicates that the Milky Way bulge was initially formed through mergers.

We present the deepest optical photometry for any dwarf elliptical (dE) galaxy based on Hubble Space Telescope Advanced Camera for Surveys (ACS) observations of the Local Group dE galaxies NGC 147 and NGC 185. Our F606W and F814W color–magnitude diagrams are the first to reach below the oldest main sequence turnoff in a dE galaxy, allowing us to determine full star formation histories in these systems. The ACS fields are located roughly ∼1.5 effective radii from the galaxy center to avoid photometric crowding. While both ACS fields show unambiguous evidence for old and intermediate age stars, the mean age of NGC 147 is ∼4–5 Gyr younger as compared to NGC 185. In NGC 147, only 40% of stars were in place 12.5 Gyr ago (z ∼ 5), with the bulk of the remaining stellar population forming between 5 to 7 Gyr. In contrast, 70% of stars were formed in NGC 185 prior to 12.5 Gyr ago with the majority of the remaining population forming between 8 to 10 Gyr ago. Star formation has ceased in both ACS fields for at least 3 Gyr. Previous observations in the central regions of NGC 185 show evidence for star formation as recent as 100 Myr ago, and a strong metallicity gradient with radius. This implies a lack of radial mixing between the center of NGC 185 and our ACS field. The lack of radial gradients in NGC 147 suggests that our inferred SFHs are more representative of its global history. We interpret the inferred differences in star formation histories to imply an earlier infall time into the M31 environment for NGC 185 as compared to NGC 147.

The central image of a strongly lensed background source places constraints on the foreground lens galaxy’s inner mass profile slope, core radius, and mass of its nuclear supermassive black hole. Using high-resolution long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations and archival Hubble Space Telescope (HST) imaging, we model the gravitational lens H-ATLAS J090311.6+003906 (also known as SDP.81) and search for the demagnified central image. There is central continuum emission from the lens galaxy’s active galactic nucleus (AGN) but no evidence of the central lensed image in any molecular line. We use the CO maps to determine the flux limit of the central image excluding the AGN continuum. We predict the flux density of the central image and use the limits from the ALMA data to constrain the innermost mass distribution of the lens. For a power-law profile with a core radius of 015 measured from HST photometry of the lens galaxy assuming that the central flux is attributed to the AGN, we find that a black hole mass of is preferred. Deeper observations with a detection of the central image will significantly improve the constraints of the innermost mass distribution of the lens galaxy.

Extremely large surveys with future experiments like Euclid and the SKA will soon allow us to access perturbation modes close to the Hubble scale, with wavenumbers If a modified gravity (MG) theory is responsible for cosmic acceleration, then the Hubble scale is a natural regime for deviations from General Relativity (GR) to become manifest. However, the majority of studies to date have concentrated on the consequences of alternative gravity theories for the subhorizon, quasi-static regime. In this paper, we investigate how modifications to the gravitational field equations affect perturbations around the Hubble scale. We choose functional forms to represent the generic scale-dependent behavior of gravity theories that modify GR at long wavelengths, and study the resulting deviations of ultra-large-scale relativistic observables from their GR behavior. We find that these are small unless modifications to the field equations are drastic. The angular dependence and redshift evolution of the deviations is highly parameterization- and survey-dependent, however, and so they are possibly a rich source of MG phenomenology if they can be measured.

The progenitor stars of several Type IIb supernovae (SNe) show indications of extended hydrogen envelopes. These envelopes might be the outcome of luminous energetic pre-explosion events, so-called precursor eruptions. We use the Palomar Transient Factory (PTF) pre-explosion observations of a sample of 27 nearby SNe IIb to look for such precursors during the final years prior to the SN explosion. No precursors are found when combining the observations in 15-day bins, and we calculate the absolute-magnitude-dependent upper limit on the precursor rate. At the 90% confidence level, SNe IIb have on average <0.86 precursors as bright as an absolute R-band magnitude of −14 in the final 3.5 years before the explosion and <0.56 events over the final year. In contrast, precursors among SNe IIn have a ≳5 times higher rate. The kinetic energy required to unbind a low-mass stellar envelope is comparable to the radiated energy of a few-weeks-long precursor that would be detectable for the closest SNe in our sample. Therefore, mass ejections, if they are common in such SNe, are radiatively inefficient or have durations longer than months. Indeed, when using 60-day bins, a faint precursor candidate is detected prior to SN 2012cs (∼2% false-alarm probability). We also report the detection of the progenitor of SN 2011dh that does not show detectable variability over the final two years before the explosion. The suggested progenitor of SN 2012P is still present, and hence is likely a compact star cluster or an unrelated object.

Using the Planck far-infrared and Arecibo GALFA 21 cm line surveys, we identified a set of isolated interstellar clouds (approximately degree-sized on the sky and comprising 100 solar masses) and assessed the ratio of gas mass to dust mass. Significant variations of the gas/dust ratio are found both from cloud to cloud and within regions of individual clouds; within the clouds, the atomic gas per unit dust decreases by more than a factor of 3 compared with the standard gas/dust ratio. Three hypotheses are considered. First, the apparently low gas/dust ratio could be due to molecular gas. Comparing to Planck CO maps, the brightest clouds have a H₂/CO ratio comparable to Galactic plane clouds, but a strong lower limit is placed on the ratio for other clouds, such that the required amount of molecular gas is far higher than would be expected based on the CO upper limits. Second, we consider self-absorbed 21 cm lines and find that the optical depth must be ∼3, significantly higher than found from surveys of radio sources. Third, grain properties may change within the clouds: they become more emissive when they are colder, while not utilizing heavy elements that already have their cosmic abundance fully locked into grains. It is possible that all three processes are active, and follow-up studies will be required to disentangle them and measure the true total gas and dust content of interstellar clouds.

The molecular carriers of the ubiquitous absorption features called the diffuse interstellar bands (DIBs) have eluded identification for many decades, in part because of the enormous parameter space spanned by the candidates and the limited set of empirical constraints afforded by observations in the diffuse interstellar medium. Detection of these features in circumstellar regions, where the environmental properties are more easily measured, is thus a promising approach to understanding the chemical nature of the carriers themselves. Here, using high-resolution spectra from the Apache Point Observatory Galactic Evolution Experiment survey, we present an analysis of the unusually asymmetric 1.5272 μm DIB feature along the sightline to the Red Square Nebula (RSN) and demonstrate the likely circumstellar origin of about half of the DIB absorption in this line of sight. This interpretation is supported both by the velocities of the feature components and by the ratio of foreground to total reddening along the line of sight. The RSN sightline offers the unique opportunity to study the behavior of DIB carriers in a constrained environment and thus to shed new light on the carriers themselves.

To understand the effects of cosmic-ray (CR) impacts on interstellar icy grains immersed in H₂ gas, we have irradiated porous water-ice films loaded with H₂ with 100 keV H⁺. The ice films were exposed to H₂ gas at different pressures following deposition and during irradiation. A net H₂ loss is observed during irradiation due to competition between ion-induced sputtering and gas adsorption. The initial H₂ loss cross-section, 4(1) × 10⁻¹⁴ cm², was independent of film thickness, H₂, and proton fluxes. In addition to sputtering, irradiation also closes nanopores, trapping H₂ in the film with binding that exceeds physical absorption energies. As a result, 2%–7% H₂ is retained in the ice following irradiation to high fluences. We find that the trapped H₂ concentration increases with decreasing Φ, the ratio of ion to H₂ fluxes, suggesting that as high as 8% solid H₂ can be trapped in interstellar ice by CR or stellar wind impacts.

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

We present full-orbit phase curve observations of the eccentric (e ∼ 0.08) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 μm bands using the Spitzer Space Telescope. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of 0.1882% ± 0.0048% and 0.2247% ± 0.0086% at 3.6 and 4.5 μm, respectively, are both consistent with a blackbody temperature of 2402 ± 35 K. We place a 2σ upper limit on the nightside flux at 3.6 μm and find it to be 9% ± 1% of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 μm, the minimum planet flux is 30% ± 5% of the maximum flux, corresponding to a brightness temperature of 1380 ± 65 K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b’s measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to overpredict the nightside emission at 3.6 μm, while underpredicting the nightside emission at 4.5 μm. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b (7.8 M_Jup) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

Palomar 5 (Pal 5) is a faint halo globular cluster associated with narrow tidal tails. It is a useful system to understand the process of tidal dissolution, as well as to constrain the potential of the Milky Way. A well-determined orbit for Pal 5 would enable detailed study of these open questions. We present here the first CCD-based proper motion measurement of Pal 5 obtained using SDSS as a first epoch and new Large Binocular Telescope/Large Binocular Camera (LBC) images as a second, giving a baseline of 15 years. We perform relative astrometry, using SDSS as a distortion-free reference, and images of the cluster and also of the Pal 5 stream for the derivation of the distortion correction for LBC. The reference frame is made up of background galaxies. We correct for differential chromatic refraction using relations obtained from SDSS colors as well as from flux-calibrated spectra, finding that the correction relations for stars and for galaxies are different. We obtain μ_α = −2.296 ± 0.186 mas yr⁻¹ and μ_δ = −2.257 ± 0.181 mas yr⁻¹ for the proper motion of Pal 5. We use this motion, and the publicly available code galpy, to model the disruption of Pal 5 in different Milky Way models consisting of a bulge, a disk, and a spherical dark matter halo. Our fits to the observed stream properties (streak and radial velocity gradient) result in a preference for a relatively large Pal 5 distance of around 24 kpc. A slightly larger absolute proper motion than what we measure also results in better matches but the best solutions need a change in distance. We find that a spherical Milky Way model, with V₀ = 220 km s⁻¹ and V_{20 kpc}, i.e., approximately at the apocenter of Pal 5, of 218 km s⁻¹, can match the data well, at least for our choice of disk and bulge parametrization.

We present spatially resolved imaging obtained with the Australia Telescope Compact Array (ATCA) of three CO lines in two high-redshift gravitationally lensed dusty star-forming galaxies, discovered by the South Pole Telescope. Strong lensing allows us to probe the structure and dynamics of the molecular gas in these two objects, at z = 2.78 and z = 5.66, with effective source-plane resolution of less than 1 kpc. We model the lensed emission from multiple CO transitions and the dust continuum in a consistent manner, finding that the cold molecular gas as traced by low-J CO always has a larger half-light radius than the 870 μm dust continuum emission. This size difference leads to up to 50% differences in the magnification factor for the cold gas compared to dust. In the z = 2.78 galaxy, these CO observations confirm that the background source is undergoing a major merger, while the velocity field of the other source is more complex. We use the ATCA CO observations and comparable resolution Atacama Large Millimeter/submillimeter Array dust continuum imaging of the same objects to constrain the CO–H₂ conversion factor with three different procedures, finding good agreement between the methods and values consistent with those found for rapidly star-forming systems. We discuss these galaxies in the context of the star formation—gas mass surface density relation, noting that the change in emitting area with observed CO transition must be accounted for when comparing high-redshift galaxies to their lower redshift counterparts.

We present the cross-correlation between the far-infrared (far-IR) background fluctuations as measured with the Herschel Space Observatory at 250, 350, and 500 μm and the near-infrared (near-IR) background fluctuations with the Spitzer Space Telescope at 3.6 and 4.5 μm. The cross-correlation between the FIR and NIR background anisotropies is detected such that the correlation coefficient at a few to 10 arcminute angular scale decreases from 0.3 to 0.1 when the FIR wavelength increases from 250 to 500 μm. We model the cross-correlation using a halo model with three components: (a) FIR bright or dusty star-forming galaxies below the masking depth in Herschel maps, (b) NIR faint galaxies below the masking depth, and (c) intra-halo light (IHL), or diffuse stars in dark matter halos, that is likely dominating the large-scale NIR fluctuations. The model is able to reasonably reproduce the auto-correlations at each of the FIR wavelengths and at 3.6 μm and their corresponding cross-correlations. While the FIR and NIR auto-correlations are dominated by faint, dusty, star-forming galaxies and IHL, respectively, we find that roughly half of the cross-correlation between the NIR and FIR backgrounds is due to the same dusty galaxies that remain unmasked at 3.6 μm. The remaining signal in the cross-correlation is due to IHL present in the same dark matter halos as those hosting the same faint and unmasked galaxies.

The Keck Array is a system of cosmic microwave background polarimeters, each similar to the Bicep2 experiment. In this paper we report results from the 2012 to 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as Bicep2. We again find an excess of B-mode power over the lensed-ΛCDM expectation of >5σ in the range 30 < ℓ < 150 and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with Bicep2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 μK arcmin) over an effective area of 400 deg² for an equivalent survey weight of 250,000 μK⁻². The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of >6σ.

The Mg ii h&k doublet are two of the primary spectral lines observed by the Sun-pointing Interface Region Imaging Spectrograph (IRIS). These lines are tracers of the magnetic and thermal environment that spans from the photosphere to the upper chromosphere. We use a double-Gaussian model to fit the Mg ii h profile for a full-Sun mosaic data set taken on 2014 August 24. We use the ensemble of high-quality profile fits to conduct a statistical study on the variability of the line profile as it relates the magnetic structure, dynamics, and center-to-limb viewing angle. The average internetwork profile contains a deeply reversed core and is weakly asymmetric at h2. In the internetwork, we find a strong correlation between h3 wavelength and profile asymmetry as well as h1 width and h2 width. The average reversal depth of the h3 core is inversely related to the magnetic field. Plage and sunspots exhibit many profiles that do not contain a reversal. These profiles also occur infrequently in the internetwork. We see indications of magnetically aligned structures in plage and network in statistics associated with the line core, but these structures are not clear or extended in the internetwork. The center-to-limb variations are compared to predictions of semi-empirical model atmospheres. We measure a pronounced limb darkening in the line core that is not predicted by the model. The aim of this work is to provide a comprehensive measurement baseline and preliminary analysis on the observed structure and formation of the Mg ii profiles observed by IRIS.

The Atmospheric Imaging Assembly (AIA) telescope on board the Solar Dynamics Observatory provides coronal extreme ultraviolet imaging over a broader temperature sensitivity range than the previous generations of instruments (Extreme Ultraviolet Imager; EUVI, EIT, and TRACE). Differential emission measure tomography (DEMT) of the solar corona based on AIA data is presented here for the first time. The main product of DEMT is the three-dimensional distribution of the local differential emission measure (LDEM). While in previous studies, based on EIT or EUVI data, there were three available EUV bands, the present study is based on the four cooler AIA bands (aimed at studying the quiet sun). The AIA filters allow exploration of new parametric LDEM models. Since DEMT is better suited for lower activity periods, we use data from Carrington Rotation 2099, when the Sun was in its most quiescent state during the AIA mission. Also, we validate the parametric LDEM models by using them to perform a bi-dimensional differential emission measure (DEM) analysis on sets of simultaneous AIA images, and comparing results with those obtained using other methods. Our study reveals a ubiquitous bimodal LDEM distribution in the quiet diffuse corona, characterized by two well-defined average centroid temperatures and We argue that the nanoflare heating scenario is less likely to explain these results, and that alternative mechanisms, such as wave dissipation, appear better supported by our results.

The million degree plasma of the solar corona must be supplied by the underlying layers of the atmosphere. The mechanism and location of energy release, and the precise source of coronal plasma, remain unresolved. In earlier work, we pursued the idea that warm plasma is supplied to the corona via direct heating of the chromosphere by nanoflares, contrary to the prevailing belief that the corona is heated in situ and the chromosphere is subsequently energized and ablated by thermal conduction. We found that single (low-frequency) chromospheric nanoflares could not explain the observed intensities, Doppler-shifts, and red/blue asymmetries in Fe xii and xiv emission lines. In the present work, we follow up on another suggestion that the corona could be powered by chromospheric nanoflares that repeat on a timescale substantially shorter than the cooling/draining timescale. That is, a single magnetic strand is re-supplied with coronal plasma before the existing plasma has time to cool and drain. We perform a series of hydrodynamic experiments and predict the Fe xii and xiv line intensities, Doppler-shifts, and red/blue asymmetries. We find that our predicted quantities disagree dramatically with observations and fully developed loop structures cannot be created by intermediate- or high-frequency chromospheric nanoflares. We conclude that the mechanism ultimately responsible for producing coronal plasma operates above the chromosphere, but this does not preclude the possibility of a similar mechanism powering the chromosphere, extreme examples of which may be responsible for heating chromospheric plasma to transition region temperatures (e.g., type II spicules).

The iron Kα line commonly observed in the X-ray spectrum of both stellar-mass and supermassive black hole (BH) candidates originates from X-ray fluorescence of the inner accretion disk. Accordingly, it can be used to map the spacetime geometry around these objects. In this paper, we extend previous work using the iron Kα line to test the Kerr BH hypothesis. We adopt the Cardoso–Pani–Rico parametrization and we test the possibility of constraining possible deviations from the Kerr solution that can be obtained from observations across the range of BH spins and inclination angles. We confirm previous claims that the iron Kα line is potentially a quite powerful probe for testing the Kerr metric given sufficiently high quality data and with systematics under control, especially in the case of fast-rotating BHs and high inclination angles since both conditions serve to maximize relativistic effects. We find that some geometric perturbations from Kerr geometry manifest more strongly in the iron line profile than others. While the perturbation parameter can be well constrained by the iron line profile, an orthogonal data set is necessary to constrain departures from Kerr geometry in

The transfer of matter between a circumbinary disk and a young binary system remains poorly understood, obscuring the interpretation of accretion indicators. To explore the behavior of these indicators in multiple systems, we have performed the first systematic time-domain study of young binaries in the ultraviolet. We obtained far- and near-ultraviolet HST/COS spectra of the young spectroscopic binaries DQ Tau and UZ Tau E. Here we focus on the continuum from 2800 to 3200 Å and on the C iv doublet (λλ1548.19, 1550.77 Å) as accretion diagnostics. Each system was observed over three or four consecutive binary orbits, at phases ∼0, 0.2, 0.5, and 0.7. Those observations are complemented by ground-based U-band measurements. Contrary to model predictions, we do not detect any clear correlation between accretion luminosity and phase. Further, we do not detect any correlation between C iv flux and phase. For both stars the appearance of the C iv line is similar to that of single Classical T Tauri Stars (CTTSs), despite the lack of stable long-lived circumstellar disks. However, unlike the case in single CTTSs, the narrow and broad components of the C iv lines are uncorrelated, and we argue that the narrow component is powered by processes other than accretion, such as flares in the stellar magnetospheres and/or enhanced activity in the upper atmosphere. We find that both stars contribute equally to the narrow component C iv flux in DQ Tau, but the primary dominates the narrow component C iv emission in UZ Tau E. The C iv broad component flux is correlated with other accretion indicators, suggesting an accretion origin. However, the line is blueshifted, which is inconsistent with its origin in an infall flow close to the star. It is possible that the complicated geometry of the region, as well as turbulence in the shock region, are responsible for the blueshifted line profiles.

We present a detailed analysis of a large-scale galactic outflow in the circumgalactic medium of a massive (), star-forming ( yr⁻¹), sub-L_* () galaxy at z = 0.39853 that exhibits a wealth of metal-line absorption in the spectra of the background quasar Q 0122−003 at an impact parameter of 163 kpc. The galaxy inclination angle () and the azimuthal angle () imply that the QSO sightline is passing through the projected minor-axis of the galaxy. The absorption system shows a multiphase, multicomponent structure with ultra-strong, wide velocity spread ( 419 km s⁻¹) and ( 285 km s⁻¹) lines that are extremely rare in the literature. The highly ionized absorption components are well explained as arising in a low density ( cm⁻³), diffuse (∼10 kpc), cool (∼10⁴ K) photoionized gas with a super-solar metallicity (). From the observed narrowness of the Lyβ profile, the non-detection of absorption, and the presence of strong absorption in the low-resolution FOS spectrum, we rule out equilibrium/non-equilibrium collisional ionization models. The low-ionization photoionized gas with a density of cm⁻³ and a metallicity of is possibly tracing recycled halo gas. We estimate an outflow mass of a mass-flow rate of a kinetic luminosity of erg s⁻¹, and a mass loading factor of ∼8 for the outflowing high-ionization gas. These are consistent with the properties of “down-the-barrel” outflows from infrared-luminous starbursts as studied by Rupke et al. Such powerful, large-scale, metal-rich outflows are the primary means of sufficient mechanical and chemical feedback as invoked in theoretical models of galaxy formation and evolution.

Near-infrared (NIR) spectra that have an angular resolution of ∼0.15 arcsec are used to examine the stellar content of the central regions of the nearby elliptical galaxy Maffei 1. The spectra were recorded at the Subaru Telescope, with wavefront distortions corrected by the RAVEN Multi-object Adaptive Optics science demonstrator. The Ballick–Ramsey C₂ absorption bandhead near 1.76 μm is detected, and models in which ∼10%–20% of the light near 1.8 μm originates from stars of spectral type C5 reproduce the depth of this feature. Archival NIR and mid-infrared images are also used to probe the structural and photometric properties of the galaxy. Comparisons with models suggest that an intermediate age population dominates the spectral energy distribution between 1 and 5 μm near the galaxy center. This is consistent not only with the presence of C stars, but also with the large Hβ index that has been measured previously for Maffei 1. The J − K color is more or less constant within 15 arcsec of the galaxy center, suggesting that the brightest red stars are well-mixed in this area.

We present observations of 1.5 square degree maps of the ¹²CO, ¹³CO, and C¹⁸O (J = 1 − 0) emission toward the complex region of the supernova remnant (SNR) W41 and SNR G22.7–0.2. A massive (), large (∼84 × 15 pc), and dense (∼10³ cm⁻³) giant molecular cloud (GMC), G23.0–0.4 with 77 km s⁻¹, is found to be adjacent to the two SNRs. The GMC displays a filamentary structure approximately along the Galactic plane. The filamentary structure of the dense molecular gas, traced by C¹⁸O (J = 1 − 0) emission, is also coincident well with the distribution of the dust-continuum emission in the direction. Two dense massive MC clumps, two 6.7 GHz methanol masers, and one H ii/SNR complex, associated with the 77 km s⁻¹ GMC G23.0–0.4, are aligned along the filamentary structure, indicating the star-forming activity within the GMC. These sources have periodic projected spacing of 018–026 along the giant filament, which is consistent with the theoretical predictions of 022. This indicates that the turbulence seems to dominate the fragmentation process of the dense gaseous filament on a large scale. The established 4.4 kpc distance of the GMC and the long dense filament traced by C¹⁸O emission, together with the rich massive star-formation groups in the nearby region, suggest that G23.0–0.4 is probably located at the near side of the Scutum–Centaurus arm in the first quadrant. Considering the large scale and the elongation structure along the Galactic plane, we speculate that the dense filamentary GMC is related to the spiral density wave of the Milky Way.

The generation of a large-scale magnetic field in the kinematic regime in the absence of an α-effect is investigated by following two different approaches: the test-field method and the multiscale stability theory relying on the homogenization technique. Our computations of the magnetic eddy diffusivity tensor of the parity-invariant flow IV of G. O. Roberts and the modified Taylor–Green flow confirm the findings of previous studies and also explain some of their apparent contradictions. The two flows have large symmetry groups; this is used to considerably simplify the eddy diffusivity tensor. Finally, a new analytic result is presented: upon expressing the eddy diffusivity tensor in terms of solutions to auxiliary problems for the adjoint operator, we derive relations between the magnetic eddy diffusivity tensors that arise for mutually reverse small-scale flows and .

We investigate several key questions of plasma heating in open-field regions of the corona that connect to the solar wind. We present results for a model of Alfvén-wave-driven turbulence for three typical open magnetic field structures: a polar coronal hole, an open flux tube neighboring an equatorial streamer, and an open flux tube near a strong-field active region. We compare time-steady, one-dimensional turbulent heating models against fully time-dependent three-dimensional reduced-magnetohydrodynamic modeling of BRAID. We find that the time-steady results agree well with time-averaged results from BRAID. The time dependence allows us to investigate the variability of the magnetic fluctuations and of the heating in the corona. The high-frequency tail of the power spectrum of fluctuations forms a power law whose exponent varies with height, and we discuss the possible physical explanation for this behavior. The variability in the heating rate is bursty and nanoflare-like in nature, and we analyze the amount of energy lost via dissipative heating in transient events throughout the simulation. The average energy in these events is 10^21.91 erg, within the “picoflare” range, and many events reach classical “nanoflare” energies. We also estimated the multithermal distribution of temperatures that would result from the heating-rate variability, and found good agreement with observed widths of coronal differential emission measure distributions. The results of the modeling presented in this paper provide compelling evidence that turbulent heating in the solar atmosphere by Alfvén waves accelerates the solar wind in open flux tubes.

Light bridges, the bright structures that divide the umbra of sunspots and pores into smaller pieces, are known to produce a wide variety of activity events in solar active regions (ARs). It is also known that the light bridges appear in the assembling process of nascent sunspots. The ultimate goal of this series of papers is to reveal the nature of light bridges in developing ARs and the occurrence of activity events associated with the light bridge structures from both observational and numerical approaches. In this first paper, exploiting the observational data obtained by Hinode, the Interface Region Imaging Spectrograph, and the Solar Dynamics Observatory, we investigate the detailed structure of the light bridge in NOAA AR 11974 and its dynamic activity phenomena. As a result, we find that the light bridge has a weak, horizontal magnetic field, which is transported from the interior by a large-scale convective upflow and is surrounded by strong, vertical fields of adjacent pores. In the chromosphere above the bridge, a transient brightening occurs repeatedly and intermittently, followed by a recurrent dark surge ejection into higher altitudes. Our analysis indicates that the brightening is the plasma heating due to magnetic reconnection at lower altitudes, while the dark surge is the cool, dense plasma ejected from the reconnection region. From the observational results, we conclude that the dynamic activity observed in a light bridge structure such as chromospheric brightenings and dark surge ejections are driven by magnetoconvective evolution within the light bridge and its interaction with the surrounding magnetic fields.

Light bridges, the bright structure dividing umbrae in sunspot regions, show various activity events. In Paper I, we reported on an analysis of multi-wavelength observations of a light bridge in a developing active region (AR) and concluded that the activity events are caused by magnetic reconnection driven by magnetconvective evolution. The aim of this second paper is to investigate the detailed magnetic and velocity structures and the formation mechanism of light bridges. For this purpose, we analyze numerical simulation data from a radiative magnetohydrodynamics model of an emerging AR. We find that a weakly magnetized plasma upflow in the near-surface layers of the convection zone is entrained between the emerging magnetic bundles that appear as pores at the solar surface. This convective upflow continuously transports horizontal fields to the surface layer and creates a light bridge structure. Due to the magnetic shear between the horizontal fields of the bridge and the vertical fields of the ambient pores, an elongated cusp-shaped current layer is formed above the bridge, which may be favorable for magnetic reconnection. The striking correspondence between the observational results of Paper I and the numerical results of this paper provides a consistent physical picture of light bridges. The dynamic activity phenomena occur as a natural result of the bridge formation and its convective nature, which has much in common with those of umbral dots and penumbral filaments.

With observations from the Interface Region Imaging Spectrograph, we track the complete evolution of ∼11 MK evaporation flows in an M1.1 flare on 2014 September 6 and an X1.6 flare on 2014 September 10. These hot flows, as indicated by the blueshifted Fe xxi 1354.08 Å line, evolve smoothly with a velocity decreasing exponentially from ∼200 km s⁻¹ to almost stationary within a few minutes. We find a good correlation between the flow velocity and energy deposition rate as represented by the hard X-ray flux observed with the Reuven Ramaty High Energy Solar Spectroscopic Imager, or time derivative of the soft X-ray flux observed with the Geostationary Operational Environmental Satellites and the HINODE X-ray Telescope, which is in general agreement with models of nonthermal electron heating. The maximum blueshift of Fe xxi appears approximately at the same time as or slightly after the impulsive enhancement of the ultraviolet continuum and the Mg ii 2798.8 Å line emission, demonstrating that the evaporation flow is closely related to heating of the lower chromosphere. Finally, while the hot Fe xxi 1354.08 Å line is entirely blueshifted with no obvious rest component, cool chromospheric and transition region lines like Si iv 1402.77 Å are often not entirely redshifted but just reveal an obvious red wing enhancement at the ribbons, suggesting that the speed of chromospheric condensation might be larger than previously thought.

Thomson optical depth τ measurements from Planck provide new insights into the reionization of the universe. In pursuit of model-independent constraints on the properties of the ionizing sources, we determine the empirical evolution of the cosmic ionizing emissivity. We use a simple two-parameter model to map out the evolution in the emissivity at z ≳ 6 from the new Planck optical depth τ measurements, from the constraints provided by quasar absorption spectra and from the prevalence of Lyα emission in z ∼ 7–8 galaxies. We find the redshift evolution in the emissivity required by the observations to be ( for a flat prior), largely independent of the assumed clumping factor C_{H ii} and entirely independent of the nature of the ionizing sources. The trend in is well-matched by the evolution of the galaxy UV-luminosity density () to a magnitude limit ≳−13 mag, suggesting that galaxies are the sources that drive the reionization of the universe. The role of galaxies is further strengthened by the conversion from the UV luminosity density ρ_UV to being possible for physically plausible values of the escape fraction f_esc, the Lyman-continuum photon production efficiency ξ_ion, and faint-end cut-off M_lim to the luminosity function. Quasars/active galactic nuclei appear to match neither the redshift evolution nor normalization of the ionizing emissivity. Based on the inferred evolution in the ionizing emissivity, we estimate that the z ∼ 10 UV-Iuminosity density is 8₋₄⁺¹⁵× lower than at z ∼ 6, consistent with the observations. The present approach of contrasting the inferred evolution of the ionizing emissivity with that of the galaxy UV luminosity density adds to the growing observational evidence that faint, star-forming galaxies drive the reionization of the universe.

Upcoming space-based surveys such as Euclid and WFIRST-AFTA plan to measure baryonic acoustic oscillations in order to study dark energy. These surveys will use IR slitless grism spectroscopy to measure redshifts of a large number of galaxies over a significant redshift range. In this paper, we use the Wide Field Camera 3 Infrared Spectroscopic Parallel Survey (WISP) to estimate the expected number of Hα emitters observable by these future surveys. WISP is an ongoing Hubble Space Telescope slitless spectroscopic survey, covering the 0.8–1.65 μm wavelength range and allowing the detection of Hα emitters up to z ∼ 1.5 and [O iii] emitters to z ∼ 2.3. We derive the Hα–[O iii] bivariate line luminosity function (LLF) for WISP galaxies at z ∼ 1 using a maximum likelihood estimator that properly accounts for uncertainties in line luminosity measurements and we demonstrate how it can be used to derive the Hα luminosity function by exclusively fitting [O iii] data. Using the [O iii] LLF and assuming that the relation between Hα and [O iii] luminosity does not change significantly over the redshift range, we predict the Hα number counts at —the upper end of the redshift range of interest for future surveys. For the redshift range we expect ∼3000 galaxies deg⁻² for a flux limit of 3 × 10⁻¹⁶ erg s⁻¹ cm⁻² (the proposed depth of the Euclid galaxy redshift survey) and ∼20,000 galaxies deg⁻² for a flux limit of ∼10⁻¹⁶ erg s⁻¹ cm⁻² (the baseline depth of the WFIRST galaxy redshift survey).

We report new observations of multiple transitions of the CS molecular lines in the SgrA region of the Galactic center, at an angular resolution of 40″ (=1.5 pc). The objective of this paper is to study the polar arc (PA), which is a molecular ridge near the SgrA region with apparent non-coplanar motions and a velocity gradient perpendicular to the Galactic plane. With our high resolution dense-gas maps, we search for the base and the origin of the PA, which is expected to be embedded in the Galactic disk. We find that the PA is connected to a continuous structure from one of the disk rings/arms in both the spatial and velocity domains. This structure near SgrA* has high CS(J = 4–3)/CS(J = 2–1) ratios The fact that this structure has eluded detection in previous observations is likely due to the combination of high excitation and low surface brightness temperature. We call this new structure the connecting ridge. We discuss the possible mechanism for forming this structure and for lifting the gas above the Galactic plane.

The well-studied blazar Markarian 421 (Mrk 421, z = 0.031) was the subject of an intensive multi-wavelength campaign when it flared in 2013 April. The recorded X-ray and very high-energy (E > 100 GeV) γ-ray fluxes are the highest ever measured from this object. At the peak of the activity, it was monitored by the hard X-ray focusing telescope Nuclear Spectroscopic Telescope Array (NuSTAR) and the Swift X-Ray Telescope (XRT). In this work, we present a detailed variability analysis of NuSTAR and Swift-XRT observations of Mrk 421 during this flaring episode. We obtained the shortest flux doubling time of 14.01 ± 5.03 minutes, which is the shortest hard X-ray (3–79 keV) variability ever recorded from Mrk 421, and is on the order of the light-crossing time of the black hole's event horizon. A pattern of extremely fast variability events superposed on slowly varying flares is found in most of the NuSTAR observations. We suggest that these peculiar variability patterns may be explained by magnetic energy dissipation and reconnection in a fast-moving compact emission region within the jet. Based on the fast variability, we derive a lower limit on the magnetic field strength of G, where δ₁ is the Doppler factor in units of 10, and ν₁₉ is the characteristic X-ray synchrotron frequency in units of 10¹⁹ Hz.

We calculate the effects of spot size on pulse profiles of moderately rotating neutron stars. Specifically, we quantify the bias introduced in radius measurements from the common assumption that spots are infinitesimally small. We find that this assumption is reasonable for spots smaller than 10°–18° and leads to errors that are ≤10% in the radius measurement, depending on the location of the spot and the inclination of the observer. We consider the implications of our results for neutron star radius measurements with the upcoming and planned X-ray missions NICER and LOFT. We calculate the expected spot size for different classes of sources and investigate the circumstances under which the assumption of a small spot is justified.

Using observations from the Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE) survey of the Magellanic Clouds (MC), we have found 35 evolved stars and stellar end products that are bright in the far-infrared. These 28 (LMC) and 7 (SMC) sources were selected from the 529 evolved star candidates in the HERITAGE far-infrared point source catalogs. Our source identification method is based on spectral confirmation, spectral energy distribution characteristics, careful examination of the multiwavelength images and includes constraints on the luminosity, resulting in a thoroughly vetted list of evolved stars. These sources span a wide range in luminosity and hence initial mass. We found 13 low- to intermediate-mass evolved stars, including asymptotic giant branch (AGB) stars, post-AGB stars, planetary nebulae, and a symbiotic star. We also identify 10 high mass stars, including 4 of the 15 known B[e] stars in the MC, 3 extreme red supergiants that are highly enshrouded by dust, a Luminous Blue Variable, a Wolf–Rayet star, and two supernova remnants. Further, we report the detection of 9 probable evolved objects which were previously undescribed in the literature. These sources are likely to be among the dustiest evolved objects in the MC. The Herschel emission may either be due to dust produced by the evolved star or it may arise from swept-up interstellar medium material.

Observations of nearby molecular clouds detect “shells,” which are likely caused by winds from young main sequence stars. However, the progenitors of these observed features are not well characterized and the mass-loss rates inferred from the gas kinematics are several orders of magnitude greater than those predicted by atomic line-driven stellar wind models. We use magnetohydrodynamic simulations to model winds launching within turbulent molecular clouds and explore the impact of wind properties on cloud morphology and turbulence. We find that winds do not produce clear features in turbulent statistics such as the Fourier spectra of density and momentum but do impact the Fourier velocity spectrum. The density and velocity distribution functions, especially as probed by CO spectral lines, strongly indicate the presence and influence of winds. We show that stellar mass-loss rates for individual stars must be yr⁻¹, similar to those estimated from observations, to reproduce shell properties. Consequently, we conclude that B and A-type main sequence stars have mass-loss rates several orders of magnitude larger that those predicted by models or that young stars are more variable than expected due to magnetic activity or accretion.

Hubble Space Telescope observations of the site of the supernova (SN) SN 2008ax obtained in 2011 and 2013 reveal that the possible progenitor object detected in pre-explosion images was in fact multiple. Four point sources are resolved in the new, higher-resolution images. We identify one of the sources with the fading SN. The other three objects are consistent with single supergiant stars. We conclude that their light contaminated the previously identified progenitor candidate. After subtraction of these stars, the progenitor appears to be significantly fainter and bluer than previously measured. Post-explosion photometry at the SN location indicates that the progenitor object has disappeared. If single, the progenitor is compatible with a supergiant star of B to mid-A spectral type, while a Wolf–Rayet (W-R) star would be too luminous in the ultraviolet to account for the observations. Moreover, our hydrodynamical modeling shows that the pre-explosion mass was 4–5 M_⊙ and the radius was 30–50 R_⊙, which is incompatible with a W-R progenitor. We present a possible interacting binary progenitor computed with our evolutionary models that reproduces all the observational evidence. A companion star as luminous as an O9–B0 main-sequence star may have remained after the explosion.

We explore the connections between the evolving galaxy and active galactic nucleus (AGN) populations. We present a simple phenomenological model that links the evolving galaxy mass function and the evolving quasar luminosity function, which makes specific and testable predictions for the distribution of host galaxy masses for AGNs of different luminosities. We show that the ϕ^∗ normalizations of the galaxy mass function and of the AGN luminosity function closely track each other over a wide range of redshifts, implying a constant “duty cycle” of AGN activity. The strong redshift evolution in the AGN L^∗ can be produced by either an evolution in the distribution of Eddington ratios, or in the m_bh/m_∗ mass ratio, or both. To try to break this degeneracy we look at the distribution of AGNs in the Sloan Digital Sky Survey (m_bh, L) plane, showing that an evolving ratio reproduces the observed data and also reproduces the local relations that connect the black hole population with the host galaxies for both quenched and star-forming populations. We stress that observational studies that compare the masses of black holes in active galaxies at high redshift with those in quiescent galaxies locally will always see much weaker evolution. Evolution of this form would produce, or could be produced by, a redshift-independent m_bh–σ relation and could explain why the local m_bh–σ relation is tighter than m_bh–m_* even if σ is not directly linked to black hole growth. Irrespective of the evolution of m_bh/m_*, the model reproduces both the appearance of “downsizing” and the so-called “sub-Eddington boundary” without any mass-dependence in the evolution of black hole growth rates.

We report on a sample of 48 nearby, star-forming galaxies observed with the Cosmic Origin Spectrograph on the Hubble Space Telescope. We measure the kinematics of warm gas in galactic outflows using a combination of four Si ii absorption lines. We use multi-wavelength ancillary data to estimate stellar masses (M_*), star formation rates (SFR), circular velocities (v_circ), and morphologies. The galaxies cover four orders of magnitude in M_* and SFR, and sample a wide range of morphologies from starbursting mergers to normal star-forming galaxies. We derive 3.0–3.5σ relations between outflow velocity and SFR, M_*, and v_circ. The outflow velocities scale as SFR^0.08−0.22, and with the range depending on whether we use a maximum or a central velocity to quantify the outflow velocity. After accounting for their increased SFR, mergers drive 32% faster outflows than non-merging galaxies, with all of the highest velocity outflows arising from mergers. Low-mass galaxies (log(M_*/ M_⊙) < 10.5) lose some low-ionization gas through galactic outflows, while more massive galaxies retain all of their low-ionization gas, unless they undergo a merger.

The Energetic Neutral Atom (ENA) full-sky maps obtained with the Interstellar Boundary Explorer (IBEX) show an unexpected bright narrow band of increased intensity. This so-called ENA ribbon results from charge exchange of interstellar neutral atoms with protons in the outer heliosphere or beyond. Among other hypotheses it has been argued that this ribbon may be related to a neutral density enhancement, or H-wave, in the local interstellar medium. Here we quantitatively demonstrate, on the basis of an analytical model of the principal large-scale heliospheric structure, that this scenario for the ribbon formation leads to results that are fully consistent with the observed location of the ribbon in the full-sky maps at all energies detected with high-energy sensor IBEX-Hi.

We investigate the effects produced mainly by broadband soft X-rays up to 2 keV (plus fast (∼keV) photoelectrons and low-energy (∼eV) induced secondary electrons) in the ice mixtures containing H₂O:CO₂:NH₃:SO₂ (10:1:1:1) at two different temperatures (50 and 90 K). The experiments are an attempt to simulate the photochemical processes induced by energetic photons in SO₂-containing ices present in cold environments in the ices surrounding young stellar objects (YSO) and in molecular clouds in the vicinity of star-forming regions, which are largely illuminated by soft X-rays. The measurements were performed using a high-vacuum portable chamber from the Laboratório de Astroquímica e Astrobiologia (LASA/UNIVAP) coupled to the spherical grating monochromator beamline at the Brazilian Synchrotron Light Source (LNLS) in Campinas, Brazil. In situ analyses were performed by a Fourier transform infrared spectrometer. Sample processing revealed the formation of several organic molecules, including nitriles, acids, and other compounds such as H₂O₂, H₃O⁺, SO₃, CO, and OCN⁻. The dissociation cross section of parental species was on the order of (2–7) × 10⁻¹⁸ cm². The ice temperature does not seem to affect the stability of SO₂ in the presence of X-rays. Formation cross sections of new species  produced were also determined. Molecular half-lives at ices toward YSOs due to the presence of incoming soft X-rays were estimated. The low values obtained employing two different models of the radiation field of YSOs (TW Hydra and typical T-Tauri star) reinforce that soft X-rays are indeed a very efficient source of molecular dissociation in such environments.

We present and test Tessellation-based Recovery of Amorphous halo Concentrations (TesseRACt), a non-parametric technique for recovering the concentration of simulated dark matter halos using Voronoi tessellation. TesseRACt is tested on idealized N-body halos that are axisymmetric, triaxial, or contain substructure and is compared to traditional least-squares fitting as well as to two non-parametric techniques that assume spherical symmetry. TesseRACt recovers halo concentrations within 3% of the true value regardless of whether the halo is spherical, axisymmetric, or triaxial. Traditional fitting and non-parametric techniques that assume spherical symmetry can return concentrations for non-spherical halos that are systematically off by as much as 10% from the true value. TesseRACt also performs significantly better when there is substructure present outside 0.5 R₂₀₀. Given that cosmological halos are rarely spherical and often contain substructure, we discuss implications for studies of halo concentration in cosmological N-body simulations including how choice of technique for measuring concentration might bias scaling relations.

We report new correlations between ratios of band intensities of the 15–20 μm emission bands of polycyclic aromatic hydrocarbons (PAHs) in a sample of 57 sources observed with the Spitzer/Infrared Spectrograph. This sample includes Large Magellanic Cloud point sources from the SAGE-Spec survey, nearby galaxies from the Spitzer Infrared Nearby Galaxies Survey survey, two Galactic interstellar medium cirrus sources, and the spectral maps of the Galactic reflection nebulae NGC 2023 and NGC 7023. We find that the 16.4, 17.4, and 17.8 μm band intensities are inter-correlated in all environments. In NGC 2023 and NGC 7023 these bands also correlate with the 11.0 and 12.7 μm band intensities. The 15.8 μm band correlates only with the 15–18 μm plateau and the 11.2 μm emission. We examine the spatial morphology of these bands and introduce radial cuts. We find that these bands can be spatially organized into three sets: the 12.7, 16.4, and 17.8 μm bands; the 11.2, 15.8 μm bands and the 15–18 μm plateau; and the 11.0 and 17.4 μm bands. We also find that the spatial distribution of the 12.7, 16.4, and 17.8 μm bands can be reconstructed by averaging the spatial distributions of the cationic 11.0 μm and neutral 11.2 μm bands. We conclude that the 17.4 μm band is dominated by cations, the 15.8 μm band by neutral species, and the 12.7, 16.4, and 17.8 μm bands by a combination of the two. These results highlight the importance of PAH ionization for spatially differentiating sub-populations by their 15–20 μm emission variability.

Recently, AMS-02 reported their results of cosmic ray (CR) observations. In addition to the AMS-02 data, we add HESS data to estimate the spectra of CR electrons and the diffuse gamma rays above TeV. In the conventional diffusion model, a global analysis is performed on the spectral features of CR electrons and the diffuse gamma rays by the GALRPOP package. The results show that the spectrum structure of the primary component of CR electrons cannot be fully reproduced by a simple power law and that the relevant break is around 100 GeV. At the 99% confidence level (C.L.) the injection indices above the break decrease from 2.54 to 2.35, but the ones below the break are only in the range of 2.746–2.751. The spectrum of CR electrons does not need to add TeV cutoff to also match the features of the HESS data. Based on the difference between the fluxes of CR electrons and their primary components, the predicted excess of CR positrons is consistent with the interpretation that these positrons originate from a pulsar or dark matter. In the analysis of the Galactic diffuse gamma rays with the indirect constraint of AMS-02 and HESS data, it is found that the fluxes of Galactic diffuse gamma rays are consistent with the GeV data of the Fermi-Large Area Telescope (LAT) in the high-latitude regions. The results indicate that inverse Compton scattering is the dominant component in the range of hundreds of GeV to tens of TeV, respectively from the high-latitude regions to the low ones, and in all of the regions of the Galaxy the flux of diffuse gamma rays is less than that of CR electrons at the energy scale of 20 TeV.

The dusty, ionized gas cloud G2 is currently passing the massive black hole in the Galactic Center at a distance of roughly 2400 Schwarzschild radii. We explore the possibility of a starting point of the cloud within the disks of young stars. We make use of the large amount of new observations in order to put constraints on G2's origin. Interpreting the observations as a diffuse cloud of gas, we employ three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations with the PLUTO code and do a detailed comparison with observational data. The simulations presented in this work update our previously obtained results in multiple ways: (1) high resolution three-dimensional hydrodynamical AMR simulations are used, (2) the cloud follows the updated orbit based on the Brackett-γ data, (3) a detailed comparison to the observed high-quality position–velocity (PV) diagrams and the evolution of the total Brackett-γ luminosity is done. We concentrate on two unsolved problems of the diffuse cloud scenario: the unphysical formation epoch only shortly before the first detection and the too steep Brackett-γ light curve obtained in simulations, whereas the observations indicate a constant Brackett-γ luminosity between 2004 and 2013. For a given atmosphere and cloud mass, we find a consistent model that can explain both, the observed Brackett-γ light curve and the PV diagrams of all epochs. Assuming initial pressure equilibrium with the atmosphere, this can be reached for a starting date earlier than roughly 1900, which is close to apo-center and well within the disks of young stars.

We calculate the abundances of electrons and ions in the hot (≳500 K), dusty parts of protoplanetary disks, treating for the first time the effects of thermionic and ion emission from the dust grains. High-temperature ionization modeling has involved simply assuming that alkali elements such as potassium occur as gas-phase atoms and are collisionally ionized following the Saha equation. We show that the Saha equation often does not hold, because free charges are produced by thermionic and ion emission and destroyed when they stick to grain surfaces. This means the ionization state depends not on the first ionization potential of the alkali atoms, but rather on the grains’ work functions. The charged species’ abundances typically rise abruptly above about 800 K, with little qualitative dependence on the work function, gas density, or dust-to-gas mass ratio. Applying our results, we find that protoplanetary disks’ dead zone, where high diffusivities stifle magnetorotational turbulence, has its inner edge located where the temperature exceeds a threshold value ≈1000 K. The threshold is set by ambipolar diffusion except at the highest densities, where it is set by Ohmic resistivity. We find that the disk gas can be diffusively loaded onto the stellar magnetosphere at temperatures below a similar threshold. We investigate whether the “short-circuit” instability of current sheets can operate in disks and find that it cannot, or works only in a narrow range of conditions; it appears not to be the chondrule formation mechanism. We also suggest that thermionic emission is important for determining the rate of Ohmic heating in hot Jupiters.

Binaries are typically excluded from direct imaging exoplanet surveys. However, the recent findings of Kepler and radial velocity programs show that planets can and do form in binary systems. Here, we suggest that visual binaries offer unique advantages for direct imaging. We show that Binary Differential Imaging (BDI), whereby two stars are imaged simultaneously at the same wavelength within the isoplanatic patch at a high Strehl ratio, offers improved point spread function (PSF) subtraction that can result in increased sensitivity to planets close to each star. We demonstrate this by observing a young visual binary separated by 4″ with MagAO/Clio-2 at 3.9 μm, where the Strehl ratio is high, the isoplanatic patch is large, and giant planets are bright. Comparing BDI to angular differential imaging (ADI), we find that BDI’s 5σ contrast is ∼0.5 mag better than ADI’s within ∼1″ for the particular binary we observed. Because planets typically reside close to their host stars, BDI is a promising technique for discovering exoplanets in stellar systems that are often ignored. BDI is also 2–4× more efficient than ADI and classical reference PSF subtraction, since planets can be detected around both the target and PSF reference simultaneously. We are currently exploiting this technique in a new MagAO survey for giant planets in 140 young nearby visual binaries. BDI on a space-based telescope would not be limited by isoplanatism effects and would therefore be an even more powerful tool for imaging and discovering planets.