The 100 most-cited articles in The Astrophysical Journal, 2021

Abstract Measurement of the distances to nearby galaxies has improved rapidly in recent decades. The ever-present challenge is to reduce systematic effects, especially as greater distances are probed and the uncertainties become larger. In this paper, we combine several recent calibrations of the tip of the red giant branch (TRGB) method. These calibrations are internally self-consistent at the 1% level. New Gaia Early Data Release 3 data provide an additional consistency check at a (lower) 5% level of accuracy, a result of the well-documented Gaia angular covariance bias. The updated TRGB calibration applied to a sample of Type Ia supernovae from the Carnegie Supernova Project results in a value of the Hubble constant of H 0 = 69.8 ± 0.6 (stat) ± 1.6 (sys) km s−1 Mpc−1. No statistically significant difference is found between the value of H 0 based on the TRGB and that determined from the cosmic microwave background. The TRGB results are also consistent to within 2σ with the SHoES and Spitzer plus Hubble Space Telescope (HST) Key Project Cepheid calibrations. The TRGB results alone…

Abstract While tidal disruption events (TDEs) have long been heralded as laboratories for the study of quiescent black holes, the small number of known TDEs and uncertainties in their emission mechanism have hindered progress toward this promise. Here we present 17 new TDEs that have been detected recently by the Zwicky Transient Facility along with Swift UV and X-ray follow-up observations. Our homogeneous analysis of the optical/UV light curves, including 22 previously known TDEs from the literature, reveals a clean separation of light-curve properties with spectroscopic class. The TDEs with Bowen fluorescence features in their optical spectra have smaller blackbody radii, lower optical luminosities, and higher disruption rates compared to the rest of the sample. The small subset of TDEs that show only helium emission lines in their spectra have the longest rise times, the highest luminosities, and the lowest rates. A high detection rate of Bowen lines in TDEs with small photometric radii could be explained by the high density that is required for this fluorescence mechanism. The…

Abstract We present version 10 of the CHIANTI package. In this release, we provide updated atomic models for several helium-like ions and for all the ions of the beryllium, carbon, and magnesium isoelectronic sequences that are abundant in astrophysical plasmas. We include rates from large-scale atomic structure and scattering calculations that are in many cases a significant improvement over the previous version, especially for the Be-like sequence, which has useful line diagnostics to measure the electron density and temperature. We have also added new ions and updated several of them with new atomic rates and line identifications. Also, we have added several improvements to the IDL software, to speed up the calculations and to estimate the suppression of dielectronic recombination.

Abstract The Hubble constant (H 0) tension between Type Ia supernovae (SNe Ia) and Planck measurements ranges from 4 to 6σ. To investigate this tension, we estimate H 0 in the ΛCDM and w 0 w a CDM (cold dark matter) models by dividing the Pantheon sample, the largest compilation of SNe Ia, into 3, 4, 20, and 40 bins. We fit the extracted H 0 values with a function mimicking the redshift evolution: g ( z ) = H 0 ( z ) = H ˜ 0 / ( 1 + z ) α , where α indicates an evolutionary parameter and H ˜ 0 = H 0 at z = 0. We set the absolute magnitude of SNe Ia so that H 0 = 73.5 km s − 1 Mpc − 1 , and we fix fiducial values for Ω 0 m Λ CDM = 0.298 and Ω 0 m w 0 w

Abstract We present a new generation of substellar atmosphere and evolution models, appropriate for application to studies of L-, T-, and Y-type brown dwarfs and self-luminous extrasolar planets. The models describe the expected temperature-pressure profiles and emergent spectra of atmospheres in radiative-convective equilibrium with effective temperatures and gravities within the ranges 200 ≤ T eff ≤ 2400 K and 2.5 ≤ log g ≤ 5.5 . These ranges encompass masses from about 0.5 to 85 Jupiter masses for a set of metallicities ([M/H] = − 0.5 to + 0.5), C/O ratios (from 0.5 to 1.5 times that of solar), and ages. These models expand the diversity of model atmospheres currently available, notably to cooler effective temperatures and greater ranges in C/O. Notable improvements from past such models include updated opacities and atmospheric chemistry. Here we describe our modeling approach and present our initial tranche of models for cloudless, chemical equilibrium atmospheres. We compare the modeled spectra, photometry, and evolution to various data sets.

Abstract The second LIGO–Virgo catalog of gravitational-wave (GW) transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes that impact the observational properties of mergers. Given our set of formation models, we find significant differences between the branching fractions of the underlying and detectable populations, and the diversity of detections suggests that multiple formation channels are at play. A mixture of channels is strongly preferred over any single channel dominating the detected population: an individual channel does not contribute to more than ≃70% of the observational sample of binary black holes. We calculate the preference between the natal spin assumptions and common-envelope efficiencies in our models, favoring natal spins of isolated black holes of ≲0.1 and marginally…

Abstract We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. CAMELS is a suite of 4233 cosmological simulations of 25 h − 1 Mpc 3 volume each: 2184 state-of-the-art (magneto)hydrodynamic simulations run with the AREPO and GIZMO codes, employing the same baryonic subgrid physics as the IllustrisTNG and SIMBA simulations, and 2049 N-body simulations. The goal of the CAMELS project is to provide theory predictions for different observables as a function of cosmology and astrophysics, and it is the largest suite of cosmological (magneto)hydrodynamic simulations designed to train machine-learning algorithms. CAMELS contains thousands of different cosmological and astrophysical models by way of varying Ω m , σ 8, and four parameters controlling stellar and active galactic nucleus feedback, following the evolution of more than 100 billion particles and fluid elements over a combined volume of ( 400 h − 1 Mpc ) 3 . We describe the simulations in detail and characterize the large range of conditions represented in terms of the matter power spectrum,…

Abstract This paper presents the gravitational-wave measurement of the Hubble constant ( H 0 ) using the detections from the first and second observing runs of the Advanced LIGO and Virgo detector network. The presence of the transient electromagnetic counterpart of the binary neutron star GW170817 led to the first standard-siren measurement of H 0 . Here we additionally use binary black hole detections in conjunction with galaxy catalogs and report a joint measurement. Our updated measurement is H 0 = 69 − 8 + 16 km s −1 Mpc −1 (68.3% of the highest density posterior interval with a flat-in-log prior) which is an improvement by a factor of 1.04 (about 4%) over the GW170817-only value of 69 − 8 + 17 km s −1 Mpc −1 . A significant additional contribution currently comes from GW170814, a loud and well-localized detection from a part of the sky thoroughly covered by the Dark Energy Survey. With numerous detections anticipated over the upcoming years, an exhaustive understanding of other systematic effects are also going to become increasingly important. These results establish the path…

Abstract We investigate the evolution of galaxy gas-phase metallicity (O/H) over the range z = 0–3.3 using samples of ∼300 galaxies at z ∼ 2.3 and ∼150 galaxies at z ∼ 3.3 from the MOSDEF survey. This analysis crucially utilizes different metallicity calibrations at z ∼ 0 and z > 1 to account for evolving interstellar medium (ISM) conditions. We find significant correlations between O/H and stellar mass ( M * ) at z ∼ 2.3 and z ∼ 3.3. The low-mass power-law slope of the mass–metallicity relation (MZR) is remarkably invariant over z = 0–3.3, such that O/H ∝ at all redshifts in this range. At fixed M * , O/H decreases with increasing redshift as d log(O/H)/ dz = −0.11 ± 0.02. We find no evidence that the fundamental metallicity relation between M * , O/H, and star formation rate evolves out to z ∼ 3.3. We employ analytic chemical evolution models to place constraints on the mass and metal loading factors of galactic outflows. The efficiency of metal removal increases toward lower M * at fixed redshift and toward higher redshift at fixed M * . These models suggest that the slope of the…

Abstract We present a synthesis of fast radio burst (FRB) morphology (the change in flux as a function of time and frequency) as detected in the 400–800 MHz octave by the FRB project on the Canadian Hydrogen Intensity Mapping Experiment (CHIME/FRB), using events from the first CHIME/FRB catalog. The catalog consists of 62 bursts from 18 repeating sources, plus 474 one-off FRBs, detected between 2018 July 25 and 2019 July 2. We identify four observed archetypes of burst morphology (“simple broadband,” “simple narrowband,” “temporally complex,” and “downward drifting”) and describe relevant instrumental biases that are essential for interpreting the observed morphologies. Using the catalog properties of the FRBs, we confirm that bursts from repeating sources, on average, have larger widths, and we show, for the first time, that bursts from repeating sources, on average, are narrower in bandwidth. This difference could be due to beaming or propagation effects, or it could be intrinsic to the populations. We discuss potential implications of these morphological differences for using…

Abstract The question of whether cosmic microwave background (CMB) temperature and polarization data from Planck favor a spatially closed universe with curvature parameter Ω K < 0 has been the subject of recent intense discussions. Attempts to break the geometrical degeneracy combining Planck data with external data sets such as baryon acoustic oscillation (BAO) measurements all point toward a spatially flat universe at the cost of significant tensions with Planck, which makes the resulting data set combination problematic. Settling this issue requires identifying a data set that can break the geometrical degeneracy while avoiding these tensions. We argue that cosmic chronometers (CCs), measurements of the expansion rate H ( z ) from the relative ages of massive early-type passively evolving galaxies, are the data set we are after. Furthermore, CCs come with the additional advantage of being virtually free of cosmological model assumptions. Combining Planck 2018 CMB temperature and polarization data with the latest CCs, we break the geometrical degeneracy and find Ω K = −0.0054 ±…

Abstract We present the Stage II results from the ongoing Satellites Around Galactic Analogs (SAGA) Survey. Upon completion, the SAGA Survey will spectroscopically identify satellite galaxies brighter than M r , o = −12.3 around 100 Milky Way (MW) analogs at z ∼ 0.01. In Stage II, we have more than quadrupled the sample size of Stage I, delivering results from 127 satellites around 36 MW analogs with an improved target selection strategy and deep photometric imaging catalogs from the Dark Energy Survey and the Legacy Surveys. We have obtained 25,372 galaxy redshifts, peaking around z = 0.2. These data significantly increase spectroscopic coverage for very low redshift objects in 17 < r o < 20.75 around SAGA hosts, creating a unique data set that places the Local Group in a wider context. The number of confirmed satellites per system ranges from zero to nine and correlates with host galaxy and brightest satellite luminosities. We find that the number and luminosities of MW satellites are consistent with being drawn from the same underlying distribution as SAGA systems. The majority…

Abstract Simulated images of a black hole surrounded by optically thin emission typically display two main features: a central brightness depression and a narrow “photon ring” consisting of strongly lensed images superimposed over the direct emission. The photon ring closely tracks a theoretical curve on the image plane corresponding to light rays that asymptote to bound photon orbits. The size and shape of this critical curve are purely governed by the Kerr geometry; in contrast, the size, shape, and depth of the observed brightness depression depend on the details of the emission region. For instance, images of spherical accretion models display a distinctive dark region—the “black hole shadow”—that completely fills the photon ring. By contrast, in models of equatorial disks extending to the event horizon, the darkest region in the image is restricted to a much smaller area—an inner shadow —whose edge lies near the direct lensed image of the equatorial horizon. Using both general relativistic MHD simulations and semi-analytic models, we demonstrate that the photon ring and inner…

Abstract This is the second paper of a series devoted to presenting an updated release of the BaSTI (a Bag of Stellar Tracks and Isochrones) stellar model and isochrone library. Following the publication of the updated solar-scaled library, here we present the library for an α -enhanced heavy element distribution. These new α -enhanced models account for all improvements and updates in the reference solar metal distribution and physics inputs, as in the new solar-scaled library. The models cover a mass range between 0.1 and 15 M ⊙ , 18 metallicities between [Fe/H] = −3.20 and +0.06 with [ α /Fe] = +0.4, and a He-to-metal enrichment ratio Δ Y /Δ Z = 1.31. For each metallicity, He-enhanced stellar models are also provided. The isochrones cover (typically) an age range between 20 Myr and 14.5 Gyr, including consistently the pre-main-sequence phase. The asteroseismic properties of the theoretical models have also been calculated. Models and isochrones have been compared with results from independent calculations, with the previous BaSTI release, and also with selected observations, to…

Abstract We present a measurement of the Hubble constant H 0 from surface brightness fluctuation (SBF) distances for 63 bright, mainly early-type galaxies out to 100 Mpc observed with the WFC3/IR on the Hubble Space Telescope (HST). The sample is drawn from several independent HST imaging programs using the F110W bandpass, with the majority of the galaxies being selected from the MASSIVE survey. The distances reach the Hubble flow with a median statistical uncertainty per measurement of 4%. We construct the Hubble diagram with these IR SBF distances and constrain H 0 using four different treatments of the galaxy velocities. For the SBF zero-point calibration, we use both the existing tie to Cepheid variables, updated for consistency with the latest determination of the distance to the Large Magellanic Cloud from detached eclipsing binaries, and a new tie to the tip of the red giant branch (TRGB) calibrated from the maser distance to NGC 4258. These two SBF calibrations are consistent with each other and with theoretical predictions from stellar population models. From a weighted…

Abstract The composition of giant planets is imprinted by their migration history and the compositional structure of their hosting disks. Studies in recent literature have investigated how the abundances of C and O can constrain the formation pathways of giant planets forming within few tens of au from a star. New ALMA observations, however, suggest planet-forming regions possibly extending to hundreds of au. We explore the implications of these wider formation environments through n -body simulations of growing and migrating giant planets embedded in planetesimal disks, coupled with a compositional model of the protoplanetary disk where volatiles are inherited from the molecular cloud and refractories are calibrated against extrasolar and Solar System data. We find that the C/O ratio provides limited insight on the formation pathways of giant planets that undergo large-scale migration. This limitation can be overcome, however, thanks to nitrogen and sulfur. Jointly using the C/N, N/O, and C/O ratios breaks any degeneracy in the formation and migration tracks of giant planets. The…

Abstract We present the third open gravitational-wave catalog (3-OGC) of compact-binary coalescences, based on the analysis of the public LIGO and Virgo data from 2015 through 2019 (O1, O2, O3a). Our updated catalog includes a population of 57 observations, including 4 binary black hole mergers that had not been previously reported. This consists of 55 binary black hole mergers and the 2 binary neutron star mergers, GW170817 and GW190425. We find no additional significant binary neutron star or neutron star–black hole merger events. The most confident new detection is the binary black hole merger GW190925_232845, which was observed by the LIGO–Hanford and Virgo observatories with  astro > 0.99 ; its primary and secondary component masses are 20.2 − 2.5 + 3.9 M ⊙ and 15.6 − 2.6 + 2.1 M ⊙ , respectively. We estimate the parameters of all binary black hole events using an up-to-date waveform model that includes both subdominant harmonics and precession effects. To enable deep follow up as our understanding of the underlying populations evolves, we make available our comprehensive…

Abstract We present the first results from the Mapping Obscuration to Reionization with ALMA (MORA) survey, the largest Atacama Large Millimeter/submillimeter Array (ALMA) blank-field contiguous survey to date (184 arcmin 2 ) and the only at 2 mm to search for dusty star-forming galaxies (DSFGs). We use the 13 sources detected above 5 σ to estimate the first ALMA galaxy number counts at this wavelength. These number counts are then combined with the state-of-the-art galaxy number counts at 1.2 and 3 mm and with a backward evolution model to place constraints on the evolution of the IR luminosity function and dust-obscured star formation in the past 13 billion years. Our results suggest a steep redshift evolution on the space density of DSFGs and confirm the flattening of the IR luminosity function at faint luminosities, with a slope of . We conclude that the dust-obscured component, which peaks at z ≈ 2–2.5, has dominated the cosmic history of star formation for the past ∼12 billion years, back to z ∼ 4. At z = 5, the dust-obscured star formation is estimated to be ∼35% of the total…

Abstract TauREx 3 is the next generation of the TauREx exoplanet atmospheric retrieval framework for Windows, Mac, and Linux. It is a complete rewrite with a full Python stack that makes it easy-to-use, high-performance, dynamic, and flexible. The new main TauREx program is built with modularity in mind, allowing the user to augment its functionalities with custom code and efficiently perform retrievals on custom parameters. We achieve this result by dynamic determination of fitting parameters, whereby TauREx 3 can detect new parameters for retrieval from user code through a simple interface. TauREx 3 can act as a library with a simple import taurex command, providing a rich set of classes and functions related to atmospheric modeling. A 10× speedup in forward model computations is achieved as compared to the previous version with a sixfold reduction in retrieval times while maintaining robust results. TauREx 3 is intended as a standalone, all-in-one package for retrievals while the TauREx 3 Python library can build or augment a user’s custom data pipeline easily.

Abstract The use of Type Ia supernovae (SNe Ia) as cosmological tools has motivated significant effort to understand what drives the intrinsic scatter of SN Ia distance modulus residuals after standardization, characterize the distribution of SN Ia colors, and explain why properties of the host galaxies of the SNe correlate with SN Ia distance modulus residuals. We use a compiled sample of ∼1450 spectroscopically confirmed photometric light curves of SNe Ia and propose a solution to these three problems simultaneously that also explains an empirical 11 σ detection of the dependence of Hubble residual scatter on SN Ia color. We introduce a physical model of color where intrinsic SN Ia colors with a relatively weak correlation with luminosity are combined with extrinsic dust-like colors ( E ( B − V )) with a wide range of extinction parameter values ( R V ). This model captures the observed trends of Hubble residual scatter and indicates that the dominant component of SN Ia intrinsic scatter is variation in R V . We also find that the recovered E ( B − V ) and R V distributions differ…

Abstract To date, modern three-dimensional (3D) supernova (SN) simulations have not demonstrated that explosion energies of 10 51 erg (=1 bethe=1 B) or more are possible for neutrino-driven SNe of non/slow-rotating M < 20 M ⊙ progenitors. We present the first such model, considering a nonrotating, solar-metallicity 18.88 M ⊙ progenitor, whose final 7 minutes of convective oxygen-shell burning were simulated in 3D and showed a violent oxygen–neon shell merger prior to collapse. A large set of 3D SN models was computed with the Prometheus-Vertex code, whose improved convergence of the two-moment equations with Boltzmann closure allows now to fully exploit the implicit neutrino-transport treatment. Nuclear burning is treated with a 23-species network. We vary the angular grid resolution and consider different nuclear equations of state and muon formation in the proto-neutron star (PNS), which requires six-species transport with coupling of all neutrino flavors across all energy–momentum groups. Elaborate neutrino transport was applied until ∼2 s after bounce. In one case, the…

Abstract Over the past few decades, the measurement precision of some pulsar timing experiments has advanced from ∼10 μ s to ∼10 ns, revealing many subtle phenomena. Such high precision demands both careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present PINT ( P INT I s N ot T empo3 ), a high-precision Python pulsar timing data analysis package, which is hosted on GitHub and available on the Python Package Index (PyPI) as pint-pulsar . PINT is well tested, validated, object oriented, and modular, enabling interactive data analysis and providing an extensible and flexible development platform for timing applications. It utilizes well-debugged public Python packages (e.g., the N um P y and A stropy libraries) and modern software development schemes (e.g., version control and efficient development with git and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility. PINT is developed and implemented without referring to, copying, or transcribing the code from other traditional…

Abstract We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range 6.3 < z ≤ 7.64, including 32 quasars at z > 6.5, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the Mg ii emission lines are in the range (0.3–3.6) × 109 M ⊙, which requires massive seed black holes with masses ≳103–104 M ⊙, assuming Eddington accretion since z = 30. The Eddington ratio distribution peaks at λ Edd ∼ 0.8 and has a mean of 1.08, suggesting high accretion rates for these quasars. The C iv–Mg ii emission-line velocity differences in our sample show an increase of C iv blueshift toward higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The Fe ii/Mg ii flux ratios derived for these quasars up to z = 7.6, compared with previous measurements at…

Abstract We present maps of the stellar streams detected in the Gaia Data Release 2 (DR2) and Early Data Release 3 (EDR3) catalogs using the STREAMFINDER algorithm. We also report the spectroscopic follow-up of the brighter DR2 stream members obtained with the high-resolution CFHT/ESPaDOnS and VLT/UVES spectrographs as well as with the medium-resolution NTT/EFOSC2 spectrograph. Two new stellar streams that do not have a clear progenitor are detected in DR2 (named Hríd and Gunnthrá), and seven are detected in EDR3 (named Gaia-6 to Gaia-12). Several candidate streams are also identified. The software also finds very long tidal tails associated with the 15 globular clusters: NGC 288, NGC 1261, NGC 1851, NGC 2298, NGC 2808, NGC 3201, M68, ω Cen, NGC 5466, Palomar 5, M5, NGC 6101, M92, NGC 6397, and NGC 7089. These stellar streams will be used in subsequent contributions in this series to chart the properties of the Galactic acceleration field on ∼100 pc to ∼100 kpc scales.

Abstract Spectroscopic detection of narrow emission lines traces the presence of circumstellar mass distributions around massive stars exploding as core-collapse supernovae. Transient emission lines disappearing shortly after the supernova explosion suggest that the material spatial extent is compact and implies an increased mass loss shortly prior to explosion. Here, we present a systematic survey for such transient emission lines (Flash Spectroscopy) among Type II supernovae detected in the first year of the Zwicky Transient Facility survey. We find that at least six out of ten events for which a spectrum was obtained within two days of the estimated explosion time show evidence for such transient flash lines. Our measured flash event fraction (>30% at 95% confidence level) indicates that elevated mass loss is a common process occurring in massive stars that are about to explode as supernovae.

Abstract We investigate a new class of habitable planets composed of water-rich interiors with massive oceans underlying H 2 -rich atmospheres, referred to here as Hycean worlds. With densities between those of rocky super-Earths and more extended mini-Neptunes, Hycean planets can be optimal candidates in the search for exoplanetary habitability and may be abundant in the exoplanet population. We investigate the bulk properties (masses, radii, and temperatures), potential for habitability, and observable biosignatures of Hycean planets. We show that Hycean planets can be significantly larger compared to previous considerations for habitable planets, with radii as large as 2.6 R ⊕ (2.3 R ⊕ ) for a mass of 10 M ⊕ (5 M ⊕ ). We construct the Hycean habitable zone (HZ), considering stellar hosts from late M to Sun-like stars, and find it to be significantly wider than the terrestrial-like HZ. While the inner boundary of the Hycean HZ corresponds to equilibrium temperatures as high as ∼500 K for late M dwarfs, the outer boundary is unrestricted to arbitrarily large orbital separations.…

Abstract We present new stellar mass functions at z ∼ 6, z ∼ 7, z ∼ 8, z ∼ 9 and, for the first time, z ∼ 10, constructed from ∼800 Lyman-break galaxies previously identified over the eXtreme Deep Field and Hubble Ultra-Deep Field parallel fields and the five Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields. Our study is distinctive due to (1) the much deeper (∼200 hr) wide-area Spitzer/Infrared Array Camera (IRAC) imaging at 3.6 μ m and 4.5 μ m from the Great Observatories Origins Deep Survey Re-ionization Era Wide-area Treasury from Spitzer program (GREATS) and (2) consideration of z ∼ 6–10 sources over a 3× larger area than those of previous Hubble Space Telescope+Spitzer studies. The Spitzer/IRAC data enable ≥2 σ rest-frame optical detections for an unprecedented 50% of galaxies down to a stellar mass limit of ∼ 10 8  ⊙ across all redshifts. Schechter fits to our volume densities suggest a combined evolution in the characteristic mass  * and normalization factor ϕ * between z ∼ 6 and z ∼ 8. The stellar mass density (SMD) increases by ∼1000× in the ∼500 Myr…

Abstract Several lines of evidence suggest that the Milky Way underwent a major merger at z ∼ 2 with the Gaia-Sausage-Enceladus (GSE) galaxy. Here we use H3 Survey data to argue that GSE entered the Galaxy on a retrograde orbit based on a population of highly retrograde stars with chemistry similar to the largely radial GSE debris. We present the first tailored N-body simulations of the merger. From a grid of ≈500 simulations we find that a GSE with M ⋆ = 5 × 108 M ⊙, M DM = 2 × 1011 M ⊙ best matches the H3 data. This simulation shows that the retrograde stars are stripped from GSE’s outer disk early in the merger. Despite being selected purely on angular momenta and radial distributions, this simulation reproduces and explains the following phenomena: (i) the triaxial shape of the inner halo, whose major axis is at ≈35° to the plane and connects GSE’s apocenters; (ii) the Hercules-Aquila Cloud and the Virgo Overdensity, which arise due to apocenter pileup; and (iii) the 2 Gyr lag between the quenching of GSE and the truncation of the age distribution of the in situ halo, which…

Abstract The streaming instability (SI) is a mechanism to aerodynamically concentrate solids in protoplanetary disks and trigger the formation of planetesimals. The SI produces strong particle clumping if the ratio of solid to gas surface density—an effective metallicity—exceeds a critical value. This critical value depends on particle sizes and disk conditions such as radial drift-inducing pressure gradients and levels of turbulence. To quantify these thresholds, we perform a suite of vertically stratified SI simulations over a range of dust sizes and metallicities. We find a critical metallicity as low as 0.4% for the optimum particle sizes and standard radial pressure gradients (normalized value of Π = 0.05). This subsolar metallicity is lower than previous results, due to improved numerical methods and computational effort. We discover a sharp increase in the critical metallicity for small solids, when the dimensionless stopping time (Stokes number) is ≤0.01. We provide simple fits to the size-dependent SI clumping threshold, including generalizations to different disk models…

Abstract We present a systematic numerical relativity study of the dynamical ejecta, winds, and nucleosynthesis in neutron star (NS) merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOSs) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamic turbulence up to 100 ms postmerger. Spiral density waves propagating from the NS remnant to the disk trigger a wind with mass flux ∼0.1–0.5 M ⊙ s −1 , which persists for the entire simulation as long as the remnant does not collapse to a black hole. This wind has average electron fraction ≳0.3 and average velocity ∼0.1–0.17 c and thus is a site for the production of weak r -process elements (mass number A < 195). Disks around long-lived remnants have masses ∼0.1–0.2 M ⊙ , temperatures peaking at ≲10 MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not…

Abstract The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf Galaxy (Sgr), Fornax (Fnx), and the now fully disrupted Gaia Sausage/Enceladus (GSE) system. We present and analyze the APOGEE chemical abundance patterns of each galaxy to draw robust conclusions about their star formation histories, by quantifying the relative abundance trends of multiple elements (C, N, O, Mg, Al, Si, Ca, Fe, Ni, and Ce), as well as by fitting chemical evolution models to the [ α /Fe]–[Fe/H] abundance plane for each galaxy. Results show that the chemical signatures of the starburst in the Magellanic Clouds (MCs) observed by Nidever et al. in the α -element abundances extend to C+N, Al, and Ni, with the major burst in the SMC occurring some 3–4 Gyr before the burst in the LMC. We find that Sgr and Fnx also exhibit chemical abundance patterns suggestive of secondary…

Abstract A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models—including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is 2.5× larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range 2000–11,000 Å (1800 Å redder) and reduced uncertainties compared to SALT2, enabling accurate use of low-z I and iz photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low-z Foundation…

Abstract We present results from the “Mint” resolution DC Justice League suite of Milky Way–like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultrafaint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to z = 0, with initial star (dark matter) particle masses of 994 (17900) M ⊙, and a force resolution of 87 pc. We study the surrounding dwarfs and UFDs, and find that the simulations match the observed dynamical properties of galaxies with −3 > M V > −19, and reproduce the scatter seen in the size–luminosity plane for r h ≳ 200 pc. We predict the vast majority of nearby galaxies will be observable by the Vera Rubin Observatory’s coadded Legacy Survey of Space and Time. We additionally show that faint dwarfs with velocity dispersions ≲5 km s−1 result from severe tidal stripping of the host halo. We investigate the quenching of UFDs in a hydrodynamical Milky Way context and find that the majority of UFDs are quenched prior to…

Abstract We present an update of the open-source photochemical kinetics code VULCAN to include C–H–N–O–S networks and photochemistry. The additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al., Venot et al., and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b,…

Abstract The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGNs) disks, where binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Here we compare the properties of these observed events to the theoretical models of mergers in AGN disks, which are obtained by performing one-dimensional N -body simulations combined with semi-analytical prescriptions. The high BH masses in GW190521 are consistent with mergers of high-generation (high-g) BHs where the initial progenitor stars had high metallicity, 2g BHs if the original progenitors were metal-poor, or 1g BHs that had gained mass via super-Eddington accretion. Other measured properties related to spin parameters in GW190521 are also consistent with mergers in AGN disks. Furthermore, mergers in the lower mass gap or those…

Abstract Using physical models, we study the sensitivity of polycyclic aromatic hydrocarbon (PAH) emission spectra to the character of the illuminating starlight, to the PAH size distribution, and to the PAH charge distribution. The starlight models considered range from the emission from a 3 Myr old starburst, rich in far-ultraviolet (FUV) radiation, to the FUV-poor spectrum of the very old population of the M31 bulge. A wide range of starlight intensities is considered. The effects of reddening in dusty clouds are investigated for different starlight spectra. For a fixed PAH abundance parameter q PAH (the fraction of the total dust mass in PAHs with <10 3 C atoms), the fraction of the infrared power appearing in the PAH emission features can vary by a factor of two as the starlight spectrum varies from FUV-poor (M31 bulge) to FUV-rich (young starburst). We show how q PAH can be measured from the strength of the 7.7 μ m emission. The fractional power in the 17 μ m feature can be suppressed by high starlight intensities.

Abstract High-resolution UV spectra of stellar H i Lyα lines from the Hubble Space Telescope (HST) provide observational constraints on the winds of coronal main-sequence stars, thanks to an astrospheric absorption signature created by the interaction between the stellar winds and the interstellar medium. We report the results of a new HST survey of M dwarf stars, yielding six new detections of astrospheric absorption. We estimate mass-loss rates for these detections and upper limits for nondetections. These new constraints allow us to characterize the nature of M dwarf winds and their dependence on coronal activity for the first time. For a clear majority of the M dwarfs, we find winds that are weaker than or comparable in strength to that of the Sun, i.e., M ̇ ≤ 1 M ̇ ⊙ . However, two of the M dwarfs have much stronger winds: YZ CMi (M4 Ve; M ̇ = 30 M ̇ ⊙ ) and GJ 15AB (M2 V+M3.5 V; M ̇ = 10 M ̇ ⊙ ). Even these winds are much weaker than expectations if the solar relation between flare energy and coronal mass ejection (CME) mass extended to M dwarfs. Thus, the solar flare/CME…

Abstract Active galactic nuclei (AGNs) are powered by geometrically thin accretion disks surrounding a central supermassive black hole. Here we explore the evolution of stars embedded in these extreme astrophysical environments (AGN stars). Because AGN disks are much hotter and denser than most components of the interstellar medium, AGN stars are subject to very different boundary conditions than normal stars. They are also strongly affected by both mass accretion, which can run away given the vast mass of the disk, and mass loss due to super-Eddington winds. Moreover, chemical mixing plays a critical role in the evolution of these stars by allowing fresh hydrogen accreted from the disk to mix into their cores. We find that, depending on the local AGN density and sound speed and the duration of the AGN phase, AGN stars can rapidly become very massive (M > 100 M ⊙). These stars undergo core collapse, leave behind compact remnants, and contribute to polluting the disk with heavy elements. We show that the evolution of AGN stars can have a profound impact on the evolution of AGN…

Abstract To date, about two dozen low-mass embedded protostars exhibit rich spectra with lines of complex organic molecules (COMs). These protostars seem to possess a different enrichment in COMs. However, the statistics of COM abundance in low-mass protostars are limited by the scarcity of observations. This study introduces the Perseus ALMA Chemistry Survey (PEACHES), which aims at unbiasedly characterizing the chemistry of COMs toward the embedded (Class 0/I) protostars in the Perseus molecular cloud. Of the 50 embedded protostars surveyed, 58% of them have emission from COMs. 56%, 32%, and 40% of the protostars have CH3OH, CH3OCHO, and N-bearing COMs, respectively. The detectability of COMs depends neither on the averaged continuum brightness temperature, a proxy of the H2 column density, nor on the bolometric luminosity and the bolometric temperature. For the protostars with detected COMs, CH3OH has a tight correlation with CH3CN, spanning more than two orders of magnitude in column densities normalized by the continuum brightness temperature, suggesting a chemical relation…

Abstract Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), M hm , using the semianalytic model Galacticus to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields M hm < 10 7.0 M ⊙ (WDM particle mass m WDM > 9.7 keV) at 95% confidence and disfavors M hm = 10 7.4 M ⊙ ( m WDM = 7.4 keV) with a 20:1 marginal likelihood…

Abstract The gas content of the complete compilation of Local Group dwarf galaxies (119 within 2 Mpc) is presented using H i survey data. Within the virial radius of the Milky Way (224 kpc here), 53 of 55 dwarf galaxies are devoid of gas to limits of M H i < 104 M ⊙. Within the virial radius of M31 (266 kpc), 27 of 30 dwarf galaxies are devoid of gas (with limits typically <105 M ⊙). Beyond the virial radii of the Milky Way and M31, the majority of the dwarf galaxies have detected H i gas and H i masses higher than the limits. When the relationship between gas content and distance is investigated using a Local Group virial radius, more of the nondetected dwarf galaxies are within this radius (85 ± 1 of the 93 nondetected dwarf galaxies) than within the virial radii of the Milky Way and M31. Using the Gaia proper-motion measurements available for 38 dwarf galaxies, the minimum gas density required to completely strip them of gas is calculated. Halo densities between 10−5 and 5 × 10−4 cm−3 are typically required for instantaneous stripping at perigalacticon. When compared to halo…

Abstract We present the systematic spectral analyses of gamma-ray bursts (GRBs) detected by the Fermi Gamma-Ray Burst Monitor during its first ten years of operation. This catalog contains two types of spectra: time-integrated spectral fits and spectral fits at the brightest time bin, from 2297 GRBs, resulting in a compendium of over 18,000 spectra. The four different spectral models used for fitting the spectra were selected based on their empirical importance to the shape of many GRBs. We describe in detail our procedure and criteria for the analyses, and present the bulk results in the form of parameter distributions both in the observer frame and in the GRB rest frame. 941 GRBs from the first four years have been refitted using the same methodology as that of the 1356 GRBs in years five through ten. The data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center.

Abstract Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, nonlinear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in the MESA stellar evolution software…

Abstract We present the discovery that ASASSN-14ko is a periodically flaring active galactic nucleus at the center of the galaxy ESO 253-G003. At the time of its discovery by the All-Sky Automated Survey for Supernovae (ASAS-SN), it was classified as a supernova close to the nucleus. The subsequent 6 yr of V - and g -band ASAS-SN observations revealed that ASASSN-14ko has nuclear flares occurring at regular intervals. The 17 observed outbursts show evidence of a decreasing period over time, with a mean period of P 0 = 114.2 ± 0.4 days and a period derivative of . The most recent outburst in 2020 May, which took place as predicted, exhibited spectroscopic changes during the rise and had a UV bright, blackbody spectral energy distribution similar to tidal disruption events (TDEs). The X-ray flux decreased by a factor of 4 at the beginning of the outburst and then returned to its quiescent flux after ∼8 days. The Transiting Exoplanet Survey Satellite observed an outburst during Sectors 4–6, revealing a rise time of 5.60 ± 0.05 days in the optical and a decline that is best fit with an…

Abstract We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0 < z < 5 and how the formation of this virialized phase affects the evolution of galactic disks. We demonstrate that when the halo mass crosses ∼10 12 M ⊙ , the cooling time of shocked gas in the inner CGM (∼0.1 R vir , where R vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average subvirial temperatures ( T ≪ T vir ), large pressure fluctuations, and supersonic inflow/outflow velocities to virial temperatures ( T ∼ T vir ), uniform pressures, and subsonic velocities. This transition occurs when the outer CGM (∼0.5 R vir ) is already subsonic and has a temperature ∼ T vir , indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We…

Abstract We compile observations of molecular gas contents and infrared-based star formation rates (SFRs) for 112 circumnuclear star-forming regions, in order to reinvestigate the form of the disk-averaged Schmidt surface density star-formation law in starbursts. We then combine these results with total gas and SFR surface densities for 153 nearby nonstarbursting disk galaxies from de los Reyes & Kennicutt (2019), to investigate the properties of the combined star formation law, following Kennicutt (1998). We confirm that the combined Schmidt law can be fitted with a single power law with slope n = 1.5 ± 0.05 (including fitting method uncertainties), somewhat steeper than the value n = 1.4 ± 0.15 found by Kennicutt. Fitting separate power laws to the nonstarbursting and starburst galaxies, however, produces very different slopes (n = 1.34 ± 0.07 and 0.98 ± 0.07, respectively), with a pronounced offset in the zero-point (∼0.6 dex) of the starburst relation to higher SFR surface densities. This offset is seen even when a common conversion factor between CO intensity and molecular…

Abstract We explore the kinematics of 27 z ≳ 6 quasar host galaxies observed in [C ii] 158 μm ([C ii]) emission with the Atacama Large Millimeter/submillimeter Array at a resolution of ≈025. We find that nine of the galaxies show disturbed [C ii] emission, due to either a close companion galaxy or a recent merger. Ten galaxies have smooth velocity gradients consistent with the emission arising from a gaseous disk. The remaining eight quasar host galaxies show no velocity gradient, suggesting that the gas in these systems is dispersion dominated. All galaxies show high velocity dispersions with a mean of 129 ± 10 km s−1. To provide an estimate of the dynamical mass within twice the half-light radius of the quasar host galaxy, we model the kinematics of the [C ii] emission line using our publicly available kinematic fitting code, qubefit. This results in a mean dynamical mass of 5.0 ± 0.8( ± 3.5) × 1010 M ⊙. Comparison between the dynamical mass and the mass of the supermassive black hole reveals that the sample falls above the locally derived bulge mass–black hole mass relation at…

Abstract Interaction-powered supernovae (SNe) explode within an optically thick circumstellar medium (CSM) that could be ejected during eruptive events. To identify and characterize such pre-explosion outbursts, we produce forced-photometry light curves for 196 interacting SNe, mostly of Type IIn, detected by the Zwicky Transient Facility between early 2018 and 2020 June. Extensive tests demonstrate that we only expect a few false detections among the 70,000 analyzed pre-explosion images after applying quality cuts and bias corrections. We detect precursor eruptions prior to 18 Type IIn SNe and prior to the Type Ibn SN 2019uo. Precursors become brighter and more frequent in the last months before the SN and month-long outbursts brighter than magnitude −13 occur prior to 25% (5–69%, 95% confidence range) of all Type IIn SNe within the final three months before the explosion. With radiative energies of up to 1049 erg, precursors could eject ∼1 M ⊙ of material. Nevertheless, SNe with detected precursors are not significantly more luminous than other SNe IIn, and the characteristic…

Abstract Massive and water-rich planets should be ubiquitous in the universe. Many of these worlds are expected to be subject to important irradiation from their host star, and display supercritical water layers surrounded by extended steam atmospheres. Irradiated ocean planets with such inflated hydrospheres have been recently shown to be good candidates for matching the mass–radius distribution of sub-Neptunes. Here we describe a model that computes a realistic structure for water-rich planets by combining an interior model with an updated equation of state (EOS) for water, and an atmospheric model that takes into account radiative transfer. We find that the use of inappropriate EOSs can lead to the overestimation of the planetary radius by up to ∼10%, depending on the planet size and composition. Our model has been applied to the GJ 9827 system as a test case and indicates Earth- or Venus-like interiors for planets b and c, respectively. Planet d could be an irradiated ocean planet with a water mass fraction (WMF) of ∼20% ± 10%. We also provide fits for the mass–radius…

Abstract We present a synthesis of the astronomical observations constraining the wavelength-dependent extinction, emission, and polarization from interstellar dust from UV to microwave wavelengths on diffuse Galactic sight lines. Representative solid-phase abundances for those sight lines are also derived. Given the sensitive new observations of polarized dust emission provided by the Planck satellite, we place particular emphasis on dust polarimetry, including continuum polarized extinction, polarization in the carbonaceous and silicate spectroscopic features, the wavelength-dependent polarization fraction of the dust emission, and the connection between optical polarized extinction and far-infrared polarized emission. Together, these constitute a set of constraints that should be reproduced by models of dust in the diffuse interstellar medium.

Abstract Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud life cycle and star formation efficiency (SFE) remains an open question. We perform radiation magnetohydrodynamic simulations of star-forming giant molecular clouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter (1 ≤ α vir,0 ≤ 5) or dimensionless mass-to-magnetic flux ratio (0.5 ≤ μ Φ,0 ≤ 8 and ∞ ); the initial mass 10 5 M ⊙ and radius 20 pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4–8 Myr. Both the final SFE ( ε * ) and time-averaged SFE per freefall time ( ε ff ) are reduced by strong turbulence and magnetic fields. The median ε * ranges between 2.1% and 9.5%. The median ε ff ranges between 1.0% and 8.0%, and anticorrelates with α vir,0 , in…

Abstract Motivated by the recent discovery of a compact object with mass in the range 2.5–2.67 M ⊙ in the binary merger GW190814, we revisit the question of the maximum mass of neutron stars (NSs). We use a Markov Chain Monte Carlo approach to generate about two million phenomenological equations of state with and without first-order phase transitions. We fix the crust equation of state and only assume causality at higher densities. We show how a strict upper bound on the maximum NS mass can be inferred from upcoming observation of NS radii and masses. The derived upper bounds depend only on relativity and causality, so it is not affected by nuclear physics uncertainties. We show how a lower limit on the maximum mass of NSs, in combination with upcoming measurements of NS radii by LIGO/Virgo and NICER, would constrain the equation of state of dense matter. Finally, we discuss the implications for GW190814.

Abstract Indications of disequilibrium throughout the Milky Way (MW) highlight the need for compact, flexible, non-parametric descriptions of phase–space distributions of galaxies. We present a new representation of the current dark matter (DM) distribution and potential derived from N -body simulations of the MW and Large Magellanic Cloud (LMC) system using basis function expansions (BFEs). We incorporate methods to maximize the physical signal in the representation. As a result, the simulations of 10 8 DM particles representing the distorted MW(MW+LMC) system can be described by ∼236(2067) coefficients. We find that the LMC induces asymmetric perturbations (odd l, m ) to the MW’s halo, which are inconsistent with oblate, prolate, or triaxial halos. Furthermore, the energy in high order even modes ( l , m > 2) is similar to average triaxial halos found in cosmological simulations. As such, the response of the MW’s halo to the LMC must be accounted for in order to recover the imprints of its assembly history. The LMC causes the outer halo (>30 kpc) to shift from the disk center of…

Abstract We report new measurements of millimeter-wave power spectra in the angular multipole range 2000 ≤ ℓ ≤ 11,000 (angular scales ). By adding 95 and 150 GHz data from the low-noise 500 deg2 SPTpol survey to the SPT-SZ three-frequency 2540 deg2 survey, we substantially reduce the uncertainties in these bands. These power spectra include contributions from the primary cosmic microwave background, cosmic infrared background, radio galaxies, and thermal and kinematic Sunyaev–Zel’dovich (SZ) effects. The data favor a thermal SZ (tSZ) power at 143 GHz of and a kinematic SZ (kSZ) power of . This is the first measurement of kSZ power at ≥3σ. However, different assumptions about the CIB or SZ models can reduce the significance down to 2.4σ in the worst case. We study the implications of the measured kSZ power for the epoch of reionization under the Calabrese et al. model for the kSZ power spectrum and find the duration of reionization to be ( at 95% confidence), when combined with our previously published tSZ bispectrum measurement. The upper limit tightens to if the assumed homogeneous…

Abstract We present a novel method to detect variable astrophysical objects and transient phenomena using anomalous excess scatter in repeated measurements from public catalogs of Gaia DR2 and Zwicky Transient Facility (ZTF) DR3 photometry. We first provide a generalized, all-sky proxy for variability using only Gaia DR2 photometry, calibrated to white dwarf stars. To ensure more robust candidate detection, we further employ a method combining Gaia with ZTF photometry and alerts. To demonstrate its efficacy, we apply this latter technique to a sample of roughly 12,100 white dwarfs within 200 pc centered on the ZZ Ceti instability strip, where hydrogen-atmosphere white dwarfs are known to pulsate. By inspecting the top 1% of the samples ranked by these methods, we demonstrate that both the Gaia-only and ZTF-informed techniques are highly effective at identifying known and new variable white dwarfs, which we verify using follow-up, high-speed photometry. We confirm variability in all 33 out of 33 (100%) observed white dwarfs within our top 1% highest-ranked candidates, both inside and…

Abstract One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed switchbacks. These δ B R / B ∼  ( 1 ) fluctuations occur over a range of timescales and in patches separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma β and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small (∼1°) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to ∼85 keV that we argue are the…

Abstract We present the first results from the ongoing, intensive, multiwavelength monitoring program of the luminous Seyfert 1 galaxy Mrk 817. While this active galactic nucleus was, in part, selected for its historically unobscured nature, we discovered that the X-ray spectrum is highly absorbed, and there are new blueshifted, broad, and narrow UV absorption lines, which suggest that a dust-free, ionized obscurer located at the inner broad-line region partially covers the central source. Despite the obscuration, we measure UV and optical continuum reverberation lags consistent with a centrally illuminated Shakura–Sunyaev thin accretion disk, and measure reverberation lags associated with the optical broad-line region, as expected. However, in the first 55 days of the campaign, when the obscuration was becoming most extreme, we observe a de-coupling of the UV continuum and the UV broad emission-line variability. The correlation recovered in the next 42 days of the campaign, as Mrk 817 entered a less obscured state. The short C iv and Lyα lags suggest that the accretion disk extends…

Abstract Stellar streams are excellent probes of the underlying gravitational potential in which they evolve. In this work, we fit dynamical models to five streams in the Southern Galactic hemisphere, combining observations from the Southern Stellar Stream Spectroscopic Survey ( S 5 ), Gaia EDR3, and the Dark Energy Survey, to measure the mass of the Large Magellanic Cloud (LMC). With an ensemble of streams, we find a mass of the LMC ranging from ∼14–19 × 10 10 M ⊙ , probed over a range of closest approach times and distances. With the most constraining stream (Orphan–Chenab), we measure an LMC mass of 18.8 − 4.0 + 3.5 × 10 10 M ⊙ , probed at a closest approach time of 310 Myr and a closest approach distance of 25.4 kpc. This mass is compatible with previous measurements, showing that a consistent picture is emerging of the LMC’s influence on structures in the Milky Way. Using this sample of streams, we find that the LMC’s effect depends on the relative orientation of the stream and LMC at their point of closest approach. To better understand this, we present a simple model based on…

Abstract We report the latest statistical analyses of superflares on solar-type (G-type main-sequence; effective temperature is 5100–6000 K) stars using all of the Kepler primary mission data and Gaia Data Release 2 catalog. We updated the flare detection method from our previous studies by using a high-pass filter to remove rotational variations caused by starspots. We also examined the sample biases on the frequency of superflares, taking into account gyrochronology and flare detection completeness. The sample sizes of solar-type and Sun-like stars (effective temperature is 5600–6000 K and rotation period is over 20 days in solar-type stars) are ∼4 and ∼12 times, respectively, compared with Notsu et al. As a result, we found 2341 superflares on 265 solar-type stars and 26 superflares on 15 Sun-like stars; the former increased from 527 to 2341 and the latter from three to 26 events compared with our previous study. This enabled us to have a more well-established view on the statistical properties of superflares. The observed upper limit of the flare energy decreases as the rotation…

Abstract We describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Lyα emitting galaxies between 1.88 < z < 3.52, in a 540 deg2 area encompassing a comoving volume of 10.9 Gpc3. No preselection of targets is involved; instead the HETDEX measurements are accomplished via a spectroscopic survey using a suite of wide-field integral field units distributed over the focal plane of the telescope. This survey measures the Hubble expansion parameter and angular diameter distance, with a final expected accuracy of better than 1%. We detail the project’s observational strategy, reduction pipeline, source detection, and catalog generation, and present initial results for science verification in the Cosmological Evolution Survey, Extended Groth Strip, and Great Observatories Origins Deep Survey North fields. We demonstrate that our data reach the required specifications in throughput, astrometric accuracy, flux limit, and object…

Abstract We have analyzed the parsec-scale jet kinematics of 447 bright radio-loud active active galactic nuclei (AGN), based on 15 GHz Very Long Baseline Array (VLBA) data obtained between 1994 August 31 and 2019 August 4. We present new total intensity and linear polarization maps obtained between 2017 January 1 and 2019 August 4 for 143 of these AGN. We tracked 1923 bright features for five or more epochs in 419 jets. The majority (60%) of the well-sampled jet features show either accelerated or nonradial motion. In 47 jets there is at least one nonaccelerating feature with an unusually slow apparent speed. Most of the jets show variations of 10°–50° in their inner jet position angle (PA) over time, although the overall distribution has a continuous tail out to 200°. AGN with spectral energy distributions peaked at lower frequencies tend to have more variable PAs, with BL Lac objects being less variable than quasars. The Fermi Large Area Telescope (LAT) gamma-ray-associated AGN also tend to have more variable PAs than the non-LAT AGN in our sample. We attribute these trends to…

Abstract Young stellar associations hold a star formation record that can persist for millions of years, revealing the progression of star formation long after the dispersal of the natal cloud. To identify nearby young stellar populations that trace this progression, we have designed a comprehensive framework for the identification of young stars and use it to identify ∼3 × 104 candidate young stars within a distance of 333 pc using Gaia DR2. Applying the HDBSCAN clustering algorithm to this sample, we identify 27 top-level groups, nearly half of which have little to no presence in previous literature. Ten of these groups have visible substructure, including notable young associations such as Orion, Perseus, Taurus, and Sco-Cen. We provide a complete subclustering analysis of all groups with substructure, using age estimates to reveal each region’s star formation history. The patterns we reveal include an apparent star formation origin for Sco-Cen along a semicircular arc, as well as clear evidence for sequential star formation moving away from that arc with a propagation speed of…

Abstract Winds from massive stars have velocities of 1000 km s −1 or more and produce hot, high-pressure gas when they shock. We develop a theory for the evolution of bubbles driven by the collective winds from star clusters early in their lifetimes, which involves interaction with the turbulent, dense interstellar medium of the surrounding natal molecular cloud. A key feature is the fractal nature of the hot bubble’s surface. The large area of this interface with surrounding denser gas strongly enhances energy losses from the hot interior, enabled by turbulent mixing and subsequent cooling at temperatures T ∼ 10 4 –10 5 K, where radiation is maximally efficient. Due to the extreme cooling, the bubble radius scales differently ( ) from the classical Weaver et al. solution and has expansion velocity and momentum lower by factors of 10–10 2 at given , with pressure lower by factors of 10 2 –10 3 . Our theory explains the weak X-ray emission and low shell expansion velocities of observed sources. We discuss further implications of our theory for observations of the hot bubbles and…

Abstract Time-domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled griz light curves for thousands of transient events up to z ≈ 0.2. This large sample of transients with four-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as…

Abstract Recently we have shown that high-energy neutrinos above 200 TeV detected by IceCube are produced within several parsecs in the central regions of radio-bright blazars, that is active galactic nuclei with jets pointing toward us. To independently test this result and extend the analysis to a wider energy range, we use public data for all neutrino energies from seven years of IceCube observations. The IceCube point-source likelihood map is analyzed against the positions of blazars from a statistically complete sample selected according to their compact radio flux density. The latter analysis delivers a 3.0 σ significance, with the combined post-trial significance of both studies being 4.1 σ . The correlation is driven by a large number of blazars. Together with fainter but physically similar sources not included in the sample, they may explain the entire IceCube astrophysical neutrino flux as derived from muon-track analyses. The neutrinos can be produced in interactions of relativistic protons with X-ray self-Compton photons in parsec-scale blazar jets.

Abstract A primary goal of exoplanet characterization is to use a planet’s current composition to understand how that planet formed. For example, the C/O ratio has long been recognized as carrying important information on the chemistry of volatile species. Refractory elements, like Fe, Mg, and Si, are usually not considered in this conversation because they condense into solids like Fe(s) or MgSiO 3 and would be removed from the observable, gaseous atmosphere in exoplanets cooler than about 2000 K. However, planets hotter than about 2000 K, called ultra-hot Jupiters (UHJs), are warm enough to largely avoid the condensation of refractory species. In this paper, we explore the insight that the measurement of refractory abundances can provide into a planet’s origins. Through refractory-to-volatile elemental abundance ratios, we can estimate a planet’s atmospheric rock-to-ice fraction and constrain planet formation and migration scenarios. We first relate a planet’s present-day refractory-to-volatile ratio to its rock-to-ice ratio from formation using various compositional models for…

Abstract Monitoring of vibrational eigenmodes of an elastic body excited by gravitational waves was one of the first concepts proposed for the detection of gravitational waves. At laboratory scale, these experiments became known as resonant bar detectors first developed by Joseph Weber in the 1960s. Due to the dimensions of these bars, the targeted signal frequencies were in the kHz range. Weber also pointed out that monitoring of vibrations of Earth or the Moon could reveal gravitational waves in the mHz band. His Lunar Surface Gravimeter experiment deployed on the Moon by the Apollo 17 crew had a technical failure, which greatly reduced the science scope of the experiment. In this article, we revisit the idea and propose a Lunar Gravitational-Wave Antenna (LGWA). We find that LGWA could become an important partner observatory for joint observations with the space-borne, laser-interferometric detector LISA and at the same time contribute an independent science case due to LGWA’s unique features. Technical challenges need to be overcome for the deployment of the experiment, and…

Abstract Rapid binary population synthesis codes are often used to investigate the evolution of compact-object binaries. They typically rely on analytical fits of single-star evolutionary tracks and parameterized models for interactive phases of evolution (e.g., mass transfer on a thermal timescale, determination of dynamical instability, and common envelope) that are crucial to predict the fate of binaries. These processes can be more carefully implemented in stellar structure and evolution codes such as MESA. To assess the impact of such improvements, we compare binary black hole mergers as predicted in models with the rapid binary population synthesis code COSMIC to models ran with MESA simulations through mass transfer and common-envelope treatment. We find that results significantly differ in terms of formation paths, the orbital periods and mass ratios of merging binary black holes, and consequently merger rates. While common-envelope evolution is the dominant formation channel in COSMIC, stable mass transfer dominates in our MESA models. Depending upon the black hole donor…

Abstract Merging neutron stars produce “kilonovae”—electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their…

Abstract The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide…

Abstract We measured the mid-infrared (MIR) extinction using Spitzer photometry and spectroscopy (3.6–37 μ m) for a sample of Milky Way sight lines (mostly) having measured ultraviolet extinction curves. We used the pair method to determine the MIR extinction that we then fit with a power law for the continuum and modified Drude profiles for the silicate features. We derived 16 extinction curves having a range of A ( V ) (1.8–5.5) and R ( V ) values (2.4–4.3). Our sample includes two dense sight lines that have 3 μ m ice feature detections and weak 2175 Å bumps. The average A ( λ )/ A ( V ) diffuse sight-line extinction curve we calculate is lower than most previous literature measurements. This agrees better with literature diffuse dust grain models, though it is somewhat higher. The 10 μ m silicate feature does not correlate with the 2175 Å bump, for the first time providing direct observational confirmation that these two features arise from different grain populations. The strength of the 10 μ m silicate feature varies by ∼2.5 and is not correlated with A ( V ) or R ( V ). It is…

Abstract We examine Lyman continuum (LyC) leakage through H ii regions regulated by turbulence and radiative feedback in a giant molecular cloud in the context of fully coupled radiation hydrodynamics (RHD). The physical relations of the LyC escape with H i covering fraction, kinematics, ionizing photon production efficiency, and emergent Lyα line profiles are studied using a series of RHD turbulence simulations performed with ramses-rt. The turbulence-regulated mechanism allows ionizing photons to leak out at early times before the onset of supernova feedback. The LyC photons escape through turbulence-generated low column density channels that are evacuated efficiently by radiative feedback via photoheating-induced shocks across the D-type ionization fronts. The Lyα photons funnel through the photoionized channels along the paths of LyC escape, resulting in a diverse Lyα spectral morphology including narrow double-peaked profiles. The Lyα peak separation is controlled by the residual H i column density of the channels, and the line asymmetry correlates with the porosity and…

Abstract The discovery of periodicity in the arrival times of the fast radio bursts (FRBs) poses a challenge to the oft-studied magnetar scenarios. However, models that postulate that FRBs result from magnetized shocks or magnetic reconnection in a relativistic outflow are not specific to magnetar engines; instead, they require only the impulsive injection of relativistic energy into a dense magnetized medium. Motivated thus, we outline a new scenario in which FRBs are powered by short-lived relativistic outflows (“flares”) from accreting black holes or neutron stars, which propagate into the cavity of the pre-existing (“quiescent”) jet. In order to reproduce FRB luminosities and rates, we are driven to consider binaries of stellar-mass compact objects undergoing super-Eddington mass transfer, similar to ultraluminous X-ray (ULX) sources. Indeed, the host galaxies of FRBs, and their spatial offsets within their hosts, show broad similarities with ULXs. Periodicity on timescales of days to years could be attributed to precession (e.g., Lens-Thirring) of the polar accretion funnel,…

Abstract Stellar population models produce radiation fields that ionize oxygen up to O +2 , defining the limit of standard H ii region models (<54.9 eV). Yet, some extreme emission-line galaxies, or EELGs, have surprisingly strong emission originating from much higher ionization potentials. We present UV HST/COS and optical LBT/MODS spectra of two nearby EELGs that have very high-ionization emission lines (e.g., He ii λλ 1640,4686 C iv λλ 1548,1550, [Fe v ] λ 4227, [Ar iv ] λλ 4711,4740). We define a four-zone ionization model that is augmented by a very high-ionization zone, as characterized by He +2 (>54.4 eV). The four-zone model has little to no effect on the measured total nebular abundances, but does change the interpretation of other EELG properties: we measure steeper central ionization gradients; higher volume-averaged ionization parameters; and higher central T e , n e , and log U values. Traditional three-zone estimates of the ionization parameter can underestimate the average log U by up to 0.5 dex. Additionally, we find a model-independent dichotomy in the abundance…

Abstract We present the merger rate density of Population III binary black holes (BHs) by means of a widely used binary population synthesis code BSE with extensions to very massive and extreme metal-poor stars. We consider not only low-mass BHs (lBHs: 5–50M ⊙) but also high-mass BHs (hBHs: 130–200M ⊙), where lBHs and hBHs are below and above the pair-instability mass gap (50–130M ⊙), respectively. Population III BH–BHs can be categorized into three subpopulations: BH–BHs without hBHs (hBH0s: m tot ≲ 100M ⊙), with one hBH (hBH1s: m tot ∼ 130–260M ⊙), and with two hBHs (hBH2s: m tot ∼ 270–400M ⊙), where m tot is the total mass of a BH–BH. Their merger rate densities at the current universe are ∼0.1 yr−1 Gpc−3 for hBH0s, and ∼0.01 yr−1 Gpc−3 for the sum of hBH1s and hBH2s, provided that the mass density of Population III stars is ∼1013 M ⊙ Gpc−3. These rates are modestly insensitive to initial conditions and single star models. The hBH1 and hBH2 mergers can dominate BH–BHs with hBHs discovered in the near future. They have low effective spins ≲0.2 in the current universe. The number…

Abstract Particles may be accelerated in magnetized coronae via magnetic reconnections and/or plasma turbulence, leading to high-energy neutrinos and soft γ-rays. We evaluate the detectability of neutrinos from nearby bright Seyfert galaxies identified by X-ray measurements. In the disk-corona model, we find that NGC 1068 is the most promising Seyfert galaxy in the Northern sky, where IceCube is the most sensitive, and show prospects for the identification of aggregated neutrino signals from Seyfert galaxies bright in X-rays. Moreover, we demonstrate that nearby Seyfert galaxies are promising targets for the next generation of neutrino telescopes such as KM3NeT and IceCube-Gen2. For KM3NeT, Cen A can be the most promising source in the Southern sky if a significant fraction of the observed X-rays come from the corona, and it could be identified in few years of KM3NeT operation. Our results reinforce the idea that hidden cores of supermassive black holes are the dominant sources of the high-energy neutrino emission and underlines the necessity of better sensitivity to medium-energy…

Abstract A distinct visual signature occurs in black holes that are surrounded by optically thin and geometrically thick emission regions. This signature is a sharp-edged dip in brightness that is coincident with the black hole’s shadow, which is the projection of the black hole’s unstable-photon region on the observer’s sky. We highlight two key mechanisms that are responsible for producing the sharp-edged dip: (i) the reduction of intensity observed in rays that intersect the unstable-photon region, and thus the perfectly absorbing event horizon, versus rays that do not (blocking); and (ii) the increase of intensity observed in rays that travel along extended, horizon-circling paths near the boundary of the unstable-photon region (path-lengthening). We demonstrate that the black hole shadow is a distinct phenomenon from the photon ring, and that models exist in which the former may be observed but not the latter. Additionally, we show that the black hole shadow and its associated visual signature differ from the more model-dependent brightness depressions associated with thin-disk…

Abstract The mergers of binary neutron stars, as well as black hole–neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining 32 different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at 1 day post-merger), the bolometric luminosity (one order of magnitude at 5 days post-merger), and the inferred mass of material from the bolometric luminosity (factor of 8 when considering the 8–10 day region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei…

Abstract We present multiepoch, parsec-scale core brightness temperature observations of 447 active galactic nucleus (AGN) jets from the MOJAVE and 2 cm Survey programs at 15 GHz from 1994 to 2019. The brightness temperature of each jet over time is characterized by its median value and variability. We find that the range of median brightness temperatures for AGN jets in our sample is much larger than the variations within individual jets, consistent with Doppler boosting being the primary difference between the brightness temperatures of jets in their median state. We combine the observed median brightness temperatures with apparent jet speed measurements to find the typical intrinsic Gaussian brightness temperature of 4.1( ± 0.6) × 10 10 K, suggesting that jet cores are at or below equipartition between particle and magnetic field energy in their median state. We use this value to derive estimates for the Doppler factor for every source in our sample. For the 309 jets with both apparent speed and brightness temperature data, we estimate their Lorentz factors and viewing angles to…

Abstract Accretion disks around black holes power some of the most luminous objects in the universe. Disks that are misaligned to the black hole spin can become warped over time by Lense–Thirring precession. Recent work has shown that strongly warped disks can become unstable, causing the disk to break into discrete rings producing a more dynamic and variable accretion flow. In a companion paper, we present numerical simulations of this instability and the resulting dynamics. In this paper, we discuss the implications of this dynamics for accreting black hole systems, with particular focus on the variability of active galactic nuclei (AGN). We discuss the timescales on which variability might manifest, as well as the impact of the observer orientation with respect to the black hole spin axis. When the disk warp is unstable near the inner edge of the disk, we find quasi-periodic behavior of the inner disk, which may explain the recent quasi-periodic eruptions observed in, for example, the Seyfert 2 galaxy GSN 069 and in the galactic nucleus of RX J1301.9+2747. These eruptions are…

Abstract In a companion paper, we develop a theory for the evolution of stellar wind-driven bubbles in dense, turbulent clouds. This theory proposes that turbulent mixing at a fractal bubble/shell interface leads to highly efficient cooling, in which the vast majority of the input wind energy is radiated away. This energy loss renders the majority of the bubble evolution momentum driven rather than energy driven, with expansion velocities and pressures orders of magnitude lower than in the classical Weaver et al. solution. In this paper, we validate our theory with three-dimensional, hydrodynamic simulations. We show that extreme cooling is not only possible, but is generic to star formation in turbulent clouds over more than three orders of magnitude in density. We quantify the few free parameters in our theory, and show that the momentum exceeds the wind input rate by only a factor . We verify that the bubble/cloud interface is a fractal with dimension . The measured turbulent amplitude ( ) in the hot gas near the interface is shown to be consistent with theoretical requirements…

Abstract We provide a simple geometric explanation for the source of switchbacks and associated large and one-sided transverse flows in the solar wind observed by the Parker Solar Probe (PSP). The more radial, sub-Parker spiral structure of the heliospheric magnetic field observed previously by Ulysses, ACE, and STEREO is created within rarefaction regions where footpoint motion from the source of fast into slow wind at the Sun creates a magnetic fieldline connection across solar wind speed shear. Conversely, when footpoints move from the source of slow wind into faster wind, a super-Parker spiral field structure is formed: below the Alfvén critical point, one-sided transverse field-aligned flows develop; above the Alfvén critical point, the field structure contracts between adjacent solar wind flows, and the radial field component decreases in magnitude with distance from the Sun, eventually reversing into a switchback. The sub-Parker and super-Parker spirals behave functionally as opposites. Observations from PSP confirm the paucity of switchbacks within rarefaction regions and…

Abstract The third observing run of Advanced LIGO and Advanced Virgo took place between 2019 April and 2020 March and resulted in dozens of gravitational-wave candidates, many of which are now published as confident detections. A crucial requirement of the third observing run was the rapid identification and public reporting of compact binary mergers, which enabled massive follow-up observation campaigns with electromagnetic and neutrino observatories. PyCBC Live is a low-latency search for compact binary mergers based on frequency-domain matched filtering, which was used during the second and third observing runs, together with other low-latency analyses, to generate these rapid alerts from the data acquired by LIGO and Virgo. This paper describes and evaluates the improvements made to PyCBC Live after the second observing run, which defined its operation and performance during the third observing run.

Abstract We perform a series of high-resolution 2D hydrodynamical simulations of equal-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) accretion disks to study whether these binaries can be driven to merger by the surrounding gas. We find that the gravitational softening adopted for the BBH has a profound impact on this result. When the softening is less than 10% of the binary separation, we show that, in agreement with recent simulations of isolated equal-mass binaries, prograde BBHs expand in time rather than contract. Eventually, however, the binary separation becomes large enough that the tidal force of the central AGN disrupts them. Only when the softening is relatively large do we find that prograde BBHs harden. We determine through detailed analysis of the binary torque, that this dichotomy is due to a loss of spiral structure in the circum-single disks orbiting each black hole when the softening is a significant fraction of the binary separation. Properly resolving these spirals—both with high resolution and small softening—results in a significant…

Abstract We extend the halo-based group finder developed by Yang et al. (2005c) to use data simultaneously with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N -body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude z ≤ 21 and redshift 0 < z ≤ 1.0 in the DESI legacy imaging surveys (the Legacy Surveys), our group finder successfully identifies more than 60% of the members in about 90% of halos with mass ≳10 12.5 h −1 M ⊙ . Detected groups with mass ≳10 12.0 h −1 M ⊙ have a purity (the fraction of true groups) greater than 90%. The halo mass assigned to each group has an uncertainty of about 0.2 dex at the high-mass end ≳10 13.5 h −1 M ⊙ and 0.45 dex at the low-mass end. Groups with more than 10 members have a redshift accuracy of ∼0.008. We apply this group finder to the Legacy Surveys DR8 and find 6.4 million groups with at least three members. About 500,000 of these groups have at least 10 members. The resulting catalog containing 3D coordinates, richness, halo…

Abstract The structure of a dwarf galaxy is an important probe of the effects of stellar feedback and environment. Using an unprecedented sample of 223 low-mass satellites from the ongoing Exploration of Local Volume Satellites survey, we explore the structures of dwarf satellites in the mass range 105.5 < M ⋆ < 108.5 M ⊙. We survey satellites around 80% of the massive, M K < − 22.4 mag, hosts in the Local Volume (LV). Our sample of dwarf satellites is complete to luminosities of M V <−9 mag and surface brightness μ 0,V < 26.5 mag arcsec−2 within at least ∼200 projected kpc of the hosts. For this sample, we find a median satellite luminosity of M V = −12.4 mag, median size of r e = 560 pc, median ellipticity of ϵ = 0.30, and median Sérsic index of n = 0.72. We separate the satellites into late- and early-type (29.6% and 70.4%, respectively). The mass–size relations are very similar between them within ∼5%, which indicates that the quenching and transformation of a late-type dwarf into an early-type one involves only very mild size evolution. Considering the distribution of apparent…

Abstract In this paper, we investigate the effects of a dark matter (DM) spike on the neighborhood of Sgr A*, the black hole (BH) in the center of the Milky Way. Our main goal is to investigate whether current and future astronomical observations of Sgr A* could detect the presence of such a DM spike. At first, we construct the spacetime metric around a static and spherically symmetric BH with a DM spike, and later, this solution is generalized for a rotating BH using the Newman–Janis–Azreg-Aïnou algorithm. For the static BH metric, we use the data of the S2 star orbiting Sgr A* to determine and analyze the constraints on the two free parameters characterizing the density and innermost boundary of the DM halo surrounding the BH. Furthermore, by making use of the available observational data for the DM spike density ρ sp and the DM spike radius R sp in the Milky Way, we consider a geometrically thick accretion disk model around the Sgr A* BH and demonstrate that the effect of DM distribution on the shadow radius and the image of the BH is considerably weak for realistic DM densities,…

Abstract We have conducted a search for new strong gravitational lensing systems in the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys’ Data Release 8. We use deep residual neural networks, building on previous work presented by Huang et al. These surveys together cover approximately one-third of the sky visible from the Northern Hemisphere, reaching a z -band AB magnitude of ∼22.5. We compile a training sample that consists of known lensing systems as well as non-lenses in the Legacy Surveys and the Dark Energy Survey. After applying our trained neural networks to the survey data, we visually inspect and rank images with probabilities above a threshold. Here we present 1210 new strong lens candidates.

Abstract Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the universe. Key aspects to these processes are the gas heating and cooling mechanisms, and although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Only a few detailed radiative transfer studies have been carried out owing to a lack of multiple line detections per galaxy. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies identified by the Planck satellite (LPs) at z ∼ 1.1–3.5. We analyze 162 CO rotational transitions (ranging from J up = 1 to 12) and 37 atomic carbon fine-structure lines ([C i]) in order to characterize the physical conditions of the gas in the sample of LPs. We simultaneously fit the CO and [C i] lines and the dust continuum emission, using two different non-LTE, radiative transfer models. The first model represents a two-component gas density, while the second assumes a turbulence-driven lognormal gas density distribution. These LPs are among the most…

Abstract We construct evolutionary models of the remnant of the merger of two carbon-oxygen (CO) core white dwarfs (WDs). With total masses in the range 1–2 M ⊙ , these remnants may either leave behind a single massive WD or undergo a merger-induced collapse to a neutron star (NS). On the way to their final fate, these objects generally experience a ∼10 kyr luminous giant phase, which may be extended if sufficient helium remains to set up a stable shell-burning configuration. The uncertain, but likely significant, mass-loss rate during this phase influences the final remnant mass and fate (WD or NS). We find that the initial CO core composition of the WD is converted to oxygen-neon (ONe) in remnants with final masses ≳1.05 M ⊙ . This implies that the CO core/ONe core transition in single WDs formed via mergers occurs at a similar mass as in WDs descended from single stars and thus that WD–WD mergers do not naturally provide a route to producing ultramassive CO-core WDs. As the remnant contracts toward a compact configuration, it experiences a “bottleneck” that sets the…

Abstract We use the Panoramic Survey Telescope and Rapid Response System 1 Survey (Pan-STARRS1, PS1) data to extend the Sloan Digital Sky Survey (SDSS) Stripe 82 quasar light curves. Combining PS1 and SDSS light curves provides a 15 yr baseline for 9248 quasars—5 yr longer than prior studies that used only SDSS. We fit the light curves with the damped random walk (DRW) model—a statistical description of their variability. We correlate the resulting DRW model parameters: asymptotic variability amplitude SF ∞ , and characteristic timescale τ , with quasar physical properties—black hole mass, bolometric luminosity, and redshift. Using simulated light curves, we find that a longer baseline allows us to better constrain the DRW parameters. After adding PS1 data, the variability amplitude is a stronger function of the black hole mass and has a weaker dependence on quasar luminosity. In addition, the characteristic timescale τ dependence on quasar luminosity is marginally weaker. We also make predictions for the fidelity of DRW model parameter retrieval when light curves will be further…

Abstract Magnetic braking (MB) likely plays a vital role in the evolution of low-mass X-ray binaries (LMXBs). However, there is still uncertainty around the physics of MB, and there are various proposed scenarios for MB in the literature. To examine and discriminate the efficiency of MB, we investigate the LMXB evolution with five proposed MB laws. Combining detailed binary evolution calculation with binary population synthesis, we obtain the expected properties of LMXBs and their descendant binary millisecond pulsars. We then discuss the strength and weakness of each MB law by comparing the calculated results with observations. We conclude that the τ-boosted MB law seems to best match the observational characteristics.

Abstract Using a general circulation model (GCM), we investigate trends in simulated hot Jupiter atmospheres for a range of irradiation temperatures (1500–4000 K), surface gravities (10 and 40 m s −2 ), and cloud conditions. Our models include simplified temperature-dependent clouds with radiative feedback and show how different cloud compositions, vertical thicknesses, and opacities shape hot Jupiter atmospheres by potentially increasing planetary albedos, decreasing photospheric pressures and nightside temperatures, and in some cases producing strong dayside thermal inversions. With decreasing irradiation, clouds progressively form on the nightside and cooler western limb, followed by the eastern limb and central dayside. We find that clouds significantly modify the radiative transport and affect the observable properties of planets colder than T irr ≈ 3000 K ( T eq ≈ 2100 K) depending on the clouds’ vertical extent. The precise strength of expected effects depends on the assumed parameters, but trends in predicted phase curves emerge from an ensemble of simulations. Clouds lead…

Abstract We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S 5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity–metallicity relation and is likely only now experiencing substantial stellar mass loss. Current dynamical models of Antlia 2 require it to have lost over 90% of its stars to tides, in tension…

Abstract Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen…

Abstract We leverage the 1 pc spatial resolution of the Leike et al. three-dimensional (3D) dust map to characterize the 3D structure of nearby molecular clouds ( d ≲ 400 pc). We start by “skeletonizing” the clouds in 3D volume density space to determine their “spines,” which we project on the sky to constrain cloud distances with ≈1% uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is ≈3:1. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, n , falling off like n −2 at large radii. Using Bayesian model…

Abstract We present Hubble Space Telescope ultraviolet and infrared observations of eight fast radio burst (FRB) host galaxies with subarcsecond localizations, including the hosts of three known repeating FRBs. We quantify their spatial distributions and locations with respect to their host galaxy light distributions, finding that they occur at moderate host-normalized offsets of 1.4r e ([0.6, 2.1]r e ; 68% interval) and on fainter regions of their hosts in terms of IR light but overall trace the radial distribution of IR light in their galaxies. The FRBs in our tested distribution do not clearly trace the distributions of any other transient population with known progenitors and are statistically distinct from the locations of LGRBs, H-poor SLSNe, SGRBs, and Ca-rich transients. We further find that most FRBs are not in regions of elevated local star formation rates and stellar mass surface densities in comparison to the mean global values of their hosts. We also place upper limits on the IR flux at the FRB positions of m IR ≳ 24.8–27.6 AB mag, constraining both satellite and…

Abstract We study the evolution of the binary black hole (BBH) mass distribution across cosmic time. The second gravitational-wave transient catalog (GWTC-2) from LIGO/Virgo contains BBH events out to redshifts z ∼ 1, with component masses in the range ∼5–80 M ⊙. In this catalog, the biggest BBHs, with m 1 ≳ 45 M ⊙, are only found at the highest redshifts, z ≳ 0.4. We ask whether the absence of high-mass observations at low redshift indicates that the mass distribution evolves: the biggest BBHs only merge at high redshift, and cease merging at low redshift. Modeling the BBH primary-mass spectrum as a power law with a sharp maximum mass cutoff (Truncated model), we find that the cutoff increases with redshift (> 99.9% credibility). An abrupt cutoff in the mass spectrum is expected from (pulsational) pair-instability supernova simulations; however, GWTC-2 is only consistent with a Truncated mass model if the location of the cutoff increases from at z < 0.4 to at z > 0.4. Alternatively, if the primary-mass spectrum has a break in the power law (Broken Power Law) at , rather than a…

Abstract The observed radii distribution of Kepler exoplanets reveals two distinct populations: those that are more likely to be terrestrials (≲1.7R ⊕) and those that are more likely to be gas-enveloped (≳2R ⊕). There exists a clear gap in the distribution of radii that separates these two kinds of planets. Mass-loss processes like photoevaporation by high-energy photons from the host star have been proposed as natural mechanisms to carve out this radius valley. These models favor underlying core mass function of sub-Neptunes that is sharply peaked at ∼4–8M ⊕, but the radial-velocity follow-up of these small planets hints at a more bottom-heavy mass function. By taking into account the initial gas accretion in gas-poor (but not gas-empty) nebula, we demonstrate that (1) the observed radius valley is a robust feature that is initially carved out at formation during late-time gas accretion; and (2) that it can be reconciled with core mass functions that are broad extending well into the sub-Earth regime. The maximally cooled isothermal limit prohibits cores lighter than ∼1–2M ⊕ from…

Abstract The close neighborhood of a supermassive black hole contains not only the accreting gas and dust but also stellar-sized objects, such as late-type and early-type stars and compact remnants that belong to the nuclear star cluster. When passing through the accretion flow, these objects perturb it by the direct action of stellar winds, as well as their magnetic and gravitational effects. By performing general-relativistic magnetohydrodynamic simulations, we investigate how the passages of a star can influence the supermassive black hole gaseous environment. We focus on the changes in the accretion rate and the emergence of blobs of plasma in the funnel of an accretion torus. We compare results from 2D and 3D numerical computations that have been started with comparable initial conditions. We find that a quasi-stationary inflow can be temporarily inhibited by a transiting star, and the plasmoids can be ejected along the magnetic field lines near the rotation axis. We observe the characteristic signatures of the perturbing motion in the power spectrum of the accretion…