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

Abstract We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from z = 0.001 to 2.26. This work features an increased sample size from the addition of multiple cross-calibrated photometric systems of SNe covering an increased redshift span, and improved treatments of systematic uncertainties in comparison to the original Pantheon analysis, which together result in a factor of 2 improvement in cosmological constraining power. For a flat ΛCDM model, we find Ω M = 0.334 ± 0.018 from SNe Ia alone. For a flat w 0 CDM model, we measure w 0 = −0.90 ± 0.14 from SNe Ia alone, H 0 = 73.5 ± 1.1 km s −1 Mpc −1 when including the Cepheid host distances and covariance (SH0ES), and w 0 = − 0.978 − 0.031 + 0.024 when combining the SN likelihood with Planck constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both w 0 values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark…

Abstract Here we present 1701 light curves of 1550 unique, spectroscopically confirmed Type Ia supernovae (SNe Ia) that will be used to infer cosmological parameters as part of the Pantheon+ SN analysis and the Supernovae and H 0 for the Equation of State of dark energy distance-ladder analysis. This effort is one part of a series of works that perform an extensive review of redshifts, peculiar velocities, photometric calibration, and intrinsic-scatter models of SNe Ia. The total number of light curves, which are compiled across 18 different surveys, is a significant increase from the first Pantheon analysis (1048 SNe), particularly at low redshift ( z ). Furthermore, unlike in the Pantheon analysis, we include light curves for SNe with z < 0.01 such that SN systematic covariance can be included in a joint measurement of the Hubble constant ( H 0 ) and the dark energy equation-of-state parameter ( w ). We use the large sample to compare properties of 151 SNe Ia observed by multiple surveys and 12 pairs/triplets of “SN siblings”—SNe found in the same host galaxy. Distance…

Abstract The GLASS-JWST Early Release Science (hereafter GLASS-JWST-ERS) Program will obtain and make publicly available the deepest extragalactic data of the ERS campaign. It is primarily designed to address two key science questions, namely, “what sources ionized the universe and when?” and “how do baryons cycle through galaxies?”, while also enabling a broad variety of first look scientific investigations. In primary mode, it will obtain NIRISS and NIRSpec spectroscopy of galaxies lensed by the foreground Hubble Frontier Field cluster, Abell 2744. In parallel, it will use NIRCam to observe two fields that are offset from the cluster center, where lensing magnification is negligible, and which can thus be effectively considered blank fields. In order to prepare the community for access to this unprecedented data, we describe the scientific rationale, the survey design (including target selection and observational setups), and present pre-commissioning estimates of the expected sensitivity. In addition, we describe the planned public releases of high-level data products, for use by…

Abstract CMB-S4—the next-generation ground-based cosmic microwave background (CMB) experiment—is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the universe. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semianalytic projection tool, targeted explicitly toward optimizing constraints on the tensor-to-scalar ratio, r , in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2–3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments, given a desired scientific goal. To form a closed-loop process, we…

Abstract We present a measurement of the high-energy astrophysical muon–neutrino flux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of 650k neutrino-induced muon tracks from the northern celestial hemisphere, corresponding to 9.5 yr of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and the extended model testing beyond simple power-law hypotheses. An updated treatment of systematic uncertainties and atmospheric background fluxes has been implemented based on recent models. The best-fit single power-law parameterization for the astrophysical energy spectrum results in a normalization of ϕ @ 100 TeV ν μ + ν ¯ μ = 1.44 − 0.26 + 0.25 × 10 − 18 GeV − 1 cm − 2 s − 1 sr − 1 and a spectral index γ SPL = 2.37 − 0.09 + 0.09 , constrained in the energy range from 15 TeV to 5 PeV. The model tests include a single power law with a spectral cutoff at high energies, a log-parabola model, several source-class-specific flux predictions from the literature, and a model-independent…

Abstract The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z > 6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7 deg 2 area for bright [C ii ] 158 μ m and [O iii ] 88 μ m lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z > 6.5 by ∼4–5× to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the first-year observations, 18 highly significant ≥ 7 σ [C ii ] 158 μ m lines have already been…

Abstract Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phase I survey operations with follow-up Swift UV and X-ray observations. Through our investigation into correlations between light-curve properties, we recover a shallow positive correlation between the peak bolometric luminosity and decay timescales. We introduce a new spectroscopic class of TDE, TDE-featureless, which are characterized by featureless optical spectra. The new TDE-featureless class shows larger peak bolometric luminosities, peak blackbody temperatures, and peak blackbody radii. We examine the differences between the X-ray bright and X-ray faint populations of TDEs in this sample, finding that X-ray bright TDEs show higher peak blackbody luminosities than the X-ray…

Abstract We present a detailed stellar population analysis of 11 bright (H < 26.6) galaxies at z = 9–11 (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star formation histories (SFHs), a flexible dust attenuation law, and a self-consistent model of emission lines. This approach allows us to assess how different priors affect our results and how well we can break degeneracies between dust attenuation, stellar ages, metallicity, and emission lines using data that probe only the rest-frame ultraviolet (UV) to optical wavelengths. We measure a median observed UV spectral slope β = − 1.87 − 0.43 + 0.35 for relatively massive star-forming galaxies ( 9 < log ( M ⋆ / M ⊙ ) < 10 ), consistent with no change from z = 4 to z = 9–10 at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sublinear slope (0.7 ± 0.2) and specific star formation rates…

Abstract Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission ( cigale ) is an efficient python code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of cigale (named x-cigale ) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply x-cigale to different samples, including Cosmological Evolution Survey (COSMOS) spectroscopic type 2 AGNs, Chandra Deep Field-South X-ray detected normal galaxies, Sloan Digital Sky Survey quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of x-cigale and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality…

Abstract We report upper limits on the Epoch of Reionization 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data (∼36 hr of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of ∼109 with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of Δ 21 2 ≤ ( 30.76 ) 2 mK 2 at k = 0.192 h Mpc−1 at z = 7.9, and also Δ 21 2 ≤ ( 95.74 ) 2 mK 2 at k = 0.256 h Mpc−1 at z = 10.4. At z = 7.9, these limits are the most sensitive to date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier k ∥ modes, at high k ∥ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable subreflections, or residual instrumental cross-coupling, and warrant further study. This analysis…

Abstract The Solar Probe ANalyzer for Ions (SPAN-I) onboard NASA’s Parker Solar Probe spacecraft is an electrostatic analyzer with time-of-flight capabilities that measures the ion composition and three-dimensional distribution function of the thermal corona and solar-wind plasma. SPAN-I measures the energy per charge of ions in the solar wind from 2 eV to 30 keV with a field of view of 247.°5 × 120° while simultaneously separating H + from He ++ to develop 3D velocity distribution functions of individual ion species. These observations, combined with reduced distribution functions measured by the Sun-pointed Solar Probe Cup, will help us further our understanding of the solar-wind acceleration and formation, the heating of the corona, and the acceleration of particles in the inner heliosphere. This paper describes the instrument hardware, including several innovative improvements over previous time-of-flight sensors, the data products generated by the experiment, and the ground calibrations of the sensor.

Abstract The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams, and satellite galaxies in action ( J ) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ( E ) to enhance the contrast between the groupings. For this endeavor, we use Gaia EDR3‒based measurements of 170 globular clusters, 41 streams, and 46 satellites to derive their J and E . To detect groups, we use the ENLINK software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of N = 6 groups, including the previously known mergers Sagittarius, Cetus, Gaia‒Sausage/Enceladus, LMS-1/Wukong, Arjuna/Sequoia/I’itoi, and one new merger that…

Abstract The Lyman continuum (LyC) cannot be observed at the epoch of reionization ( z ≳ 6) owing to intergalactic H i absorption. To identify LyC emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Ly α , such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O iii ]/[O ii ] flux ratio and star formation rate surface density, predict whether a galaxy is an LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H i column density, ionization parameter, and stellar feedback. We find that LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

Abstract We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at 〈 z obs 〉 = 0.8 with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector , we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including nonparametric star formation histories, emission lines, variable dust attenuation law, and dust and active galactic nucleus emission), together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust–age–metallicity degeneracy. We find a large diversity of star formation histories: although the most massive ( M ⋆ > 2 × 10 11 M ⊙ ) galaxies formed the earliest (formation redshift of z f ≈ 5–10 with a short star formation timescale of τ SF ≲ 1 Gyr), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of z f with M ⋆ . Interestingly, several low-mass galaxies have formation redshifts of z f ≈ 5–8. Star-forming galaxies evolve about the star-forming main…

Abstract We present a recalibration of the photometric systems in the Pantheon+ sample of Type Ia supernovae (SNe Ia) including those in the SH0ES distance-ladder measurement of H 0 . We utilize the large and uniform sky coverage of the public Pan-STARRS stellar photometry catalog to cross calibrate against tertiary standards released by individual SN Ia surveys. The most significant updates over the “SuperCal” cross calibration used for the previous Pantheon and SH0ES analyses are: (1) expansion of the number of photometric systems (now 25) and filters (now 105), (2) solving for all filter offsets in all systems simultaneously to produce a calibration uncertainty covariance matrix for cosmological-model constraints, and (3) accounting for the change in the fundamental flux calibration of the Hubble Space Telescope CALSPEC standards from previous versions on the order of 1.5% over a Δ λ of 4000 Å. We retrain the SALT2 model and find that our new model coupled with the new calibration of the light curves themselves causes a net distance modulus change ( d μ / dz ) of 0.04 mag over…

Abstract We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration (t E ≃ 270 days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of 6 yr, reveals a clear relativistic astrometric deflection of the background star’s apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of Earth’s motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of 7.1 ± 1.3 M ⊙ and a distance of 1.58 ± 0.18 kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion…

Abstract Recently, the Hydrogen Epoch of Reionization Array (HERA) has produced the experiment’s first upper limits on the power spectrum of 21 cm fluctuations at z ∼ 8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization from these limits. We find that the IGM must have been heated above the adiabatic-cooling threshold by z ∼ 8, independent of uncertainties about IGM ionization and the radio background. Combining HERA limits with complementary observations constrains the spin temperature of the z ∼ 8 neutral IGM to 27 K 〈 T ¯ S 〉 630 K (2.3 K 〈 T ¯ S 〉 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the cosmic microwave background dominates the z ∼ 8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones. The z ∼ 10 limits require even earlier heating if dark-matter interactions cool the hydrogen gas. If…

Abstract Gravitational-wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, R BBH ( z ). We make predictions for R BBH ( z ) as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations, COMPAS . We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterized by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about 30 M ⊙ and short delay times ( t delay ≲ 1 Gyr), while the stable RLOF channel primarily forms systems with BH masses above 30 M ⊙ and long delay times ( t delay ≳ 1 Gyr). We provide a new fit for the metallicity-dependent specific star formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction of R BBH ( z ). This leads to a distinct…

Abstract We use the panchromatic spectral energy distribution (SED)-fitting code Prospector to measure the galaxy log M *–logSFR relationship (the star-forming sequence ) across 0.2 < z < 3.0 using the COSMOS-2015 and 3D-HST UV-IR photometric catalogs. We demonstrate that the chosen method of identifying star-forming galaxies introduces a systematic uncertainty in the inferred normalization and width of the star-forming sequence, peaking for massive galaxies at ∼0.5 and ∼0.2 dex, respectively. To avoid this systematic, we instead parameterize the density of the full galaxy population in the log M *–logSFR–redshift plane using a flexible neural network known as a normalizing flow. The resulting star-forming sequence has a low-mass slope near unity and a much flatter slope at higher masses, with a normalization 0.2–0.5 dex lower than typical inferences in the literature. We show this difference is due to the sophistication of the Prospector stellar populations modeling: the nonparametric star formation histories naturally produce higher masses while the combination of individualized…

Abstract High dust density in the midplane of protoplanetary disks is favorable for efficient grain growth and can allow fast formation of planetesimals and planets, before disks dissipate. Vertical settling and dust trapping in pressure maxima are two mechanisms allowing dust to concentrate in geometrically thin and high-density regions. In this work, we aim to study these mechanisms in the highly inclined protoplanetary disk SSTC2D J163131.2-242627 (Oph 163131, i ∼ 84°). We present new high-angular-resolution continuum and 12 CO ALMA observations of Oph 163131. The gas emission appears significantly more extended in the vertical and radial direction compared to the dust emission, consistent with vertical settling and possibly radial drift. In addition, the new continuum observations reveal two clear rings. The outer ring, located at ∼100 au, is well-resolved in the observations, allowing us to put stringent constraints on the vertical extent of millimeter dust particles. We model the disk using radiative transfer and find that the scale height of millimeter-sized grains is 0.5 au…

Abstract We present a novel analytic framework to model the steady-state structure of multiphase galactic winds comprised of a hot, volume-filling component and a cold, clumpy component. We first derive general expressions for the structure of the hot phase for arbitrary mass, momentum, and energy source terms. Next, informed by recent simulations, we parameterize the cloud–wind mass transfer rates, which are set by the competition between turbulent mixing and radiative cooling. This enables us to cast the cloud–wind interaction as a source term for the hot phase and thereby simultaneously solve for the evolution of both phases, fully accounting for their bidirectional influence. With this model, we explore the nature of galactic winds over a broad range of conditions. We find that (i) with realistic parameter choices, we naturally produce a hot, low-density wind that transports energy while entraining a significant flux of cold clouds, (ii) mixing dominates the cold cloud acceleration and decelerates the hot wind, (iii) during mixing thermalization of relative kinetic energy…

Abstract We present Hubble Space Telescope (HST) photometry of 17 Cepheids in open clusters and their cluster mean parallaxes from Gaia EDR3. These parallaxes are more precise than those from individual Cepheids (G < 8 mag) previously used to measure the Hubble constant because they are derived from an average of >300 stars per cluster. Cluster parallaxes also have smaller systematic uncertainty because their stars lie in the range (G > 13 mag) where the Gaia parallax calibration is the most comprehensive. Cepheid photometry employed in the period–luminosity relation was measured using the same HST instrument (WFC3) and filters (F555W, F814W, F160W) as extragalactic Cepheids in Type Ia supernova hosts. We find no evidence of residual parallax offset in this magnitude range, zp = −3 ± 4 μas, consistent with the results from Lindegren et al. and most studies. The Cepheid luminosity (at P = 10 d and solar metallicity) in the HST near-infrared, Wesenheit magnitude system derived from the cluster sample is M H , 1 W = − 5.902 ± 0.025 mag and −5.890 ± 0.018 mag with or without…

Abstract The inner ∼200 pc region of the Galaxy contains a 4 million M⊙ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic-ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation-driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total-intensity and spectral-index mosaics of the region from 20 pointings (144 hr on-target in total), covering 6.5 square degrees with an angular resolution of 4″ at a central frequency of…

Abstract We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal Type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey. Pre-SN activity was detected in riz -bands at −130 days and persisted at relatively constant flux until first light. Soon after discovery, “flash” spectroscopy of SN 2020tlf revealed narrow, symmetric emission lines that resulted from the photoionization of circumstellar material (CSM) shed in progenitor mass-loss episodes before explosion. Surprisingly, this novel display of pre-SN emission and associated mass loss occurred in a red supergiant (RSG) progenitor with zero-age main-sequence mass of only 10–12 M ⊙ , as inferred from nebular spectra. Modeling of the light curve and multi-epoch spectra with the non-LTE radiative-transfer code CMFGEN and radiation-hydrodynamical code HERACLES suggests a dense CSM limited to r ≈ 10 15 cm, and mass-loss rate of 10 −2 M ⊙ yr −1 . The luminous light-curve plateau and persistent blue excess indicates an extended progenitor,…

Abstract We introduce GWFAST, a novel Fisher-matrix code for gravitational-wave studies, tuned toward third-generation gravitational-wave detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE). We use it to perform a comprehensive study of the capabilities of ET alone, and of a network made by ET and two CE detectors, as well as to provide forecasts for the forthcoming O4 run of the LIGO-Virgo-KAGRA (LVK) collaboration. We consider binary neutron stars, binary black holes, and neutron star–black hole binaries, and compute basic metrics such as the distribution of signal-to-noise ratio (S/N), the accuracy in the reconstruction of various parameters (including distance, sky localization, masses, spins, and, for neutron stars, tidal deformabilities), and the redshift distribution of the detections for different thresholds in S/N and different levels of accuracy in localization and distance measurement. We examine the expected distribution and properties of golden events, with especially large values of the S/N. We also pay special attention to the dependence of the results…

Abstract Various profiles of matter distribution in galactic halos (such as the Navarro–Frenk–White, Burkert, Hernquist, Moore, Taylor–Silk models, and others) are considered here as the source term for the Einstein equations. We solve these equations and find exact solutions that represent the metric of a central black hole immersed in a galactic halo. Even though in the general case the solution is numerical, very accurate general analytical metrics, which include all the particular models, are found in the astrophysically relevant regime, when the mass of the galaxy is much smaller than the characteristic scale in the halo.

Abstract We combine Gaia early data release 3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we use a Monte Carlo sample for the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare…

Abstract We present a comprehensive optical and near-infrared census of the fields of 90 short gamma-ray bursts (GRBs) discovered in 2005–2021, constituting all short GRBs for which host galaxy associations are feasible (≈60% of the total Swift short GRB population). We contribute 274 new multi-band imaging observations across 58 distinct GRBs and 26 spectra of their host galaxies. Supplemented by literature and archival survey data, the catalog contains 542 photometric and 42 spectroscopic data sets. The photometric catalog reaches 3 σ depths of ≳24–27 mag and ≳23–26 mag for the optical and near-infrared bands, respectively. We identify host galaxies for 84 bursts, in which the most robust associations make up 56% (50/90) of events, while only a small fraction, 6.7%, have inconclusive host associations. Based on new spectroscopy, we determine 18 host spectroscopic redshifts with a range of z ≈ 0.15–1.5 and find that ≈23%–41% of Swift short GRBs originate from z > 1. We also present the galactocentric offset catalog for 84 short GRBs. Taking into account the large range of…

Abstract Our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population, which reflects the star formation and enrichment in the few most massive progenitors that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet provide a comprehensive picture of such a metal-poor “heart” of the Milky Way. We use information from Gaia Data Release 3, especially the XP spectra, to construct a sample of 2 million bright ( G BP < 15.5 mag) giant stars within 30° of the Galactic center (GC) with robust [M/H] estimates, δ [M/H] ≲ 0.1. For ∼1.25 million stars we calculate orbits from Gaia Radial Velocity Spectrometer velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes ∼18,000 stars with −2.7 < [M/H] < −1.5, representing a stellar mass of ≳5 × 10 7 M ⊙ . The spatial distribution of these [M/H] < −1.5 stars has a Gaussian extent of only σ R GC ∼ 2.7 kpc around the GC, with most orbits confined to the…

Abstract We present homogeneous abundance analysis of over 20 elements for 385 very metal-poor (VMP) stars based on the LAMOST survey and follow-up observations with the Subaru Telescope. It is the largest high-resolution VMP sample (including 363 new objects) studied by a single program, and the first attempt to accurately determine evolutionary stages for such a large sample based on Gaia parallaxes. The sample covers a wide metallicity range from [Fe/H] ≲ −1.7 down to [Fe/H] ∼ −4.3, including over 110 objects with [Fe/H] ≤ −3.0. The expanded coverage in evolutionary status makes it possible to define the abundance trends respectively for giants and turnoff stars. The newly obtained abundance data confirm most abundance trends found by previous studies, but also provide useful updates and new samples of outliers. The Li plateau is seen in main-sequence turnoff stars with −2.5 < [Fe/H] < −1.7 in our sample, whereas the average Li abundance is clearly lower at lower metallicity. Mg, Si, and Ca are overabundant with respect to Fe, showing decreasing trend with increasing metallicity.…

Abstract The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0–1.93 μ m, R = 435–865) and the 0.28 deg 2 wide-field camera. We present a reference HLSS that maps 2000 deg 2 and achieves an emission-line flux limit of 10 −16 erg s −1 cm −2 at 6.5 σ , requiring ∼0.6 yr of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semianalytic galaxy formation model to a cosmological N -body simulation and use this mock survey to create pixel-level simulations of 4 deg 2 of HLSS grism spectroscopy. We find that the reference HLSS would measure ∼10 million H α galaxy redshifts that densely map large-scale structure at z = 1–2 and 2 million [O iii ] galaxy redshifts that sparsely map structures at z = 2–3. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the…

Abstract We present the results from a new search for candidate galaxies at z ≈ 8.5–11 discovered over the 850 arcmin2 area probed by the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). We use a photometric-redshift selection including both Hubble and Spitzer Space Telescope photometry to robustly identify galaxies in this epoch at H 160 < 26.6. We use a detailed vetting procedure, including screening against persistence and stellar contamination, and the inclusion of ground-based imaging and follow-up Hubble Space Telescope imaging to build a robust sample of 11 candidate galaxies, three presented here for the first time. The inclusion of Spitzer/IRAC photometry in the selection process reduces contamination, and yields more robust redshift estimates than Hubble alone. We constrain the evolution of the rest-frame ultraviolet luminosity function via a new method of calculating the observed number densities without choosing a prior magnitude bin size. We find that the abundance at our brightest probed luminosities (M UV = − 22.3) is consistent with…

Abstract The growing population of compact binary mergers detected with gravitational waves contains multiple events that are challenging to explain through isolated binary evolution. Such events have higher masses than are expected in isolated binaries, component spin tilt angles that are misaligned, and/or nonnegligible orbital eccentricities. We investigate the orbital eccentricities of 62 binary black hole candidates from the third gravitational-wave transient catalog of the LIGO–Virgo–KAGRA Collaboration with an aligned-spin, moderate-eccentricity waveform model. Within this framework, we find that at least four of these events show significant support for eccentricity e 10 ≥ 0.1 at a gravitational-wave frequency of 10 Hz (>60% credibility, under a log-uniform eccentricity prior that spans the range 10−4 < e 10 < 0.2). Two of these events are new additions to the population: GW191109 and GW200208_22. If the four eccentric candidates are truly eccentric, our results suggest that densely populated star clusters may produce 100% of the observed mergers. However, it remains likely…

Abstract The ultraviolet (UV) continuum slope ( β , where f λ ∝ λ β ) of galaxies is sensitive to a variety of properties, from the metallicity and age of the stellar population to dust attenuation throughout the galaxy. Considerable attention has focused on identifying reionization-era galaxies with very blue UV slopes ( β < −3). Not only do such systems provide a signpost of low-metallicity stars, but they also identify galaxies likely to leak ionizing photons from their H ii regions as such blue UV slopes require the reddening effect of nebular continuum to be diminished. In this paper we present a search for reionization-era galaxies with very blue UV colors in recent JWST/NIRCam imaging of the Extended Groth Strip field. We characterize UV slopes for a large sample of z ≃ 7–11 galaxies, finding a median of β = −2.0. Two lower luminosity ( M UV ≃ −19.5) and lower stellar mass (6–10 × 10 7 M ⊙ ) systems exhibit extremely blue UV slopes ( β = −2.9 to −3.1) and rest-optical photometry indicating weak nebular line emission. Each system is very compact ( r e ≲ 260 pc) with very high…

Abstract We present new determinations of the rest-UV luminosity functions (LFs) at z = 2–9 to extremely low luminosities (>−14 mag) from a sample of >2500 lensed galaxies found behind the Hubble Frontier Fields (HFF) clusters. For the first time, we present faint-end slope results from lensed samples that are fully consistent with blank-field results over the redshift range z = 2–9, while reaching to much lower luminosities than possible from the blank-field studies. Combining the deep lensed sample with the large blank-field samples allows us to set tight constraints on the faint-end slope α of the z = 2–9 UV LFs and its evolution. We find a smooth flattening in α from −2.28 ± 0.10 (z = 9) to −1.53 ± 0.03 (z = 2) with cosmic time (dα/dz = −0.11 ± 0.01), fully consistent with dark matter halo buildup. We utilize these new results to present new measurements of the evolution in the UV luminosity density ρ UV brighter than −13 mag from z ∼ 9 to z ∼ 2. Accounting for the star formation rate (SFR) densities to faint luminosities implied by our LF results, we find that unobscured star…

Abstract We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot…

Abstract We measure the low- J CO line ratios R 21 ≡ CO (2–1)/CO (1–0), R 32 ≡ CO (3–2)/CO (2–1), and R 31 ≡CO (3–2)/CO (1–0) using whole-disk CO maps of nearby galaxies. We draw CO (2–1) from PHANGS-ALMA, HERACLES, and follow-up IRAM surveys; CO (1–0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO (3–2) from the James Clerk Maxwell Telescope Nearby Galaxy Legacy Survey and Atacama Pathfinder Experiment Large APEX Sub-Millimetre Array mapping. All together, this yields 76, 47, and 29 maps of R 21 , R 32 , and R 31 at 20″ ∼ 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log ( M ⋆ / M ⊙ ) = 10.25 – 11 , and star formation rate (SFR) = 1–5 M ⊙ yr −1 , dominate the sample. We find galaxy-integrated mean values and a 16%–84% range of R 21 = 0.65 (0.50–0.83), R 32 = 0.50 (0.23–0.59), and R 31 = 0.31 (0.20–0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/ M ⋆ and SFR/ L CO . Within galaxies, we measure central enhancements with…

Abstract We present measurements of black hole masses and Eddington ratios (λ Edd) for a sample of 38 bright (M 1450 < −24.4 mag) quasars at 5.8 ≲ z ≲ 7.5, derived from Very Large Telescope/X–shooter near–IR spectroscopy of their broad C iv and Mg ii emission lines. The black hole masses (on average, M BH ∼ 4.6 × 109 M ⊙) and accretion rates (0.1 ≲ λ Edd ≲ 1.0) are broadly consistent with that of similarly luminous 0.3 ≲ z ≲ 2.3 quasars, but there is evidence for a mild increase in the Eddington ratio above z ≳ 6. Combined with deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [C II] 158 μm line from the host galaxies and VLT/MUSE investigations of the extended Lyα halos, this study provides fundamental clues to models of the formation and growth of the first massive galaxies and black holes. Compared to local scaling relations, z ≳ 5.7 black holes appear to be over-massive relative to their hosts, with accretion properties that do not change with host galaxy morphologies. Assuming that the kinematics of the T ∼ 104 K gas, traced by the extended Lyα halos,…

Abstract Separating the components of redshift due to expansion and peculiar motion in the nearby universe (z < 0.1) is critical for using Type Ia Supernovae (SNe Ia) to measure the Hubble constant (H 0) and the equation-of-state parameter of dark energy (w). Here, we study the two dominant “motions” contributing to nearby peculiar velocities: large-scale, coherent-flow (CF) motions and small-scale motions due to gravitationally associated galaxies deemed to be in a galaxy group. We use a set of 584 low-z SNe from the Pantheon+ sample, and evaluate the efficacy of corrections to these motions by measuring the improvement of SN distance residuals. We study multiple methods for modeling the large and small-scale motions and show that, while group assignments and CF corrections individually contribute to small improvements in Hubble residual scatter, the greatest improvement comes from the combination of the two (relative standard deviation of the Hubble residuals, Rel. SD, improves from 0.167 to 0.157 mag). We find the optimal flow corrections derived from various local density maps…

Abstract Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo’s third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and…

Abstract The tip of the red giant branch has been used to measure distances to 500 nearby galaxies with the Hubble Space Telescope (HST) which are available in the Color–Magnitude Diagrams and Tip of the Red Giant Branch (CMDs/TRGB) catalog on the Extragalactic Distance Database (EDD). Our established methods are employed to perform an independent reduction of the targets presented by the Carnegie-Chicago Hubble Program (CCHP) in the series of papers culminating in Freedman (2021). Our distinct methodology involves modeling the observed luminosity function of red giant branch and asymptotic giant branch stars, which differs from the edge-detection algorithms employed by the CCHP. We find excellent agreement between distances for 11 hosts with new imaging, all at D < 20 Mpc. However, we are unable to measure the TRGB for four hosts that use archival data designed to measure distances with Cepheids, all at D > 23 Mpc. With two new HST observations taken in the halo of the megamaser host NGC 4258, the first with the same ACS F606W and F814W filters and state of the electronics used for…

Abstract Ground-based gravitational-wave detectors like Cosmic Explorer (CE) can be tuned to improve their sensitivity at high or low frequencies by tuning the response of the signal extraction cavity. Enhanced sensitivity above 2 kHz enables measurements of the post-merger gravitational-wave spectrum from binary neutron star mergers, which depends critically on the unknown equation of state of hot, ultra-dense matter. Improved sensitivity below 500 Hz favors precision tests of extreme gravity with black hole ringdown signals and improves the detection prospects while facilitating an improved measurement of source properties for compact binary inspirals at cosmological distances. At intermediate frequencies, a more sensitive detector can better measure the tidal properties of neutron stars. We present and characterize the performance of tuned CE configurations that are designed to optimize detections across different astrophysical source populations. These tuning options give CE the flexibility to target a diverse set of science goals with the same detector infrastructure. We find…

Abstract The star formation rate (SFR) in galactic disks depends on both the quantity of the available interstellar medium (ISM) gas and its physical state. Conversely, the ISM’s physical state depends on the SFR, because the “feedback” energy and momentum injected by recently formed massive stars is crucial to offsetting losses from turbulent dissipation and radiative cooling. The ISM’s physical state also responds to the gravitational field that confines it, with increased weight driving higher pressure. In a quasi-steady state, it is expected that the mean total pressure of different thermal phases will match each other, that the component pressures and total pressure will satisfy thermal and dynamical equilibrium requirements, and that the SFR will adjust as needed to provide the requisite stellar radiation and supernova feedback. The pressure-regulated, feedback-modulated (PRFM) theory of the star-forming ISM formalizes these ideas, leading to a prediction that the SFR per unit area, Σ SFR , will scale nearly linearly with ISM weight  . In terms of the large-scale gas surface…

Abstract We present two bright galaxy candidates at z ∼ 12–13 identified in our H-dropout Lyman break selection with 2.3 deg2 near-infrared deep imaging data. These galaxy candidates, selected after careful screening of foreground interlopers, have spectral energy distributions showing a sharp discontinuity around 1.7 μm, a flat continuum at 2–5 μm, and nondetections at <1.2 μm in the available photometric data sets, all of which are consistent with a z > 12 galaxy. An ALMA program targeting one of the candidates shows a tentative 4σ [O iii] 88 μm line at z = 13.27, in agreement with its photometric redshift estimate. The number density of the z ∼ 12–13 candidates is comparable to that of bright z ∼ 10 galaxies and is consistent with a recently proposed double-power-law luminosity function rather than the Schechter function, indicating little evolution in the abundance of bright galaxies from z ∼ 4 to 13. Comparisons with theoretical models show that the models cannot reproduce the bright end of rest-frame ultraviolet luminosity functions at z ∼ 10–13. Combined with recent studies…

Abstract The primordial He abundance Y P is a powerful probe of cosmology. Currently, Y P is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor (<0.1Z ⊙) galaxies (EMPGs) having reliable He/H measurements with He i λ10830 near-infrared (NIR) emission. Here we present deep Subaru NIR spectroscopy for 10 EMPGs. Combining the existing optical data, He/H values of 5 out of the 10 EMPGs are reliably derived by the Markov chain Monte Carlo algorithm. Adding the existing 3 EMPGs and 51 moderately metal-poor (0.1–0.4Z ⊙) galaxies with reliable He/H estimates, we obtain Y P = 0.2370 − 0.0034 + 0.0033 by linear regression in the (He/H) − (O/H) plane, where we increase the number of EMPGs from three to eight anchoring He/H of the most metal-poor gas in galaxies. Although our Y P measurement and previous measurements are consistent, our result is slightly (∼1σ) smaller due to our EMPGs. Including the existing primordial deuterium D P measurement, we constrain the effective number of neutrino species N eff and the baryon-to-photon…

Abstract The Event Horizon Telescope (EHT) observation unveiled the first image of supermassive black hole Sgr A* showing a shadow of diameter θ sh = 48.7 ± 7 μ as with fractional deviation from the Schwarzschild black hole shadow diameter δ = − 0.08 − 0.09 + 0.09 ( VLTI ) , − 0.04 − 0.10 + 0.09 ( Keck ) . The Sgr A* shadow size is within 10% of the Kerr predictions, providing us with another tool to investigate the nature of strong-field gravity. We use the Sgr A* shadow observables to constrain metrics of four independent and well-motivated, parametrically different from Kerr spacetime, rotating regular spacetimes and the corresponding no-horizon spacetimes. We present constraints on the deviation parameter g of rotating regular black holes. The shadow angular diameter θ sh within the 1 σ region places bounds on the parameters a and g . Together with EHT bounds on the θ sh and δ of Sgr A*, our analysis concludes that the three rotating regular black holes, i.e., Bardeen, Hayward, and Simpson–Visser black holes, and corresponding no-horizon spacetimes agree with the EHT results of…

Abstract A star that approaches a supermassive black hole (SMBH) on a circular extreme mass ratio inspiral (EMRI) can undergo Roche lobe overflow (RLOF), resulting in a phase of long-lived mass transfer onto the SMBH. If the interval separating consecutive EMRIs is less than the mass-transfer timescale driven by gravitational wave emission (typically ∼1–10 Myr), the semimajor axes of the two stars will approach each another on scales of ≲ hundreds to thousands of gravitational radii. Close flybys tidally strip gas from one or both RLOFing stars, briefly enhancing the mass-transfer rate onto the SMBH and giving rise to a flare of transient X-ray emission. If both stars reside in a common orbital plane, these close interactions will repeat on a timescale as short as hours, generating a periodic series of flares with properties (amplitudes, timescales, sources lifetimes) remarkably similar to the “quasi-periodic eruptions” (QPEs) recently observed from galactic nuclei hosting low-mass SMBHs. A cessation of QPE activity is predicted on a timescale of months to years, due to nodal…

Abstract Gamma-ray bursts (GRBs) are classified into long and short events. Long GRBs (LGRBs) are associated with the end states of very massive stars, while short GRBs (SGRBs) are linked to the merger of compact objects. GRB 200826A was a peculiar event, because by definition it was an SGRB, with a rest-frame duration of ∼0.5 s. However, this event was energetic and soft, which is consistent with LGRBs. The relatively low redshift ( z = 0.7486) motivated a comprehensive, multiwavelength follow-up campaign to characterize its host, search for a possible associated supernova (SN), and thus understand the origin of this burst. To this aim we obtained a combination of deep near-infrared (NIR) and optical imaging together with spectroscopy. Our analysis reveals an optical and NIR bump in the light curve whose luminosity and evolution are in agreement with several SNe associated to LGRBs. Analysis of the prompt GRB shows that this event follows the E p,i – E iso relation found for LGRBs. The host galaxy is a low-mass star-forming galaxy, typical of LGRBs, but with one of the highest star…

Abstract The recent theoretical advance known as the minimal geometric deformation (MGD) method has initiated renewed interest in investigating higher-curvature gravitational effects in relativistic astrophysics. In this work, we model a strange star within the context of Einstein–Gauss–Bonnet gravity with the help of the MGD technique. Starting off with the Tolman metric ansatz, together with the MIT bag model equation of state applicable to hadronic matter, anisotropy is introduced via the superposition of the seed source and the decoupled energy-momentum tensor. The solution of the governing systems of equations bifurcates into two distinct models, namely, the mimicking of the θ sector to the seed radial pressure and energy density and a regular fluid model. Each of these models can be interpreted as self-gravitating static, compact objects with the exterior described by the vacuum Boulware–Deser solution. Utilizing observational data for three stellar candidates, namely PSR J1614–2230, PSR J1903+317, and LMC X-4, we subject our solutions to rigorous viability tests based on…

Abstract The Event Horizon Telescope (EHT) collaboration recently unveiled the first image of the supermassive black hole M87*, which exhibited a ring of angular diameter θ d = 42 ± 3 μ as, a circularity deviation of Δ C ≤ 0.1, and also inferred a black hole mass of M = (6.5 ± 0.7) × 10 9 M ⊙ . This provides a new window onto tests of theories of gravity in the strong-field regime, including probes of violations of the no-hair theorem. It is widely believed that the Kerr metric describes astrophysical black holes, as encapsulated in the critical but untested no-hair theorem. Modeling Horndeski gravity black holes—with an additional hair parameter h besides the mass M and spin a of the Kerr black hole—as the supermassive black hole M87*, we observe that to be a viable astrophysical black hole candidate, the EHT result constrains the ( a , h ) parameter space. However, a systematic bias analysis indicates that rotating Horndeski black hole shadows may or may not capture Kerr black hole shadows, depending on the parameter values; the latter is the case over a substantial part of the…

Abstract The observations of GW170817/AT2017gfo have provided us with evidence that binary neutron star mergers are sites of r-process nucleosynthesis. However, the observed signatures in the spectra of GW170817/AT2017gfo have not been fully decoded, especially in the near-infrared (NIR) wavelengths. In this paper, we investigate the kilonova spectra over the entire wavelength range with the aim of elemental identification. We systematically calculate the strength of bound–bound transitions by constructing a hybrid line list that is accurate for important strong transitions and complete for weak transitions. We find that the elements on the left side of the periodic table, such as Ca, Sr, Y, Zr, Ba, La, and Ce, tend to produce prominent absorption lines in the spectra. This is because such elements have a small number of valence electrons and low-lying energy levels, resulting in strong transitions. By performing self-consistent radiative transfer simulations for the entire ejecta, we find that La iii and Ce iii appear in the NIR spectra, which can explain the absorption features at…

Abstract A number of studies based on the data collected by the Hubble Space Telescope (HST) GO-13297 program “HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation” have investigated the photometric properties of a large sample of Galactic globular clusters and revolutionized our understanding of their stellar populations. In this paper, we expand upon previous studies by focusing our attention on the stellar clusters’ internal kinematics. We computed proper motions for stars in 56 globular clusters and one open cluster by combining the GO-13297 images with archival HST data. The astrophotometric catalogs released with this paper represent the most complete and homogeneous collection of proper motions of stars in the cores of stellar clusters to date, and expand the information provided by the current (and future) Gaia data releases to much fainter stars and into the crowded central regions. We also census the general kinematic properties of stellar clusters by computing the velocity dispersion and anisotropy radial profiles of their…

Abstract We report the kinematic, orbital, and chemical properties of 12 stellar streams with no evident progenitors using line-of-sight velocities and metallicities from the Southern Stellar Stream Spectroscopic Survey (S 5), proper motions from Gaia EDR3, and distances derived from distance tracers or the literature. This data set provides the largest homogeneously analyzed set of streams with full 6D kinematics and metallicities. All streams have heliocentric distances between ∼10 and 50 kpc. The velocity and metallicity dispersions show that half of the stream progenitors were disrupted dwarf galaxies (DGs), while the other half originated from disrupted globular clusters (GCs), hereafter referred to as DG and GC streams. Based on the mean metallicities of the streams and the mass–metallicity relation, the luminosities of the progenitors of the DG streams range between those of Carina and Ursa Major I (−9.5 ≲ M V ≲ −5.5). Four of the six GC streams have mean metallicities of [Fe/H] < −2, more metal poor than typical Milky Way (MW) GCs at similar distances. Interestingly, the…

Abstract We report the discovery of MAGAZ3NE J095924+022537, a spectroscopically confirmed protocluster at z = 3.3665 − 0.0012 + 0.0009 around a spectroscopically confirmed UVJ -quiescent ultramassive galaxy (UMG; M ⋆ = 2.34 − 0.34 + 0.23 × 10 11 M ⊙ ) in the COSMOS UltraVISTA field. We present a total of 38 protocluster members (14 spectroscopic and 24 photometric), including the UMG. Notably, and in marked contrast to protoclusters previously reported at this epoch that have been found to contain predominantly star-forming members, we measure an elevated fraction of quiescent galaxies relative to the coeval field ( 73.3 − 16.9 + 26.7 % versus 11.6 − 4.9 + 7.1 % for galaxies with stellar mass M ⋆ ≥ 10 11 M ⊙ ). This high quenched fraction provides a striking and important counterexample to the seeming ubiquitousness of star-forming galaxies in protoclusters at z > 2 and suggests, rather, that protoclusters exist in a diversity of evolutionary states in the early universe. We discuss the possibility that we might be observing either “early mass quenching” or nonclassical…

Abstract We combine spectroscopic, photometric, and astrometric information from APOGEE data release 17 and Gaia early data release 3 to perform a self-consistent characterization of Gaia-Sausage/Enceladus (GSE), the remnant of the last major merger experienced by the Milky Way, considering stars and globular clusters (GCs) altogether. Our novel set of chemodynamical criteria to select genuine stars of GSE yields a metallicity distribution function with a median [Fe/H] of −1.22 and 0.23 dex dispersion. Stars from GSE present an excess of [Al/Fe] and [Mg/Mn] (also [Mg/Fe]) in comparison to surviving Milky Way dwarf satellites, which can be explained by differences in star formation efficiencies and timescales between these systems. However, stars from Sequoia, another proposed accreted halo substructure, essentially overlap the GSE footprint in all analyzed chemical-abundance spaces, but present lower metallicities. Among the probable GCs of GSE with APOGEE observations available, we find no evidence for atypical [Fe/H] spreads with the exception of ω Centauri (ωCen). Under the…

Abstract We present a new investigation of the intergalactic medium near reionization using dark gaps in the Lyβ forest. With its lower optical depth, Lyβ offers a potentially more sensitive probe to any remaining neutral gas compared to the commonly used Lyα line. We identify dark gaps in the Lyβ forest using spectra of 42 QSOs at z em > 5.5, including new data from the XQR-30 VLT Large Programme. Approximately 40% of these QSO spectra exhibit dark gaps longer than 10 h −1 Mpc at z ≃ 5.8. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Lyα forest are caused solely by ionizing ultraviolet background fluctuations and with models that include large neutral hydrogen patches at z < 6 due to a late end to reionization. Of particular interest is a very long (L = 28 h −1 Mpc) and dark (τ eff ≳ 6) gap persisting down to z ≃ 5.5 in the Lyβ forest of the z = 5.85 QSO PSO J025−11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of 〈x H I〉 ≳ 5%…

Abstract We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly modeled sample to date. Using the Prospector stellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of z = 0.64 − 0.32 + 0.83 (68% confidence), including nine photometric redshifts at z ≳ 1. We further find a median mass-weighted age of t m = 0.8 − 0.53 + 2.71 Gyr, stellar mass of log( M * / M ⊙ ) = 9.69 − 0.65 + 0.75 , star formation rate of SFR = 1.44 − 1.35 + 9.37 M ⊙ yr −1 , stellar metallicity of log( Z * / Z ⊙ ) = − 0.38 − 0.42 + 0.44 , and dust attenuation of A V = 0.43 − 0.36 + 0.85 mag (68% confidence). Overall, the majority of short GRB hosts are star-forming (≈84%), with small fractions that are either transitioning (≈6%) or quiescent (≈10%); however, we observe a much larger fraction (≈40%) of quiescent and transitioning hosts at z ≲ 0.25, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal…

Abstract Giant planets can interact with multiple and chemically diverse environments in protoplanetary disks while they form and migrate to their final orbits. The way this interaction affects the accretion of gas and solids shapes the chemical composition of the planets and of their atmospheres. Here we investigate the effects of different chemical structures of the host protoplanetary disk on the planetary composition. We consider both scenarios of molecular (inheritance from the prestellar cloud) and atomic (complete chemical reset) initial abundances in the disk. We focus on four elemental tracers of different volatility: C, O, N, and S. We explore the entire extension of possible formation regions suggested by observations by coupling the disk chemical scenarios with N-body simulations of forming and migrating giant planets. The planet formation process produces giant planets with chemical compositions significantly deviating from that of the host disk. We find that the C/N, N/O, and S/N ratios follow monotonic trends with the extent of migration. The C/O ratio shows a more…

Abstract Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly 50 M ⊙ and 100 M ⊙, while above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational-wave emission with the current LIGO/Virgo/Kagra interferometers and future detectors such as LISA or the Einstein Telescope, can form both mass-gap BHs and IMBHs. Here we explore the possibility that mass-gap BHs and IMBHs are born as a result of successive BH mergers in dense star clusters. In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products after they receive significant recoil kicks due to anisotropic emission of gravitational radiation. Using for the first time simulations that include full stellar evolution, we show that a massive stellar BH seed can easily grow to ∼103–104 M ⊙ as a result of repeated mergers with other smaller BHs. We find that lowering the cluster metallicity leads to…

Abstract Strong gravitational lensing of gravitational wave sources offers a novel probe of both the lens galaxy and the binary source population. In particular, the strong lensing event rate and the time-delay distribution of multiply imaged gravitational-wave binary coalescence events can be used to constrain the mass distribution of the lenses as well as the intrinsic properties of the source population. We calculate the strong lensing event rate for a range of second- (2G) and third-generation (3G) detectors, including Advanced LIGO/Virgo, A+, Einstein Telescope (ET), and Cosmic Explorer (CE). For 3G detectors, we find that ∼0.1% of observed events are expected to be strongly lensed. We predict detections of ∼1 lensing pair per year with A+, and ∼50 pairs per year with ET/CE. These rates are highly sensitive to the characteristic galaxy velocity dispersion, σ * , implying that observations of the rates will be a sensitive probe of lens properties. We explore using the time-delay distribution between multiply imaged gravitational-wave sources to constrain properties of the…

Abstract We present observations of SN 2021csp, the second example of a newly identified type of supernova (SN) hallmarked by strong, narrow, P Cygni carbon features at early times (Type Icn). The SN appears as a fast and luminous blue transient at early times, reaching a peak absolute magnitude of −20 within 3 days due to strong interaction between fast SN ejecta (v ≈ 30,000 km s−1) and a massive, dense, fast-moving C/O wind shed by the WC-like progenitor months before explosion. The narrow-line features disappear from the spectrum 10–20 days after explosion and are replaced by a blue continuum dominated by broad Fe features, reminiscent of Type Ibn and IIn supernovae and indicative of weaker interaction with more extended H/He-poor material. The transient then abruptly fades ∼60 days post-explosion when interaction ceases. Deep limits at later phases suggest minimal heavy-element nucleosynthesis, a low ejecta mass, or both, and imply an origin distinct from that of classical Type Ic SNe. We place SN 2021csp in context with other fast-evolving interacting transients, and discuss…

Abstract The filamentary nature of accretion streams found around embedded sources suggests that protostellar disks experience heterogenous infall from the star-forming environment, consistent with the accretion behavior onto star-forming cores in top-down star-cluster formation simulations. This may produce disk substructures in the form of rings, gaps, and spirals that continue to be identified by high-resolution imaging surveys in both embedded Class 0/I and later Class II sources. We present a parameter study of anisotropic infall, informed by the properties of accretion flows onto protostellar cores in numerical simulations, and varying the relative specific angular momentum of incoming flows as well as their flow geometry. Our results show that anisotropic infall perturbs the disk and readily launches the Rossby wave instability. It forms vortices at the inner and outer edges of the infall zone where material is deposited. These vortices drive spiral waves and angular momentum transport, with some models able to drive stresses corresponding to a viscosity parameter on the…

Abstract A promising energy range to look for angular correlations between cosmic rays of extragalactic origin and their sources is at the highest energies, above a few tens of EeV (1 EeV ≡ 1018 eV). Despite the flux of these particles being extremely low, the area of ∼3000 km2 covered at the Pierre Auger Observatory, and the 17 yr data-taking period of the Phase 1 of its operations, have enabled us to measure the arrival directions of more than 2600 ultra-high-energy cosmic rays above 32 EeV. We publish this data set, the largest available at such energies from an integrated exposure of 122,000 km2 sr yr, and search it for anisotropies over the 3.4π steradians covered with the Observatory. Evidence for a deviation in excess of isotropy at intermediate angular scales, with ∼15° Gaussian spread or ∼25° top-hat radius, is obtained at the 4σ significance level for cosmic-ray energies above ∼40 EeV.

Abstract Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap at ∼60–130 M ⊙. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing run (O4), we perform Bayesian analysis on quasi-circular nonprecessing, spinning IMBH binaries (IMBHBs) with total masses 50–500 M ⊙, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at 2σ, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of ∼10%–40% at a signal-to-noise ratio of 20. Focusing on the stellar-mass gap with new tabulations of the 12C(α, γ)16O reaction rate and its uncertainties, we evolve massive helium core stars using MESA to establish the lower and upper edges of the mass gap as ≃ 59 − 13 + 34 M ⊙ and ≃ 139 − 14 + 30 M ⊙ respectively, where the error bars give the mass range that follows from the ±3σ uncertainty in the 12C(α, γ)16O nuclear reaction…

Abstract Multiwavelength observations indicate that some starburst galaxies show a dominant nonthermal contribution from their central region. These active galactic nuclei (AGN)-starburst composites are of special interest, as both phenomena on their own are potential sources of highly energetic cosmic rays and associated γ -ray and neutrino emission. In this work, a homogeneous, steady-state two-zone multimessenger model of the nonthermal emission from the AGN corona as well as the circumnuclear starburst region is developed and subsequently applied to the case of NGC 1068, which has recently shown some first indications of high-energy neutrino emission. Here, we show that the entire spectrum of multimessenger data—from radio to γ -rays including the neutrino constraint—can be described very well if both, starburst and AGN corona, are taken into account. Using only a single emission region is not sufficient.

Abstract The Cosmic Origins Spectrograph (COS) Legacy Archive Spectroscopic SurveY (CLASSY) provides the first high-resolution spectral catalog of 45 local high-z analogs in the ultraviolet (UV; 1200–2000 Å) to investigate their stellar and gas properties. Here we present a toolkit of UV interstellar medium (ISM) diagnostics, analyzing the main emission lines of CLASSY spectra (N iv] λ λ1483,87, C iv λλ1548,51, He ii λ1640, O iii]λ λ1661,6, Si iii] λλ1883,92, C iii] λ1907,9). Specifically, our aim is to provide accurate diagnostics for the reddening E(B − V), electron density n e , electron temperature T e , metallicity 12+log(O/H), and ionization parameter log(U), taking the different ISM ionization zones into account. We calibrate our UV toolkit using well-known optical diagnostics, analyzing archival optical spectra for all CLASSY targets. We find that UV density diagnostics estimate n e values that are ∼1–2 dex higher (e.g., n e (C iii]λ λ1907,9) ∼ 104 cm−3) than those inferred from their optical counterparts (e.g., n e ([S ii]λ λ6717,31) ∼ 102 cm−3; n e ([Ar iv]λ λ4714,41) ∼…

Abstract Natal kicks and spins are characteristic properties of neutron stars (NSs) and black holes (BHs). Both offer valuable clues to dynamical processes during stellar core collapse and explosion. Moreover, they influence the evolution of stellar multiple systems and the gravitational-wave signals from their inspiral and merger. Observational evidence of a possibly generic spin-kick alignment has been interpreted as an indication that NS spins are either induced with the NS kicks or inherited from the progenitor rotation, which thus might play a dynamically important role during stellar collapse. Current three-dimensional supernova simulations suggest that NS kicks are transferred in the first seconds of the explosion, mainly by anisotropic mass ejection and, on a secondary level, anisotropic neutrino emission. By contrast, the NS spins are only determined minutes to hours later by the angular momentum associated with the fallback of matter that does not become gravitationally unbound in the supernova. Here, we propose a novel scenario to explain spin-kick alignment as a…

Abstract Population III stars ended the cosmic dark ages and began early cosmological reionization and chemical enrichment. However, in spite of their importance to the evolution of the early universe, their properties remain uncertain because of the limitations on previous numerical simulations and the lack of any observational constraints. Here, we investigate Population III star formation in five primordial halos using 3D radiation-hydrodynamical cosmological simulations. We find that multiple stars form in each minihalo and that their numbers increase over time, with up to 23 stars forming in one of the halos. Radiative feedback from the stars generates strong outflows, deforms the surrounding protostellar disk, and delays star formation for a few thousand years. Star formation rates vary with halo, and depend on the mass accretion onto the disk, the halo spin number, and the fraction of massive stars in the halo. The stellar masses in our models range from 0.1–37 M ⊙, and of the 55 stars that form in our models, 12 are >10 M ⊙ and most of the others are 1–10 M ⊙. Our…

Abstract The general behavior of the nuclear equation of state (EOS), relevant for the description of neutron stars (NSs), is studied within a Bayesian approach applied to a set of models based on a density-dependent relativistic mean-field description of nuclear matter. The EOS is subjected to a minimal number of constraints based on nuclear saturation properties and the low-density pure neutron matter EOS obtained from a precise next-to-next-to-next-to-leading order (N 3 LO) calculation in chiral effective field theory ( χ EFT). The posterior distributions of the model parameters obtained under these minimal constraints are employed to construct the distributions of various nuclear matter properties and NS properties such as radii, tidal deformabilities, central energy densities, speeds of sound, etc. We found that a 90% confidence interval for the allowed NS mass–radius relationship and tidal deformabilities is compatible with GW170817 and recent Neutron star Interior Composition ExploreR observations, without invoking the exotic degrees of freedom. A central speed of sound of…

Abstract We explore atmospheric escape from close-in exoplanets with the highest mass-loss rates. First, we locate the transition from stellar X-ray and UV-driven escape to rapid Roche lobe overflow, which occurs once the 10–100 nbar pressure level in the atmosphere reaches the Roche lobe. Planets enter this regime when the ratio of the substellar radius to the polar radius along the visible surface pressure level, which aligns with a surface of constant Roche potential, is X/Z ≳ 1.2 for Jovian planets ( M p ≳ 100 M ⊕ ) and X/Z ≳ 1.02 for sub-Jovian planets ( M p ≈ 10–100 M ⊕ ). Around a Sun-like star, this regime applies to orbital periods of less than two days for planets with radii of about 3–14R ⊕ . Our results agree with the properties of known transiting planets and can explain parts of the sub-Jovian desert in the population of known exoplanets. Second, we present detailed numerical simulations of atmospheric escape from a planet like Uranus or Neptune orbiting close to a Sun-like star that support the results above and point to interesting qualitative differences between hot…

Abstract We present a survey for photometric variability in young, low-mass brown dwarfs with the Spitzer Space Telescope. The 23 objects in our sample show robust signatures of youth and share properties with directly imaged exoplanets. We present three new young objects: 2MASS J03492367+0635078, 2MASS J09512690−8023553, and 2MASS J07180871−6415310. We detect variability in 13 young objects, and find that young brown dwarfs are highly likely to display variability across the L2–T4 spectral type range. In contrast, the field dwarf variability occurrence rate drops for spectral types >L9. We examine the variability amplitudes of young objects and find an enhancement in maximum amplitudes compared to field dwarfs. We speculate that the observed range of amplitudes within a spectral type may be influenced by secondary effects such as viewing inclination and/or rotation period. We combine our new rotation periods with the literature to investigate the effects of mass on angular momentum evolution. While high-mass brown dwarfs (>30 M Jup ) spin up over time, the same trend is not…

Abstract We present a comprehensive census of the active galactic nuclei (AGNs) in the GOODS-S/HUDF region from the X-ray to the radio, covering both the obscured and unobscured populations. This work includes a robust analysis of the source optical to mid-IR spectral energy distributions (SEDs) featuring (semi)empirical AGN and galaxy dust emission models and Bayesian fitting techniques, ultra-deep Very Large Array 3 and 6 GHz observations, and an integrated analysis of various AGN selection techniques, including X-ray properties, UV to mid-IR SED analysis, optical spectral features, mid-IR colors, radio-loudness and spectral slope, and AGN variability. In total, we report ∼900 AGNs over the ∼170 arcmin2 3D-HST GOODS-S footprint, which has doubled the AGN number identified in the previous X-ray sample, with ∼26% of our sample undetected in the deepest Chandra image. With a summary of AGN demographics from different selection methods, we find that no one single band or technique comes close to selecting a complete AGN sample despite the great depth of the data in GOODS-S/HUDF. We…

Abstract A 20 s cadence Transiting Exoplanet Survey Satellite monitoring campaign of 226 low-mass flare stars during Cycle 3 recorded 3792 stellar flares of ≥1032 erg. We explore the time-resolved emission and substructure in 440 of the largest flares observed at high signal-to-noise, 97% of which released energies of ≥1033 erg. We discover degeneracy present at 2 minute cadence between sharply peaked and weakly peaked flares is common, although 20 s cadence breaks these degeneracies. We better resolve the rise phases and find 46% of large flares exhibit substructure during the rise phase. We observe 49 candidate quasi-periodic pulsations (QPP) and confirm 17 at ≥3σ. Most of our QPPs have periods less than 10 minutes, suggesting short-period optical QPPs are common. We find QPPs in both the rise and decay phases of flares, including a rise-phase QPP in a large flare from Proxima Cen. We confirm that the Davenport et al. template provides a good fit to most classical flares observed at high cadence, although 9% favor Gaussian peaks instead. We characterize the properties of complex…

Abstract Gravitational-wave (GW) detections are starting to reveal features in the mass distribution of double compact objects. The lower end of the black hole (BH) mass distribution is especially interesting as few formation channels contribute here and because it is more robust against variations in the cosmic star formation than the high-mass end. In this work we explore the stable mass transfer channel for the formation of GW sources with a focus on the low-mass end of the mass distribution. We conduct an extensive exploration of the uncertain physical processes that impact this channel. We note that, for fiducial assumptions, this channel reproduces the peak at ∼9 M ☉ in the GW-observed binary BH mass distribution remarkably well and predicts a cutoff mass that coincides with the upper edge of the purported neutron star–black hole (NS–BH) mass gap. The peak and cutoff mass are a consequence of the unique properties of this channel; namely (1) the requirement of stability during the mass transfer phases, and (2) the complex way in which the final compact object masses scale with…

Abstract Interstellar dust extinction curves provide valuable information about dust properties, including the composition and size of the dust grains, and are essential to correct observations for the effects of interstellar dust. In this work, we measure a representative sample of near-infrared (NIR; 0.8–5.5 μ m) spectroscopic extinction curves for the first time, enabling us to investigate the extinction at wavelengths where it is usually only measured in broad photometric bands. We use IRTF/SpeX spectra of a sample of reddened and comparison stars to measure 15 extinction curves with the pair method. Our sample spans A ( V ) values from 0.78 to 5.65 and R ( V ) values from 2.43 to 5.33. We confirm that the NIR extinction curves are well fit by a power law, with indices and amplitudes differing from sight line to sight line. Our average diffuse NIR extinction curve can be represented by a single power law with index α = 1.7, but because of the sight line-to-sight line variations, the shape of any average curve will depend on the parental sample. We find that most of the variation…

Abstract We use the anelastic spherical harmonic code to model the convective dynamo of solar-type stars. Based on a series of 15 3D MHD simulations spanning four bins in rotation and mass, we show what mechanisms are at work in these stellar dynamos with and without magnetic cycles and how global stellar parameters affect the outcome. We also derive scaling laws for the differential rotation and magnetic field based on these simulations. We find a weaker trend between differential rotation and stellar rotation rate, ( Δ Ω ∝ ( ∣ Ω ∣ / Ω ⊙ ) 0.46 ) in the MHD solutions than in their HD counterpart ∣ Ω ∣ / Ω ⊙ 0.66 ), yielding a better agreement with the observational trends based on power laws. We find that for a fluid Rossby number between 0.15 ≲ Ro f ≲ 0.65, the solutions possess long magnetic cycle, if Ro f ≲ 0.42 a short cycle and if Ro f ≳ 1 (antisolar-like differential rotation), a statistically steady state. We show that short-cycle dynamos follow the classical Parker–Yoshimura rule whereas the long-cycle period ones do not. We also find efficient energy transfer between…

Abstract The extremely high brightness temperature of fast radio bursts (FRBs) requires that their emission mechanism must be “coherent,” either through concerted particle emission by bunches or through the exponential growth of a plasma wave mode or radiation amplitude via certain maser mechanisms. The bunching mechanism has been mostly discussed within the context of curvature radiation or cyclotron/synchrotron radiation. Here we propose a family of models invoking the coherent inverse Compton scattering (ICS) of bunched particles that may operate within or just outside of the magnetosphere of a flaring magnetar. Crustal oscillations during the flaring event may excite low-frequency electromagnetic waves near the magnetar surface. The X-mode of these waves could penetrate through the magnetosphere. Bunched relativistic particles in the charge-starved region inside the magnetosphere or in the current sheet outside the magnetosphere would upscatter these low-frequency waves to produce gigahertz emission to power FRBs. The ICS mechanism has a much larger emission power for individual…

Abstract We consider the flare prediction problem that distinguishes flare-imminent active regions that produce an M- or X-class flare in the succeeding 24 hr, from quiet active regions that do not produce any flares within ±24 hr. Using line-of-sight magnetograms and parameters of active regions in two data products covering Solar Cycles 23 and 24, we train and evaluate two deep learning algorithms—a convolutional neural network (CNN) and a long short-term memory (LSTM)—and their stacking ensembles. The decisions of CNN are explained using visual attribution methods. We have the following three main findings. (1) LSTM trained on data from two solar cycles achieves significantly higher true skill scores (TSSs) than that trained on data from a single solar cycle with a confidence level of at least 0.95. (2) On data from Solar Cycle 23, a stacking ensemble that combines predictions from LSTM and CNN using the TSS criterion achieves a significantly higher TSS than the “select-best” strategy with a confidence level of at least 0.95. (3) A visual attribution method called “integrated…

Abstract Prestellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures (<10 K) allow the formation of thick icy dust mantles, which will be partially preserved in future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in prestellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH 3 and, in particular, its deuterated isotopologues appear to maintain high abundances where CO molecules are mainly in the solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH 2 D freeze-out toward the L1544 prestellar core, suggestive of the presence of a “complete depletion zone” within a ≃1800 au radius, in agreement with astrochemical prestellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH 2 D becomes…

Abstract The Multi-slit Solar Explorer (MUSE) is a proposed mission composed of a multislit extreme ultraviolet (EUV) spectrograph (in three spectral bands around 171 Å, 284 Å, and 108 Å) and an EUV context imager (in two passbands around 195 Å and 304 Å). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (≤0.″5) and temporal resolution (down to ∼0.5 s for sit-and-stare observations), thanks to its innovative multislit design. By obtaining spectra in four bright EUV lines (Fe ix 171 Å, Fe xv 284 Å, Fe xix –Fe xxi 108 Å) covering a wide range of transition regions and coronal temperatures along 37 slits simultaneously, MUSE will, for the first time, “freeze” (at a cadence as short as 10 s) with a spectroscopic raster the evolution of the dynamic coronal plasma over a wide range of scales: from the spatial scales on which energy is released (≤0.″5) to the large-scale (∼170″ × 170″) atmospheric response. We use numerical modeling to showcase how MUSE will constrain the properties of the solar atmosphere on spatiotemporal scales (≤0.″5,…

Abstract We present a new derivation of the Milky Way’s current star formation rate (SFR) based on the data of the Herschel InfraRed Galactic Plane Survey (Hi-GAL). We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to 2.0 ± 0.7 M ⊙ yr −1 , of which 1.7 ± 0.6 M ⊙ yr −1 coming from clumps with reliable heliocentric distance assignment. This value is in general agreement with estimates found in the literature of last decades. The profile of SFR density averaged in Galactocentric rings is found to be qualitatively similar to others previously computed, with a peak corresponding to the Central Molecular Zone and another one around Galactocentric radius R gal ∼ 5 kpc, followed by an exponential decrease as log ( Σ SFR / [ M ⊙ yr − 1 kpc − 2 ] ) = a R gal / [ kpc ] + b , with a = −0.28 ± 0.01. In this regard, the fraction of SFR produced within and outside the solar circle is 84% and 16%, respectively; the…

Abstract We search for features in the mass distribution of detected compact binary coalescences which signify the transition between neutron stars (NSs) and black holes (BHs). We analyze all gravitational-wave (GW) detections by the LIGO Scientific Collaboration, the Virgo Collaboration, and the KAGRA Collaboration (LVK) made through the end of the first half of the third observing run, and find clear evidence for two different populations of compact objects based solely on GW data. We confidently (99.3%) find a steepening relative to a single power law describing NSs and low-mass BHs below 2.4 − 0.5 + 0.5 M ⊙ , which is consistent with many predictions for the maximum NS mass. We find suggestions of the purported lower mass gap between the most massive NSs and the least massive BHs, but are unable to conclusively resolve it with current data. If it exists, we find the lower mass gap’s edges to lie at 2.2 − 0.5 + 0.7 M ⊙ and 6.0 − 1.4 + 2.4 M ⊙ . We reexamine events that have been deemed “exceptional” by the LVK collaborations in the context of these features. We analyze GW190814…

Abstract We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10 m Keck telescope to search for Lyα emission from candidate galaxies at z ∼ 9–10 in four of the CANDELS fields (GOODS-N, EGS, UDS, and COSMOS). We observed 11 target galaxies, detecting Lyα from one object in ∼8.1 hr of integration, at z = 8.665 ± 0.001 with an integrated signal-to-noise ratio > 7. This galaxy is in the CANDELS Extended Groth Strip (EGS) field and lies physically close (3.5 physical Mpc [pMpc]) to another confirmed galaxy in this field with Lyα detected at z = 8.683. The detection of Lyα suggests the existence of large (∼1 pMpc) ionized bubbles fairly early in the reionization process. We explore the ionizing output needed to create bubbles of this size at this epoch and find that such a bubble requires more than the ionizing power provided by the full expected population of galaxies (by integrating the UV luminosity function down to M UV = −13). The Lyα we detect would be able to escape the predominantly neutral intergalactic medium at this epoch if our…

Abstract Motivated by the large observed diversity in the properties of extragalactic extinction by dust, we reanalyze the Cepheid calibration used to infer the Hubble constant, H 0, from Type Ia supernovae, using Cepheid data in 19 Type Ia supernova host galaxies from Riess et al. and anchor data from Riess et al. Unlike the SH0ES team, we do not enforce a fixed universal color–luminosity relation to correct the Cepheid magnitudes. Instead, we focus on a data-driven method, where the optical colors and near-infrared magnitudes of the Cepheids are used to derive individual color–luminosity relations for each Type Ia supernova host and anchor galaxy. We present two different analyses, one based on Wesenheit magnitudes, resulting in H 0 = 73.2 ± 1.3 km s−1 Mpc−1, a 4.2σ tension with the value inferred from the cosmic microwave background. In the second approach, we calibrate an individual extinction law for each galaxy, with noninformative priors using color excesses, yielding H 0 = 73.9 ± 1.8 km s−1 Mpc−1, in 3.4σ tension with the Planck value. Although the two methods yield similar…

Abstract The repeating FRB 20201124A was first discovered by CHIME/FRB in November of 2020, after which it was seen to repeat a few times over several months. It entered a period of high activity in April of 2021, at which time several observatories recorded tens to hundreds more bursts from the source. These follow-up observations enabled precise localization and host-galaxy identification. In this paper, we report on the CHIME/FRB-detected bursts from FRB 20201124A, including their best-fit morphologies, fluences, and arrival times. The large exposure time of the CHIME/FRB telescope toward the location of this source allows us to constrain its rates of activity. We analyze the repetition rates over different spans of time, constraining the rate prior to discovery to <3.4 day−1 (at 3σ), and demonstrate a significant change in the event rate following initial detection. Lastly, we perform a maximum-likelihood estimation of a power-law luminosity function, finding a best-fit index α = −4.6 ± 1.3 ± 0.6, with a break at a fluence threshold of F min ∼ 16.6 Jy ms, consistent with the…

Abstract Coherent curvature radiation as the radiation mechanism for fast radio bursts (FRBs) has been discussed since FRBs were discovered. We study the spectral and polarization properties of repeating FRBs within the framework of coherent curvature radiation by charged bunches in the magnetosphere of a highly magnetized neutron star. The spectra can be generally characterized by multisegmented broken power laws, and evolve as bunches move and the line of sight sweeps. Emitted waves are highly linear polarized and polarization angles are flat across the burst envelopes, if the line of sight is confined to the beam within an angle of 1/γ, while a circular polarization fraction becomes strong for off-beam cases. The spectro-temporal pulse-to-pulse properties can be a natural consequence due to the magnetospheric geometry. We investigate the relationship between drift rate, central frequency, and temporal duration. The radius-to-frequency mapping is derived and simulated within the assumptions of both dipolar and quadrupolar magnetic configurations. The geometric results show that…

Abstract We report the robust detection of coherent, localized deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disk around HD 163296. The analysis is performed using the discminer channel map modeling framework on 12CO J = 2–1 DSHARP data. Not only orbital radius but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at R = 94 ± 6 au, ϕ = 50° ± 3° (P94), is near the center of one of the gaps in dust continuum emission and is consistent with a planet mass of 1 M Jup. The other planet, located at R = 261 ± 4 au, ϕ = 57° ± 1° (P261), is in the region where a velocity kink was previously observed in 12CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disk and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles, we detect gas gaps at R = 38, 88, and 136 au, closely matching the location of their dust and kinematical…

Abstract We measure homogeneous distances to M31 and 38 associated stellar systems (−16.8 ≤ M V ≤ −6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From >700 orbits of new/archival Advanced Camera for Surveys imaging, we identify >4700 RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 day and 0.04 mag. Based on period–Wesenheit–metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of ∼20 kpc (3%) and ∼10 kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that ∼80% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms 7–23 kpc) planar “arc” of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess the physical proximity of notable…

Abstract We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with Cosmic Origins Spectrograph (COS) data. The outflows are traced by blueshifted absorption lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static interstellar medium (ISM) from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow’s mean velocity and velocity width, and find that both correlate in a highly significant way with the star formation rate, galaxy mass, and circular velocity over ranges of four orders of magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon…

Abstract We perform joint modeling of the composite rest-frame far-UV and optical spectra of redshift 1.85 ≤ z ≤ 3.49 star-forming galaxies to deduce key properties of the massive stars, ionized interstellar medium (ISM), and neutral ISM, with the aim of investigating the principal factors affecting the production and escape of Lyα photons. Our sample consists of 136 galaxies with deep Keck/LRIS and MOSFIRE spectra covering, respectively, Lyβ through C iii] λλ1907, 1909 and [O ii], [Ne iii], Hβ, [O iii], Hα, [N ii], and [S ii]. Spectral and photoionization modeling indicates that the galaxies are uniformly consistent with stellar population synthesis models that include the effects of stellar binarity. Over the dynamic range of our sample, there is little variation in stellar and nebular abundance with Lyα equivalent width, W λ (Lyα), and only a marginal anticorrelation between age and W λ (Lyα). The inferred range of ionizing spectral shapes is insufficient to solely account for the variation in W λ (Lyα); rather, the covering fraction of optically thick H i appears to be the…

Abstract The Milky Way (MW) stellar halo contains relics of ancient mergers that tell the story of our galaxy’s formation. Some of them are identified due to their similarity in energy, actions, and chemistry, referred to as the “chemodynamical space,” and are often attributed to distinct merger events. It is also known that our galaxy went through a significant merger event that shaped the local stellar halo during its first billion years. Previous studies using N-body only and cosmological hydrodynamical simulations have shown that such a single massive merger can produce several “signatures” in the chemodynamical space, which can potentially be misinterpreted as distinct merger events. Motivated by these, in this work we use a subset of the GASTRO library, which consists of several smoothed particle hydrodynamics+N-body models of a single accretion event in a MW-like galaxy. Here, we study models with orbital properties similar to the main merger event of our galaxy and explore the implications to known stellar halo substructures. We find that (i) supernova feedback efficiency…

Abstract Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully coherent search for such signals from eighteen pulsars in data from LIGO and Virgo’s third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow both the frequency and the time derivative of the frequency of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search, we look in O3 data for long-duration (hours–months) transient gravitational waves in the…

Abstract We evolve two high-resolution general relativistic magnetohydrodynamic simulations of advection-dominated accretion flows around nonspinning black holes (BHs), each over a duration ∼3 × 10 5 GM BH / c 3 . One model captures the evolution of a weakly magnetized (SANE) disk and the other that of a magnetically arrested disk (MAD). Magnetic flux eruptions in the MAD model push out gas from the disk and launch strong winds with outflow efficiencies at times reaching 10% of the incoming accretion power. Despite the substantial power in these winds, average mass outflow rates remain low out to a radius ∼100 GM BH / c 2 , only reaching ∼60%–80% of the horizon accretion rate. The average outward angular momentum transport is primarily radial in both modes of accretion, but with a clear distinction: magnetic flux eruption–driven disk winds cause a strong vertical flow of angular momentum in the MAD model, while for the SANE model, the magnetorotational instability (MRI) moves angular momentum mostly equatorially through the disk. Further, we find that the MAD state is highly…

Abstract The Kepler and Transiting Exoplanet Survey Satellite (TESS) missions have generated over 100,000 potential transit signals that must be processed in order to create a catalog of planet candidates. During the past few years, there has been a growing interest in using machine learning to analyze these data in search of new exoplanets. Different from the existing machine learning works, ExoMiner , the proposed deep learning classifier in this work, mimics how domain experts examine diagnostic tests to vet a transit signal. ExoMiner is a highly accurate, explainable, and robust classifier that (1) allows us to validate 301 new exoplanets from the MAST Kepler Archive and (2) is general enough to be applied across missions such as the ongoing TESS mission. We perform an extensive experimental study to verify that ExoMiner is more reliable and accurate than the existing transit signal classifiers in terms of different classification and ranking metrics. For example, for a fixed precision value of 99%, ExoMiner retrieves 93.6% of all exoplanets in the test set (i.e., recall =…

Abstract In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolving approximate waveform degeneracies, which may result in multimodal posterior distributions. The family of waveform models we use includes a new fast time-domain model (IMRPhenomTPHM), which allows us to extensively test different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper, Nitz & Capano, associated with possible support for a high-mass-ratio source, but confirm their finding of a multimodal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the pair…

Abstract Quasiperiodic erupters are a remarkable class of objects exhibiting very-large-amplitude quasiperiodic X-ray flares. Although numerous dynamical models have been proposed to explain them, relatively little attention has been given to using the properties of their radiation to constrain their dynamics. Here we show that the observed luminosity, spectrum, repetition period, duty cycle, and fluctuations in the latter two quantities point toward a model in which a main-sequence star on a moderately eccentric orbit around a supermassive black hole periodically transfers mass to the Roche lobe of the black hole; orbital dynamics lead to mildly relativistic shocks near the black hole; and thermal X-rays at the observed temperature are emitted by the gas as it flows away from the shock. Strong X-ray irradiation of the star by the flare itself augments the mass transfer, creates fluctuations in flare timing, and stirs turbulence in the stellar atmosphere that amplifies the magnetic field to a level at which magnetic stresses can accelerate infall of the transferred mass toward the…

Abstract Direct imaging studies have mainly used low-resolution spectroscopy (R ∼ 20–100) to study the atmospheres of giant exoplanets and brown dwarf companions, but the presence of clouds has often led to degeneracies in the retrieved atmospheric abundances (e.g., carbon-to-oxygen ratio, metallicity). This precludes clear insights into the formation mechanisms of these companions. The Keck Planet Imager and Characterizer (KPIC) uses adaptive optics and single-mode fibers to transport light into NIRSPEC (R ∼ 35,000 in the K band), and aims to address these challenges with high-resolution spectroscopy. Using an atmospheric retrieval framework based on petitRADTRANS, we analyze the KPIC high-resolution spectrum (2.29–2.49 μm) and the archival low-resolution spectrum (1–2.2 μm) of the benchmark brown dwarf HD 4747 B (m = 67.2 ± 1.8 M Jup, a = 10.0 ± 0.2 au, T eff ≈ 1400 K). We find that our measured C/O and metallicity for the companion from the KPIC high-resolution spectrum agree with those of its host star within 1σ–2σ. The retrieved parameters from the K-band high-resolution…

Abstract The distribution of the dark matter (DM) in DM-admixed neutron stars (DANSs) is supposed to result in either a dense dark core or an extended dark halo, subject to the DM fraction of the DANS (f χ ) and the DM properties, such as the mass (m χ ) and the strength of the self-interaction (y). In this paper, we perform an in-depth analysis of the formation criterion for dark cores/dark halos, and point out that the relative distribution of these two components is essentially determined by the ratio of the central enthalpy of the DM component to that of the baryonic matter component inside the DANSs. For the critical case where the radii of the DM and the baryonic matter are the same, we further derive an analytical formula to describe the dependence of f χ crit on m χ and y for a given DANS mass. The relative distribution of the two components in DANSs can lead to different observational effects. We here focus on the modification of the pulsar pulse profile, due to the extra light-bending effect in the case of a dark halo existence, and conduct the first investigation into the…

Abstract The recent NICER measurement of the radius of the neutron star PSR J0740+6620, and the inferred small variation of radii from 1.4 to 2.1 M ⊙ , reveal key features of the equation of state of neutron star matter. The pressure rises rapidly in the regime of baryon density n ∼ 2–4 times nuclear saturation density, n 0 —the region where we expect hadronic matter to be undergoing transformation into quark matter—and the pressure in the nuclear regime is greater than predicted by microscopic many-body variational calculations of nuclear matter. To incorporate these insights into the microscopic physics from the nuclear to the quark matter regimes, we construct an equation of state, QHC21, within the framework of quark–hadron crossover. We include nuclear matter results primarily based on the state-of-the-art chiral effective field theory, but also note results of using nuclear matter variational calculations based on empirical nuclear forces. We employ explicit nuclear degrees of freedom only up to n ∼ 1.5 n 0 , in order to explore the possibility of further physical degrees of…

Abstract We explore the three-dimensional properties of convective, luminous (L ≈ 104.5–105 L ⊙), hydrogen-rich envelopes of red supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using Athena++. These computations comprise ≈30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈3M ⊙ of the simulated M ≈ 15M ⊙ stars. This work reveals a radius, R corr, around which the nature of the convection changes. For r > R corr, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well studied in less luminous stars, but in RSGs, the near- (or above-)Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outward rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This “halo” of material will impact predictions for both shock breakout and…