Table of contents

We present proof of the axion as a cold dark matter (CDM) candidate to the fully nonlinear order perturbations based on Einstein's gravity. We consider the axion as a coherently oscillating massive classical scalar field without interaction. We present the fully nonlinear and exact, except for ignoring the transverse-tracefree tensor-type perturbation, hydrodynamic equations for an axion fluid in Einstein's gravity. We show that the axion has the characteristic pressure and anisotropic stress; the latter starts to appear from the second-order perturbation. But these terms do not directly affect the hydrodynamic equations in our axion treatment. Instead, what behaves as the effective pressure term in relativistic hydrodynamic equations is the perturbed lapse function and the relativistic result coincides exactly with the one known in the previous non-relativistic studies. The effective pressure term leads to a Jeans scale that is of the solar-system scale for conventional axion mass. As the fully nonlinear and relativistic hydrodynamic equations for an axion fluid coincide exactly with the ones of a zero-pressure fluid in the super-Jeans scale, we have proved the CDM nature of such an axion in that scale.

As the disk formation mechanism(s) in Be stars is(are) as yet unknown, we investigate the role of rapidly rotating radiation-driven winds in this process. We implemented the effects of high stellar rotation on m-CAK models accounting for the shape of the star, the oblate finite disk correction factor, and gravity darkening. For a fast rotating star, we obtain a two-component wind model, i.e., a fast, thin wind in the polar latitudes and an Ω-slow, dense wind in the equatorial regions. We use the equatorial mass densities to explore Hα emission profiles for the following scenarios: (1) a spherically symmetric star, (2) an oblate star with constant temperature, and (3) an oblate star with gravity darkening. One result of this work is that we have developed a novel method for solving the gravity-darkened, oblate m-CAK equation of motion. Furthermore, from our modeling we find that (a) the oblate finite disk correction factor, for the scenario considering the gravity darkening, can vary by at least a factor of two between the equatorial and polar directions, influencing the velocity profile and mass-loss rate accordingly, (b) the Hα profiles predicted by our model are in agreement with those predicted by a standard power-law model for following values of the line-force parameters: $1.5\lesssim k\lesssim 3,\alpha \sim 0.6$, and $\,\delta \gtrsim 0.1$, and (c) the contribution of the fast wind component to the Hα emission line profile is negligible; therefore, the line profiles arise mainly from the equatorial disks of Be stars.

Evidence for dust around supermassive black holes (SMBHs) in the early universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black holes (IMBHs) of $\sim {10}^{4}\mbox{--}{10}^{6}\,{M}_{\odot }$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time-averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from $\sim {10}^{-4}$ to $\sim {10}^{-2}$ in clouds with initial gas densities of ${n}_{{\rm{H}}}=10\mbox{--}1000\,{\mathrm{cm}}^{-3}$. Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between $\lambda \lesssim 20\,\mu {\rm{m}}$ and $\gtrsim 100\,\mu {\rm{m}}$ is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the phase of high accretion, producing higher flux density at $\lesssim 20\,\mu {\rm{m}}$. Therefore, we suggest that a combination of mid-IR observations by the James Webb Space Telescope and far-IR observations by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs. We also extend our simple modeling to SMBHs of ${10}^{8}\mbox{--}{10}^{9}\,{M}_{\odot }$ and show that ALMA can detect SMBHs of $\sim {10}^{9}\,{M}_{\odot }$ at $z\gtrsim 5$.

We study the mean absorption spectrum of the Damped Lyα (DLA) population at z ∼ 2.6 by stacking normalized, rest-frame-shifted spectra of ∼27,000 DLA systems from the DR12 of the Baryon Oscillation Spectroscopic Survey (BOSS)/SDSS-III. We measure the equivalent widths of 50 individual metal absorption lines in five intervals of DLA hydrogen column density, five intervals of DLA redshift, and overall mean equivalent widths for an additional 13 absorption features from groups of strongly blended lines. The mean equivalent width of low-ionization lines increases with N_{H i}, whereas for high-ionization lines the increase is much weaker. The mean metal line equivalent widths decrease by a factor ∼1.1–1.5 from z ∼ 2.1 to z ∼ 3.5, with small or no differences between low- and high-ionization species. We develop a theoretical model, inspired by the presence of multiple absorption components observed in high-resolution spectra, to infer mean metal column densities from the equivalent widths of partially saturated metal lines. We apply this model to 14 low-ionization species and to Al iii, S iii, Si iii, C iv, Si iv, N v, and O vi. We use an approximate derivation for separating the equivalent width contributions of several lines to blended absorption features, and infer mean equivalent widths and column densities from lines of the additional species N i, Zn ii, C ii*, Fe iii, and S iv. Several of these mean column densities of metal lines in DLAs are obtained for the first time; their values generally agree with measurements of individual DLAs from high-resolution, high signal-to-noise ratio spectra when they are available.

We present a systematic survey of multiple velocity-resolved H₂O spectra using Herschel/Heterodyne Instrument for the Far Infrared (HIFI) toward nine nearby actively star-forming galaxies. The ground-state and low-excitation lines (E_up ≤ 130 K) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines (130 K ≤ E_up ≤ 350 K) typically display line shapes similar to CO. We analyze the HIFI observation together with archival SPIRE/PACS H₂O data using a state-of-the-art 3D radiative transfer code that includes the interaction between continuum and line emission. The water excitation models are combined with information on the dust and CO spectral line energy distribution to determine the physical structure of the interstellar medium (ISM). We identify two ISM components that are common to all galaxies: a warm (${T}_{\mathrm{dust}}\sim 40\mbox{--}70$ K), dense ($n({\rm{H}})\sim {10}^{5}\mbox{--}{10}^{6}\,{\mathrm{cm}}^{-3}$) phase that dominates the emission of medium-excitation H₂O lines. This gas phase also dominates the far-IR emission and the CO intensities for ${J}_{\mathrm{up}}\gt 8$. In addition, a cold (${T}_{\mathrm{dust}}\sim 20\mbox{--}30$ K), dense ($n({\rm{H}})\sim {10}^{4}\mbox{--}{10}^{5}\,{\mathrm{cm}}^{-3}$), more extended phase is present. It outputs the emission in the low-excitation H₂O lines and typically also produces the prominent line absorption features. For the two ULIRGs in our sample (Arp 220 and Mrk 231) an even hotter and more compact (R_s ≤ 100 pc) region is present, which is possibly linked to AGN activity. We find that collisions dominate the water excitation in the cold gas and for lines with ${E}_{\mathrm{up}}\leqslant 300$ K and ${E}_{\mathrm{up}}\leqslant 800$ K in the warm and hot component, respectively. Higher-energy levels are mainly excited by IR pumping.

The environment around protoplanetary disks (PPDs) regulates processes that drive the chemical and structural evolution of circumstellar material. We perform a detailed empirical survey of warm molecular hydrogen (H₂) absorption observed against H i-Lyα (Lyα: λ1215.67) emission profiles for 22 PPDs, using archival Hubble Space Telescope ultraviolet (UV) spectra to identify H₂ absorption signatures and quantify the column densities of H₂ ground states in each sightline. We compare thermal equilibrium models of H₂ to the observed H₂ rovibrational level distributions. We find that, for the majority of targets, there is a clear deviation in high-energy states (T_exc ≳ 20,000 K) away from thermal equilibrium populations (T(H₂) ≳ 3500 K). We create a metric to estimate the total column density of non-thermal H₂ (N(H₂)_nLTE) and find that the total column densities of thermal (N(H₂)) and N(H₂)_nLTE correlate for transition disks and targets with detectable C iv-pumped H₂ fluorescence. We compare N(H₂) and N(H₂)_nLTE to circumstellar observables and find that N(H₂)_nLTE correlates with X-ray and far-UV luminosities, but no correlations are observed with the luminosities of discrete emission features (e.g., Lyα, C iv). Additionally, N(H₂) and N(H₂)_nLTE are too low to account for the H₂ fluorescence observed in PPDs, so we speculate that this H₂ may instead be associated with a diffuse, hot, atomic halo surrounding the planet-forming disk. We create a simple photon-pumping model for each target to test this hypothesis and find that Lyα efficiently pumps H₂ levels with T_exc ≥ 10,000 K out of thermal equilibrium.

We investigate the formation of protoplanetary disks around nine solar-mass stars formed in the context of a (40 pc)³ Giant Molecular Cloud model, using ramses adaptive mesh refinement simulations extending over a scale range of about 4 million, from an outer scale of 40 pc down to cell sizes of 2 au. Our most important result is that the accretion process is heterogeneous in multiple ways: in time, in space, and among protostars of otherwise similar mass. Accretion is heterogeneous in time, in the sense that accretion rates vary during the evolution, with generally decreasing profiles, whose slopes vary over a wide range, and where accretion can increase again if a protostar enters a region with increased density and low speed. Accretion is heterogeneous in space, because of the mass distribution, with mass approaching the accreting star–disk system in filaments and sheets. Finally, accretion is heterogeneous among stars, since the detailed conditions and dynamics in the neighborhood of each star can vary widely. We also investigate the sensitivity of disk formation to physical conditions and test their robustness by varying numerical parameters. We find that disk formation is robust even when choosing the least favorable sink particle parameters, and that turbulence cascading from larger scales is a decisive factor in disk formation. We also investigate the transport of angular momentum, finding that the net inward mechanical transport is compensated for mainly by an outward-directed magnetic transport, with a contribution from gravitational torques usually subordinate to the magnetic transport.

During stellar evolution, especially in the pre-main-sequence phase, stellar structure and rotation evolve significantly, causing major changes in the dynamics and global flows of the star. We wish to assess the consequences of these changes on stellar dynamo, internal magnetic field topology, and activity level. To do so, we have performed a series of 3D HD and MHD simulations with the ASH code. We choose five different models characterized by the radius of their radiative zone following an evolutionary track computed by a 1D stellar evolution code. These models characterized stellar evolution from 1 to 50 Myr. By introducing a seed magnetic field in the fully convective model and spreading its evolved state through all four remaining cases, we observe systematic variations in the dynamical properties and magnetic field amplitude and topology of the models. The five MHD simulations develop a strong dynamo field that can reach an equipartition state between the kinetic and magnetic energies and even superequipartition levels in the faster-rotating cases. We find that the magnetic field amplitude increases as it evolves toward the zero-age main sequence. Moreover, the magnetic field topology becomes more complex, with a decreasing axisymmetric component and a nonaxisymmetric one becoming predominant. The dipolar components decrease as the rotation rate and the size of the radiative core increase. The magnetic fields possess a mixed poloidal-toroidal topology with no obvious dominant component. Moreover, the relaxation of the vestige dynamo magnetic field within the radiative core is found to satisfy MHD stability criteria. Hence, it does not experience a global reconfiguration but slowly relaxes by retaining its mixed stable poloidal-toroidal topology.

We study the C8.4-class solar flare SOL2016-05-14T11:34 UT using high-resolution spectral imaging in the Ca ii 8542 Å line obtained with the CRISP imaging spectropolarimeter on the Swedish 1 m Solar Telescope. Spectroscopic inversions of the Ca ii 8542 Å line using the non-LTE code NICOLE are used to investigate the evolution of the temperature and velocity structure in the flaring chromosphere. A comparison of the temperature stratification in flaring and non-flaring areas reveals strong footpoint heating during the flare peak in the lower atmosphere. The temperature of the flaring footpoints between $\mathrm{log}\,{\tau }_{500}\,\approx -2.5\,\mathrm{and}\,-3.5$, where τ₅₀₀ is the continuum optical depth at 500 nm, is $\sim 5\mbox{--}6.5\,\mathrm{kK}$ close to the flare peak, reducing gradually to $\sim 5\,\mathrm{kK}$. The temperature in the middle and upper chromosphere, between $\mathrm{log}\,{\tau }_{500}\approx -3.5$ and −5.5, is estimated to be ∼6.5–20 kK, decreasing to preflare temperatures, ∼5–10 kK, after approximately 15 minutes. However, the temperature stratification of the non-flaring areas is unchanged. The inverted velocity fields show that the flaring chromosphere is dominated by weak downflowing condensations at the formation height of Ca ii 8542 Å.

We present a sample of 1148 ab-type RR Lyrae (RRLab) variables identified from Catalina Surveys Data Release 1, combined with SDSS DR8 and LAMOST DR4 spectral data. We first use a large sample of 860 Galactic halo RRLab stars and derive the circular velocity distributions for the stellar halo. With the precise distances and carefully determined radial velocities (the center-of-mass radial velocities) and by considering the pulsation of the RRLab stars in our sample, we can obtain a reliable and comparable stellar halo circular velocity curve. We follow two different prescriptions for the velocity anisotropy parameter β in the Jeans equation to study the circular velocity curve and mass profile. Additionally, we test two different solar peculiar motions in our calculation. The best result we obtained with the adopted solar peculiar motion 1 of (U, V, W) = (11.1, 12, 7.2) km s⁻¹ is that the enclosed mass of the Milky Way within 50 kpc is (3.75 ± 1.33) × 10¹¹M_⊙ based on β = 0 and the circular velocity 180 ± 31.92 (km s⁻¹) at 50 kpc. This result is consistent with dynamical model results, and it is also comparable to the results of previous similar works.

We present a new method to quantify the value of the escape fraction of ionizing photons, and the existence of ultra-faint galaxies clustered around brighter objects during the epoch of cosmic reionization, using the diffuse Lyα, continuum, and Hα emission observed around galaxies at $z\sim 6$. We model the surface brightness profiles of the diffuse halos, considering the fluorescent emission powered by ionizing photons escaping from the central galaxies, and the nebular emission from satellite star-forming sources, by extending the formalisms developed in Mas-Ribas & Dijkstra and Mas-Ribas et al. The comparison between our predicted profiles and Lyα observations at z = 5.7 and z = 6.6 favors a low ionizing escape fraction, ${f}_{\mathrm{esc}}^{\mathrm{ion}}\sim 5 \%$, for galaxies in the range $-19\gtrsim {M}_{\mathrm{UV}}\gtrsim -21.5$. However, uncertainties and possible systematics in the observations do not allow for firm conclusions. We predict Hα and rest-frame visible continuum observations with the James Webb Space Telescope (JWST), and show that it will be able to detect extended (a few tens of kiloparsecs) fluorescent Hα emission powered by ionizing photons escaping from a bright, $L\gtrsim 5{L}^{* }$, galaxy. Such observations could differentiate fluorescent emission from nebular emission by satellite sources. We discuss how observations and stacking several objects may provide unique constraints on the escape fraction for faint galaxies and/or the abundance of ultra-faint radiation sources.

We present the results from 1.5 and 5 GHz phase-referenced VLBA and 1.5 GHz Karl G. Jansky Very Large Array (VLA) observations of the Seyfert 2 galaxy KISSR 1219, which exhibits double-peaked emission lines in its optical spectrum. The VLA and VLBA data reveal a one-sided core-jet structure at roughly the same position angles, providing evidence of an active galactic nucleus outflow. The absence of dual parsec-scale radio cores puts the binary black-hole picture in doubt for the case of KISSR 1219. The high brightness temperatures of the parsec-scale core and jet components (>10⁶ K) are consistent with this interpretation. Doppler boosting with jet speeds of ≳0.55c to ≳0.25c, going from parsec to kiloparsec scales, at a jet inclination ≳50° can explain the jet one-sidedness in this Seyfert 2 galaxy. A blueshifted broad emission line component in [O iii] is also indicative of an outflow in the emission line gas at a velocity of ∼350 km s⁻¹, while the [O i] doublet lines suggest the presence of shock-heated gas. A detailed line ratio study using the MAPPINGS III code further suggests that a shock+precursor model can explain the line ionization data well. Overall, our data suggest that the radio outflow in KISSR 1219 is pushing the emission line clouds, both ahead of the jet and in a lateral direction, giving rise to the double peak emission line spectra.

Of the 30 or so Galactic magnetars, about 8 are in supernova remnants (SNRs). One of the most extreme magnetars, 1E 1841−045, is at the center of the SNR Kes 73 (G27.4+0.0), whose age is uncertain. We measure its expansion using three Chandra observations over 15 years, obtaining a mean rate of $0.023 \% \pm 0.002 \%$ yr⁻¹. For a distance of 8.5 kpc, we obtain a shell velocity of 1100 km s⁻¹ and infer a blast wave speed of 1400 km s⁻¹. For Sedov expansion into a uniform medium, this gives an age of 1800 years. Derived emission measures imply an ambient density of about 2 cm⁻³ and an upper limit on the swept-up mass of about $70\,{M}_{\odot }$, with lower limits of tens of ${M}_{\odot }$, confirming that Kes 73 is in an advanced evolutionary stage. Our spectral analysis shows no evidence for enhanced abundances as would be expected from a massive progenitor. Our derived total energy is $1.9\times {10}^{51}$ erg, giving a very conservative lower limit to the magnetar's initial period of about 3 ms, unless its energy was lost by non-electromagnetic means. We see no evidence of a wind-blown bubble as would be produced by a massive progenitor, or any evidence that the progenitor of Kes 73/1E 1841−045 was anything but a normal red supergiant producing a Type IIP supernova, though a short-lived stripped-envelope progenitor cannot be absolutely excluded. Kes 73's magnetar thus joins SGR 1900+14 as magnetars resulting from relatively low-mass progenitors.

We constrain the possible presence of a central black hole (BH) in the center of the Large Magellanic Cloud. This requires spectroscopic measurements over an area of the order of a square degree, due to the poorly known position of the kinematic center. Such measurements are now possible with the impressive field of view of the Multi Unit Spectroscopic Explorer (MUSE) on the ESO Very Large Telescope. We used the Calcium Triplet (∼850 nm) spectral lines in many short-exposure MUSE pointings to create a two-dimensional integrated-light line-of-sight velocity map from the $\sim {10}^{8}$ individual spectra, taking care to identify and remove Galactic foreground populations. The data reveal a clear velocity gradient at an unprecedented spatial resolution of 1 arcmin². We fit kinematic models to arrive at a $3\sigma$ upper-mass limit of ${10}^{7.1}\,{M}_{\odot }$ for any central BH—consistent with the known scaling relations for supermassive black holes and their host systems. This adds to the growing body of knowledge on the presence of BHs in low-mass and dwarf galaxies, and their scaling relations with host-galaxy properties, which can shed light on theories of BH growth and host system interaction.

Type Ia supernovae (SNe Ia) are generally agreed to arise from thermonuclear explosions of carbon–oxygen white dwarfs. The actual path to explosion, however, remains elusive, with numerous plausible parent systems and explosion mechanisms suggested. Observationally, SNe Ia have multiple subclasses, distinguished by their light curves and spectra. This raises the question of whether these indicate that multiple mechanisms occur in nature or that explosions have a large but continuous range of physical properties. We revisit the idea that normal and 91bg-like SNe can be understood as part of a spectral sequence in which changes in temperature dominate. Specifically, we find that a single ejecta structure is sufficient to provide reasonable fits of both the normal SN Ia SN 2011fe and the 91bg-like SN 2005bl, provided that the luminosity and thus temperature of the ejecta are adjusted appropriately. This suggests that the outer layers of the ejecta are similar, thus providing some support for a common explosion mechanism. Our spectral sequence also helps to shed light on the conditions under which carbon can be detected in premaximum SN Ia spectra—we find that emission from iron can "fill in" the carbon trough in cool SNe Ia. This may indicate that the outer layers of the ejecta of events in which carbon is detected are relatively metal-poor compared to events in which carbon is not detected.

We present new Submillimeter Array (SMA) observations of CO(2–1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from Herschel observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of ∼3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

Recently, several ultraluminous X-ray (ULX) sources were shown to host a neutron star (NS) accretor. We perform a suite of evolutionary calculations, which show that, in fact, NSs are the dominant type of ULX accretor. Although black holes (BH) dominate early epochs after the star-formation burst, NSs outweigh them after a few 100 Myr and may appear as late as a few gigayears after the end of the star-formation episode. If star formation is a prolonged and continuous event (i.e., not a relatively short burst), NS accretors dominate the ULX population at any time in the solar metallicity environment, whereas BH accretors dominate when the metallicity is sub-solar. Our results show a very clear (and testable) relation between the companion/donor evolutionary stage and the age of the system. A typical NSULX consists of a $\sim 1.3\,{M}_{\odot }$ NS and $\sim 1.0\,{M}_{\odot }$ Red Giant. A typical BH ULX consists of a $\sim 8\,{M}_{\odot }$ BH and $\sim 6\,{M}_{\odot }$ main-sequence star. Additionally, we find that the very luminous ULXs (${L}_{X}\gtrsim {10}^{41}$ erg s⁻¹) are predominantly BH systems ($\sim 9\,{M}_{\odot }$) with Hertzsprung-gap donors ($\sim 2\,{M}_{\odot }$). Nevertheless, some NSULX systems may also reach extremely high X-ray luminosities (≳10⁴¹ erg s⁻¹).

Since the field-line mixing model of Giacalone et al. suggests that ion dropouts cannot happen in the "gradual" solar energetic particle (SEP) event because of the large size of the particle source region in the event, the observational evidence of ion dropouts in the gradual SEP event should challenge the model. We have searched for the presence of ion dropouts in the gradual SEP event during solar cycle 23. From 10 SEP events the synchronized occurrence of ion and electron dropouts is identified in 12 periods. Our main observational facts, including the mean width of electron–ion dropout periods being consistent with the solar wind correlation scale, during the dropout period the dominance of the slab turbulence component and the enhanced turbulence power parallel to the mean magnetic field, and the ion gyroradius dependence of the edge steepness in dropout periods, are all in support of the solar wind turbulence origin of dropout events. Also, our observation indicates that a wide longitude distribution of SEP events could be due to the increase of slab turbulence fraction with the increased longitude distance from the flare-associated active region.

We present multi-wavelength radio observations obtained with the VLA of the protoplanetary disk surrounding the young brown dwarf 2MASS J04442713+2512164 (2M0444) in the Taurus star-forming region. 2M0444 is the brightest known brown dwarf disk at millimeter wavelengths, making this an ideal target to probe radio emission from a young brown dwarf. Thermal emission from dust in the disk is detected at 6.8 and 9.1 mm, whereas the 1.36 cm measured flux is dominated by ionized gas emission. We combine these data with previous observations at shorter sub-mm and mm wavelengths to test the predictions of dust evolution models in gas-rich disks after adapting their parameters to the case of 2M0444. These models show that the radial drift mechanism affecting solids in a gaseous environment has to be either completely made inefficient, or significantly slowed down by very strong gas pressure bumps in order to explain the presence of mm/cm-sized grains in the outer regions of the 2M0444 disk. We also discuss the possible mechanisms for the origin of the ionized gas emission detected at 1.36 cm. The inferred radio luminosity for this emission is in line with the relation between radio and bolometric luminosity valid for for more massive and luminous young stellar objects, and extrapolated down to the very low luminosity of the 2M0444 brown dwarf.

We identify sources with extremely hard X-ray spectra (i.e., with photon indices of ${\rm{\Gamma }}\lesssim 0.6$) in the 13 deg²NuSTAR serendipitous survey, to search for the most highly obscured active galactic nuclei (AGNs) detected at $\gt 10\,\mathrm{keV}$. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower-energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broadband (0.5–24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; ${N}_{{\rm{H}}}\gt 1.5\times {10}^{24}$ cm⁻²) AGNs at low redshift ($z\lt 0.1$) and a likely CT AGN at higher redshift (z = 0.16). Most of the extreme sources would not have been identified as highly obscured based on the low-energy ($\lt 10$ keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray–mid-IR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10–40 keV luminosities of the extreme sources cover a broad range, from $\approx 5\times {10}^{42}$ to 10⁴⁵ erg s⁻¹. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts ($z\lt 0.07$), where we measure a high CT fraction of ${f}_{\mathrm{CT}}^{\mathrm{obs}}={30}_{-12}^{+16} \%$. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50% ± 33%) compared to control samples of "normal" AGNs.

We present a halo-model-based approach to calculate the cross-correlation between $21\,\mathrm{cm}$ H i intensity fluctuations and $\mathrm{Ly}\alpha$ emitters (LAE) during the epoch of reionization (EoR). Ionizing radiation around dark matter halos are modeled as bubbles with the size and growth determined based on the reionization photon production, among other physical parameters. The cross-correlation shows a clear negative-to-positive transition, associated with transition from ionized to neutral hydrogen in the intergalactic medium during EoR. The cross-correlation is subject to several foreground contaminants, including foreground radio point sources important for $21\,\mathrm{cm}$ experiments and low-z interloper emission lines, such as ${\rm{H}}\alpha$, O iii, and O ii, for $\mathrm{Ly}\alpha$ experiments. Our calculations show that by masking out high fluxes in the $\mathrm{Ly}\alpha$ measurement, the correlated foreground contamination on the $21\,\mathrm{cm}$–$\mathrm{Ly}\alpha$ cross-correlation can be dramatically reduced. We forecast the detectability of $21\,\mathrm{cm}$–$\mathrm{Ly}\alpha$ cross-correlation at different redshifts and adopt a Fisher matrix approach to estimate uncertainties on the key EoR parameters that have not been well constrained by other observations of reionization. This halo-model-based approach enables us to explore the EoR parameter space rapidly for different $21\,\mathrm{cm}$ and $\mathrm{Ly}\alpha$ experiments.

Extended main-sequence turn-off (eMSTO) regions are a common feature in color–magnitude diagrams of young- and intermediate-age star clusters in the Magellanic Clouds. The nature of eMSTOs remains debated in the literature. The currently most popular scenarios are extended star formation activity and ranges of stellar rotation rates. Here we study details of differences in main-sequence turn-off (MSTO) morphology expected from spreads in age versus spreads in rotation rates, using Monte Carlo simulations with the Geneva syclist isochrone models that include the effects of stellar rotation. We confirm a recent finding of Niederhofer et al. that a distribution of stellar rotation velocities yields an MSTO extent that is proportional to the cluster age, as observed. However, we find that stellar rotation yields MSTO crosscut widths that are generally smaller than observed ones at a given age. We compare the simulations with high-quality Hubble Space Telescope data of NGC 1987 and NGC 2249, which are the two only relatively massive star clusters with an age of ∼1 Gyr for which such data is available. We find that the distribution of stars across the eMSTOs of these clusters cannot be explained solely by a distribution of stellar rotation velocities, unless the orientations of rapidly rotating stars are heavily biased toward an equator-on configuration. Under the assumption of random viewing angles, stellar rotation can account for ∼60% and ∼40% of the observed FWHM widths of the eMSTOs of NGC 1987 and NGC 2249, respectively. In contrast, a combination of distributions of stellar rotation velocities and stellar ages fits the observed eMSTO morphologies very well.

We obtain high-resolution spectra of nine red giant branch stars in NGC 6681 and perform the first detailed abundance analysis of stars in this cluster. We confirm cluster membership for these stars based on consistent radial velocities of 214.5 ± 3.7 km s⁻¹ and find a mean [Fe/H] = −1.63 ± 0.07 dex and [α/Fe] = 0.42 ± 0.11 dex. Additionally, we confirm the existence of a Na–O anti-correlation in NGC 6681 and identify two populations of stars with unique abundance trends. With the use of HST photometry from Sarajedini et al. and Piotto et al. we are able to identify these two populations as discrete sequences in the cluster CMD. Although we cannot confirm the nature of the polluter stars responsible for the abundance differences in these populations, these results do help put constraints on possible polluter candidates.

Using standard 1D-LTE model atmosphere analysis, we provide an in-depth investigation of iron abundance as derived from neutral and singly ionization iron lines (Fe i, ii) in nearby star clusters. Specifically, we replicate the discrepancy regarding Δ[Fe/H], wherein the difference of Fe ii–Fe i increases for stars of the same cluster with decreasing T_eff, reaching an astonishing 1.0 dex at T_eff ∼ 4000 K. Previous studies have investigated this anomaly in the Pleiades and Hyades clusters with no concrete solution. In this analysis, we probe two samples: 63 wide binary field stars where the primary star is of Sun-like temperatures and the secondary is a K-dwarf, ranging from 4231 K ≤ T_eff ≤ 6453 K, and 33 Hyades stars of temperatures 4268 K ≤ T_eff ≤ 6072 K. Previous studies have found discrepancies on the order of 1.0 dex. However, we find that these studies have neglected line-blending effects of certain Fe ii lines, namely λ = {4508.29 Å, 4993.34 Å, 5197.58 Å, 5325.55 Å, 5425.26 Å, 6456.38 Å}. When these lines are removed from the line-list, we find Δ[Fe/H] decreases to ∼0.6 dex in the field binaries and ∼0.3 dex in the Hyades. The reason for this remaining trend is investigated by probing NLTE effects, as well as age and activity considerations using Ca ii H+K emission and Li absorption, but these results appear to be small to negligible.

In this work we investigate the thermal structure of an off-limb active region (AR) in various non-flaring areas, as it provides key information on the way these structures are heated. In particular, we concentrate on the very hot component ($\gt 3\,\mathrm{MK}$) as it is a crucial element to distinguish between different heating mechanisms. We present an analysis using Fe and Ca emission lines from both the Solar Ultraviolet Measurement of Emitted Radiation (SUMER) on board the Solar and Heliospheric Observatory (SOHO) and the EUV Imaging Spectrometer (EIS) on board Hinode. A data set covering all ionization stages from Fe x to Fe xix has been used for the thermal analysis (both differential emission measure and emission measure, EM). Ca xiv is used for the SUMER-EIS radiometric cross calibration. We show that the very hot plasma is present and persistent almost everywhere in the core of the limb AR. The off-limb AR is clearly structured in Fe xviii. Almost everywhere, the EM analysis reveals plasma at 10 MK (visible in Fe xix emission), which is down to 0.1% of EM of the main $3\,\mathrm{MK}$ plasma. We estimate the power-law index of the hot tail of the EM to be between −8.5 and −4.4. However, the question about the possible existence of a small minor peak at around $10\,\mathrm{MK}$ remains open. The absence in some part of the AR of the Fe xix and Fe xxiii lines (which fall into our spectral range) enables us to determine an upper limit on the EM at these temperatures. Our results include a new Ca xiv 943.59 Å atomic model.

We present a detection of 89 candidates of ultra-diffuse galaxies (UDGs) in a 4.9 degree² field centered on the Hickson Compact Group 95 (HCG 95) using deep g- and r-band images taken with the Chinese Near Object Survey Telescope. This field contains one rich galaxy cluster (Abell 2588 at z = 0.199) and two poor clusters (Pegasus I at z = 0.013 and Pegasus II at z = 0.040). The 89 candidates are likely associated with the two poor clusters, giving about 50–60 true UDGs with a half-light radius ${r}_{{\rm{e}}}\gt 1.5\,\mathrm{kpc}$ and a central surface brightness $\mu (g,0)\gt 24.0$ mag arcsec⁻². Deep $z^{\prime}$-band images are available for 84 of the 89 galaxies from the Dark Energy Camera Legacy Survey (DECaLS), confirming that these galaxies have an extremely low central surface brightness. Moreover, our UDG candidates are spread over a wide range in g − r color, and ∼26% are as blue as normal star-forming galaxies, which is suggestive of young UDGs that are still in formation. Interestingly, we find that one UDG linked with HCG 95 is a gas-rich galaxy with H i mass $1.1\times {10}^{9}$M_⊙ detected by the Very Large Array, and has a stellar mass of ${M}_{\star }\sim 1.8\times {10}^{8}$M_⊙. This indicates that UDGs at least partially overlap with the population of nearly dark galaxies found in deep H i surveys. Our results show that the high abundance of blue UDGs in the HCG 95 field is favored by the environment of poor galaxy clusters residing in H i-rich large-scale structures.

In the context of the GAs Stripping Phenomena in galaxies with Muse (GASP) survey, we present the characterization of JO204, a jellyfish galaxy in A957, a relatively low-mass cluster with $M=4.4\times {10}^{14}\,{M}_{\odot }$. This galaxy shows a tail of ionized gas that extends up to 30 kpc from the main body in the opposite direction of the cluster center. No gas emission is detected in the galaxy outer disk, suggesting that gas-stripping is proceeding outside-in. The stellar component is distributed as a regular disk galaxy; the stellar kinematics shows a symmetric rotation curve with a maximum radial velocity of 200 km s⁻¹ out to 20 kpc from the galaxy center. The radial velocity of the gas component in the central part of the disk follows the distribution of the stellar component; the gas kinematics in the tail retains the rotation of the galaxy disk, indicating that JO204 is moving at high speed in the intracluster medium. Both the emission and radial-velocity maps of the gas and stellar components indicate ram-pressure as the most likely primary mechanism for gas-stripping, as expected given that JO204 is close to the cluster center and it is likely at the first infall in the cluster. The spatially resolved star formation history of JO204 provides evidence that the onset of ram-pressure-stripping occurred in the last 500 Myr, quenching the star formation activity in the outer disk, where the gas has been already completely stripped. Our conclusions are supported by a set of hydrodynamic simulations.

Planetary obliquity determines the meridional distribution of the annual mean insolation. For obliquity exceeding 55°, the weakest insolation occurs at the equator. Stable partial snow and ice cover on such a planet would be in the form of a belt about the equator rather than polar caps. An analytical model of planetary climate is used to investigate the stability of ice caps and ice belts over the widest possible range of parameters. The model is a non-dimensional diffusive Energy Balance Model, representing insolation, heat transport, and ice−albedo feedback on a spherical planet. A complete analytical solution for any obliquity is given and validated against numerical solutions of a seasonal model in the "deep-water" regime of weak seasonal ice line migration. Multiple equilibria and unstable transitions between climate states (ice-free, Snowball, or ice cap/belt) are found over wide swaths of parameter space, including a "Large Ice-Belt Instability" and "Small Ice-Belt Instability" at high obliquity. The Snowball catastrophe is avoided at weak radiative forcing in two different scenarios: weak albedo feedback and inefficient heat transport (favoring stable partial ice cover), or efficient transport at high obliquity (favoring ice-free conditions). From speculative assumptions about distributions of planetary parameters, three-fourths to four-fifths of all planets with stable partial ice cover should be in the form of Earth-like polar caps.

The goal of the Event Horizon Telescope (EHT) is to provide spatially resolved images of Sgr A*, the source associated with the Galactic Center black hole. Because Sgr A* varies on timescales that are short compared to an EHT observing campaign, it is interesting to ask whether variability contains information about the structure and dynamics of the accretion flow. In this paper, we introduce "time-domain filtering," a technique to filter time fluctuating images with specific temporal frequency ranges and to demonstrate the power and usage of the technique by applying it to mock millimeter wavelength images of Sgr A*. The mock image data is generated from the General Relativistic Magnetohydrodynamic (GRMHD) simulation and the general relativistic ray-tracing method. We show that the variability on each line of sight is tightly correlated with a typical radius of emission. This is because disk emissivity fluctuates on a timescale of the order of the local orbital period. Time-domain filtered images therefore reflect the model dependent emission radius distribution, which is not accessible in time-averaged images. We show that, in principle, filtered data have the power to distinguish between models with different black-hole spins, different disk viewing angles, and different disk orientations in the sky.

Recent observations have detected galaxies at high-redshift $z\sim 6\mbox{--}11$, and revealed the diversity of their physical properties, from normal star-forming galaxies to starburst galaxies. To understand the properties of these observed galaxies, it is crucial to understand the star formation (SF) history of high-redshift galaxies under the influence of stellar feedback. In this work, we present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution toward observable galaxies at $z\sim 6$. We focus on three different galaxies that end up in halos with masses ${M}_{{\rm{h}}}=2.4\times {10}^{10}\,{h}^{-1}\ {M}_{\odot }$ (Halo-10), $1.6\times {10}^{11}\,{h}^{-1}\ {M}_{\odot }$ (Halo-11), and $0.7\times {10}^{12}\,{h}^{-1}\,{M}_{\odot }$ (Halo-12) at z = 6. Our simulations also probe the impacts of different subgrid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that SF occurs intermittently due to supernova (SN) feedback at $z\gtrsim 10$, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no feedback or lower SF efficiency models can sustain a galactic disk for long periods $\gtrsim 10\,\mathrm{Myr}$. The expulsion of gas at the galactic center also affects the inner dark matter density profile for a short period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the SF rates and stellar masses of observed Lyα emitters at $z\sim 7\mbox{--}8$ fairly well given the observational uncertainties.

Hot-Jupiters are subject to extreme radiation and plasma flows coming from their host stars. Past ultraviolet Hubble Space Telescope observations, supported by hydrodynamic models, confirmed that these factors lead to the formation of an extended envelope, part of which lies beyond the Roche lobe. We use gas-dynamic simulations to study the impact of time variations in the parameters of the stellar wind, namely that of coronal mass ejections (CMEs), on the envelope of the typical hot-Jupiter HD 209458b. We consider three CMEs characterized by different velocities and densities, taking their parameters from typical CMEs observed for the Sun. The perturbations in the ram-pressure of the stellar wind during the passage of each CME tear off most of the envelope that is located beyond the Roche lobe. This leads to a substantial increase of the mass-loss rates during the interaction with the CME. We find that the mass lost by the planet during the whole crossing of a CME is of ≈10¹⁵ g, regardless of the CME taken into consideration. We also find that over the course of 1 Gyr, the mass lost by the planet because of CME impacts is comparable to that lost because of high-energy stellar irradiation.

We present an analysis of ${[{\rm{O}}{\rm{I}}]}_{63}$, [O iii]₈₈, [N ii]₁₂₂, and ${[{\rm{C}}{\rm{II}}]}_{158}$ far-infrared (FIR) fine-structure line observations obtained with Herschel/PACS, for ∼240 local luminous infrared galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey. We find pronounced declines ("deficits") of line-to-FIR continuum emission for [N ii]₁₂₂, ${[{\rm{O}}{\rm{I}}]}_{63}$, and ${[{\rm{C}}{\rm{II}}]}_{158}$ as a function of FIR color and infrared luminosity surface density, ${{\rm{\Sigma }}}_{\mathrm{IR}}$. The median electron density of the ionized gas in LIRGs, based on the [N ii]₁₂₂/[N ii]₂₀₅ ratio, is ${n}_{{\rm{e}}}$ = 41 cm⁻³. We find that the dispersion in the ${[{\rm{C}}{\rm{II}}]}_{158}$ deficit of LIRGs is attributed to a varying fractional contribution of photodissociation regions (PDRs) to the observed ${[{\rm{C}}{\rm{II}}]}_{158}$ emission, f($[{\rm{C}}\,{\rm{II}}{]}_{158}^{\mathrm{PDR}}$) = $[{\rm{C}}\,{\rm{II}}{]}_{158}^{\mathrm{PDR}}$/${[{\rm{C}}{\rm{II}}]}_{158}$, which increases from ∼60% to ∼95% in the warmest LIRGs. The ${[{\rm{O}}{\rm{I}}]}_{63}$/$[{\rm{C}}\,{\rm{II}}{]}_{158}^{\mathrm{PDR}}$ ratio is tightly correlated with the PDR gas kinetic temperature in sources where ${[{\rm{O}}{\rm{I}}]}_{63}$ is not optically thick or self-absorbed. For each galaxy, we derive the average PDR hydrogen density, ${n}_{{\rm{H}}}$, and intensity of the interstellar radiation field, G, in units of ${G}_{0}$ and find G/${n}_{{\rm{H}}}$ ratios of ∼0.1–50 ${G}_{0}$ cm³, with ULIRGs populating the upper end of the distribution. There is a relation between G/${n}_{{\rm{H}}}$ and ${{\rm{\Sigma }}}_{\mathrm{IR}}$, showing a critical break at ${{\rm{\Sigma }}}_{\mathrm{IR}}^{* }$ ≃ 5 × 10¹⁰L_⊙ kpc⁻². Below ${{\rm{\Sigma }}}_{\mathrm{IR}}^{* }$, G/${n}_{{\rm{H}}}$ remains constant, ≃0.32 ${G}_{0}$ cm³, and variations in ${{\rm{\Sigma }}}_{\mathrm{IR}}$ are driven by the number density of star-forming regions within a galaxy, with no change in their PDR properties. Above ${{\rm{\Sigma }}}_{\mathrm{IR}}^{* }$, G/${n}_{{\rm{H}}}$ increases rapidly with ${{\rm{\Sigma }}}_{\mathrm{IR}}$, signaling a departure from the typical PDR conditions found in normal star-forming galaxies toward more intense/harder radiation fields and compact geometries typical of starbursting sources.

We compare analytic predictions of supernova light curves with recent high-quality data from SN2011fe (Ia), KSN2011b (Ia), and the Palomar Transient Factory and the La Silla-QUEST variability survey (LSQ) (Ia). Because of the steady, fast cadence of observations, KSN2011b provides unique new information on SNe Ia: the smoothness of the light curve, which is consistent with significant large-scale mixing during the explosion, possibly due to 3D effects (e.g., Rayleigh–Taylor instabilities), and provides support for a slowly varying leakage (mean opacity). For a more complex light curve (SN2008D, SN Ib), we separate the luminosity due to multiple causes and indicate the possibility of a radioactive plume. The early rise in luminosity is shown to be affected by the opacity (leakage rate) for thermal and non-thermal radiation. A general derivation of Arnett's rule again shows that it depends upon all processes heating the plasma, not just radioactive ones, so that SNe Ia will differ from SNe Ibc if the latter have multiple heating processes.

We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis to the first 7.4 years of Fermi observations, and in two separate energy bands 0.1–0.8 GeV and 0.8–300 GeV, a total of 4547 flares were detected with significance greater than $6\sigma$ (before trials), on the timescale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources were identified. Based on positional coincidence, likely counterparts have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of freshly accelerated electrons is never harder than $p\sim 2$.

We present molecular gas-mass estimates for a sample of 13 local galaxies whose kinematic and star-forming properties closely resemble those observed in z ≈ 1.5 main-sequence galaxies. Plateau de Bure observations of the CO[1-0] emission line and Herschel Space Observatory observations of the dust emission both suggest molecular gas-mass fractions of ∼20%. Moreover, dust emission modeling finds T_dust < 30 K, suggesting a cold dust distribution compared to their high infrared luminosity. The gas-mass estimates argue that z ∼ 0.1 DYNAMO galaxies not only share similar kinematic properties with high-z disks, but they are also similarly rich in molecular material. Pairing the gas-mass fractions with existing kinematics reveals a linear relationship between f_gas and σ/v_c, consistent with predictions from stability theory of a self-gravitating disk. It thus follows that high gas-velocity dispersions are a natural consequence of large gas fractions. We also find that the systems with the lowest t_dep (∼0.5 Gyr) have the highest ratios of σ/v_c and more pronounced clumps, even at the same high molecular gas fraction.

We perform the first systematic study of how dynamical stellar tides and general relativistic (GR) effects affect the dynamics and outcomes of binary-single interactions. For this, we have constructed an N-body code that includes tides in the affine approximation, where stars are modeled as self-similar ellipsoidal polytropes, and GR corrections using the commonly used post-Newtonian formalism. Using this numerical formalism, we are able resolve the leading effect from tides and GR across several orders of magnitude in both stellar radius and initial target binary separation. We find that the main effect from tides is the formation of two-body tidal captures that form during the chaotic and resonant evolution of the triple system. The two stars undergoing the capture spiral in and merge. The inclusion of tides can thus lead to an increase in the stellar coalescence rate. We also develop an analytical framework for calculating the cross section of tidal inspirals between any pair of objects with similar mass. From our analytical and numerical estimates, we find that the rate of tidal inspirals relative to collisions increases as the initial semimajor axis of the target binary increases and the radius of the interacting tidal objects decreases. The largest effect is therefore found for triple systems hosting white dwarfs and neutron stars (NSs). In this case, we find the rate of highly eccentric white dwarf—NS mergers to likely be dominated by tidal inspirals. While tidal inspirals occur rarely, we note that they can give rise to a plethora of thermonuclear transients, such as Ca-rich transients.

Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the properties of the SN outburst for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code Prometheus, using a presupernova model of a 15 M_⊙ star as initial data. On approaching homologous expansion, the hydrodynamic and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian, radiation hydrodynamics) code Crab to follow the evolution of the blast wave during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about $0.5\times {10}^{51}$ erg produces ⁵⁶Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of ⁵⁶Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic decline in luminosity from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface.

We present a new map of interstellar reddening, covering the 39% of the sky with low H i column densities (${N}_{{\rm{H}}{\rm{I}}}\lt 4\times {10}^{20}\,$ cm⁻² or $E(B-V)\approx 45$ mmag) at $16\buildrel{\,\prime}\over{.} 1$ resolution, based on all-sky observations of Galactic H i emission by the HI4PI Survey. In this low-column-density regime, we derive a characteristic value of ${N}_{{\rm{H}}{\rm{I}}}/E(B-V)=8.8\,\times \,{10}^{21}\,{\mathrm{cm}}^{2}\,{\mathrm{mag}}^{-1}$ for gas with $| {v}_{\mathrm{LSR}}| \lt 90$ km s⁻¹ and find no significant reddening associated with gas at higher velocities. We compare our H i-based reddening map with the Schlegel et al. (SFD) reddening map and find them consistent to within a scatter of $\simeq 5$ mmag. Further, the differences between our map and the SFD map are in excellent agreement with the low-resolution ($4\buildrel{\circ}\over{.} 5$) corrections to the SFD map derived by Peek and Graves based on observed reddening toward passive galaxies. We therefore argue that our H i-based map provides the most accurate interstellar reddening estimates in the low-column-density regime to date. Our reddening map is made publicly available at doi.org/10.7910/DVN/AFJNWJ.

The spatial distribution of oxygen in the interstellar medium of galaxies is the key to understanding how efficiently metals that are synthesized in massive stars can be redistributed across a galaxy. We present here a case study in the nearby spiral galaxy NGC 1365 using 3D optical data obtained in the TYPHOON Program. We find systematic azimuthal variations of the H ii region oxygen abundance imprinted on a negative radial gradient. The 0.2 dex azimuthal variations occur over a wide radial range of 0.3–0.7 R₂₅ and peak at the two spiral arms in NGC 1365. We show that the azimuthal variations can be explained by two physical processes: gas undergoes localized, sub-kiloparsec-scale self-enrichment when orbiting in the inter-arm region, and experiences efficient, kiloparsec-scale mixing-induced dilution when spiral density waves pass through. We construct a simple chemical evolution model to quantitatively test this picture and find that our toy model can reproduce the observations. This result suggests that the observed abundance variations in NGC 1365 are a snapshot of the dynamical local enrichment of oxygen modulated by spiral-driven, periodic mixing and dilution.

The C i 135.58 nm line is located in the wavelength range of NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission. We study the formation and diagnostic potential of this line by means of non local-thermodynamic-equilibrium modeling, employing both 1D and 3D radiation-magnetohydrodynamic models. The C i/C ii ionization balance is strongly influenced by photoionization by Lyα emission. The emission in the C i 135.58 nm line is dominated by a recombination cascade and the line forming region is optically thick. The Doppler shift of the line correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.5 Mm height. With IRIS, the C i 135.58 nm line is usually observed together with the O i 135.56 nm line, and from the Doppler shift of both lines, we obtain the velocity difference between the line forming regions of the two lines. From the ratio of the C i/O i line core intensity, we can determine the distance between the C i and the O i forming layers. Combined with the velocity difference, the velocity gradient at mid-chromospheric heights can be derived. The C i/O i total intensity line ratio is correlated with the inverse of the electron density in the mid-chromosphere. We conclude that the C i 135.58 nm line is an excellent probe of the middle chromosphere by itself, and together with the O i 135.56 nm line the two lines provide even more information, which complements other powerful chromospheric diagnostics of IRIS such as the Mg ii h and k lines and the C ii lines around 133.5 nm.

Asteroseismology is a useful tool that is usually used to probe stellar interiors and to determine stellar fundamental parameters, such as stellar mass, radius, and surface gravity. In order to probe stellar interiors, making comparisons between observations and models is usually used with the ${\chi }^{2}$-minimization method. The work of Wu & Li reported that the best parameter determined by the ${\chi }^{2}$-matching process is the acoustic radius for pure p-mode oscillations. In the present work, based on the theoretical calculations of Wu & Li, we will independently analyze the seismic observations of KIC 6225718 to determine its fundamental parameters and to investigate its interior properties. First, in order to test the method, we use it in the Sun to determine its fundamental parameters and to investigate interiors. Second, we independently determine the fundamental parameters of KIC 6225718 without any other non-seismic constraint. Therefore, those determined fundamental parameters are independent of those determined by other methods. They can be regarded as independent references in other analyses. Finally, we analyze the stellar internal structure and find that KIC 6225718 has a convective core with the size of 0.078–0.092 ${R}_{\odot }$. Its overshooting parameter ${f}_{\mathrm{ov}}$ in the core is around 0.010. In addition, its center hydrogen ${X}_{{\rm{c}}}$ is about 0.264–0.355.

We compare the optical properties of the host galaxies of radio-quiet (RQ) and radio-loud (RL) Type 2 active galactic nuclei (AGNs) to infer whether the jet production efficiency depends on the host properties or is determined just by intrinsic properties of the accretion flows. We carefully select galaxies from SDSS, FIRST, and NVSS catalogs. We confirm previous findings that the fraction of RL AGNs depends on the black-hole (BH) masses, and on the Eddington ratio. The comparison of the nature of the hosts of RL and RQ AGNs, therefore, requires pair-matching techniques. Matching in BH mass and Eddington ratio allows us to study the differences between galaxies hosting RL and RQ AGNs that have the same basic accretion parameters. We show that these two samples differ predominantly in the host-galaxy concentration index, morphological type (in the RL sample the frequency of elliptical galaxies becoming larger with increasing radio loudness), and nebular extinction (galaxies with highest radio loudness showing only low nebular extinction). Contrary to some previous studies, we find no significant difference between our radio-loud and radio-quiet samples regarding merger/interaction features.

On the surface of icy dust grains in the dense regions of the interstellar medium, a rich chemistry can take place. Due to the low temperature, reactions that proceed via a barrier can only take place through tunneling. The reaction ${\rm{H}}+{{\rm{H}}}_{2}{{\rm{O}}}_{2}\longrightarrow {{\rm{H}}}_{2}{\rm{O}}+\mathrm{OH}$ is such a case with a gas-phase barrier of ∼26.5 kJ mol⁻¹. Still, the reaction is known to be involved in water formation on interstellar grains. Here, we investigate the influence of a water ice surface and of bulk ice on the reaction rate constant. Rate constants are calculated using instanton theory down to 74 K. The ice is taken into account via multiscale modeling, describing the reactants and the direct surrounding at the quantum mechanical level with density functional theory (DFT), while the rest of the ice is modeled on the molecular mechanical level with a force field. We find that H₂O₂ binding energies cannot be captured by a single value, but rather they depend on the number of hydrogen bonds with surface molecules. In highly amorphous surroundings, the binding site can block the routes of attack and impede the reaction. Furthermore, the activation energies do not correlate with the binding energies of the same sites. The unimolecular rate constants related to the Langmuir–Hinshelwood mechanism increase as the activation energy decreases. Thus, we provide a lower limit for the rate constant and argue that rate constants can have values up to two order of magnitude larger than this limit.

The localization of the repeating fast radio burst (FRB), FRB 121102, suggests that it is associated with a persistent radio-luminous compact source in the FRB host galaxy. Using the FIRST radio catalog, I present a search for luminous persistent sources in nearby galaxies, with radio luminosities $\gt 10 \%$ of the FRB 121102 persistent source luminosity. The galaxy sample contains about 30% of the total galaxy g-band luminosity within $\lt 108$ Mpc, in a footprint of 10,600 deg². After rejecting sources likely due to active galactic nuclei activity or background sources, I am left with 11 candidates that are presumably associated with galactic disks or star-formation regions. At least some of these candidates are likely to be due to chance alignment. In addition, I find 85 sources within $1^{\prime\prime}$ of galactic nuclei. Assuming that the radio persistent sources are not related to galactic nuclei and that they follow the galaxy g-band light, the 11 sources imply a 95% confidence upper limit on the space density of luminous persistent sources of $\lesssim 5\times {10}^{-5}$ Mpc⁻³, and that at any given time only a small fraction of galaxies host a radio-luminous persistent source ($\lesssim {10}^{-3}$${L}_{* }^{-1}$). Assuming a persistent source lifetime of 100 years, this implies a birth rate of $\lesssim 5\times {10}^{-7}$ yr⁻¹ Mpc⁻³. Given the FRB volumetric rate, and assuming that all FRBs repeat and are associated with persistent radio sources, this sets a lower limit on the rate of FRB events per persistent source of $\gtrsim 0.8$ yr⁻¹. I argue that these 11 candidates are good targets for FRB searches and I estimate the FRB event rate from these candidates.

We present a demonstration of delensing the observed cosmic microwave background (CMB) B-mode polarization anisotropy. This process of reducing the gravitational-lensing-generated B-mode component will become increasingly important for improving searches for the B modes produced by primordial gravitational waves. In this work, we delens B-mode maps constructed from multi-frequency SPTpol observations of a 90 deg² patch of sky by subtracting a B-mode template constructed from two inputs: SPTpol E-mode maps and a lensing potential map estimated from the Herschel 500 μm map of the cosmic infrared background. We find that our delensing procedure reduces the measured B-mode power spectrum by $28$% in the multipole range $300\lt {\ell }\lt 2300;$ this is shown to be consistent with expectations from simulations and to be robust against systematics. The null hypothesis of no delensing is rejected at $6.9\sigma$. Furthermore, we build and use a suite of realistic simulations to study the general properties of the delensing process and find that the delensing efficiency achieved in this work is limited primarily by the noise in the lensing potential map. We demonstrate the importance of including realistic experimental non-idealities in the delensing forecasts used to inform instrument and survey-strategy planning of upcoming lower-noise experiments, such as CMB-S4.

We observed the $J=5-4$ rotational lines of the normal species and three ¹³C isotopologues of HC₃N at the 45 GHz band toward two low-mass starless cores, L1521B and L134N (L183), using the Nobeyama 45 m radio telescope in order to study the main formation pathways of HC₃N in each core. The abundance ratios of the three ¹³C isotopologues in L1521B are derived to be [H¹³CCCN]:[HC¹³CCN]:[HCC¹³CN] = 0.98 (±0.14):1.00:1.52 (±0.16) ($1\sigma$). The fractionation pattern is consistent with that at the cyanopolyyne peak in Taurus Molecular Cloud-1. This fractionation pattern suggests that the main formation pathway of HC₃N is the neutral–neutral reaction between C₂H₂ and CN. On the other hand, their abundance ratios in L134N are found to be [H¹³CCCN]:[HC¹³CCN]:[HCC¹³CN] = 1.5 (±0.2):1.0:2.1 (±0.4) ($1\sigma$), which are different from those in L1521B. From this fractionation pattern, we propose that the reaction between HNC and CCH is a possible main formation pathway of HC₃N in L134N. We find out that the main formation pathways of the same molecule are not common even in similar physical conditions. We discuss the possible factors that could make a difference in the fractionation pattern between L134N and L1521B/TMC-1.

In the presence of a magnetic field and weakly ionizing winds, ohmic dissipation is expected to take place in the envelopes of Jovian and lower-mass planets alike. While the process has been investigated on the former, there have been no studies done on mini-Neptunes so far. From structure and thermal evolution models, we determine that the required energy deposition for halting the contraction of mini-Neptunes increases with planetary mass and envelope fraction. Scaled to the insolation power, the ohmic heating needed is small: $\sim {10}^{-5}$ orders of magnitude lower than for exo-Jupiters $\sim {10}^{-2}$. Conversely, from solving the magnetic induction equation, we find that ohmic energy is dissipated more readily for lower-mass planets and those with larger envelope fractions. Combining these two trends, we find that ohmic dissipation in hot mini-Neptunes is strong enough to inflate their radii ($\sim {10}^{15}$ W for ${T}_{\mathrm{eq}}=1400\,{\rm{K}}$). The implication is that the radii of hot mini-Neptunes may be attributed in part to ohmic heating. Thus, there is a trade-off between ohmic dissipation and H/He content for hot mini-Neptunes, adding a new degeneracy for the interpretation of the composition of such planets. In addition, ohmic dissipation would make mini-Neptunes more vulnerable to atmospheric evaporation.

Dust drifting inward in protoplanetary disks is subject to increasing temperatures. In laboratory experiments, we tempered basaltic dust between 873 K and 1273 K and find that the dust grains change in size and composition. These modifications influence the outcome of self-consistent low speed aggregation experiments showing a transition temperature of 1000 K. Dust tempered at lower temperatures grows to a maximum aggregate size of 2.02 ± 0.06 mm, which is 1.49 ± 0.08 times the value for dust tempered at higher temperatures. A similar size ratio of 1.75 ± 0.16 results for a different set of collision velocities. This transition temperature is in agreement with orbit temperatures deduced for observed extrasolar planets. Most terrestrial planets are observed at positions equivalent to less than 1000 K. Dust aggregation on the millimeter-scale at elevated temperatures might therefore be a key factor for terrestrial planet formation.

Solar wind fluctuations reveal the ubiquity of intermittency, which is believed to affect the spectral signatures of turbulence. In this work, based on simulation of driven compressible MHD turbulence, we apply the wavelet technique to the magnetic field and velocity to identify intermittency, and we analyze the influence of the intermittency on the quasi-perpendicular scaling in the inertial range. The numerical results show that the original magnetic and velocity fluctuations are anisotropic, and have a power anisotropy with a spectral index approaching the Iroshnikov–Kraichnan $-3/2$ scaling in the direction quasi-perpendicular to the local mean magnetic field. As in observations of the solar wind fluctuations, as the scale decreases in the simulation, the calculated probability distribution functions (pdfs) of the wavelet coefficients become extended on both tails of the non-Gaussian distribution, with a rapid increase in flatness. After intermittency has been removed from the driven turbulence, at each scale, the pdfs approach a Gaussian distribution, with the flatness being ∼3. Meanwhile, the quasi-perpendicular scaling for both fluctuations becomes steeper and close to a Kolmogorov $-5/3$ scaling, which may be a result of the stronger intermittency in the quasi-perpendicular direction and at the smaller scales. These results suggest that there is intermittency superposed on the "background" turbulence that seems to have the Kolmogorov scaling, whereby the overall slope is getting flatter with the involvement of intermittency.

The progenitor systems of the class of "Ca-rich transients" is a key open issue in time domain astrophysics. These intriguing objects exhibit unusually strong calcium line emissions months after explosion, fall within an intermediate luminosity range, are often found at large projected distances from their host galaxies, and may play a vital role in enriching galaxies and the intergalactic medium. Here we present multiwavelength observations of iPTF15eqv in NGC 3430, which exhibits a unique combination of properties that bridge those observed in Ca-rich transients and SNe Ib/c. iPTF15eqv has among the highest [Ca ii]/[O i] emission line ratios observed to date, yet is more luminous and decays more slowly than other Ca-rich transients. Optical and near-infrared photometry and spectroscopy reveal signatures consistent with the supernova explosion of a $\lesssim 10\,{M}_{\odot }$ star that was stripped of its H-rich envelope via binary interaction. Distinct chemical abundances and ejecta kinematics suggest that the core collapse occurred through electron-capture processes. Deep limits on possible radio emission made with the Jansky Very Large Array imply a clean environment (n ≲ 0.1 cm⁻³) within a radius of $\sim {10}^{17}$ cm. Chandra X-ray Observatory observations rule out alternative scenarios involving the tidal disruption of a white dwarf (WD) by a black hole, for masses >100 M_⊙. Our results challenge the notion that spectroscopically classified Ca-rich transients only originate from WD progenitor systems, complicate the view that they are all associated with large ejection velocities, and indicate that their chemical abundances may vary widely between events.

X-ray observations show that galaxy clusters have a very large range of morphologies. The most disturbed systems, which are good to study how clusters form and grow and to test physical models, may potentially complicate cosmological studies because the cluster mass determination becomes more challenging. Thus, we need to understand the cluster properties of our samples to reduce possible biases. This is complicated by the fact that different experiments may detect different cluster populations. For example, Sunyaev–Zeldovich (SZ) selected cluster samples have been found to include a greater fraction of disturbed systems than X-ray selected samples. In this paper we determine eight morphological parameters for the Planck Early Sunyaev–Zeldovich (ESZ) objects observed with XMM-Newton. We found that two parameters, concentration and centroid shift, are the best to distinguish between relaxed and disturbed systems. For each parameter we provide the values that allow selecting the most relaxed or most disturbed objects from a sample. We found that there is no mass dependence on the cluster dynamical state. By comparing our results with what was obtained with REXCESS clusters, we also confirm that the ESZ clusters indeed tend to be more disturbed, as found by previous studies.

The standard disk is often inadequate to model disk-dominated cataclysmic variables (CVs) and generates a spectrum that is bluer than the observed UV spectra. X-ray observations of these systems reveal an optically thin boundary layer (BL) expected to appear as an inner hole in the disk. Consequently, we truncate the inner disk. However, instead of removing the inner disk, we impose the no-shear boundary condition at the truncation radius, thereby lowering the disk temperature and generating a spectrum that better fits the UV data. With our modified disk, we analyze the archival UV spectra of three novalikes that cannot be fitted with standard disks. For the VY Scl systems MV Lyr and BZ Cam, we fit a hot inflated white dwarf (WD) with a cold modified disk ($\dot{M}\,\sim$ a few 10⁻⁹M_⊙ yr⁻¹). For V592 Cas, the slightly modified disk ($\dot{M}\sim 6\times {10}^{-9}\,{M}_{\odot }\,{\mathrm{yr}}^{-1}$) completely dominates the UV. These results are consistent with Swift X-ray observations of these systems, revealing BLs merged with ADAF-like flows and/or hot coronae, where the advection of energy is likely launching an outflow and heating the WD, thereby explaining the high WD temperature in VY Scl systems. This is further supported by the fact that the X-ray hardness ratio increases with the shallowness of the UV slope in a small CV sample we examine. Furthermore, for 105 disk-dominated systems, the International Ultraviolet Explorer spectra UV slope decreases in the same order as the ratio of the X-ray flux to optical/UV flux: from SU UMa's, to U Gem's, Z Cam's, UX UMa's, and VY Scl's.

The protons in large solar energetic particle events are accelerated in the inner heliosphere by fast shocks produced by coronal mass ejections. Unless there are other sources, the protons these shocks act upon would be those of the solar wind (SW). The efficiency of the acceleration depends on the kinetic energy of the protons. For a 2000 km s⁻¹ shock, the most effective proton energies would be 30–100 keV; i.e., within the suprathermal tail component of the SW. We investigate one possible additional source of such protons: those resulting from the decay of solar-flare-produced neutrons that escape from the Sun into the low corona. The neutrons are produced by interactions of flare-accelerated ions with the solar atmosphere. We discuss the production of low-energy neutrons in flares and their decay on a interplanetary magnetic field line near the Sun. We find that even when the flaring conditions are optimal, the 30–100 keV neutron-decay proton density produced by even a very large solar flare would be only about 10% of that of the 30–100 keV SW suprathermal tail. We discuss the implication of a seed-particle source of more frequent, small flares.

Internal shocks between propagating plasma shells, originally ejected at different times with different velocities, are believed to play a major role in dissipating the kinetic energy, thereby explaining the observed light curves and spectra in a large range of transient objects. Even if initially the colliding plasmas are cold, following the first collision, the plasma shells are substantially heated, implying that in a scenario of multiple collisions, most collisions take place between plasmas of non-zero temperatures. Here, we calculate the dynamical properties of plasmas resulting from a collision between arbitrarily hot plasma shells, moving at arbitrary speeds. We provide simple analytical expressions valid for both ultrarelativistic and Newtonian velocities for both hot and cold plasmas. We derive the minimum criteria required for the formation of the two-shock wave system, and show that in the relativistic limit, the minimum Lorentz factor is proportional to the square root of the ratio of the initial plasmas enthalpies. We provide basic scaling laws of synchrotron emission from both the forward and reverse-shock waves, and show how these can be used to deduce the properties of the colliding shells. Finally, we discuss the implications of these results in the study of several astronomical transients, such as X-ray binaries, radio-loud quasars, and gamma-ray bursts.

During the Space Telescope and Optical Reverberation Mapping Project observations of NGC 5548, the continuum and emission-line variability became decorrelated during the second half of the six-month-long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with the Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all three observations: the UV emission-line flux decrease, the soft-excess increase, and the emission-line anomaly.

We develop a numerical model to study the time-dependent modulation of galactic cosmic rays in the inner heliosphere. In the model, a time-delayed modified Parker heliospheric magnetic field (HMF) and a new diffusion coefficient model, NLGCE-F, from Qin & Zhang, are adopted. In addition, the latitudinal dependence of magnetic turbulence magnitude is assumed to be $\sim (1+{\sin }^{2}\theta )/2$ from the observations of Ulysses, and the radial dependence is assumed to be $\sim {r}^{S}$, where we choose an expression of S as a function of the heliospheric current sheet tilt angle. We show that the analytical expression used to describe the spatial variation of HMF turbulence magnitude agrees well with the Ulysses, Voyager 1, and Voyager 2 observations. By numerically calculating the modulation code, we get the proton energy spectra as a function of time during the recent solar minimum, it is shown that the modulation results are consistent with the Payload for Antimatter-Matter Exploration and Light-nuclei Astrophysics measurements.

We present the result of our spectroscopic follow-up observation for faint quasar candidates at z ∼ 5 in part of the Canada–France–Hawaii Telescope Legacy Survey wide field. We select nine photometric candidates and identify three z ∼ 5 faint quasars, one z ∼ 4 faint quasar, and a late-type star. Since two faint quasar spectra show the C iv emission line without suffering from a heavy atmospheric absorption, we estimate their black hole masses (${M}_{\mathrm{BH}}$) and Eddington ratios ($L/{L}_{\mathrm{Edd}}$). The inferred $\mathrm{log}{M}_{\mathrm{BH}}$ are 9.04 ± 0.14 and 8.53 ± 0.20, respectively. In addition, the inferred $\mathrm{log}(L/{L}_{\mathrm{Edd}})$ are −1.00 ± 0.15 and −0.42 ± 0.22, respectively. If we adopt that $L/{L}_{\mathrm{Edd}}=\mathrm{constant}\ \mathrm{or}\propto {(1+z)}^{2}$, the seed black hole masses (${M}_{\mathrm{seed}}$) of our z ∼ 5 faint quasars are expected to be $\gt {10}^{5}\,{M}_{\odot }$ in most cases. We also compare the observational results with a mass accretion model, where angular momentum is lost due to supernova explosions. Accordingly, ${M}_{\mathrm{BH}}$ of the z ∼ 5 faint quasars in our sample can be explained even if ${M}_{\mathrm{seed}}$ is $\sim {10}^{3}\,{M}_{\odot }$. Since z ∼ 6 luminous qusars and our z ∼ 5 faint quasars are not on the same evolutionary track, z ∼ 6 luminous quasars and our z ∼ 5 quasars are not the same populations but different populations, due to the difference of a period of the mass supply from host galaxies. Furthermore, we confirm that one can explain ${M}_{\mathrm{BH}}$ of z ∼ 6 luminous quasars and our z ∼ 5 faint quasars even if their seed black holes are formed at z ∼ 7.

We study the optical light curve (LC) relations of Type Ia supernovae (SNe Ia) for their use in cosmology using high-quality photometry published by the Carnegie Supernova Project (CSP-I). We revisit the classical luminosity decline rate (Δm₁₅) relation and the Lira relation, as well as investigate the time evolution of the (B − V) color and B(B − V), which serves as the basis of the color–stretch relation and Color–MAgnitude Intercept Calibrations (CMAGIC). Our analysis is based on explosion and radiation transport simulations for spherically symmetric delayed-detonation models (DDT) producing normal-bright and subluminous SNe Ia. Empirical LC relations can be understood as having the same physical underpinnings, i.e., opacities, ionization balances in the photosphere, and radioactive energy deposition changing with time from below to above the photosphere. Some three to four weeks past maximum, the photosphere recedes to ⁵⁶Ni-rich layers of similar density structure, leading to a similar color evolution. An important secondary parameter is the central density ρ_c of the WD because at higher densities, more electron-capture elements are produced at the expense of ⁵⁶Ni production. This results in a Δm₁₅ spread of 0.1 mag in normal-bright and 0.7 mag in subluminous SNe Ia and ≈0.2 mag in the Lira relation. We show why color–magnitude diagrams emphasize the transition between physical regimes and enable the construction of templates that depend mostly on Δm₁₅ with little dispersion in both the CSP-I sample and our DDT models. This allows intrinsic SN Ia variations to be separated from the interstellar reddening characterized by E(B − V) and R_B. Invoking different scenarios causes a wide spread in empirical relations, which may suggest one dominant scenario.

Photospheric velocities and stellar activity features such as spots and faculae produce measurable radial velocity signals that currently obscure the detection of sub-meter-per-second planetary signals. However, photospheric velocities are imprinted differently in a high-resolution spectrum than are Keplerian Doppler shifts. Photospheric activity produces subtle differences in the shapes of absorption lines due to differences in how temperature or pressure affects the atomic transitions. In contrast, Keplerian Doppler shifts affect every spectral line in the same way. With a high enough signal-to-noise (S/N) and resolution, statistical techniques can exploit differences in spectra to disentangle the photospheric velocities and detect lower-amplitude exoplanet signals. We use simulated disk-integrated time-series spectra and principal component analysis (PCA) to show that photospheric signals introduce spectral line variability that is distinct from that of Doppler shifts. We quantify the impact of instrumental resolution and S/N for this work.

Line-intensity mapping surveys probe large-scale structure through spatial variations in molecular line emission from a population of unresolved cosmological sources. Future such surveys of carbon monoxide line emission, specifically the CO(1-0) line, face potential contamination from a disjointed population of sources emitting in a hydrogen cyanide emission line, HCN(1-0). This paper explores the potential range of the strength of HCN emission and its effect on the CO auto power spectrum, using simulations with an empirical model of the CO/HCN–halo connection. We find that effects on the observed CO power spectrum depend on modeling assumptions but are very small for our fiducial model, which is based on current understanding of the galaxy–halo connection. Given the fiducial model, we expect the bias in overall CO detection significance due to HCN to be less than 1%.

We study the dependence of galaxy clustering on H i mass using ∼16,000 galaxies with redshift in the range of $0.0025\lt z\lt 0.05$ and H i mass of ${M}_{{\rm{H}}{\rm{I}}}\gt {10}^{8}\,{M}_{\odot }$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with ${M}_{{\rm{H}}{\rm{I}}}$ above different thresholds and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak H i mass dependence, we find both the clustering amplitudes on scales above a few megaparsecs and the bias factors to increase significantly with increasing H i mass for ${M}_{{\rm{H}}{\rm{I}}}\gt {10}^{9}\,{M}_{\odot }$. For H i mass thresholds below $\sim {10}^{9}\,{M}_{\odot }$, the inferred galaxy bias factors are systematically lower than the minimum halo bias from mass-selected halo samples. We extend the simple halo model, in which the galaxy content is only determined by halo mass, by including the halo formation time as an additional parameter. A model that puts H i-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of H i mass with halo mass and formation time, as well as the halo occupation distributions and H i mass functions for central and satellite galaxies. These results are compared with the predictions from semianalytic galaxy formation models and hydrodynamic galaxy formation simulations.

The electromagnetic (EM) follow-up of a gravitational-wave (GW) event requires scanning a wide sky region, defined by the so-called "skymap," to detect and identify a transient counterpart. We propose a novel method that exploits the information encoded in the GW signal to construct a "detectability map," which represents the time-dependent ("when") probability of detecting the transient at each position of the skymap ("where"). Focusing on the case of a neutron star binary inspiral, we model the associated short gamma-ray burst afterglow and macronova emission using the probability distributions of binary parameters (sky position, distance, orbit inclination, mass ratio) extracted from the GW signal as inputs. The resulting family of possible light curves is the basis for constructing the detectability map. As a practical example, we apply the method to a simulated GW signal produced by a neutron star merger at 75 Mpc whose localization uncertainty is very large (∼1500 deg²). We construct observing strategies for optical, infrared, and radio facilities based on the detectability maps, taking VST, VISTA, and MeerKAT as prototypes. Assuming limiting fluxes of $r\sim 24.5$, $J\sim 22.4$ (AB magnitudes), and $500\,\mu \mathrm{Jy}$ ($1.4\,\mathrm{GHz}$) for ∼1000 s of exposure each, the afterglow and macronova emissions are successfully detected with a minimum observing time of 7, 15, and 5 hr respectively.

With 7 years of Interstellar Boundary Explorer (IBEX) measurements of energetic neutral atoms (ENAs), IBEX has shown a clear correlation between dynamic changes in the solar wind and the heliosphere's response in the formation of ENAs. In this paper, we investigate temporal variations in the latitudinal-dependent ENA spectrum from IBEX and their relationship to the solar wind speed observed at 1 au. We find that the variation in latitude of the transition in ENA spectral indices between low (≲1.8) and high (≳1.8) values, as well as the distribution of ENA spectral indices at high and low latitudes, correlates well with the evolution of the fast and slow solar wind latitudinal structure observed near 1 au. This correlation includes a delay due to the time it takes the solar wind to propagate to the termination shock and into the inner heliosheath, and for ENAs to be generated via charge-exchange and travel back toward 1 au. Moreover, we observe a temporal asymmetry in the steepening of the ENA spectrum in the northern and southern hemispheres, consistent with asymmetries observed in the solar wind and polar coronal holes. While this asymmetry is observed near the upwind direction of the heliosphere, it is not yet observed in the tail direction, suggesting a longer line-of-sight integration distance or different processing of the solar wind plasma downstream of the termination shock.

We derive infrared and radio flux densities of all ∼1000 known Galactic H ii regions in the Galactic longitude range $17\buildrel{\circ}\over{.} 5\lt {\ell }\lt 65^\circ$. Our sample comes from the Wide-Field Infrared Survey Explorer (WISE) catalog of Galactic H ii regions. We compute flux densities at six wavelengths in the infrared (Spitzer GLIMPSE 8 μm, WISE 12 μm and 22 μm, Spitzer MIPSGAL 24 μm, and Herschel Hi-GAL 70 μm and 160 μm) and two in the radio (MAGPIS 20 cm and VGPS 21 cm). All H ii region infrared flux densities are strongly correlated with their ∼20 cm flux densities. All H ii regions used here, regardless of physical size or Galactocentric radius, have similar infrared to radio flux density ratios and similar infrared colors, although the smallest regions (r < 1 pc), have slightly elevated IR to radio ratios. The colors ${\mathrm{log}}_{10}({F}_{24\mu {\rm{m}}}/{F}_{12\mu {\rm{m}}})\geqslant 0$ and ${\mathrm{log}}_{10}({F}_{70\mu {\rm{m}}}/{F}_{12\mu {\rm{m}}})\geqslant 1.2$, and ${\mathrm{log}}_{10}({F}_{24\mu {\rm{m}}}/{F}_{12\mu {\rm{m}}})\geqslant 0$ and ${\mathrm{log}}_{10}({F}_{160\mu {\rm{m}}}/{F}_{70\mu {\rm{m}}})\leqslant 0.67$ reliably select H ii regions, independent of size. The infrared colors of ∼22% of H ii regions, spanning a large range of physical sizes, satisfy the IRAS color criteria of Wood & Churchwell for H ii regions, after adjusting the criteria to the wavelengths used here. Because these color criteria are commonly thought to select only ultra-compact H ii regions, this result indicates that the true ultra-compact H ii region population is uncertain. Compared to a sample of IR color indices from star-forming galaxies, H ii regions show higher ${\mathrm{log}}_{10}({F}_{70\mu {\rm{m}}}/{F}_{12\mu {\rm{m}}})$ ratios. We find a weak trend of decreasing infrared to ∼20 cm flux density ratios with increasing R_gal, in agreement with previous extragalactic results, possibly indicating a decreased dust abundance in the outer Galaxy.

Crucial information on nova nucleosynthesis can be potentially inferred from γ-ray signals powered by ¹⁸F decay. Therefore, the reaction network producing and destroying this radioactive isotope has been extensively studied in the last years. Among those reactions, the ¹⁸F(p, α)¹⁵O cross-section has been measured by means of several dedicated experiments, both using direct and indirect methods. The presence of interfering resonances in the energy region of astrophysical interest has been reported by many authors including the recent applications of the Trojan Horse Method. In this work, we evaluate what changes are introduced by the Trojan Horse data in the ¹⁸F(p, α)¹⁵O astrophysical factor recommended in a recent R-matrix analysis, accounting for existing direct and indirect measurements. Then the updated reaction rate is calculated and parameterized and implications of the new results on nova nucleosynthesis are thoroughly discussed.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way (MW) as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the MW. Using Markov Chain Monte Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass–metallicity relation of the MW satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass–metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g., "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

High-energy γ-rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are well understood and are reproduced by existing CR propagation models employing 2D galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degrees scales, indicating the need to include the details of the Galactic spiral structure and thus requiring 3D spatial modeling. In this paper, the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disk or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy γ-ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the nonthermal interstellar emissions from the Galaxy.

We have compiled a sample of 26 metal-poor galaxies with 12 + log(O/H) < 8.1 with both infrared continuum and 1.4 GHz radio continuum data. By comparing to galaxies at higher metallicity, we have investigated the IR–radio relationship's dependence on metallicity at the 24, 70, 100, and 160 μm bands, as well as the integrated FIR luminosity. It is found that metal-poor galaxies have on average lower ${q}_{\mathrm{IR}}$ than metal-rich ones, with larger offsets at longer IR wavelengths, from −0.06 dex in ${q}_{24\mu {\rm{m}}}$ to −0.6 dex in ${q}_{160\mu {\rm{m}}}$. The ${q}_{\mathrm{IR}}$ of all galaxies as a whole at 160 μm show positive trends with the metallicity and IR-to-FUV ratio and negative trends with the IR color, while those at lower IR wavelengths show weaker correlations. We propose a mechanism that invokes the combined effects of low obscured-SFR-to-total-SFR fraction and warm dust temperature at low metallicity to interpret the above behavior of ${q}_{\mathrm{IR}}$, with the former reducing the IR radiation and the latter further reducing the IR emission at longer IR wavelengths. Other mechanisms that are related to the radio emission, including the enhanced magnetic field strength and increased thermal radio contribution, are unable to reconcile the IR-wavelength-dependent differences of ${q}_{\mathrm{IR}}$ between metal-poor and metal-rich galaxies. In contrast to ${q}_{\mathrm{IR}}$, the mean total-SFR-to-radio ratio of metal-poor galaxies is the same as that for metal-rich galaxies, indicating the 1.4 GHz radio emission is still an effective tracer of SFRs at low metallicity.

We present three-dimensional atmospheric circulation models of a hypothetical "warm Jupiter" planet, for a range of possible obliquities from 0° to 90°. We model a Jupiter-mass planet on a 10 day orbit around a Sun-like star, since this hypothetical planet sits at the boundary between planets for which we expect that tidal forces should have aligned their rotation axes with their orbital axes (i.e., ones with zero obliquity) and planets whose timescale for tidal alignment is longer than the typical age of an exoplanet system. In line with observational progress, which is pushing atmospheric characterization for planets on longer orbital periods, we calculate the observable signatures of obliquity for a transiting warm Jupiter: in orbital phase curves of thermal emission and in the hemispheric flux gradients that could be measured by eclipse mapping. For both of these predicted measurements, the signal that we would see depends strongly on our viewing geometry relative to the orientation of the planet's rotation axis, and we thoroughly identify the degeneracies that result. We compare these signals to the predicted sensitivities of current and future instruments and determine that the James Webb Space Telescope should be able to constrain the obliquities of nearby warm Jupiters to be small (if $\leqslant 10^\circ$) or to directly measure them if significantly non-zero ($\geqslant 30^\circ$) using the technique of eclipse mapping. For a bright target and assuming photon-limited precision, this could be done with a single secondary eclipse observation.

Disks around brown dwarfs (BDs) are excellent laboratories to study the first steps of planet formation in cold and low-mass disk conditions. The radial-drift velocities of dust particles in BD disks higher than in disks around more massive stars. Therefore, BD disks are expected to be more depleted in millimeter-sized grains compared to disks around T Tauri or Herbig Ae/Be stars. However, recent millimeter observations of BD disks revealed low millimeter spectral indices, indicating the presence of large grains in these disks and challenging models of dust evolution. We present 3 mm photometric observations carried out with the IRAM/Plateau de Bure Interferometer (PdBI) of three BD disks in the Taurus star-forming region, which have been observed with ALMA at 0.89 mm. The disks were not resolved and only one was detected with enough confidence (∼3.5σ) with PdBI. Based on these observations, we obtain the values and lower limits of the spectral index and find low values (α_mm ≲ 3.0). We compare these observations in the context of particle trapping by an embedded planet, a promising mechanism to explain the observational signatures in more massive and warmer disks. We find, however, that this model cannot reproduce the current millimeter observations for BD disks, and multiple-strong pressure bumps globally distributed in the disk remain as a favorable scenario to explain observations. Alternative possibilities are that the gas masses in the BD disk are very low (∼2 × 10⁻³M_Jup) such that the millimeter grains are decoupled and do not drift, or fast growth of fluffy aggregates.

We compare the structure of molecular gas at 40 pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into 370 pc and 1.1 kpc resolution elements, and within each we estimate the molecular gas depletion time (${\tau }_{\mathrm{Dep}}^{\mathrm{mol}}$), the star-formation efficiency per free-fall time (${\epsilon }_{\mathrm{ff}}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, Σ, line width, σ, and $b\equiv {\rm{\Sigma }}/{\sigma }^{2}\propto {\alpha }_{\mathrm{vir}}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star-formation efficiency per free-fall time, ${\epsilon }_{\mathrm{ff}}(\langle {{\rm{\Sigma }}}_{40\mathrm{pc}}\rangle )\sim 0.3 \% \mbox{--}0.36 \%$, lower than adopted in many models and found for local clouds. Furthermore, the efficiency per free-fall time anti-correlates with both Σ and σ, in some tension with turbulent star-formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b\equiv {\rm{\Sigma }}/{\sigma }^{2}$, tracing the strength of self-gravity, with ${\tau }_{\mathrm{Dep}}^{\mathrm{mol}}\propto {b}^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher b) forms stars at a higher rate (low ${\tau }_{\mathrm{Dep}}^{\mathrm{mol}}$). The different regions of the galaxy mostly overlap in ${\tau }_{\mathrm{Dep}}^{\mathrm{mol}}$ as a function of b, so that low b explains the surprisingly high ${\tau }_{\mathrm{Dep}}^{\mathrm{mol}}$ found toward the inner spiral arms found by Meidt et al. (2013).

Recent determinations of the radial distributions of mono-metallicity populations (MMPs, i.e., stars in narrow bins in [Fe/H] within wider [α/Fe] ranges) by the SDSS-III/APOGEE DR12 survey cast doubts on the classical thin- and thick-disk dichotomy. The analysis of these observations led to the non-$[\alpha$/Fe] enhanced populations splitting into MMPs with different surface densities according to their [Fe/H]. By contrast, $[\alpha$/Fe] enhanced (i.e., old) populations show a homogeneous behavior. We analyze these results in the wider context of disk formation within non-isolated halos embedded in the Cosmic Web, resulting in a two-phase mass assembly. By performing hydrodynamical simulations in the context of the ΛCDM model, we have found that the two phases of halo mass assembly (an early fast phase, followed by a slow phase with low mass-assembly rates) are very relevant to determine the radial structure of MMP distributions, while radial mixing only plays a secondary role, depending on the coeval dynamical and/or destabilizing events. Indeed, while the frequent dynamical violent events occuring at high redshift remove metallicity gradients and imply efficient stellar mixing, the relatively quiescent dynamics after the transition keeps [Fe/H] gaseous gradients and prevents newly formed stars from suffering strong radial mixing. By linking the two-component disk concept with the two-phase halo mass-assembly scenario, our results set halo virialization (the event marking the transition from the fast to the slow phases) as the separating event that marks periods that are characterized by different physical conditions under which thick- and thin-disk stars were born.

We present observations of CO(3–2) and ¹³CO(3–2) emission near the supernebula in the dwarf galaxy NGC 5253, which contains one of the best examples of a potential globular cluster in formation. The 03 resolution images reveal an unusual molecular cloud, "Cloud D1," that is coincident with the radio-infrared supernebula. The ∼6 pc diameter cloud has a linewidth, Δ v = 21.7 $\mathrm{km}\,{{\rm{s}}}^{-1}$, that reflects only the gravitational potential of the star cluster residing within it. The corresponding virial mass is 2.5 × 10⁵${M}_{\odot }$. The cluster appears to have a top-heavy initial mass function, with M_* ≳ 1–2 ${M}_{\odot }$. Cloud D1 is optically thin in CO(3–2), probably because the gas is hot. Molecular gas mass is very uncertain but constitutes <35% of the dynamical mass within the cloud boundaries. In spite of the presence of an estimated ∼1500–2000 O stars within the small cloud, the CO appears relatively undisturbed. We propose that Cloud D1 consists of molecular clumps or cores, possibly star-forming, orbiting with more evolved stars in the core of the giant cluster.

Interactions between pairs of isolated dwarf galaxies provide a critical window into low-mass hierarchical, gas-dominated galaxy assembly and the build-up of stellar mass in low-metallicity systems. We present the first Very Large Telescope/Multi Unit Spectroscopic Explorer (VLT/MUSE) optical integral field unit (IFU) observations of the interacting dwarf pair dm1647+21 selected from the TiNy Titans survey. The Hα emission is widespread and corresponds to a total unobscured star formation rate (SFR) of 0.44 M_⊙ yr⁻¹, which is 2.7 times higher than the SFR inferred from Sloan Digital Sky Survey (SDSS) data. The implied specific SFR (sSFR) for the system is elevated by more than an order of magnitude above non-interacting dwarfs in the same mass range. This increase is dominated by the lower-mass galaxy, which has a sSFR enhancement of >50. Examining the spatially resolved maps of classic optical line diagnostics, we find that the interstellar medium (ISM) excitation can be fully explained by star formation. The velocity field of the ionized gas is not consistent with simple rotation. Dynamical simulations indicate that the irregular velocity field and the stellar structure is consistent with the identification of this system as an ongoing interaction between two dwarf galaxies. The widespread, clumpy enhancements in the star formation in this system point to important differences in the effect of mergers on dwarf galaxies, compared to massive galaxies; rather than the funneling of gas to the nucleus and giving rise to a nuclear starburst, starbursts in low-mass galaxy mergers may be triggered by large-scale ISM compression, and thus may be more distributed.

Stellar coronal activity has been shown to persist into the low-mass star regime, down to late M-dwarf spectral types. However, there is now an accumulation of evidence suggesting that at the end of the main sequence, there is a transition in the nature of the magnetic activity from chromospheric and coronal to planet-like and auroral, from local impulsive heating via flares and MHD wave dissipation to energy dissipation from strong large-scale magnetospheric current systems. We examine this transition and the prevalence of auroral activity in brown dwarfs through a compilation of multiwavelength surveys of magnetic activity, including radio, X-ray, and optical. We compile the results of those surveys and place their conclusions in the context of auroral emission as a consequence of large-scale magnetospheric current systems that accelerate energetic electron beams and drive the particles to impact the cool atmospheric gas. We explore the different manifestations of auroral phenomena, like Hα, in brown dwarf atmospheres and define their distinguishing characteristics. We conclude that large-amplitude photometric variability in the near-infrared is most likely a consequence of clouds in brown dwarf atmospheres, but that auroral activity may be responsible for long-lived stable surface features. We report a connection between auroral Hα emission and quiescent radio emission in electron cyclotron maser instability pulsing brown dwarfs, suggesting a potential underlying physical connection between quiescent and auroral emissions. We also discuss the electrodynamic engines powering brown dwarf aurorae and the possible role of satellites around these systems both to power the aurorae and seed the magnetosphere with plasma.

The Wolf–Rayet (WR) nebula NGC 3199 has been suggested to be a bow shock around its central star, WR 18, which is presumably a runaway star, because optical images of the nebula show a dominating arc of emission southwest of the star. We present the XMM-Newton detection of extended X-ray emission from NGC 3199, unveiling the powerful effect of the fast wind from WR 18. The X-ray emission is brighter in the region southeast of the star and an analysis of the spectral properties of the X-ray emission reveals abundance variations: (i) regions close to the optical arc present nitrogen-rich gas enhanced by the stellar wind from WR 18 and (ii) gas at the eastern region exhibits abundances close to those reported for the nebular abundances derived from optical studies, which is a signature of an efficient mixing of the nebular material with the stellar wind. The dominant plasma temperature and electron density are estimated to be T ≈ 1.2 × 10⁶ K and n_e = 0.3 cm⁻³ with an X-ray luminosity in the 0.3–3.0 keV energy range of L_X = 2.6 × 10³⁴ erg s⁻¹. Combined with information derived from Herschel and the recent Gaia first data release, we conclude that WR 18 is not a runaway star and that the formation, chemical variations, and the shape of NGC 3199 depend on the initial configuration of the interstellar medium.

The mixed morphology class of supernova remnants has centrally peaked X-ray emission along with a shell-like morphology in radio emission. White & Long proposed that these remnants are evolving in a cloudy medium wherein the clouds are evaporated via thermal conduction once being overrun by the expanding shock. Their analytical model made detailed predictions regarding temperature, density, and emission profiles as well as shock evolution. We present numerical hydrodynamical models in 2D and 3D including thermal conduction, testing the White & Long model and presenting results for the evolution and emission from remnants evolving in a cloudy medium. We find that, while certain general results of the White & Long model hold, such as the way the remnants expand and the flattening of the X-ray surface brightness distribution, in detail there are substantial differences. In particular we find that the X-ray luminosity is dominated by emission from shocked cloud gas early on, leading to a bright peak, which then declines and flattens as evaporation becomes more important. In addition, the effects of thermal conduction on the intercloud gas, which is not included in the White & Long model, are important and lead to further flattening of the X-ray brightness profile as well as lower X-ray emission temperatures.

Following previous work, we further confirm that the cosmic evolution of steep-spectrum radio-loud AGNs (active galactic nuclei) can be reproduced by a simple combination of density evolution (DE) and luminosity evolution (LE). This mixture evolution scenario can naturally explain the luminosity-dependent evolution of radio-loud AGNs. Our models successfully fitted a large amount of data on radio luminosity functions of steep-spectrum sources and multi-frequency source counts. The modeling indicates that the DE slowly increases as ${(1+z)}^{0.3\sim 1.3}$ out to $z\sim 0.8$, and then rapidly decreases as ${(1+z)}^{-6.8\sim -5.7}$, while the LE rapidly increases as ${(1+z)}^{4.8}$ out to a higher redshift (at least $z\gt 3.5$). We find a high-redshift decline (i.e., redshift cutoff) in the number density of steep-spectrum radio sources, but we cannot conclude whether such a decline is sharp or shallow. We believe that whether a redshift cutoff occurs or not depends mainly on DE, while its steepness is decided by LE, which, however, cannot be well constrained due to the lack of high-redshift samples. Most intriguingly, according to our mixture evolution scenario, there appears to be no need for different evolution for the low- and high-power radio-loud AGNs. Both types of sources experience the same combined evolution of DE and LE.

We present composite broad-line region (BLR) reverberation mapping lag measurements for Hα, Hβ, He iiλ4686, and Mg ii for a sample of 144, z ≲ 1 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Using only the 32-epoch spectroscopic light curves in the first six-month season of SDSS-RM observations, we compile correlation function measurements for individual objects and then coadd them to allow the measurement of the average lags for our sample at mean redshifts of 0.4 (for Hα) and ∼0.65 (for the other lines). At similar quasar luminosities and redshifts, the sample-averaged lag decreases in the order of Mg ii, Hα, Hβ, and He ii. This decrease in lags is accompanied by an increase in the mean line width of the four lines, and is roughly consistent with the virialized motion for BLR gas in photoionization equilibrium. These are among the first RM measurements of stratified BLR structure at z > 0.3. Dividing our sample by luminosity, Hα shows clear evidence of increasing lags with luminosity, consistent with the expectation from the measured BLR size–luminosity relation based on Hβ. The other three lines do not show a clear luminosity trend in their average lags due to the limited dynamic range of luminosity probed and the poor average correlation signals in the divided samples, a situation that will be improved with the incorporation of additional photometric and spectroscopic data from SDSS-RM. We discuss the utility and caveats of composite lag measurements for large statistical quasar samples with reverberation mapping data.

We undertook coordinated campaigns with the Green Bank, Effelsberg, and Arecibo radio telescopes during Chandra X-ray Observatory and XMM-Newton observations of the repeating fast radio burst FRB 121102 to search for simultaneous radio and X-ray bursts. We find 12 radio bursts from FRB 121102 during 70 ks total of X-ray observations. We detect no X-ray photons at the times of radio bursts from FRB 121102 and further detect no X-ray bursts above the measured background at any time. We place a 5σ upper limit of 3 × 10⁻¹¹ erg cm⁻² on the 0.5–10 keV fluence for X-ray bursts at the time of radio bursts for durations $\lt 700$ ms, which corresponds to a burst energy of 4 × 10⁴⁵ erg at the measured distance of FRB 121102. We also place limits on the 0.5–10 keV fluence of 5 × 10⁻¹⁰ and 1 × 10⁻⁹ erg cm⁻² for bursts emitted at any time during the XMM-Newton and Chandra observations, respectively, assuming a typical X-ray burst duration of 5 ms. We analyze data from the Fermi Gamma-ray Space Telescope Gamma-ray Burst Monitor and place a 5σ upper limit on the 10–100 keV fluence of 4 × 10⁻⁹ erg cm⁻² (5 × 10⁴⁷ erg at the distance of FRB 121102) for gamma-ray bursts at the time of radio bursts. We also present a deep search for a persistent X-ray source using all of the X-ray observations taken to date and place a 5σ upper limit on the 0.5–10 keV flux of 4 × 10⁻¹⁵ erg s⁻¹ cm⁻² (3 × 10⁴¹ erg s⁻¹ at the distance of FRB 121102). We discuss these non-detections in the context of the host environment of FRB 121102 and of possible sources of fast radio bursts in general.

The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models have focused on the formation of the satellites starting from satellitesimals, the question of how satellitesimals themselves form from smaller dust particles has not yet been addressed. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust-to-gas accretion rates onto the circumplanetary disk, ${\dot{M}}_{{\rm{d}}}/{\dot{M}}_{{\rm{g}}}$, is high and the strength of turbulence, α, is not too low. The collision velocity is lower than the critical velocity of fragmentation when α is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by ${\dot{M}}_{{\rm{d}}}/{\dot{M}}_{{\rm{g}}}\geqslant 1$ and ${10}^{-4}\leqslant \alpha \lt {10}^{-2}$. The former condition is generally difficult to achieve, suggesting that the in situ satellitesimal formation via particle sticking is viable only under extreme conditions. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as ${\dot{M}}_{{\rm{d}}}/{\dot{M}}_{{\rm{g}}}$ is low.

Gravitational waves (GWs) from binary black hole (BBH) mergers provide a new probe of massive-star evolution and the formation channels of binary compact objects. By coupling the growing sample of BBH systems with population synthesis models, we can begin to constrain the parameters of such models and glean unprecedented knowledge about the inherent physical processes that underpin binary stellar evolution. In this study, we apply a hierarchical Bayesian model to mass measurements from a synthetic GW sample to constrain the physical prescriptions in population models and the relative fraction of systems generated from various channels. We employ population models of two canonical formation scenarios in our analysis—isolated binary evolution involving a common-envelope phase and dynamical formation within globular clusters—with model variations for different black hole natal kick prescriptions. We show that solely with chirp mass measurements, it is possible to constrain natal kick prescriptions and the relative fraction of systems originating from each formation channel with ${ \mathcal O }(100)$ of confident detections. This framework can be extended to include additional formation scenarios, model parameters, and measured properties of the compact binary.

The discovery of long-lived electrostatic coherent structures with large-amplitude electric fields ($1\leqslant E\,\leqslant 500$ mV/m) by the Van Allen Probes has revealed alternative routes through which planetary radiation belts' acceleration can take place. Following previous reports showing that small phase-space holes, with $q\phi /{T}_{e}^{c}\simeq {10}^{-2}\mbox{--}{10}^{-3}$, could result from electron interaction with large-amplitude whistlers, we demonstrate one possible mechanism through which holes can grow nonlinearly (i.e., $\gamma \propto \sqrt{\phi }$) and subcritically as a result of momentum exchange between hot and cold electron populations. Our results provide an explanation for the common occurrence and fast growth of large-amplitude electron phase-space holes in the Earth's radiation belts.

The Faint Infrared Grism Survey (FIGS) is a deep Hubble Space Telescope (HST) WFC3/IR (Wide Field Camera 3 Infrared) slitless spectroscopic survey of four deep fields. Two fields are located in the Great Observatories Origins Deep Survey-North (GOODS-N) area and two fields are located in the Great Observatories Origins Deep Survey-South (GOODS-S) area. One of the southern fields selected is the Hubble Ultra Deep Field. Each of these four fields were observed using the WFC3/G102 grism (0.8 μm–1.15 μm continuous coverage) with a total exposure time of 40 orbits (≈100 kilo-seconds) per field. This reaches a $3\sigma$ continuum depth of $\approx 26$ AB magnitudes and probes emission lines to $\sim {10}^{-17}\,\mathrm{erg}\,{{\rm{s}}}^{-1}\,{\mathrm{cm}}^{-2}$. This paper details the four FIGS fields and the overall observational strategy of the project. A detailed description of the Simulation Based Extraction (SBE) method used to extract and combine over 10,000 spectra of over 2000 distinct sources brighter than ${m}_{F105W}=26.5$ mag is provided. High fidelity simulations of the observations is shown to significantly improve the background subtraction process, the spectral contamination estimates, and the final flux calibration. This allows for the combination of multiple spectra to produce a final high quality, deep, 1D spectra for each object in the survey.

We present the experimental phase functions of three types of millimeter-sized dust grains consisting of enstatite, quartz, and volcanic material from Mount Etna, respectively. The three grains present similar sizes but different absorbing properties. The measurements are performed at 527 nm covering the scattering angle range from 3° to 170°. The measured phase functions show two well-defined regions: (i) soft forward peaks and (ii) a continuous increase with the scattering angle at side- and back-scattering regions. This behavior at side- and back-scattering regions is in agreement with the observed phase functions of the Fomalhaut and HR 4796A dust rings. Further computations and measurements (including polarization) for millimeter-sized grains are needed to draw some conclusions about the fluffy or compact structure of the dust grains.

This work combined coronagraphic visible light (VL) and UV data to provide with an unprecedented view of the inner corona where the nascent solar wind is accelerated. The UV (H i Lyα) and VL (polarized brightness) images (reconstructed with SOHO/UVCS, LASCO, and Mauna Loa data) have been analyzed with the Doppler dimming technique to provide for the first time daily 2D images of the radial wind speed between 1 and 6 R_⊙ over 1 month of observations. Results show that both polar and equatorial regions are characterized at the base of the corona by plasma outflows at speeds $\gt 100$ km s⁻¹. The plasma is then decelerated within ∼1.5 R_⊙ at the poles and ∼2.0 R_⊙ at the equator, where local minima of the expansion speeds are reached, and gently reaccelerated higher up, reaching speeds typical of fast and slow wind components. The mass flux is highly variable with latitude and time at the equator and more uniform and stable over the poles. The polar flow is asymmetric, with speeds above the south pole lower than those above the north pole. A correlation (anticorrelation) between the wind speed and its density is found below (above) ∼1.8 R_⊙. The 2D distribution of forces responsible for deceleration and reacceleration of solar wind is provided and interpreted in terms of Alfvén waves. These results provide a possible connection between small-scale outflows reported with other instruments at the base of the corona and bulk wind flows measured higher up.