Showing 1–50 of 78 results for author: Brooks, A M

The Simulated Oxygen Shortage (SOS): Mapping the Missing OVI in Simulated Dwarf Galaxies to Subgrid Physics

Authors: Daniel R. Piacitelli, Alyson M. Brooks, N. Nicole Sanchez, Hetvi Khatri, Charlotte Christensen, Cameron Hummels, Nishant Mishra, Akaxia Cruz, Ben Keller, Thomas R. Quinn, Sijing Shen, James Wadsley

Abstract: Observations reveal extended \OVI\, reservoirs in the circumgalactic medium (CGM) of dwarf galaxies, yet current simulations systematically underpredict \OVI\, column densities. Utilizing two suites run with different simulation codes, the \MM\, simulations (Marvelous Massive Dwarfs and Marvel-ous Dwarfs) and the publicly available FIRE-2 simulations, we explore the role of subgrid models and the… ▽ More Observations reveal extended \OVI\, reservoirs in the circumgalactic medium (CGM) of dwarf galaxies, yet current simulations systematically underpredict \OVI\, column densities. Utilizing two suites run with different simulation codes, the \MM\, simulations (Marvelous Massive Dwarfs and Marvel-ous Dwarfs) and the publicly available FIRE-2 simulations, we explore the role of subgrid models and the resulting CGM phase in shaping \OVI\, production. By comparing observationally derived \OVI\, masses to the mass of oxygen produced over the galaxies' star formation history, we find evidence for an underproduction of oxygen for low-mass simulated galaxies. Despite clear differences in feedback implementation, CGM structure, and metal mixing, we find that \OVI\, in both suites generally self-selects cool/warm ($\rm log\, T\,/K \sim 4.5$), diffuse ($\rm log\,n_{gas}\,/cm^{-3} \sim -5.0 $), and moderately metal-enriched ($\rm log\, Z/Z_{\odot} \sim -1 $) material at large radii from the galaxy. We show that neither the choice of ultraviolet background nor plausible variations in CGM thermal structure can close the gap with observations. Taken together, our results point to a possible underproduction of oxygen in low-mass galaxies. Feedback prescriptions contribute via insufficient metal transport to large radii. Hence, the \OVI\, deficit may motivate an investigation of current modeling choices for supernova yields, star formation, and feedback in low-metallicity environments. △ Less

Submitted 14 May, 2026; originally announced May 2026.

Comments: 28 pages, 16 Figures

$V/σ$ Trends with Mass for Dwarf Galaxies from the Marvelous Massive Dwarfs Suite

Authors: Dilys Ruan, Alyson M. Brooks, Leonardo A. Barba, Mithi A. C. de los Reyes, Akaxia Cruz, Robel Geda, Annika H. G. Peter, Benjamin W. Keller, Thomas Quinn, James W. Wadsley

Abstract: Galaxy formation scenarios can be interpreted through galaxy morphology and the level of rotational versus pressure support, quantified through the ratio of a galaxy's rotation speed to its velocity dispersion: $V/σ$. Observational studies of dwarf galaxies find that $V/σ$ does not strongly depend on environment, and may weakly depend on galaxy mass, which could shift our understanding of how dwar… ▽ More Galaxy formation scenarios can be interpreted through galaxy morphology and the level of rotational versus pressure support, quantified through the ratio of a galaxy's rotation speed to its velocity dispersion: $V/σ$. Observational studies of dwarf galaxies find that $V/σ$ does not strongly depend on environment, and may weakly depend on galaxy mass, which could shift our understanding of how dwarf galaxies form. We utilize the Marvelous Massive Dwarfs suite to examine whether $V/σ$ depends on mass in simulations, and understand how this varies for different baryonic components of the galaxy: HI gas, young stars ($<$ 1 Gyr) and old stars ($>$ 1 Gyr). We use a simulation sample of 67 isolated dwarf galaxies with M$_\star=10^6-10^9$ M$_\odot$ and produce line-of-sight maps for rotation speed and dispersion for different viewing angles of each galaxy. We find that $V/σ$ increases with mass, and that HI gas and young stars are more rotation-supported ($V/σ\approx 1-13$) while old stars are more dispersion-supported ($V/σ\approx 0.2-5$). This result is consistent with the scenario where young stars are born from dynamically cold gas in the interstellar medium and undergo dynamical heating over time. We quantify the effects of spatial resolution in observational determinations of $V/σ$ and find that existing observations using old stars may underestimate the intrinsic $V/σ$. We find a correlation between $V/σ_\mathrm{HI,global}$ and HI line profile shape that is qualitatively similar to previous simulation results, but we find higher $V/σ_\mathrm{HI,global}$ compared to prior work which found values $\lesssim 2$ for most galaxies in this mass range. Our results motivate future work to examine $V/σ$ and dwarf galaxy formation with different kinematic tracers of the galaxy. △ Less

Submitted 7 May, 2026; originally announced May 2026.

Comments: 26 pages, 12 figures, submitted to ApJ

MARVELously Dark: the density profile evolution of dwarf halos in velocity-dependent SIDM

Authors: Anna Engelhardt, Ferah Munshi, Annika H. G. Peter, Ethan O. Nadler, Akaxia Cruz, Alyson M. Brooks, Zhichao Carton Zeng, Thomas R. Quinn, Blake Keith

Abstract: Self-interacting dark matter (SIDM) with a sufficiently large cross section has been shown to naturally produce constant dark matter (DM) cores, as well as core-collapse, at the centers of dwarf halos on cosmic timescales, potentially reducing tensions with observation. Here, we present halos from a new dark matter only (DMO) cosmological (SIDM) simulation: Ms.Marvel DMO with a velocity-dependent… ▽ More Self-interacting dark matter (SIDM) with a sufficiently large cross section has been shown to naturally produce constant dark matter (DM) cores, as well as core-collapse, at the centers of dwarf halos on cosmic timescales, potentially reducing tensions with observation. Here, we present halos from a new dark matter only (DMO) cosmological (SIDM) simulation: Ms.Marvel DMO with a velocity-dependent self-interaction cross section with $σ/m_\text{max} = 50$ cm$^2$/g at $v_\text{max} = 35$ km/s. We compare these to the CDM suite of Storm simulations including both DMO and dark matter + hydrodynamics runs, in order to test core-formation (and core-collapse) across different dark matter models. We show that Ms.Marvel DMO can reproduce core slopes consistent with observations of isolated dwarf galaxies and more massive ($\text{M}_{vir} \gtrsim 10^{10} M_{\odot}$) CDM dwarf halos that include stellar feedback from the matched CDM run (Storm CDM+baryons). We identify nine Ms.Marvel SIDM DMO halos in the core-collapse phase of gravothermal evolution with halo masses below $2\times 10^9 M_{\odot}$. We find that using core slope to measure the core-collapse timescales of Ms.Marvel DMO halos agrees well with predicted collapse times estimated with the parametric model for SIDM halos introduced by Yang et al.(2023). Additionally, compared to central density, core slope is less sensitive to both the radius of measurement and halo merger history. These results indicate that the slope of the inner DM density profile more cleanly differentiates core-collapsed versus core-forming halos than central density amplitude. △ Less

Submitted 11 June, 2026; v1 submitted 30 January, 2026; originally announced January 2026.

Comments: 47 pages, 15 figures

CASCO: Cosmological and AStrophysical parameters from Cosmological simulations and Observations IV. Testing warm dark matter cosmologies with galaxy scaling relations: A joint simulation-observation study using DREAMS simulations

Authors: M. Silvestrini, C. Tortora, V. Busillo, Alyson M. Brooks, A. Farahi, A. M. Garcia, N. Kallivayalil, N. R. Napolitano, J. C. Rose, P. Torrey, F. Villaescusa-Navarro, M. Vogelsberger

Abstract: Small-scale discrepancies in the standard Lambda cold dark matter paradigm have motivated the exploration of alternative dark matter (DM) models, such as warm dark matter (WDM). We investigate the constraining power of galaxy scaling relations on cosmological, astrophysical, and WDM parameters through a joint analysis of hydrodynamic simulations and observational data. Our study is based on the DR… ▽ More Small-scale discrepancies in the standard Lambda cold dark matter paradigm have motivated the exploration of alternative dark matter (DM) models, such as warm dark matter (WDM). We investigate the constraining power of galaxy scaling relations on cosmological, astrophysical, and WDM parameters through a joint analysis of hydrodynamic simulations and observational data. Our study is based on the DREAMS project and combines large-volume uniform-box simulations with high-resolution Milky Way zoom-in runs in a $Λ$WDM cosmology. To ensure consistency between the different simulation sets, we apply calibrations to account for resolution effects, allowing us to exploit the complementary strengths of the two suites. We compare simulated relations, including stellar size, DM mass and fraction within the stellar half-mass radius, and the total-to-stellar mass ratio, with two complementary galaxy samples: the SPARC catalog of nearby spirals and the LVDB catalog of dwarf galaxies in the Local Volume. Using a bootstrap-based fitting procedure, we show that key cosmological parameters ($Ω_m$, $σ_8$) and supernova feedback strength can be recovered with good accuracy, particularly from the uniform-box simulations. While the WDM particle mass remains unconstrained, the zoom-in simulations reveal subtle WDM-induced trends at low stellar masses in both the DM mass and total-to-stellar mass ratio. We also find that the galaxy stellar mass function exhibits a measurable dependence on the WDM particle mass below log10(M_*/Msun) <~ 8, which appears separable from the impact of feedback, suggesting it as a promising complementary probe. Our results highlight the importance of combining multi-resolution simulations with diverse observational datasets to jointly constrain baryonic processes and DM properties. △ Less

Submitted 12 January, 2026; originally announced January 2026.

Journal ref: A&A 706, A382 (2026)

The DREAMS Project: Disentangling the Impact of Halo-to-Halo Variance and Baryonic Feedback on Milky Way Dark Matter Speed Distributions

Authors: Ethan Lilie, Jonah C. Rose, Mariangela Lisanti, Alex M. Garcia, Paul Torrey, Kassidy E. Kollmann, Jiaxuan Li, Olivia Mostow, Bonny Y. Wang, Stephanie O'Neil, Xuejian Shen, Alyson M. Brooks, Arya Farahi, Nitya Kallivayalil, Lina Necib, Andrew B. Pace, Mark Vogelsberger

Abstract: Direct detection experiments require information about the local dark matter speed distribution to produce constraints on dark matter candidates, or infer their properties in the event of a discovery. In this paper, we analyze how the uncertainty in the dark matter speed distribution near the Sun is affected by baryonic feedback, halo-to-halo variance, and halo mass. To do so, we harness the stati… ▽ More Direct detection experiments require information about the local dark matter speed distribution to produce constraints on dark matter candidates, or infer their properties in the event of a discovery. In this paper, we analyze how the uncertainty in the dark matter speed distribution near the Sun is affected by baryonic feedback, halo-to-halo variance, and halo mass. To do so, we harness the statistical power of the new DREAMS Cold Dark Matter simulation suite, which is comprised of 1024 zoom-in Milky Way-mass halos with varied initial conditions as well as cosmological and astrophysical parameters. Applying a normalizing flows emulator to these simulations, we find that the uncertainty in the local DM speed distribution is dominated by halo-to-halo variance and, to a lesser extent, uncertainty in host halo mass. Uncertainties in supernova and black hole feedback (from the IllustrisTNG model in this case) are negligible in comparison. Using the DREAMS suite, we present a state-of-the-art prediction for the DM speed distribution in the Milky Way. Although the Standard Halo Model is contained within the uncertainty of this prediction, individual galaxies may have distributions that differ from it. Lastly, we apply our DREAMS results to the XENON1T experiment and demonstrate that the astrophysical uncertainties are comparable to the experimental ones, solidifying previous results in the literature obtained with a smaller sample of simulated Milky Way-mass halos. △ Less

Submitted 3 December, 2025; originally announced December 2025.

Comments: 26 pages, 7 figures, 3 appendices, 7 additional figures. See also arxiv.2512.00148 and arXiv:2512.02095

The DREAMS Project: Disentangling the Impact of Halo-to-Halo Variance and Baryonic Feedback on Milky Way Satellite Galaxies

Authors: Jonah C. Rose, Mariangela Lisanti, Paul Torrey, Francisco Villaescusa-Navarro, Alex M. Garcia, Arya Farahi, Carrie Filion, Alyson M. Brooks, Nitya Kallivayalil, Kassidy E. Kollmann, Ethan Lilie, Jiaxuan Li, Olivia Mostow, Akaxia Cruz, Tri Nguyen, Sandip Roy, Andrew B. Pace, Niusha Ahvazi, Stephanie O'Neil, Xuejian Shen, Francis-Yan Cyr-Racine, Adrian M. Price-Whelan, Marla Geha, Lina Necib, Mark Vogelsberger , et al. (2 additional authors not shown)

Abstract: We analyze the properties of satellite galaxies around 1,024 Milky Way-mass hosts from the DREAMS Project, simulated within a $Λ$CDM cosmology. Utilizing the TNG galaxy-formation model, the DREAMS simulations incorporate both baryonic physics and cosmological uncertainties for a large sample of galaxies with diverse environments and formation histories. We investigate the relative impact of the ph… ▽ More We analyze the properties of satellite galaxies around 1,024 Milky Way-mass hosts from the DREAMS Project, simulated within a $Λ$CDM cosmology. Utilizing the TNG galaxy-formation model, the DREAMS simulations incorporate both baryonic physics and cosmological uncertainties for a large sample of galaxies with diverse environments and formation histories. We investigate the relative impact of the physical uncertainty from the galaxy-formation model on predicted satellite properties using four metrics: the satellite stellar mass function, radial distribution, inner slope of dark matter density profile, and stellar half-light radius. We compare these predictions to observations from the SAGA Survey and the DREAMS N-body simulations and find that uncertainties from baryonic physics modeling are subdominant to the scatter arising from halo-to-halo variance. Where baryonic modeling does affect satellites, the supernova wind energy has the largest effect on the satellite properties that we investigate. Specifically, increased supernova wind energy suppresses the stellar mass of satellites and results in more extended stellar half-light radii. The adopted wind speed has only a minor impact, and other astrophysical and cosmological parameters show no measurable effect. Our findings highlight the robustness of satellite properties against uncertainties in baryonic physics modeling. △ Less

Submitted 1 December, 2025; originally announced December 2025.

Comments: 16 pages, 5 figures, 2 tables, and 2 appendices with an additional 4 figures

The DREAMS Project: A New Suite of 1,024 Simulations to Contextualize the Milky Way and Assess Physics Uncertainties

Authors: Jonah C. Rose, Mariangela Lisanti, Paul Torrey, Francisco Villaescusa-Navarro, Alex M. Garcia, Arya Farahi, Carrie Filion, Alyson M. Brooks, Nitya Kallivayalil, Kassidy E. Kollmann, Ethan Lilie, Bonny Y. Wang, Akaxia Cruz, Sandip Roy, Andrew B. Pace, Niusha Ahvazi, Stephanie O'Neil, Cian Roche, Xuejian Shen, Mark Vogelsberger

Abstract: We introduce a new suite of 1,024 cosmological and hydrodynamical zoom-in simulations of Milky Way-mass halos, run with Cold Dark Matter, as part of the DREAMS Project. Each simulation in the suite has a unique set of initial conditions and combination of cosmological and astrophysical parameters. The suite is designed to quantify theoretical uncertainties from halo-to-halo variance, as well as st… ▽ More We introduce a new suite of 1,024 cosmological and hydrodynamical zoom-in simulations of Milky Way-mass halos, run with Cold Dark Matter, as part of the DREAMS Project. Each simulation in the suite has a unique set of initial conditions and combination of cosmological and astrophysical parameters. The suite is designed to quantify theoretical uncertainties from halo-to-halo variance, as well as stellar and black hole feedback. We develop a novel weighting scheme that prioritizes regions of the input parameter space, yielding galaxies consistent with the observed present-day stellar mass--halo mass relation. The resulting galaxy population exhibits a wide diversity in structural properties that encompasses those of the actual Milky Way, providing a powerful statistical sample for galactic archaeology. To demonstrate the suite's scientific utility, we investigate the connection between a galaxy's merger history, focusing on Gaia-Sausage-Enceladus~(GSE) analogs, and its present-day properties. We find that galaxies with a GSE analog have lower star formation rates, more compact disks, and more spherical stellar halos. Crucially, significant halo-to-halo scatter remains, demonstrating that matching more than the most significant events in the Milky Way's past is necessary to recover its present-day properties. Our results highlight the necessity for large statistical samples to disentangle the stochastic nature of galaxy formation and robustly model the Milky Way's unique history. △ Less

Submitted 28 November, 2025; originally announced December 2025.

Comments: 19 pages, 5 figures, 2 tables, 2 appendices with an additional 11 figures

The Merger-Driven Formation of Classical Low Surface Brightness Galaxies in Romulus25

Authors: Anna C. Wright, Alyson M. Brooks, Michael Tremmel, Jason E. Young, Ferah Munshi, Tom Quinn

Abstract: We use the Romulus25 cosmological simulation volume to study a large sample of late-type gas-rich galaxies with low central surface brightnesses known as classical low surface brightness (LSB) galaxies and compare them to a mass-matched sample of high surface brightness (HSB) galaxies. We find that classical LSB galaxies make up a substantial fraction of the galaxy population, accounting for ~60%… ▽ More We use the Romulus25 cosmological simulation volume to study a large sample of late-type gas-rich galaxies with low central surface brightnesses known as classical low surface brightness (LSB) galaxies and compare them to a mass-matched sample of high surface brightness (HSB) galaxies. We find that classical LSB galaxies make up a substantial fraction of the galaxy population, accounting for ~60% of all central galaxies with 8$\leq$log$_\mathrm{10}$(M$_\star$/M$_\odot$)$\leq$10. In Romulus25, classical LSB galaxies are predominantly formed through major mergers in which the secondary galaxy is co-rotating and aligned with the primary galaxy's gas disk and/or has above average orbital angular momentum at infall. The merger product is a high spin galaxy in which star formation is spread out and inefficient, allowing the galaxy to build up a large supply of relatively unenriched gas. The star formation rates of LSB galaxies are nearly constant over time, leading to stellar populations that are, on average, slightly older and therefore optically redder than those of similar HSB galaxies. However, because LSB galaxies are diffuse and metal-poor, they have very little internal reddening, causing them to appear bluer than HSB galaxies. We also find that, when compared to the bulges of HSB galaxies, the bulges of LSB galaxies are similar in mass, but are lower surface brightness, redder, and more diffuse on average. Despite these differences, classical LSB galaxies are part of the continuum of the galaxy population in Romulus25, constituting one of many evolutionary paths. △ Less

Submitted 23 July, 2026; v1 submitted 28 July, 2025; originally announced July 2025.

Comments: main text is 27 pages, 17 figures; accepted to ApJ

Survivors and Zombies: The Quenching and Disruption of Satellites around Milky Way Analogs

Authors: Debosmita Pathak, Charlotte R. Christensen, Alyson M. Brooks, Ferah Munshi, Anna C. Wright, Courtney Carter

Abstract: It is necessary to understand the full accretion history of the Milky Way in order to contextualize the properties of observed Milky Way satellite galaxies and the stellar halo. This paper compares the dynamical properties and star-formation histories of surviving and disrupted satellites around Milky Way-like galaxies using the DC Justice League suite of very high-resolution cosmological zoom-in… ▽ More It is necessary to understand the full accretion history of the Milky Way in order to contextualize the properties of observed Milky Way satellite galaxies and the stellar halo. This paper compares the dynamical properties and star-formation histories of surviving and disrupted satellites around Milky Way-like galaxies using the DC Justice League suite of very high-resolution cosmological zoom-in simulations of Milky Way analogs and their halo environments. We analyze the full census of galaxies accreted within the past 12 Gyrs, which including both surviving satellites at $z=0$, and dwarf galaxies that disrupted and merged with the host prior to $z=0$. Our simulations successfully reproduce the trends in $M_*$-[Fe/H]-[$α$/Fe] observed in surviving Milky Way satellites and disrupted stellar streams, indicating earlier star-formation for disrupted progenitors. We find the likelihood and timescales for quenching and disruption are strongly correlated with the mass and time of infall. In particular, none of the galaxies accreted more than 12 Gyrs ago survived, and only 20% of all accreted galaxies with $M_*>10^8M_\odot$ survive. Additionally, satellites with highly radial trajectories are more likely to quench and disrupt. Disruption proceeds quickly for $\geq10^6M_\odot$ satellites accreted $10{-}12$ Gyr ago, often on timescales similar to the $\sim300$ Myr snapshot spacing. For high-mass satellites, the disruption timescale is faster than the quenching timescale. As a result, 92% of disrupted galaxies remain star-forming up until disruption. In contrast, Ultra Faint Dwarfs (UFDs) tend to quench prior to accretion, and 94% of UFDs accreted up to 12 Gyr ago survive at $z=0$. △ Less

Submitted 13 August, 2025; v1 submitted 28 May, 2025; originally announced May 2025.

Comments: 18 pages, 7 figures, 2 tables, published in ApJ

Marvelous Metals: Surveying the Circumgalactic Medium of Simulated Dwarf Galaxies

Authors: Daniel R. Piacitelli, Alyson M. Brooks, Charlotte Christensen, N. Nicole Sanchez, Yakov Faerman, Sijing Shen, Akaxia Cruz, Ben Keller, Thomas R. Quinn, James Wadsley

Abstract: Dwarf galaxies are uniquely sensitive to energetic feedback processes and are known to experience substantial mass and metal loss from their disk. Here, we investigate the circumgalactic medium (CGM) of 64 isolated dwarf galaxies ($6.0<$log(M$_*/M_{\odot}$)$<9.5$) at $z=0$ from the Marvel-ous Dwarfs and Marvelous Massive Dwarfs hydrodynamic simulations. Our galaxies produce column densities broadl… ▽ More Dwarf galaxies are uniquely sensitive to energetic feedback processes and are known to experience substantial mass and metal loss from their disk. Here, we investigate the circumgalactic medium (CGM) of 64 isolated dwarf galaxies ($6.0<$log(M$_*/M_{\odot}$)$<9.5$) at $z=0$ from the Marvel-ous Dwarfs and Marvelous Massive Dwarfs hydrodynamic simulations. Our galaxies produce column densities broadly consistent with current observations. We investigate these column densities in the context of mass and metal retention rates and the physical properties of the CGM. We find $48\pm11\%$ of all baryons within $R_{200c}$ reside in the CGM, with $\sim70\%$ of CGM mass existing in a warm gas phase, $10^{4.5}<T<10^{5.5}$ K that dominates beyond $r/R_{200c}\sim0.5$. Further, the warm and cool ($10^{4.0}<T<10^{4.5}$ K) gas phases each retain $5-10\%$ of metals formed by the dwarf galaxy. The significant fraction of mass and metals residing in the warm CGM phase provides an interpretation for the lack of observed low ion detections beyond $b/R_{200c}\sim0.5$ at $z\sim0$, as the majority of mass in this region exists in higher ions. We find a weak correlation between galaxy mass and total CGM metal retention despite the fraction of metals lost from the halo increasing from $\sim10\%$ to $>40\%$ towards lower masses. Our findings highlight the CGM (primarily its warm component) as a key reservoir of mass and metals for dwarf galaxies across stellar masses and underscore its importance in understanding the baryon cycle in the low-mass regime. Finally, we provide individual galaxy properties of our full sample and quantify the fraction of ultraviolet observable mass to support future observational programs, particularly those aimed at performing a metal budget around dwarf galaxies. △ Less

Submitted 3 December, 2025; v1 submitted 13 May, 2025; originally announced May 2025.

Comments: 31 pages, 16 figures

Predictions for Detecting a Turndown in the Baryonic Tully Fisher Relation

Authors: Dilys Ruan, Alyson M. Brooks, Akaxia Cruz, Annika H. G. Peter, Benjamin W. Keller, Thomas Quinn, James Wadsley, Elizabeth A. K. Adams

Abstract: The baryonic Tully Fisher relation (bTFR) provides an empirical connection between baryonic mass and dynamical mass (measured by the maximum rotation velocity) for galaxies. Due to the impact of baryonic feedback in the shallower potential wells of dwarf galaxies, the bTFR is predicted to turn down at low masses from the extrapolated power-law relation at high masses. The low-mass end of the bTFR… ▽ More The baryonic Tully Fisher relation (bTFR) provides an empirical connection between baryonic mass and dynamical mass (measured by the maximum rotation velocity) for galaxies. Due to the impact of baryonic feedback in the shallower potential wells of dwarf galaxies, the bTFR is predicted to turn down at low masses from the extrapolated power-law relation at high masses. The low-mass end of the bTFR is poorly constrained due to small samples and difficulty in connecting the galaxy's gas kinematics to its dark matter halo. Simulations can help us understand this connection and interpret observations. We measure the bTFR with 75 dwarf galaxies from the Marvel-ous and Marvelous Massive Dwarfs hydrodynamic simulations. Our sample has M$_\star = 10^6-10^9$ M$_\odot$, and is mostly gas dominated. We compare five velocity methods: V$_\mathrm{out,circ}$ (spatially resolved mass-enclosed), V$_\mathrm{out,mid}$ (spatially resolved midplane gravitational potential), and unresolved HI linewidths at different percentages of the peak flux (W$_\mathrm{10}$, W$_\mathrm{20}$, and W$_\mathrm{50}$). We find an intrinsic turndown in the bTFR for maximum halo speeds $\lesssim 50$ km s$^{-1}$, or total baryonic mass M$_\mathrm{bary}\lesssim 10^{8.5}$ M$_\odot$. We find that observing HI in lower-mass galaxies to the conventional surface density limit of 1 M$_\odot$ pc$^{-2}$ is not enough to detect a turndown in the bTFR; none of the HI velocity methods, spatially resolved or unresolved, recover the turndown, and we find bTFR slopes consistent with observations of higher-mass galaxies. However, we predict that the turndown can be recovered by resolved rotation curves if the HI limit is $\lesssim 0.08$ M$_\odot$ pc$^{-2}$, which is within the sensitivity of current HI surveys like FEASTS and MHONGOOSE. △ Less

Submitted 18 July, 2025; v1 submitted 20 March, 2025; originally announced March 2025.

Comments: 16 pages, 11 figures, Accepted to MNRAS

Journal ref: Mon Not R Astron Soc (2025) 2180-2196

The JWST Resolved Stellar Populations Early Release Science Program. VIII. The Spatially Resolved Star Formation History of WLM

Authors: Roger E. Cohen, Kristen B. W. McQuinn, Alessandro Savino, Max J. B. Newman, Daniel R. Weisz, Andrew E. Dolphin, Martha L. Boyer, Matteo Correnti, Marla C. Geha, Mario Gennaro, Karoline M. Gilbert, Nitya Kallivayalil, Jack T. Warfield, Benjamin F. Williams, Alyson M. Brooks, Andrew A. Cole, Evan D. Skillman, Christopher T. Garling, Jason S. Kalirai, Jay Anderson

Abstract: We measure radial stellar age gradients in the relatively isolated gas-rich dwarf irregular WLM, combining JWST NIRCam and NIRISS imaging with six archival Hubble fields over semi-major axis equivalent distances of 0$\lesssim$R$_{SMA}$$\lesssim$4 kpc ($\lesssim$3R$_{hl}$). Fitting lifetime star formation histories (SFHs) to resolved color-magnitude diagrams (CMDs), radial age gradients are quantif… ▽ More We measure radial stellar age gradients in the relatively isolated gas-rich dwarf irregular WLM, combining JWST NIRCam and NIRISS imaging with six archival Hubble fields over semi-major axis equivalent distances of 0$\lesssim$R$_{SMA}$$\lesssim$4 kpc ($\lesssim$3R$_{hl}$). Fitting lifetime star formation histories (SFHs) to resolved color-magnitude diagrams (CMDs), radial age gradients are quantified using $τ_{90}$ and $τ_{50}$, the lookback times to form 90\% and 50\% of the cumulative stellar mass. We find that globally, the outskirts of WLM are older on average, with ($δ$$τ_{90}$, $δ$$τ_{50}$)/$δ$R$_{SMA}=$(0.82$^{+0.10}_{-0.10}$, 1.60$^{+0.23}_{-0.22}$) Gyr/kpc (stat.), in good agreement with simulations. However, we also detect an azimuthal dependence of radial stellar age gradients, finding that stars on the leading edge of WLM (relative to its proper motion) are both younger and have a flatter age gradient compared to the trailing edge. This difference persists over 0.6$\lesssim$R$_{SMA}$$\lesssim$3.2 kpc ($\sim$0.5$-$2.5R$_{hl}$) and lookback times up to $\sim$8 Gyr, and is robust to assumed stellar evolutionary model. Our results are consistent with star formation triggered by ram pressure stripping from a circumgalactic and/or intergalactic medium, suggested by recent HI observations. If confirmed, processes typifying dense environments, such as ram pressure stripping, may be more relevant to the evolution of isolated galaxies than previously thought. △ Less

Submitted 19 February, 2025; originally announced February 2025.

Comments: ApJ in press. 23 pages, 8 figures, 3 tables

A MARVEL-ous study of how well galaxy shapes reflect Dark Matter halo shapes in Cold Dark Matter Simulations

Authors: Blake Keith, Ferah Munshi, Alyson M. Brooks, Jordan Van Nest, Anna Engelhardt, Akaxia Cruz, Ben Keller, Thomas Quinn, James Wadsley

Abstract: We present a 3D shape analysis of both dark matter (DM) and stellar matter (SM) in simulated dwarf galaxies to determine whether stellar shape traces DM shape. Using 80 central and satellite galaxies from three simulation suites (Marvelous Massive Dwarfs, Marvelous Dwarfs, and DC Justice League) spanning stellar masses of $10^6$--$10^{10}$ $M_\odot$, we measure 3D shapes through the moment of iner… ▽ More We present a 3D shape analysis of both dark matter (DM) and stellar matter (SM) in simulated dwarf galaxies to determine whether stellar shape traces DM shape. Using 80 central and satellite galaxies from three simulation suites (Marvelous Massive Dwarfs, Marvelous Dwarfs, and DC Justice League) spanning stellar masses of $10^6$--$10^{10}$ $M_\odot$, we measure 3D shapes through the moment of inertia tensor at two times the effective radius to derive axis ratios ($C/A$, $B/A$) and triaxiality. We find that stellar shape does indeed follow DM halo shape for our dwarf galaxies. However, the presence of a stellar disk in more massive dwarfs ($M_* \gtrsim 10^{7.5}$ $M_\odot$) pulls the distribution of stellar $C/A$ ratios to lower values, while in lower mass galaxies the gravitational potential remains predominantly shaped by DM. Similarly, stellar triaxiality generally tracks dark matter halo triaxiality, with this relationship being particularly strong for non-disky galaxies though weaker in disky systems. This correlation is reinforced by strong alignment between SM and DM axes, particularly in disk galaxies. Further, we find no detectable difference in either SM or DM shape comparing two different SNe feedback implementations, demonstrating that shape measurements may be robust to different implementations of baryonic feedback in dwarf galaxies. We also observe that a dwarf galaxy's shape is largely unperturbed by recent mergers (with merger ratios $>4$). This comprehensive study demonstrates that stellar shape measurements can serve as a reliable tool for inferring DM shapes in dwarf galaxies. △ Less

Submitted 12 June, 2025; v1 submitted 27 January, 2025; originally announced January 2025.

Comments: 18 pages, 16 figures, Published in ApJ

Journal ref: Journal-ref: Astrophysical Journal 986 (2025) 138

Testable predictions of outside-in age gradients in dwarf galaxies of all types

Authors: Claire L. Riggs, Alyson M. Brooks, Ferah Munshi, Charlotte R. Christensen, Roger E. Cohen, Thomas R. Quinn, James Wadsley

Abstract: We use a sample of 73 simulated satellite and central dwarf galaxies spanning a stellar mass range of $10^{5.3}-10^{9.1} M_\odot$ to investigate the origin of their stellar age gradients. We find that dwarf galaxies often form their stars "inside-out," i.e., the stars form at successively larger radii over time. However, the oldest stars get reshuffled beyond the star forming radius by fluctuation… ▽ More We use a sample of 73 simulated satellite and central dwarf galaxies spanning a stellar mass range of $10^{5.3}-10^{9.1} M_\odot$ to investigate the origin of their stellar age gradients. We find that dwarf galaxies often form their stars "inside-out," i.e., the stars form at successively larger radii over time. However, the oldest stars get reshuffled beyond the star forming radius by fluctuations in the gravitational potential well caused by stellar feedback (the same mechanisms that cause dwarfs to form dark matter cores). The result is that many dwarfs appear to have an "outside-in" age gradient at $z=0$, with younger stellar populations more centrally concentrated. However, for the reshuffled galaxies with the most extended star formation, young stars can form out to the large radii to which the old stars have been reshuffled, erasing the age gradient. We find that major mergers do not play a significant role in setting the age gradients of dwarfs. We find similar age gradient trends in satellites and field dwarfs, suggesting environment plays only a minor role, if any. Finally, we find that the age gradient trends are imprinted on the galaxies at later times, suggesting that the stellar reshuffling dominates after the galaxies have formed 50% of their stellar mass. The later reshuffling is at odds with results from the FIRE-2 simulations. Hence, age gradients offer a test of current star formation and feedback models that can be probed via observations of resolved stellar populations. △ Less

Submitted 10 December, 2024; v1 submitted 19 August, 2024; originally announced August 2024.

Comments: 25 pages, 9 figures. Accepted to ApJ

Journal ref: ApJ. 977, 1 (2024) 20

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

Bursting with Feedback: The Relationship between Feedback Model and Bursty Star Formation Histories in Dwarf Galaxies

Authors: Bianca Azartash-Namin, Anna Engelhardt, Ferah Munshi, B. W. Keller, Alyson M. Brooks, Jordan Van Nest, Charlotte R. Christensen, Tom Quinn, James Wadsley

Abstract: We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capa… ▽ More We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capable of producing galaxies that are cored and reproduce observed scaling relations for metallicity, luminosity, mass, and size. We show that continuous bursty star formation and the resulting stellar feedback are able to sustain dark matter cores in the higher dwarf galaxy mass regime, while the majority of ultra-faint and classical dwarfs retain cuspy central dark matter density profiles. We find that both subgrid SN models are able to create bursty star formation histories. We find that effective core formation peaks at M$_*$/M$_\odot$ $\simeq 5 \times 10^{-3}$ for both feedback models. Galaxies simulated with superbubble feedback peak at lower mean burstiness values relative to blastwave feedback, indicating that core formation in the superbubble sample may be less motivated by the burstiness of star formation. △ Less

Submitted 17 July, 2024; v1 submitted 11 January, 2024; originally announced January 2024.

Comments: 18 pages, 7 figures

Closing the Gap between Observed Low-Mass Galaxy HI Kinematics and CDM Predictions

Authors: Amy Sardone, Annika H. G. Peter, Alyson M. Brooks, Jane Kaczmarek

Abstract: Testing the standard cosmological model ($Λ$CDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matte… ▽ More Testing the standard cosmological model ($Λ$CDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matter halos predicted to host them, obscuring our understanding of the low-mass end of the galaxy-halo connection. We use deep HI observations of low-mass galaxies at high spectral resolution in combination with cosmological simulations of dwarf galaxies to better understand the connection between dwarf galaxy kinematics and low-mass halos. Specifically, we use HI line widths to directly compare to the maximum velocities in a dark matter halo, and find that each deeper measurement approaches the expected one-to-one relationship between the observed kinematics and the predicted kinematics in $Λ$CDM. We also measure baryonic masses and place these on the Baryonic Tully-Fisher relation (BTFR). Again, our deepest measurements approach the theoretical predictions for the low-mass end of this relation, a significant improvement on similar measurements based on line widths measured at 50\% and 20\% of the peak. Our data also hints at the rollover in the BTFR predicted by hydrodynamical simulations of $Λ$CDM for low-mass galaxies. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 13 pages of text, 4 figures, submitted to AAS Journals

The Role of Mass and Environment on Satellite distributions around Milky Way analogs in the Romulus25 simulation

Authors: Jordan Van Nest, Ferah Munshi, Charlotte Christensen, Alyson M. Brooks, Michael Tremmel, Thomas R. Quinn

Abstract: We study satellite counts and quenched fractions for satellites of Milky Way analogs in Romulus25, a large-volume cosmological hydrodynamic simulation. Depending on the definition of a Milky Way analog, we have between 66 and 97 Milky Way analogs in Romulus25, a 25 Mpc per-side uniform volume simulation. We use these analogs to quantify the effect of environment and host properties on satellite po… ▽ More We study satellite counts and quenched fractions for satellites of Milky Way analogs in Romulus25, a large-volume cosmological hydrodynamic simulation. Depending on the definition of a Milky Way analog, we have between 66 and 97 Milky Way analogs in Romulus25, a 25 Mpc per-side uniform volume simulation. We use these analogs to quantify the effect of environment and host properties on satellite populations. We find that the number of satellites hosted by a Milky Way analog increases predominantly with host stellar mass, while environment, as measured by the distance to a Milky Way-mass or larger halo, may have a notable impact in high isolation. Similarly, we find that the satellite quenched fraction for our analogs also increases with host stellar mass, and potentially in higher-density environments. These results are robust for analogs within 3 Mpc of another Milky Way-mass or larger halo, the environmental parameter space where the bulk of our sample resides. We place these results in the context of observations through comparisons to the Exploration of Local VolumE Satellites and Satellites Around Galactic Analogs surveys. Our results are robust to changes in Milky Way analog selection criteria, including those that mimic observations. Finally, as our samples naturally include Milky Way-Andromeda pairs, we examine quenched fractions in pairs vs isolated systems. We find potential evidence, though not conclusive, that pairs, defined as being within 1 Mpc of another Milky Way-mass or larger halo, may have higher satellite quenched fractions. △ Less

Submitted 26 September, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

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

AGN quenching in simulated dwarf galaxies

Authors: Ray S. Sharma, Alyson M. Brooks, Michael Tremmel, Jillian Bellovary, Thomas R. Quinn

Abstract: We examine the quenching characteristics of $328$ isolated dwarf galaxies $\left(10^{8} < M_{\rm star}/M_\odot < 10^{10} \right)$ within the \Rom{} cosmological hydrodynamic simulation. Using mock observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological sim… ▽ More We examine the quenching characteristics of $328$ isolated dwarf galaxies $\left(10^{8} < M_{\rm star}/M_\odot < 10^{10} \right)$ within the \Rom{} cosmological hydrodynamic simulation. Using mock observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies below $M_{\rm star} < 10^{9} M_\odot$ not detected within the NSA. We find that the presence of massive black holes (MBHs) in \Rom{} is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs between $z=0.5-1$, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses, $M_{\rm star} \lesssim 10^{9.3} M_\odot$, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses, $M_{\rm star} \gtrsim 10^{9.3} M_\odot$, quenching occurs rapidly within $1$ Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent over-quenching of dwarf galaxies due to MBH suggests higher resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 23 pages, 11 figures; submitted to ApJ

Emergent Microrobotic Oscillators via Asymmetry-Induced Order

Authors: Jing Fan Yang, Thomas A. Berrueta, Allan M. Brooks, Albert Tianxiang Liu, Ge Zhang, David Gonzalez-Medrano, Sungyun Yang, Volodymyr B. Koman, Pavel Chvykov, Lexy N. LeMar, Marc Z. Miskin, Todd D. Murphey, Michael S. Strano

Abstract: Spontaneous low-frequency oscillations on the order of several hertz are the drivers of many crucial processes in nature. From bacterial swimming to mammal gaits, the conversion of static energy inputs into slowly oscillating electrical and mechanical power is key to the autonomy of organisms across scales. However, the fabrication of slow artificial oscillators at micrometre scales remains a majo… ▽ More Spontaneous low-frequency oscillations on the order of several hertz are the drivers of many crucial processes in nature. From bacterial swimming to mammal gaits, the conversion of static energy inputs into slowly oscillating electrical and mechanical power is key to the autonomy of organisms across scales. However, the fabrication of slow artificial oscillators at micrometre scales remains a major roadblock towards the development of fully-autonomous microrobots. Here, we report the emergence of a low-frequency relaxation oscillator from a simple collective of active microparticles interacting at the air-liquid interface of a peroxide drop. Their collective oscillations form chemomechanical and electrochemical limit cycles that enable the transduction of ambient chemical energy into periodic mechanical motion and on-board electrical currents. Surprisingly, the collective can oscillate robustly even as more particles are introduced, but only when we add a single particle with modified reactivity to intentionally break the system's permutation symmetry. We explain such emergent order through a novel thermodynamic mechanism for asymmetry-induced order. The energy harvested from the stabilized system oscillations enables the use of on-board electronic components, which we demonstrate by cyclically and synchronously driving microrobotic arms. This work highlights a new strategy for achieving low-frequency oscillations at the microscale that are otherwise difficult to observe outside of natural systems, paving the way for future microrobotic autonomy. △ Less

Submitted 26 September, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: Main text contains 13 pages and 4 figures. Supplementary information contains 21 pages and 16 supplementary figures. For associated supplementary videos, see https://www.dropbox.com/sh/2bwenfiifqnkx3i/AABcLH2mVQ_8uPxnnbzu4rGWa?dl=0

Journal ref: Nat.Commun. 13 (2022) 5734

A hidden population of massive black holes in simulated dwarf galaxies

Authors: Ray S. Sharma, Alyson M. Brooks, Michael Tremmel, Jillian Bellovary, Angelo Ricarte, Thomas R. Quinn

Abstract: We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consisten… ▽ More We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity at $M_{\rm star} < 3\times10^{10} M_{\odot}$, $M_{\rm 200} < 3\times10^{11} M_\odot$. Local dwarf AGN in \textsc{Romulus25} follow observed scaling relations between AGN x-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable x-ray observations. Since $z=2$, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result, the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low x-ray luminosities and high contamination from x-ray binaries and the hot interstellar medium that are undetectable by current x-ray surveys. These hidden MBHs make up $76\%$ of all MBHs in local dwarf galaxies, and include many MBHs that are undermassive relative to their host galaxy's stellar mass. Their detection relies not only on greater instrument sensitivity but on better modeling of x-ray contaminants or multi-wavelength surveys. Our results indicate dwarf AGN were substantially more active in the past despite being low-luminosity today, and indicate future deep x-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at least $z=2$. △ Less

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

What's in a name? Quantifying the interplay between the Definition, Orientation, and Shape of Ultra-diffuse Galaxies Using the Romulus Simulations

Authors: Jordan D. Van Nest, F. Munshi, A. C. Wright, M. Tremmel, A. M. Brooks, D. Nagai, T. Quinn

Abstract: We explore populations of ultra-diffuse galaxies (UDGs) in isolated, satellite, and cluster environments using the Romulus25 and RomulusC simulations, including how the populations vary with UDG definition and viewing orientation. Using a fiducial definition of UDGs, we find that isolated UDGs have notably larger semi-major (b/a) and smaller semi-minor (c/a) axis ratios than their non-UDG counterp… ▽ More We explore populations of ultra-diffuse galaxies (UDGs) in isolated, satellite, and cluster environments using the Romulus25 and RomulusC simulations, including how the populations vary with UDG definition and viewing orientation. Using a fiducial definition of UDGs, we find that isolated UDGs have notably larger semi-major (b/a) and smaller semi-minor (c/a) axis ratios than their non-UDG counterparts, i.e., they are more oblate, or diskier. This is in line with previous results that adopted the same UDG definition and showed that isolated UDGs form via early, high-spin mergers. However, the choice of UDG definition can drastically affect what subsets of a dwarf population are classified as UDGs, changing the number of UDGs by up to approximately 45% of the dwarf population. We also find that a galaxy's classification as a UDG is dependent on its viewing orientation, and this dependence decreases as environmental density increases. Overall, we conclude that some definitions for UDGs used in the literature manage to isolate a specific formation mechanism for isolated dwarfs, while less restrictive definitions erase a link to the formation mechanism. Thus, how we define UDG populations must be considered if we want to understand the formation and evolution of UDGs. △ Less

Submitted 25 September, 2023; v1 submitted 30 August, 2021; originally announced August 2021.

Joint Gas and Stellar Dynamical Models of WLM: An isolated dwarf galaxy within a cored, prolate DM halo

Authors: Gigi Y. C. Leung, Ryan Leaman, Giuseppina Battaglia, Glenn van de Ven, Alyson M. Brooks, Jorge Peñarrubia, Kim A. Venn

Abstract: We present multi-tracer dynamical models of the low mass ($M_{*} \sim 10^{7}$), isolated dwarf irregular galaxy WLM in order to simultaneously constrain the inner slope of the dark matter (DM) halo density profile ($γ$) and flattening ($q_\mathrm{DM}$), and the stellar orbital anisotropy ($β_{z}, β_{r}$). For the first time, we show how jointly constraining the mass distribution from the HI gas ro… ▽ More We present multi-tracer dynamical models of the low mass ($M_{*} \sim 10^{7}$), isolated dwarf irregular galaxy WLM in order to simultaneously constrain the inner slope of the dark matter (DM) halo density profile ($γ$) and flattening ($q_\mathrm{DM}$), and the stellar orbital anisotropy ($β_{z}, β_{r}$). For the first time, we show how jointly constraining the mass distribution from the HI gas rotation curve and solving the Jeans' equations with discrete stellar kinematics leads to a factor of $\sim2$ reduction in the uncertainties on $γ$. The mass-anisotropy degeneracy is also partially broken, leading to reductions on uncertainty by $\sim 30\%$ on $M_\mathrm{vir}$ (and $\sim 70\%$ at the half-light radius) and $\sim 25\%$ on anisotropy. Our inferred value of $γ= 0.3 \pm 0.1$ is robust to the halo geometry, and in excellent agreement with predictions of stellar feedback driven DM core creation. The derived prolate geometry of the DM halo with $q_\mathrm{DM} = 2 \pm 1$ is consistent with $Λ$CDM simulations of dwarf galaxy halos. While self-interacting DM (SIDM) models with $σ/m_{X} \sim 0.6$ can reproduce this cored DM profile, the interaction events may sphericalise the halo. The simultaneously cored and prolate DM halo may therefore present a challenge for SIDM. Finally we find that the radial profile of stellar anisotropy in WLM ($β_{r}$) follows a nearly identical trend of increasing tangential anisotropy to the classical dSphs, Fornax and Sculptor. Given WLM's orbital history, this result may call into question whether such anisotropy is a consequence of tidal stripping in only one pericentric passage or if it instead is a feature of the largely self-similar formation and evolutionary pathways for some dwarf galaxies. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: arxiv version, published in MNRAS. 23 pages, 14 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2021, Volume 500, Issue 1, pp.410-429

Stellar Migration and Chemical Enrichment in the Milky Way Disc: A Hybrid Model

Authors: James W. Johnson, David H. Weinberg, Fiorenzo Vincenzo, Jonathan C. Bird, Sarah R. Loebman, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen, Emily J. Griffith

Abstract: We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradien… ▽ More We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc's abundance structure: (i) the dependence of the [O/Fe]-[Fe/H] distribution on radius $R_\text{gal}$ and midplane distance $|z|$; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with $R_\text{gal}$ and $|z|$; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with $R_\text{gal}$, $|z|$, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] $>$ 0.1; (v) a population of young and intermediate-age $α$-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age-[O/H] and age-[Fe/H] relations, with large scatter and a median age of $\sim$4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate $\sim$2 Gyr ago improve agreement with the observed age-[Fe/H] and age-[O/H] relations, but worsen agreement with the observed age-[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice). △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 24 pages; 18 figures; to be submitted to MNRAS; comments welcome

Ultra-faint dwarfs in a Milky Way context: Introducing the Mint Condition DC Justice League Simulations

Authors: Elaad Applebaum, Alyson M. Brooks, Charlotte R. Christensen, Ferah Munshi, Thomas R. Quinn, Sijing Shen, Michael Tremmel

Abstract: We present results from the "Mint" resolution DC Justice League suite of Milky Way-like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultra-faint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to $z=0$, with initial star (dark matter) particl… ▽ More We present results from the "Mint" resolution DC Justice League suite of Milky Way-like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultra-faint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to $z=0$, with initial star (dark matter) particle masses of 994 (17900) M$_\odot$, and a force resolution of 87 pc. We study the surrounding dwarfs and UFDs, and find the simulations match the observed dynamical properties of galaxies with $-3 < M_V < -19$, and reproduce the scatter seen in the size-luminosity plane for $r_h\gtrsim200$ pc. We predict the vast majority of nearby galaxies will be observable by the Vera Rubin Observatory's co-added Legacy Survey of Space and Time (LSST). We additionally show that faint dwarfs with velocity dispersions $\lesssim5$ km/s result from severe tidal stripping of the host halo. We investigate the quenching of UFDs in a hydrodynamical Milky Way context, and find that the majority of UFDs are quenched prior to interactions with the Milky Way, though some of the quenched UFDs retain their gas until infall. Additionally these simulations yield some unique dwarfs that are the first of their kind to be simulated, e.g., an HI-rich field UFD, a late-forming UFD that has structural properties similar to Crater 2, as well as a compact dwarf satellite that has no dark matter at $z=0$. △ Less

Submitted 17 November, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: 33 pages, 17 figures, ApJ accepted version

Quenching timescales of dwarf satellites around Milky Way-mass hosts

Authors: Hollis B. Akins, Charlotte R. Christensen, Alyson M. Brooks, Ferah Munshi, Elaad Applebaum, Anna Engelhardt, Lucas Chamberland

Abstract: Observations of the low-mass satellites in the Local Group have shown high fractions of gas-poor, quiescent galaxies relative to isolated dwarfs, implying that the host halo environment plays an important role in the quenching of dwarf galaxies. In this work, we present measurements of the quenched fractions and quenching timescales of dwarf satellite galaxies in the DC Justice League suite of 4 h… ▽ More Observations of the low-mass satellites in the Local Group have shown high fractions of gas-poor, quiescent galaxies relative to isolated dwarfs, implying that the host halo environment plays an important role in the quenching of dwarf galaxies. In this work, we present measurements of the quenched fractions and quenching timescales of dwarf satellite galaxies in the DC Justice League suite of 4 high-resolution cosmological zoom-in simulations of Milky Way-mass halos. We show that these simulations accurately reproduce the satellite luminosity functions of observed nearby galaxies, as well as the variation in satellite quenched fractions from $M_* \sim 10^5$ solar masses to $10^{10}$ solar masses. We then trace the histories of satellite galaxies back to $z \sim 15$, and find that many satellites with $M_* \sim 10^{6-8}$ solar masses quench within 2 Gyr of infall into the host halo, while others in the same mass range remain star-forming for as long as 5 Gyr. We show that this scatter can be explained by the satellite's gas mass and the ram pressure it feels at infall. Finally, we identify a characteristic stellar mass scale of $10^8$ solar masses above which infalling satellites are largely resistant to rapid environmental quenching. △ Less

Submitted 2 February, 2021; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 16 pages, 7 figures, accepted to ApJ

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Authors: Kartheik G. Iyer, Sandro Tacchella, Shy Genel, Christopher C. Hayward, Lars Hernquist, Alyson M. Brooks, Neven Caplar, Romeel Davé, Benedikt Diemer, John C. Forbes, Eric Gawiser, Rachel S. Somerville, Tjitske K. Starkenburg

Abstract: Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIR… ▽ More Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 31 pages, 17 figures (+ appendix). Resubmitted to MNRAS after responding to referee's comments. Comments are welcome!

The Formation of Isolated Ultra-Diffuse Galaxies in Romulus25

Authors: Anna C. Wright, Michael Tremmel, Alyson M. Brooks, Ferah Munshi, Daisuke Nagai, Ray S. Sharma, Thomas R. Quinn

Abstract: We use the \textsc{Romulus25} cosmological simulation volume to identify the largest-ever simulated sample of {\it field} ultra-diffuse galaxies (UDGs). At $z=0$, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. UDGs have moderately elevated HI masses, being 70\% (300\%) more HI-rich than typical isolated dwarf galaxi… ▽ More We use the \textsc{Romulus25} cosmological simulation volume to identify the largest-ever simulated sample of {\it field} ultra-diffuse galaxies (UDGs). At $z=0$, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. UDGs have moderately elevated HI masses, being 70\% (300\%) more HI-rich than typical isolated dwarf galaxies at luminosities brighter (fainter) than M$_\mathrm{B}$=-14. However, UDGs are consistent with the general isolated dwarf galaxy population and make up $\sim$20\% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The HI masses, effective radii, and overall appearances of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at $z=0$, the UDGs in our sample are the products of major mergers. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier - almost never occurring after $z=1$, produce a temporary boost in spin, and cause star formation to be redistributed to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii, producing steeper negative g-r color gradients. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, demonstrating that UDGs can potentially be formed through multiple mechanisms. △ Less

Submitted 7 January, 2021; v1 submitted 15 May, 2020; originally announced May 2020.

Comments: 20 pages, 16 figures, 2 movies, accepted to MNRAS

Self-Interacting Dark Matter and the Delay of Super-Massive Black Hole Growth

Authors: Akaxia Cruz, Andrew Pontzen, Marta Volonteri, Thomas R. Quinn, Michael Tremmel, Alyson M. Brooks, Nicole N. Sanchez, Ferah Munshi, Arianna Di Cintio

Abstract: Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find t… ▽ More Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model. △ Less

Submitted 17 April, 2020; originally announced April 2020.

The Formation of Ultra-Diffuse Galaxies in the RomulusC Galaxy Cluster Simulation

Authors: Michael Tremmel, Anna C. Wright, Alyson M. Brooks, Ferah Munshi, Daisuke Nagai, Thomas R. Quinn

Abstract: We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate… ▽ More We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate from the overall cluster dwarf population. UDGs in cluster environments form primarily from dwarf galaxies that experienced early cluster in-fall and subsequent quenching due to ram pressure. The ensuing dimming of these dwarf galaxies due to passive stellar evolution results in a population of very low surface brightness galaxies that are otherwise typical dwarfs. UDGs and non-UDGs alike are affected by tidal interactions with the cluster potential. Tidal stripping of dark matter, as well as mass loss from stellar evolution, results in the adiabatic expansion of stars, particularly in the lowest mass dwarfs. High mass dwarf galaxies show signatures of tidal heating while low mass dwarfs that survive until $z=0$ typically have not experienced such impulsive interactions. There is little difference between UDGs and non-UDGs in terms of their dark matter halos, stellar morphology, colors, and location within the cluster. In most respects cluster UDG and non-UDGs alike are similar to isolated dwarf galaxies, except for the fact that they are typically quenched. △ Less

Submitted 13 July, 2020; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 26 Pages, 22 Figures. Accepted for publication in MNRAS

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope

Authors: Alex Drlica-Wagner, Yao-Yuan Mao, Susmita Adhikari, Robert Armstrong, Arka Banerjee, Nilanjan Banik, Keith Bechtol, Simeon Bird, Kimberly K. Boddy, Ana Bonaca, Jo Bovy, Matthew R. Buckley, Esra Bulbul, Chihway Chang, George Chapline, Johann Cohen-Tanugi, Alessandro Cuoco, Francis-Yan Cyr-Racine, William A. Dawson, Ana Díaz Rivero, Cora Dvorkin, Denis Erkal, Christopher D. Fassnacht, Juan García-Bellido, Maurizio Giannotti , et al. (75 additional authors not shown)

Abstract: Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d… ▽ More Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ . △ Less

Submitted 24 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 96 pages, 22 figures, 1 table

Report number: FERMILAB-PUB-19-048-A-AE

Understanding Dwarf Galaxies in order to Understand Dark Matter

Authors: Alyson M. Brooks

Abstract: Much progress has been made in recent years by the galaxy simulation community in making realistic galaxies, mostly by more accurately capturing the effects of baryons on the structural evolution of dark matter halos at high resolutions. This progress has altered theoretical expectations for galaxy evolution within a Cold Dark Matter (CDM) model, reconciling many earlier discrepancies between theo… ▽ More Much progress has been made in recent years by the galaxy simulation community in making realistic galaxies, mostly by more accurately capturing the effects of baryons on the structural evolution of dark matter halos at high resolutions. This progress has altered theoretical expectations for galaxy evolution within a Cold Dark Matter (CDM) model, reconciling many earlier discrepancies between theory and observations. Despite this reconciliation, CDM may not be an accurate model for our Universe. Much more work must be done to understand the predictions for galaxy formation within alternative dark matter models. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: Refereed contribution to the Proceedings of the Simons Symposium on Illuminating Dark Matter, to be published by Springer

Roll of the Dice: A Stochastically Sampled IMF Alters the Stellar Content of Simulated Dwarf Galaxies

Authors: Elaad Applebaum, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen

Abstract: Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-fain… ▽ More Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to $z=6$, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in halos with mass $\lesssim10^{8.5}$ M$_\odot$, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ${\sim}$10 per cent higher star formation and ${\sim}$30 per cent higher cold gas content. Going forward, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses. △ Less

Submitted 20 February, 2020; v1 submitted 31 October, 2018; originally announced November 2018.

Comments: 14 pages, 11 figures, Updated to reflect accepted version, MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2020, 492, 8-21

Dancing in the Dark: Uncertainty in ultra-faint dwarf galaxy predictions from cosmological simulations

Authors: Ferah Munshi, Alyson M. Brooks, Charlotte Christensen, Elaad Applebaum, Kelly Holley-Bockelmann, Thomas R. Quinn, James Wadsley

Abstract: The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation… ▽ More The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation model. The total number of dwarf galaxies that forms is different by a factor of 2 between the two runs, but most of these are satellites, leading to a factor of 5 difference in the number of luminous UFD companions around more massive, isolated dwarfs. The first run yields a 47\% chance of finding a satellite around a M$_{halo}$ $\sim 10^{10}$ M$_{\odot}$ host, while the H$_2$ run predicts only a 16\% chance. Metallicity is the primary physical parameter that creates this difference. As metallicity decreases, the formation of H$_2$ is slowed and relegated to higher-density material. Thus, our H$_2$ run is unable to form many (and often, any) stars before reionization removes gas. These results emphasize that predictions for UFD properties made using hydrodynamic simulations, in particular regarding the frequency of satellites around dwarf galaxies, the slope of the stellar mass function at low masses, as well as the properties of ultra-faint galaxies occupying the smallest halos, are extremely sensitive to the subgrid physics of star formation contained within the simulation. However, upcoming discoveries of ultra-faint dwarfs will provide invaluable constraining power on the physics of the first star formation. △ Less

Submitted 8 February, 2019; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: 18 pages, 7 Figures. Accepted version, ApJ

Reignition of Star Formation in Dwarf Galaxies

Authors: Anna C. Wright, Alyson M. Brooks, Daniel R. Weisz, Charlotte R. Christensen

Abstract: The Local Group hosts a number of star-forming dwarf galaxies that show evidence of periods of little to no star formation. We use a suite of cosmological simulations to study how star formation is reignited in such galaxies. We focus on isolated galaxies at $z=0$ with halo masses between 9.2$\times$10$^8$ M$_\odot$ and 8.4$\times$10$^9$ M$_\odot$, where star formation is typically shut off by rei… ▽ More The Local Group hosts a number of star-forming dwarf galaxies that show evidence of periods of little to no star formation. We use a suite of cosmological simulations to study how star formation is reignited in such galaxies. We focus on isolated galaxies at $z=0$ with halo masses between 9.2$\times$10$^8$ M$_\odot$ and 8.4$\times$10$^9$ M$_\odot$, where star formation is typically shut off by reionization or by supernova feedback. Nearly 20% of these simulated galaxies later restart star formation, due to interactions with streams of gas in the intergalactic medium, indicating that this mechanism is relatively common in this mass range and that many isolated dwarfs at $z=0$ may not have been isolated throughout their histories. The source of this gas is not necessarily cosmic filaments. Rather, the dwarfs interact with gas thrown off by nearby galaxy mergers or streams extending from neighboring galaxies. While high ram pressure interactions of this nature lead to stripping, the encounters that reignite star formation are low density and/or low velocity and thus low ram pressure, resulting in compression of the hot gas in the halos of our dwarfs. The gas mass bound up in hot halos can be substantial -- at least an order of magnitude greater than the mass contained in HI. Consequently, we find that dwarfs that have experienced reignition tend to be more HI-rich and have a higher M$_{HI}$/M$_{*}$ ratio at $z=0$ than galaxies with continuous star formation. Using this fact, we identify galaxies in the Local Volume that might have "gappy" star formation histories, and can be studied by the Hubble Space Telescope or the James Webb Space Telescope. △ Less

Submitted 11 October, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: 16 pages, 11 figures, 1 table, 2 movies (online), accepted to MNRAS

A unified model for age-velocity dispersion relations in Local Group galaxies: Disentangling ISM turbulence and latent dynamical heating

Authors: Ryan Leaman, J. Trevor Mendel, Emily Wisnioski, Alyson M. Brooks, Michael A. Beasley, Else Starkenburg, Marie Martig, Giuseppina Battaglia, Charlotte Christensen, Andrew A. Cole, T. J. L. de Boer, Drew Wills

Abstract: We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $σ$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born o… ▽ More We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $σ$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born out of a more turbulent ISM at high redshift and retain that larger velocity dispersion till present day - consistent with recent IFU studies. To ascertain the dominant mechanism and better understand the impact of instabilities and feedback, we develop models based on observed SFHs of these Local Group galaxies in order to create an evolutionary formalism which describes the ISM velocity dispersion due to a galaxy's evolving gas fraction. These empirical models relax the common assumption that the stars are born from gas which has constant velocity dispersion at all redshifts. Using only the observed SFHs as input, the ISM velocity dispersion and a mid-plane scattering model fits the observed AVRs of low mass galaxies without fine tuning. Higher mass galaxies above $M_{vir} > 10^{11}$ need a larger contribution from latent dynamical heating processes (for example minor mergers), in excess of the ISM model. Using the SFHs we also find that supernovae feedback does not appear to be a dominant driver of the gas velocity dispersion compared to gravitational instabilities - at least for dispersions $σ\gtrsim 25$ km/s. Together our results point to stars being born with a velocity dispersion close to that of the gas at the time of their formation, with latent dynamical heating operating with a galaxy mass-dependent efficiency. These semi-empirical relations may help constrain the efficiency of feedback and its impact on the physics of disk settling in galaxy formation simulations. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: 21 pages, 10 figures. accepted for publication in MNRAS

Going, going, gone dark: Quantifying the scatter in the faintest dwarf galaxies

Authors: Ferah Munshi, Alyson M. Brooks, Elaad Applebaum, Daniel R. Weisz, Fabio Governato, Thomas R. Quinn

Abstract: We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies,… ▽ More We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies, run with the N-body + SPH code, ChaNGA, at a high resolution of 60 pc. For M$_{halo} > 10^9$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter. However, the scatter in the SMHM relation increases dramatically for lower-mass halos. We find that some of this scatter is due to {\em dark dwarfs}, halos devoid of stars. However, even when only considering well-resolved halos that contain a stellar population, the scatter in stellar mass reaches nearly 1 dex for M$_{halo}$($z=0$) 10$^7$ M$_{\odot}$. Much of this scatter is due to including satellites of the dwarf galaxies that have had their halo masses reduced through tidal stripping. The fraction of dark dwarfs (those that contain no stars) increases substantially with decreasing halo mass. When these dark halos are considered, the true scatter in the SMHM at low masses is even larger. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. We provide a formula to stochastically populate low-mass halos following our results. Our predicted large scatter at low halo masses increases the slope of the resulting stellar mass function on the ultra-faint dwarf galaxy scales currently being probed by such surveys as the Dark Energy Survey or the Hyper-Suprime Cam Subaru Strategic Program, and in the future by the Large Synoptic Survey Telescope. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 13 pages, 5 figures

How to Reconcile the Observed Velocity Function of Galaxies with Theory

Authors: Alyson M. Brooks, Emmanouil Papastergis, Charlotte R. Christensen, Fabio Governato, Adrienne Stilp, Thomas R. Quinn, James Wadsley

Abstract: Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-t… ▽ More Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-to-apples comparison yields an "observed" velocity function in agreement with observations, reconciling the large number of low-mass halos expected in a LCDM cosmological model with the low number of observed dwarfs at a given velocity. We then explore the source of the discrepancy between observations and theory, and conclude that the dearth of observed dwarf galaxies is primarily explained by two effects. The first effect is that galactic rotational velocities derived from the HI linewidth severely underestimate the maximum halo velocity. The second effect is that a large fraction of halos at the lowest masses are too faint to be detected by current galaxy surveys. We find that cored dark matter density profiles can contribute to the lower observed velocity of galaxies, but only for galaxies in which the velocity is measured interior to the size of the core (~3 kpc). △ Less

Submitted 25 October, 2017; v1 submitted 26 January, 2017; originally announced January 2017.

Comments: updated after acceptance for publication in ApJ

The Role of Baryons in Creating Statistically Significant Planes of Satellites around Milky Way-Mass Galaxies

Authors: Sheehan H. Ahmed, Alyson M. Brooks, Charlotte R. Christensen

Abstract: We investigate whether the inclusion of baryonic physics influences the formation of thin, coherently rotating planes of satellites such as those seen around the Milky Way and Andromeda. For four Milky Way-mass simulations, each run both as dark matter-only and with baryons included, we are able to identify a planar configuration that significantly maximizes the number of plane satellite members.… ▽ More We investigate whether the inclusion of baryonic physics influences the formation of thin, coherently rotating planes of satellites such as those seen around the Milky Way and Andromeda. For four Milky Way-mass simulations, each run both as dark matter-only and with baryons included, we are able to identify a planar configuration that significantly maximizes the number of plane satellite members. The maximum plane member satellites are consistently different between the dark matter-only and baryonic versions of the same run due to the fact that satellites are both more likely to be destroyed and to infall later in the baryonic runs. Hence, studying satellite planes in dark matter-only simulations is misleading, because they will be composed of different satellite members than those that would exist if baryons were included. Additionally, the destruction of satellites in the baryonic runs leads to less radially concentrated satellite distributions, a result that is critical to making planes that are statistically significant compared to a random distribution. Since all planes pass through the centre of the galaxy, it is much harder to create a plane from a random distribution if the satellites have a low radial concentration. We identify Andromeda's low radial satellite concentration as a key reason why the plane in Andromeda is highly significant. Despite this, when co-rotation is considered, none of the satellite planes identified for the simulated galaxies are as statistically significant as the observed planes around the Milky Way and Andromeda, even in the baryonic runs. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Comments: 16 pages, 9 figures, submitted to MNRAS

Search for Gamma-ray Emission from Dark Matter Annihilation in the Small Magellanic Cloud with the Fermi Large Area Telescope

Authors: Regina Caputo, Matthew R. Buckley, Pierrick Martin, Eric Charles, Alyson M. Brooks, Alex Drlica-Wagner, Jennifer M. Gaskins, Matthew Wood

Abstract: The Small Magellanic Cloud (SMC) is the second-largest satellite galaxy of the Milky Way and is only 60 kpc away. As a nearby, massive, and dense object with relatively low astrophysical backgrounds, it is a natural target for dark matter indirect detection searches. In this work, we use six years of Pass 8 data from the Fermi Large Area Telescope to search for gamma-ray signals of dark matter ann… ▽ More The Small Magellanic Cloud (SMC) is the second-largest satellite galaxy of the Milky Way and is only 60 kpc away. As a nearby, massive, and dense object with relatively low astrophysical backgrounds, it is a natural target for dark matter indirect detection searches. In this work, we use six years of Pass 8 data from the Fermi Large Area Telescope to search for gamma-ray signals of dark matter annihilation in the SMC. Using data-driven fits to the gamma-ray backgrounds, and a combination of N-body simulations and direct measurements of rotation curves to estimate the SMC DM density profile, we found that the SMC was well described by standard astrophysical sources, and no signal from dark matter annihilation was detected. We set conservative upper limits on the dark matter annihilation cross section. These constraints are in agreement with stronger constraints set by searches in the Large Magellanic Cloud and approach the canonical thermal relic cross section at dark matter masses lower than 10 GeV in the $b\bar{b}$ and $τ^+τ^-$ channels. △ Less

Submitted 29 March, 2016; v1 submitted 2 March, 2016; originally announced March 2016.

Comments: 17 pages, 11 figures. Accepted by PRD

Journal ref: Phys. Rev. D 93, 062004 (2016)

Assessing Astrophysical Uncertainties in Direct Detection with Galaxy Simulations

Authors: Jonathan D. Sloane, Matthew R. Buckley, Alyson M. Brooks, Fabio Governato

Abstract: We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the Solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence… ▽ More We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the Solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions which are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Further, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulation. We show the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way's estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models which are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the Solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame. △ Less

Submitted 29 July, 2016; v1 submitted 20 January, 2016; originally announced January 2016.

Comments: 10 pages, 4 figures Updated after referee comments

Bulge Formation via Mergers in Cosmological Simulations

Authors: Alyson M. Brooks, Charlotte R. Christensen

Abstract: The latest generation of cosmological simulations are on the verge of being able to resolve the structure of bulges for the first time. Hence, we review the current state of bulge formation in cosmological simulations, and discuss open questions that can be addressed in the near future by simulators, with a particular focus on merger-driven bulge growth. Galaxy mergers have long been assumed to pr… ▽ More The latest generation of cosmological simulations are on the verge of being able to resolve the structure of bulges for the first time. Hence, we review the current state of bulge formation in cosmological simulations, and discuss open questions that can be addressed in the near future by simulators, with a particular focus on merger-driven bulge growth. Galaxy mergers have long been assumed to produce classical bulges in disk galaxies. Under this bulge-formation model, though, the high rates of mergers in Cold Dark Matter galaxy formation theory predict many more classical bulges than are observed. Furthermore, simulations of galaxy formation continue to generally produce too massive of bulges. Feedback offers a promising avenue for reducing merger-driven bulge growth by maintaining high gas fractions in galaxies and ejecting low-angular momentum gas driven to the centers of galaxies. After reviewing the results of relevant research that has been published to date, we use cosmological simulations to explore the ability of feedback to reduce or even prevent bulge growth during mergers. In dwarf galaxies, mergers actually reduce the central concentration of galaxies as the induced burst of star formation drives out low angular momentum material. This result shows the potential for feedback to reduce central mass growth. However, we also demonstrate that it is very difficult for current stellar feedback models to reproduce the small bulges observed in more massive disk galaxies like the Milky Way. We argue that feedback models need to be improved, or an additional source of feedback such as AGN is necessary to generate the required outflows. △ Less

Submitted 12 November, 2015; originally announced November 2015.

Comments: The first 5 sections are review, but Section 6 introduces new analysis. We demonstrate that mergers prevent bulge formation in low mass galaxies. Approximately 15 pages in MNRAS format

Journal ref: Review article written for "Galactic Bulges," published by Springer, 2015, edited by E. Laurikainen, R.F. Peletier, D.A. Gadotti

Particle tagging and its implications for stellar population dynamics

Authors: Theo Le Bret, Andrew Pontzen, Andrew P. Cooper, Carlos Frenk, Adi Zolotov, Alyson M. Brooks, Fabio Governato, Owen H. Parry

Abstract: We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our stud… ▽ More We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: 10 pages, 6 figures

Search for Gamma-ray Emission from Dark Matter Annihilation in the Large Magellanic Cloud with the Fermi Large Area Telescope

Authors: Matthew R. Buckley, Eric Charles, Jennifer M. Gaskins, Alyson M. Brooks, Alex Drlica-Wagner, Pierrick Martin, Geng Zhao

Abstract: At a distance of 50 kpc and with a dark matter mass of $\sim10^{10}$ M$_{\odot}$, the Large Magellanic Cloud (LMC) is a natural target for indirect dark matter searches. We use five years of data from the Fermi Large Area Telescope (LAT) and updated models of the gamma-ray emission from standard astrophysical components to search for a dark matter annihilation signal from the LMC. We perform a rot… ▽ More At a distance of 50 kpc and with a dark matter mass of $\sim10^{10}$ M$_{\odot}$, the Large Magellanic Cloud (LMC) is a natural target for indirect dark matter searches. We use five years of data from the Fermi Large Area Telescope (LAT) and updated models of the gamma-ray emission from standard astrophysical components to search for a dark matter annihilation signal from the LMC. We perform a rotation curve analysis to determine the dark matter distribution, setting a robust minimum on the amount of dark matter in the LMC, which we use to set conservative bounds on the annihilation cross section. The LMC emission is generally very well described by the standard astrophysical sources, with at most a $1-2σ$ excess identified near the kinematic center of the LMC once systematic uncertainties are taken into account. We place competitive bounds on the dark matter annihilation cross section as a function of dark matter particle mass and annihilation channel. △ Less

Submitted 15 April, 2015; v1 submitted 3 February, 2015; originally announced February 2015.

Comments: 33 pages, 22 figures Version 2: minor corrections and clarifications after journal peer review process

All about baryons: revisiting SIDM predictions at small halo masses

Authors: A. Bastidas Fry, F. Governato, A. Pontzen, T. Quinn, M. Tremmel, L. Anderson, H. Menon, A. M. Brooks, J. Wadsley

Abstract: We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $Λ$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that natu… ▽ More We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $Λ$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that naturally drives potential fluctuations able to create dark matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density at halo peak velocities V$_{max}$ $<$ 30 Km/s. This is due to the fact that the central regions of very low--mass field halos have relatively low central velocity dispersion and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM halos with V$_{max}$ $<$ 30 km/s have inefficient star formation, and hence weak supernova feedback. At a fixed 2 cm2/g SIDM cross section, the DM content of very low mass CDM and SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses ($\sim$ 10$^{10}$ solar masses), the introduction of baryonic processes creates field dwarf galaxies with dark matter cores and central DM$+$baryon distributions that are effectively indistinguishable between CDM and SIDM. Both models are in broad agreement with observed Local Group field galaxies across the range of masses explored. To significantly differentiate SIDM from CDM at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases to values larger than 2 $cm^{2}/g$ for halos with V$_{max}$ < 25-30 Km/s needs to be introduced. △ Less

Submitted 16 June, 2015; v1 submitted 2 January, 2015; originally announced January 2015.

Comments: 13 pages, 12 Figures, one Appendix, MNRAS in press. Two figures added in final version showing the simulated systems on the stellar mass halo mass relation and the SIDM vs CDM mass deficiency at a function of radius

Journal ref: Monthly Notices of the Royal Astronomical Society 2015, 452, 1468

Delayed Star Formation in Isolated Dwarf Galaxies: HST Star Formation History of the Aquarius Dwarf Irregular

Authors: Andrew A. Cole, Daniel R. Weisz, Andrew E. Dolphin, Evan D. Skillman, Alan W. McConnachie, Alyson M. Brooks, Ryan Leaman

Abstract: We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing… ▽ More We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ~10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ~10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ~2). The star formation rate increased dramatically ~6-8 Gyr ago (z ~ 0.7-1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar M(HI)/M(stellar), dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies. △ Less

Submitted 4 September, 2014; originally announced September 2014.

Comments: 13 pages, 9 figures (4 in color), accepted by ApJ

Consequences of bursty star formation on galaxy observables at high redshifts

Authors: Alberto Domínguez, Brian Siana, Alyson M. Brooks, Charlotte R. Christensen, Gustavo Bruzual, Daniel P. Stark, Anahita Alavi

Abstract: The star formation histories (SFHs) of dwarf galaxies are thought to be \emph{bursty}, with large -- order of magnitude -- changes in the star formation rate on timescales similar to O-star lifetimes. As a result, the standard interpretations of many galaxy observables (which assume a slowly varying SFH) are often incorrect. Here, we use the SFHs from hydro-dynamical simulations to investigate the… ▽ More The star formation histories (SFHs) of dwarf galaxies are thought to be \emph{bursty}, with large -- order of magnitude -- changes in the star formation rate on timescales similar to O-star lifetimes. As a result, the standard interpretations of many galaxy observables (which assume a slowly varying SFH) are often incorrect. Here, we use the SFHs from hydro-dynamical simulations to investigate the effects of bursty SFHs on sample selection and interpretation of observables and make predictions to confirm such SFHs in future surveys. First, because dwarf galaxies' star formation rates change rapidly, the mass-to-light ratio is also changing rapidly in both the ionizing continuum and, to a lesser extent, the non-ionizing UV continuum. Therefore, flux limited surveys are highly biased toward selecting galaxies in the \emph{burst} phase and very deep observations are required to detect all dwarf galaxies at a given stellar mass. Second, we show that a $\log_{10}[νL_ν(1500{\rm Å})/L_{{\rm H}α}]>2.5$ implies a very recent quenching of star formation and can be used as evidence of stellar feedback regulating star formation. Third, we show that the ionizing continuum can be significantly higher than when assuming a constant SFH, which can affect the interpretation of nebular emission line equivalent widths and direct ionizing continuum detections. Finally, we show that a star formation rate estimate based on continuum measurements only (and not on nebular tracers such as the hydrogen Balmer lines) will not trace the rapid changes in star formation and will give the false impression of a star-forming main sequence with low dispersion. △ Less

Submitted 15 May, 2015; v1 submitted 25 August, 2014; originally announced August 2014.

Comments: 12 pages, 10 figures, and 2 tables. Version accepted by MNRAS

The Milky Way Tomography with SDSS. V. Mapping the Dark Matter Halo

Authors: Sarah R. Loebman, Zeljko Ivezic, Thomas R. Quinn, Jo Bovy, Charlotte R. Christensen, Mario Juric, Rok Roskar, Alyson M. Brooks, Fabio Governato

Abstract: We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to genera… ▽ More We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior. △ Less

Submitted 24 November, 2014; v1 submitted 22 August, 2014; originally announced August 2014.

Comments: 15 pages, 21 figures, accepted to ApJ. If viewing using Preview, some figures may appear blurry. If possible, please view using Adobe PDF Viewer