Table of contents

Protogalactic environments are typically identified using quasar absorption lines and can manifest as Damped Lyman-alpha Absorbers (DLAs) and Lyman Limit Systems (LLSs). We use radio observations of Faraday effects to test whether these galactic building blocks host a magnetized medium, by combining DLA and LLS detections with 1.4 GHz polarization data from the NRAO VLA Sky Survey (NVSS). We obtain a control, a DLA, and an LLS sample consisting of 114, 19, and 27 lines of sight, respectively. Using a Bayesian framework and weakly informative priors, we are unable to detect either coherent or random magnetic fields in DLAs: the regular coherent fields must be $\leqslant 2.8$μG, and the lack of depolarization suggests the weakly magnetized gas in DLAs is non-turbulent and quiescent. However, we find a mild suggestive indication that LLSs have coherent magnetic fields, with a 71.5% probability that LLSs have higher $| \mathrm{RM}|$ than a control, although this is sensitive to the redshift distribution. We also find a strong indication that LLSs host random magnetic fields, with a 95.5% probability that LLS lines of sight have lower polarized fractions than a control. The regular coherent fields within the LLSs must be $\leqslant 2.4$μG, and the magnetized gas must be highly turbulent with a typical turbulent length scale on the order of ≈5–20 pc. Our results are consistent with the standard dynamo paradigm, whereby magnetism in protogalaxies increases in coherence over cosmic time, and with a hierarchical galaxy formation scenario, with the DLAs and LLSs exploring different stages of magnetic field evolution in galaxies.

There is good evidence that the centers of massive early-type galaxies have a bottom-heavy stellar initial mass function (IMF) compared to that of the Milky Way. Here we study the radial variation of the IMF within such galaxies, using a combination of high-quality Keck spectroscopy and a new suite of stellar population synthesis models that cover a wide range in metallicity. As in the previous studies in this series, the models are fitted directly to the spectra and treat all elemental abundance ratios as free parameters. Using newly obtained spectroscopy for six galaxies, including deep data extending to $\sim 1{R}_{{\rm{e}}}$ for the galaxies NGC 1407, NGC 1600, and NGC 2695, we find that the IMF varies strongly with galactocentric radius. For all six galaxies the IMF is bottom-heavy in the central regions, with average mass-to-light ratio "mismatch" parameter $\alpha \equiv {({\text{}}M/L)/({\text{}}M/L)}_{\mathrm{MW}}\approx 2.5$ at R = 0. The IMF rapidly becomes more bottom-light with increasing radius, flattening off near the Milky Way value ($\alpha \approx 1.1$) at $R\gt 0.4{R}_{{\rm{e}}}$. A consequence is that the luminosity-weighted average IMF depends on the measurement aperture: within $R={R}_{{\rm{e}}}$ we find $\langle \alpha {\rangle }_{L}=1.3\mbox{--}1.5$, consistent with recent lensing and dynamical results from SLACS and ${\mathrm{ATLAS}}^{3{\rm{D}}}$. Our results are also consistent with several earlier studies that were based on analyses of radial gradients of line indices. The observed IMF gradients support galaxy formation models in which the central regions of massive galaxies had a different formation history than their outer parts. Finally, we make use of the high signal-to-noise central spectra of NGC 1407 and NGC 2695 to demonstrate how we can disentangle IMF effects and abundance effects.

We present a framework to simultaneously constrain the values and uncertainties of the strength of convective core overshooting, metallicity, extinction, distance, and age in stellar populations. We then apply the framework to archival Hubble Space Telescope observations of six stellar clusters in the Large Magellanic Cloud that have reported ages between $\sim 1\mbox{--}2.5\,\mathrm{Gyr}$. Assuming a canonical value of the strength of core convective overshooting, we recover the well-known age–metallicity correlation, and additional correlations between metallicity and extinction and metallicity and distance. If we allow the strength of core overshooting to vary, we find that for intermediate-aged stellar clusters, the measured values of distance and extinction are negligibly effected by uncertainties of core overshooting strength. However, cluster age and metallicity may have disconcertingly large systematic shifts when ${{\rm{\Lambda }}}_{{\rm{c}}}$ is allowed to vary by more than $\pm 0.05\ {H}_{p}$. Using the six stellar clusters, we combine their posterior distribution functions to obtain the most probable core overshooting value, ${0.500}_{-0.134}^{+0.016}{H}_{p}$, which is in line with canonical values.

The Ca ii K spectroheliograms spanning over a century (1907–2007) from Kodaikanal Solar Observatory, India, have recently been digitized and calibrated. Applying a fully automated algorithm (which includes contrast enhancement and the "Watershed method") to these data, we have identified the supergranules and calculated the associated parameters, such as scale, circularity, and fractal dimension. We have segregated the quiet and active regions and obtained the supergranule parameters separately for these two domains. In this way, we have isolated the effect of large-scale and small-scale magnetic fields on these structures and find a significantly different behavior of the supergranule parameters over solar cycles. These differences indicate intrinsic changes in the physical mechanism behind the generation and evolution of supergranules in the presence of small-scale and large-scale magnetic fields. This also highlights the need for further studies using solar dynamo theory along with magneto-convection models.

We investigate the occurrence of radio minihalos—diffuse radio sources of unknown origin observed in the cores of some galaxy clusters—in a statistical sample of 58 clusters drawn from the Planck Sunyaev–Zel'dovich cluster catalog using a mass cut (M₅₀₀ > 6 × 10¹⁴M_⊙). We supplement our statistical sample with a similarly sized nonstatistical sample mostly consisting of clusters in the ACCEPT X-ray catalog with suitable X-ray and radio data, which includes lower-mass clusters. Where necessary (for nine clusters), we reanalyzed the Very Large Array archival radio data to determine whether a minihalo is present. Our total sample includes all 28 currently known and recently discovered radio minihalos, including six candidates. We classify clusters as cool-core or non-cool-core according to the value of the specific entropy floor in the cluster center, rederived or newly derived from the Chandra X-ray density and temperature profiles where necessary (for 27 clusters). Contrary to the common wisdom that minihalos are rare, we find that almost all cool cores—at least 12 out of 15 (80%)—in our complete sample of massive clusters exhibit minihalos. The supplementary sample shows that the occurrence of minihalos may be lower in lower-mass cool-core clusters. No minihalos are found in non-cool cores or "warm cores." These findings will help test theories of the origin of minihalos and provide information on the physical processes and energetics of the cluster cores.

Observations and modeling suggest that dust abundance (gas-to-dust ratio, G/D) depends on (surface) density. Variations of the G/D provide timescale constraints for the different processes involved in the life cycle of metals in galaxies. Recent G/D measurements based on Herschel data suggest a factor of 5–10 decrease in dust abundance between the dense and diffuse interstellar media (ISM) in the Magellanic Clouds. However, the relative nature of the Herschel measurements precludes definitive conclusions as to the magnitude of those variations. We investigate variations of the dust abundance in the LMC and SMC using all-sky far-infrared surveys, which do not suffer from the limitations of Herschel on their zero-point calibration. We stack the dust spectral energy distribution (SED) at 100, 350, 550, and 850 microns from IRAS and Planck in intervals of gas surface density, model the stacked SEDs to derive the dust surface density, and constrain the relation between G/D and gas surface density in the range 10–100 M_⊙ pc⁻² on ∼80 pc scales. We find that G/D decreases by factors of 3 (from 1500 to 500) in the LMC and 7 (from $1.5\times {10}^{4}$ to 2000) in the SMC between the diffuse and dense ISM. The surface-density-dependence of G/D is consistent with elemental depletions, and with simple modeling of the accretion of gas-phase metals onto dust grains. This result has important implications for the sub-grid modeling of galaxy evolution, and for the calibration of dust-based gas-mass estimates, both locally and at high redshift.

The Galex Nearby Young Star Survey (GALNYSS) has yielded a sample of ∼2000 UV-selected objects that are candidate nearby ($D\lesssim 150\,\mathrm{pc}$), young (age ∼ 10–100 Myr), late-type stars. Here, we evaluate the distances and ages of the subsample of (19) GALNYSS stars with Gaia Data Release 1 (DR1) parallax distances $D\leqslant 120\,\mathrm{pc}$. The overall youth of these 19 mid-K to early-M stars is readily apparent from their positions relative to the loci of main-sequence stars and giants in Gaia-based color-magnitude and color-color diagrams constructed for all stars detected by Galex and the Wide-field Infrared Space Explorer for which parallax measurements are included in DR1. The isochronal ages of all 19 stars lie in the range ∼10–100 Myr. Comparison with Li-based age estimates indicates a handful of these stars may be young main-sequence binaries rather than pre-main sequence stars. Nine of the 19 objects have not previously been considered as nearby, young stars, and all but one of these are found at declinations north of +30°. The Gaia DR1 results presented here indicate that the GALNYSS sample includes several hundred nearby, young stars, a substantial fraction of which have not been previously recognized as having ages $\lesssim 100\,\mathrm{Myr}$.

We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions (MEs) that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid ME model, with two different basis sets and a thick-disk correction, accurately reproduces the overall structures of the Milky Way. Through non-axisymmetric Galactic structures of an elongated bar and spiral arms, gas clouds in the disk inflow to the nuclear region and form a central molecular zone-like nuclear ring. We find that the size of the nuclear ring evolves into $\sim 240\,\mathrm{pc}$ at $T\sim 1500\,\mathrm{Myr}$, regardless of the initial size. For most simulation runs, the rate of gas inflow to the nuclear region is equilibrated to $\sim 0.02\,{M}_{\odot }\,{\mathrm{yr}}^{-1}$. The nuclear ring is off-centered, relative to the Galactic center, by the lopsided central mass distribution of the Galaxy model, and thus an asymmetric mass distribution of the nuclear ring arises accordingly. The vertical asymmetry of the Galaxy model also causes the nuclear ring to be tilted along the Galactic plane. During the first ∼100 Myr, the vertical frequency of the gas motion is twice that of the orbital frequency, thus the projected nuclear ring shows a twisted, $\infty$-like shape.

We report the analysis of the first resolved caustic-crossing binary-source microlensing event OGLE-2016-BLG-1003. The event is densely covered by round-the-clock observations of three surveys. The light curve is characterized by two nested caustic-crossing features, which is unusual for typical caustic-crossing perturbations. From the modeling of the light curve, we find that the anomaly is produced by a binary source passing over a caustic formed by a binary lens. The result proves the importance of high-cadence and continuous observations, and the capability of second-generation microlensing experiments to identify such complex perturbations that are previously unknown. However, the result also raises the issues of the limitations of current analysis techniques for understanding lens systems beyond two masses and of determining the appropriate multiband observing strategy of survey experiments.

The range of currently proposed active galactic nucleus (AGN) far-infrared templates results in uncertainties in retrieving host galaxy information from infrared observations and also undermines constraints on the outer part of the AGN torus. We discuss how to test and reconcile these templates. Physically, the fraction of the intrinsic AGN IR-processed luminosity compared with that from the central engine should be consistent with the dust-covering factor. In addition, besides reproducing the composite spectral energy distributions (SEDs) of quasars, a correct AGN IR template combined with an accurate library of star-forming galaxy templates should be able to reproduce the IR properties of the host galaxies, such as the luminosity-dependent SED shapes and aromatic feature strengths. We develop tests based on these expected behaviors and find that the shape of the AGN intrinsic far-IR emission drops off rapidly starting at ∼20 μm and can be matched by an Elvis et al.-like template with a minor modification. Despite the variations in the near- to mid-IR bands, AGNs in quasars and Seyfert galaxies have remarkably similar intrinsic far-IR SEDs at λ ∼ 20–100 μm, suggesting a similar emission character of the outermost region of the circumnuclear torus. The variations of the intrinsic AGN IR SEDs among the type-1 quasar population can be explained by the changing relative strengths of four major dust components with similar characteristic temperatures, and there is evidence for compact AGN-heated dusty structures at sub-kiloparsec scales in the far-IR.

We consider the formation of binary black hole (BH) mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a BH binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov–Kozai cycles) with stellar evolution. After a BH triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most BH mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a BH merger rate in the range (0.3–1.3) Gpc⁻³ yr⁻¹, or up to ≈2.5 Gpc⁻³ yr⁻¹ if the BH orbital planes have initially random orientation. Finally, we show that BH mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ≈10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.

The broadband emission of pulsar wind nebulae (PWNe) is well described by non-thermal emissions from accelerated electrons and positrons. However, the standard shock acceleration model of PWNe does not account for the hard spectrum in radio wavelengths. The origin of the radio-emitting particles is also important to determine the pair production efficiency in the pulsar magnetosphere. Here, we propose a possible resolution for the particle energy distribution in PWNe; the radio-emitting particles are not accelerated at the pulsar wind termination shock but are stochastically accelerated by turbulence inside PWNe. We upgrade our past one-zone spectral evolution model to include the energy diffusion, i.e., the stochastic acceleration, and apply the model to the Crab Nebula. A fairly simple form of the energy diffusion coefficient is assumed for this demonstrative study. For a particle injection to the stochastic acceleration process, we consider the continuous injection from the supernova ejecta or the impulsive injection associated with supernova explosion. The observed broadband spectrum and the decay of the radio flux are reproduced by tuning the amount of the particle injected to the stochastic acceleration process. The acceleration timescale and the duration of the acceleration are required to be a few decades and a few hundred years, respectively. Our results imply that some unveiled mechanisms, such as back reaction to the turbulence, are required to make the energies of stochastically and shock-accelerated particles comparable.

SS 433 is an X-ray binary and the source of sub-relativistic, precessing, baryonic jets. We present high-resolution spectrograms of SS 433 in the infrared H and K bands. The spectrum is dominated by hydrogen and helium emission lines. The precession phase of the emission lines from the jet continues to be described by a constant period, ${P}_{\mathrm{jet}}=162.375\,\mathrm{days}$. The limit on any secularly changing period is $| \dot{P}| \lesssim {10}^{-5}$. The He i $\lambda 2.0587\,\mu {\rm{m}}$ line has complex and variable P-Cygni absorption features produced by an inhomogeneous wind with a maximum outflow velocity near 900 km s⁻¹. The He ii emission lines in the spectrum also arise in this wind. The higher members of the hydrogen Brackett lines show a double-peaked profile with symmetric wings extending more than ±1500 km s⁻¹ from the line center. The lines display radial velocity variations in phase with the radial velocity variation expected of the compact star, and they show a distortion during disk eclipse that we interpret as a rotational distortion. We fit the line profiles with a model in which the emission comes from the surface of a symmetric, Keplerian accretion disk around the compact object. The outer edge of the disk has velocities that vary from 110 to 190 km s⁻¹. These comparatively low velocities place an important constraint on the mass of the compact star: its mass must be less than $2.2\,{M}_{\odot }$ and is probably less than $1.6\,{M}_{\odot }$.

The relativistic jets created by some active galactic nuclei are important agents of AGN feedback. In spite of this, our understanding of what produces these jets is still incomplete. X-ray observations, which can probe the processes operating in the central regions in the immediate vicinity of the supermassive black hole, the presumed jet launching point, are potentially particularly valuable in illuminating the jet formation process. Here, we present the hard X-ray NuSTAR observations of the radio-loud quasar 4C 74.26 in a joint analysis with quasi-simultaneous, soft X-ray Swift observations. Our spectral analysis reveals a high-energy cutoff of ${183}_{-35}^{+51}$ keV and confirms the presence of ionized reflection in the source. From the average spectrum we detect that the accretion disk is mildly recessed, with an inner radius of R_in = 4–180 R_g. However, no significant evolution of the inner radius is seen during the three months covered by our NuSTAR campaign. This lack of variation could mean that the jet formation in this radio-loud quasar differs from what is observed in broad-line radio galaxies.

We report the results of three Chandra observations covering most of the extent of the TeV γ-ray source HESS J1616–508 and a search for a lower-energy counterpart to this source. We detect 56 X-ray sources, 37 of which have counterparts at lower frequencies, including a young massive star cluster, but none of them appear to be a particularly promising counterpart to the TeV source. The brightest X-ray source, CXOU J161423.4–505738, with a flux F_{0.5–7 keV} ≈ 5 × 10⁻¹³ erg cm⁻² s⁻¹, has a hard spectrum that is well fit by a power-law model with a photon index Γ = 0.2 ± 0.3 and is a likely intermediate polar CV candidate. No counterparts of this source were detected at other wavelengths. CVs are not known to produce extended TeV emission, and the source is also largely offset (19') from HESS J1616–508, making them unlikely to be associated. We have also set an upper limit on the X-ray flux of PSR J1614–5048 in the 0.5–8 keV band (F_{0.5–8 keV} < 5 × 10⁻¹⁵ erg cm⁻² s⁻¹ at a 90% confidence level). This makes PSR J1614–5048 one of the least X-ray-efficient pulsars known, with an X-ray efficiency ${\eta }_{0.5\mbox{--}8\mathrm{keV}}={L}_{0.5\mbox{--}8\mathrm{keV}}/\dot{E}\lt 2\times {10}^{-5}$. We find no evidence supporting the association between the pulsar and the TeV source. We rule out a number of X-ray sources as possible counterparts to the TeV emission and do not find a plausible counterpart among the other sources. Lastly, we discuss the possible relation of PSR J1617–5055 to HESS J1616–508 in light of the new observations.

Star formation from the interstellar medium of galactic disks is a basic process controlling the evolution of galaxies. Understanding the star formation rate (SFR) in a local patch of a disk with a given gas mass is thus an important challenge for theoretical models. Here we simulate a kiloparsec region of a disk, following the evolution of self-gravitating molecular clouds down to subparsec scales, as they form stars that then inject feedback energy by dissociating and ionizing UV photons and supernova explosions. We assess the relative importance of each feedback mechanism. We find that H₂-dissociating feedback results in the largest absolute reduction in star formation compared to the run with no feedback. Subsequently adding photoionization feedback produces a more modest reduction. Our fiducial models that combine all three feedback mechanisms yield, without fine-tuning, SFRs that are in excellent agreement with observations, with H₂-dissociating photons playing a crucial role. Models that only include supernova feedback—a common method in galaxy evolution simulations—settle to similar SFRs, but with very different temperatures and chemical states of the gas, and with very different spatial distributions of young stars.

We study winds in 12 X-ray AGN host galaxies at $z\sim 1$. We find, using the low-ionization Fe iiλ2586 absorption in the stacked spectra, that the probability distribution function (PDF) of the centroid-velocity shift in AGNs has 50th (median), 16th, and 84th percentiles of (−87, −251, +86) km s⁻¹ respectively. The PDF of the velocity dispersion in AGNs has 50th (median), 84th, and 16th percentiles of (139, 253, 52) km s⁻¹ respectively. The centroid velocity and the velocity dispersions are obtained from a two-component (ISM+wind) absorption-line model. The equivalent width PDF of the outflow in AGNs has 50th (median), 84th, and 16th percentiles of (0.4, 0.8, 0.1) Å. There is a strong ISM component in Fe ii$\lambda 2586$ absorption (with (1.2, 1.5, 0.8) Å, implying the presence of a substantial amount cold gas in the host galaxies. For comparison, star-forming and X-ray undetected galaxies at a similar redshift, matched roughly in stellar mass and galaxy inclination, have a centroid-velocity PDF with percentiles of (−74, −258, +90) km s⁻¹, and a velocity-dispersion PDF with percentiles of (150, 259, 57) km s⁻¹. Thus, winds in the AGN are similar to star formation-driven winds, and are too weak to escape and expel substantial cool gas from galaxies. Our sample doubles the previous sample of AGNs studied at $z\sim 0.5$ and extends the analysis to $z\sim 1$. A joint reanalysis of the $z\sim 0.5$ AGN sample and our sample yields consistent results to the measurements above.

We present new theoretical period–luminosity–metallicity (PLZ) relations for RR Lyræ stars (RRLs) at Spitzer and WISE wavelengths. The PLZ relations were derived using nonlinear, time-dependent convective hydrodynamical models for a broad range of metal abundances (Z = 0.0001–0.0198). In deriving the light curves, we tested two sets of atmospheric models and found no significant difference between the resulting mean magnitudes. We also compare our theoretical relations to empirical relations derived from RRLs in both the field and in the globular cluster M4. Our theoretical PLZ relations were combined with multi-wavelength observations to simultaneously fit the distance modulus, μ₀, and extinction, A_V, of both the individual Galactic RRL and of the cluster M4. The results for the Galactic RRL are consistent with trigonometric parallax measurements from Gaia's first data release. For M4, we find a distance modulus of μ₀ = 11.257 ± 0.035 mag with A_V = 1.45 ± 0.12 mag, which is consistent with measurements from other distance indicators. This analysis has shown that, when considering a sample covering a range of iron abundances, the metallicity spread introduces a dispersion in the PL relation on the order of 0.13 mag. However, if this metallicity component is accounted for in a PLZ relation, the dispersion is reduced to ∼0.02 mag at mid-infrared wavelengths.

The thermal plasma beta in the solar wind and the solar corona is of the order of $\beta \sim 1$ and $\beta \ll 1$. Zank et al. developed 2D and slab turbulence transport model equations of the order of $\beta \sim 1$ and $\beta \ll 1$ using nearly incompressible (NI) theory. We solve the Zank et al. NI MHD coupled turbulence transport equations for the inhomogeneous solar wind from 1 to 75 au, and compare the numerical solutions to Voyager 2 observations. We find that (1) the 2D turbulent energies are larger than the slab energies throughout the heliosphere; (2) the 2D turbulent energies decrease more slowly than the slab turbulent energies within ∼4 au, while the slab energies increase and the 2D energies flatten in the outer heliosphere; (3) the 2D normalized cross-helicity decreases faster than the slab normalized cross-helicity within ∼4 au; (4) the 2D normalized residual energy is more magnetically dominated than the slab; (5) the variance of density fluctuations decreases more rapidly than ${r}^{-4}$ within ∼10 au, and more slowly in the outer heliosphere; and (6) the observed variance in magnetic field fluctuations as a function of the thermal plasma beta is described by the two-component turbulence transport model. In summary, the NI MHD two-component Zank et al. turbulence transport model captures the behavior of the forward, backward, and total energies in the fluctuating Elsässer variables, the variance in the magnetic field, kinetic energy, and density fluctuations, the cross-helicities and residual energies, the thermal temperature and plasma beta, and the various correlation lengths.

Planets are thought to form via accretion from a remnant disk of gas and solids around a newly formed star. During this process, material in the disk either remains bound to the star as part of either a planet, a smaller celestial body, or makes up part of the the interplanetary medium; falls into the star; or is ejected from the system. Herein we use dynamical models to probe the abundance and properties of ejected material during late-stage planet formation and estimate their contribution to the free-floating planet population. We present 300 N-body simulations of terrestrial planet formation around a solar-type star, with and without giant planets present, using a model that accounts for collisional fragmentation. In simulations with Jupiter and Saturn analogs, about one-third of the initial (∼5 M_⊕) disk mass is ejected, about half in planets more massive than Mercury but with a mass lower than 0.3 M_⊕, and the remainder in smaller bodies. Most ejections occur within 25 Myr, which is shorter than the timescale typically required for Earth-mass planets to grow (30–100 Myr). When giant planets are omitted from our simulations, almost no material is ejected within 200 Myr and only about 1% of the initial disk is ejected by 2 Gyr. We show that about 2.5 terrestrial-mass planets are ejected per star in the Galaxy. We predict that the space-borne microlensing search for free-floating planets from the Wide-Field Infra-Red Space Telescope will discover up to 15 Mars-mass planets, but few free-floating Earth-mass planets.

We have run a new suite of simulations that solve hydrodynamics and radiative transfer simultaneously to study helium ii reionization. Our suite of simulations employs various models for populating quasars inside of dark matter halos, which affect the He ii reionization history. In particular, we are able to explore the impact that differences in the timing and duration of reionization have on observables. We examine the thermal signature that reionization leaves on the intergalactic medium (IGM), and measure the temperature-density relation. As previous studies have shown, we confirm that the photoheating feedback from helium ii reionization raises the temperature of the IGM by several thousand kelvin. To compare against observations, we generate synthetic Lyα forest sightlines on-the-fly and match the observed effective optical depth ${\tau }_{\mathrm{eff}}(z)$ of hydrogen to recent observations. We show that when the simulations have been normalized to have the same values of ${\tau }_{\mathrm{eff}}$, the effect that helium ii reionization has on observations of the hydrogen Lyα forest is minimal. Specifically, the flux PDF and the one-dimensional power spectrum are sensitive to the thermal state of the IGM, but do not show direct evidence for the ionization state of helium. We show that the peak temperature of the IGM typically corresponds to the time of 90%–95% helium ionization by volume, and is a relatively robust indicator of the timing of reionization. Future observations of helium reionization from the hydrogen Lyα forest should thus focus on measuring the temperature of the IGM, especially at mean density. Detecting the peak in the IGM temperature would provide valuable information about the timing of the end of helium ii reionization.

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three-dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation subgrid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold and (2) "Magnetically Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between the non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial substructure and more disturbed kinematics.

We present the results of the search for gravitational waves (GWs) associated with γ-ray bursts detected during the first observing run of the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO). We find no evidence of a GW signal for any of the 41 γ-ray bursts for which LIGO data are available with sufficient duration. For all γ-ray bursts, we place lower bounds on the distance to the source using the optimistic assumption that GWs with an energy of ${10}^{-2}{M}_{\odot }{c}^{2}$ were emitted within the $16$–$500$ Hz band, and we find a median 90% confidence limit of 71 Mpc at 150 Hz. For the subset of 19 short/hard γ-ray bursts, we place lower bounds on distance with a median 90% confidence limit of 90 Mpc for binary neutron star (BNS) coalescences, and 150 and 139 Mpc for neutron star–black hole coalescences with spins aligned to the orbital angular momentum and in a generic configuration, respectively. These are the highest distance limits ever achieved by GW searches. We also discuss in detail the results of the search for GWs associated with GRB 150906B, an event that was localized by the InterPlanetary Network near the local galaxy NGC 3313, which is at a luminosity distance of $54$ Mpc (z = 0.0124). Assuming the γ-ray emission is beamed with a jet half-opening angle $\leqslant 30^\circ$, we exclude a BNS and a neutron star–black hole in NGC 3313 as the progenitor of this event with confidence >99%. Further, we exclude such progenitors up to a distance of 102 Mpc and 170 Mpc, respectively.

Proper motions of collisionless, pointlike objects in a spherically symmetric system—for example, stars in a galaxy—can be used to test whether that system is in equilibrium, with no assumptions regarding isotropy. In particular, the fourth-order spherical Jeans equation yields expressions for two observable quantities characterizing the departure from equilibrium, both of which can be expressed in terms of time derivatives of first and third moments of the velocities. As illustrations, we compute these quantities for tracer distributions drawn from an exact equilibrium configuration, and also from near-equilibrium configurations generated using the N-body code GALIC.

Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the three-dimensional velocity dispersions $({\sigma }_{r},{\sigma }_{\theta },{\sigma }_{\phi })$ or by the velocity anisotropy parameter $\beta =1-({\sigma }_{\theta }^{2}+{\sigma }_{\phi }^{2})/(2{\sigma }_{r}^{2})$. Direct measurements of $({\sigma }_{r},{\sigma }_{\theta },{\sigma }_{\phi })$ consistently suggest β = 0.5–0.7 for nearby halo stars. In contrast, the value of β at large Galactocentric radius r is still controversial, since reliable proper motion data are available for only a handful of stars. In the last decade, several authors have tried to estimate β for distant halo stars by fitting the observed line-of-sight velocities at each radius with simple velocity distribution models (local fitting methods). Some results of local fitting methods imply $\beta \lt 0$ at $r\gtrsim 20\,\mathrm{kpc}$, which is inconsistent with recent predictions from cosmological simulations. Here we perform mock-catalog analyses to show that the estimates of β based on local fitting methods are reliable only at $r\leqslant 15\,\mathrm{kpc}$ with the current sample size (∼10³ stars at a given radius). As r increases, the line-of-sight velocity (corrected for the solar reflex motion) becomes increasingly closer to the Galactocentric radial velocity, so it becomes increasingly more difficult to estimate the tangential velocity dispersion $({\sigma }_{\theta },{\sigma }_{\phi })$ from the line-of-sight velocity distribution. Our results suggest that the forthcoming Gaia data will be crucial for understanding the velocity distribution of halo stars at $r\geqslant 20\,\mathrm{kpc}$.

We present a study of the effective (half-light) radii and other structural properties of a systematically selected sample of young, massive star clusters (≥5 × 10³${M}_{\odot }$ and ≤200 Myr) in two nearby spiral galaxies, NGC 628 and NGC 1313. We use Hubble Space Telescope (HST) WFC3/UVIS and archival ACS/WFC data obtained by the Legacy Extragalactic UV Survey (LEGUS), an HST Treasury Program. We measure effective radii with GALFIT, a two-dimensional image-fitting package, and with a new technique to estimate effective radii from the concentration index of observed clusters. The distribution of effective radii from both techniques spans ∼0.5–10 pc and peaks at 2–3 pc for both galaxies. We find slight positive correlations between effective radius and cluster age in both galaxies, but no significant relationship between effective radius and galactocentric distance. Clusters in NGC 1313 display a mild increase in effective radius with cluster mass, but the trend disappears when the sample is divided into age bins. We show that the vast majority of the clusters in both galaxies are much older than their dynamical times, suggesting they are gravitationally bound objects. We find that about half of the clusters in NGC 628 are underfilling their Roche lobes, based on their Jacobi radii. Our results suggest that the young, massive clusters in NGC 628 and NGC 1313 are expanding, due to stellar mass loss or two-body relaxation, and are not significantly influenced by the tidal fields of their host galaxies.

We conduct a systematic search for galaxies at $z=0.1\mbox{--}1.5$ with [O ii]$\lambda 3727$, [O iii]$\lambda 5007$, or Hα$\lambda 6563$ emission lines extended over at least 30 kpc by using deep narrowband and broadband imaging in the Subaru-XMM Deep Survey field. These extended emission-line galaxies are dubbed [O ii], [O iii], or Hα blobs. Based on a new selection method that securely selects extended emission-line galaxies, we find 77 blobs at $z=0.40\mbox{--}1.46$ with the isophotal area of emission lines down to $1.2\times {10}^{-18}$ erg s⁻¹ cm⁻² kpc⁻². Four of them are spectroscopically confirmed to be [O iii] blobs at z = 0.83. We identify AGN activities in eight blobs with X-ray and radio data, and find that the fraction of AGN contribution increases with increasing isophotal area of the extended emission. With the Kolmogorov–Smirnov (KS) and Anderson–Darling tests, we confirm that the stellar-mass distributions of Hα and [O ii] blobs are not drawn from those of the emitters at the $\gt 90$% confidence level in that Hα and [O ii] blobs are located at the massive end of the distributions, but cannot reject a null hypothesis of being the same distributions in terms of the specific star formation rates. It is suggested that galactic-scale outflows tend to be more prominent in more massive star-forming galaxies. Exploiting our sample homogeneously selected over the large area, we derive the number densities of blobs at each epoch. The number densities of blobs decrease drastically with redshifts at a rate that is larger than that of the decrease of cosmic star formation densities.

We discuss the nature of the small areas of rapidly diverging, open magnetic flux that form in the strong unipolar fields at the peripheries of active regions (ARs), according to coronal extrapolations of photospheric field measurements. Because such regions usually have dark counterparts in extreme-ultraviolet (EUV) images, we refer to them as coronal holes, even when they appear as narrow lanes or contain sunspots. Revisiting previously identified "AR sources" of slow solar wind from 1998 and 1999, we find that they are all associated with EUV coronal holes; the absence of well-defined He i 1083.0 nm counterparts to some of these holes is attributed to the large flux of photoionizing radiation from neighboring AR loops. Examining a number of AR-associated EUV holes during the 2014 activity maximum, we confirm that they are characterized by wind speeds of ∼300–450 km s⁻¹, O⁷⁺/O⁶⁺ ratios of ∼0.05–0.4, and footpoint field strengths typically of order 30 G. The close spacing between ARs at sunspot maximum limits the widths of unipolar regions and their embedded holes, while the continual emergence of new flux leads to rapid changes in the hole boundaries. Because of the highly nonradial nature of AR fields, the smaller EUV holes are often masked by the overlying canopy of loops, and may be more visible toward one solar limb than at central meridian. As sunspot activity declines, the AR remnants merge to form much larger, weaker, and longer-lived unipolar regions, which harbor the "classical" coronal holes that produce recurrent high-speed streams.

The ages of the components in very short period pre-main-sequence (PMS) binaries are essential to an understanding of their formation. We considered a sample of seven PMS eclipsing binaries (EBs) with ages 1–6.3 MY and component masses 0.2–1.4 ${M}_{\odot }$. The very high precision with which their masses and radii have been measured and the capability provided by the Modules for Experiments in Stellar Astrophysics to calculate their evolutionary tracks at exactly the measured masses allows the determination of age differences of the components independent of their luminosities and effective temperatures. We found that the components of five EBs, ASAS J052821+0338.5, Parenago 1802, JW 380, CoRoT 223992193, and UScoCTIO 5, formed within 0.3 MY of each other. The parameters for the components of V1174 Ori imply an implausible large age difference of 2.7 MY and should be reconsidered. The seventh EB in our sample, RX J0529.4+0041 fell outside the applicability of our analysis.

The structural isomers ethanol (CH₃CH₂OH) and dimethyl ether (CH₃OCH₃) were detected in several low-, intermediate-, and high-mass star-forming regions, including Sgr B2, Orion, and W33A, with the relative abundance ratios of ethanol/dimethyl ether varying from about 0.03 to 3.4. Until now, no experimental data regarding the formation mechanisms and branching ratios of these two species in laboratory simulation experiments could be provided. Here, we exploit tunable photoionization reflectron time-of-flight mass spectrometry (PI-ReTOF-MS) to detect and analyze the production of complex organic molecules (COMs) resulting from the exposure of water/methane (H₂O/CH₄) ices to energetic electrons. The main goal is to understand the formation mechanisms in star-forming regions of two C₂H₆O isomers: ethanol (CH₃CH₂OH) and dimethyl ether (CH₃OCH₃). The results show that the experimental branching ratios favor the synthesis of ethanol versus dimethyl ether (31 ± 11:1). This finding diverges from the abundances observed toward most star-forming regions, suggesting that production routes on interstellar grains to form dimethyl ether might be missing; alternatively, ethanol can be overproduced in the present simulation experiments, such as via radical–radical recombination pathways involving ethyl and hydroxyl radicals. Finally, the PI-ReTOF-MS data suggest the formation of methylacetylene (C₃H₄), ketene (CH₂CO), propene (C₃H₆), vinyl alcohol (CH₂CHOH), acetaldehyde (CH₃CHO), and methyl hydroperoxide (CH₃OOH), in addition to ethane (C₂H₆), methanol (CH₃OH), and CO₂ detected from infrared spectroscopy. The yield of all the confirmed species is also determined.

The infrared dark cloud (IRDC) G028.23-00.19 hosts a massive (1500 M_⊙), cold (12 K), and 3.6–70 μm IR dark clump (MM1) that has the potential to form high-mass stars. We observed this prestellar clump candidate with the Submillimeter Array (∼35 resolution) and Jansky Very Large Array (∼21 resolution) in order to characterize the early stages of high-mass star formation and to constrain theoretical models. Dust emission at 1.3 mm wavelength reveals five cores with masses ≤15 M_⊙. None of the cores currently have the mass reservoir to form a high-mass star in the prestellar phase. If the MM1 clump will ultimately form high-mass stars, its embedded cores must gather a significant amount of additional mass over time. No molecular outflows are detected in the CO (2-1) and SiO (5-4) transitions, suggesting that the SMA cores are starless. By using the NH₃ (1, 1) line, the velocity dispersion of the gas is determined to be transonic or mildly supersonic (ΔV_nt/ΔV_th ∼ 1.1–1.8). The cores are not highly supersonic as some theories of high-mass star formation predict. The embedded cores are four to seven times more massive than the clump thermal Jeans mass and the most massive core (SMA1) is nine times less massive than the clump turbulent Jeans mass. These values indicate that neither thermal pressure nor turbulent pressure dominates the fragmentation of MM1. The low virial parameters of the cores (0.1–0.5) suggest that they are not in virial equilibrium, unless strong magnetic fields of ∼1–2 mG are present. We discuss high-mass star formation scenarios in a context based on IRDC G028.23-00.19, a study case believed to represent the initial fragmentation of molecular clouds that will form high-mass stars.

The accreting millisecond X-ray pulsar SAX J1808.4−3658 shows a peculiar orbital evolution that proceeds at a very fast pace. It is important to identify the underlying mechanism responsible for this behavior because it can help to understand how this system evolves and which physical processes (such as mass loss or spin–orbit coupling) are occurring in the binary. It has also been suggested that, when in quiescence, SAX J1808.4−3658 turns on as a radio pulsar, a circumstance that might provide a link between accreting millisecond pulsars and black-widow (BW) radio pulsars. In this work, we report the results of a deep radio pulsation search at 2 GHz using the Green Bank Telescope in 2014 August and an X-ray study of the 2015 outburst with Chandra, Swift XRT, and INTEGRAL. In quiescence, we detect no radio pulsations and place the strongest limit to date on the pulsed radio flux density of any accreting millisecond pulsar. We also find that the orbit of SAX J1808.4−3658 continues evolving at a fast pace. We compare the orbital evolution of SAX J1808.4−3658 to that of several other accreting and nonaccreting binaries, including BWs, redbacks, cataclysmic variables, black holes, and neutron stars in low-mass X-ray binaries. We discuss two possible scenarios: either the neutron star has a large moment of inertia and is ablating the donor, generating mass loss with an efficiency of 40%, or the donor star has a strong magnetic field of at least 1 kG and is undergoing quasi-cyclic variations due to spin–orbit coupling.

We obtain a new set of analytical solutions for the evolution of a self-gravitating accretion disk by holding the Toomre parameter close to its threshold and obtaining the stress parameter from the cooling rate. In agreement with the previous numerical solutions, furthermore, the accretion rate is assumed to be independent of the disk radius. Extreme situations where the entire disk is either optically thick or optically thin are studied independently, and the obtained solutions can be used for exploring the early or the final phases of a protoplanetary disk evolution. Our solutions exhibit decay of the accretion rate as a power-law function of the age of the system, with exponents −0.75 and −1.04 for optically thick and thin cases, respectively. Our calculations permit us to explore the evolution of the snow line analytically. The location of the snow line in the optically thick regime evolves as a power-law function of time with the exponent −0.16; however, when the disk is optically thin, the location of the snow line as a function of time with the exponent −0.7 has a stronger dependence on time. This means that in an optically thin disk inward migration of the snow line is faster than an optically thick disk.

We present results from daily monitoring of gamma-rays in the energy range from ∼0.5 to ∼100 TeV with the first 17 months of data from the High Altitude Water Cherenkov (HAWC) Observatory. Its wide field of view of 2 steradians and duty cycle of $\gt 95$% are unique features compared to other TeV observatories that allow us to observe every source that transits over HAWC for up to ∼6 hr each sidereal day. This regular sampling yields unprecedented light curves from unbiased measurements that are independent of seasons or weather conditions. For the Crab Nebula as a reference source, we find no variability in the TeV band. Our main focus is the study of the TeV blazars Markarian (Mrk) 421 and Mrk 501. A spectral fit for Mrk 421 yields a power-law index ${\rm{\Gamma }}=2.21\pm {0.14}_{\mathrm{stat}}\pm {0.20}_{\mathrm{sys}}$ and an exponential cut-off ${E}_{0}=5.4\pm {1.1}_{\mathrm{stat}}\pm {1.0}_{\mathrm{sys}}$ TeV. For Mrk 501, we find an index ${\rm{\Gamma }}=1.60\pm {0.30}_{\mathrm{stat}}\pm {0.20}_{\mathrm{sys}}$ and exponential cut-off ${E}_{0}=5.7\pm {1.6}_{\mathrm{stat}}\pm {1.0}_{\mathrm{sys}}$ TeV. The light curves for both sources show clear variability and a Bayesian analysis is applied to identify changes between flux states. The highest per-transit fluxes observed from Mrk 421 exceed the Crab Nebula flux by a factor of approximately five. For Mrk 501, several transits show fluxes in excess of three times the Crab Nebula flux. In a comparison to lower energy gamma-ray and X-ray monitoring data with comparable sampling, we cannot identify clear counterparts for the most significant flaring features observed by HAWC.

Galactic outflows are observed everywhere in star-forming disk galaxies and are critical for galaxy formation. Supernovae (SNe) play the key role in driving the outflows, but there is no consensus as to how much energy, mass, and metal they can launch out of the disk. We perform 3D, high-resolution hydrodynamic simulations to study SNe-driven outflows from stratified media. Assuming the SN rate scales with gas surface density Σ_gas as in the Kennicutt–Schmidt relation, we find that the mass loading factor, η_m, defined as the mass outflow flux divided by the star formation surface density, decreases with increasing Σ_gas as ${\eta }_{{\rm{m}}}\propto {{\rm{\Sigma }}}_{\mathrm{gas}}^{-0.61}$. Approximately Σ_gas ≲ 50 M_⊙ pc⁻² marks when η_m ≳ 1. About 10%–50% of the energy and 40%–80% of the metals produced by SNe end up in the outflows. The tenuous hot phase (T > 3 × 10⁵ K), which fills 60%–80% of the volume at the midplane, carries the majority of the energy and metals in the outflows. We discuss how various physical processes, including the vertical distribution of SNe, photoelectric heating, external gravitational field, and SN rate, affect the loading efficiencies. The relative scale height of gas and SNe is a very important factor in determining the loading efficiencies.

We investigate the star formation properties of a large sample of ∼2300 X-ray-selected Type 2 Active Galactic Nuclei (AGNs) host galaxies out to $z\sim 3$ in the Chandra COSMOS Legacy Survey in order to understand the connection between the star formation and nuclear activity. Making use of the existing multi-wavelength photometric data available in the COSMOS field, we perform a multi-component modeling from far-infrared to near-ultraviolet using a nuclear dust torus model, a stellar population model and a starburst model of the spectral energy distributions (SEDs). Through detailed analyses of SEDs, we derive the stellar masses and the star formation rates (SFRs) of Type 2 AGN host galaxies. The stellar mass of our sample is in the range of $9\lt \mathrm{log}\,{M}_{\mathrm{stellar}}/{M}_{\odot }\lt 12$ with uncertainties of ∼0.19 dex. We find that Type 2 AGN host galaxies have, on average, similar SFRs compared to the normal star-forming galaxies with similar M_stellar and redshift ranges, suggesting no significant evidence for enhancement or quenching of star formation. This could be interpreted in a scenario, where the relative massive galaxies have already experienced substantial growth at higher redshift ($z\gt 3$), and grow slowly through secular fueling processes hosting moderate-luminosity AGNs.

We present measurements of the proper motion of the sub-parsec scale jet at 22 GHz in the nearby FR I galaxy 3C 66B. Observations were made using Very Long Baseline Array (VLBA) at six epochs over four years. A phase-referencing technique was used to improve the image quality of the weak and diffuse jet components. We find that the inner knots are almost stationary, although one of them was expected to be detected with an apparent speed of 0.2 mas yr⁻¹, according to 8 GHz monitoring at the same observation epochs. Clear flux variations are not observed in the core at 22 GHz; in contrast, clear flux enhancement is observed in the core at 8 GHz. We discuss a possible explanation: if the jet has helical structure, the viewing angles of the jet at 8 and 22 GHz differ by a few degrees, if the jet direction is almost along our line of sight. Although these results may imply the existence of a two-zone jet, which has been suggested in certain radio galaxies, it cannot explain the fact that the jet at the higher frequency is slower than that at the lower frequency.

This article reports the results of X-ray studies of the extended TeV γ-ray source VER J2019+368. Suzaku observations conducted to examine properties of the X-ray pulsar wind nebula (PWN) around PSR J2021+3651 revealed that the western region of the X-ray PWN has a source extent of $15^{\prime} \times 10^{\prime}$ with the major axis oriented to that of the TeV emission. The PWN-west spectrum was closely fitted by a power law for absorption at $N({\rm{H}})=({8.2}_{-1.1}^{+1.3})\times {10}^{21}\,{\mathrm{cm}}^{-2}$ and a photon index of ${\rm{\Gamma }}=2.05\pm 0.12$, with no obvious change in the index within the X-ray PWN. The measured X-ray absorption indicates that the distance to the source is much less than the $10\,\mathrm{kpc}$ inferred by radio data. Aside from the PWN, no extended emission was observed around PSR J2021+3651 even by Suzaku. Archival data from the XMM-Newton were also analyzed to complement the Suzaku observations, indicating that the eastern region of the X-ray PWN has a similar spectrum ($N({\rm{H}})=(7.5\pm 0.9)\times {10}^{21}\,{\mathrm{cm}}^{-2}$ and ${\rm{\Gamma }}=2.03\pm 0.10$) and source extent up to at least $12^{\prime}$ along the major axis. The lack of significant change in the photon index and the source extent in X-ray are used to constrain the advection velocity or the diffusion coefficient for accelerated X-ray-producing electrons. A mean magnetic field of $\sim 3\,\mu {\rm{G}}$ is required to account for the measured X-ray spectrum and reported TeV γ-ray spectrum. A model calculation of synchrotron radiation and inverse Compton scattering was able to explain $\sim 80 \%$ of the reported TeV flux, indicating that the X-ray PWN is a major contributor of VER J2019+368.

Pseudo-Newtonian potentials are a tool often used in theoretical astrophysics to capture some key features of a black hole space-time in a Newtonian framework. As a result, one can use Newtonian numerical codes, and Newtonian formalism, in general, in an effective description of important astrophysical processes such as accretion onto black holes. In this paper, we develop a general pseudo-Newtonian formalism, which pertains to the motion of particles, light, and fluids in stationary space-times. In return, we are able to assess the applicability of the pseudo-Newtonian scheme. The simplest and most elegant formulas are obtained in space-times without gravitomagnetic effects, such as the Schwarzschild rather than the Kerr space-time; the quantitative errors are smallest for motion with low binding energy. Included is a ready-to-use set of fluid equations in Schwarzschild space-time in Cartesian and radial coordinates.

In this paper, we study how a flux rope (FR) is formed and evolves into the corresponding structure of a coronal mass ejection (CME) numerically driven by photospheric converging motion. A two-and-a-half-dimensional magnetohydrodynamics simulation is conducted in a chromosphere-transition-corona setup. The initial arcade-like linear force-free magnetic field is driven by an imposed slow motion converging toward the magnetic inversion line at the bottom boundary. The convergence brings opposite-polarity magnetic flux to the polarity inversion, giving rise to the formation of an FR by magnetic reconnection and eventually to the eruption of a CME. During the FR formation, an embedded prominence gets formed by the levitation of chromospheric material. We confirm that the converging flow is a potential mechanism for the formation of FRs and a possible triggering mechanism for CMEs. We investigate the thermal, dynamical, and magnetic properties of the FR and its embedded prominence by tracking their thermal evolution, analyzing their force balance, and measuring their kinematic quantities. The phase transition from the initiation phase to the acceleration phase of the kinematic evolution of the FR was observed in our simulation. The FR undergoes a series of quasi-static equilibrium states in the initiation phase; while in the acceleration phase the FR is driven by Lorentz force and the impulsive acceleration occurs. The underlying physical reason for the phase transition is the change of the reconnection mechanism from the Sweet–Parker to the unsteady bursty regime of reconnection in the evolving current sheet underneath the FR.

Drift effects play a significant role in the transport of charged particles in the heliosphere. A turbulent magnetic field is also known to reduce the effects of particle drifts. The exact nature of this reduction, however, is not clear. This study aims to provide some insight into this reduction and proposes a relatively simple, tractable means of modeling it that provides results in reasonable agreement with numerical simulations of the drift coefficient in a turbulent magnetic field.

We measured phosphorus abundances in 22 FGK dwarfs and giants that span −0.55 < [Fe/H] $\lt$ 0.2 using spectra obtained with the Phoenix high-resolution infrared spectrometer on the Kitt Peak National Observatory Mayall 4 m telescope, the Gemini South Telescope, and the Arcturus spectral atlas. We fit synthetic spectra to the P i feature at 10581 Å to determine abundances for our sample. Our results are consistent with previously measured phosphorus abundances; the average [P/Fe] ratio measured in [Fe/H] bins of 0.2 dex for our stars are within ∼1σ compared to averages from other IR phosphorus studies. Our study provides more evidence that models of chemical evolution using the results of theoretical yields are underproducing phosphorus compared to the observed abundances. Our data better fit a chemical evolution model with phosphorus yields increased by a factor of 2.75 compared to models with unadjusted yields. We also found average [P/Si] = 0.02 ± 0.07 and [P/S] = 0.15 ± 0.15 for our sample, showing no significant deviations from the solar ratios for [P/Si] and [P/S] ratios.

Galactic star formation scaling relations show increased scatter from kpc to sub-kpc scales. Investigating this scatter may hold important clues to how the star formation process evolves in time and space. Here, we combine different molecular gas tracers, different star formation indicators probing distinct populations of massive stars, and knowledge of the evolutionary state of each star-forming region to derive the star formation properties of ∼150 star-forming complexes over the face of the Large Magellanic Cloud (LMC). We find that the rate of massive star formation ramps up when stellar clusters emerge and boost the formation of subsequent generations of massive stars. In addition, we reveal that the star formation efficiency of individual giant molecular clouds (GMCs) declines with increasing cloud gas mass (${M}_{\mathrm{cloud}}$). This trend persists in Galactic star-forming regions and implies higher molecular gas depletion times for larger GMCs. We compare the star formation efficiency per freefall time (${\epsilon }_{\mathrm{ff}}$) with predictions from various widely used analytical star formation models. While these models can produce large dispersions in ${\epsilon }_{\mathrm{ff}}$ similar to those in observations, the origin of the model-predicted scatter is inconsistent with observations. Moreover, all models fail to reproduce the observed decline of ${\epsilon }_{\mathrm{ff}}$ with increasing ${M}_{\mathrm{cloud}}$ in the LMC and the Milky Way. We conclude that analytical star formation models idealizing global turbulence levels and cloud densities and assuming a stationary star formation rate (SFR) are inconsistent with observations from modern data sets tracing massive star formation on individual cloud scales. Instead, we reiterate the importance of local stellar feedback in shaping the properties of GMCs and setting their massive SFR.

Our Cycle 0 ALMA observations confirmed that the Boomerang Nebula is the coldest known object in the universe, with a massive high-speed outflow that has cooled significantly below the cosmic background temperature. Our new CO 1–0 data reveal heretofore unseen distant regions of this ultra-cold outflow, out to ≳120,000 au. We find that in the ultra-cold outflow, the mass-loss rate ($\dot{M}$) increases with radius, similar to its expansion velocity (V)—taking $V\propto r$, we find $\dot{M}\propto {r}^{0.9\mbox{--}2.2}$. The mass in the ultra-cold outflow is $\gtrsim 3.3$M_⊙, and the Boomerang's main-sequence progenitor mass is $\gtrsim 4$M_⊙. Our high angular resolution ($\sim 0\buildrel{\prime\prime}\over{.} 3$) CO J = 3–2 map shows the inner bipolar nebula's precise, highly collimated shape, and a dense central waist of size (FWHM) ∼1740 au × 275 au. The molecular gas and the dust as seen in scattered light via optical Hubble Space Telescope imaging show a detailed correspondence. The waist shows a compact core in thermal dust emission at 0.87–3.3 mm, which harbors $(4\mbox{--}7)\times {10}^{-4}$M_⊙ of very large (∼millimeter-to-centimeter sized), cold ($\sim 20\mbox{--}30$ K) grains. The central waist (assuming its outer regions to be expanding) and fast bipolar outflow have expansion ages of $\lesssim 1925\,\mathrm{years}$ and $\leqslant 1050\,\mathrm{years}$: the "jet-lag" (i.e., torus age minus the fast-outflow age) in the Boomerang supports models in which the primary star interacts directly with a binary companion. We argue that this interaction resulted in a common-envelope configuration, while the Boomerang's primary was an RGB or early-AGB star, with the companion finally merging into the primary's core, and ejecting the primary's envelope that now forms the ultra-cold outflow.

A key goal of the Stage IV dark energy experiments Euclid, LSST, and WFIRST is to measure the growth of structure with cosmic time from weak lensing analysis over large regions of the sky. Weak lensing cosmology will be challenging: in addition to highly accurate galaxy shape measurements, statistically robust and accurate photometric redshift (photo-z) estimates for billions of faint galaxies will be needed in order to reconstruct the three-dimensional matter distribution. Here we present an overview of and initial results from the Complete Calibration of the Color–Redshift Relation (C3R2) survey, which is designed specifically to calibrate the empirical galaxy color–redshift relation to the Euclid depth. These redshifts will also be important for the calibrations of LSST and WFIRST. The C3R2 survey is obtaining multiplexed observations with Keck (DEIMOS, LRIS, and MOSFIRE), the Gran Telescopio Canarias (GTC; OSIRIS), and the Very Large Telescope (VLT; FORS2 and KMOS) of a targeted sample of galaxies that are most important for the redshift calibration. We focus spectroscopic efforts on undersampled regions of galaxy color space identified in previous work in order to minimize the number of spectroscopic redshifts needed to map the color–redshift relation to the required accuracy. We present the C3R2 survey strategy and initial results, including the 1283 high-confidence redshifts obtained in the 2016A semester and released as Data Release 1.

Studying the dynamics of filaments at the pre-eruption phase can shed light on the precursor of eruptive events. Such high-resolution studies (of the order of 01) are highly desirable yet very rare. In this work, we present a detailed observation of a pre-eruption evolution of a filament obtained by the 1.6 m New Solar Telescope (NST) at the Big Bear Solar Observatory (BBSO). One end of the filament is anchored at the sunspot in the NOAA active region (AR) 11515, which is well observed by NST Hα off-bands from four hours before to one hour after the filament eruption. A M1.6 flare is associated with the eruption. We observed persistent downflowing materials along the Hα multi-threaded component of the loop toward the AR end during the pre-eruption phase. We traced the trajectories of plasma blobs along the Hα threads and obtained a plane-of-sky velocity of 45 km s⁻¹ on average. Furthermore, we estimated the real velocities of the downflows and the altitude of the filament by matching the observed Hα threads with magnetic field lines extrapolated from a nonlinear force-free field model. Observations of chromospheric brightenings at the footpoints of the falling plasma blobs are also presented. The lower limit of the kinetic energy per second of the downflows through the brightenings is found to be ∼10²¹ erg. Larger FOV observations from BBSO full-disk Hα images show that the AR end of the filament started ascending four hours before the flare. We attribute the observed downflows at the AR end of the filament to the draining effect of the filament rising prior to its eruption. During the slow-rise phase, the downflows continuously drained away ∼10¹⁵g mass from the filament over a few hours, which is believed to be essential for the instability, and could be an important precursor of eruptive events.

We investigate a sample of 622 blazars with measured fluxes at 12 wavebands across the radio-to-gamma-ray spectrum but without spectroscopic or photometric redshifts. This sample includes hundreds of sources with newly analyzed X-ray spectra reported here. From the synchrotron peak frequencies, estimated by fitting the broadband spectral energy distributions (SEDs), we find that the fraction of high-synchrotron-peaked blazars in these 622 sources is roughly the same as in larger samples of blazars that do have redshifts. We characterize the no-redshift blazars using their infrared colors, which lie in the distinct locus called the WISE blazar strip, then estimate their redshifts using a KNN regression based on the redshifts of the closest blazars in the WISE color–color plot. Finally, using randomly drawn values from plausible redshift distributions, we simulate the SEDs of these blazars and compare them to known blazar SEDs. Based on all these considerations, we conclude that blazars without redshift estimates are unlikely to be high-luminosity, high-synchrotron-peaked objects, which had been suggested in order to explain the "blazar sequence"—an observed trend of SED shape with luminosity—as a selection effect. Instead, the observed properties of no-redshift blazars are compatible with a causal connection between jet power and electron cooling, i.e., a true blazar sequence.

We present an analysis of the internal bulk rotation in the metal-poor globular cluster (GC) NGC 5024 (M53) using radial velocities (RVs) of individual cluster members. We use RV measurements from a previous abundance study of M53 done using the Hydra multi-object spectrograph on the WIYN 3.5 m telescope. The Hydra sample greatly increases the number of RVs available in the central regions of the cluster where the internal rotation is the strongest. The sample of cluster members is further increased through two previous kinematic studies of M53. The combined total sample contains 245 cluster members. With our sample, we are able to create a velocity dispersion profile of the cluster and derive a central velocity dispersion ${\sigma }_{0}=4.0\pm 0.3\ \mathrm{km}\,{{\rm{s}}}^{-1};$ we find that M53 inner regions are characterized by a peak amplitude of rotation equal to $1.4\pm 0.1\ \mathrm{km}\,{{\rm{s}}}^{-1}$ corresponding to a relatively high value of the ratio of the rotation speed to central velocity dispersion (${V}_{\mathrm{rot}}/{\sigma }_{0}=0.35\pm 0.04$). Our data also reveal a radial variation in the orientation of the projected rotation axis suggesting complex internal kinematics.

Present-day semi-empirical models of solar irradiance (SI) variations reconstruct SI changes measured on timescales greater than a day by using spectra computed in one dimensional atmosphere models (1D models), which are representative of various solar surface features. Various recent studies have pointed out, however, that the spectra synthesized in 1D models do not reflect the radiative emission of the inhomogenous atmosphere revealed by high-resolution solar observations. We aimed to derive observation-based atmospheres from such observations and test their accuracy for SI estimates. We analyzed spectropolarimetric data of the Fe i 630 nm line pair in photospheric regions that are representative of the granular quiet-Sun pattern (QS) and of small- and large-scale magnetic features, both bright and dark with respect to the QS. The data were taken on 2011 August 6, with the CRisp Imaging Spectropolarimeter at the Swedish Solar Telescope, under excellent seeing conditions. We derived atmosphere models of the observed regions from data inversion with the SIR code. We studied the sensitivity of results to spatial resolution and temporal evolution, and discuss the obtained atmospheres with respect to several 1D models. The atmospheres derived from our study agree well with most of the 1D models we compare our results with, both qualitatively and quantitatively (within 10%), except for pore regions. Spectral synthesis computations of the atmosphere obtained from the QS observations return an SI between 400 and 2400 nm that agrees, on average, within 2.2% with standard reference measurements, and within −0.14% with the SI computed on the QS atmosphere employed by the most advanced semi-empirical model of SI variations.

The properties of disks around brown dwarfs and very low mass stars (hereafter VLMOs) provide important boundary conditions on the process of planet formation and inform us about the numbers and masses of planets than can form in this regime. We use the Herschel Space Observatory PACS spectrometer to measure the continuum and [O i] 63 μm line emission toward 11 VLMOs with known disks in the Taurus and Chamaeleon I star-forming regions. We fit radiative transfer models to the spectral energy distributions of these sources. Additionally, we carry out a grid of radiative transfer models run in a regime that connects the luminosity of our sources with brighter T Tauri stars. We find that VLMO disks with sizes 1.3–78 au, smaller than typical T Tauri disks, fit well the spectral energy distributions assuming that disk geometry and dust properties are stellar mass independent. Reducing the disk size increases the disk temperature, and we show that VLMOs do not follow previously derived disk temperature–stellar luminosity relationships if the disk outer radius scales with stellar mass. Only 2 out of 11 sources are detected in [O i] despite a better sensitivity than was achieved for T Tauri stars, suggesting that VLMO disks are underluminous. Using thermochemical models, we show that smaller disks can lead to the unexpected [O i] 63 μm nondetections in our sample. The disk outer radius is an important factor in determining the gas and dust observables. Hence, spatially resolved observations with ALMA—to establish if and how disk radii scale with stellar mass—should be pursued further.

We report the appearance of a new radio source at a projected offset of 460 pc from the nucleus of Cygnus A. The flux density of the source (which we designate Cygnus A-2) rose from an upper limit of <0.5 mJy in 1989 to 4 mJy in 2016 (ν = 8.5 GHz), but is currently not varying by more than a few percent per year. The radio luminosity of the source is comparable to the most luminous known supernovae, it is compact in Very Long Baseline Array observations down to a scale of 4 pc, and it is coincident with a near-infrared point source seen in pre-existing adaptive optics and HST observations. The most likely interpretation of this source is that it represents a secondary supermassive black hole in a close orbit around the Cygnus A primary, though an exotic supernova model cannot be ruled out. The gravitational influence of a secondary SMBH at this location may have played an important role in triggering the rapid accretion that has powered the Cygnus A radio jet over the past 10⁷ years.

PHL 6625 is a luminous quasi-stellar object (QSO) at z = 0.3954 located behind the nearby galaxy NGC 247 (z = 0.0005). Hubble Space Telescope observations revealed an arc structure associated with it. We report on spectroscopic observations with the Very Large Telescope and multiwavelength observations from the radio to the X-ray band for the system, suggesting that PHL 6625 and the arc are a close pair of merging galaxies, instead of a strong gravitational lens system. The QSO host galaxy is estimated to be (4–28) × 10¹⁰M_☉ and the mass of the companion galaxy is estimated to be M_* = (6.8 ± 2.4) × 10⁹M_☉, suggesting that this is a minor merger system. The QSO displays typical broad emission lines, from which a black hole mass of about (2–5) × 10⁸M_☉ and an Eddington ratio of about 0.01–0.05 can be inferred. The system represents an interesting and rare case where a QSO is associated with an ongoing minor merger, analogous to Arp 142.

We present the results of very long baseline interferometry (VLBI) observations of gamma-ray bright blazar S5 0716+714 using the Korean VLBI Network (KVN) at the 22, 43, 86, and 129 GHz bands, as part of the Interferometric Monitoring of Gamma-ray Bright active galactic nuclei (iMOGABA) KVN key science program. Observations were conducted in 29 sessions from 2013 January 16 to 2016 March 1, with the source being detected and imaged at all available frequencies. In all epochs, the source was compact on the milliarcsecond scale, yielding a compact VLBI core dominating the synchrotron emission on these scales. Based on the multiwavelength data between 15 GHz (Owens Valley Radio Observatory) and 230 GHz (Submillimeter Array), we found that the source shows multiple prominent enhancements of the flux density at the centimeter (cm) and millimeter (mm) wavelengths, with mm enhancements leading cm enhancements by −16 ± 8 days. The turnover frequency was found to vary between 21 and 69 GHz during our observations. By assuming a synchrotron self-absorption model for the relativistic jet emission in S5 0716+714, we found the magnetic field strength in the mas emission region to be ≤5 mG during the observing period, yielding a weighted mean of 1.0 ± 0.6 mG for higher turnover frequencies (e.g., >45 GHz).

Complex organic molecules (COMs) have been observed toward several low-mass young stellar objects (LYSOs). Small and heterogeneous samples have so far precluded conclusions on typical COM abundances, as well as the origin(s) of abundance variations between sources. We present observations toward 16 deeply embedded (Class 0/I) low-mass protostars using the IRAM 30 m telescope. We detect CH₂CO, CH₃CHO, CH₃OCH₃, CH₃OCHO, CH₃CN, HNCO, and HC₃N toward 67%, 37%, 13%, 13%, 44%, 81%, and 75% of sources, respectively. Median column densities derived using survival analysis range between 6.0 × 10¹⁰ cm⁻² (CH₃CN) and 2.4 × 10¹² cm⁻² (CH₃OCH₃), and median abundances range between 0.48% (CH₃CN) and 16% (HNCO) with respect to CH₃OH. Column densities for each molecule vary by about one order of magnitude across the sample. Abundances with respect to CH₃OH are more narrowly distributed, especially for oxygen-bearing species. We compare observed median abundances with a chemical model for low-mass protostars and find fair agreement, although some modeling work remains to bring abundances higher with respect to CH₃OH. Median abundances with respect to CH₃OH in LYSOs are also found to be generally comparable to observed abundances in hot cores, hot corinos, and massive YSOs. Compared with comets, our sample is comparable for all molecules except HC₃N and CH₂CO, which likely become depleted at later evolutionary stages.

The Outer Scutum–Centaurus (OSC) spiral arm is the most distant molecular spiral arm in the Milky Way, but until recently little was known about this structure. Discovered by Dame and Thaddeus, the OSC lies ∼15 kpc from the Galactic Center. Due to the Galactic warp, it rises to nearly 4° above the Galactic Plane in the first Galactic quadrant, leaving it unsampled by most Galactic plane surveys. Here we observe H ii region candidates spatially coincident with the OSC using the Very Large Array to image radio continuum emission from 65 targets and the Green Bank Telescope to search for ammonia and water maser emission from 75 targets. This sample, drawn from the Wide-field Infrared Survey Explorer Catalog of Galactic H ii Regions, represents every H ii region candidate near the longitude–latitude $({\ell },b)$ locus of the OSC. Coupled with their characteristic mid-infrared morphologies, detection of radio continuum emission strongly suggests that a target is a bona fide H ii region. Detections of associated ammonia or water maser emission allow us to derive a kinematic distance and determine if the velocity of the region is consistent with that of the OSC. Nearly 60% of the observed sources were detected in radio continuum, and more than 20% have ammonia or water maser detections. The velocities of these sources mainly place them beyond the Solar orbit. These very distant high-mass stars have stellar spectral types as early as O4. We associate high-mass star formation at 2 new locations with the OSC, increasing the total number of detected H ii regions in the OSC to 12.

Using archival Rossi X-ray Timing Explorer (RXTE) data, we studied the low-frequency quasi-periodic oscillations (LFQPOs) in the neutron star low-mass X-ray binary (LMXB) Cir X-1 and examined their contribution to frequency–frequency correlations for Z sources. We also studied the orbital phase effects on the LFQPO properties and found them to be phase independent. Comparing LFQPO frequencies in different classes of LMXBs, we found that systems that show both Z and atoll states form a common track with atoll/BH sources in the so-called WK correlation, while persistent Z systems are offset by a factor of about two. We found that neither source luminosity nor mass accretion rate is related to the shift of persistent Z systems. We discuss the possibility of a misidentification of fundamental frequency for horizontal branch oscillations from persistent Z systems and interpreted the oscillations in terms of models based on relativistic precession.

We have examined 40 Nuclear Spectroscopic Telescope Array (NuSTAR) light curves (LCs) of five TeV emitting high synchrotron peaked blazars: 1ES 0229+200, Mrk 421, Mrk 501, 1ES 1959+650, and PKS 2155−304. Four of the blazars showed intraday variability in the NuSTAR energy range of 3–79 keV. Using an autocorrelation function analysis we searched for intraday variability timescales in these LCs and found indications of several between 2.5 and 32.8 ks in eight LCs of Mrk 421, a timescale around 8.0 ks for one LC of Mrk 501, and timescales of 29.6 and 57.4 ks in two LCs of PKS 2155-304. The other two blazars' LCs do not show any evidence for intraday variability timescales shorter than the lengths of those observations; however, the data were both sparser and noisier for them. We found positive correlations with zero lag between soft (3–10 keV) and hard (10–79 keV) bands for most of the LCs, indicating that their emissions originate from the same electron population. We examined spectral variability using a hardness ratio analysis and noticed a general "harder-when-brighter" behavior. The 22 LCs of Mrk 421 observed between 2012 July and 2013 April show that this source was in a quiescent state for an extended period of time and then underwent an unprecedented double-peaked outburst while monitored on a daily basis during 2013 April 10–16. We briefly discuss models capable of explaining these blazar emissions.

We analyze the K2 light curve of the TRAPPIST-1 system. The Fourier analysis of the data suggests P_rot = 3.295 ± 0.003 days. The light curve shows several flares, of which we analyzed 42 events with integrated flare energies of 1.26 × 10³⁰–1.24 × 10³³ erg. Approximately 12% of the flares were complex, multi-peaked eruptions. The flaring and the possible rotational modulation shows no obvious correlation. The flaring activity of TRAPPIST-1 probably continuously alters the atmospheres of the orbiting exoplanets, which makes these less favorable for hosting life.

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

We report on 2.4 yr of radio timing measurements of the magnetar PSR J1622−4950 using the Parkes Observatory, between 2011 November and 2014 March. During this period the torque on the neutron star (inferred from the rotational frequency derivative) varied greatly, though much less erratically than during the 2 yr following its discovery in 2009. During the last year of our measurements the frequency derivative decreased in magnitude monotonically by 20%, to a value of −1.3 × 10⁻¹³ s⁻², a factor of 8 smaller than when it was discovered. The flux density continued to vary greatly during our monitoring through 2014 March, reaching a relatively steady low level after late 2012. The pulse profile varied secularly on a similar timescale as the flux density and torque. A relatively rapid transition in all three properties was evident in early 2013. After PSR J1622−4950 was detected in all of our 87 observations up to 2014 March, we did not detect the magnetar in our resumed monitoring starting in 2015 January and have not detected it in any of the 30 observations conducted through 2016 September.

Core-collapse supernova (SN) explosions expose the structure and environment of massive stars at the moment of their death. We use the global fitting technique of Pejcha & Prieto to estimate a set of physical parameters of 19 normal SNe II, such as their distance moduli, reddenings, ⁵⁶Ni masses ${M}_{\mathrm{Ni}}$, and explosion energies ${E}_{\exp }$ from multicolor light curves and photospheric velocity curves. We confirm and characterize known correlations between ${M}_{\mathrm{Ni}}$ and bolometric luminosity at 50 days after the explosion, and between ${M}_{\mathrm{Ni}}$ and ${E}_{\exp }$. We pay special attention to the observed distribution of ${M}_{\mathrm{Ni}}$ coming from a joint sample of 38 SNe II, which can be described as a skewed-Gaussian-like distribution between $0.005\,{M}_{\odot }$ and $0.280\,{M}_{\odot }$, with a median of $0.031\,{M}_{\odot }$, mean of $0.046\,{M}_{\odot }$, standard deviation of $0.048\,{M}_{\odot }$, and skewness of 3.050. We use a two-sample Kolmogorov–Smirnov test and two-sample Anderson–Darling test to compare the observed distribution of ${M}_{\mathrm{Ni}}$ to results from theoretical hydrodynamical codes of core-collapse explosions with the neutrino mechanism presented in the literature. Our results show that the theoretical distributions obtained from the codes tested in this work, KEPLER and Prometheus Hot Bubble, are compatible with the observations irrespective of different pre-SN calibrations and different maximum mass of the progenitors.

Active galactic nuclei (AGNs) in the high-redshift universe are thought to reside in overdense environments. However, recent works provide controversial results, partly due to the use of different techniques and possible suppression of nearby galaxy formation by AGN feedback. We conducted deep and wide-field imaging observations with the Suprime-Cam on the Subaru Telescope and searched for Lyα emitters (LAEs) around two quasi-stellar objects (QSOs) at z ∼ 4.9 and a radio galaxy at z ∼ 4.5 by using narrowband filters to address these issues more robustly. In the QSO fields, we obtained additional broadband images to select Lyman break galaxies (LBGs) at z ∼ 5 for comparison. We constructed a photometric sample of 301 LAEs and 170 LBGs in total. A wide field of view (34' × 27', corresponding to 80 × 60 comoving Mpc²) of the Suprime-Cam enabled us to probe galaxies in the immediate vicinities of the AGNs and in the blank fields simultaneously and compare various properties of them in a consistent manner. The two QSOs are located near local density peaks (<2σ), and one of the QSOs has a close companion LAE with projected separation of 80 physical kpc. The radio galaxy is found to be near a void of LAEs. The number densities of LAEs/LGBs in a larger spatial scale around the AGNs are not significantly different from those in blank fields. No sign of feedback is found down to ${L}_{\mathrm{Ly}\alpha }\sim {10}^{41.8}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$. Our results suggest that high-redshift AGNs are not associated with extreme galaxy overdensity and that this cannot be attributed to the effect of AGN feedback.

The polarization of the light that is scattered by the coronal ions is influenced by the anisotropic illumination from the photosphere and the magnetic field structuring in the solar corona. The properties of the coronal magnetic fields can be well studied by understanding the polarization properties of coronal forbidden emission lines that arise from magnetic dipole (M1) transitions in the highly ionized atoms that are present in the corona. We present the classical scattering theory of the forbidden lines for a more general case of arbitrary-strength magnetic fields. We derive the scattering matrix for M1 transitions using the classical magnetic dipole model of Casini & Lin and applying the scattering matrix approach of Stenflo. We consider a two-level atom model and neglect collisional effects. The scattering matrix so derived is used to study the Stokes profiles formed in coronal conditions in those regions where the radiative excitations dominate collisional excitations. To this end, we take into account the integration over a cone of an unpolarized radiation from the solar disk incident on the scattering atoms. Furthermore, we also integrate along the line of sight to calculate the emerging polarized line profiles. We consider radial and dipole magnetic field configurations and spherically symmetric density distributions. For our studies we adopt the atomic parameters corresponding to the [Fe xiii] 10747 Å coronal forbidden line. We also discuss the nature of the scattering matrix for M1 transitions and compare it with that for the electric dipole (E1) transitions.

We present ALMA and ATCA observations of the luminous blue variable RMC 127. The radio maps show for the first time the core of the nebula and evidence that the nebula is strongly asymmetric with a Z-pattern shape. Hints of this morphology are also visible in the archival Hubble Space Telescope${\rm{H}}\alpha$ image, which overall resembles the radio emission. The emission mechanism in the outer nebula is optically thin free–free in the radio. At high frequencies, a component of point-source emission appears at the position of the star, up to the ALMA frequencies. The rising flux density distribution (${S}_{\nu }\sim {\nu }^{0.78\pm 0.05}$) of this object suggests thermal emission from the ionized stellar wind and indicates a departure from spherical symmetry with ${n}_{e}(r)\propto {r}^{-2}$. We examine different scenarios to explain this excess of thermal emission from the wind and show that this can arise from a bipolar outflow, supporting the suggestion by other authors that the stellar wind of RMC 127 is aspherical. We fit the data with two collimated ionized wind models, and we find that the mass-loss rate can be a factor of two or more smaller than in the spherical case. We also fit the photometry obtained by IR space telescopes and deduce that the mid- to far-IR emission must arise from extended, cool ($\sim 80\,{\rm{K}}$) dust within the outer ionized nebula. Finally, we discuss two possible scenarios for the nebular morphology: the canonical single-star expanding shell geometry and a precessing jet model assuming the presence of a companion star.

We report the large effort that is producing comprehensive high-level young star cluster (YSC) catalogs for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC 628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC 628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times 10⁵${M}_{\odot }$. After their formation, YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer time frame, confirming their being potentially bound systems. Associations disappear on timescales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC 628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass (≤10⁴${M}_{\odot }$) clusters, suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC 628.

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs, and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of 10⁵–${10}^{7}\,{M}_{\odot }$—masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17\,{M}_{\odot }$ white dwarf by a ${10}^{5}\,{M}_{\odot }\,\mathrm{BH}$ reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions.

We compare several common subgrid implementations of active galactic nucleus (AGN) feedback, focusing on the effects of different triggering mechanisms and the differences between thermal and kinetic feedback. Our main result is that pure thermal feedback that is centrally injected behaves differently from feedback with even a small kinetic component. Specifically, pure thermal feedback results in excessive condensation and smothering of the AGN by cold gas because the feedback energy does not propagate to large enough radii. We do not see large differences between implementations of different triggering mechanisms, as long as the spatial resolution is sufficiently high, probably because all of the implementations tested here trigger strong AGN feedback under similar conditions. In order to assess the role of resolution, we vary the size of the "accretion zone" in which properties are measured to determine the AGN accretion rate and resulting feedback power. We find that a larger accretion zone results in steadier jets but can also allow too much cold gas condensation in simulations with a Bondi-like triggering algorithm. We also vary the opening angle of jet precession and find that a larger precession angle causes more of the jet energy to thermalize closer to the AGN, thereby producing results similar to pure thermal feedback. Our simulations confirm that AGNs can regulate the thermal state of cool-core galaxy clusters and maintain the core in a state that is marginally susceptible to thermal instability followed by precipitation.