Showing 1–50 of 194 results for author: Boylan-Kolchin, M

Stellar Metallicities and Gradients in the Faint M31 Satellites Andromeda XVI and Andromeda XXVIII

Authors: Sal Wanying Fu, Daniel R. Weisz, Else Starkenburg, Nicolas Martin, Michelle L. M. Collins, Alessandro Savino, Michael Boylan-Kolchin, Patrick Côté, Andrew E. Dolphin, Nicolas Longeard, Mario L. Mateo, Francisco J. Mercado, Nathan R. Sandford, Evan D. Skillman

Abstract: We present $\sim300$ stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI ($M_V = -7.5$) and Andromeda XXVIII ($M_V = -8.8$) derived using metallicity-sensitive Calcium H & K narrow-band Hubble Space Telescope imaging. These are the first individual stellar metallicities in And~XVI (95 stars). Our And~XXVIII sample (191 stars) is a factor of $\sim15$ increase over litera… ▽ More We present $\sim300$ stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI ($M_V = -7.5$) and Andromeda XXVIII ($M_V = -8.8$) derived using metallicity-sensitive Calcium H & K narrow-band Hubble Space Telescope imaging. These are the first individual stellar metallicities in And~XVI (95 stars). Our And~XXVIII sample (191 stars) is a factor of $\sim15$ increase over literature metallicities. For And~XVI, we measure $\langle \mbox{[Fe/H]}\rangle = -2.17^{+0.05}_{-0.05}$, $σ_{\mbox{[Fe/H]}}=0.33^{+0.07}_{-0.07}$, and $\nabla_{\mbox{[Fe/H]}} = -0.23\pm0.15$ dex $R_e^{-1}$. We find that And XVI is more metal-rich than MW UFDs of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure $\langle \mbox{[Fe/H]}\rangle = -1.95^{+0.04}_{-0.04}$, $σ_{\mbox{[Fe/H]}}=0.34^{+0.07}_{-0.07}$, and $\nabla_{\mbox{[Fe/H]}} = -0.46 \pm 0.10$~dex~$R_e^{-1}$, placing it on the dwarf galaxy mass-metallicity relation. Neither galaxy has a metallicity distribution function with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint ($L \lesssim 10^6~L_{\odot}$) galaxies outside the Milky Way halo. Both galaxies' gradients are consistent with predictions from the FIRE simulations, where an age-gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal poor ($\mbox{[Fe/H]} < -3.0$) star candidates, which make up 7% of And~XVI's MDF and 6% of And~XXVIII's. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 17 pages, 5 figures, 4 tables, ApJ submitted; comments welcome!

Elevated UV luminosity density at Cosmic Dawn explained by non-evolving, weakly-mass dependent star formation efficiency

Authors: Robert Feldmann, Michael Boylan-Kolchin, James S. Bullock, Onur Çatmabacak, Claude-André Faucher-Giguère, Christopher C. Hayward, Dušan Kereš, Alexandres Lazar, Lichen Liang, Jorge Moreno, Pascal A. Oesch, Eliot Quataert, Xuejian Shen, Guochao Sun

Abstract: Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR… ▽ More Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 26 pages, 14 figures, 5 tables, submitted to MNRAS, comments welcome

Early Galaxies and Early Dark Energy: A Unified Solution to the Hubble Tension and Puzzles of Massive Bright Galaxies revealed by JWST

Authors: Xuejian Shen, Mark Vogelsberger, Michael Boylan-Kolchin, Sandro Tacchella, Rohan P. Naidu

Abstract: JWST has revealed a large population of ultra-violet (UV)-bright galaxies at $z\gtrsim 10$ and possibly overly massive galaxies at $z\gtrsim 7$, challenging standard galaxy formation models in the $Λ$CDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tensi… ▽ More JWST has revealed a large population of ultra-violet (UV)-bright galaxies at $z\gtrsim 10$ and possibly overly massive galaxies at $z\gtrsim 7$, challenging standard galaxy formation models in the $Λ$CDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $4\lesssim z \lesssim10$ in both $Λ$CDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $z\simeq 12$ based on spectroscopically confirmed galaxies exhibit no tension with the benchmark model. Photometric constraints at $12 \lesssim z\lesssim 16$ can be fully explained within EDE via either moderately increased star formation efficiencies ($ε_{\ast}\sim 3-10\%$ at $M_{\rm halo}\sim 10^{10.5}\,{\rm M}_\odot$) or enhanced UV variabilities ($σ_{\rm UV}\sim 0.8-1.3$ mag at $M_{\rm halo}\sim 10^{10.5}\,{\rm M}_\odot$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $Λ$CDM, especially at $z\gtrsim 14$, where the required $σ_{\rm UV}$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST ($z \sim 14-16$) will be crucial for differentiating modified galaxy formation physics from new cosmological physics. △ Less

Submitted 8 July, 2024; v1 submitted 21 June, 2024; originally announced June 2024.

Comments: Minor changes to references and adding a few new lines to compare in plots. Submitted to MNRAS

Confronting the Diversity Problem: The Limits of Galaxy Rotation Curves as a tool to Understand Dark Matter Profiles

Authors: Isabel S. Sands, Philip F. Hopkins, Xuejian Shen, Michael Boylan-Kolchin, James Bullock, Claude-Andre Faucher-Giguere, Francisco J. Mercado, Jorge Moreno, Lina Necib, Xiaowei Ou, Sarah Wellons, Andrew Wetzel

Abstract: While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold… ▽ More While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold Dark Matter cosmology on galactic scales. We analyze a suite of Feedback in Realistic Environments (FIRE) simulations of $10^{10}-10^{12}$ $M_\odot$ halos with standard cold dark matter, and compare the true circular velocity to rotation curve reconstructions. We find that, for galaxies with well-ordered gaseous disks, the measured rotation curve may deviate from true circular velocity by at most 10% within the radius of the disk. However, non-equilibrium behavior, non-circular motions, and non-thermal and non-kinetic stresses may cause much larger discrepancies of 50% or more. Most rotation curve reconstructions underestimate the true circular velocity, while some reconstructions transiently over-estimate it in the central few kiloparsecs due to dynamical phenomena. We further demonstrate that the features that contribute to these failures are not always visibly obvious in HI observations. If such dwarf galaxies are included in galaxy catalogs, they may give rise to the appearance of "artificial" rotation curve diversity that does not reflect the true variation in underlying dark matter profiles. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 20 pages, 10 figures

Sliding into DM: Determining the local dark matter density and speed distribution using only the local circular speed of the Galaxy

Authors: Patrick G. Staudt, James S. Bullock, Michael Boylan-Kolchin, Andrew Wetzel, Xiaowei Ou

Abstract: We use FIRE-2 zoom simulations of Milky Way size disk galaxies to derive easy-to-use relationships between the observed circular speed of the Galaxy at the Solar location, $v_\mathrm{c}$, and dark matter properties of relevance for direct detection experiments: the dark matter density, the dark matter velocity dispersion, and the speed distribution of dark matter particles near the Solar location.… ▽ More We use FIRE-2 zoom simulations of Milky Way size disk galaxies to derive easy-to-use relationships between the observed circular speed of the Galaxy at the Solar location, $v_\mathrm{c}$, and dark matter properties of relevance for direct detection experiments: the dark matter density, the dark matter velocity dispersion, and the speed distribution of dark matter particles near the Solar location. We find that both the local dark matter density and 3D velocity dispersion follow tight power laws with $v_\mathrm{c}$. Using this relation together with the observed circular speed of the Milky Way at the Solar radius, we infer the local dark matter density and velocity dispersion near the Sun to be $ρ= 0.42\pm 0.06\,\mathrm{GeV\,cm^{-3}}$ and $σ_\mathrm{3D} = 280^{+19}_{-18}\,\mathrm{km\,s^{-1}}$. We also find that the distribution of dark matter particle speeds is well-described by a modified Maxwellian with two shape parameters, both of which correlate with the observed $v_{\rm c}$. We use that modified Maxwellian to predict the speed distribution of dark matter near the Sun and find that it peaks at a most probable speed of $250\,\mathrm{km\,s^{-1}}$ and begins to truncate sharply above $470\,\mathrm{km\,s^{-1}}$. This peak speed is somewhat higher than expected from the standard halo model, and the truncation occurs well below the formal escape speed to infinity, with fewer very-high-speed particles than assumed in the standard halo model. △ Less

Submitted 11 March, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: 22 pages, 12 figures, 2 tables; added a reference in section 1, corrected sqrt(2/3) to sqrt(3/2) in section 6

Primordial Rotating Disk Composed of $\geq$15 Dense Star-Forming Clumps at Cosmic Dawn

Authors: S. Fujimoto, M. Ouchi, K. Kohno, F. Valentino, C. Giménez-Arteaga, G. B. Brammer, L. J. Furtak, M. Kohandel, M. Oguri, A. Pallottini, J. Richard, A. Zitrin, F. E. Bauer, M. Boylan-Kolchin, M. Dessauges-Zavadsky, E. Egami, S. L. Finkelstein, Z. Ma, I. Smail, D. Watson, T. A. Hutchison, J. R. Rigby, B. D. Welch, Y. Ao, L. D. Bradley , et al. (21 additional authors not shown)

Abstract: Early galaxy formation, initiated by the dark matter and gas assembly, evolves through frequent mergers and feedback processes into dynamically hot, chaotic structures. In contrast, dynamically cold, smooth rotating disks have been observed in massive evolved galaxies merely 1.4 billion years after the Big Bang, suggesting rapid morphological and dynamical evolution in the early Universe. Probing… ▽ More Early galaxy formation, initiated by the dark matter and gas assembly, evolves through frequent mergers and feedback processes into dynamically hot, chaotic structures. In contrast, dynamically cold, smooth rotating disks have been observed in massive evolved galaxies merely 1.4 billion years after the Big Bang, suggesting rapid morphological and dynamical evolution in the early Universe. Probing this evolution mechanism necessitates studies of young galaxies, yet efforts have been hindered by observational limitations in both sensitivity and spatial resolution. Here we report high-resolution observations of a strongly lensed and quintuply imaged, low-luminosity, young galaxy at $z=6.072$ (dubbed the Cosmic Grapes), 930 million years after the Big Bang. Magnified by gravitational lensing, the galaxy is resolved into at least 15 individual star-forming clumps with effective radii of $r_{\rm e}\simeq$ 10--60 parsec (pc), which dominate $\simeq$ 70\% of the galaxy's total flux. The cool gas emission unveils a smooth, underlying rotating disk characterized by a high rotational-to-random motion ratio and a gravitationally unstable state (Toomre $Q \simeq$ 0.2--0.3), with high surface gas densities comparable to local dusty starbursts with $\simeq10^{3-5}$ $M_{\odot}$/pc$^{2}$. These gas properties suggest that the numerous star-forming clumps are formed through disk instabilities with weak feedback effects. The clumpiness of the Cosmic Grapes significantly exceeds that of galaxies at later epochs and the predictions from current simulations for early galaxies. Our findings shed new light on internal galaxy substructures and their relation to the underlying dynamics and feedback mechanisms at play during their early formation phases, potentially explaining the high abundance of bright galaxies observed in the early Universe and the dark matter core-cusp problem. △ Less

Submitted 4 March, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: Submitted. 44 pages, 14 figures, 5 tables. Comments are welcome! See also the companion papers on arXiv. Valentino+2024: arXiv:2402.17845 Giménez-Arteaga+2024: arXiv:2402.17875

An attractive model: simulating fuzzy dark matter with attractive self-interactions

Authors: Connor A. Painter, Michael Boylan-Kolchin, Philip Mocz, Mark Vogelsberger

Abstract: Fuzzy Dark Matter (FDM) comprised of ultralight ($m \sim 10^{-22}~\rm{eV}$) boson particles has received significant attention as a viable alternative to Cold Dark Matter (CDM), as it approximates CDM on large scales ($\gtrsim 1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (t… ▽ More Fuzzy Dark Matter (FDM) comprised of ultralight ($m \sim 10^{-22}~\rm{eV}$) boson particles has received significant attention as a viable alternative to Cold Dark Matter (CDM), as it approximates CDM on large scales ($\gtrsim 1$ Mpc) while potentially resolving some of its small-scale problems via kiloparsec-scale quantum interference. However, the most basic FDM model, with one free parameter (the boson mass), is subject to a tension: small boson masses yield the desired cores of dwarf galaxies but underpredict structure in the Lyman-$α$ forest, while large boson masses render FDM effectively identical to CDM. This Catch-22 problem may be alleviated by considering an axion-like particle with attractive particle self-interactions. We simulate an idealized FDM halo with self-interactions parameterized by an energy decay constant $f \sim 10^{15}~\rm{GeV}$ related to the axion symmetry-breaking conjectured to solve the strong-CP problem in particle physics. We observe solitons, a hallmark of FDM, condensing within a broader halo envelope, and find that the density profile and soliton mass depend on self-interaction strength. We propose generalized formulae to extend those from previous works to include self-interactions. We also investigate a critical mass threshold predicted for strong interactions at which the soliton collapses into a compact, unresolved state. We find that the collapse happens quickly and its effects are initially contained to the central region of the halo. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 20 pages, 16 figures; to be submitted to MNRAS

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

Stellar Metallicities and Gradients in the Isolated, Quenched Low-Mass Galaxy Tucana

Authors: Sal Wanying Fu, Daniel R. Weisz, Else Starkenburg, Nicolas Martin, Francisco J. Mercado, Alessandro Savino, Michael Boylan-Kolchin, Patrick Côté, Andrew E. Dolphin, Nicolas Longeard, Mario L. Mateo, Jenna Samuel, Nathan R. Sandford

Abstract: We measure the metallicities of 374 red giant branch (RGB) stars in the isolated, quenched dwarf galaxy Tucana using Hubble Space Telescope (HST) narrow-band (F395N) Calcium H & K (CaHK) imaging. Our sample is a factor of $\sim7$ larger than what is published. Our main findings are: (i) A global metallicity distribution function (MDF) with $\langle \mbox{[Fe/H]} \rangle = -1.55 \pm 0.04$ and… ▽ More We measure the metallicities of 374 red giant branch (RGB) stars in the isolated, quenched dwarf galaxy Tucana using Hubble Space Telescope (HST) narrow-band (F395N) Calcium H & K (CaHK) imaging. Our sample is a factor of $\sim7$ larger than what is published. Our main findings are: (i) A global metallicity distribution function (MDF) with $\langle \mbox{[Fe/H]} \rangle = -1.55 \pm 0.04$ and $σ_{\mbox{[Fe/H]}}=0.54\pm0.03$; (ii) A metallicity gradient of $-0.54 \pm 0.07$ dex $R_e^{-1}$ ($-2.1 \pm 0.3$ dex kpc$^{-1}$) over the extent of our imaging ($\sim 2.5 R_e$), which is steeper than literature measurements. Our finding is consistent with predicted gradients from the publicly-available FIRE-2 simulations, in which bursty star formation creates stellar population gradients and dark matter cores; (iii) Tucana's bifurcated RGB has distinct metallicities: a blue RGB with $\langle \mbox{[Fe/H]} \rangle = -1.78 \pm 0.06$ and $σ_{\mbox{[Fe/H]}}=0.44^{+0.07}_{-0.06}$, and a red RGB with $\langle \mbox{[Fe/H]} \rangle = -1.08 \pm 0.07$ and $σ_{\mbox{[Fe/H]}}=0.42 \pm 0.06$; (iv) At fixed stellar mass, Tucana is more MR than MW satellites by $\sim 0.4$ dex, but its blue RGB is chemically comparable to the satellites. Tucana's MDF appears consistent with star-forming isolated dwarfs, though MDFs of the latter are not as well-populated; (v) $\sim2$% of Tucana's stars have $\mbox{[Fe/H]} < -3$ and 20% $\mbox{[Fe/H]} > -1$. We provide a catalog for community spectroscopic follow-up. △ Less

Submitted 17 April, 2024; v1 submitted 10 December, 2023; originally announced December 2023.

Comments: Replaced with ApJ published version; 23 pages, 18 figures

The JWST Resolved Stellar Populations Early Release Science Program IV: The Star Formation History of the Local Group Galaxy WLM

Authors: Kristen. B. W. McQuinn, Max J. B. Newman, Alessandro Savino, Andrew E. Dolphin, Daniel R. Weisz, Benjamin F. Williams, Martha L. Boyer, Roger E. Cohen, Matteo Correnti, Andrew A. Cole, Marla C. Geha, Mario Gennaro, Nitya Kallivayalil, Karin M. Sandstrom, Evan D. Skillman, Jay Anderson, Alberto Bolatto, Michael Boylan-Kolchin, Christopher T. Garling, Karoline M. Gilbert, Leo Girardi, Jason S. Kalirai, Alessandro Mazzi, Giada Pastorelli, Hannah Richstein , et al. (1 additional authors not shown)

Abstract: We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galax… ▽ More We present the first star formation history (SFH) and age-metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color-magnitude diagram (CMD) reaches below the oldest main sequence turn-off with a SNR=10 at M_F090W=+4.6 mag; this is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering a SFH and AMR with the NIRCam F090W-F150W filter combination and provide validation of the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation re-ignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo~A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ~5 Gyr; the stellar mass then rapidly doubled in ~2.5 Gyr, followed by constant star formation over the last ~5 Gyr. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 27 pages, 10 figures, 1 table

Evidence for a Shallow Evolution in the Volume Densities of Massive Galaxies at $z=4$ to $8$ from CEERS

Authors: Katherine Chworowsky, Steven L. Finkelstein, Michael Boylan-Kolchin, Elizabeth J. McGrath, Kartheik G. Iyer, Casey Papovich, Mark Dickinson, Anthony J. Taylor, L. Y. Aaron Yung, Pablo Arrabal Haro, Micaela B. Bagley, Bren E. Backhaus, Rachana Bhatawdekar, Yingjie Cheng, Nikko J. Cleri, Justin W. Cole, M. C. Cooper, Luca Costantin, Avishai Dekel, Maximilien Franco, Seiji Fujimoto, Christopher C. Hayward, Benne W. Holwerda, Marc Huertas-Company, Michaela Hirschmann , et al. (14 additional authors not shown)

Abstract: We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we inc… ▽ More We analyze the evolution of massive (log$_{10}$ [$M_\star/M_\odot$] $>10$) galaxies at $z \sim$ 4--8 selected from the JWST Cosmic Evolution Early Release Science (CEERS) survey. We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam wide-field slitless spectroscopy. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei (AGN). We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for $z \sim$ 1--4. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 20 pages, 10 figures

Testing the near-far connection with FIRE simulations: inferring the stellar mass function of the proto-Local Group at z > 6 using the fossil record of present-day galaxies

Authors: Pratik J. Gandhi, Andrew Wetzel, Michael Boylan-Kolchin, Robyn E. Sanderson, Alessandro Savino, Daniel R. Weisz, Erik J. Tollerud, Guochao Sun, Claude-Andre Faucher-Giguere

Abstract: The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil r… ▽ More The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful `near-far' technique for studying the reionisation-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass (Mstar < 10^9 Msun) galaxies across 11 Milky Way- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterise their progenitors at z ~ 6 - 9, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with Mstar ~ 10^5 Msun (10^9 Msun) had ~1 (~30) progenitors at z ~ 7, and its main progenitor comprised ~100% (~50%) of the total stellar mass of all its progenitors at z ~ 7. We show that although only ~ 15% of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at z ~ 6 - 9. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just the low-mass galaxies at z = 0 and not the MW/M31 hosts correctly recovers the slope of the SMF down to Mstar ~ 10^4.5 Msun at z > 6. Thus, we validate the `near-far' approach as an unbiased tool for probing low-mass reionisation-era galaxies. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 6 figures, 3 appendix figures, and 1 table. Submitted to MNRAS

COSMOS-Web: Intrinsically Luminous z$\gtrsim$10 Galaxy Candidates Test Early Stellar Mass Assembly

Authors: Caitlin M. Casey, Hollis B. Akins, Marko Shuntov, Olivier Ilbert, Louise Paquereau, Maximilien Franco, Christopher C. Hayward, Steven L. Finkelstein, Michael Boylan-Kolchin, Brant E. Robertson, Natalie Allen, Malte Brinch, Olivia R. Cooper, Xuheng Ding, Nicole E. Drakos, Andreas L. Faisst, Seiji Fujimoto, Steven Gillman, Santosh Harish, Michaela Hirschmann, Shuowen Jin, Jeyhan S. Kartaltepe, Anton M. Koekemoer, Vasily Kokorev, Daizhong Liu , et al. (17 additional authors not shown)

Abstract: We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble… ▽ More We report the discovery of 15 exceptionally luminous $10\lesssim z\lesssim14$ candidate galaxies discovered in the first 0.28 deg$^2$ of JWST/NIRCam imaging from the COSMOS-Web Survey. These sources span rest-frame UV magnitudes of $-20.5>M_{\rm UV}>-22$, and thus constitute the most intrinsically luminous $z\gtrsim10$ candidates identified by JWST to-date. Selected via NIRCam imaging with Hubble ACS/F814W, deep ground-based observations corroborate their detection and help significantly constrain their photometric redshifts. We analyze their spectral energy distributions using multiple open-source codes and evaluate the probability of low-redshift solutions; we conclude that 12/15 (80%) are likely genuine $z\gtrsim10$ sources and 3/15 (20%) likely low-redshift contaminants. Three of our $z\sim12$ candidates push the limits of early stellar mass assembly: they have estimated stellar masses $\sim5\times10^{9}\,M_\odot$, implying an effective stellar baryon fraction of $ε_{\star}\sim0.2-0.5$, where $ε_{\star}\equiv M_{\star}/(f_{b}M_{halo})$. The assembly of such stellar reservoirs is made possible due to rapid, burst-driven star formation on timescales $<$100\,Myr where the star-formation rate may far outpace the growth of the underlying dark matter halos. This is supported by the similar volume densities inferred for $M_\star\sim10^{10}\,M_\odot$ galaxies relative to $M_\star\sim10^{9}\,M_\odot$ -- both about $10^{-6}$ Mpc$^{-3}$ -- implying they live in halos of comparable mass. At such high redshifts, the duty cycle for starbursts would be of order unity, which could cause the observed change in the shape of the UVLF from a double powerlaw to Schechter at $z\approx8$. Spectroscopic redshift confirmation and ensuing constraints of their masses will be critical to understanding how, and if, such early massive galaxies push the limits of galaxy formation in $Λ$CDM. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 30 pages, 9 figures; ApJ submitted

Hooks & Bends in the Radial Acceleration Relation: Discriminatory Tests for Dark Matter and MOND

Authors: Francisco J. Mercado, James S. Bullock, Jorge Moreno, Michael Boylan-Kolchin, Philip F. Hopkins, Andrew Wetzel, Claude-André Faucher-Giguère, Jenna Samuel

Abstract: The Radial Acceleration Relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, $a_{\rm tot}(r)$, with that accounted for by baryons at the same radius, $a_{\rm bar}(r)$. The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the Cosmological Constant +… ▽ More The Radial Acceleration Relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, $a_{\rm tot}(r)$, with that accounted for by baryons at the same radius, $a_{\rm bar}(r)$. The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the Cosmological Constant + Cold Dark Matter ($Λ$CDM) paradigm. We use zoom simulations of 20 galaxies with stellar masses of $M_{\star} \, \simeq 10^{7-11} \, M_{\odot}$ to study the RAR in the \texttt{FIRE-2} simulations. We highlight the existence of simulated galaxies with non-monotonic RAR tracks that ``hook'' down from the average relation. These hooks are challenging to explain in Modified Inertia theories of MOND, but naturally arise in all of our \lcdm-simulated galaxies that are dark-matter dominated at small radii and have feedback-induced cores in their dark matter haloes. We show, analytically and numerically, that downward hooks are expected in such cored haloes because they have non-monotonic acceleration profiles. We also extend the relation to accelerations below those traced by disc galaxy rotation curves. In this regime, our simulations exhibit ``bends'' off of the MOND-inspired extrapolation of the RAR, which, at large radii, approach $a_{\rm tot} \, \approx \, a_{\rm bar} \, /f_{\rm b}$, where $f_{\rm b}$ is the cosmic baryon fraction. Future efforts to search for these hooks and bends in real galaxies will provide interesting tests for MOND and $Λ$CDM. △ Less

Submitted 19 March, 2024; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 15 pages, 9 figures; Accepted for publication in MNRAS

Metallicity Distribution Functions of 13 Ultra-Faint Dwarf Galaxy Candidates from Hubble Space Telescope Narrowband Imaging

Authors: Sal Wanying Fu, Daniel R. Weisz, Else Starkenburg, Nicolas Martin, Alessandro Savino, Michael Boylan-Kolchin, Patrick Cote, Andrew E. Dolphin, Alexander P. Ji, Nicolas Longeard, Mario L. Mateo, Ekta Patel, Nathan R. Sandford

Abstract: We present uniformly measured stellar metallicities for 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; $M_V = -7.1$ to $-0.8$) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope (HST). This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the… ▽ More We present uniformly measured stellar metallicities for 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; $M_V = -7.1$ to $-0.8$) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope (HST). This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the number of metallicities in all known MW UFDs. We provide the first well-populated MDFs for all galaxies in this sample, with [Fe/H] ranging from -3.0 dex to -2.0 dex, and $σ_{\mbox{[Fe/H]}}$ ranging from 0.3 dex to 0.7 dex. We find a nearly constant [Fe/H] ~ -2.6 over 3 decades in luminosity (~$10^2 - 10^5 L_{\odot}$), suggesting that the mass-metallicity relationship does not hold for such faint systems. We find a larger fraction (24%) of EMP ([Fe/H]<-3.0) stars across our sample compared to the literature (14%), but note that uncertainties in our most metal-poor measurements make this an upper limit. We find 19% of stars in our UFD sample to be metal-rich ([Fe/H]>-2.0), consistent with the sum of literature spectroscopic studies. MW UFDs are known to be predominantly >13 Gyr old, meaning that all stars in our sample are truly ancient, unlike metal-poor stars in the MW, which have a range of possible ages. Our UFD metallicities are not well-matched to known streams in the MW, providing further evidence that known MW substructures are not related to UFDs. We include a catalog of our stars to encourage community follow-up studies, including priority targets for ELT-era observations. △ Less

Submitted 22 November, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 37 pages, 20 figures, ApJ accepted and published

Journal ref: ApJ 958 167 (2023)

The Absolute Age of M92

Authors: Jiaqi, Ying, Brian Chaboyer, Emily M. Boudreaux, Catherine Slaughter, Michael Boylan-Kolchin, Daniel Weisz

Abstract: The \textit{absolute age} of a simple stellar population is of fundamental interest for a wide range of applications but is difficult to measure in practice, as it requires an understanding of the uncertainties in a variety of stellar evolution processes as well as the uncertainty in the distance, reddening and composition. As a result, most studies focus only on the \textit{relative age} by assum… ▽ More The \textit{absolute age} of a simple stellar population is of fundamental interest for a wide range of applications but is difficult to measure in practice, as it requires an understanding of the uncertainties in a variety of stellar evolution processes as well as the uncertainty in the distance, reddening and composition. As a result, most studies focus only on the \textit{relative age} by assuming that stellar evolution calculations are accurate and using age determinations techniques that are relatively independent of distance and reddening. Here, we construct $20,000$ sets of theoretical isochrones through Monte Carlo simulation using the Dartmouth Stellar Evolution Program to measure the absolute age of the globular cluster M92. For each model, we vary a range of input physics used in the stellar evolution models, including opacities, nuclear reaction rates, diffusion coefficients, atmospheric boundary conditions, helium abundance, and treatment of convection. We also explore variations in the distance and reddening as well as its overall metallicity and $α$ enhancement. We generate simulated Hess diagrams around the main-sequence turn-off region from each set of isochrones and use a Voronoi binning method to fit the diagrams to HST ACS data. We find the age of M92 to be $13.80 \pm 0.75$ Gyr. The $5.4\%$ error in the absolute age is dominated by the uncertainty in the distance to M92 ($\sim 80\%$ of the error budget); of the remaining parameters, only the total metallicity, $α$ element abundance, and treatment of helium diffusion contribute significantly to the total error. △ Less

Submitted 3 June, 2023; originally announced June 2023.

Comments: 15 Pages, 14 Figures, 2 Tables; Accepted for Publication AJ

The Hubble Space Telescope Survey of M31 Satellite Galaxies II. The Star Formation Histories of Ultra-Faint Dwarf Galaxies

Authors: A. Savino, D. R. Weisz, E. D. Skillman, A. Dolphin, A. A. Cole, N. Kallivayalil, A. Wetzel, J. Anderson, G. Besla, M. Boylan-Kolchin, T. M. Brown, J. S. Bullock, M. L. M. Collins, M. C. Cooper, A. J. Deason, A. L. Dotter, M. Fardal, A. M. N. Ferguson, T. K. Fritz, M. C. Geha, K. M. Gilbert, P. Guhathakurta, R. Ibata, M. J. Irwin, M. Jeon , et al. (12 additional authors not shown)

Abstract: We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find th… ▽ More We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50\% of their stellar mass by $z=5$ (12.6~Gyr ago), similar to known UFDs around the MW, but that 10-40\% of their stellar mass formed at later times. We uncover one remarkable UFD, \A{XIII}, which formed only 10\% of its stellar mass by $z=5$, and 75\% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\times10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs. △ Less

Submitted 13 September, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures, 5 appendices, accepted for publication in ApJ

The impact of UV variability on the abundance of bright galaxies at $z \geq 9$

Authors: Xuejian Shen, Mark Vogelsberger, Michael Boylan-Kolchin, Sandro Tacchella, Rahul Kannan

Abstract: JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the $Λ$CDM cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function (UVLF) of galaxies at $z\geq 9$. Our study differentiates between the $\textit{median UV radiation yield}$… ▽ More JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the $Λ$CDM cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function (UVLF) of galaxies at $z\geq 9$. Our study differentiates between the $\textit{median UV radiation yield}$ and the $\textit{variability of UV luminosities}$ of galaxies at a fixed dark matter halo mass. We primarily focus on the latter effect, which depends on halo assembly and galaxy formation processes and can significantly increase the abundance of UV-bright galaxies due to the upscatter of galaxies in lower-mass haloes. We find that a relatively low level of variability, $σ_{\rm UV} \approx 0.75$ mag, matches the observational constraints at $z\approx 9$. However, increasingly larger $σ_{\rm UV}$ is necessary when moving to higher redshifts, reaching $σ_{\rm UV} \approx 2.0\,(2.5)\,{\rm mag}$ at $z\approx 12$ ($16$). This implied variability is consistent with expectations of physical processes in high-redshift galaxies such as bursty star formation and cycles of dust clearance. Photometric constraints from JWST at $z\gtrsim 9$ therefore can be reconciled with a standard $Λ$CDM-based galaxy formation model calibrated at lower redshifts without the need for adjustments to the median UV radiation yield. △ Less

Submitted 29 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 6 pages, 3 figures. Accepted and to be published on MNRAS

An analytic surface density profile for $Λ$CDM halos and gravitational lensing studies

Authors: Alexandres Lazar, James S. Bullock, Anna Nierenberg, Leonidas Moustakas, Michael Boylan-Kolchin

Abstract: We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NF… ▽ More We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NFW) profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, $R_{-1}$, where the log slope of surface density is equal to $-1$, and an associated surface density, $Σ_{-1}$, we can represent the expected lensing signal from line-of-sight halos statistically, for an ensemble of halo orientations, using a distribution of {\em projected concentration} parameters, $\mathcal{C}_{\rm vir} := r_{\rm vir}/ R_{-1}$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use. △ Less

Submitted 3 January, 2024; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: 15 pages, 13 figures, Accepted to MNRAS

Cosmological Structure Formation and Soliton Phase Transition in Fuzzy Dark Matter with Axion Self-Interactions

Authors: Philip Mocz, Anastasia Fialkov, Mark Vogelsberger, Michael Boylan-Kolchin, Pierre-Henri Chavanis, Mustafa A. Amin, Sownak Bose, Tibor Dome, Lars Hernquist, Lachlan Lancaster, Matthew Notis, Connor Painter, Victor H. Robles, Jesus Zavala

Abstract: We investigate cosmological structure formation in Fuzzy Dark Matter (FDM) with an attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We f… ▽ More We investigate cosmological structure formation in Fuzzy Dark Matter (FDM) with an attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 10 pages, 3 figures, submitted to mnras

What Causes The Formation of Disks and End of Bursty Star Formation?

Authors: Philip F. Hopkins, Alexander B. Gurvich, Xuejian Shen, Zachary Hafen, Michael Y. Grudic, Shalini Kurinchi-Vendhan, Christopher C. Hayward, Fangzhou Jiang, Matthew E. Orr, Andrew Wetzel, Dusan Keres, Jonathan Stern, Claude-Andre Faucher-Giguere, James Bullock, Coral Wheeler, Kareem El-Badry, Sarah R. Loebman, Jorge Moreno, Michael Boylan-Kolchin, Eliot Quataert

Abstract: As they grow, galaxies can transition from irregular/spheroidal with 'bursty' star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existin… ▽ More As they grow, galaxies can transition from irregular/spheroidal with 'bursty' star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of disks and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disk formation is promoted when the mass profile becomes sufficiently centrally-concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical center, ceases to support the global 'breathing modes' which can persist indefinitely in less-concentrated profiles and efficiently destroy disks, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disk. Smooth SF is promoted by the potential or escape velocity (not circular velocity) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disk 'settling'), and observations. △ Less

Submitted 28 August, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: Updated to match MNRAS accepted version. 45 pages, 32 figures

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

HETDEX Public Source Catalog 1: 220K Sources Including Over 50K Lyman Alpha Emitters from an Untargeted Wide-area Spectroscopic Survey

Authors: Erin Mentuch Cooper, Karl Gebhardt, Dustin Davis, Daniel J. Farrow, Chenxu Liu, Gregory Zeimann, Robin Ciardullo, John J. Feldmeier, Niv Drory, Donghui Jeong, Barbara Benda, William P. Bowman, Michael Boylan-Kolchin, Oscar A. Chavez Ortiz, Maya H. Debski, Mona Dentler, Maximilian Fabricius, Rameen Farooq, Steven L. Finkelstein, Eric Gawiser, Caryl Gronwall, Gary J. Hill, Ulrich Hopp, Lindsay R. House, Steven Janowiecki , et al. (21 additional authors not shown)

Abstract: We present the first publicly released catalog of sources obtained from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). HETDEX is an integral field spectroscopic survey designed to measure the Hubble expansion parameter and angular diameter distance at 1.88<z<3.52 by using the spatial distribution of more than a million Ly-alpha-emitting galaxies over a total target area of 540 deg^2.… ▽ More We present the first publicly released catalog of sources obtained from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). HETDEX is an integral field spectroscopic survey designed to measure the Hubble expansion parameter and angular diameter distance at 1.88<z<3.52 by using the spatial distribution of more than a million Ly-alpha-emitting galaxies over a total target area of 540 deg^2. The catalog comes from contiguous fiber spectra coverage of 25 deg^2 of sky from January 2017 through June 2020, where object detection is performed through two complementary detection methods: one designed to search for line emission and the other a search for continuum emission. The HETDEX public release catalog is dominated by emission-line galaxies and includes 51,863 Lyα-emitting galaxy (LAE) identifications and 123,891 OII-emitting galaxies at z<0.5. Also included in the catalog are 37,916 stars, 5274 low-redshift (z<0.5) galaxies without emission lines, and 4976 active galactic nuclei. The catalog provides sky coordinates, redshifts, line identifications, classification information, line fluxes, OII and Ly-alpha line luminosities where applicable, and spectra for all identified sources processed by the HETDEX detection pipeline. Extensive testing demonstrates that HETDEX redshifts agree to within deltaz < 0.02, 96.1% of the time to those in external spectroscopic catalogs. We measure the photometric counterpart fraction in deep ancillary Hyper Suprime-Cam imaging and find that only 55.5% of the LAE sample has an r-band continuum counterpart down to a limiting magnitude of r~26.2 mag (AB) indicating that an LAE search of similar sensitivity with photometric pre-selection would miss nearly half of the HETDEX LAE catalog sample. Data access and details about the catalog can be found online at http://hetdex.org/. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: 38 pages, 20 figures. Data access and details about the catalog can be found online at http://hetdex.org/. A copy of the catalogs presented in this work (Version 3.2) is available to download at Zenodo doi:10.5281/zenodo.7448504

A jolt to the system: ram pressure on low-mass galaxies in simulations of the Local Group

Authors: Jenna Samuel, Bhavya Pardasani, Andrew Wetzel, Isaiah Santistevan, Michael Boylan-Kolchin, Jorge Moreno, Claude-Andre Faucher-Giguere

Abstract: Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies ($M_{*}\sim10^{5-9}\,\rm{M}_{\odot}$) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mas… ▽ More Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies ($M_{*}\sim10^{5-9}\,\rm{M}_{\odot}$) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mass galaxies over time. The impulsiveness of ram pressure -- the maximum ram pressure scaled to the integrated ram pressure prior to quenching -- correlates with whether a galaxy is quiescent or star-forming. The time-scale between maximum ram pressure and quenching is anticorrelated with impulsiveness, such that high impulsiveness corresponds to quenching time-scales $<1$ Gyr. Galaxies in low-mass groups ($M_\mathrm{*,host}\sim10^{7-9}\,\rm{M}_{\odot}$) outside of MW haloes experience typical ram pressure only slightly lower than ram pressure on MW satellites, helping to explain effective quenching via group pre-processing. Ram pressure on MW satellites rises sharply with decreasing distance to the host, and, at a fixed physical distance, more recent pericentre passages are typically associated with higher ram pressure because of greater gas density in the inner host halo at late times. Furthermore, the ram pressure and gas density in the inner regions of Local Group-like paired host haloes is higher at small angles off the host galaxy disc compared to isolated hosts. The quiescent fraction of satellites within these low-latitude regions is also elevated in the simulations and observations, signaling possible anisotropic quenching via ram pressure around MW-mass hosts. △ Less

Submitted 25 August, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 13 pages, 10 figures, + appendices. Accepted for publication in MNRAS

Orientations of DM Halos in FIRE-2 Milky Way-mass Galaxies

Authors: Jay Baptista, Robyn Sanderson, Dan Huber, Andrew Wetzel, Omid Sameie, Michael Boylan-Kolchin, Jeremy Bailin, Philip F. Hopkins, Claude-Andre Faucher-Giguere, Sukanya Chakrabarti, Drona Vargya, Nondh Panithanpaisal, Arpit Arora, Emily Cunningham

Abstract: The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FI… ▽ More The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $σ= 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Submitted to ApJ, 19 pages, 12 figures, 3 tables

Dwarf galaxy formation with and without dark matter-baryon streaming velocities

Authors: Anna T. P. Schauer, Michael Boylan-Kolchin, Katelyn Colston, Omid Sameie, Volker Bromm, James S. Bullock, Andrew Wetzel

Abstract: We study how supersonic streaming velocities of baryons relative to dark matter -- a large-scale effect imprinted at recombination and coherent over $\sim 3$ Mpc scales -- affects the formation of dwarf galaxies at $z \gtrsim 5$. We perform cosmological hydrodynamic simulations, including and excluding streaming velocities, in regions centered on halos with $M_{\rm vir}(z=0) \approx 10^{10}$ M… ▽ More We study how supersonic streaming velocities of baryons relative to dark matter -- a large-scale effect imprinted at recombination and coherent over $\sim 3$ Mpc scales -- affects the formation of dwarf galaxies at $z \gtrsim 5$. We perform cosmological hydrodynamic simulations, including and excluding streaming velocities, in regions centered on halos with $M_{\rm vir}(z=0) \approx 10^{10}$ M$_{\odot}$; the simulations are part of the Feedback In Realistic Environments (FIRE) project and run with FIRE-3 physics. Our simulations comprise many thousands of systems with halo masses between $M_{\rm vir} = 2\times10^{5}$ M$_{\odot}$ and $2\times10^9$ M$_{\odot}$ in the redshift range $z=20-5$. A few hundred of these galaxies form stars and have stellar masses ranging from 100 to $10^7$ M$_{\odot}$. While star formation is globally delayed by approximately 50 Myr in the streaming relative to non-streaming simulations and the number of luminous galaxies is correspondingly suppressed at high redshift in the streaming runs, these effects decay with time. By $z=5$, the properties of the simulated galaxies are nearly identical in the streaming versus non-streaming runs, indicating that any effects of streaming velocities on the properties of galaxies at the mass scale of classical dwarfs and larger do not persist to $z=0$. △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: 12 pages, 8 figures, submitted to ApJ

Born this way: thin disc, thick disc, and isotropic spheroid formation in FIRE-2 Milky-Way-mass galaxy simulations

Authors: Sijie Yu, James S. Bullock, Alexander B. Gurvich, Zachary Hafen, Jonathan Stern, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Andrew Wetzel, Philip F. Hopkins, Jorge Moreno

Abstract: We investigate the formation of Milky-Way-mass galaxies using FIRE-2 LCDM cosmological zoom-in simulations by studying the orbital evolution of stars formed in the main progenitor of the galaxy, from birth to the present day. We classify in situ stars as isotropic spheroid, thick-disc, and thin-disc according to their orbital circularities and show that these components are assembled in a time-ord… ▽ More We investigate the formation of Milky-Way-mass galaxies using FIRE-2 LCDM cosmological zoom-in simulations by studying the orbital evolution of stars formed in the main progenitor of the galaxy, from birth to the present day. We classify in situ stars as isotropic spheroid, thick-disc, and thin-disc according to their orbital circularities and show that these components are assembled in a time-ordered sequence from early to late times, respectively. All simulated galaxies experience an early phase of bursty star formation that transitions to a late-time steady phase. This transition coincides with the time that the inner CGM virializes. During the early bursty phase, galaxies have irregular morphologies and new stars are born on radial orbits; these stars evolve into an isotropic spheroidal population today. The bulk of thick-disc stars form at intermediate times, during a clumpy-disc ``spin-up'' phase, slightly later than the peak of spheroid formation. At late times, once the CGM virializes and star formation ``cools down," stars are born on circular orbits within a narrow plane. Those stars mostly inhabit thin discs today. Broadly speaking, stars with disc-like or spheroid-like orbits today were born that way. Mergers onto discs and secular processes do affect kinematics in our simulations, but play only secondary roles in populating thick-disc and in situ spheroid populations at z=0. The age distributions of spheroid, thick disc, and thin disc populations scale self-similarly with the steady-phase transition time, which suggests that morphological age dating can be linked to the CGM virialization time in galaxies. △ Less

Submitted 19 June, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: 16 pages, 10 figures; this version was accepted to MNRAS on June 9, 2023

On the Cosmic Web Elongation in Fuzzy Dark Matter Cosmologies: Effects on Density Profiles, Shapes and Alignments of Halos

Authors: Tibor Dome, Anastasia Fialkov, Philip Mocz, Björn Malte Schäfer, Michael Boylan-Kolchin, Mark Vogelsberger

Abstract: The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($Λ$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like… ▽ More The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($Λ$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like power spectrum cutoff on small scales in the primordial density perturbations. We study halo density profiles, shapes and alignments in FDM-like cosmologies (the latter two for the first time) by providing fits and quantifying departures from $Λ$CDM as a function of the particle mass $m$. Compared to $Λ$CDM, the concentrations of FDM-like halos are lower, peaking at an $m$-dependent halo mass and thus breaking the approximate universality of density profiles in $Λ$CDM. The intermediate-to-major and minor-to-major shape parameter profiles are monotonically increasing with ellipsoidal radius in $N$-body simulations of $Λ$CDM. In FDM-like cosmologies, the monotonicity is broken, halos are more elongated around the virial radius than their $Λ$CDM counterparts and less elongated closer to the center. Finally, intrinsic alignment correlations, stemming from the deformation of initially spherically collapsing halos in an ambient gravitational tidal field, become stronger with decreasing $m$. At $z\sim 4$, we find a $6.4 σ$-significance in the fractional differences between the isotropised linear alignment magnitudes $D_{\text{iso}}$ in the $m=10^{-22}$ eV model and $Λ$CDM. Such FDM-like imprints on the internal properties of virialised halos are expected to be strikingly visible in the high-$z$ Universe. △ Less

Submitted 19 January, 2023; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 20 pages, 10 figures, 3 tables

Journal ref: MNRAS, 519, 4183-4202 (2022)

Stress Testing $Λ$CDM with High-redshift Galaxy Candidates

Authors: Michael Boylan-Kolchin

Abstract: Early data from JWST have revealed a bevy of high-redshift galaxy candidates with unexpectedly high stellar masses. An immediate concern is the consistency of these candidates with galaxy formation in the standard cosmological model. In the $Λ$CDM paradigm, the stellar mass ($M_\star$) of a galaxy is limited by the available baryonic reservoir of its host dark matter halo. The mass function of dar… ▽ More Early data from JWST have revealed a bevy of high-redshift galaxy candidates with unexpectedly high stellar masses. An immediate concern is the consistency of these candidates with galaxy formation in the standard cosmological model. In the $Λ$CDM paradigm, the stellar mass ($M_\star$) of a galaxy is limited by the available baryonic reservoir of its host dark matter halo. The mass function of dark matter halos therefore imposes an absolute upper limit on the number density $n(>M_\star,z)$ and stellar mass density $ρ_{\star}(>M_\star,z)$ of galaxies more massive than $M_\star$ at any epoch $z$. Here I show that the most massive galaxy candidates in JWST observations at $z\sim 7-10$ lie at the very edge of these limits, indicating an important unresolved issue with the properties of galaxies derived from the observations, how galaxies form at early times in $Λ$CDM, or within this standard cosmology itself. △ Less

Submitted 20 February, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 4 pages, 2 figures; author's version of paper to appear in Nature Astronomy. V2: revised with updated observational data on most massive galaxy candidates at $z \sim 7-10$

Journal ref: Nature Astronomy 7 (2023) 731-735

On the probability of the extremely lensed $z$=6.2 Earendel source being a Population~III star

Authors: Anna T. P. Schauer, Volker Bromm, Niv Drory, Michael Boylan-Kolchin

Abstract: The recent discovery of the extremely lensed Earendel object at $z=6.2$ is remarkable in that it is likely a single star or stellar multiple, observed within the first billion years of cosmic history. Depending on its mass, which is still uncertain but will soon be more tightly constrained with the James Webb Space Telescope, the Earendel star might even be a member of the first generation of star… ▽ More The recent discovery of the extremely lensed Earendel object at $z=6.2$ is remarkable in that it is likely a single star or stellar multiple, observed within the first billion years of cosmic history. Depending on its mass, which is still uncertain but will soon be more tightly constrained with the James Webb Space Telescope, the Earendel star might even be a member of the first generation of stars, the so-called Population~III (Pop~III). By combining results from detailed cosmological simulations of the assembly of the first galaxies, including the enrichment of the pristine gas with heavy chemical elements, with assumptions on key stellar parameters, we quantify the probability that Earendel has indeed a Pop~III origin. We find that this probability is non-negligible throughout the mass range inferred for Earendel, specifically ranging from a few percent at the lower-mass end to near unity for some Pop~III initial mass function (IMF) models towards the high-mass end of the allowed range. For models that extend the metal-enriched IMF to $500$\,M$_\odot$, the likelihood of Earendel being a Pop~III star stays at the few to ten percent level. We discuss the implications of such a discovery for the overall endeavor to probe the hitherto so elusive first stars in the universe. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: 5 pages, 2 figures, resubmitted to ApJL

Dissipative Dark Matter on FIRE: II. Observational signatures and constraints from local dwarf galaxies

Authors: Xuejian Shen, Philip F. Hopkins, Lina Necib, Fangzhou Jiang, Michael Boylan-Kolchin, Andrew Wetzel

Abstract: We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with… ▽ More We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with $M_{\ast} \sim 10^{5}\operatorname{-}10^{9}\,{\rm M}_{\odot}$ in both isolation and within Milky Way-mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of stellar/neutral gas disk formation in dSIDM with $(σ/m)\leq 1\,{\rm cm^2\,g^{-1}}$ but they are still consistent with observed galaxy sizes and masses. In addition, as a result of the steeper central density profiles developed in dSIDM, the sub-kpc circular velocities of isolated dwarfs in models with $(σ/m)\geq 0.1\,{\rm cm^2\,g^{-1}}$ are enhanced by about a factor of two, which are still consistent with the measured stellar velocity dispersions of Local Group dwarfs but in tension with the HI rotation curves of more massive field dwarfs. Meanwhile, for satellites of the simulated Milky Way-mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help address the missing compact dwarf satellites in CDM. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with constant cross-section $(σ/m) \gtrsim 0.1\,{\rm cm^2\,g^{-1}}$ (assuming $f_{\rm diss}=0.5$) are effectively ruled out in bright dwarfs ($M_{\rm halo}\sim 10^{11}\,{\rm M}_{\odot}$) by circular velocity constraints. However, models with lower effective cross-sections (at this halo mass/velocity scale) are still viable and can give rise to non-trivial observable signatures. △ Less

Submitted 15 June, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: To be submitted to MNRAS. 18 pages, 12 figures

The Hubble Space Telescope Survey of M31 Satellite Galaxies I. RR Lyrae-based Distances and Refined 3D Geometric Structure

Authors: Alessandro Savino, Daniel R. Weisz, Evan D. Skillman, Andrew Dolphin, Nitya Kallivayalil, Andrew Wetzel, Jay Anderson, Gurtina Besla, Michael Boylan-Kolchin, James S. Bullock, Andrew A. Cole, Michelle L. M. Collins, M. C. Cooper, Alis J. Deason, Aaron L. Dotter, Mark Fardal, Annette M. N. Ferguson, Tobias K. Fritz, Marla C. Geha, Karoline M. Gilbert, Puragra Guhathakurta, Rodrigo Ibata, Michael J. Irwin, Myoungwon Jeon, Evan Kirby , et al. (11 additional authors not shown)

Abstract: We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days… ▽ More We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system. △ Less

Submitted 12 September, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 32 pages, 15 figures, accepted for publication on ApJ

Formation of proto-globular cluster candidates in cosmological simulations of dwarf galaxies at $z>4$

Authors: Omid Sameie, Michael Boylan-Kolchin, Philip F. Hopkins, Andrew Wetzel, Xiangcheng Ma, James S. Bullock, Kareem El-Badry, Eliot Quataert, Jenna Samuel, Anna T. P. Schauer, Daniel R. Weisz

Abstract: We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at $z=0$) as part of the "Feedback in Realistic Environment" (FIRE) project. Compact ($r_{1/2}<30$ pc), relatively massive ($0.5 \times 10^5 \lesssim M_{\star}/{\rm M}_\odot \lesssim 5\times10^5$),… ▽ More We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at $z=0$) as part of the "Feedback in Realistic Environment" (FIRE) project. Compact ($r_{1/2}<30$ pc), relatively massive ($0.5 \times 10^5 \lesssim M_{\star}/{\rm M}_\odot \lesssim 5\times10^5$), self-bound stellar clusters form at $11\gtrsim z \gtrsim 5$ in progenitors with $M_{\rm vir} \approx 10^9\,{\rm M}_\odot$. Cluster formation is triggered when at least $10^7\,{\rm M}_\odot$ of dense, turbulent gas reaches $Σ_{\rm gas} \approx 10^4\, {\rm M}_\odot\, {\rm pc}^{-2}$ as a result of the compressive effects of supernova feedback or from cloud-cloud collisions. The clusters can survive for $2-3\,{\rm Gyr}$; absent numerical effects, they would likely survive substantially longer, perhaps to $z=0$. The longest-lived clusters are those that form at significant distance -- several hundreds of pc -- from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter halos as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in \textit{Hubble Space Telescope} lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with the \textit{James Webb Space Telescope}. By contrast, the star clusters' host galaxies will remain undetectable. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 14 pages, 5 figures, submitted to MNRAS

Extinguishing the FIRE: environmental quenching of satellite galaxies around Milky Way-mass hosts in simulations

Authors: Jenna Samuel, Andrew Wetzel, Isaiah Santistevan, Erik Tollerud, Jorge Moreno, Michael Boylan-Kolchin, Jeremy Bailin, Bhavya Pardasani

Abstract: The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in $z=0$ galaxies at $M*=10^{5-10}\rm{M}_{\odot}$ around 14 MW-mass hosts from the FIRE-2 simulations. Lower-mass s… ▽ More The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in $z=0$ galaxies at $M*=10^{5-10}\rm{M}_{\odot}$ around 14 MW-mass hosts from the FIRE-2 simulations. Lower-mass satellites ($M*\lesssim10^7\rm{M}_{\odot}$) are mostly quiescent and higher-mass satellites ($M*\gtrsim10^8\rm{M}_{\odot}$) are mostly star-forming, with intermediate-mass satellites ($M*\approx10^{7-8}\rm{M}_{\odot}$) split roughly equally between quiescent and star-forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M*=10^{8-9}\rm{M}_{\odot}$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and LMC-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as preprocessing in low-mass groups prior to MW infall. Lower-mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced $\geq1-2$ pericentre passages ($\approx2.5-5$ Gyr) within a MW-mass halo. Most galaxies with $M*\gtrsim10^{6.5}\rm{M}_{\odot}$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the SAGA survey's quiescent fraction at $M*\gtrsim10^8\rm{M}_{\odot}$, they do not generally reproduce SAGA's turnover at lower masses. △ Less

Submitted 17 June, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 17 pages, 13 figures, + appendices. Accepted by MNRAS

Snowmass2021 Cosmic Frontier: Synergies between dark matter searches and multiwavelength/multimessenger astrophysics

Authors: Shin'ichiro Ando, Sebastian Baum, Michael Boylan-Kolchin, Esra Bulbul, Michael Burgess, Ilias Cholis, Philip von Doetinchem, JiJi Fan, Patrick J. Harding, Shunsaku Horiuchi, Rebecca K. Leane, Oscar Macias, Katie Mack, Kohta Murase, Lina Necib, Ibles Olcina, Laura Olivera-Nieto, Jong-Chul Park, Kerstin Perez, Marco Regis, Nicholas L. Rodd, Carsten Rott, Kuver Sinha, Volodymyr Takhistov, Yun-Tse Tsai , et al. (1 additional authors not shown)

Abstract: This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dar… ▽ More This whitepaper focuses on the astrophysical systematics which are encountered in dark matter searches. Oftentimes in indirect and also in direct dark matter searches, astrophysical systematics are a major limiting factor to sensitivity to dark matter. Just as there are many forms of dark matter searches, there are many forms of backgrounds. We attempt to cover the major systematics arising in dark matter searches using photons -- radio and gamma rays -- to cosmic rays, neutrinos and gravitational waves. Examples include astrophysical sources of cosmic messengers and their interactions which can mimic dark matter signatures. In turn, these depend on commensurate studies in understanding the cosmic environment -- gas distributions, magnetic field configurations -- as well as relevant nuclear astrophysics. We also cover the astrophysics governing celestial bodies and galaxies used to probe dark matter, from black holes to dwarf galaxies. Finally, we cover astrophysical backgrounds related to probing the dark matter distribution and kinematics, which impact a wide range of dark matter studies. In the future, the rise of multi-messenger astronomy, and novel analysis methods to exploit it for dark matter, will offer various strategic ways to continue to enhance our understanding of astrophysical backgrounds to deliver improved sensitivity to dark matter. △ Less

Submitted 13 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021

Great Balls of FIRE I: The formation of star clusters across cosmic time in a Milky Way-mass galaxy

Authors: Michael Y. Grudić, Zachary Hafen, Carl L. Rodriguez, Dávid Guszejnov, Astrid Lamberts, Andrew Wetzel, Michael Boylan-Kolchin, Claude-André Faucher-Giguère

Abstract: The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium (ISM) conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modeling cluster formation in galaxy simulations: mapping giant molecular… ▽ More The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium (ISM) conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modeling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 MHD galaxy simulation onto a cluster population according to a GMC-scale cluster formation model calibrated to higher-resolution simulations, obtaining detailed properties of the galaxy's star clusters in mass, metallicity, space, and time. We find $\sim 10\%$ of all stars formed in the galaxy originate in gravitationally-bound clusters overall, and this fraction increases in regions with elevated $Σ_{\rm gas}$ and $Σ_{\rm SFR}$, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies -- no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as $7\times 10^6 M_\odot$ form in massive, dense clouds with high star formation efficiency. The initial mass-radius relation of young star clusters is consistent with an environmentally-dependent 3D density that increases with $Σ_{\rm gas}$ and $Σ_{\rm SFR}$. The model does not reproduce the age and metallicity statistics of old ($>11\rm Gyr$) globular clusters found in the Milky Way, possibly because it forms stars more slowly at $z>3$. △ Less

Submitted 1 December, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: Accepted to MNRAS

Snowmass2021 CMB-HD White Paper

Authors: The CMB-HD Collaboration, :, Simone Aiola, Yashar Akrami, Kaustuv Basu, Michael Boylan-Kolchin, Thejs Brinckmann, Sean Bryan, Caitlin M. Casey, Jens Chluba, Sebastien Clesse, Francis-Yan Cyr-Racine, Luca Di Mascolo, Simon Dicker, Thomas Essinger-Hileman, Gerrit S. Farren, Michael A. Fedderke, Simone Ferraro, George M. Fuller, Nicholas Galitzki, Vera Gluscevic, Daniel Grin, Dongwon Han, Matthew Hasselfield, Renee Hlozek , et al. (40 additional authors not shown)

Abstract: CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter… ▽ More CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. Note some text overlap with CMB-HD Astro2020 APC and RFI (arXiv:1906.10134, arXiv:2002.12714). Science case further broadened and updated

FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations

Authors: Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

Abstract: Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numeric… ▽ More Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation

Authors: Andrew Wetzel, Christopher C. Hayward, Robyn E. Sanderson, Xiangcheng Ma, Daniel Angles-Alcazar, Robert Feldmann, T. K Chan, Kareem El-Badry, Coral Wheeler, Shea Garrison-Kimmel, Farnik Nikakhtar, Nondh Panithanpaisal, Arpit Arora, Alexander B. Gurvich, Jenna Samuel, Omid Sameie, Viraj Pandya, Zachary Hafen, Cameron Hummels, Sarah Loebman, Michael Boylan-Kolchin, James S. Bullock, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert , et al. (1 additional authors not shown)

Abstract: We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, includ… ▽ More We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. △ Less

Submitted 29 March, 2023; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

Authors: Jorge Moreno, Shany Danieli, James S. Bullock, Robert Feldmann, Philip F. Hopkins, Onur Catmabacak, Alexander Gurvich, Alexandres Lazar, Courtney Klein, Cameron B. Hummels, Zachary Hafen, Francisco J. Mercado, Sijie Yu, Fangzhou Jiang, Coral Wheeler, Andrew Wetzel, Daniel Angles-Alcazar, Michael Boylan-Kolchin, Eliot Quataert, Claude-Andre Faucher-Giguere, Dusan Keres

Abstract: The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark ma… ▽ More The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 55 pages, 4 figures, 13 supplementary figures. 29 pages, 4 figures. Accepted for publication in Nature Astronomy. To appear on 14-February-2022. Published version: https://www.nature.com/articles/s41550-021-01598-4

Hot-mode accretion and the physics of thin-disk galaxy formation

Authors: Zachary Hafen, Jonathan Stern, James Bullock, Alex B. Gurvich, Sijie Yu, Claude-Andre Faucher-Giguere, Drummond B. Fielding, Daniel Angles-Alcazar, Eliot Quataert, Andrew Wetzel, Tjitske Starkenburg, Michael Boylan-Kolchin, Jorge Moreno, Robert Feldmann, Kareem El-Badry, T. K. Chan, Cameron Trapp, Dusan Keres, Philip F. Hopkins

Abstract: We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that:… ▽ More We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe. △ Less

Submitted 6 June, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 17 pages, 11 figures. Comments welcome. Accompanying interactive visualization at zhafen.github.io/rotating-cooling-flows

The galaxy-halo size relation of low-mass galaxies in FIRE

Authors: Eric Rohr, Robert Feldmann, James Bullock, Onur Çatmabacak, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Dušan Kereš, Lichen Liang, Jorge Moreno, Andrew Wetzel

Abstract: Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy-halo size relation for low-mass ($M_\star \sim 10^{7-9} {\rm M_\odot}$) central galaxies over the past 12.5 b… ▽ More Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy-halo size relation for low-mass ($M_\star \sim 10^{7-9} {\rm M_\odot}$) central galaxies over the past 12.5 billion years with the help of cosmological volume simulations (FIREbox) from the Feedback in Realistic Environments (FIRE) project. We find a nearly linear relationship between the half-stellar mass galaxy size $R_{1/2}$ and the parent dark matter halo virial radius $R_{\rm vir}$. This relation evolves only weakly since redshift $z = 5$: $R_{1/2} {\rm kpc} = (0.053\pm0.002)(R_{\rm vir}/35 {\rm kpc})^{0.934\pm0.054}$, with a nearly constant scatter $\langle σ\rangle = 0.084 [{\rm dex}]$. Whilst this ratio is similar to what is expected from models where galaxy disc sizes are set by halo angular momentum, the low-mass galaxies in our sample are not angular momentum supported, with stellar rotational to circular velocity ratios $v_{\rm rot} / v_{\rm circ} \sim 0.15$. Introducing redshift as another parameter to the GHSR does not decrease the scatter. Furthermore, this scatter does not correlate with any of the halo properties we investigate -- including spin and concentration -- suggesting that baryonic processes and feedback physics are instead critical in setting the scatter in the galaxy-halo size relation. Given the relatively small scatter and the weak dependence of the galaxy-halo size relation on redshift and halo properties for these low-mass central galaxies, we propose using galaxy sizes as an independent method from stellar masses to infer halo masses. △ Less

Submitted 16 April, 2024; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 14 pages, 9 figures, 3 tables + 3 appendices (5 pages, 5 figures, 2 tables). Accepted for publication in MNRAS. Reuploaded after proofing stage

Journal ref: MNRAS, 510, 3967-3985 (2022)

Amplified J-factors in the Galactic Center for velocity-dependent darkmatter annihilation in FIRE simulations

Authors: Daniel McKeown, James S. Bullock, Francisco J. Mercado, Zachary Hafen, Michael Boylan-Kolchin, Andrew Wetzel, Lina Necib, Philip F. Hopkins, Sijie Yu

Abstract: We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross sections $σ_v$, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross section is either either p-wave (… ▽ More We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross sections $σ_v$, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross section is either either p-wave ($\propto v^2/c^2$) or d-wave ($\propto v^4/c^4$). We use 12 pairs of simulations, each run with dark-matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of $\sim 2.5-4$. They also have a larger range of central ($\sim 400$ pc) dark matter densities than the DMO runs ($ρ_{\rm FIRE}/ρ_{\rm DMO} \simeq 0.5 - 3$) owing to the competing effects of baryonic contraction and feedback. At 3 degrees from the Galactic Center, FIRE J-factors are $5-50$ (p-wave) and $15-500$ (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 degrees is more modest and can be higher or lower ($\sim 0.3-6$), though the shape of the emission profile is flatter (less peaked towards the Galactic Center) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observe p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future. △ Less

Submitted 14 April, 2022; v1 submitted 4 November, 2021; originally announced November 2021.

Comments: 15 pages, 9 figures, submitted to MNRAS

Metallicity Distribution Function of the Eridanus~II Ultra-Faint Dwarf Galaxy from Hubble Space Telescope Narrow-band Imaging

Authors: Sal Wanying Fu, Daniel R. Weisz, Else Starkenburg, Nicolas Martin, Alexander P. Ji, Ekta Patel, Michael Boylan-Kolchin, Patrick Cote, Andrew E. Dolphin, Nicolas Longeard, Mario L. Mateo, Nathan R. Sandford

Abstract: We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more… ▽ More We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more stars than current spectroscopic MDF determinations. We find that Eri II has a mean metallicity of [Fe/H]$=$-2.50$^{+0.07}_{-0.07}$ and a dispersion of $σ_{\mbox{[Fe/H]}}=0.42^{+0.06}_{-0.06}$, which are consistent with spectroscopic MDFs, though more precisely constrained owing to a larger sample. We identify a handful of extremely metal-poor star candidates (EMP; [Fe/H] $< -3$) that are marginally bright enough for spectroscopic follow up. Eri II's MDF appears well-described by a leaky box chemical evolution model. We also compute an updated orbital history for Eri II using Gaia eDR3 proper motions, and find that it is likely on first infall into the Milky Way. Our findings suggest that Eri II underwent an evolutionary history similar to that of an isolated galaxy. Compared to MDFs for select cosmological simulations of similar mass galaxies, we find that Eri II has a lower fraction of stars with [Fe/H] $< -3$, though such comparisons should currently be treated with caution due to a paucity of simulations, selection effects, and known limitations of Ca H&K for EMPs. This study demonstrates the power of deep HST CaHK imaging for measuring the MDFs of UFDs. △ Less

Submitted 29 October, 2021; originally announced November 2021.

Comments: Resubmitted to ApJ with revisions following a positive referee report of the initial draft; 29 pages, 19 figures, and 3 tables, including appendix. Comments welcome!

From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

Authors: Mauro Bernardini, Robert Feldmann, Daniel Anglés-Alcázar, Mike Boylan-Kolchin, James Bullock, Lucio Mayer, Joachim Stadel

Abstract: Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U… ▽ More Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields. △ Less

Submitted 31 March, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 21 pages, 12 figures, accepted by MNRAS, comments welcome

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Authors: Erin Kado-Fong, Robyn E. Sanderson, Jenny E. Greene, Emily C. Cunningham, Coral Wheeler, T. K. Chan, Kareem El-Badry, Philip F. Hopkins, Andrew Wetzel, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Song Huang, Eliot Quataert, Tjitske Starkenburg

Abstract: Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosm… ▽ More Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies. △ Less

Submitted 2 May, 2022; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 28 pages, 15 figures. Accepted to ApJ

Sizing from the Smallest Scales: The Mass of the Milky Way

Authors: M. K. Rodriguez Wimberly, M. C. Cooper, D. C. Baxter, M. Boylan-Kolchin, J. S. Bullock, S. P. Fillingham, A. P. Ji, L. V. Sales, J. D. Simon

Abstract: As the Milky Way and its satellite system become more entrenched in near field cosmology efforts, the need for an accurate mass estimate of the Milky Way's dark matter halo is increasingly critical. With the second and early third data releases of stellar proper motions from {\it Gaia}, several groups calculated full $6$D phase-space information for the population of Milky Way satellite galaxies.… ▽ More As the Milky Way and its satellite system become more entrenched in near field cosmology efforts, the need for an accurate mass estimate of the Milky Way's dark matter halo is increasingly critical. With the second and early third data releases of stellar proper motions from {\it Gaia}, several groups calculated full $6$D phase-space information for the population of Milky Way satellite galaxies. Utilizing these data in comparison to subhalo properties drawn from the Phat ELVIS simulations, we constrain the Milky Way dark matter halo mass to be $\sim 1-1.2\times10^{12}~\msun$. We find that the kinematics of subhalos drawn from more- or less-massive hosts (i.e. $>1.2\times10^{12}~\msun$ or $<10^{12}~\msun$) are inconsistent, at the $3σ$ confidence level, with the observed velocities of the Milky Way satellites. The preferred host halo mass for the Milky Way is largely insensitive to the exclusion of systems associated with the Large Magellanic Cloud, changes in galaxy formation thresholds, and variations in observational completeness. As more Milky Way satellites are discovered, their velocities (radial, tangential, and total) plus Galactocentric distances will provide further insight into the mass of the Milky Way dark matter halo. △ Less

Submitted 10 May, 2022; v1 submitted 1 September, 2021; originally announced September 2021.

Comments: Key Figures: 4 and 5; 15 pages; 8 figures; accepted for publication in MNRAS 3 May 2022

Globular Clusters and Streaming Velocities: Testing the new formation channel in high-resolution cosmological simulations

Authors: Anna T. P. Schauer, Volker Bromm, Michael Boylan-Kolchin, Simon C. O. Glover, Ralf S. Klessen

Abstract: The formation of globular clusters and their relation to the distribution of dark matter have long puzzled astronomers. One of the most recently-proposed globular cluster formation channels ties ancient star clusters to the large-scale streaming velocity of baryons relative to dark matter in the early Universe. These streaming velocities affect the global infall of baryons into dark matter halos,… ▽ More The formation of globular clusters and their relation to the distribution of dark matter have long puzzled astronomers. One of the most recently-proposed globular cluster formation channels ties ancient star clusters to the large-scale streaming velocity of baryons relative to dark matter in the early Universe. These streaming velocities affect the global infall of baryons into dark matter halos, the high-redshift halo mass function, and the earliest generations of stars. In some cases, streaming velocities may result in dense regions of dark-matter-free gas that becomes Jeans unstable, potentially leading to the formation of compact star clusters. We investigate this hypothesis using cosmological hydrodynamical simulations that include a full chemical network and the formation and destruction of H$_2$, a process crucial for the formation of the first stars. We find that high-density gas in regions with significant streaming velocities -- which constitute approximately 1\% of the Universe -- is indeed somewhat offset from the centers of dark matter halos, but this offset is typically significantly smaller than the virial radius. Gas outside of dark matter halos never reaches Jeans-unstable densities in our simulations. We postulate that low-level ($Z \approx 10^{-3}\,Z_{\odot}$) metal enrichment by Population III supernovae may enable cooling in the extra-virial regions, allowing gas outside of dark matter halos to cool to the CMB temperature and become Jeans-unstable. Follow-up simulations that include both streaming velocities and metal enrichment by Population III supernovae are needed to understand if streaming velocities provide one path for the formation of globular clusters in the early Universe. △ Less

Submitted 18 August, 2021; originally announced August 2021.

Comments: 16 pages, 10 figures, submitted to ApJ

The effects of LMC-mass environments on their dwarf satellite galaxies in the FIRE simulations

Authors: Ethan D. Jahn, Laura V. Sales, Andrew Wetzel, Jenna Samuel, Kareem El-Badry, Michael Boylan-Kolchin, James S. Bullock

Abstract: Characterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass halos ($M_{200m} \sim 10^{11}$ M$_\odot$) on their populations of satellites ($M_\star \geq 10^4$ M$_\odot$) using a… ▽ More Characterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass halos ($M_{200m} \sim 10^{11}$ M$_\odot$) on their populations of satellites ($M_\star \geq 10^4$ M$_\odot$) using a suite of zoom-in simulations from the Feedback In Realistic Environments (FIRE) project. Our simulations predict significant hot coronas with $T\sim10^5$ K and $M_\text{gas}\sim10^{9.5}$ M$_\odot$. We identify signatures of environmental quenching in dwarf satellite galaxies, particularly for satellites with intermediate mass ($M_\star = 10^{6-7}$ M$_\odot$). The gas content of such objects indicates ram-pressure as the likely quenching mechanism, sometimes aided by star formation feedback. Satellites of LMC-mass hosts replicate the stellar mass dependence of the quiescent fraction found in satellites of MW mass hosts (i.e. that the quiescent fraction increases as stellar mass decreases). Satellites of LMC-mass hosts have a wider variety of quenching times when compared to the strongly bi-modal distribution of quenching times of nearby centrals. Finally, we identify significant tidal stellar structures around four of our six LMC-analogs, suggesting that stellar streams may be common. These tidal features originated from satellites on close orbits, extend to $\sim$80 kpc from the central galaxy, and contain $\sim10^{6-7}$ M$_\odot$~of stars. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 14 pages, 10 figures, submitted to MNRAS

HETDEX [OIII] Emitters I: A spectroscopically selected low-redshift population of low-mass, low-metallicity galaxies

Authors: Briana Indahl, Greg Zeimann, Gary J. Hill, William P. Bowman, Robin Ciardullo, Niv Drory, Eric Gawiser, Ulrich Hopp, Steven Janowiecki, Michael Boylan-Kolchin, Erin Mentuch Cooper, Dustin Davis, Daniel Farrow, Steven Finkelstein, Caryl Gronwall, Andreas Kelz, Kristen B. W. McQuinn, Don Schneider, Sarah E. Tuttle

Abstract: We assemble a sample of 17 low metallicity (7.45 < log(O/H)+12 < 8.12) galaxies with z < 0.1 found spectroscopically, without photometric pre-selection, in early data from the Hobby Eberly Telescope Dark Energy Experiment (HETDEX). Star forming galaxies that occupy the lowest mass and metallicity end of the mass-metallicity relation tend to be under sampled in continuum-based surveys as their spec… ▽ More We assemble a sample of 17 low metallicity (7.45 < log(O/H)+12 < 8.12) galaxies with z < 0.1 found spectroscopically, without photometric pre-selection, in early data from the Hobby Eberly Telescope Dark Energy Experiment (HETDEX). Star forming galaxies that occupy the lowest mass and metallicity end of the mass-metallicity relation tend to be under sampled in continuum-based surveys as their spectra are typically dominated by emission from newly forming stars. We search for galaxies with high [OIII]$λ$5007 / [OII]$λ$3727, implying highly ionized nebular emission often indicative of low metallicity systems. With the Second Generation Low Resolution Spectrograph on the Hobby Eberly Telescope we acquired follow-up spectra, with higher resolution and broader wavelength coverage, of each low-metallicity candidate in order to confirm the redshift, measure the H$α$ and [NII] line strengths and, in many cases, obtain deeper spectra of the blue lines. We find our galaxies are consistent with the mass-metallicity relation of typical low mass galaxies. However, galaxies in our sample tend to have similar specific star formation rates (sSFRs) as the incredibly rare "blueberry" galaxies found in (Yang et. al. 2017). We illustrate the power of spectroscopic surveys for finding low mass and metallicity galaxies and reveal that we find a sample of galaxies that are a hybrid between the properties of typical dwarf galaxies and the more extreme blueberry galaxies. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Comments: 32 pages, 12 figures, 8 tables, Accepted for publication in ApJ