Showing 1–50 of 62 results for author: Nadler, E O

Dark Matter Coupled to Radiation: Limits from the Milky Way Satellites

Authors: Wendy Crumrine, Ethan O. Nadler, Rui An, Vera Gluscevic

Abstract: Interactions between dark matter (DM) and relativistic particles at early times suppress structure formation on small scales. In particular, the scattering process transfers heat and momentum from radiation to DM, ultimately reducing the abundance of low-mass DM halos and the dwarf galaxies they host. Herein, we derive limits on DM--photon and DM--neutrino scattering cross section using the Milky… ▽ More Interactions between dark matter (DM) and relativistic particles at early times suppress structure formation on small scales. In particular, the scattering process transfers heat and momentum from radiation to DM, ultimately reducing the abundance of low-mass DM halos and the dwarf galaxies they host. Herein, we derive limits on DM--photon and DM--neutrino scattering cross section using the Milky Way (MW) satellite galaxy population. We consider temperature-independent interactions parameterized by DM mass ($m_χ$) and DM--radiation interaction cross section ($σ_{χ\text{--}i}$, where $i$ represents the target species). By requiring that the linear matter power spectra are strictly less suppressed than in the case of a thermal-relic warm DM, we derive the following $95\%$ upper limits at $m_χ=1$ MeV: $σ_{χ\text{--}γ}<1.50\times10^{-38}\text{cm}^2$ and $σ_{χ\text{--}ν}<2.41\times10^{-38}\text{cm}^2$. Our bounds on $σ_{χ\text{--}i}$ depend linearly on $m_χ$ for $m_χ\gtrsim 1~\mathrm{MeV}$ and improve upon previous limits by an order of magnitude. The mass dependence of our limit approaches $m_χ^3$ at lower masses due to the effects of DM sound speed; at $m_χ=100~\mathrm{keV}$, we arrive at an upper limit over three orders of magnitude more stringent than achieved in previous explorations. Upcoming dwarf galaxy surveys will further improve the sensitivity of similar analyses, complementing laboratory and indirect detection searches for DM--radiation interactions. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 7 pages, 2 figures, 1 table, submitted to PRD

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Authors: Daneng Yang, Ethan O. Nadler, Hai-Bo Yu

Abstract: We systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in arXiv:2305.16176 with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and suhbhalos. We focus on the accuracy of the model in predicting halo d… ▽ More We systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in arXiv:2305.16176 with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and suhbhalos. We focus on the accuracy of the model in predicting halo density profiles at $z=0$ and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within $10\%-50\%$ for most of the simulated (sub)halos, $50\%-100\%$ for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts. △ Less

Submitted 15 June, 2024; originally announced June 2024.

Comments: 20 pages, 19 figures

The SAGA Survey. V. Modeling Satellite Systems around Milky Way-mass Galaxies with Updated UniverseMachine

Authors: Yunchong Wang, Ethan O. Nadler, Yao-Yuan Mao, Risa H. Wechsler, Tom Abel, Peter Behroozi, Marla Geha, Yasmeen Asali, Mithi A. C. de los Reyes, Erin Kado-Fong, Nitya Kallivayalil, Erik J. Tollerud, Benjamin Weiner, John F. Wu

Abstract: Environment plays a critical role in shaping the assembly of low-mass galaxies. Here, we use the UniverseMachine (UM) galaxy-halo connection framework and the Data Release 3 of the Satellites Around Galactic Analogs (SAGA) Survey to place dwarf galaxy star formation and quenching into a cosmological context. UM is a data-driven forward model that flexibly parameterizes galaxy star formation rates… ▽ More Environment plays a critical role in shaping the assembly of low-mass galaxies. Here, we use the UniverseMachine (UM) galaxy-halo connection framework and the Data Release 3 of the Satellites Around Galactic Analogs (SAGA) Survey to place dwarf galaxy star formation and quenching into a cosmological context. UM is a data-driven forward model that flexibly parameterizes galaxy star formation rates (SFR) using only halo mass and assembly history. We add a new quenching model to UM, tailored for galaxies with stellar masses $\lesssim 10^9$ solar masses, and constrain the model down to a stellar mass $\gtrsim 10^7$ solar masses using new SAGA observations of 101 satellite systems around Milky Way (MW)-mass hosts and a sample of isolated field galaxies in a similar mass range from the Sloan Digital Sky Survey (SDSS). The new best-fit model, 'UM-SAGA,' reproduces the satellite stellar mass functions, average SFRs, and quenched fractions in SAGA satellites while keeping isolated dwarfs mostly star forming. The enhanced quenching in satellites relative to isolated field galaxies leads the model to maximally rely on halo assembly to explain the observed environmental quenching. Extrapolating the model down to a stellar mass $\sim 10^{6.5}$ solar masses yields a quenched fraction of $\gtrsim$ 30% for isolated field galaxies and $\gtrsim$ 80% for satellites of MW-mass hosts at this stellar mass. This specific prediction can soon be tested by spectroscopic surveys to reveal the relative importance of internal feedback, cessation of mass and gas accretion, satellite-specific gas processes, and reionization for the evolution of faint low-mass galaxies. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 33 pages, 14 figures, 2 tables. This paper is part of the SAGA Survey Data Release 3. Survey website: https://sagasurvey.org

The SAGA Survey. IV. The Star Formation Properties of 101 Satellite Systems around Milky Way-mass Galaxies

Authors: Marla Geha, Yao-Yuan Mao, Risa H. Wechsler, Yasmeen Asali, Erin Kado-Fong, Nitya Kallivayalil, Ethan O. Nadler, Erik J. Tollerud, Benjamin Weiner, Mithi A. C. de los Reyes, Yunchong Wang, John F. Wu

Abstract: We present the star-forming properties of 378 satellite galaxies around 101 Milky Way analogs in the Satellites Around Galactic Analogs (SAGA) Survey, focusing on the environmental processes that suppress or quench star formation. In the SAGA stellar mass range of 10^6 to 10^10 solar masses, we present quenched fractions, star-forming rates, gas-phase metallicities, and gas content. The fraction o… ▽ More We present the star-forming properties of 378 satellite galaxies around 101 Milky Way analogs in the Satellites Around Galactic Analogs (SAGA) Survey, focusing on the environmental processes that suppress or quench star formation. In the SAGA stellar mass range of 10^6 to 10^10 solar masses, we present quenched fractions, star-forming rates, gas-phase metallicities, and gas content. The fraction of SAGA satellites that are quenched increases with decreasing stellar mass and shows significant system-to-system scatter. SAGA satellite quenched fractions are highest in the central 100 kpc of their hosts and decline out to the virial radius. Splitting by specific star formation rate (sSFR), the least star-forming satellite quartile follows the radial trend of the quenched population. The median sSFR of star-forming satellites increases with decreasing stellar mass and is roughly constant with projected radius. Star-forming SAGA satellites are consistent with the star formation rate--stellar mass relationship determined in the Local Volume, while the median gas-phase metallicity is higher and median HI gas mass is lower at all stellar masses. We investigate the dependence of the satellite quenched fraction on host properties. Quenched fractions are higher in systems with larger host halo mass, but this trend is only seen in the inner 100 kpc; we do not see significant trends with host color or star formation rate. Our results suggest that lower mass satellites and satellites inside 100 kpc are more efficiently quenched in a Milky Way-like environment, with these processes acting sufficiently slowly to preserve a population of star-forming satellites at all stellar masses and projected radii. △ Less

Submitted 25 July, 2024; v1 submitted 22 April, 2024; originally announced April 2024.

Comments: 25 pages, 13 figures, 2 tables. Accepted by ApJ. This paper is part of the SAGA Survey Data Release 3. Data and notebooks are available on the survey website: https://sagasurvey.org

The SAGA Survey. III. A Census of 101 Satellite Systems around Milky Way-mass Galaxies

Authors: Yao-Yuan Mao, Marla Geha, Risa H. Wechsler, Yasmeen Asali, Yunchong Wang, Erin Kado-Fong, Nitya Kallivayalil, Ethan O. Nadler, Erik J. Tollerud, Benjamin Weiner, Mithi A. C. de los Reyes, John F. Wu

Abstract: We present the third Data Release (DR3) of the Satellites Around Galactic Analogs (SAGA) Survey, a spectroscopic survey characterizing satellite galaxies around Milky Way (MW)-mass galaxies. The SAGA Survey DR3 includes 378 satellites identified across 101 MW-mass systems in the distance range 25-40.75 Mpc, and an accompanying redshift catalog of background galaxies (including about 46,000 taken b… ▽ More We present the third Data Release (DR3) of the Satellites Around Galactic Analogs (SAGA) Survey, a spectroscopic survey characterizing satellite galaxies around Milky Way (MW)-mass galaxies. The SAGA Survey DR3 includes 378 satellites identified across 101 MW-mass systems in the distance range 25-40.75 Mpc, and an accompanying redshift catalog of background galaxies (including about 46,000 taken by SAGA) in the SAGA footprint of 84.7 sq. deg. The number of confirmed satellites per system ranges from zero to 13, in the stellar mass range 10^6 to 10^10 solar masses. Based on a detailed completeness model, this sample accounts for 94% of the true satellite population down to a stellar mass of 10^7.5 solar masses. We find that the mass of the most massive satellite in SAGA systems is the strongest predictor of satellite abundance; one-third of the SAGA systems contain LMC-mass satellites, and they tend to have more satellites than the MW. The SAGA satellite radial distribution is less concentrated than the MW, and the SAGA quenched fraction below 10^8.5 solar masses is lower than the MW, but in both cases, the MW is within 1 sigma of SAGA system-to-system scatter. SAGA satellites do not exhibit a clear corotating signal as has been suggested in the MW/M31 satellite systems. Although the MW differs in many respects from the typical SAGA system, these differences can be reconciled if the MW is an older, slightly less massive host with a recently accreted LMC/SMC system. △ Less

Submitted 25 July, 2024; v1 submitted 22 April, 2024; originally announced April 2024.

Comments: 38 pages, 23 figures, 7 tables. Accepted by ApJ. This paper is part of the SAGA Survey Data Release 3. Data and notebooks will be available on the survey website: https://sagasurvey.org

Milky Way-est: Cosmological Zoom-in Simulations with Large Magellanic Cloud and Gaia-Sausage-Enceladus Analogs

Authors: Deveshi Buch, Ethan O. Nadler, Risa H. Wechsler, Yao-Yuan Mao

Abstract: We present Milky Way-est, a suite of 20 cosmological cold-dark-matter-only zoom-in simulations of Milky Way (MW)-like host halos. Milky Way-est hosts are selected such that they: ($i$) are consistent with the MW's measured halo mass and concentration, ($ii$) accrete a Large Magellanic Cloud (LMC)-like ($\approx 10^{11}~M_{\odot}$) subhalo within the last $1.3~\mathrm{Gyr}$ on a realistic orbit, pl… ▽ More We present Milky Way-est, a suite of 20 cosmological cold-dark-matter-only zoom-in simulations of Milky Way (MW)-like host halos. Milky Way-est hosts are selected such that they: ($i$) are consistent with the MW's measured halo mass and concentration, ($ii$) accrete a Large Magellanic Cloud (LMC)-like ($\approx 10^{11}~M_{\odot}$) subhalo within the last $1.3~\mathrm{Gyr}$ on a realistic orbit, placing them near $50~\mathrm{kpc}$ from the host center at $z\approx 0$, and ($iii$) undergo a $>$1:5 sub-to-host halo mass ratio merger with a Gaia-Sausage-Enceladus (GSE)-like system at early times ($0.67<z<3$). Hosts satisfying these LMC and GSE constraints constitute $< 1\%$ of all halos in the MW's mass range, and their total masses grow rapidly at late times due to LMC analog accretion. Compared to hosts of a similar final halo mass that are not selected to include LMC and GSE analogs, Milky Way-est hosts contain $22\%$ more subhalos with present-day virial masses above $10^8~M_{\odot}$ throughout the virial radius, on average. This enhancement reaches $\approx 80\%$ in the inner $100~\mathrm{kpc}$ and is largely, if not entirely, due to LMC-associated subhalos. These systems also induce spatial anisotropy in Milky Way-est subhalo populations, with $\approx 60\%$ of the total subhalo population within $100~\mathrm{kpc}$ found in the current direction of the LMC. Meanwhile, we find that GSE-associated subhalos do not significantly contribute to present-day Milky Way-est subhalo populations. These results provide context for our Galaxy's dark matter structure and subhalo population and will help interpret a range of measurements that are currently only possible in the MW. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 19 pages, 10 figures, 2 tables

SASHIMI-SIDM: Semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales

Authors: Shin'ichiro Ando, Shunichi Horigome, Ethan O. Nadler, Daneng Yang, Hai-Bo Yu

Abstract: We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that… ▽ More We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses; this result is uniquely enabled by our semi-analytical approach. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 21 pages, 7 figures. The code accompanying this submission is available at https://github.com/shinichiroando/sashimi-si

Tidal evolution of cored and cuspy dark matter halos

Authors: Xiaolong Du, Andrew Benson, Zhichao Carton Zeng, Tommaso Treu, Annika H. G. Peter, Charlie Mace, Fangzhou Jiang, Shengqi Yang, Charles Gannon, Daniel Gilman, Anna. M. Nierenberg, Ethan O. Nadler

Abstract: The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first cons… ▽ More The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first consider idealized N-body simulations of subhalos within their host halo, using a generalized mass density profile that describes their properties in a variety of dark matter models at infall, including popular warm, cold, and self-interacting ones. Using these simulations we construct tidal "tracks" for the evolution of subhalos based on their conditions at infall. Second, we use the results of these simulations to build semi-analytic models for tidal effects, including stripping and heating and implement them within the code GALACTICUS. Our semi-analytic models can accurately predict the tidal evolution of both cored and cuspy subhalos, including the bound mass and density profiles, providing a powerful and efficient tool for studying the post-infall properties of subhalos in different dark matter models. △ Less

Submitted 15 July, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: 25 pages, 50 figures. Version 2: matches published version. NOTE that the definition of dynamical timescale is slightly changed in version 2. The best-fit tidal stripping parameter is updated accordingly

Journal ref: Phys. Rev. D 110, 023019 (2024)

SAGAbg I: A Near-Unity Mass Loading Factor in Low-Mass Galaxies via their Low-Redshift Evolution in Stellar Mass, Oxygen Abundance, and Star Formation Rate

Authors: Erin Kado-Fong, Marla Geha, Yao-Yuan Mao, Mithi A. C. de los Reyes, Risa H. Wechsler, Yasmeen Asali, Nitya Kallivayalil, Ethan O. Nadler, Erik J. Tollerud, Benjamin Weiner

Abstract: Measuring the relation between star formation and galactic winds is observationally difficult. In this work we make an indirect measurement of the mass loading factor (the ratio between mass outflow rate and star formation rate) in low-mass galaxies using a differential approach to modeling the low-redshift evolution of the star-forming main sequence and mass-metallicity relation. We use the SAGA… ▽ More Measuring the relation between star formation and galactic winds is observationally difficult. In this work we make an indirect measurement of the mass loading factor (the ratio between mass outflow rate and star formation rate) in low-mass galaxies using a differential approach to modeling the low-redshift evolution of the star-forming main sequence and mass-metallicity relation. We use the SAGA (Satellites Around Galactic Analogs) background galaxies, those spectra observed by the SAGA survey that are not associated with the main SAGA host galaxies, to construct a sample of 11925 spectroscopically confirmed low-mass galaxies from $0.01\lesssim z \leq 0.21$ and measure a auroral line metallicity for 120 galaxies. The crux of the method is to use the lowest redshift galaxies as the boundary condition of our model, and to infer a mass-loading factor for the sample by comparing the expected evolution of the low redshift reference sample in stellar mass, gas-phase metallicity, and star formation rate against the observed properties of the sample at higher redshift. We infer a mass-loading factor of $η_{\rm m}=0.92^{+1.76}_{-0.74}$, which is in line with direct measurements of the mass-loading factor from the literature despite the drastically different set of assumptions needed for each approach. While our estimate of the mass-loading factor is in good agreement with recent galaxy simulations that focus on resolving the dynamics of the interstellar medium, it is smaller by over an order of magnitude than the mass-loading factor produced by many contemporary cosmological simulations. △ Less

Submitted 5 March, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 37 pages, 19 figures, accepted to ApJ

Forecasts for Galaxy Formation and Dark Matter Constraints from Dwarf Galaxy Surveys

Authors: Ethan O. Nadler, Vera Gluscevic, Trey Driskell, Risa H. Wechsler, Leonidas A. Moustakas, Andrew Benson, Yao-Yuan Mao

Abstract: The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)--mass cosmological zoom-in simulations and an empiric… ▽ More The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)--mass cosmological zoom-in simulations and an empirical galaxy--halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy--halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of $M_{50}=10^8~M_{\mathrm{\odot}}$ at the $1σ$ level; however, a tail toward low $M_{50}$ prevents a $2σ$ measurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on $M_{50}$ at the $1σ$ level, but the $2σ$ tail toward low $M_{50}$ persists. We project that observations of one (two) complete satellite populations can constrain warm DM models with $m_{\mathrm{WDM}}\approx 10~\mathrm{keV}$ ($20~\mathrm{keV}$). Subhalo mass function (SHMF) suppression can be constrained to $\approx 70\%$, $60\%$, and $50\%$ that in cold dark matter (CDM) at peak virial masses of $10^8$, $10^9$, and $10^{10}~M_{\mathrm{\odot}}$, respectively; SHMF enhancement constraints are weaker ($\approx 20$, $4$, and $2$ times that in CDM, respectively) due to galaxy--halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics. △ Less

Submitted 17 May, 2024; v1 submitted 18 January, 2024; originally announced January 2024.

Comments: 38 pages, 23 figures, 2 tables. Updated to published version

Journal ref: ApJ 967, 61 (2024)

Dynamical friction in self-interacting ultralight dark matter

Authors: Noah Glennon, Nathan Musoke, Ethan O. Nadler, Chanda Prescod-Weinstein, Risa H. Wechsler

Abstract: We explore how dynamical friction in an ultralight dark matter (ULDM) background is affected by dark matter self-interactions. We calculate the force of dynamical friction on a point mass moving through a uniform ULDM background with self-interactions, finding that the force of dynamical friction vanishes for sufficiently strong repulsive self-interactions. Using the pseudospectral solver… ▽ More We explore how dynamical friction in an ultralight dark matter (ULDM) background is affected by dark matter self-interactions. We calculate the force of dynamical friction on a point mass moving through a uniform ULDM background with self-interactions, finding that the force of dynamical friction vanishes for sufficiently strong repulsive self-interactions. Using the pseudospectral solver $\texttt{UltraDark.jl}$, we show with simulations that reasonable values of the ULDM self-interaction strength and particle mass cause $\mathcal{O}(1)$ differences in the acceleration of an object like a supermassive black hole (SMBH) traveling near the center of a soliton, relative to the case with no self-interactions. For example, repulsive self-interactions with $λ= 10^{-90}$ yield a deceleration due to dynamical friction $\approx70\%$ smaller than a model with no self-interactions. We discuss the observational implications of our results for SMBHs near soliton centers and for massive satellite galaxies falling into ultralight axion halos and show that outcomes are dependent on whether a self-interaction is present or not. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 17 pages; animations available at https://doi.org/10.5281/zenodo.7927475 or https://www.youtube.com/playlist?list=PLHrf0iQS5SY6COihRVYJkz31smUyDvFwW

MSC Class: 85A40

Journal ref: Phys.Rev.D 109 (2024) 6, 063501

A search for faint resolved galaxies beyond the Milky Way in DES Year 6: A new faint, diffuse dwarf satellite of NGC 55

Authors: M. McNanna, K. Bechtol, S. Mau, E. O. Nadler, J. Medoff, A. Drlica-Wagner, W. Cerny, D. Crnojevic, B. Mutlu-Pakdil, A. K. Vivas, A. B. Pace, J. L. Carlin, M. L. M. Collins, P. S. Ferguson, D. Martinez-Delgado, C. E. Martinez-Vazquez, N. E. D. Noel, A. H. Riley, D. J. Sand, A. Smercina, E. Tollerud, R. H. Wechsler, T. M. C. Abbott, M. Aguena, O. Alves , et al. (41 additional authors not shown)

Abstract: We report results from a systematic wide-area search for faint dwarf galaxies at heliocentric distances from 0.3 to 2 Mpc using the full six years of data from the Dark Energy Survey (DES). Unlike previous searches over the DES data, this search specifically targeted a field population of faint galaxies located beyond the Milky Way virial radius. We derive our detection efficiency for faint, resol… ▽ More We report results from a systematic wide-area search for faint dwarf galaxies at heliocentric distances from 0.3 to 2 Mpc using the full six years of data from the Dark Energy Survey (DES). Unlike previous searches over the DES data, this search specifically targeted a field population of faint galaxies located beyond the Milky Way virial radius. We derive our detection efficiency for faint, resolved dwarf galaxies in the Local Volume with a set of synthetic galaxies and expect our search to be complete to $M_V$ ~ $(-7, -10)$ mag for galaxies at $D = (0.3, 2.0)$ Mpc respectively. We find no new field dwarfs in the DES footprint, but we report the discovery of one high-significance candidate dwarf galaxy at a distance of $2.2\substack{+0.05\\-0.12}$ Mpc, a potential satellite of the Local Volume galaxy NGC 55, separated by $47$ arcmin (physical separation as small as 30 kpc). We estimate this dwarf galaxy to have an absolute V-band magnitude of $-8.0\substack{+0.5\\-0.3}$ mag and an azimuthally averaged physical half-light radius of $2.2\substack{+0.5\\-0.4}$ kpc, making this one of the lowest surface brightness galaxies ever found with $μ= 32.3$ mag ${\rm arcsec}^{-2}$. This is the largest, most diffuse galaxy known at this luminosity, suggesting possible tidal interactions with its host. △ Less

Submitted 4 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 20 pages, 7 figures

Report number: FERMILAB-PUB-23-478-PPD

A comprehensive model for the formation and evolution of the faintest Milky Way dwarf satellites

Authors: Niusha Ahvazi, Andrew Benson, Laura V. Sales, Ethan O. Nadler, Sachi Weerasooriya, Xiaolong Du, Mia Sauda Bovill

Abstract: In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the… ▽ More In this study, we modify the semi-analytic model Galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass-metallicity relation for dwarf galaxies requires us to include H$_2$ cooling, an updated UV background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size$-$mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy-halo connection. We predict a total of $300 ^{+75} _{-99}$ satellites with an absolute $V$-band magnitude (M$_{V}$) less than $0$ within $300$ kpc from our Milky Way-analogs. The fraction of subhalos that host a galaxy at least this bright drops to $50\%$ by a halo peak mass of $\sim 8.9 \times 10^{7}$ M$_{\odot}$, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population. △ Less

Submitted 18 March, 2024; v1 submitted 25 August, 2023; originally announced August 2023.

Comments: 20 pages, 11 figures. Accepted for publication in MNRAS

Symfind: Addressing the Fragility of Subhalo Finders and Revealing the Durability of Subhalos

Authors: Philip Mansfield, Elise Darragh-Ford, Yunchong Wang, Ethan O. Nadler, Risa H. Wechsler

Abstract: A major question in $Λ$CDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind, a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower… ▽ More A major question in $Λ$CDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind, a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower masses than commonly used halo-finding tools, with a focus on Rockstar and consistent-trees. These longer survival mean that at a fixed peak subhalo mass, we find $\approx 15\%{-}40\%$ more subhalos within the virial radius, $R_\textrm{vir}$, and $\approx 35\%-120\%$ more subhalos within $R_\textrm{vir}/4$ in the Symphony dark-matter-only simulation suite. More subhalos are found as resolution is increased. We perform extensive numerical testing. In agreement with idealized simulations, we show that the $v_{\rm max}$ of subhalos is only resolved at high resolutions ($n_\textrm{peak}\gtrsim3\times 10^4$), but that mass loss itself can be resolved at much more modest particle counts ($n_\textrm{peak}\gtrsim4\times 10^3$). We show that Rockstar converges to false solutions for the mass function, radial distribution, and disruption masses of subhalos. We argue that our new method can trace resolved subhalos until the point of typical galaxy disruption without invoking ``orphan'' modeling. We outline a concrete set of steps for determining whether other subhalo finders meet the same criteria. We publicly release Symfind catalogs and particle data for the Symphony simulation suite at \url{http://web.stanford.edu/group/gfc/symphony}. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 45 pages, 19 figures

A Self-interacting Dark Matter Solution to the Extreme Diversity of Low-mass Halo Properties

Authors: Ethan O. Nadler, Daneng Yang, Hai-Bo Yu

Abstract: The properties of low-mass dark matter (DM) halos appear to be remarkably diverse relative to cold, collisionless DM predictions, even in the presence of baryons. We show that self-interacting DM (SIDM) can simultaneously explain observations of halo diversity at two opposite extremes$-$the inner density profile of the dense substructure perturbing the strong lens galaxy SDSSJ0946+1006 and the rot… ▽ More The properties of low-mass dark matter (DM) halos appear to be remarkably diverse relative to cold, collisionless DM predictions, even in the presence of baryons. We show that self-interacting DM (SIDM) can simultaneously explain observations of halo diversity at two opposite extremes$-$the inner density profile of the dense substructure perturbing the strong lens galaxy SDSSJ0946+1006 and the rotation curves of isolated, gas-rich ultradiffuse galaxies (UDGs). To achieve this, we present the first cosmological zoom-in simulation featuring strong DM self-interactions in a galaxy group environment centered on a $10^{13}~M_{\mathrm{\odot}}$ host halo. In our SIDM simulation, most surviving subhalos of the group-mass host are deeply core-collapsed, yielding excellent candidates for the observed dense strong-lens perturber. Self-interactions simultaneously create kiloparsec-scale cores in low-concentration isolated halos, which could host the observed UDGs. Our scenario can be further tested with observations of DM structure and galaxies over a wide mass range. △ Less

Submitted 30 November, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 10 pages, 5 figures. Updated to published version

Journal ref: ApJL 958, L39 (2023)

A Parametric Model for Self-Interacting Dark Matter Halos

Authors: Daneng Yang, Ethan O. Nadler, Hai-Bo Yu, Yi-Ming Zhong

Abstract: We propose a parametric model for studying self-interacting dark matter (SIDM) halos. The model uses an analytical density profile, calibrated using a controlled N-body SIDM simulation that covers the entire gravothermal evolution, including core-forming and -collapsing phases. By normalizing the calibrated density profile, we obtain a universal description for SIDM halos at any evolution phase. T… ▽ More We propose a parametric model for studying self-interacting dark matter (SIDM) halos. The model uses an analytical density profile, calibrated using a controlled N-body SIDM simulation that covers the entire gravothermal evolution, including core-forming and -collapsing phases. By normalizing the calibrated density profile, we obtain a universal description for SIDM halos at any evolution phase. The model allows us to infer properties of SIDM halos based on their cold dark matter (CDM) counterparts. As a basic application, we only require two characteristic parameters of an isolated CDM halo at $z=0$. We then extend the model to incorporate effects induced by halo mass changes, such as major mergers or tidal stripping, making it applicable to both isolated halos and subhalos. The parametric model is tested and validated using cosmological zoom-in SIDM simulations available in the literature. △ Less

Submitted 21 March, 2024; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 22 pages, 14 figures, typos fixed, added example rotation curves and an appendix for the model based on the CORE-NFW ("Read") profile; example scripts provided through a link

Journal ref: JCAP 02 (2024) 032

Non-thermal warm dark matter limits from small-scale structure

Authors: Arka Banerjee, Subinoy Das, Anshuman Maharana, Ethan O. Nadler, Ravi Kumar Sharma

Abstract: We present small-scale structure constraints on sterile dark matter produced from a heavy mediator particle, inspired by models of moduli decay. Dark matter particles produced through this mechanism can contribute to the entire dark matter energy density but the particles have a non-thermal phase-space distribution; however, we show that the resulting linear matter power spectra can be mapped to e… ▽ More We present small-scale structure constraints on sterile dark matter produced from a heavy mediator particle, inspired by models of moduli decay. Dark matter particles produced through this mechanism can contribute to the entire dark matter energy density but the particles have a non-thermal phase-space distribution; however, we show that the resulting linear matter power spectra can be mapped to effective thermal-relic warm dark matter models. This production mechanism is therefore subject to warm dark matter constraints from small-scale structure as probed by ultra-faint dwarf galaxy abundances and strong gravitational lensing flux ratio statistics. We use the correspondence to thermal-relic models to derive a lower bound on the non-thermal particle mass of $107\ \mathrm{keV}$, at $95\%$ confidence. These are the first and most stringent constraints derived on sterile dark matter produced via the heavy mediator decay scenario we consider. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Comments are most welcome

Journal ref: Phys. Rev. D 108, 043518 (2023)

Can Neutrino Self-interactions Save Sterile Neutrino Dark Matter?

Authors: Rui An, Vera Gluscevic, Ethan O. Nadler, Yue Zhang

Abstract: Sterile neutrinos only interact with the Standard Model through the neutrino sector, and thus represent a simple dark matter (DM) candidate with many potential astrophysical and cosmological signatures. Recently, sterile neutrinos produced through self-interactions of active neutrinos have received attention as a particle candidate that can yield the entire observed DM relic abundance without viol… ▽ More Sterile neutrinos only interact with the Standard Model through the neutrino sector, and thus represent a simple dark matter (DM) candidate with many potential astrophysical and cosmological signatures. Recently, sterile neutrinos produced through self-interactions of active neutrinos have received attention as a particle candidate that can yield the entire observed DM relic abundance without violating the most stringent constraints from X-ray observations. We examine consistency of this production mechanism with the abundance of small-scale structure in the universe, as captured by the population of ultra-faint dwarf galaxies orbiting the Milky Way, and derive a lower bound on the sterile-neutrino particle mass of $37$ keV. Combining these results with previous limits from particle physics and astrophysics excludes $100\%$ sterile neutrino DM produced by strong neutrino self-coupling, mediated by a heavy ($\gtrsim 1~\mathrm{GeV}$) scalar particle; however, data permits sterile-neutrino DM production via a light mediator. △ Less

Submitted 15 November, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 9 pages, 5 figures, comments are welcome

Journal ref: ApJL 954, L18 (2023)

Growing the First Galaxies' Merger Trees

Authors: Ethan O. Nadler, Andrew Benson, Trey Driskell, Xiaolong Du, Vera Gluscevic

Abstract: Modelling the growth histories of specific galaxies often involves generating the entire population of objects that arise in a given cosmology and selecting systems with appropriate properties. This approach is highly inefficient when targeting rare systems such as the extremely luminous high-redshift galaxy candidates detected by JWST. Here, we present a novel framework for generating merger tree… ▽ More Modelling the growth histories of specific galaxies often involves generating the entire population of objects that arise in a given cosmology and selecting systems with appropriate properties. This approach is highly inefficient when targeting rare systems such as the extremely luminous high-redshift galaxy candidates detected by JWST. Here, we present a novel framework for generating merger trees with branches that are guaranteed to achieve a desired halo mass at a chosen redshift. This method augments extended Press Schechter theory solutions with constrained random processes known as Brownian bridges and is implemented in the open-source semi-analytic model $\texttt{Galacticus}$. We generate ensembles of constrained merger trees to predict the growth histories of seven high-redshift JWST galaxy candidates, finding that these systems most likely merge $\approx 2~\mathrm{Gyr}$ after the observation epoch and occupy haloes of mass $\gtrsim 10^{14}~M_{\mathrm{\odot}}$ today. These calculations are thousands of times more efficient than existing methods, are analytically controlled, and provide physical insights into the evolution of haloes with rapid early growth. Our constrained merger tree implementation is publicly available at http://github.com/galacticusorg/galacticus. △ Less

Submitted 23 March, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Comments: 20 pages, 12 figures. Updated to published version

Journal ref: MNRAS, 521, 3201 (2023)

Target Selection and Sample Characterization for the DESI LOW-Z Secondary Target Program

Authors: Elise Darragh-Ford, John F. Wu, Yao-Yuan Mao, Risa H. Wechsler, Marla Geha, Jaime E. Forero-Romero, ChangHoon Hahn, Nitya Kallivayalil, John Moustakas, Ethan O. Nadler, Marta Nowotka, J. E. G. Peek, Erik J. Tollerud, Benjamin Weiner, J. Aguilar, S. Ahlen, D. Brooks, A. P. Cooper, A. de la Macorra, A. Dey, K. Fanning, A. Font-Ribera, S. Gontcho A Gontcho, K. Honscheid, T. Kisner , et al. (17 additional authors not shown)

Abstract: We introduce the DESI LOW-Z Secondary Target Survey, which combines the wide-area capabilities of the Dark Energy Spectroscopic Instrument (DESI) with an efficient, low-redshift target selection method. Our selection consists of a set of color and surface brightness cuts, combined with modern machine learning methods, to target low-redshift dwarf galaxies ($z$ < 0.03) between $19 < r < 21$ with hi… ▽ More We introduce the DESI LOW-Z Secondary Target Survey, which combines the wide-area capabilities of the Dark Energy Spectroscopic Instrument (DESI) with an efficient, low-redshift target selection method. Our selection consists of a set of color and surface brightness cuts, combined with modern machine learning methods, to target low-redshift dwarf galaxies ($z$ < 0.03) between $19 < r < 21$ with high completeness. We employ a convolutional neural network (CNN) to select high-priority targets. The LOW-Z survey has already obtained over 22,000 redshifts of dwarf galaxies (M$_* < 10^9$ M$_\odot$), comparable to the number of dwarf galaxies discovered in SDSS-DR8 and GAMA. As a spare fiber survey, LOW-Z currently receives fiber allocation for just ~50% of its targets. However, we estimate that our selection is highly complete: for galaxies at $z < 0.03$ within our magnitude limits, we achieve better than 95% completeness with ~1% efficiency using catalog-level photometric cuts. We also demonstrate that our CNN selections $z<0.03$ galaxies from the photometric cuts subsample at least ten times more efficiently while maintaining high completeness. The full five-year DESI program will expand the LOW-Z sample, densely mapping the low-redshift Universe, providing an unprecedented sample of dwarf galaxies, and providing critical information about how to pursue effective and efficient low-redshift surveys. △ Less

Submitted 12 June, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 24 pages, 14 figures, data to reproduce figures: https://zenodo.org/record/7422591

Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations

Authors: Philip F. Hopkins, Ethan O. Nadler, Michael Y. Grudic, Xuejian Shen, Isabel Sands, Fangzhou Jiang

Abstract: Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simula… ▽ More Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multi-physics simulations where species with different softening laws may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated timestep criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously-discussed adaptive softening schemes based on the local neighbor particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the tidal softening scheme not only represents a physically-motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method, but imposes negligible timestep or other computational penalties, ensures that pairwise two-body scattering is small compared to smooth background forces, and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code. △ Less

Submitted 28 August, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: 25 pages, 13 figures. Updated to match MNRAS accepted version. Six appendices with additional details and tests expanded

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

Authors: Daneng Yang, Ethan O. Nadler, Hai-Bo Yu

Abstract: We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)--like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches $\approx 100~\mathrm{cm}^2\ \mathrm{g}^{-1}$ at relative velocities of $\approx 10~\mathrm{km\ s}^{-1}$, motivated by the diverse and orbitally dependent ce… ▽ More We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)--like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches $\approx 100~\mathrm{cm}^2\ \mathrm{g}^{-1}$ at relative velocities of $\approx 10~\mathrm{km\ s}^{-1}$, motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos' gravothermal evolution, demonstrating that deviations from the collisionless $R_{\rm max}$--$V_{\rm max}$ relation can be used to select deeply core-collapsed halos, where $V_{\rm max}$ is a halo's maximum circular velocity, and $R_{\rm max}$ is the radius at which it occurs. We predict that a sizable fraction ($\approx 20\%$) of subhalos with masses down to $\approx 10^8~M_{\odot}$ is deeply core collapsed in our SIDM model. Core-collapsed systems form $\approx 10\%$ of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above $\approx 10^9~M_{\odot}$ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies. △ Less

Submitted 2 June, 2023; v1 submitted 24 November, 2022; originally announced November 2022.

Comments: 29 pages, 17 figures

Journal ref: 2023 ApJ 949 67

Six More Ultra-Faint Milky Way Companions Discovered in the DECam Local Volume Exploration Survey

Authors: W. Cerny, C. E. Martínez-Vázquez, A. Drlica-Wagner, A. B. Pace, B. Mutlu-Pakdil, T. S. Li, A. H. Riley, D. Crnojević, C. R. Bom, J. A. Carballo-Bello, J. L. Carlin, A. Chiti, Y. Choi, M. L. M. Collins, E Darragh-Ford, P. S. Ferguson, M. Geha, D. Martínez-Delgado, P. Massana, S. Mau, G. E. Medina, R. R. Muñoz, E. O. Nadler, K. A. G. Olsen, A. Pieres , et al. (6 additional authors not shown)

Abstract: We report the discovery of six ultra-faint Milky Way satellites discovered through matched-filter searches conducted using Dark Energy Camera (DECam) data processed as part of the second data release of the DECam Local Volume Exploration (DELVE) survey. Leveraging deep Gemini/GMOS-N imaging (for four candidates) as well as follow-up DECam imaging (for two candidates), we characterize the morpholog… ▽ More We report the discovery of six ultra-faint Milky Way satellites discovered through matched-filter searches conducted using Dark Energy Camera (DECam) data processed as part of the second data release of the DECam Local Volume Exploration (DELVE) survey. Leveraging deep Gemini/GMOS-N imaging (for four candidates) as well as follow-up DECam imaging (for two candidates), we characterize the morphologies and stellar populations of these systems. We find that these candidates all share faint absolute magnitudes ($M_{V} \geq -3.2$ mag) and old, metal-poor stellar populations ($τ> 10$ Gyr, [Fe/H] $< -1.4$ dex). Three of these systems are more extended ($r_{1/2} > 15$ pc), while the other three are compact ($r_{1/2} < 10$ pc). From these properties, we infer that the former three systems (Boötes V, Leo Minor I, and Virgo II) are consistent with ultra-faint dwarf galaxy classifications, whereas the latter three (DELVE 3, DELVE 4, and DELVE 5) are likely ultra-faint star clusters. Using data from the Gaia satellite, we confidently measure the proper motion of Boötes V, Leo Minor I, and DELVE 4, and tentatively detect a proper motion signal from DELVE 3 and DELVE 5; no signal is detected for Virgo II. We use these measurements to explore possible associations between the newly-discovered systems and the Sagittarius dwarf spheroidal, the Magellanic Clouds, and the Vast Polar Structure, finding several plausible associations. Our results offer a preview of the numerous ultra-faint stellar systems that will soon be discovered by the Vera C. Rubin Observatory and highlight the challenges of classifying the faintest stellar systems. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 30 pages, 12 Figures (including Appendix). Submitted to ApJ. We encourage the reader to also review Smith et al. 2022, "Discovery of a new Local Group Dwarf Galaxy Candidate in UNIONS: Boötes V" (arxiv: 2209.08242), who independently present the discovery of one of the candidates reported here. We are working to make code and data products available

Report number: FERMILAB-PUB-22-704-LDRD-PPD

Report of the Topical Group on Cosmic Probes of Dark Matter for Snowmass 2021

Authors: Alex Drlica-Wagner, Chanda Prescod-Weinstein, Hai-Bo Yu, Andrea Albert, Mustafa Amin, Arka Banerjee, Masha Baryakhtar, Keith Bechtol, Simeon Bird, Simon Birrer, Torsten Bringmann, Regina Caputo, Sukanya Chakrabarti, Thomas Y. Chen, Djuna Croon, Francis-Yan Cyr-Racine, William A. Dawson, Cora Dvorkin, Vera Gluscevic, Daniel Gilman, Daniel Grin, Renée Hložek, Rebecca K. Leane, Ting S. Li, Yao-Yuan Mao , et al. (15 additional authors not shown)

Abstract: Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to… ▽ More Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to inform our understanding of the fundamental nature of dark matter in the coming decade. △ Less

Submitted 13 December, 2022; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: Report of the CF3 Topical Group for Snowmass 2021; 35 pages, 10 figures, many references. V3 updates Fig 3-2 and the author list

Structure Formation and the Global 21-cm Signal in the Presence of Coulomb-like Dark Matter-Baryon Interactions

Authors: Trey Driskell, Ethan O. Nadler, Jordan Mirocha, Andrew Benson, Kimberly K. Boddy, Timothy D. Morton, Jack Lashner, Rui An, Vera Gluscevic

Abstract: Many compelling dark matter (DM) scenarios feature Coulomb-like interactions between DM particles and baryons, in which the cross section for elastic scattering scales with relative particle velocity as $v^{-4}$. Previous studies have invoked such interactions to produce heat exchange between cold DM and baryons and alter the temperature evolution of hydrogen. In this study, we present a comprehen… ▽ More Many compelling dark matter (DM) scenarios feature Coulomb-like interactions between DM particles and baryons, in which the cross section for elastic scattering scales with relative particle velocity as $v^{-4}$. Previous studies have invoked such interactions to produce heat exchange between cold DM and baryons and alter the temperature evolution of hydrogen. In this study, we present a comprehensive study of the effects of Coulomb-like scattering on structure formation, in addition to the known effects on the thermal history of hydrogen. We find that interactions which significantly alter the temperature of hydrogen at Cosmic Dawn also dramatically suppress the formation of galaxies that source the Lyman-$α$ background, further affecting the global 21-cm signal. In particular, an interaction cross section at the current observational upper limit leads to a decrease in the abundance of star-forming halos by a factor of $\sim 2$ at $z\sim 20$, relative to cold, collisionless DM. We also find that DM that is 100% millicharged cannot reproduce the depth and the timing of the reported EDGES anomaly in any part of the parameter space. These results critically inform modeling of the global 21-cm signal and structure formation in cosmologies with DM-baryon scattering, with repercussions for future and upcoming cosmological data analysis. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: 21 pages, 10 figures, submitted to PRD

Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass

Authors: Ethan O. Nadler, Philip Mansfield, Yunchong Wang, Xiaolong Du, Susmita Adhikari, Arka Banerjee, Andrew Benson, Elise Darragh-Ford, Yao-Yuan Mao, Sebastian Wagner-Carena, Risa H. Wechsler, Hao-Yi Wu

Abstract: We present Symphony, a compilation of $262$ cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from $10^{11}$$-$$10^{15}~M_{\mathrm{\odot}}$. This compilation includes three existing simulation suites at the cluster and Milky Way$-$mass scales, and two new suites: $39$ Large Magellanic Cloud-mass ($10^{11}~M_{\mathrm{\odot}}$) and $49$ strong-lens-anal… ▽ More We present Symphony, a compilation of $262$ cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from $10^{11}$$-$$10^{15}~M_{\mathrm{\odot}}$. This compilation includes three existing simulation suites at the cluster and Milky Way$-$mass scales, and two new suites: $39$ Large Magellanic Cloud-mass ($10^{11}~M_{\mathrm{\odot}}$) and $49$ strong-lens-analog ($10^{13}~M_{\mathrm{\odot}}$) group-mass hosts. Across the entire host halo mass range, the highest-resolution regions in these simulations are resolved with a dark matter particle mass of $\approx 3\times 10^{-7}$ times the host virial mass and a Plummer-equivalent gravitational softening length of $\approx 9\times 10^{-4}$ times the host virial radius, on average. We measure correlations between subhalo abundance and host concentration, formation time, and maximum subhalo mass, all of which peak at the Milky Way host halo mass scale. Subhalo abundances are $\approx 50\%$ higher in clusters than in lower-mass hosts at fixed sub-to-host halo mass ratios. Subhalo radial distributions are approximately self-similar as a function of host mass and are less concentrated than hosts' underlying dark matter distributions. We compare our results to the semianalytic model $\mathrm{\texttt{Galacticus}}$, which predicts subhalo mass functions with a higher normalization at the low-mass end and radial distributions that are slightly more concentrated than Symphony. We use $\mathrm{\texttt{UniverseMachine}}$ to model halo and subhalo star formation histories in Symphony, and we demonstrate that these predictions resolve the formation histories of the halos that host nearly all currently observable satellite galaxies in the universe. To promote open use of Symphony, data products are publicly available at http://web.stanford.edu/group/gfc/symphony. △ Less

Submitted 16 March, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 39 pages, 20 figures, 2 tables. Updated to published version

Journal ref: ApJ 945, 159 (2023)

Tidal disruption of solitons in self-interacting ultralight axion dark matter

Authors: Noah Glennon, Ethan O. Nadler, Nathan Musoke, Arka Banerjee, Chanda Prescod-Weinstein, Risa H. Wechsler

Abstract: Ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal… ▽ More Ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal evolution of solitons is affected by ULA self-interaction strength and sign. Specifically, using the pseudospectral solver UltraDark.jl, we simulate the tidal disruption of self-interacting solitonic cores as they orbit a $10^{11}~M_{\mathrm{\odot}}$ Navarro-Frenk-White CDM host halo potential for a range of orbital parameters, assuming a fiducial ULA particle mass of $10^{-22}\mathrm{eV}$. We find that repulsive (attractive) self-interactions significantly accelerate (decelerate) soliton tidal disruption. We also identify a degeneracy between the self-interaction strength and soliton mass that determines the efficiency of tidal disruption, such that disruption timescales are affected at the $\sim 50\%$ level for variations in the dimensionless ULA self-coupling from $λ=-10^{-92}$ to $λ=10^{-92}$. △ Less

Submitted 6 July, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Comments: 15 pages, 10 figures. This paper was made to be as similar to the PRD version as possible

Gravothermal solutions of SIDM halos: mapping from constant to velocity-dependent cross section

Authors: Shengqi Yang, Xiaolong Du, Zhichao Carton Zeng, Andrew Benson, Fangzhou Jiang, Ethan O. Nadler, Annika H. G. Peter

Abstract: The scale-free gravothermal fluid formalism has long proved effective in describing the evolution of self-interacting dark matter halos with a constant dark matter particle cross section. However, whether the gravothermal fluid solutions match numerical simulations for velocity-dependent cross-section scenarios remains untested. In this work, we provide a fast mapping method that relates the const… ▽ More The scale-free gravothermal fluid formalism has long proved effective in describing the evolution of self-interacting dark matter halos with a constant dark matter particle cross section. However, whether the gravothermal fluid solutions match numerical simulations for velocity-dependent cross-section scenarios remains untested. In this work, we provide a fast mapping method that relates the constant-cross-section gravothermal solution to models with arbitrary velocity dependence in the cross section. We show that the gravothermal solutions after mapping are in good agreement with Arepo N-Body simulation results. We illustrate the power of this approach by applying this fast mapping method to a halo hosting a low surface brightness galaxy UGC 128. We show that this fast mapping method can be used to constrain free parameters in a physically motivated cross-section model and illustrate parameter space favored by the rotation curve measurement. △ Less

Submitted 2 March, 2023; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 20 pages, 11 figures

Snowmass2021 Cosmic Frontier White Paper: Dark Matter Physics from Halo Measurements

Authors: Keith Bechtol, Simon Birrer, Francis-Yan Cyr-Racine, Katelin Schutz, Susmita Adhikari, Mustafa Amin, Arka Banerjee, Simeon Bird, Nikita Blinov, Kimberly K. Boddy, Celine Boehm, Kevin Bundy, Malte Buschmann, Sukanya Chakrabarti, David Curtin, Liang Dai, Alex Drlica-Wagner, Cora Dvorkin, Adrienne L. Erickcek, Daniel Gilman, Saniya Heeba, Stacy Kim, Vid Iršič, Alexie Leauthaud, Mark Lovell , et al. (19 additional authors not shown)

Abstract: The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determin… ▽ More The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics. △ Less

Submitted 24 April, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: White paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). 88 pages, 9 figures. Comments welcome

Snowmass2021: Vera C. Rubin Observatory as a Flagship Dark Matter Experiment

Authors: Yao-Yuan Mao, Annika H. G. Peter, Susmita Adhikari, Keith Bechtol, Simeon Bird, Simon Birrer, Jonathan Blazek, Jeffrey L. Carlin, Nushkia Chamba, Johann Cohen-Tanugi, Francis-Yan Cyr-Racine, Tansu Daylan, Birendra Dhanasingham, Alex Drlica-Wagner, Cora Dvorkin, Christopher Fassnacht, Eric Gawiser, Maurizio Giannotti, Vera Gluscevic, Alma Gonzalez-Morales, Renee Hlozek, M. James Jee, Stacy Kim, Akhtar Mahmood, Rachel Mandelbaum , et al. (8 additional authors not shown)

Abstract: Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requir… ▽ More Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Snowmass2021 Cosmic Frontier White Paper: Cosmological Simulations for Dark Matter Physics

Authors: Arka Banerjee, Kimberly K. Boddy, Francis-Yan Cyr-Racine, Adrienne L. Erickcek, Daniel Gilman, Vera Gluscevic, Stacy Kim, Benjamin V. Lehmann, Yao-Yuan Mao, Philip Mocz, Ferah Munshi, Ethan O. Nadler, Lina Necib, Aditya Parikh, Annika H. G. Peter, Laura Sales, Mark Vogelsberger, Anna C. Wright

Abstract: Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CM… ▽ More Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics. △ Less

Submitted 21 March, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Astrophysical and Cosmological Probes of Dark Matter

Authors: Kimberly K. Boddy, Mariangela Lisanti, Samuel D. McDermott, Nicholas L. Rodd, Christoph Weniger, Yacine Ali-Haïmoud, Malte Buschmann, Ilias Cholis, Djuna Croon, Adrienne L. Erickcek, Vera Gluscevic, Rebecca K. Leane, Siddharth Mishra-Sharma, Julian B. Muñoz, Ethan O. Nadler, Priyamvada Natarajan, Adrian Price-Whelan, Simona Vegetti, Samuel J. Witte

Abstract: While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for t… ▽ More While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade. △ Less

Submitted 12 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021, 24 pages, 4 figures, comments welcome

From Images to Dark Matter: End-To-End Inference of Substructure From Hundreds of Strong Gravitational Lenses

Authors: Sebastian Wagner-Carena, Jelle Aalbers, Simon Birrer, Ethan O. Nadler, Elise Darragh-Ford, Philip J. Marshall, Risa H. Wechsler

Abstract: Constraining the distribution of small-scale structure in our universe allows us to probe alternatives to the cold dark matter paradigm. Strong gravitational lensing offers a unique window into small dark matter halos ($<10^{10} M_\odot$) because these halos impart a gravitational lensing signal even if they do not host luminous galaxies. We create large datasets of strong lensing images with real… ▽ More Constraining the distribution of small-scale structure in our universe allows us to probe alternatives to the cold dark matter paradigm. Strong gravitational lensing offers a unique window into small dark matter halos ($<10^{10} M_\odot$) because these halos impart a gravitational lensing signal even if they do not host luminous galaxies. We create large datasets of strong lensing images with realistic low-mass halos, Hubble Space Telescope (HST) observational effects, and galaxy light from HST's COSMOS field. Using a simulation-based inference pipeline, we train a neural posterior estimator of the subhalo mass function (SHMF) and place constraints on populations of lenses generated using a separate set of galaxy sources. We find that by combining our network with a hierarchical inference framework, we can both reliably infer the SHMF across a variety of configurations and scale efficiently to populations with hundreds of lenses. By conducting precise inference on large and complex simulated datasets, our method lays a foundation for extracting dark matter constraints from the next generation of wide-field optical imaging surveys. △ Less

Submitted 6 March, 2023; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: Code available at https://github.com/swagnercarena/paltas

Journal ref: ApJ 942 75 (2023)

R2-D2: Roman and Rubin -- From Data to Discovery

Authors: Suvi Gezari, Misty Bentz, Kishalay De, K. Decker French, Aaron Meisner, Michelle Ntampaka, Robert Jedicke, Ekta Patel, Daniel Perley, Robyn Sanderson, Christian Aganze, Igor Andreoni, Eric F. Bell, Edo Berger, Ian Dell'Antonio, Ryan Foley, Henry Hsieh, Mansi Kasliwal, Joel Kastner, Charles D. Kilpatrick, J. Davy Kirkpatrick, Casey Lam, Karen Meech, Dante Minniti, Ethan O. Nadler , et al. (6 additional authors not shown)

Abstract: The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to… ▽ More The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

Milky Way Satellite Census. IV. Constraints on Decaying Dark Matter from Observations of Milky Way Satellite Galaxies

Authors: S. Mau, E. O. Nadler, R. H. Wechsler, A. Drlica-Wagner, K. Bechtol, G. Green, D. Huterer, T. S. Li, Y. -Y. Mao, C. E. Martínez-Vázquez, M. McNanna, B. Mutlu-Pakdil, A. B. Pace, A. Peter, A. H. Riley, L. Strigari, M. -Y. Wang, M. Aguena, S. Allam, J. Annis, D. Bacon, E. Bertin, S. Bocquet, D. Brooks, D. L. Burke , et al. (44 additional authors not shown)

Abstract: We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime $τ$ comparable to the age of the Universe to a lighter DM particle (with mass splitting $ε$) and to a dark radiation species. These decays impart a characteristic "kick velo… ▽ More We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime $τ$ comparable to the age of the Universe to a lighter DM particle (with mass splitting $ε$) and to a dark radiation species. These decays impart a characteristic "kick velocity," $V_{\mathrm{kick}}=εc$, on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function of $τ$ and $V_{\mathrm{kick}}$ using a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy--halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models with $τ< 18\ \mathrm{Gyr}$ ($29\ \mathrm{Gyr}$) for $V_{\mathrm{kick}}=20\ \mathrm{km}\, \mathrm{s}^{-1}$ ($40\ \mathrm{km}\, \mathrm{s}^{-1}$) at $95\%$ confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble and $S_8$ tensions. △ Less

Submitted 27 June, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: 26 pages, 11 figures, 1 table. Updated to published version

Report number: FERMILAB-PUB-21-690-PPD

Journal ref: ApJ 932, 128 (2022)

Extending the SAGA Survey (xSAGA) I: Satellite Radial Profiles as a Function of Host Galaxy Properties

Authors: John F. Wu, J. E. G. Peek, Erik J. Tollerud, Yao-Yuan Mao, Ethan O. Nadler, Marla Geha, Risa H. Wechsler, Nitya Kallivayalil, Benjamin J. Weiner

Abstract: We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ gal… ▽ More We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ galaxies from more distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of $> 100,000$ CNN-selected low-$z$ galaxies, we identify $>20,000$ probable satellites located between 36-300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar mass range $9.5 < \log(M_\star/M_\odot) < 11$. We characterize the incompleteness and contamination for CNN-selected samples, and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36-300 kpc does not depend strongly on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations, and highlight the promise of using machine learning for extending galaxy samples of wide-area surveys. △ Less

Submitted 25 January, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: 27 pages, 18 figures, 1 table. Accepted to ApJ. Code available at https://github.com/jwuphysics/xsaga

Journal ref: ApJ (2022), 927, 121

The Effects of Dark Matter and Baryonic Physics on the Milky Way Subhalo Population in the Presence of the Large Magellanic Cloud

Authors: Ethan O. Nadler, Arka Banerjee, Susmita Adhikari, Yao-Yuan Mao, Risa H. Wechsler

Abstract: Given recent developments in our understanding of the Large Magellanic Cloud's (LMC) impact on the Milky Way's (MW) dark matter subhalo population, we compare the signatures of dark matter and baryonic physics on subhalos in MW systems with realistic LMC analogs. In particular, we study the effects of self-interacting dark matter (SIDM), warm dark matter (WDM), and the Galactic disk on the peak ma… ▽ More Given recent developments in our understanding of the Large Magellanic Cloud's (LMC) impact on the Milky Way's (MW) dark matter subhalo population, we compare the signatures of dark matter and baryonic physics on subhalos in MW systems with realistic LMC analogs. In particular, we study the effects of self-interacting dark matter (SIDM), warm dark matter (WDM), and the Galactic disk on the peak maximum circular velocity ($V_{\mathrm{peak}}$) function, radial distribution, and spatial distribution of MW and LMC-associated subhalos using cosmological dark matter-only zoom-in simulations of MW+LMC systems. For a fixed abundance of subhalos expected to host dwarf galaxies ($V_{\mathrm{peak}}\gtrsim 20\ \mathrm{km\ s}^{-1}$), SIDM and WDM can produce a similar mass-dependent suppression of the subhalo $V_{\mathrm{peak}}$ function, while disk disruption is mass independent. Subhalos in the inner regions of the MW are preferentially disrupted by both self-interactions and the disk, while suppression in WDM is radially independent. The relative abundance of LMC-associated subhalos is not strongly affected by disk disruption or WDM, but is significantly suppressed in SIDM due to self-interactions with the LMC at early times and with the MW during LMC infall at late times, erasing spatial anisotropy in the MW subhalo population. These results provide avenues to distinguish dark matter and baryonic physics by combining properties of the MW and LMC subhalo populations probed by upcoming observations of satellite galaxies and stellar streams. △ Less

Submitted 12 October, 2021; v1 submitted 24 September, 2021; originally announced September 2021.

Comments: 9 pages, 2 figures, 1 table. Updated to published version

Report number: FERMILAB-PUB-21-318-T

Journal ref: ApJL 920, 11 (2021)

Observational constraints on dark matter scattering with electrons

Authors: David Nguyen, Dimple Sarnaaik, Kimberly K. Boddy, Ethan O. Nadler, Vera Gluscevic

Abstract: We present new observational constraints on the elastic scattering of dark matter with electrons for dark matter masses between 10 keV and 1 TeV. We consider scenarios in which the momentum-transfer cross section has a power-law dependence on the relative particle velocity, with a power-law index $n \in \{-4,-2,0,2,4,6\}$. We search for evidence of dark matter scattering through its suppression of… ▽ More We present new observational constraints on the elastic scattering of dark matter with electrons for dark matter masses between 10 keV and 1 TeV. We consider scenarios in which the momentum-transfer cross section has a power-law dependence on the relative particle velocity, with a power-law index $n \in \{-4,-2,0,2,4,6\}$. We search for evidence of dark matter scattering through its suppression of structure formation. Measurements of the cosmic microwave background temperature, polarization, and lensing anisotropy from \textit{Planck} 2018 data and of the Milky Way satellite abundance measurements from the Dark Energy Survey and Pan-STARRS1 show no evidence of interactions. We use these data sets to obtain upper limits on the scattering cross section, comparing them with exclusion bounds from electronic recoil data in direct detection experiments. Our results provide the strongest bounds available for dark matter--electron scattering derived from the distribution of matter in the Universe, extending down to sub-MeV dark matter masses, where current direct detection experiments lose sensitivity. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: modified CLASS code available at https://github.com/kboddy/class_public/tree/dmeff

Report number: UTTG-04-2021

The Signatures of Self-Interacting Dark Matter and Subhalo Disruption on Cluster Substructure

Authors: Joy Bhattacharyya, Susmita Adhikari, Arka Banerjee, Surhud More, Amit Kumar, Ethan O. Nadler, Suchetana Chatterjee

Abstract: The abundance, distribution and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross-section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simula… ▽ More The abundance, distribution and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross-section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simulations are suppressed relative to their cold dark matter (CDM) counterparts. Once the population of disrupted subhalos -- which may host orphan galaxies -- are taken into account, satellite galaxy populations in CDM and SIDM models can be reconciled. However, even in this case, the inner structure of subhalos are significantly different in the two dark matter models. We study the feasibility of using the weak lensing signal from the subhalo density profiles to distinguish between the cold and self-interacting dark matter while accounting for the potential contribution of orphan galaxies. We find that the effects of self-interactions on the density profile of subhalos can appear degenerate with subhalo disruption in CDM, when orphans are accounted for. With current error bars from the Subaru Hyper Suprime-Cam Strategic Program, we find that subhalos in the outskirts of clusters (where disruption is less prevalent) can be used to constrain dark matter physics. In the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time will give precise measurements of the weak lensing profile and can be used to constrain $σ_T/m$ at the $\sim 1$ cm$^2$ g$^{-1}$ level at $v\sim 2000$ km s$^{-1}$. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 19 pages, 12 figures, to be submitted to ApJ

Journal ref: ApJ 932 (2022) 30

The DECam Local Volume Exploration Survey: Overview and First Data Release

Authors: A. Drlica-Wagner, J. L. Carlin, D. L. Nidever, P. S. Ferguson, N. Kuropatkin, M. Adamów, W. Cerny, Y. Choi, J. H. Esteves, C. E. Martínez-Vázquez, S. Mau, A. E. Miller, B. Mutlu-Pakdil, E. H. Neilsen, K. A. G. Olsen, A. B. Pace, A. H. Riley, J. D. Sakowska, D. J. Sand, L. Santana-Silva, E. J. Tollerud, D. L. Tucker, A. K. Vivas, E. Zaborowski, A. Zenteno , et al. (45 additional authors not shown)

Abstract: The DECam Local Volume Exploration survey (DELVE) is a 126-night survey program on the 4-m Blanco Telescope at the Cerro Tololo Inter-American Observatory in Chile. DELVE seeks to understand the characteristics of faint satellite galaxies and other resolved stellar substructures over a range of environments in the Local Volume. DELVE will combine new DECam observations with archival DECam data to… ▽ More The DECam Local Volume Exploration survey (DELVE) is a 126-night survey program on the 4-m Blanco Telescope at the Cerro Tololo Inter-American Observatory in Chile. DELVE seeks to understand the characteristics of faint satellite galaxies and other resolved stellar substructures over a range of environments in the Local Volume. DELVE will combine new DECam observations with archival DECam data to cover ~15000 deg$^2$ of high-Galactic-latitude (|b| > 10 deg) southern sky to a 5$σ$ depth of g,r,i,z ~ 23.5 mag. In addition, DELVE will cover a region of ~2200 deg$^2$ around the Magellanic Clouds to a depth of g,r,i ~ 24.5 mag and an area of ~135 deg$^2$ around four Magellanic analogs to a depth of g,i ~ 25.5 mag. Here, we present an overview of the DELVE program and progress to date. We also summarize the first DELVE public data release (DELVE DR1), which provides point-source and automatic aperture photometry for ~520 million astronomical sources covering ~5000 deg$^2$ of the southern sky to a 5$σ$ point-source depth of g=24.3, r=23.9, i=23.3, and z=22.8 mag. DELVE DR1 is publicly available via the NOIRLab Astro Data Lab science platform. △ Less

Submitted 2 September, 2021; v1 submitted 12 March, 2021; originally announced March 2021.

Comments: 30 pages, 14 figures, 6 tables; updated to match published version; public data release at https://datalab.noirlab.edu/delve

Report number: FERMILAB-PUB-21-075-AE-LDRD

Journal ref: ApJS 256, 2 (2021)

UniverseMachine: Predicting Galaxy Star Formation over Seven Decades of Halo Mass with Zoom-in Simulations

Authors: Yunchong Wang, Ethan O. Nadler, Yao-Yuan Mao, Susmita Adhikari, Risa H. Wechsler, Peter Behroozi

Abstract: We apply the empirical galaxy--halo connection model UniverseMachine to dark matter-only zoom-in simulations of isolated Milky Way (MW)--mass halos along with their parent cosmological simulations. This application extends \textsc{UniverseMachine} predictions into the ultra-faint dwarf galaxy regime ($ 10^{2}\,\mathrm{M_{\odot}} \leqslant M_{\ast} \leqslant 10^{5}\,\mathrm{M_{\odot}}$) and yields… ▽ More We apply the empirical galaxy--halo connection model UniverseMachine to dark matter-only zoom-in simulations of isolated Milky Way (MW)--mass halos along with their parent cosmological simulations. This application extends \textsc{UniverseMachine} predictions into the ultra-faint dwarf galaxy regime ($ 10^{2}\,\mathrm{M_{\odot}} \leqslant M_{\ast} \leqslant 10^{5}\,\mathrm{M_{\odot}}$) and yields a well-resolved stellar mass--halo mass (SMHM) relation over the peak halo mass range $10^8\,\mathrm{M_{\odot}}$ to $10^{15}\,\mathrm{M_{\odot}}$. The extensive dynamic range provided by the zoom-in simulations allows us to assess specific aspects of dwarf galaxy evolution predicted by \textsc{UniverseMachine}. In particular, although UniverseMachine is not constrained for dwarf galaxies with $M_* \lesssim 10^{8}\,\mathrm{M_{\odot}}$, our predicted SMHM relation is consistent with that inferred for MW satellite galaxies at $z=0$ using abundance matching. However, UniverseMachine predicts that nearly all galaxies are actively star forming below $M_{\ast}\sim 10^{7}\,\mathrm{M_{\odot}}$ and that these systems typically form more than half of their stars at $z\lesssim 4$, which is discrepant with the star formation histories of Local Group dwarf galaxies that favor early quenching. This indicates that the current UniverseMachine model does not fully capture galaxy quenching physics at the low-mass end. We highlight specific improvements necessary to incorporate environmental and reionization-driven quenching for dwarf galaxies, and provide a new tool to connect dark matter accretion to star formation over the full dynamic range that hosts galaxies. △ Less

Submitted 2 June, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: Accepted for publication in ApJ. 21 pages, 8 figures, 1 table. Minor edits to v1. Added discussion of reionization contribution of ultra-faint dwarfs in sections 4.4 and 5

Journal ref: ApJ, 915, 116 (2021)

Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies

Authors: Ethan O. Nadler, Simon Birrer, Daniel Gilman, Risa H. Wechsler, Xiaolong Du, Andrew Benson, Anna M. Nierenberg, Tommaso Treu

Abstract: Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and… ▽ More Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semianalytic model $\mathrm{\texttt{Galacticus}}$ to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ M_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ M_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts. △ Less

Submitted 9 August, 2021; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: 23 pages, 9 figures, 1 table. Updated to published version

Journal ref: ApJ 917, 7 (2021)

Searching for Dwarf Galaxies in ${\it Gaia}$ DR2 Phase-Space Data Using Wavelet Transforms

Authors: Elise Darragh-Ford, Ethan O. Nadler, Sean McLaughlin, Risa H. Wechsler

Abstract: We present a wavelet-based algorithm to identify dwarf galaxies in the Milky Way in ${\it Gaia}$ DR2 data. Our algorithm detects overdensities in 4D position--proper motion space, making it the first search to explicitly use velocity information to search for dwarf galaxy candidates. We optimize our algorithm and quantify its performance by searching for mock dwarfs injected into ${\it Gaia}$ DR2… ▽ More We present a wavelet-based algorithm to identify dwarf galaxies in the Milky Way in ${\it Gaia}$ DR2 data. Our algorithm detects overdensities in 4D position--proper motion space, making it the first search to explicitly use velocity information to search for dwarf galaxy candidates. We optimize our algorithm and quantify its performance by searching for mock dwarfs injected into ${\it Gaia}$ DR2 data and for known Milky Way satellite galaxies. Comparing our results with previous photometric searches, we find that our search is sensitive to undiscovered systems at Galactic latitudes~$\lvert b\rvert>20^{\circ}$ and with half-light radii larger than the 50% detection efficiency threshold for Pan-STARRS1 (PS1) at (${\it i}$) absolute magnitudes of =$-7<M_V<-3$ and distances of $32$ kpc $< D < 64$ kpc, and (${\it ii}$) $M_V< -4$ and $64$ kpc $< D < 128$ kpc. Based on these results, we predict that our search is expected to discover $5 \pm 2$ new satellite galaxies: four in the PS1 footprint and one outside the Dark Energy Survey and PS1 footprints. We apply our algorithm to the ${\it Gaia}$ DR2 dataset and recover $\sim 830$ high-significance candidates, out of which we identify a "gold standard" list of $\sim 200$ candidates based on cross-matching with potential candidates identified in a preliminary search using ${\it Gaia}$ EDR3 data. All of our candidate lists are publicly distributed for future follow-up studies. We show that improvements in astrometric measurements provided by ${\it Gaia}$ EDR3 increase the sensitivity of this technique; we plan to continue to refine our candidate list using future data releases. △ Less

Submitted 14 July, 2021; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 17 pages, 10 figures, 3 tables; submitted to ApJ. Candidate lists available at https://dwarfswaves.github.io. Updated to published version

Report number: SLAC-PUB-17614

Journal ref: ApJ 915 48 (2021)

Bounds on velocity-dependent dark matter-proton scattering from Milky Way satellite abundance

Authors: Karime Maamari, Vera Gluscevic, Kimberly K. Boddy, Ethan O. Nadler, Risa H. Wechsler

Abstract: We use the latest measurements of the Milky Way satellite population from the Dark Energy Survey and Pan-STARRS1 to infer the most stringent astrophysical bound to date on velocity-dependent interactions between dark matter particles and protons. We model the momentum-transfer cross section as a power law of the relative particle velocity $v$ with a free normalizing amplitude,… ▽ More We use the latest measurements of the Milky Way satellite population from the Dark Energy Survey and Pan-STARRS1 to infer the most stringent astrophysical bound to date on velocity-dependent interactions between dark matter particles and protons. We model the momentum-transfer cross section as a power law of the relative particle velocity $v$ with a free normalizing amplitude, $σ_\text{MT}=σ_0 v^n$, to broadly capture the interactions arising within the non-relativistic effective theory of dark matter-proton scattering. The scattering leads to a momentum and heat transfer between the baryon and dark matter fluids in the early Universe, ultimately erasing structure on small physical scales and reducing the abundance of low-mass halos that host dwarf galaxies today. From the consistency of observations with the cold collisionless dark matter paradigm, using a new method that relies on the most robust predictions of the linear perturbation theory, we infer an upper limit on $σ_0$ of $1.4\times10^{-23}$, $2.1\times10^{-19}$, and $1.0\times10^{-12}\ \mathrm{cm}^2$, for interaction models with $n=2,4,6$, respectively, for a dark matter particle mass of $10\ \mathrm{MeV}$. These results improve observational limits on dark matter--proton scattering by orders of magnitude and thus provide an important guide for viable sub-GeV dark matter candidates. △ Less

Submitted 8 March, 2022; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 8 pages, 3 figures, 1 table. Updated to correct typo in Table 1. Associated code available at https://github.com/maamari/VelocityDependentLimits

Report number: UTTG-13-2020

Journal ref: ApJL 907, 46 (2021)

Constraints on the epoch of dark matter formation from Milky Way satellites

Authors: Subinoy Das, Ethan O. Nadler

Abstract: A small fraction of thermalized dark radiation that transitions into cold dark matter (CDM) between big bang nucleosynthesis and matter-radiation equality can account for the entire dark matter relic density. Because of its transition from dark radiation, "late-forming dark matter" (LFDM) suppresses the growth of linear matter perturbations and imprints the oscillatory signatures of dark radiation… ▽ More A small fraction of thermalized dark radiation that transitions into cold dark matter (CDM) between big bang nucleosynthesis and matter-radiation equality can account for the entire dark matter relic density. Because of its transition from dark radiation, "late-forming dark matter" (LFDM) suppresses the growth of linear matter perturbations and imprints the oscillatory signatures of dark radiation perturbations on small scales. The cutoff scale in the linear matter power spectrum is set by the redshift $z_T$ of the phase transition; tracers of small-scale structure can therefore be used to infer the LFDM formation epoch. Here, we use a forward model of the Milky Way (MW) satellite galaxy population to address the question: How late can dark matter form? For dark radiation with strong self-interactions, which arises in theories of neutrinolike LFDM, we report $z_{T}>5.5\times 10^6$ at $95\%$ confidence based on the abundance of known MW satellite galaxies. This limit rigorously accounts for observational incompleteness corrections, marginalizes over uncertainties in the connection between dwarf galaxies and dark matter halos, and improves upon galaxy clustering and Lyman-$α$ forest constraints by nearly an order of magnitude. We show that this limit can also be interpreted as a lower bound on $z_T$ for LFDM that free-streams prior to its phase transition, although dedicated simulations will be needed to analyze this case in detail. Thus, dark matter created by a transition from dark radiation must form no later than one week after the big bang. △ Less

Submitted 20 January, 2021; v1 submitted 2 October, 2020; originally announced October 2020.

Comments: 10 pages, 3 figures. Updated to published version

Journal ref: Phys. Rev. D 103, 043517 (2021)

The SAGA Survey. II. Building a Statistical Sample of Satellite Systems around Milky Way-like Galaxies

Authors: Yao-Yuan Mao, Marla Geha, Risa H. Wechsler, Benjamin Weiner, Erik J. Tollerud, Ethan O. Nadler, Nitya Kallivayalil

Abstract: We present the Stage II results from the ongoing Satellites Around Galactic Analogs (SAGA) Survey. Upon completion, the SAGA Survey will spectroscopically identify satellite galaxies brighter than $ M_{r,o} = -12.3 $ around 100 Milky Way (MW) analogs at $ z \sim 0.01 $. In Stage II, we have more than quadrupled the sample size of Stage I, delivering results from 127 satellites around 36 MW analogs… ▽ More We present the Stage II results from the ongoing Satellites Around Galactic Analogs (SAGA) Survey. Upon completion, the SAGA Survey will spectroscopically identify satellite galaxies brighter than $ M_{r,o} = -12.3 $ around 100 Milky Way (MW) analogs at $ z \sim 0.01 $. In Stage II, we have more than quadrupled the sample size of Stage I, delivering results from 127 satellites around 36 MW analogs with an improved target selection strategy and deep photometric imaging catalogs from the Dark Energy Survey and the Legacy Surveys. We have obtained 25,372 galaxy redshifts, peaking around $ z = 0.2 $. These data significantly increase spectroscopic coverage for very low redshift objects in $ 17 < r_o < 20.75 $ around SAGA hosts, creating a unique data set that places the Local Group in a wider context. The number of confirmed satellites per system ranges from zero to nine, and correlates with host galaxy and brightest satellite luminosities. We find that the number and the luminosities of MW satellites are consistent with being drawn from the same underlying distribution as SAGA systems. The majority of confirmed SAGA satellites are star forming, and the quenched fraction increases as satellite stellar mass and projected radius from the host galaxy decrease. Overall, the satellite quenched fraction among SAGA systems is lower than that in the Local Group. We compare the luminosity functions and radial distributions of SAGA satellites with theoretical predictions based on cold dark matter simulations and an empirical galaxy-halo connection model and find that the results are broadly in agreement. △ Less

Submitted 17 December, 2020; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: 42 pages, 22 figures, 3 tables. See Sec. 9 (p.28) for summary and figure index. Main results are shown in Sec. 6-8 (Fig. 8-18). Accepted by ApJ. Data available on survey website: https://sagasurvey.org

Journal ref: Astrophysical Journal 907 (2021) 85

Milky Way Satellite Census. III. Constraints on Dark Matter Properties from Observations of Milky Way Satellite Galaxies

Authors: E. O. Nadler, A. Drlica-Wagner, K. Bechtol, S. Mau, R. H. Wechsler, V. Gluscevic, K. Boddy, A. B. Pace, T. S. Li, M. McNanna, A. H. Riley, J. García-Bellido, Y. -Y. Mao, G. Green, D. L. Burke, A. Peter, B. Jain, T. M. C. Abbott, M. Aguena, S. Allam, J. Annis, S. Avila, D. Brooks, M. Carrasco Kind, J. Carretero , et al. (43 additional authors not shown)

Abstract: We perform a comprehensive study of Milky Way (MW) satellite galaxies to constrain the fundamental properties of dark matter (DM). This analysis fully incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and marginalizes over uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW d… ▽ More We perform a comprehensive study of Milky Way (MW) satellite galaxies to constrain the fundamental properties of dark matter (DM). This analysis fully incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and marginalizes over uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk. Our results are consistent with the cold, collisionless DM paradigm and yield the strongest cosmological constraints to date on particle models of warm, interacting, and fuzzy dark matter. At $95\%$ confidence, we report limits on (i) the mass of thermal relic warm DM, $m_{\rm WDM} > 6.5\ \mathrm{keV}$ (free-streaming length, $λ_{\rm{fs}} \lesssim 10\,h^{-1}\ \mathrm{kpc}$), (ii) the velocity-independent DM-proton scattering cross section, $σ_{0} < 8.8\times 10^{-29}\ \mathrm{cm}^{2}$ for a $100\ \mathrm{MeV}$ DM particle mass (DM-proton coupling, $c_p \lesssim (0.3\ \mathrm{GeV})^{-2}$), and (iii) the mass of fuzzy DM, $m_φ> 2.9 \times 10^{-21}\ \mathrm{eV}$ (de Broglie wavelength, $λ_{\rm{dB}} \lesssim 0.5\ \mathrm{kpc}$). These constraints are complementary to other observational and laboratory constraints on DM properties. △ Less

Submitted 1 March, 2021; v1 submitted 31 July, 2020; originally announced August 2020.

Comments: 15 pages, 5 figures, 2 tables. Updated to published version

Report number: FERMILAB-PUB-20-277-AE

Journal ref: Phys. Rev. Lett. 126, 091101 (2021)

Signatures of Velocity-Dependent Dark Matter Self-Interactions in Milky Way-mass Halos

Authors: Ethan O. Nadler, Arka Banerjee, Susmita Adhikari, Yao-Yuan Mao, Risa H. Wechsler

Abstract: We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-in… ▽ More We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of $w$ that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of $w$. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of $\sim 200\ \rm{km\ s}^{-1}$, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models. △ Less

Submitted 18 June, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

Comments: 19 pages, 10 figures, 1 table. Updated to published version

Journal ref: ApJ 896, 112 (2020)

Milky Way Satellite Census. II. Galaxy--Halo Connection Constraints Including the Impact of the Large Magellanic Cloud

Authors: E. O. Nadler, R. H. Wechsler, K. Bechtol, Y. -Y. Mao, G. Green, A. Drlica-Wagner, M. McNanna, S. Mau, A. B. Pace, J. D. Simon, A. Kravtsov, S. Dodelson, T. S. Li, A. H. Riley, M. Y. Wang, T. M. C. Abbott, M. Aguena, S. Allam, J. Annis, S. Avila, G. M. Bernstein, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell , et al. (41 additional authors not shown)

Abstract: The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, w… ▽ More The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches. △ Less

Submitted 16 April, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 26 pages, 10 figures, 3 tables. Updated to published version

Report number: FERMILAB-PUB-19-606-AE

Journal ref: ApJ 893, 48 (2020)

Milky Way Satellite Census. I. The Observational Selection Function for Milky Way Satellites in DES Y3 and Pan-STARRS DR1

Authors: A. Drlica-Wagner, K. Bechtol, S. Mau, M. McNanna, E. O. Nadler, A. B. Pace, T. S. Li, A. Pieres, E. Rozo, J. D. Simon, A. R. Walker, R. H. Wechsler, T. M. C. Abbott, S. Allam, J. Annis, E. Bertin, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, M. Costanzi, L. N. da Costa, J. De Vicente, S. Desai , et al. (33 additional authors not shown)

Abstract: We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ~80% of the sky. Our search for satellite galaxies targets ~25,000 deg$^2$ of the high-Galactic-latitude sky reaching a 10$σ$ point-source depth of… ▽ More We report the results of a systematic search for ultra-faint Milky Way satellite galaxies using data from the Dark Energy Survey (DES) and Pan-STARRS1 (PS1). Together, DES and PS1 provide multi-band photometry in optical/near-infrared wavelengths over ~80% of the sky. Our search for satellite galaxies targets ~25,000 deg$^2$ of the high-Galactic-latitude sky reaching a 10$σ$ point-source depth of $\gtrsim$ 22.5 mag in the $g$ and $r$ bands. While satellite galaxy searches have been performed independently on DES and PS1 before, this is the first time that a self-consistent search is performed across both data sets. We do not detect any new high-significance satellite galaxy candidates, while recovering the majority of satellites previously detected in surveys of comparable depth. We characterize the sensitivity of our search using a large set of simulated satellites injected into the survey data. We use these simulations to derive both analytic and machine-learning models that accurately predict the detectability of Milky Way satellites as a function of their distance, size, luminosity, and location on the sky. To demonstrate the utility of this observational selection function, we calculate the luminosity function of Milky Way satellite galaxies, assuming that the known population of satellite galaxies is representative of the underlying distribution. We provide access to our observational selection function to facilitate comparisons with cosmological models of galaxy formation and evolution. △ Less

Submitted 16 April, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Comments: 34 pages, 11 figures, 6 tables; updated to match published version. Selection functions available at: https://github.com/des-science/mw-sats

Report number: FERMILAB-PUB-19-604-AE

Journal ref: ApJ 893, 1 (2020)