Showing 1–50 of 288 results for author: Conroy, C

Modeling the Ages and Chemical Abundances of Elliptical Galaxies

Authors: Nicole Marcelina Gountanis, David H. Weinberg, Aliza G. Beverage, Nathan R. Sandford, Charlie Conroy, Mariska Kriek

Abstract: Spectroscopic studies of elliptical galaxies show that their stellar population ages, mean metallicity, and $α$-enhancement traced by [Mg/Fe] all increase with galaxy stellar mass or velocity dispersion. We use one-zone galactic chemical evolution (GCE) models with a flexible star formation history (SFH) to model the age, [Mg/H], and [Mg/Fe] inferred from simple stellar population (SSP) fits to ob… ▽ More Spectroscopic studies of elliptical galaxies show that their stellar population ages, mean metallicity, and $α$-enhancement traced by [Mg/Fe] all increase with galaxy stellar mass or velocity dispersion. We use one-zone galactic chemical evolution (GCE) models with a flexible star formation history (SFH) to model the age, [Mg/H], and [Mg/Fe] inferred from simple stellar population (SSP) fits to observed ellipticals at $z \sim 0$ and $z \sim 0.7$. We show that an SSP fit to the spectrum computed from a full GCE model gives ages and abundances close to the light-weighted, logarithmically averaged values of the composite stellar population, <age>, <[Mg/H]>, and <[Mg/Fe]>. With supernova Mg and Fe yields fixed to values motivated by Milky Way stellar populations, we find that predicted <[Mg/H]>-<age> and <[Mg/Fe]>-<age> relations are surprisingly insensitive to SFH parameters: older galaxies have higher <[Mg/Fe]>, but the detailed form of the SFH has limited impact. The star formation efficiency and outflow efficiency affect the early and late evolution of <[Mg/H]>, respectively; explaining observed trends requires higher star formation efficiency and lower outflows in more massive galaxies. With core collapse supernova yields calibrated to the plateau [Mg/Fe]$_{\rm cc} \approx0.45$ observed in many Milky Way studies, our models underpredict the observed <[Mg/Fe]> ratios of ellipticals by 0.05-0.1 dex. Increasing the core collapse yield ratio to [Mg/Fe]$_{\rm cc} = 0.55$ improves the agreement, though the models still lie below the data. We discuss potential resolutions of this discrepancy, including the possibility that many ellipticals terminate their star formation with a self-enriching, terminating burst that reduces the light-weighted age and boosts <[Mg/Fe]>. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Reconciling M/L Ratios Across Cosmic Time: a Concordance IMF for Massive Galaxies

Authors: Pieter van Dokkum, Charlie Conroy

Abstract: The stellar initial mass function (IMF) is thought to be bottom-heavy in the cores of the most massive galaxies, with an excess of low mass stars compared to the Milky Way. However, studies of the kinematics of quiescent galaxies at 2<z<5 find M/L ratios that are better explained with lighter IMFs. Light IMFs have also been proposed for the unexpected populations of luminous galaxies that JWST has… ▽ More The stellar initial mass function (IMF) is thought to be bottom-heavy in the cores of the most massive galaxies, with an excess of low mass stars compared to the Milky Way. However, studies of the kinematics of quiescent galaxies at 2<z<5 find M/L ratios that are better explained with lighter IMFs. Light IMFs have also been proposed for the unexpected populations of luminous galaxies that JWST has uncovered at z>7, to reduce tensions with galaxy formation models. Here we explore IMFs that are simultaneously bottom-heavy, with a steep slope at low stellar masses, and top-heavy, with a shallow slope at high masses. We derive a form of the IMF for massive galaxies that is consistent with measurements in the local universe and yet produces relatively low M/L ratios at high redshift. This `concordance' IMF has slopes $γ_1=2.40\pm0.09$, $γ_2=2.00\pm0.14$, and $γ_3=1.85\pm0.11$ in the regimes 0.1-0.5 Msun, 0.5-1 Msun, and >1 Msun respectively. The IMF parameter $α$, the mass excess compared to a Milky Way IMF, ranges from $\log(α)\approx+0.3$ for present-day galaxies to $\log(α)\approx-0.1$ for their star forming progenitors. The concordance IMF applies only to the central regions of the most massive galaxies, with velocity dispersions ~300 km/s, and their progenitors. However, it can be generalized using a previously-measured relation between $α$ and $σ$. We arrive at the following modification to the Kroupa (2001) IMF for galaxies with $σ\gtrsim 160$ km/s: $γ_1\approx1.3+4.3\logσ_{160}$; $γ_2\approx2.3-1.2\logσ_{160}$; and $γ_3\approx2.3-1.7\logσ_{160}$, with $σ_{160}=σ/160$ km/s. If galaxies grow primarily inside-out, so that velocity dispersions are relatively stable, these relations should also hold at high redshift. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ Letters. Figure 3 shows the proposed IMF for massive galaxies. For those not liking 3D figures: the information in Fig. 2 is mostly illustrative

Metals in Star-Forming Galaxies with KCWI. I. Methodology and First Results on the Abundances of Iron, Magnesium, and Oxygen

Authors: Zhuyun Zhuang, Evan N. Kirby, Charles C. Steidel, Mithi A. C. de los Reyes, Nikolaus Z. Prusinski, N. Leethochawalit, Minjung Park, Charlie Conroy, Evan H. Nuñez

Abstract: Understanding the chemical enrichment of different elements is crucial to gaining a complete picture of galaxy chemical evolution. In this study, we present a new sample of 46 low-redshift, low-mass star-forming galaxies at $M_*\sim 10^{8-10}M_{\odot}$ along with two quiescent galaxies at $M_*\sim 10^{8.8}M_{\odot}$ observed with the Keck Cosmic Web Imager (KCWI), aiming to investigate the chemica… ▽ More Understanding the chemical enrichment of different elements is crucial to gaining a complete picture of galaxy chemical evolution. In this study, we present a new sample of 46 low-redshift, low-mass star-forming galaxies at $M_*\sim 10^{8-10}M_{\odot}$ along with two quiescent galaxies at $M_*\sim 10^{8.8}M_{\odot}$ observed with the Keck Cosmic Web Imager (KCWI), aiming to investigate the chemical evolution of galaxies in the transition zone between Local Group satellites and massive field galaxies. We develop a novel method to simultaneously determine stellar abundances of iron and magnesium in star-forming galaxies. With the gas-phase oxygen abundance (O/H)$_{\rm g}$ measured using the strong line method, we are able to make the first-ever apples-to-apples comparison of $α$ elements in the stars and the ISM. We find that the [Mg/H]$_*$-[O/H]$_{\rm g}$ relation is much tighter than the [Fe/H]$_*$-[O/H]$_{\rm g}$ relation, which can be explained by the similar production processes of $α$ elements. Most galaxies in our sample exhibit higher [O/H]$_{\rm g}$ than [Fe/H]$_*$ and [Mg/H]$_*$. In addition, we construct mass-metallicity relations (MZRs) measured as three different elements (Fe$_*$, Mg$_*$, O$_{\rm g}$). Compared to the gas O-MZR, the stellar Fe- and Mg-MZRs show larger scatter driven by variations in specific star formation rates (sSFR), with star-forming galaxies exhibiting higher sSFR and lower stellar abundances at fixed mass. The excess of [O/H]$_{\rm g}$ compared to stellar abundances as well as the anti-correlation between sSFR and stellar abundance suggests that galaxy quenching of intermediate-mass galaxies at $M_*\sim 10^{8-10}M_{\odot}$ is primarily driven by starvation. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 36 pages, 16 figures, 4 tables; accepted for publication in ApJ. Main results in Figure 7, 10 and 11

Carbon and Iron Deficiencies in Quiescent Galaxies at z=1-3 from JWST-SUSPENSE: Implications for the Formation Histories of Massive Galaxies

Authors: Aliza G. Beverage, Martje Slob, Mariska Kriek, Charlie Conroy, Guillermo Barro, Rachel Bezanson, Gabriel Brammer, Chloe M. Cheng, Anna de Graaff, Natascha M. Förster Schreiber, Marijn Franx, Brian Lorenz, Pavel E. Mancera Piña, Danilo Marchesini, Adam Muzzin, Andrew B. Newman, Sedona H. Price, Alice E. Shapley, Mauro Stefanon, Katherine A. Suess, Pieter van Dokkum, David Weinberg, Daniel R. Weisz

Abstract: We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log M/M$_\odot$=10.2-11.2) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 dex in [Mg/H], implying rapid formation and quenching befor… ▽ More We present the stellar metallicities and multi-element abundances (C, Mg, Si, Ca, Ti, Cr, and Fe) of 15 massive (log M/M$_\odot$=10.2-11.2) quiescent galaxies at z=1-3, derived from ultradeep JWST-SUSPENSE spectra. Compared to quiescent galaxies at z~0, these galaxies exhibit a deficiency of 0.25 dex in [C/H], 0.16 dex in [Fe/H], and 0.07 dex in [Mg/H], implying rapid formation and quenching before significant enrichment from asymptotic giant branch stars and Type Ia supernovae. Additionally, we find that galaxies that form at higher redshift have higher [Mg/Fe] and lower [Fe/H] and [Mg/H], irrespective of their observed redshift. The evolution in [Fe/H] and [C/H] is therefore primarily explained by lower redshift samples naturally including galaxies with longer star-formation timescales. On the other hand, the lower [Mg/H] can be explained by galaxies forming at earlier epochs expelling larger gas reservoirs during their quenching phase. Consequently, the mass-metallicity relation, primarily reflecting [Mg/H], is also lower at z=1-3 compared to the lower redshift relation, though the slopes are similar. Finally, we compare our results to standard stellar population modeling approaches employing solar abundance patterns and non-parametric star-formation histories (using Prospector). Our SSP-equivalent ages agree with the mass-weighted ages from Prospector, while the metallicities disagree significantly. Nonetheless, the metallicities better reflect [Fe/H] than total [Z/H]. We also find that star-formation timescales inferred from elemental abundances are significantly shorter than those from Prospector, and we discuss the resulting implications for the early formation of massive galaxies. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ; 18 pages, 6 figures, 1 table

Our Halo of Ice and Fire: Strong Kinematic Asymmetries in the Galactic Halo

Authors: Jiwon Jesse Han, Charlie Conroy, Dennis Zaritsky, Ana Bonaca, Nelson Caldwell, Vedant Chandra, Yuan-Sen Ting

Abstract: The kinematics of the stellar halo hold important clues to the assembly history and mass distribution of the Galaxy. In this study, we map the kinematics of stars across the Galactic halo with the H3 Survey. We find a complex distribution that breaks both azimuthal symmetry about the $Z$-axis and mirror symmetry about the Galactic plane. This asymmetry manifests as large variations in the radial v… ▽ More The kinematics of the stellar halo hold important clues to the assembly history and mass distribution of the Galaxy. In this study, we map the kinematics of stars across the Galactic halo with the H3 Survey. We find a complex distribution that breaks both azimuthal symmetry about the $Z$-axis and mirror symmetry about the Galactic plane. This asymmetry manifests as large variations in the radial velocity dispersion $σ_r$ from as ``cold'' as 70 $\text{km}\text{ s}^{-1}$ to as ``hot'' as 160 $\text{km}\text{ s}^{-1}$. We use stellar chemistry to distinguish accreted stars from in-situ stars in the halo, and find that the accreted population has higher $σ_r$ and radially biased orbits, while the in-situ population has lower $σ_r$ and isotropic orbits. As a result, the Galactic halo kinematics are highly heterogeneous and poorly approximated as being spherical or axisymmetric. We measure radial profiles of $σ_r$ and the anisotropy parameter $β$ over Galactocentric radii $10-80\text{ kpc}$, and find that discrepancies in the literature are due to the nonspherical geometry and heterogeneous nature of the halo. Investigating the effect of strongly asymmetric $σ_r$ and $β$ on equilibrium models is a path forward to accurately constraining the Galactic gravitational field, including its total mass. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: submitted to ApJ; comments welcome

AGN Feedback in Quiescent Galaxies at Cosmic Noon Traced by Ionized Gas Emission

Authors: Letizia Bugiani, Sirio Belli, Minjung Park, Rebecca L. Davies, J. Trevor Mendel, Benjamin D. Johnson, Amir H. Khoram, Chloë Benton, Andrea Cimatti, Charlie Conroy, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7<z<3.5$ observed with JWST/NIRSpec by the Blue Jay survey. Robust detection of emission lines in $75\%$ of the sample indicates the presence of ongoing ionizing sources in this passive population. The H$α$ line luminosities confirm that the population is quiescent, with star formation… ▽ More We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7<z<3.5$ observed with JWST/NIRSpec by the Blue Jay survey. Robust detection of emission lines in $75\%$ of the sample indicates the presence of ongoing ionizing sources in this passive population. The H$α$ line luminosities confirm that the population is quiescent, with star formation rates that are at least ten times lower than the main sequence of star formation. The quiescent sample is clearly separate from the star-forming population in line diagnostic diagrams, and occupies a region usually populated by active galactic nuclei (AGN). Analysis of the observed line ratios, equivalent widths, and velocity dispersions leads us to conclude that in most cases the gas is ionized by AGN activity, despite the lack of X-ray detections. A subset of the sample also hosts ionized and/or neutral outflows. Our results show, for the first time using a representative sample, that low luminosity AGN are extremely common among quiescent galaxies at high redshift. These low luminosity AGN may play a key role in quenching star formation and in maintaining massive galaxies quiescent from Cosmic Noon to $z\sim0$. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 20 pages, 18 figures. Submitted to ApJ

On the Origin of High Velocity Clouds in the Galaxy

Authors: Scott Lucchini, Jiwon Jesse Han, Lars Hernquist, Charlie Conroy

Abstract: The origin of our Galaxy's high-velocity clouds (HVCs) remains a mystery after many decades of effort. In this paper, we use the TNG50 simulation of the IllustrisTNG project to identify cool, dense clouds that match observations of Galactic HVCs. We track these clouds back in time to determine their origin. We find that only 17% of HVCs (by mass) can be tracked directly to the disk, and 21% to sat… ▽ More The origin of our Galaxy's high-velocity clouds (HVCs) remains a mystery after many decades of effort. In this paper, we use the TNG50 simulation of the IllustrisTNG project to identify cool, dense clouds that match observations of Galactic HVCs. We track these clouds back in time to determine their origin. We find that only 17% of HVCs (by mass) can be tracked directly to the disk, and 21% to satellite accretion. The majority of HVCs (62%) arise from warm and hot circumgalactic gas that cools through thermal instability. They then obtain their anomalous velocities through interactions with the turbulent circumgalactic medium. At TNG50 resolution, we do not see evidence for HVCs forming out of very low metallicity intergalactic material. Instead, low metallicity HVCs are most likely associated with satellite accretion. These results suggest that Galactic HVCs are highly heterogeneous in their origin, and can provide insight into the physical processes that shape the circumgalactic medium such as disk outflows, satellite accretion, and thermal instabilities. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 15 pages, 12 figures, submitted. Comments welcome

All-Sky Kinematics of the Distant Halo: The Reflex Response to the LMC

Authors: Vedant Chandra, Rohan P. Naidu, Charlie Conroy, Nicolas Garavito-Camargo, Chervin Laporte, Ana Bonaca, Phillip A. Cargile, Emily Cunningham, Jiwon Jesse Han, Benjamin D. Johnson, Hans-Walter Rix, Yuan-Sen Ting, Turner Woody, Dennis Zaritsky

Abstract: The infall of the Large Magellanic Cloud (LMC) is predicted to displace the inner Milky Way (MW), imprinting an apparent 'reflex motion' on the observed velocities of distant halo stars. We construct the largest all-sky spectroscopic dataset of luminous red giant stars from $50-160$ kpc, including a new survey of the southern celestial hemisphere. We fit the full 6D kinematics of our data to measu… ▽ More The infall of the Large Magellanic Cloud (LMC) is predicted to displace the inner Milky Way (MW), imprinting an apparent 'reflex motion' on the observed velocities of distant halo stars. We construct the largest all-sky spectroscopic dataset of luminous red giant stars from $50-160$ kpc, including a new survey of the southern celestial hemisphere. We fit the full 6D kinematics of our data to measure the amplitude and direction of the inner MW's motion towards the outer halo. The observed velocity grows with distance such that, relative to halo stars at $100$ kpc, the inner MW is lurching at $\approx 40$ km s$^{-1}$ towards a recent location along the LMC's past orbit. Our measurements align with N-body simulations of the halo's response to a $1.8 \times 10^{11} M_\odot$ LMC on first infall, suggesting that the LMC is at least 15% as massive as the MW. Our findings highlight the dramatic disequilibrium of the MW outskirts, and will enable more accurate measurements of the total mass of our Galaxy. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 25 pages, 15 figures. Submitted to ApJ

A Model for Eruptive Mass Loss in Massive Stars

Authors: Shelley J. Cheng, Jared A. Goldberg, Matteo Cantiello, Evan B. Bauer, Mathieu Renzo, Charlie Conroy

Abstract: Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well-understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in the MESA stellar evolution code. The model identifies regions in the star where the energy associated wit… ▽ More Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well-understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in the MESA stellar evolution code. The model identifies regions in the star where the energy associated with a locally super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with masses $10-100~M_\odot$ at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to $\sim10^{-2}~M_\odot \mathrm{yr}^{-1}$ can be driven by this mechanism, and occurs in a vertical band on the HR diagram between $3.5 \lesssim \log(T_\mathrm{eff}/\mathrm{K}) \lesssim 4.0$. This predicted eruptive mass loss prevents stars of initial masses $\gtrsim20~M_\odot$ from evolving to become red supergiants, with the stars instead ending their lives as blue supergiants, and therefore offers a possible explanation for the observed lack of red supergiants in that mass regime. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 16 pages, 9 figures, submitted to ApJ

Widespread rapid quenching at cosmic noon revealed by JWST deep spectroscopy

Authors: Minjung Park, Sirio Belli, Charlie Conroy, Benjamin D. Johnson, Rebecca L. Davies, Joel Leja, Sandro Tacchella, J. Trevor Mendel, Chloë Benton, Letizia Bugiani, Razieh Emami, Amirhossein Khoram, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Bryan A. Terrazas, Rainer Weinberger

Abstract: Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation historie… ▽ More Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation histories by fitting spectral energy distribution models to the JWST/NIRSpec $R\sim1000$ spectra. We find that massive quiescent galaxies can be split into three categories with roughly equal numbers of galaxies according to their SFHs: 1) Relatively old galaxies quenched at early epochs; 2) Galaxies that are rapidly and recently quenched after a flat or bursty formation history (depending on the assumed prior); 3) Galaxies that are rapidly and recently quenched after a major starburst. Most recently quenched galaxies show neutral gas outflows, probed by blueshifted $\rm Na\,I\,D$ absorption, and ionized gas emission, with line ratios consistent with active galactic nucleus (AGN) diagnostics. This suggests that AGN activity drives multi-phase gas outflows, leading to rapid quenching. By tracing back the SFHs of the entire sample, we predict the number density of massive quiescent galaxies at $z=4-6$: $n=3.0\pm1.4\times10^{-5}\,\rm Mpc^{-3}$. The two oldest massive quiescent galaxies in our sample appear to have extremely early formation and quenching ($z\gtrsim6$), possibly descendants of early post-starbursts at $z>3$. These galaxies still show neutral gas reservoirs and low-level star formation, consistent with weak H$α$ emission, perhaps because the ejective AGN feedback that caused rapid quenching has weakened over time. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ; 21 pages, 9 figures, 3 tables (excluding appendix)

The JWST-SUSPENSE Ultradeep Spectroscopic Program: Survey Overview and Star-Formation Histories of Quiescent Galaxies at 1 < z < 3

Authors: Martje Slob, Mariska Kriek, Aliza G. Beverage, Katherine A. Suess, Guillermo Barro, Rachel Bezanson, Gabriel Brammer, Chloe M. Cheng, Charlie Conroy, Anna de Graaff, Natascha M. Förster Schreiber, Marijn Franx, Brian Lorenz, Pavel E. Mancera Piña, Danilo Marchesini, Adam Muzzin, Andrew B. Newman, Sedona H. Price, Alice E. Shapley, Mauro Stefanon, Pieter van Dokkum, Daniel R. Weisz

Abstract: We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quie… ▽ More We present an overview and first results from the Spectroscopic Ultradeep Survey Probing Extragalactic Near-infrared Stellar Emission (SUSPENSE), executed with NIRSpec on JWST. The primary goal of the SUSPENSE program is to characterize the stellar, chemical, and kinematic properties of massive quiescent galaxies at cosmic noon. In a single deep NIRSpec/MSA configuration, we target 20 distant quiescent galaxy candidates ($z=1-3$, $H_{AB}\le23$), as well as 53 star-forming galaxies at $z=1-4$. With 16~hr of integration and the G140M-F100LP dispersion-filter combination, we observe numerous Balmer and metal absorption lines for all quiescent candidates. We derive stellar masses (log$M_*/M_{\odot}\sim10.2-11.5$) and detailed star-formation histories (SFHs) and show that all 20 candidate quiescent galaxies indeed have quenched stellar populations. These galaxies show a variety of mass-weighted ages ($0.8-3.3$~Gyr) and star formation timescales ($\sim0.5-4$~Gyr), and four out of 20 galaxies were already quiescent by $z=3$. On average, the $z>1.75$ $[z<1.75]$ galaxies formed 50\% of their stellar mass before $z=4$ $[z=3]$. Furthermore, the typical SFHs of galaxies in these two redshift bins ($z_{\text{mean}}=2.2~[1.3]$) indicate that galaxies at higher redshift formed earlier and over shorter star-formation timescales compared to lower redshifts. Although this evolution is naturally explained by the growth of the quiescent galaxy population over cosmic time, number density calculations imply that mergers and/or late-time star formation also contribute to the evolution. In future work, we will further unravel the early formation, quenching, and late-time evolution of these galaxies by extending this work with studies on their chemical abundances, resolved stellar populations and kinematics. △ Less

Submitted 18 July, 2024; v1 submitted 18 April, 2024; originally announced April 2024.

Comments: Accepted in ApJ; 25 pages, 14 figures, 2 tables (excluding appendices)

Anatomy of an ionized bubble: NIRCam grism spectroscopy of the $z=6.6$ double-peaked Lyman-$α$ emitter COLA1 and its environment

Authors: Alberto Torralba-Torregrosa, Jorryt Matthee, Rohan P. Naidu, Ruari Mackenzie, Gabriele Pezzulli, Anne Hutter, Pablo Arnalte-Mur, Siddhartha Gurung-López, Sandro Tacchella, Pascal Oesch, Daichi Kashino, Charlie Conroy, David Sobral

Abstract: The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$α$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$α$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected thro… ▽ More The increasingly neutral intergalactic gas at $z>6$ impacts the Lyman-$α$ flux observed from galaxies. One luminous galaxy, COLA1, stands out because of its unique double-peaked Ly$α$ line at $z=6.6$, unseen in any simulation of reionization. Here we present JWST/NIRCam wide-field slitless spectroscopy in a 21 arcmin$^2$ field centered on COLA1. We find 141 galaxies spectroscopically-selected through the [OIII]($\lambda4969,5008$) doublet at $5.35<z<6.95$, with 40 of these sources showing H$β$. For COLA1 we additionally detect [OIII]$_{4363}$ and H$γ$. We measure a systemic redshift of $z=6.5917$ for COLA1, confirming the double-peak nature of the Ly$α$ profile. This implies that it resides in a highly ionized bubble and that it is leaking ionizing photons with a high escape fraction $f_{\rm esc}{\rm (LyC)}=20$-$50$%, making it a prime laboratory to study Lyman continuum escape in the Epoch of Reionization. COLA1 shows all the signs of a prolific ionizer with a Ly$α$ escape fraction of $81\pm5\%$, Balmer decrement indicating no dust, a steep UV slope ($β_{\rm UV}=-3.2\pm 0.4$), and a star-formation surface density $\gtrsim 10\times$ that of typical galaxies at similar redshift. We detect 5 galaxies in COLA1's close environment ($Δz<0.02$). Exploiting the high spectroscopic completeness inherent to grism surveys, and using mock simulations that mimic the selection function, we show the that number of detected companions is very typical for a similarly UV-bright ($M_{\rm{UV}}\sim-21.3$) galaxy; that is, the ionized bubble around COLA1 is unlikely due to an excessively large over-density. Instead, the measured ionizing properties suggest that COLA1 by itself might be powering the bubble required to explain its double-peaked Ly$α$ profile ($R_{\rm ion}\approx0.7$ pMpc), with minor contribution from detected neighbours ($-17.5>M_{\rm UV}>-19.5$). △ Less

Submitted 26 June, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: 15 pages, 13 figures. Accepted for publication in A&A. Main figs. 1, 2, 10 (SFR surface density), 11 (spectra and positions of galaxies in the bubble). Corrected error in escape fraction estimate from SFR surface density

The JWST Resolved Stellar Populations Early Release Science Program V. DOLPHOT Stellar Photometry for NIRCam and NIRISS

Authors: Daniel R. Weisz, Andrew E. Dolphin, Alessandro Savino, Kristen B. W. McQuinn, Max J. B. Newman, Benjamin F. Williams, Nitya Kallivayalil, Jay Anderson, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Karin M. Sandstrom, Andrew A. Cole, Jack T. Warfield, Evan D. Skillman, Roger E. Cohen, Rachael Beaton, Alessandro Bressan, Alberto Bolatto, Michael Boylan-Kolchin, Alyson M. Brooks, James S. Bullock, Charlie Conroy, Michael C. Cooper, Julianne J. Dalcanton , et al. (16 additional authors not shown)

Abstract: We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy),… ▽ More We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT's photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT's photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ($\lesssim0.05$ mag) chip-to-chip variations in NIRCam's zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 30 pages, 17 figures. Accepted to ApJS. Data products to be hosted on MAST. For DOLPHOT/JWST tutorials, see https://dolphot-jwst.readthedocs.io/en/latest/ . For more program and DOLPHOT info, see https://ers-stars.github.io

Spectacular nucleosynthesis from early massive stars

Authors: Alexander P. Ji, Sanjana Curtis, Nicholas Storm, Vedant Chandra, Kevin C. Schlaufman, Keivan G. Stassun, Alexander Heger, Marco Pignatari, Adrian M. Price-Whelan, Maria Bergemann, Guy S. Stringfellow, Carla Frohlich, Henrique Reggiani, Erika M. Holmbeck, Jamie Tayar, Shivani P. Shah, Emily J. Griffith, Chervin F. P. Laporte, Andrew R. Casey, Keith Hawkins, Danny Horta, William Cerny, Pierre Thibodeaux, Sam A. Usman, Joao A. S. Amarante , et al. (17 additional authors not shown)

Abstract: Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931… ▽ More Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 11 pages + 22 page appendix, accepted to ApJL

The Heavy Metal Survey: The Evolution of Stellar Metallicities, Abundance Ratios, and Ages of Massive Quiescent Galaxies Since z~2

Authors: Aliza G. Beverage, Mariska Kriek, Katherine A. Suess, Charlie Conroy, Sedona H. Price, Guillermo Barro, Rachel Bezanson, Marijn Franx, Brian Lorenz, Yilun Ma, Lamiya Mowla, Imad Pasha, Pieter van Dokkum, Daniel Weisz

Abstract: We present the elemental abundances and ages of 19 massive quiescent galaxies at $z\sim1.4$ and $z\sim2.1$ from the Keck Heavy Metal Survey. The ultra-deep LRIS and MOSFIRE spectra were modeled using a full-spectrum stellar population fitting code with variable abundance patterns. The galaxies have iron abundances between [Fe/H] = -0.5 and -0.1 dex, with typical values of $-0.2$ [$-0.3$] at… ▽ More We present the elemental abundances and ages of 19 massive quiescent galaxies at $z\sim1.4$ and $z\sim2.1$ from the Keck Heavy Metal Survey. The ultra-deep LRIS and MOSFIRE spectra were modeled using a full-spectrum stellar population fitting code with variable abundance patterns. The galaxies have iron abundances between [Fe/H] = -0.5 and -0.1 dex, with typical values of $-0.2$ [$-0.3$] at $z\sim1.4$ [$z\sim2.1$]. We also find a tentative $\logσ_v$-[Fe/H] relation at $z\sim1.4$. The magnesium-to-iron ratios span [Mg/Fe] = 0.1--0.6 dex, with typical values of $0.3$ [$0.5$] dex at $z\sim1.4$ [$z\sim2.1$]. The ages imply formation redshifts of $z_{\rm form}=2-8$. Compared to quiescent galaxies at lower redshifts, we find [Fe/H] was $\sim0.2$ dex lower at $z=1.4-2.1$. We find no evolution in [Mg/Fe] out to $z\sim1.4$, though the $z\sim2.1$ galaxies are $0.2$ dex enhanced compared to $z=0-0.7$. A comparison of these results to a chemical evolution model indicates that galaxies at higher redshift form at progressively earlier epochs and over shorter star-formation timescales, with the $z\sim2.1$ galaxies forming the bulk of their stars over 150 Myr at $z_{\rm form}\sim4$. This evolution cannot be solely attributed to an increased number of quiescent galaxies at later times; several Heavy Metal galaxies have extreme chemical properties not found in massive galaxies at $z\sim0.0-0.7$. Thus, the chemical properties of individual galaxies must evolve over time. Minor mergers also cannot fully account for this evolution as they cannot increase [Fe/H], particularly in galaxy centers. Consequently, the build-up of massive quiescent galaxies since $z\sim2.1$ may require further mechanisms such as major mergers and/or central star formation. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Submitted to ApJ (21 pages, 11 figures)

The Heavy Metal Survey: Star Formation Constraints and Dynamical Masses of 21 Massive Quiescent Galaxies at $z=1.3-2.3$

Authors: Mariska Kriek, Aliza G. Beverage, Sedona H. Price, Katherine A. Suess, Guillermo Barro, Rachel S. Bezanson, Charlie Conroy, Sam E. Cutler, Marijn Franx, Jamie Lin, Brian Lorenz, Yilun Ma, Ivelina G. Momcheva, Lamiya A. Mowla, Imad Pasha, Pieter van Dokkum, Katherine E. Whitaker

Abstract: In this paper, we present the Heavy Metal Survey, which obtained ultradeep medium-resolution spectra of 21 massive quiescent galaxies at $1.3<z<2.3$ with Keck/LRIS and MOSFIRE. With integration times of up to 16\,hr per band per galaxy, we observe numerous Balmer and metal absorption lines in atmospheric windows. We successfully derive spectroscopic redshifts for all 21 galaxies and for 19 we also… ▽ More In this paper, we present the Heavy Metal Survey, which obtained ultradeep medium-resolution spectra of 21 massive quiescent galaxies at $1.3<z<2.3$ with Keck/LRIS and MOSFIRE. With integration times of up to 16\,hr per band per galaxy, we observe numerous Balmer and metal absorption lines in atmospheric windows. We successfully derive spectroscopic redshifts for all 21 galaxies and for 19 we also measure stellar velocity dispersions ($σ_v$), ages, and elemental abundances, as detailed in an accompanying paper. Except for one emission-line active galactic nucleus, all galaxies are confirmed as quiescent through their faint or absent H$α$ emission and evolved stellar spectra. For most galaxies exhibiting faint H$α$, elevated [NII]/H$α$ suggests a non-star-forming origin. We calculate dynamical masses ($M_{\rm dyn}$) by combining $σ_v$ with structural parameters obtained from HST/COSMOS(-DASH), and compare them with stellar masses ($M_*$) derived using spectrophotometric modeling, considering various assumptions. For a fixed initial mass function (IMF), we observe a strong correlation between $M_{\rm dyn}/M_*$ and $σ_v$. This correlation may suggest that a varying IMF, with high-$σ_v$ galaxies being more bottom heavy, was already in place at $z\sim2$. When implementing the $σ_v$-dependent IMF found in the cores of nearby early-type galaxies \textit{and} correcting for biases in our stellar mass and size measurements, we find a low scatter in $M_{\rm dyn}/M_*$ of 0.14 dex. However, these assumptions result in unphysical stellar masses, which exceed the dynamical masses by 34%. This tension suggests that distant quiescent galaxies do not simply grow inside-out into today's massive early-type galaxies and the evolution is more complicated. △ Less

Submitted 18 July, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: Published in ApJ

Journal ref: ApJ (2024), 966, 36

JWST Reveals Widespread AGN-Driven Neutral Gas Outflows in Massive z ~ 2 Galaxies

Authors: Rebecca L. Davies, Sirio Belli, Minjung Park, J. Trevor Mendel, Benjamin D. Johnson, Charlie Conroy, Chloë Benton, Letizia Bugiani, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: We use deep JWST/NIRSpec R~1000 slit spectra of 113 galaxies at 1.7 < z < 3.5, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na I D absorption (beyond the stellar contribution) in 46% of massive galaxies ($\log$ M$_*$/M$_\odot >$ 10), with similar incidence rates in star-forming and que… ▽ More We use deep JWST/NIRSpec R~1000 slit spectra of 113 galaxies at 1.7 < z < 3.5, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na I D absorption (beyond the stellar contribution) in 46% of massive galaxies ($\log$ M$_*$/M$_\odot >$ 10), with similar incidence rates in star-forming and quenching systems. Half of the absorption profiles are blueshifted by at least 100 km/s, providing unambiguous evidence for neutral gas outflows. Galaxies with strong Na I D absorption are distinguished by enhanced emission line ratios consistent with AGN ionization. We conservatively measure mass outflow rates of 3 - 100 $M_\odot$ yr$^{-1}$; comparable to or exceeding ionized gas outflow rates measured for galaxies at similar stellar mass and redshift. The outflows from the quenching systems (log(sSFR)[yr$^{-1}$] $\lesssim$ -10) have mass loading factors of 4 - 360, and the energy and momentum outflow rates exceed the expected injection rates from supernova explosions, suggesting that these galaxies could possibly be caught in a rapid blowout phase powered by the AGN. Our findings suggest that AGN-driven ejection of cold gas may be a dominant mechanism for fast quenching of star formation at z~2. △ Less

Submitted 30 January, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 17 pages, 8 figures, accepted for publication in MNRAS

The Three-Phase Evolution of the Milky Way

Authors: Vedant Chandra, Vadim A. Semenov, Hans-Walter Rix, Charlie Conroy, Ana Bonaca, Rohan P. Naidu, Rene Andrae, Jiadong Li, Lars Hernquist

Abstract: We illustrate the formation and evolution of the Milky Way over cosmic time, utilizing a sample of 10 million red giant stars with full chemodynamical information, including metallicities and $α$-abundances from low-resolution Gaia XP spectra. The evolution of angular momentum as a function of metallicity - a rough proxy for stellar age, particularly for high-[$α$/Fe] stars - displays three distin… ▽ More We illustrate the formation and evolution of the Milky Way over cosmic time, utilizing a sample of 10 million red giant stars with full chemodynamical information, including metallicities and $α$-abundances from low-resolution Gaia XP spectra. The evolution of angular momentum as a function of metallicity - a rough proxy for stellar age, particularly for high-[$α$/Fe] stars - displays three distinct phases: the disordered and chaotic protogalaxy, the kinematically-hot old disk, and the kinematically-cold young disk. The old high-$α$ disk starts at [Fe/H] $\approx -1.0$, 'spinning up' from the nascent protogalaxy, and then exhibits a smooth 'cooldown' toward more ordered and circular orbits at higher metallicities. The young low-$α$ disk is kinematically cold throughout its metallicity range, with its observed properties modulated by a strong radial gradient. We interpret these trends using Milky Way analogs from the TNG50 cosmological simulation, identifying one that closely matches the kinematic evolution of our Galaxy. This halo's protogalaxy spins up into a relatively thin and misaligned high-$α$ disk at early times, which is subsequently heated and torqued by a major gas-rich merger. The merger contributes a large amount of low-metallicity gas and angular momentum, from which the kinematically cold low-$α$ stellar disk is subsequently born. This simulated history parallels several observed features of the Milky Way, particularly the decisive 'GSE' merger that likely occurred at $z \approx 2$. Our results provide an all-sky perspective on the emerging picture of our Galaxy's three-phase formation, impelled by the three physical mechanisms of spinup, merger, and cooldown. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 21 pages, 12 figures, submitted to ApJ

Detection of Accretion Shelves Out to the Virial Radius of a Low-Mass Galaxy with JWST

Authors: Charlie Conroy, Benjamin D. Johnson, Pieter van Dokkum, Alis Deason, Sandro Tacchella, Sirio Belli, William P. Bowman, Rohan P. Naidu, Minjung Park, Roberto Abraham, Razieh Emami

Abstract: We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emis… ▽ More We report the serendipitous discovery of an extended stellar halo surrounding the low-mass galaxy Ark 227 ($M_\ast=5\times10^9 M_\odot$; d=35 Mpc) in deep JWST NIRCam imaging from the Blue Jay Survey. The F200W-F444W color provides robust star-galaxy separation, enabling the identification of stars at very low density. By combining resolved stars at large galactocentric distances with diffuse emission from NIRCam and Dragonfly imaging at smaller distances, we trace the surface brightness and color profiles of this galaxy over the entire extent of its predicted dark matter halo, from 0.1-100 kpc. Controlled N-body simulations have predicted that minor mergers create "accretion shelves" in the surface brightness profile at large radius. We observe such a feature in Ark 227 at 10-20 kpc, which, according to models, could be caused by a merger with total mass ratio 1:10. The metallicity declines over this radial range, further supporting the minor merger scenario. There is tentative evidence of a second shelf at $μ_V\approx 35$ mag arcsec$^{-2}$ extending from 50-100 kpc, along with a corresponding drop in metallicity. The stellar mass in this outermost envelope is $\approx10^7M_\odot$. These results suggest that Ark 227 experienced multiple mergers with a spectrum of lower-mass galaxies -- a scenario that is broadly consistent with the hierarchical growth of structure in a cold dark matter-dominated universe. Finally, we identify an ultra-faint dwarf associated with Ark 227 with $M_\ast\approx10^5 M_\odot$ and $μ_{V,e}=28.1$ mag arcsec$^{-2}$, demonstrating that JWST is capable of detecting very low-mass dwarfs to distances of at least ~30 Mpc. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 10 pages, 9 figures, submitted to ApJ

A massive compact quiescent galaxy at z=2 with a complete Einstein ring in JWST imaging

Authors: Pieter van Dokkum, Gabriel Brammer, Bingjie Wang, Joel Leja, Charlie Conroy

Abstract: One of the surprising results from HST was the discovery that many of the most massive galaxies at z~2 are very compact, having half-light radii of only 1-2 kpc. The interpretation is that massive galaxies formed inside-out, with their cores largely in place by z~2 and approximately half of their present-day mass added later through minor mergers. Here we present a compact, massive, quiescent gala… ▽ More One of the surprising results from HST was the discovery that many of the most massive galaxies at z~2 are very compact, having half-light radii of only 1-2 kpc. The interpretation is that massive galaxies formed inside-out, with their cores largely in place by z~2 and approximately half of their present-day mass added later through minor mergers. Here we present a compact, massive, quiescent galaxy at $z_{\rm phot}=1.94^{+0.13}_{-0.17}$ with a complete Einstein ring. The ring was found in the JWST COSMOS-Web survey and is produced by a background galaxy at $z_{\rm phot}=2.98^{+0.42}_{-0.47}$. Its 1.54" diameter provides a direct measurement of the mass of the "pristine" core of a massive galaxy, observed before mixing and dilution of its stellar population during the 10 Gyr of galaxy evolution between z=2 and z=0. We find a mass of $M_{\rm lens}=6.5^{+3.7}_{-1.5} \times 10^{11}$ Msun within a radius of 6.6 kpc. The stellar mass within the same radius is $M_{\rm stars}= 1.1^{+0.2}_{-0.3} \times 10^{11}$ Msun for a Chabrier initial mass function (IMF), and the fiducial dark matter mass is $M_{\rm dm} = 2.6^{+1.6}_{-0.7} \times 10^{11}$ Msun. Additional mass is needed to explain the lensing results, either in the form of a higher-than-expected dark matter density or a bottom-heavy IMF. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Nature Astronomy, in press. As it happens this is the highest redshift lens currently known. Found by eyeballing the wonderful COSMOS-Web data

A Tilted Dark Halo Origin of the Galactic Disk Warp and Flare

Authors: Jiwon Jesse Han, Charlie Conroy, Lars Hernquist

Abstract: The outer disk of the Milky Way Galaxy is warped and flared. Several mechanisms have been proposed to explain these phenomena, but none have quantitatively reproduced both features. Recent work has demonstrated that the Galactic stellar halo is tilted with respect to the disk plane, suggesting that at least some component of the dark matter halo may also be tilted. Here we show that a dark halo ti… ▽ More The outer disk of the Milky Way Galaxy is warped and flared. Several mechanisms have been proposed to explain these phenomena, but none have quantitatively reproduced both features. Recent work has demonstrated that the Galactic stellar halo is tilted with respect to the disk plane, suggesting that at least some component of the dark matter halo may also be tilted. Here we show that a dark halo tilted in the same direction as the stellar halo can induce a warp and flare in the Galactic disk at the same amplitude and orientation as the data. In our model the warp is visible in both the gas and stars of all ages, which is consistent with the breadth of observational tracers of the warp. These results, in combination with data in the stellar halo, provide compelling evidence that our Galaxy is embedded in a tilted dark matter halo. This misalignment of the dark halo and the disk holds clue to the formation history of the Galaxy, and represents the next step in the dynamical modeling of the Galactic potential. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: Published in Nature Astronomy 14 September 2023

Tilted Dark Halos are Common, Long-Lived, and can Warp Galactic Disks

Authors: Jiwon Jesse Han, Vadim Semenov, Charlie Conroy, Lars Hernquist

Abstract: In the $Λ$-CDM paradigm, the dark halo governs the gravitational potential within which a galaxy can form and evolve. In this Letter we show that the present-day inner ($r<50\text{ kpc}$) dark halo can be significantly misaligned with the stellar disk. To this end, we use the TNG50 run from the cosmological magneto-hydrodynamic IllustrisTNG simulation suite. Such "tilted" dark halos can arise from… ▽ More In the $Λ$-CDM paradigm, the dark halo governs the gravitational potential within which a galaxy can form and evolve. In this Letter we show that the present-day inner ($r<50\text{ kpc}$) dark halo can be significantly misaligned with the stellar disk. To this end, we use the TNG50 run from the cosmological magneto-hydrodynamic IllustrisTNG simulation suite. Such "tilted" dark halos can arise from a variety of processes including major mergers, massive fly-bys, or interactions with satellite companions. Furthermore, we show that tilted dark halos: (1) are well traced by tilted stellar halos, (2) can maintain their tilt for $>$ 5 Gyr in isolated evolution, and (3) can generate warps in the outer disks that are stable over many Gyr. A tilted dark halo holds clues to important events in the formation history of a galaxy, and could help explain the abundance of warped disks in galaxy observations, including the Milky Way. △ Less

Submitted 14 September, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: Submitted to ApJL (previous version had wrong abstract)

Extending the Chemical Reach of the H3 Survey: Detailed Abundances of the Dwarf-galaxy Stellar Stream Wukong/LMS-1

Authors: Guilherme Limberg, Alexander P. Ji, Rohan P. Naidu, Anirudh Chiti, Silvia Rossi, Sam A. Usman, Yuan-Sen Ting, Dennis Zaritsky, Ana Bonaca, Lais Borbolato, Joshua S. Speagle, Vedant Chandra, Charlie Conroy

Abstract: We present the first detailed chemical-abundance analysis of stars from the dwarf-galaxy stellar stream Wukong/LMS-1 covering a wide metallicity range ($-3.5 < \rm[Fe/H] \lesssim -1.3$). We find abundance patterns that are effectively indistinguishable from the bulk of Indus and Jhelum, a pair of smaller stellar streams proposed to be dynamically associated with Wukong/LMS-1. We confirmed a carbon… ▽ More We present the first detailed chemical-abundance analysis of stars from the dwarf-galaxy stellar stream Wukong/LMS-1 covering a wide metallicity range ($-3.5 < \rm[Fe/H] \lesssim -1.3$). We find abundance patterns that are effectively indistinguishable from the bulk of Indus and Jhelum, a pair of smaller stellar streams proposed to be dynamically associated with Wukong/LMS-1. We confirmed a carbon-enhanced metal-poor star ($\rm[C/Fe] > +0.7$ and $\rm[Fe/H] \sim -2.9$) in Wukong/LMS-1 with strong enhancements in Sr, Y, and Zr, which is peculiar given its solar-level [Ba/Fe]. Wukong/LMS-1 stars have high abundances of $α$ elements up to $\rm[Fe/H] \gtrsim -2$, which is expected for relatively massive dwarfs. Towards the high-metallicity end, Wukong/LMS-1 becomes $α$-poor, revealing that it probably experienced fairly standard chemical evolution. We identified a pair of N- and Na-rich stars in Wukong/LMS-1, reminiscent of multiple populations in globular clusters. This indicates that this dwarf galaxy contained at least one globular cluster that was completely disrupted in addition to two intact ones previously known to be associated with Wukong/LMS-1, which is possibly connected to similar evidence found in Indus. From these $\geq$3 globular clusters, we estimate the total mass of Wukong/LMS-1 to be ${\approx}10^{10} M_\odot$, representing ${\sim}1$% of the present-day Milky Way. Finally, the [Eu/Mg] ratio in Wukong/LMS-1 continuously increases with metallicity, making this the first example of a dwarf galaxy where the production of $r$-process elements is clearly dominated by delayed sources, presumably neutron-star mergers. △ Less

Submitted 5 April, 2024; v1 submitted 25 August, 2023; originally announced August 2023.

Comments: Accepted to MNRAS. New version fixes abundance uncertainties, which significantly affects elements like Al and N. Use new version of abundance files for correct error bars

Star Formation Shut Down by Multiphase Gas Outflow in a Galaxy at a Redshift of 2.45

Authors: Sirio Belli, Minjung Park, Rebecca L. Davies, J. Trevor Mendel, Benjamin D. Johnson, Charlie Conroy, Chloë Benton, Letizia Bugiani, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: Large-scale outflows driven by supermassive black holes are thought to play a fundamental role in suppressing star formation in massive galaxies. However, direct observational evidence for this hypothesis is still lacking, particularly in the young universe where star formation quenching is remarkably rapid, thus requiring effective removal of gas as opposed to slow gas heating. While outflows of… ▽ More Large-scale outflows driven by supermassive black holes are thought to play a fundamental role in suppressing star formation in massive galaxies. However, direct observational evidence for this hypothesis is still lacking, particularly in the young universe where star formation quenching is remarkably rapid, thus requiring effective removal of gas as opposed to slow gas heating. While outflows of ionized gas are commonly detected in massive distant galaxies, the amount of ejected mass is too small to be able to suppress star formation. Gas ejection is expected to be more efficient in the neutral and molecular phases, but at high redshift these have only been observed in starbursts and quasars. Here we report JWST spectroscopy of a massive galaxy experiencing rapid quenching at redshift z=2.445. We detect a weak outflow of ionized gas and a powerful outflow of neutral gas, with a mass outflow rate that is sufficient to quench the star formation. Neither X-ray or radio activity are detected; however, the presence of a supermassive black hole is suggested by the properties of the ionized gas emission lines. We thus conclude that supermassive black holes are able to rapidly suppress star formation in massive galaxies by efficiently ejecting neutral gas. △ Less

Submitted 10 May, 2024; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 21 pages, 9 figures. Published in Nature

Discovery of the Magellanic Stellar Stream Out to 100 Kiloparsecs

Authors: Vedant Chandra, Rohan P. Naidu, Charlie Conroy, Ana Bonaca, Dennis Zaritsky, Phillip A. Cargile, Nelson Caldwell, Benjamin D. Johnson, Jiwon Jesse Han, Yuan-Sen Ting

Abstract: The Magellanic Stream (MS) - an enormous ribbon of gas spanning $140^\circ$ of the southern sky trailing the Magellanic Clouds - has been exquisitely mapped in the five decades since its discovery. However, despite concerted efforts, no stellar counterpart to the MS has been conclusively identified. This stellar stream would reveal the distance and 6D kinematics of the MS, constraining its formati… ▽ More The Magellanic Stream (MS) - an enormous ribbon of gas spanning $140^\circ$ of the southern sky trailing the Magellanic Clouds - has been exquisitely mapped in the five decades since its discovery. However, despite concerted efforts, no stellar counterpart to the MS has been conclusively identified. This stellar stream would reveal the distance and 6D kinematics of the MS, constraining its formation and the past orbital history of the Clouds. We have been conducting a spectroscopic survey of the most distant and luminous red giant stars in the Galactic outskirts. From this dataset, we have discovered a prominent population of 13 stars matching the extreme angular momentum of the Clouds, spanning up to $100^\circ$ along the MS at distances of $60-120$ kpc. Furthermore, these kinemetically-selected stars lie along a [$α$/Fe]-deficient track in chemical space from $-2.5 < \mathrm{[Fe/H]} < -0.5$, consistent with their formation in the Clouds themselves. We identify these stars as high-confidence members of the Magellanic Stellar Stream. Half of these stars are metal-rich and closely follow the gaseous MS, whereas the other half are more scattered and metal-poor. We argue that the metal-rich stream is the recently-formed tidal counterpart to the MS, and speculate that the metal-poor population was thrown out of the SMC outskirts during an earlier interaction between the Clouds. The Magellanic Stellar Stream provides a strong set of constraints - distances, 6D kinematics, and birth locations - that will guide future simulations towards unveiling the detailed history of the Clouds. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 21 pages, 12 figures. Submitted to ApJ

Formation of Galactic Disks II: the Physical Drivers of Disk Spin-up

Authors: Vadim A. Semenov, Charlie Conroy, Vedant Chandra, Lars Hernquist, Dylan Nelson

Abstract: Using a representative sample of Milky Way (MW)-like galaxies from the TNG50 cosmological-volume simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum… ▽ More Using a representative sample of Milky Way (MW)-like galaxies from the TNG50 cosmological-volume simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum, suggesting an exchange of angular momentum and/or outflow recycling, leading to continuous feeding of pre-aligned material from the co-rotating circumgalactic medium. We show that disk formation in TNG50 is correlated with stellar bulge formation, in qualitative agreement with a recent theoretical model of disk formation facilitated by steep gravitational potentials. Disk formation is also correlated with the formation of a hot circumgalactic halo with a significant fraction of the inflow occurring at sub- and transonic velocities. In the context of recent theoretical works connecting disk settling and hot halo formation, our results imply that the subsonic part of the inflow may settle into a disk while the remaining supersonic inflow will perturb this disk via the chaotic cold accretion. We find that disks tend to form when the host halos become more massive than $\sim (1-2) \times 10^{11} M_\odot$, consistent with previous theoretical findings and observational estimates of the pre-disk protogalaxy remnant in the MW. Our results do not prove that either co-rotating outflow recycling, gravitational potential steepening, or hot halo formation cause disk formation but they show that all these processes occur concurrently and may play an important role in disk growth. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 22 pages, 15 figures; submitted to ApJ; comments are welcome

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Formation of Galactic Disks I: Why Did the Milky Way's Disk Form Unusually Early?

Authors: Vadim A. Semenov, Charlie Conroy, Vedant Chandra, Lars Hernquist, Dylan Nelson

Abstract: Recent results from spectroscopic and astrometric surveys of nearby stars suggest that the stellar disk of our Milky Way (MW) was formed quite early, within the first few billion years of its evolution. Chemokinematic signatures of disk formation in cosmological zoom-in simulations appear to be in tension with these data, implying that MW-like disk formation is delayed in simulations. We investiga… ▽ More Recent results from spectroscopic and astrometric surveys of nearby stars suggest that the stellar disk of our Milky Way (MW) was formed quite early, within the first few billion years of its evolution. Chemokinematic signatures of disk formation in cosmological zoom-in simulations appear to be in tension with these data, implying that MW-like disk formation is delayed in simulations. We investigate the formation of galactic disks using a representative sample of MW-like galaxies from the cosmological-volume simulation TNG50. We find that on average MW-mass disks indeed form later than the local data suggest. However, their formation time and metallicity exhibit a substantial scatter, such that $\sim$10% of MW-mass galaxies form disks early, similar to the MW. Thus, although the MW is unusual, it is consistent with the overall population of MW-mass disk galaxies. The direct MW analogs assemble most of their mass early, $\gtrsim 10$ Gyr ago, and are not affected by destructive mergers after that. In addition, these galaxies form their disks during the early enrichment stage when the interstellar medium metallicity increases rapidly, with only $\sim$25% of early-forming disks being as metal-poor as the MW was at the onset of disk formation, [Fe/H] $\approx -1.0$. In contrast, most MW-mass galaxies either form disks from already enriched material or experience late destructive mergers that reset the signatures of galactic disk formation to later times and higher metallicities. Finally, we also show that earlier disk formation leads to more dominant rotationally supported stellar disks at redshift zero. △ Less

Submitted 22 January, 2024; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 19 pages, 13 figures + appendix; accepted for publication in ApJ

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

From Carbon to Cobalt: Chemical compositions and ages of $z\sim0.7$ quiescent galaxies

Authors: Aliza G. Beverage, Mariska Kriek, Charlie Conroy, Nathan R. Sandford, Rachel Bezanson, Marijn Franx, Arjen van der Wel, Daniel R. Weisz

Abstract: We present elemental abundance patterns (C, N, Mg, Si, Ca, Ti, V, Cr, Fe, Co, and Ni) for a population of 135 massive quiescent galaxies at $z\sim0.7$ with ultra-deep rest-frame optical spectroscopy drawn from the LEGA-C survey. We derive average ages and elemental abundances in four bins of stellar velocity dispersion ($σ_v$) ranging from 150$~$km$\,$s$^{-1}$ to 250$~$km$\,$s$^{-1}$ using a full-… ▽ More We present elemental abundance patterns (C, N, Mg, Si, Ca, Ti, V, Cr, Fe, Co, and Ni) for a population of 135 massive quiescent galaxies at $z\sim0.7$ with ultra-deep rest-frame optical spectroscopy drawn from the LEGA-C survey. We derive average ages and elemental abundances in four bins of stellar velocity dispersion ($σ_v$) ranging from 150$~$km$\,$s$^{-1}$ to 250$~$km$\,$s$^{-1}$ using a full-spectrum hierarchical Bayesian model. The resulting elemental abundance measurements are precise to 0.05$\,$dex. The majority of elements, as well as the total metallicity and stellar age, show a positive correlation with $σ_v$. Thus, the highest dispersion galaxies formed the earliest and are the most metal-rich. We find only mild or non-significant trends between [X/Fe] and $σ_v$, suggesting that the average star-formation timescale does not strongly depend on velocity dispersion. To first order, the abundance patterns of the $z\sim0.7$ quiescent galaxies are strikingly similar to those at $z\sim0$. However, at the lowest velocity dispersions the $z\sim0.7$ galaxies have slightly enhanced N, Mg, Ti, and Ni abundance ratios and earlier formation redshifts than their $z\sim0$ counterparts. Thus, while the higher-mass quiescent galaxy population shows little evolution, the low-mass quiescent galaxies population has grown significantly over the past six billion years. Finally, the abundance patterns of both $z\sim0$ and $z\sim0.7$ quiescent galaxies differ considerably from theoretical prediction based on a chemical evolution model, indicating that our understanding of the enrichment histories of these galaxies is still very limited. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 24 pages, 12 figures, accepted for publication in ApJ

A red giant orbiting a black hole

Authors: Kareem El-Badry, Hans-Walter Rix, Yvette Cendes, Antonio C. Rodriguez, Charlie Conroy, Eliot Quataert, Keith Hawkins, Eleonora Zari, Melissa Hobson, Katelyn Breivik, Arne Rau, Edo Berger, Sahar Shahaf, Rhys Seeburger, Kevin B. Burdge, David W. Latham, Lars A. Buchhave, Allyson Bieryla, Dolev Bashi, Tsevi Mazeh, Simchon Faigler

Abstract: We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. The system, which we call Gaia BH2, contains a $\sim 1M_{\odot}$ red giant and a dark companion with mass $M_2 = 8.9\pm 0.3\,M_{\odot}$ that is very likely a BH. The orbital period, $P_{\rm orb} = 1277$ days, is much longer than that of any previously studied BH binary. Our radial velocity (RV)… ▽ More We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. The system, which we call Gaia BH2, contains a $\sim 1M_{\odot}$ red giant and a dark companion with mass $M_2 = 8.9\pm 0.3\,M_{\odot}$ that is very likely a BH. The orbital period, $P_{\rm orb} = 1277$ days, is much longer than that of any previously studied BH binary. Our radial velocity (RV) follow-up over a 7-month period spans more than 90% of the orbit's dynamic range in RV and is in excellent agreement with predictions of the Gaia solution. UV imaging and high-resolution optical spectra rule out all plausible luminous companions that could explain the orbit. The star is a bright ($G=12.3$), slightly metal-poor ($\rm [Fe/H]=-0.22$) low-luminosity giant ($T_{\rm eff}=4600\,\rm K$; $R = 7.8\,R_{\odot}$; $\log\left[g/\left({\rm cm\,s^{-2}}\right)\right] = 2.6$). The binary's orbit is moderately eccentric ($e=0.52$). The giant is strongly enhanced in $α-$elements, with $\rm [α/Fe] = +0.26$, but the system's Galactocentric orbit is typical of the thin disk. We obtained X-ray and radio nondetections of the source near periastron, which support BH accretion models in which the net accretion rate at the horizon is much lower than the Bondi-Hoyle-Lyttleton rate. At a distance of 1.16 kpc, Gaia BH2 is the second-nearest known BH, after Gaia BH1. Its orbit -- like that of Gaia BH1 -- seems too wide to have formed through common envelope evolution. Gaia BH1 and BH2 have orbital periods at opposite edges of the Gaia DR3 sensitivity curve, perhaps hinting at a bimodal intrinsic period distribution for wide BH binaries. Dormant BH binaries like Gaia BH1 and Gaia BH2 likely significantly outnumber their close, X-ray bright cousins, but their formation pathways remain uncertain. △ Less

Submitted 19 March, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: Accepted to MNRAS

A candidate runaway supermassive black hole identified by shocks and star formation in its wake

Authors: Pieter van Dokkum, Imad Pasha, Maria Luisa Buzzo, Stephanie LaMassa, Zili Shen, Michael A. Keim, Roberto Abraham, Charlie Conroy, Shany Danieli, Kaustav Mitra, Daisuke Nagai, Priyamvada Natarajan, Aaron J. Romanowsky, Grant Tremblay, C. Megan Urry, Frank C. van den Bosch

Abstract: The interaction of a runaway supermassive black hole (SMBH) with the circumgalactic medium (CGM) can lead to the formation of a wake of shocked gas and young stars behind it. Here we report the serendipitous discovery of an extremely narrow linear feature in HST/ACS images that may be an example of such a wake. The feature extends 62 kpc from the nucleus of a compact star-forming galaxy at z=0.964… ▽ More The interaction of a runaway supermassive black hole (SMBH) with the circumgalactic medium (CGM) can lead to the formation of a wake of shocked gas and young stars behind it. Here we report the serendipitous discovery of an extremely narrow linear feature in HST/ACS images that may be an example of such a wake. The feature extends 62 kpc from the nucleus of a compact star-forming galaxy at z=0.964. Keck LRIS spectra show that the [OIII]/H$β$ ratio varies from ~1 to ~10 along the feature, indicating a mixture of star formation and fast shocks. The feature terminates in a bright [OIII] knot with a luminosity of 1.9x10$^{41}$ ergs/s. The stellar continuum colors vary along the feature, and are well-fit by a simple model that has a monotonically increasing age with distance from the tip. The line ratios, colors, and the overall morphology are consistent with an ejected SMBH moving through the CGM at high speed while triggering star formation. The best-fit time since ejection is ~39 Myr and the implied velocity is v~1600 km/s. The feature is not perfectly straight in the HST images, and we show that the amplitude of the observed spatial variations is consistent with the runaway SMBH interpretation. Opposite the primary wake is a fainter and shorter feature, marginally detected in [OIII] and the rest-frame far-ultraviolet. This feature may be shocked gas behind a binary SMBH that was ejected at the same time as the SMBH that produced the primary wake. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJ Letters. The key data are in Figure 1: a really odd thin streak in HST images, with a complex emission line spectrum. Figure 7 is an illustration of our proposed interpretation

Journal ref: Astrophysical Journal Letters, 946, L50 (2023)

The Eighteenth Data Release of the Sloan Digital Sky Surveys: Targeting and First Spectra from SDSS-V

Authors: Andrés Almeida, Scott F. Anderson, Maria Argudo-Fernández, Carles Badenes, Kat Barger, Jorge K. Barrera-Ballesteros, Chad F. Bender, Erika Benitez, Felipe Besser, Dmitry Bizyaev, Michael R. Blanton, John Bochanski, Jo Bovy, William Nielsen Brandt, Joel R. Brownstein, Johannes Buchner, Esra Bulbul, Joseph N. Burchett, Mariana Cano Díaz, Joleen K. Carlberg, Andrew R. Casey, Vedant Chandra, Brian Cherinka, Cristina Chiappini, Abigail A. Coker , et al. (129 additional authors not shown)

Abstract: The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM),… ▽ More The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM), including input catalogs and selection functions for their numerous scientific objectives. We describe the production of the targeting databases and their calibration- and scientifically-focused components. DR18 also includes ~25,000 new SDSS spectra and supplemental information for X-ray sources identified by eROSITA in its eFEDS field. We present updates to some of the SDSS software pipelines and preview changes anticipated for DR19. We also describe three value-added catalogs (VACs) based on SDSS-IV data that have been published since DR17, and one VAC based on the SDSS-V data in the eFEDS field. △ Less

Submitted 6 July, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: Accepted to ApJS

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

Authors: Daniel R. Weisz, Kristen B. W. McQuinn, Alessandro Savino, Nitya Kallivayalil, Jay Anderson, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Andrew E. Dolphin, Karin M. Sandstrom, Andrew A. Cole, Benjamin F. Williams, Evan D. Skillman, Roger E. Cohen, Max J. B. Newman, Rachael Beaton, Alessandro Bressan, Alberto Bolatto, Michael Boylan-Kolchin, Alyson M. Brooks, James S. Bullock, Charlie Conroy, M. C. Cooper, Julianne J. Dalcanton, Aaron L. Dotter , et al. (17 additional authors not shown)

Abstract: We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We descr… ▽ More We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 23 pages, 7 Figures, 3 Tables. Submitted to AAS Journals. Comments welcome

Distant Echoes of the Milky Way's Last Major Merger

Authors: Vedant Chandra, Rohan P. Naidu, Charlie Conroy, Alexander P. Ji, Hans-Walter Rix, Ana Bonaca, Phillip Cargile, Jiwon Jesse Han, Benjamin D. Johnson, Yuan-Sen Ting, Turner Woody, Dennis Zaritsky

Abstract: The majority of the Milky Way's stellar halo consists of debris from our Galaxy's last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from GSE have been kinematically and chemically studied in the inner $30$ kpc of our Galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metal… ▽ More The majority of the Milky Way's stellar halo consists of debris from our Galaxy's last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from GSE have been kinematically and chemically studied in the inner $30$ kpc of our Galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to $100$ kpc. We obtain follow-up spectra of stars in two strong overdensities - including the previously identified Outer Virgo Overdensity - and find them to be relatively metal-rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming $60-90$ kpc counterparts to the $15-20$ kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from $50-100$ kpc, in the same plane as the Sagittarius stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our Galaxy's stellar halo. △ Less

Submitted 30 June, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 22 pages, 10 figures. Accepted to ApJ

Journal ref: The Astrophysical Journal, 951, 26 (2023)

What is Missing from the Local Stellar Halo?

Authors: Katherine Sharpe, Rohan P. Naidu, Charlie Conroy

Abstract: The Milky Way's stellar halo, which extends to $>100$ kpc, encodes the evolutionary history of our Galaxy. However, most studies of the halo to date have been limited to within a few kpc of the Sun. Here, we characterize differences between this local halo and the stellar halo in its entirety. We construct a composite stellar halo model by combining observationally motivated N-body simulations of… ▽ More The Milky Way's stellar halo, which extends to $>100$ kpc, encodes the evolutionary history of our Galaxy. However, most studies of the halo to date have been limited to within a few kpc of the Sun. Here, we characterize differences between this local halo and the stellar halo in its entirety. We construct a composite stellar halo model by combining observationally motivated N-body simulations of the Milky Way's nine most massive disrupted dwarf galaxies that account for almost all of the mass in the halo. We find that (1) the representation by mass of different dwarf galaxies in the local halo compared to the whole halo can be significantly overestimated (e.g., the Helmi Streams) or underestimated (e.g., Cetus) and (2) properties of the overall halo (e.g., net rotation) inferred via orbit integration of local halo stars are significantly biased, because e.g., highly retrograde debris from Gaia-Sausage-Enceladus is missing from the local halo. Therefore, extrapolations from the local to the global halo should be treated with caution. From analysis of a sample of 11 MW-like simulated halos, we identify a population of recently accreted ($\lesssim5$ Gyrs) and disrupted galaxies on high angular momenta orbits that are entirely missing from local samples, and awaiting discovery in the outer halo. Our results motivate the need for surveys of halo stars extending to the Galaxy's virial radius. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 12 pages, 6 figures. Submitted to ApJ. Comments warmly welcomed

Novae in M51: a New, Much Higher Rate from Multi-epoch HST Data

Authors: Shifra Mandel, Michael M. Shara, David Zurek, Charlie Conroy, Pieter van Dokkum

Abstract: Accurate determination of the rates of nova eruptions in different kinds of galaxies give us strong constraints on those galaxies' underlying white dwarf and binary populations, and those stars' spatial distributions. Until 2016, limitations inherent in ground-based surveys of external galaxies - and dust extinction in the Milky Way - significantly hampered the determination of those rates and how… ▽ More Accurate determination of the rates of nova eruptions in different kinds of galaxies give us strong constraints on those galaxies' underlying white dwarf and binary populations, and those stars' spatial distributions. Until 2016, limitations inherent in ground-based surveys of external galaxies - and dust extinction in the Milky Way - significantly hampered the determination of those rates and how much they differ between different types of galaxies. Infrared Galactic surveys and dense cadence Hubble Space Telescope (HST)-based surveys are overcoming these limitations, leading to sharply increased nova-in-galaxy rates relative to those previously claimed. Here we present 14 nova candidates that were serendipitously observed during a year-long HST survey of the massive spiral galaxy M51 (the "Whirlpool Galaxy"). We use simulations based on observed nova light curves to model the incompleteness of the HST survey in unprecedented detail, determining a nova detection efficiency $ε= 20.3$ percent. The survey's M51 area coverage, combined with $ε$, indicates a conservative M51 nova rate of $172^{+46}_{-37}$ novae yr$^{-1}$, corresponding to a luminosity-specific nova rate (LSNR) of $\sim10.4^{+2.8}_{-2.2}$ novae yr$^{-1}$/$10^{10} L_{\odot,K}$. Both these rates are approximately an order of magnitude higher than those estimated by ground-based studies, contradicting claims of universal low nova rates in all types of galaxies determined by low cadence, ground-based surveys. They demonstrate that, contrary to theoretical models, the HST-determined LSNR in a giant elliptical galaxy (M87) and a giant spiral galaxy (M51) likely do not differ by an order of magnitude or more, and may in fact be quite similar. △ Less

Submitted 17 October, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Rapid Quenching of Galaxies at Cosmic Noon

Authors: Minjung Park, Sirio Belli, Charlie Conroy, Sandro Tacchella, Joel Leja, Sam E. Cutler, Benjamin D. Johnson, Erica J. Nelson, Razieh Emami

Abstract: The existence of massive quiescent galaxies at high redshift seems to require rapid quenching, but it is unclear whether all quiescent galaxies have gone through this phase and what physical mechanisms are involved. To study rapid quenching, we use rest-frame colors to select 12 young quiescent galaxies at $z \sim 1.5$. From spectral energy distribution fitting, we find that they all experienced i… ▽ More The existence of massive quiescent galaxies at high redshift seems to require rapid quenching, but it is unclear whether all quiescent galaxies have gone through this phase and what physical mechanisms are involved. To study rapid quenching, we use rest-frame colors to select 12 young quiescent galaxies at $z \sim 1.5$. From spectral energy distribution fitting, we find that they all experienced intense starbursts prior to rapid quenching. We confirm this with deep Magellan/FIRE spectroscopic observations for a subset of seven galaxies. Broad emission lines are detected for two galaxies and are most likely caused by AGN activity. The other five galaxies do not show any emission features, suggesting that gas has already been removed or depleted. Most of the rapidly quenched galaxies are more compact than normal quiescent galaxies, providing evidence for a central starburst in the recent past. We estimate an average transition time of $300\,\rm Myr$ for the rapid quenching phase. Approximately $4\%$ of quiescent galaxies at $z=1.5$ have gone through rapid quenching; this fraction increases to $23\%$ at $z=2.2$. We identify analogs in the TNG100 simulation and find that rapid quenching for these galaxies is driven by AGN, and for half of the cases, gas-rich major mergers seem to trigger the starburst. We conclude that these massive quiescent galaxies are not just rapidly quenched but also rapidly formed through a major starburst. We speculate that mergers drive gas inflow towards the central regions and grow supermassive black holes, leading to rapid quenching by AGN feedback. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: Submitted to ApJ. Comments are welcome

Dwarf galaxy archaeology from chemical abundances and star formation histories

Authors: James W. Johnson, Charlie Conroy, Benjamin D. Johnson, Annika H. G. Peter, Phillip A. Cargile, Ana Bonaca, Rohan P. Naidu, Turner Woody, Yuan-Sen Ting, Jiwon Jesse Han, Joshua S. Speagle

Abstract: We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding timescales of $τ_\text{in} = 1.01 \pm 0.13$ Gyr for GSE and… ▽ More We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding timescales of $τ_\text{in} = 1.01 \pm 0.13$ Gyr for GSE and $τ_\text{in} = 3.08^{+3.19}_{-1.16}$ Gyr for Wukong/LMS-1. GSE formed stars for $τ_\text{tot} = 5.40^{+0.32}_{-0.31}$ Gyr, sustaining star formation for $\sim$$1.5 - 2$ Gyr after its first infall into the Milky Way $\sim$10 Gyr ago. Our fit suggests that star formation lasted for $τ_\text{tot} = 3.36^{+0.55}_{-0.47}$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass $M_\star$ predicted by simulations and semi-analytic models of galaxy formation. Our fitting method is based only on poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes ($20 \leq N \leq 2000$) and measurement uncertainties ($0.01 \leq σ_\text{[$α$/Fe]}, σ_\text{[Fe/H]} \leq 0.5$; $0.02 \leq σ_{\log_{10}(\text{age})} \leq 1$). Our inferred values of the outflow mass-loading factor reasonably match $η\propto M_\star^{-1/3}$ as predicted by galactic wind models. Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space. △ Less

Submitted 4 October, 2022; originally announced October 2022.

Comments: 25 pages; 13 figures; submitted to MNRAS; comments welcome

No peaks without valleys: The stable mass transfer channel for gravitational-wave sources in light of the neutron star-black hole mass gap

Authors: L. A. C. van Son, S. E. de Mink, M. Renzo, S. Justham, E. Zapartas, K. Breivik, T. Callister, W. M. Farr, C. Conroy

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… ▽ More 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 $\sim9 \mathrm{M_{\odot}}$ 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 BH mass gap. The peak and cutoff mass are a consequence of 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 the initial mass. We provide an analytical expression for the cutoff in the primary component mass and show that this adequately matches our numerical results. Our results imply that selection effects resulting from the formation channel alone can provide an explanation for the purported neutron star--BH mass gap in GW detections. This provides an alternative to the commonly adopted view that the gap emerges during BH formation. △ Less

Submitted 10 November, 2022; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ associated code is available at https://github.com/LiekeVanSon/LowMBH_and_StableChannel

A Sun-like star orbiting a black hole

Authors: Kareem El-Badry, Hans-Walter Rix, Eliot Quataert, Andrew W. Howard, Howard Isaacson, Jim Fuller, Keith Hawkins, Katelyn Breivik, Kaze W. K. Wong, Antonio C. Rodriguez, Charlie Conroy, Sahar Shahaf, Tsevi Mazeh, Frédéric Arenou, Kevin B. Burdge, Dolev Bashi, Simchon Faigler, Daniel R. Weisz, Rhys Seeburger, Silvia Almada Monter, Jennifer Wojno

Abstract: We report discovery of a bright, nearby ($G = 13.8;\,\,d = 480\,\rm pc$) Sun-like star orbiting a dark object. We identified the system as a black hole candidate via its astrometric orbital solution from the Gaia mission. Radial velocities validated and refined the Gaia solution, and spectroscopy ruled out significant light contributions from another star. Joint modeling of radial velocities and a… ▽ More We report discovery of a bright, nearby ($G = 13.8;\,\,d = 480\,\rm pc$) Sun-like star orbiting a dark object. We identified the system as a black hole candidate via its astrometric orbital solution from the Gaia mission. Radial velocities validated and refined the Gaia solution, and spectroscopy ruled out significant light contributions from another star. Joint modeling of radial velocities and astrometry constrains the companion mass to $M_2 = 9.62\pm 0.18\,M_{\odot}$. The spectroscopic orbit alone sets a minimum companion mass of $M_2>5\,M_{\odot}$; if the companion were a $5\,M_{\odot}$ star, it would be $500$ times more luminous than the entire system. These constraints are insensitive to the mass of the luminous star, which appears as a slowly-rotating G dwarf ($T_{\rm eff}=5850\,\rm K$, $\log g = 4.5$, $M=0.93\,M_{\odot}$), with near-solar metallicity ($\rm [Fe/H] = -0.2$) and an unremarkable abundance pattern. We find no plausible astrophysical scenario that can explain the orbit and does not involve a black hole. The orbital period, $P_{\rm orb}=185.6$ days, is longer than that of any known stellar-mass black hole binary. The system's modest eccentricity ($e=0.45$), high metallicity, and thin-disk Galactic orbit suggest that it was born in the Milky Way disk with at most a weak natal kick. How the system formed is uncertain. Common envelope evolution can only produce the system's wide orbit under extreme and likely unphysical assumptions. Formation models involving triples or dynamical assembly in an open cluster may be more promising. This is the nearest known black hole by a factor of 3, and its discovery suggests the existence of a sizable population of dormant black holes in binaries. Future Gaia releases will likely facilitate the discovery of dozens more. △ Less

Submitted 28 February, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 30 pages, 15 figures, 1-2 black holes. Accepted to MNRAS

The locations of features in the mass distribution of merging binary black holes are robust against uncertainties in the metallicity-dependent cosmic star formation history

Authors: L. A. C. van Son, S. E. de Mink, M. Chruslinska, C. Conroy, R. Pakmor, L. Hernquist

Abstract: New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation h… ▽ More New observational facilities are probing astrophysical transients such as stellar explosions and gravitational wave (GW) sources at ever increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history ($\mathcal{S}(Z,z)$) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for $\mathcal{S}(Z,z)$. Variations of this function can be easily interpreted, because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behaviour well. As an example application, we investigate the effect of systematic variations in the $\mathcal{S}(Z,z)$ parameters on the predicted mass distribution of locally merging binary black holes (BBH). Our main findings are: I) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and II) the low mass end is least affected by these variations. This is promising as it increases our chances to constrain the physics that governs the formation of these objects. △ Less

Submitted 27 February, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: Accepted for publication in ApJ, 6 Figures, code available at https://github.com/LiekeVanSon/SFRD_fit, made using showyourwork

Stellar Bars in Isolated Gas-Rich Spiral Galaxies Do Not Slow Down

Authors: Angus Beane, Lars Hernquist, Elena D'Onghia, Federico Marinacci, Charlie Conroy, Jia Qi, Laura V. Sales, Paul Torrey, Mark Vogelsberger

Abstract: Elongated bar-like features are ubiquitous in galaxies, occurring at the centers of approximately two-thirds of spiral disks in the nearby Universe. Due to gravitational interactions between the bar and the other components of galaxies, it is expected that angular momentum and matter will redistribute over long (Gyr) timescales in barred galaxies. Previous work ignoring the gas phase of galaxies h… ▽ More Elongated bar-like features are ubiquitous in galaxies, occurring at the centers of approximately two-thirds of spiral disks in the nearby Universe. Due to gravitational interactions between the bar and the other components of galaxies, it is expected that angular momentum and matter will redistribute over long (Gyr) timescales in barred galaxies. Previous work ignoring the gas phase of galaxies has conclusively demonstrated that bars should slow their rotation over time due to their interaction with dark matter halos. We have performed a simulation of a Milky Way-like galactic disk hosting a strong bar which includes a state-of-the-art model of the interstellar medium and a live dark matter halo. In this simulation the bar pattern does not slow down over time, and instead remains at a stable, constant rate of rotation. This behavior has been observed in previous simulations using more simplified models for the interstellar gas, but the apparent lack of secular evolution has remained unexplained. We find that the presence of the gas phase arrests the process by which the dark matter halo slows down a bar, a phenomenon we term bar locking. This locking is responsible for stabilizing the bar pattern speed. We find that in a Milky Way-like disk, a gas fraction of only about 5\% is necessary for this mechanism to operate. Our result naturally explains why nearly all observed bars rotate rapidly and is especially relevant for our understanding of how the Milky Way arrived at its present state. △ Less

Submitted 4 June, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 20 pages, 15 figures, submitted to ApJ; updated in response to referee comments; comments welcome

The Poor Old Heart of the Milky Way

Authors: Hans-Walter Rix, Vedant Chandra, René Andrae, Adrian M. Price-Whelan, David H. Weinberg, Charlie Conroy, Morgan Fouesneau, David W. Hogg, Francesca De Angeli, Rohan P. Naidu, Maosheng Xiang, Daniela Ruz-Mieres

Abstract: Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components 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 pr… ▽ More Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components 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 DR3, especially the XP spectra, to construct a sample of 2 million bright (BP $<15.5$ mag) giant stars within $30^\circ$ of the Galactic Center with robust [M/H] estimates, $δ$ [M/H] $\lesssim 0.1$. For most sample members we can calculate orbits based on Gaia RVS velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes $\sim 18,000$ stars with $-2.7<$ [M/H] $<-1.5$, representing a stellar mass of $\gtrsim 5\times 10^7$ M$_\odot$. The spatial distribution of these [M/H] $<-1.5$ stars has a Gaussian extent of only $σ_{\mathrm{R_{GC}}} \sim 2.7$ kpc around the Galactic center, with most of these orbits being confined to the inner Galaxy. At high orbital eccentricities, there is clear evidence for accreted halo stars in their pericentral orbit phase. Stars with [M/H] $< -2$ show no net rotation, whereas those with [M/H] $\sim -1$ are rotation dominated. Most of the tightly bound stars show $[α/\text{Fe}]$-enhancement and [Al/Fe]-[Mn/Fe] abundance patterns expected for an origin in the more massive portions of the proto-Galaxy. These central, metal-poor stars most likely predate the oldest part of the disk ($τ_{\text{age}}\approx 12.5$ Gyrs), which implies that they formed at $z\gtrsim 5$, forging the proto-Milky Way. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 21 pages, 9 figures, submitted to ApJ

Magnetic braking saturates: evidence from the orbital period distribution of low-mass detached eclipsing binaries from ZTF

Authors: Kareem El-Badry, Charlie Conroy, Jim Fuller, Rocio Kiman, Jan van Roestel, Antonio C. Rodriguez, Kevin B. Burdge

Abstract: We constrain the orbital period ($P_{\rm orb}$) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods ($P_{\rm orb}\lesssim 2$ days), binaries are predicted to evolve significantly due to magnetic braking (MB), which shr… ▽ More We constrain the orbital period ($P_{\rm orb}$) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods ($P_{\rm orb}\lesssim 2$ days), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass ($0.1\lesssim M_1/M_{\odot} < 0.9$) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$), from $P_{\rm orb}=10$ days down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods, with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirm that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer leads to significant additional angular momentum loss in CVs. △ Less

Submitted 11 October, 2022; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: Table 1 included in source files. Accepted to MNRAS

The Stellar Halo of the Galaxy is Tilted & Doubly Broken

Authors: Jiwon Jesse Han, Charlie Conroy, Benjamin D. Johnson, Joshua S. Speagle, Ana Bonaca, Vedant Chandra, Rohan P. Naidu, Yuan-Sen Ting, Turner Woody, Dennis Zaritsky

Abstract: Modern Galactic surveys have revealed an ancient merger that dominates the stellar halo of our Galaxy (\textit{Gaia}-Sausage-Enceladus, GSE). Using chemical abundances and kinematics from the H3 Survey, we identify 5559 halo stars from this merger in the radial range $r_{\text{Gal}}=6-60\text{ kpc}$. We forward model the full selection function of H3 to infer the density profile of this accreted c… ▽ More Modern Galactic surveys have revealed an ancient merger that dominates the stellar halo of our Galaxy (\textit{Gaia}-Sausage-Enceladus, GSE). Using chemical abundances and kinematics from the H3 Survey, we identify 5559 halo stars from this merger in the radial range $r_{\text{Gal}}=6-60\text{ kpc}$. We forward model the full selection function of H3 to infer the density profile of this accreted component of the stellar halo. We consider a general ellipsoid with principal axes allowed to rotate with respect to the Galactocentric axes, coupled with a multiply-broken power law. The best-fit model is a triaxial ellipsoid (axes ratios 10:8:7) tilted $25^\circ$ above the Galactic plane towards the Sun and a doubly-broken power law with breaking radii at 12 kpc and 28 kpc. This result resolves the long-standing dichotomy in literature values of the halo breaking radius, being at either $\sim15\text{ kpc}$ or $\sim30\text{ kpc}$ assuming a singly-broken power law. N-body simulations suggest that the breaking radii are connected to apocenter pile-ups of stellar orbits, and so the observed double-break provides new insight into the initial conditions and evolution of the GSE merger. Furthermore, the tilt and triaxiality of the stellar halo could imply that a fraction of the underlying dark matter halo is also tilted and triaxial. This has important implications for dynamical mass modeling of the Galaxy as well as direct dark matter detection experiments. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: Submitted to ApJ; comments warmly welcome

Schrodinger's Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

Authors: Rohan P. Naidu, Pascal A. Oesch, David J. Setton, Jorryt Matthee, Charlie Conroy, Benjamin D. Johnson, John R. Weaver, Rychard J. Bouwens, Gabriel B. Brammer, Pratika Dayal, Garth D. Illingworth, Laia Barrufet, Sirio Belli, Rachel Bezanson, Sownak Bose, Kasper E. Heintz, Joel Leja, Ecaterina Leonova, Rui Marques-Chaves, Mauro Stefanon, Sune Toft, Arjen van der Wel, Pieter van Dokkum, Andrea Weibel, Katherine E. Whitaker

Abstract: $JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}… ▽ More $JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $Λ$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2μ$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($Δz\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

A Ghost in Boötes: The Least Luminous Disrupted Dwarf Galaxy

Authors: Vedant Chandra, Charlie Conroy, Nelson Caldwell, Ana Bonaca, Rohan P. Naidu, Dennis Zaritsky, Phillip A. Cargile, Jiwon Jesse Han, Benjamin D. Johnson, Joshua S. Speagle, Yuan-Sen Ting, Turner Woody

Abstract: We report the discovery of Specter, a disrupted ultrafaint dwarf galaxy revealed by the H3 Spectroscopic Survey. We detected this structure via a pair of comoving metal-poor stars at a distance of 12.5 kpc, and further characterized it with Gaia astrometry and follow-up spectroscopy. Specter is a $25^\circ \times 1^\circ$ stream of stars that is entirely invisible until strict kinematic cuts are a… ▽ More We report the discovery of Specter, a disrupted ultrafaint dwarf galaxy revealed by the H3 Spectroscopic Survey. We detected this structure via a pair of comoving metal-poor stars at a distance of 12.5 kpc, and further characterized it with Gaia astrometry and follow-up spectroscopy. Specter is a $25^\circ \times 1^\circ$ stream of stars that is entirely invisible until strict kinematic cuts are applied to remove the Galactic foreground. The spectroscopic members suggest a stellar age $τ\gtrsim 12$ Gyr and a mean metallicity $\langle\text{[Fe/H]}\rangle = -1.84_{-0.18}^{+0.16}$, with a significant intrinsic metallicity dispersion $σ_{ \text{[Fe/H]}} = 0.37_{-0.13}^{+0.21}$. We therefore argue that Specter is the disrupted remnant of an ancient dwarf galaxy. With an integrated luminosity $M_{\text{V}} \approx -2.6$, Specter is by far the least-luminous dwarf galaxy stream known. We estimate that dozens of similar streams are lurking below the detection threshold of current search techniques, and conclude that spectroscopic surveys offer a novel means to identify extremely low surface brightness structures. △ Less

Submitted 29 November, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 17 pages, 10 figures. Accepted to ApJ

Two Remarkably Luminous Galaxy Candidates at $z\approx10-12$ Revealed by JWST

Authors: Rohan P. Naidu, Pascal A. Oesch, Pieter van Dokkum, Erica J. Nelson, Katherine A. Suess, Gabriel Brammer, Katherine E. Whitaker, Garth Illingworth, Rychard Bouwens, Sandro Tacchella, Jorryt Matthee, Natalie Allen, Rachel Bezanson, Charlie Conroy, Ivo Labbe, Joel Leja, Ecaterina Leonova, Dan Magee, Sedona H. Price, David J. Setton, Victoria Strait, Mauro Stefanon, Sune Toft, John R. Weaver, Andrea Weibel

Abstract: The first few hundred Myrs at $z>10$ mark the last major uncharted epoch in the history of the Universe, where only a single galaxy (GNz11 at $z\approx11$) is currently spectroscopically confirmed. Here we present a search for luminous $z>10$ galaxies with $JWST$/NIRCam photometry spanning $\approx1-5μ$m and covering 49 arcmin$^{2}$ from the public Early Release Science programs (CEERS and GLASS).… ▽ More The first few hundred Myrs at $z>10$ mark the last major uncharted epoch in the history of the Universe, where only a single galaxy (GNz11 at $z\approx11$) is currently spectroscopically confirmed. Here we present a search for luminous $z>10$ galaxies with $JWST$/NIRCam photometry spanning $\approx1-5μ$m and covering 49 arcmin$^{2}$ from the public Early Release Science programs (CEERS and GLASS). Our most secure candidates are two $M_{\rm{UV}}\approx-21$ systems: GLASS-z12 and GLASS-z10. These galaxies display abrupt $\gtrsim1.8$ mag breaks in their spectral energy distributions, consistent with complete absorption of flux bluewards of Lyman-$α$ that is redshifted to $z=12.4^{+0.1}_{-0.3}$ and $z=10.4^{+0.4}_{-0.5}$. Lower redshift interlopers such as quiescent galaxies with strong Balmer breaks would be comfortably detected at $>5σ$ in multiple bands where instead we find no flux. From SED modeling we infer that these galaxies have already built up $\sim 10^9$ solar masses in stars over the $\lesssim300-400$ Myrs after the Big Bang. The brightness of these sources enable morphological constraints. Tantalizingly, GLASS-z10 shows a clearly extended exponential light profile, potentially consistent with a disk galaxy of $r_{\rm{50}}\approx0.7$ kpc. These sources, if confirmed, join GNz11 in defying number density forecasts for luminous galaxies based on Schechter UV luminosity functions, which require a survey area $>10\times$ larger than we have studied here to find such luminous sources at such high redshifts. They extend evidence from lower redshifts for little or no evolution in the bright end of the UV luminosity function into the cosmic dawn epoch, with implications for just how early these galaxies began forming. This, in turn, suggests that future deep $JWST$ observations may identify relatively bright galaxies to much earlier epochs than might have been anticipated. △ Less

Submitted 25 October, 2022; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: Accepted for publication in ApJL. Analysis updated with in-flight NIRCam chip-wise zero-points -- Figure 6 shows a comparison with v1

Monochromatic globular clusters as a critical test of formation models for the dark matter deficient galaxies NGC1052-DF2 and NGC1052-DF4

Authors: Pieter van Dokkum, Zili Shen, Aaron J. Romanowsky, Roberto Abraham, Charlie Conroy, Shany Danieli, Dhruba Dutta Chowdhury, Michael A. Keim, J. M. Diederik Kruijssen, Joel Leja, Sebastian Trujillo-Gomez

Abstract: It was recently proposed that the dark matter-deficient ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group could be the products of a "bullet dwarf" collision between two gas-rich progenitor galaxies. In this model DF2 and DF4 formed at the same time in the immediate aftermath of the collision, and a strong prediction is that their globular clusters should have nearly identical stellar popula… ▽ More It was recently proposed that the dark matter-deficient ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group could be the products of a "bullet dwarf" collision between two gas-rich progenitor galaxies. In this model DF2 and DF4 formed at the same time in the immediate aftermath of the collision, and a strong prediction is that their globular clusters should have nearly identical stellar populations. Here we test this prediction by measuring accurate F606W-F814W colors from deep HST/ACS imaging. We find that the clusters are extremely homogeneous. The mean color difference between the globular clusters in DF2 and DF4 is $-0.003\pm 0.005$ mag and the observed scatter for the combined sample of 18 clusters with $M_V<-8.6$ in both galaxies is $0.015 \pm 0.002$ mag. After accounting for observational uncertainties and stochastic cluster-to-cluster variation in the number of red giants, the remaining scatter is $0.008^{+0.005}_{-0.006}$ mag. Both the color difference and the scatter are an order of magnitude smaller than in other dwarf galaxies, and we infer that the bullet scenario passes an important test that could have falsified it. No other formation models have predicted this extreme uniformity of the globular clusters in the two galaxies. We find that the galaxies themselves are slightly redder than the clusters, consistent with a previously-measured metallicity difference. Numerical simulations have shown that such differences are expected in the bullet scenario, as the galaxies continued to self-enrich after the formation of the globular clusters. △ Less

Submitted 21 December, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: Published in ApJ Letters. It is difficult to visualize the uniformity of these strange clusters. Our best attempt is Fig. 2, where we compare them to globular clusters in Virgo galaxies

Journal ref: ApJ, 940, L9 (2022)