-
Unraveling the role of merger histories in the population of Insitu stars: linking IllustrisTNG cosmological simulation to H3 survey
Authors:
Razieh Emami,
Lars Hernquist,
Randall Smith,
James F. Steiner,
Grant Tremblay,
Douglas Finkbeiner,
Mark Vogelsberger,
Josh Grindlay,
Federico Marinacci,
Kung-Yi Su,
Cecilia Garraffo,
Yuan-Sen Ting,
Phillip A. Cargile,
Rebecca L. Davies,
Chloë E. Benton,
Yijia Li,
Letizia Bugiani,
Amir H. Khoram,
Sownak Bose
Abstract:
We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to…
▽ More
We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to define insitu stars. To uncover the underlying factors contributing to deviations between TNG50 and H3 distributions, we scrutinize correlation coefficients among internal drivers(e.g., virial radius, star formation rate [SFR]) and merger-related parameters (such as the effective mass-ratio, mean distance, average redshift, total number of mergers, average spin-ratio and maximum spin alignment between merging galaxies). Notably, we identify significant correlations between deviations from observational data and key parameters such as the median slope of virial radius, mean SFR values, and the rate of SFR change across different redshift scans. Furthermore, positive correlations emerge between deviations from observational data and parameters related to galaxy mergers. We validate these correlations using the Random Forest Regression method. Our findings underscore the invaluable insights provided by the H3 survey in unraveling the cosmic history of galaxies akin to the Milky Way, thereby advancing our understanding of galactic evolution and shedding light on the formation and evolution of Milky Way-like galaxies in cosmological simulations.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Brightest Cluster Galaxy Offsets in Cold Dark Matter
Authors:
Cian Roche,
Michael McDonald,
Josh Borrow,
Mark Vogelsberger,
Xuejian Shen,
Volker Springel,
Lars Hernquist,
Ruediger Pakmor,
Sownak Bose,
Rahul Kannan
Abstract:
The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the univers…
▽ More
The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above $10^{14}\rm{M_\odot}$, we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of $\sim 1\rm{kpc}$ in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of $1.48\rm{kpc}$. We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the "light traces mass" approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of $\sim10$ and $\sim30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section.
△ Less
Submitted 26 April, 2024; v1 submitted 1 February, 2024;
originally announced February 2024.
-
Varying primordial state fractions in exo- and endothermic SIDM simulations of Milky Way-mass haloes
Authors:
Aidan Leonard,
Stephanie O'Neil,
Xuejian Shen,
Mark Vogelsberger,
Olivia Rosenstein,
Hoatian Shangguan,
Yuanhong Teng,
Jiayi Hu
Abstract:
Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering occur with equal probability at high relative velocities (…
▽ More
Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering occur with equal probability at high relative velocities ($v_{\rm rel}\gtrsim400~{\rm km/s})$. In a suite of cosmological zoom-in simulation of Milky Way-size haloes, we vary the primordial state fractions to understand the impact of inelastic dark matter self-interactions on halo structure and evolution. In particular, we test how the initial conditions impact the present-day properties of dark matter haloes. Depending on the primordial state fraction, scattering reactions will be dominated by either exothermic or endothermic effects for high and low initial excited state fractions respectively. We find that increasing the initial excited fraction reduces the mass of the main halo, as well as the number of subhaloes on all mass scales. The main haloes are cored, with lower inner densities and higher outer densities compared with CDM. Additionally, we find that the shape of the main halo becomes more spherical the higher the initial excited state fraction is. Finally, we show that the number of satellites steadily decreases with initial excited state fraction across all satellite masses.
△ Less
Submitted 28 May, 2024; v1 submitted 24 January, 2024;
originally announced January 2024.
-
The EB-correlation in Resolved Polarized Images: Connections to Astrophysics of Black Holes
Authors:
Razieh Emami,
Sheperd S. Doeleman,
Maciek Wielgus,
Dominic Chang,
Koushik Chatterjee,
Randall Smith,
Matthew Liska,
James F. Steiner,
Angelo Ricarte,
Ramesh Narayan,
Grant Tremblay,
Douglas Finkbeiner,
Lars Hernquist,
Chi-Kwan Chan,
Lindy Blackburn,
Ben S. Prather,
Paul Tiede,
Avery E. Broderick,
Mark Vogelsberger,
Charles Alcock,
Freek Roelofs
Abstract:
We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the…
▽ More
We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the time-averaged EB-correlation contains novel information that might be linked to the BH spin, accretion state and the electron temperature. A detailed comparison with a semi-analytic approach with different azimuthal expansion modes shows that to recover the morphology of the real/imaginary part of the correlation function and its phase, we require higher orders of these azimuthal modes. To extract the phase features, we propose to use the Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of different geometrical ring models with various magnetic and velocity field morphologies and show that both the image space and visibility based EB-correlation morphologies in MAD simulations can be explained with simple fluid and magnetic field geometries as used in ring models. SANEs by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday Rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of EB-correlation might be well explained by the current models for BH spins as well as electron temperatures, while others may require a larger theoretical surveys.
△ Less
Submitted 3 May, 2023; v1 submitted 30 April, 2023;
originally announced May 2023.
-
THESAN-HR: Galaxies in the Epoch of Reionization in warm dark matter, fuzzy dark matter and interacting dark matter
Authors:
Xuejian Shen,
Josh Borrow,
Mark Vogelsberger,
Enrico Garaldi,
Aaron Smith,
Rahul Kannan,
Sandro Tacchella,
Jesús Zavala,
Lars Hernquist,
Jessica Y. -C. Yeh,
Chunyuan Zheng
Abstract:
Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (THESAN-HR), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matt…
▽ More
Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (THESAN-HR), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the UV luminosity functions at $z\gtrsim 6$ are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to "catch up" with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope (JWST) imaging surveys of deep, lensed fields hold potential for uncovering the faint, low-mass galaxy population, which could provide constraints on altDM models.
△ Less
Submitted 13 April, 2023;
originally announced April 2023.
-
Probing plasma composition with the next generation Event Horizon Telescope (ngEHT)
Authors:
Razieh Emami,
Richard Anantua,
Angelo Ricarte,
Sheperd S. Doeleman,
Avery Broderick,
George Wong,
Lindy Blackburn,
Maciek Wielgus,
Ramesh Narayan,
Grant Tremblay,
Charles Alcock,
Lars Hernquist,
Randall Smith,
Matthew Liska,
Priyamvada Natarajan,
Mark Vogelsberger,
Brandon Curd,
Joana A. Kramer
Abstract:
We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population…
▽ More
We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population of non-thermal electrons as well as a R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations get enhanced per adding pairs. On the contrary, we show that at higher frequencies a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
Tracing the hot spot motion using the next generation Event Horizon Telescope (ngEHT)
Authors:
Razieh Emami,
Paul Tiede,
Sheperd S. Doeleman,
Freek Roelofs,
Maciek Wielgus,
Lindy Blackburn,
Matthew Liska,
Koushik Chatterjee,
Bart Ripperda,
Antonio Fuentes,
Avery Broderick,
Lars Hernquist,
Charles Alcock,
Ramesh Narayan,
Randall Smith,
Grant Tremblay,
Angelo Ricarte,
He Sun,
Richard Anantua,
Yuri Y. Kovalev,
Priyamvada Natarajan,
Mark Vogelsberger
Abstract:
We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its bo…
▽ More
We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its boundary during its orbit, resulting in a distinct signature. We subdivide the motion to two distinct phases; the first phase refers to the appearance of the hot spot modelled as a bright blob, followed by a subsequent shearing phase simulated as a stretched ellipse. We employ different observational arrays, including EHT(2017,2022) and the next generation event horizon telescope (ngEHTp1, ngEHT) arrays, in which few new additional sites are added to the observational array. We make dynamical image reconstructions for each of these arrays. Subsequently, we infer the hot spot phase in the first phase followed by the axes ratio and the ellipse area in the second phase. We focus on the direct observability of the orbiting hot spot in the sub-mm wavelength. Our analysis demonstrates that newly added dishes may easily trace the first phase as well as part of the second phase, before the flux is reduced substantially. The algorithm used in this work can be extended to any other types of the dynamical motion. Consequently, we conclude that the ngEHT is a key to directly observe the dynamical motions near variable sources, such as SgrA*.
△ Less
Submitted 15 November, 2022; v1 submitted 12 November, 2022;
originally announced November 2022.
-
Kinematic signatures of impulsive supernova feedback in dwarf galaxies
Authors:
Jan D. Burger,
Jesús Zavala,
Laura V. Sales,
Mark Vogelsberger,
Federico Marinacci,
Paul Torrey
Abstract:
Impulsive supernova feedback and non-standard dark matter models, such as self-interacting dark matter (SIDM), are the two main contenders for the role of the dominant core formation mechanism at the dwarf galaxy scale. Here we show that the impulsive supernova cycles that follow episodes of bursty star formation leave distinct features in the distribution function of stars: groups of stars with s…
▽ More
Impulsive supernova feedback and non-standard dark matter models, such as self-interacting dark matter (SIDM), are the two main contenders for the role of the dominant core formation mechanism at the dwarf galaxy scale. Here we show that the impulsive supernova cycles that follow episodes of bursty star formation leave distinct features in the distribution function of stars: groups of stars with similar ages and metallicities develop overdense shells in phase space. If cores are formed through supernova feedback, we predict the presence of such features in star-forming dwarf galaxies with cored host halos. Their systematic absence would favor alternative dark matter models, such as SIDM, as the dominant core formation mechanism.
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
Endothermic self-interacting dark matter in Milky Way-like dark matter haloes
Authors:
Stephanie O'Neil,
Mark Vogelsberger,
Saniya Heeba,
Katelin Schutz,
Jonah C. Rose,
Paul Torrey,
Josh Borrow,
Ryan Low,
Rakshak Adhikari,
Mikhail V. Medvedev,
Tracy R. Slatyer,
Jesús Zavala
Abstract:
Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state,…
▽ More
Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state, which allows for both elastic and inelastic scattering. In particular, the model includes a significant probability for particles to up-scatter from the ground state to the excited state. We simulate a suite of zoom-in Milky Way-sized N-body haloes with six models with different scattering cross sections to study the effects of up-scattering in SIDM models. We find that the up-scattering reaction greatly increases the central densities of the main halo through the loss of kinetic energy. However, the physical model still results in significant coring due to the presence of elastic scattering and down-scattering. These effects are not as apparent in the subhalo population compared to the main halo, but the number of subhaloes is reduced compared to CDM.
△ Less
Submitted 28 May, 2024; v1 submitted 28 October, 2022;
originally announced October 2022.
-
Snowmass Theory Frontier: Astrophysics and Cosmology
Authors:
Daniel Green,
Joshua T. Ruderman,
Benjamin R. Safdi,
Jessie Shelton,
Ana Achúcarro,
Peter Adshead,
Yashar Akrami,
Masha Baryakhtar,
Daniel Baumann,
Asher Berlin,
Nikita Blinov,
Kimberly K. Boddy,
Malte Buschmann,
Giovanni Cabass,
Robert Caldwell,
Emanuele Castorina,
Thomas Y. Chen,
Xingang Chen,
William Coulton,
Djuna Croon,
Yanou Cui,
David Curtin,
Francis-Yan Cyr-Racine,
Christopher Dessert,
Keith R. Dienes
, et al. (62 additional authors not shown)
Abstract:
We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process.
We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process.
△ Less
Submitted 14 September, 2022;
originally announced September 2022.
-
An Effective Bias Expansion for 21 cm Cosmology in Redshift Space
Authors:
Wenzer Qin,
Katelin Schutz,
Aaron Smith,
Enrico Garaldi,
Rahul Kannan,
Tracy R. Slatyer,
Mark Vogelsberger
Abstract:
A near-future detection of the 21cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21cm brightness temperature field during the epoch of reionization, incorporating renormali…
▽ More
A near-future detection of the 21cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21cm brightness temperature field during the epoch of reionization, incorporating renormalized bias and a treatment of redshift space distortions. Notably, we confirm that in redshift space, measures of the 21cm brightness, e.g the power spectrum, should have irreducible contributions that lack a bias coefficient and therefore contain direct, astrophysics-free information about the cosmological density field; in this work, we study this effect beyond linear order. To validate our theoretical treatment, we fit the predicted EFT Fourier-space shapes to the THESAN suite of hydrodynamical simulations of reionization at the field level, where the considerable number of modes prevents overfitting. We find agreement at the level of a few percent between the 21cm power spectrum from the EFT fits and simulations over the wavenumber range $k \lesssim 0.8$ h/Mpc and neutral fraction $x_\mathrm{HI} \gtrsim 0.4$, which is imminently measurable by the Hydrogen Epoch of Reionization Array (HERA) and future experiments. The ability of the EFT to describe the 21cm signal extends to simulations that have different astrophysical prescriptions for reionization as well as simulations with interacting dark matter.
△ Less
Submitted 28 November, 2022; v1 submitted 12 May, 2022;
originally announced May 2022.
-
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.
-
X-ray morphology of cluster-mass haloes in self-interacting dark matter
Authors:
Xuejian Shen,
Thejs Brinckmann,
David Rapetti,
Mark Vogelsberger,
Adam Mantz,
Jesús Zavala,
Steven W. Allen
Abstract:
We perform cosmological zoom-in simulations of $19$ relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq 1\operatorname{-}3\times 10^{15}\,{\rm M}_{\odot}$. Both the dark matter and the intracluster ga…
▽ More
We perform cosmological zoom-in simulations of $19$ relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with $M_{\rm 200} \simeq 1\operatorname{-}3\times 10^{15}\,{\rm M}_{\odot}$. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of $84$ relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of $r\simeq 0.1\operatorname{-}0.2\,R_{\rm 200}$. The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as a systematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-core clusters agrees better with SIDM models with cross-section $(σ/m)= 0.5\operatorname{-}1~{\rm cm}^2/{\rm g}$ than CDM and SIDM with $(σ/m)=0.1\,{\rm cm}^2/{\rm g}$. Our statistical analysis indicates that the latter two models are disfavored at the $68\%$ confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.
△ Less
Submitted 1 November, 2022; v1 submitted 31 January, 2022;
originally announced February 2022.
-
Core-collapse, evaporation and tidal effects: the life story of a self-interacting dark matter subhalo
Authors:
Zhichao Carton Zeng,
Annika H. G. Peter,
Xiaolong Du,
Andrew Benson,
Stacy Kim,
Fangzhou Jiang,
Francis-Yan Cyr-Racine,
Mark Vogelsberger
Abstract:
Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, follow…
▽ More
Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of $\lesssim 10^8 M_\odot$ subhalos in lens-mass systems ($\sim 10^{13}M_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhalos with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhalos ($\lesssim 1/1000$ host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhalos to core-collapse, finding that it is nearly impossible to drive core-collapse in subhalos in SIDM models with constant cross sections. Any discovery of ultra-compact dark substructures with future substructure lensing observations favors additional degrees of freedom, such as velocity-dependence, in the cross section.
△ Less
Submitted 22 August, 2022; v1 submitted 1 October, 2021;
originally announced October 2021.
-
Inferring the Morphology of Stellar Distribution in TNG50: Twisted and Twisted-Stretched shapes
Authors:
Razieh Emami,
Lars Hernquist,
Charles Alcock,
Shy Genel,
Sownak Bose,
Rainer Weinberger,
Mark Vogelsberger,
Xuejian Shen,
Joshua S. Speagle,
Federico Marinacci,
John C. Forbes,
Paul Torrey
Abstract:
We investigate the morphology of the stellar distribution in a sample of Milky Way (MW) like galaxies in the TNG50 simulation. Using a local in shell iterative method (LSIM) as the main approach, we explicitly show evidence of twisting (in about 52% of halos) and stretching (in 48% of them) in the real space. This is matched with the re-orientation observed in the eigenvectors of the inertia tenso…
▽ More
We investigate the morphology of the stellar distribution in a sample of Milky Way (MW) like galaxies in the TNG50 simulation. Using a local in shell iterative method (LSIM) as the main approach, we explicitly show evidence of twisting (in about 52% of halos) and stretching (in 48% of them) in the real space. This is matched with the re-orientation observed in the eigenvectors of the inertia tensor and gives us a clear picture of having a re-oriented stellar distribution. We make a comparison between the shape profile of dark matter (DM) halo and stellar distribution and quite remarkably see that their radial profiles are fairly close, especially at small galactocentric radii where the stellar disk is located. This implies that the DM halo is somewhat aligned with stars in response to the baryonic potential. The level of alignment mostly decreases away from the center. We study the impact of substructures in the orbital circularity parameter. It is demonstrated that in some cases, far away substructures are counter-rotating compared with the central stars and may flip the sign of total angular momentum and thus the orbital circularity parameter. Truncating them above 150 kpc, however, retains the disky structure of the galaxy as per initial selection. Including the impact of substructures in the shape of stars, we explicitly show that their contribution is subdominant. Overlaying our theoretical results to the observational constraints from previous literature, we establish fair agreement.
△ Less
Submitted 5 June, 2021; v1 submitted 22 December, 2020;
originally announced December 2020.
-
ETHOS -- An Effective Theory of Structure Formation: Impact of Dark Acoustic Oscillations on Cosmic Dawn
Authors:
Julian B. Muñoz,
Sebastian Bohr,
Francis-Yan Cyr-Racine,
Jesús Zavala,
Mark Vogelsberger
Abstract:
Upcoming data of the 21-cm hydrogen line during cosmic dawn ($z\sim 10-30$) will revolutionize our understanding of the astrophysics of the first galaxies. Here we present a case study on how to exploit those same measurements to learn about the nature of dark matter (DM) at small scales. Focusing on the Effective Theory of Structure Formation (ETHOS) paradigm, we run a suite of simulations coveri…
▽ More
Upcoming data of the 21-cm hydrogen line during cosmic dawn ($z\sim 10-30$) will revolutionize our understanding of the astrophysics of the first galaxies. Here we present a case study on how to exploit those same measurements to learn about the nature of dark matter (DM) at small scales. Focusing on the Effective Theory of Structure Formation (ETHOS) paradigm, we run a suite of simulations covering a broad range of DM microphysics, connecting the output of $N$-body simulations to dedicated 21-cm simulations to predict the evolution of the 21-cm signal across the entire cosmic dawn. We find that observatories targeting both the global signal and the 21-cm power spectrum are sensitive to all ETHOS models we study, and can distinguish them from CDM if the suppression wavenumber is smaller than $k\approx 300\, h/$Mpc, even when accounting for feedback with a phenomenological model. This is an order of magnitude smaller comoving scales than currently constrained by other data sets, including the Lyman-$α$ forest. Moreover, if a prospective 21-cm detection confirmed a deficiency of power at small scales, we show that ETHOS models with strong dark acoustic oscillations can be discriminated from the pure suppression of warm dark matter, showing the power of 21-cm data to understand the behavior of DM at the smallest physical scales.
△ Less
Submitted 10 November, 2020;
originally announced November 2020.
-
DM halo morphological types of MW-like galaxies in the TNG50 simulation: Simple, Twisted, or Stretched
Authors:
Razieh Emami,
Shy Genel,
Lars Hernquist,
Charles Alcock,
Sownak Bose,
Rainer Weinberger,
Mark Vogelsberger,
Federico Marinacci,
Abraham Loeb,
Paul Torrey,
John C. Forbes
Abstract:
We present a comprehensive analysis of the shape of dark matter (DM) halos in a sample of 25 Milky Way-like galaxies in TNG50 simulation. Using an Enclosed Volume Iterative Method (EVIM), we infer an oblate-to-triaxial shape for the DM halo with the median $T \simeq 0.24 $. We group DM halos in 3 different categories. Simple halos (32% of population) establish principal axes whose ordering in magn…
▽ More
We present a comprehensive analysis of the shape of dark matter (DM) halos in a sample of 25 Milky Way-like galaxies in TNG50 simulation. Using an Enclosed Volume Iterative Method (EVIM), we infer an oblate-to-triaxial shape for the DM halo with the median $T \simeq 0.24 $. We group DM halos in 3 different categories. Simple halos (32% of population) establish principal axes whose ordering in magnitude does not change with radius and whose orientations are almost fixed throughout the halo. Twisted halos (32% of population), experience levels of gradual rotations throughout their radial profiles. Finally, stretched halos (36% of population) demonstrate a stretching in their principal axes lengths where the ordering of different eigenvalues change with radius. Subsequently, the halo experiences a "rotation" of $\sim$90 deg where the stretching occurs. Visualizing the 3D ellipsoid of each halo, for the first time, we report signs of re-orienting ellipsoid in twisted and stretched halos. We examine the impact of baryonic physics on DM halo shape through a comparison to dark matter only (DMO) simulations. This suggests a triaxial (prolate) halo. We analyze the impact of substructure on DM halo shape in both hydro and DMO simulations and confirm that their impacts are subdominant. We study the distribution of satellites in our sample. In simple and twisted halos, the angle of satellites' angular momentum with galaxy's angular momentum grows with radius. However, stretched halos show a flat distribution of angles. Overlaying our theoretical outcome on the observational results presented in the literature establishes a fair agreement.
△ Less
Submitted 23 March, 2021; v1 submitted 19 September, 2020;
originally announced September 2020.
-
Simulating the "hidden giant" in cold and self-interacting dark matter models
Authors:
Omid Sameie,
Sukanya Chakrabarti,
Hai-Bo Yu,
Michael Boylan-Kolchin,
Mark Vogelsberger,
Jesus Zavala,
Lars Hernquist
Abstract:
We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial o…
▽ More
We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way's tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.
△ Less
Submitted 11 June, 2020;
originally announced June 2020.
-
ETHOS -- An effective parametrization and classification for structure formation: the non-linear regime at $z\gtrsim5$
Authors:
Sebastian Bohr,
Jesús Zavala,
Francis-Yan Cyr-Racine,
Mark Vogelsberger,
Torsten Bringmann,
Christoph Pfrommer
Abstract:
We propose two effective parameters that fully characterise galactic-scale structure formation at high redshifts ($z\gtrsim5$) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal Warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and…
▽ More
We propose two effective parameters that fully characterise galactic-scale structure formation at high redshifts ($z\gtrsim5$) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal Warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and explore this parameter space, we use high-redshift zoom-in simulations that cover a wide range of non-linear scales from those where DM should behave as CDM ($k\sim10\,h\,{\rm Mpc}^{-1}$), down to those characterised by the onset of galaxy formation ($k\sim500\,h\,{\rm Mpc}^{-1}$). We show that the two physically motivated parameters $h_{\rm peak}$ and $k_{\rm peak}$, the amplitude and scale of the first DAO peak, respectively, are sufficient to parametrize the linear matter power spectrum and classify the DM models as belonging to effective non-linear structure formation regions. These are defined by their relative departure from Cold DM ($k_{\rm peak}\rightarrow\infty$) and WDM ($h_{\rm peak}=0$) according to the non-linear matter power spectrum and halo mass function. We identify a region where the DAOs still leave a distinct signature from WDM down to $z=5$, while a large part of the DAO parameter space is shown to be degenerate with WDM. Our framework can then be used to seamlessly connect a broad class of particle DM models to their structure formation properties at high redshift without the need of additional $N$-body simulations.
△ Less
Submitted 7 September, 2020; v1 submitted 2 June, 2020;
originally announced June 2020.
-
The diverse dark matter density at sub-kiloparsec scales in Milky Way satellites:implications for the nature of dark matter
Authors:
Jesús Zavala,
Mark R. Lovell,
Mark Vogelsberger,
Jan D. Burger
Abstract:
Milky Way (MW) satellites reside within dark matter (DM) subhalos with a broad distribution of circular velocity profiles. This diversity is enhanced with the inclusion of ultra-faint satellites, which seemingly have very high DM densities, albeit with large systematic uncertainties. We argue that if confirmed, this large diversity in the MW satellite population poses a serious test for the struct…
▽ More
Milky Way (MW) satellites reside within dark matter (DM) subhalos with a broad distribution of circular velocity profiles. This diversity is enhanced with the inclusion of ultra-faint satellites, which seemingly have very high DM densities, albeit with large systematic uncertainties. We argue that if confirmed, this large diversity in the MW satellite population poses a serious test for the structure formation theory with possible implications for the DM nature. For the Cold Dark Matter model, the diversity might be a signature of the combined effects of subhalo tidal disruption by the MW disk and strong supernova feedback. For models with a dwarf-scale cutoff in the power spectrum, the diversity is a consequence of the lower abundance of dwarf-scale halos. This diversity is most challenging for Self-Interacting Dark Matter (SIDM) models with cross sections $σ/m_χ\gtrsim1~$cm$^2$g$^{-1}$ where subhalos have too low densities to explain the ultra-faint galaxies. We propose a novel solution to explain the diversity of MW satellites based on the gravothermal collapse of SIDM haloes. This solution requires a velocity-dependent cross section that predicts a bimodal distribution of cuspy dense (collapsed) subhaloes consistent with the ultra-faint satellites, and cored lower density subhaloes consistent with the brighter satellites.
△ Less
Submitted 30 August, 2019; v1 submitted 22 April, 2019;
originally announced April 2019.
-
Self-Interacting Dark Matter Subhalos in the Milky Way's Tides
Authors:
Omid Sameie,
Hai-Bo Yu,
Laura V. Sales,
Mark Vogelsberger,
Jesus Zavala
Abstract:
We study evolution of self-interacting dark matter (SIDM) subhalos in the Milky Way (MW) tidal field. The interaction between the subhalos and the MW's tides lead to more diverse dark matter distribution in the inner region, compared to their cold dark matter counterparts. We test this scenario with two MW satellite galaxies, Draco and Fornax, opposite extremes in the inner dark matter content, an…
▽ More
We study evolution of self-interacting dark matter (SIDM) subhalos in the Milky Way (MW) tidal field. The interaction between the subhalos and the MW's tides lead to more diverse dark matter distribution in the inner region, compared to their cold dark matter counterparts. We test this scenario with two MW satellite galaxies, Draco and Fornax, opposite extremes in the inner dark matter content, and find that they can be accommodated within the SIDM model proposed to explain the diverse rotation curves of spiral galaxies in the field.
△ Less
Submitted 6 April, 2020; v1 submitted 16 April, 2019;
originally announced April 2019.
-
ETHOS - an effective theory of structure formation: formation of the first haloes and their stars
Authors:
Mark R. Lovell,
Jesús Zavala,
Mark Vogelsberger
Abstract:
A cutoff in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes and their galaxies at high and low redshift. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cutoff, and that has been shown to be an alternative to the cold dark matter (CDM) mo…
▽ More
A cutoff in the linear matter power spectrum at dwarf galaxy scales has been shown to affect the abundance, formation mechanism and age of dwarf haloes and their galaxies at high and low redshift. We use hydrodynamical simulations of galaxy formation within the ETHOS framework in a benchmark model that has such a cutoff, and that has been shown to be an alternative to the cold dark matter (CDM) model that alleviates its dwarf-scale challenges. We show how galaxies in this model form differently to CDM on a halo-by-halo basis, at redshifts $z\ge6$. We show that ETHOS haloes at the half-mode mass scale form with 50~per~cent less mass than their CDM counterparts due to their later formation times, yet they retain more of their gas reservoir due to the different behaviour of gas and dark matter during the monolithic collapse of the first haloes in models with a galactic-scale cutoff. As a result, galaxies in ETHOS haloes near the cutoff scale grow rapidly between $z=10-6$ and by $z=6$ end up having very similar stellar masses, higher gas fractions and higher star formation rates relative to their CDM counterparts. We highlight these differences by making predictions for how the number of galaxies with old stellar populations is suppressed in ETHOS for both $z=6$ galaxies and for gas-poor Local Group fossil galaxies. Interestingly, we find an age gradient in ETHOS between galaxies that form in high and low density environments.
△ Less
Submitted 11 December, 2018;
originally announced December 2018.
-
Gravitational Lensing and the Power Spectrum of Dark Matter Substructure: Insights from the ETHOS N-body Simulations
Authors:
Ana Díaz Rivero,
Cora Dvorkin,
Francis-Yan Cyr-Racine,
Jesús Zavala,
Mark Vogelsberger
Abstract:
Strong gravitational lensing has been identified as a promising astrophysical probe to study the particle nature of dark matter. In this paper we present a detailed study of the power spectrum of the projected mass density (convergence) field of substructure in a Milky Way-sized halo. This power spectrum has been suggested as a key observable that can be extracted from strongly lensed images and y…
▽ More
Strong gravitational lensing has been identified as a promising astrophysical probe to study the particle nature of dark matter. In this paper we present a detailed study of the power spectrum of the projected mass density (convergence) field of substructure in a Milky Way-sized halo. This power spectrum has been suggested as a key observable that can be extracted from strongly lensed images and yield important clues about the matter distribution within the lens galaxy. We use two different $N$-body simulations from the ETHOS framework: one with cold dark matter and another with self-interacting dark matter and a cutoff in the initial power spectrum. Despite earlier works that identified $ k \gtrsim 100$ kpc$^{-1}$ as the most promising scales to learn about the particle nature of dark matter we find that even at lower wavenumbers - which are actually within reach of observations in the near future - we can gain important information about dark matter. Comparing the amplitude and slope of the power spectrum on scales $0.1 \lesssim k/$kpc$^{-1} \lesssim 10$ from lenses at different redshifts can help us distinguish between cold dark matter and other exotic dark matter scenarios that alter the abundance and central densities of subhalos. Furthermore, by considering the contribution of different mass bins to the power spectrum we find that subhalos in the mass range $10^7 - 10^8$ M$_{\odot}$ are on average the largest contributors to the power spectrum signal on scales $2 \lesssim k/$kpc$^{-1} \lesssim 15$, despite the numerous subhalos with masses $> 10^8$ M$_{\odot}$ in a typical lens galaxy. Finally, by comparing the power spectra obtained from the subhalo catalogs to those from the particle data in the simulation snapshots we find that the seemingly-too-simple halo model is in fact a fairly good approximation to the much more complex array of substructure in the lens.
△ Less
Submitted 31 August, 2018;
originally announced September 2018.
-
Towards an improved model of self-interacting dark matter haloes
Authors:
Anastasia Sokolenko,
Kyrylo Bondarenko,
Thejs Brinckmann,
Jesús Zavala,
Mark Vogelsberger,
Torsten Bringmann,
Alexey Boyarsky
Abstract:
In this work, we discuss the relation between the strength of the self-interaction of dark matter particles and the predicted properties of the inner density distributions of dark matter haloes. We present the results of $N$-body simulations for 28 haloes performed with the same initial conditions for cold dark matter and for self-interacting dark matter, for a range of cross-sections. We provide…
▽ More
In this work, we discuss the relation between the strength of the self-interaction of dark matter particles and the predicted properties of the inner density distributions of dark matter haloes. We present the results of $N$-body simulations for 28 haloes performed with the same initial conditions for cold dark matter and for self-interacting dark matter, for a range of cross-sections. We provide a simple phenomenological description of these results and compare with the semi-analytical model typically used in the literature. Using these results, we then predict how the inner dark matter surface density and core radius should depend on the self-interaction cross-section for observed haloes.
△ Less
Submitted 30 January, 2019; v1 submitted 29 June, 2018;
originally announced June 2018.
-
Evaporating the Milky Way halo and its satellites with inelastic self-interacting dark matter
Authors:
Mark Vogelsberger,
Jesús Zavala,
Katelin Schutz,
Tracy R. Slatyer
Abstract:
Self-interacting dark matter provides a promising alternative for the cold dark matter paradigm to solve potential small-scale galaxy formation problems. Nearly all self-interacting dark matter simulations so far have considered only elastic collisions. Here we present simulations of a galactic halo within a generic inelastic model using a novel numerical implementation in the Arepo code to study…
▽ More
Self-interacting dark matter provides a promising alternative for the cold dark matter paradigm to solve potential small-scale galaxy formation problems. Nearly all self-interacting dark matter simulations so far have considered only elastic collisions. Here we present simulations of a galactic halo within a generic inelastic model using a novel numerical implementation in the Arepo code to study arbitrary multi-state inelastic dark matter scenarios. For this model we find that inelastic self-interactions can: (i) create larger subhalo density cores compared to elastic models for the same cross section normalisation; (ii) lower the abundance of satellites without the need for a power spectrum cutoff; (iii) reduce the total halo mass by about 10%; (iv) inject the energy equivalent of O(100) million Type II supernovae in galactic haloes through level de-excitation; (v) avoid the gravothermal catastrophe due to removal of particles from halo centers. We conclude that a ~5 times larger elastic cross section is required to achieve the same central density reduction as the inelastic model. This implies that well-established constraints on self-interacting cross sections have to be revised if inelastic collisions are the dominant mode. In this case significantly smaller cross sections can achieve the same core density reduction thereby increasing the parameter space of allowed models considerably.
△ Less
Submitted 16 March, 2019; v1 submitted 8 May, 2018;
originally announced May 2018.
-
ETHOS -- an effective theory of structure formation: predictions for the high-redshift Universe -- abundance of galaxies and reionization
Authors:
Mark R. Lovell,
Jesús Zavala,
Mark Vogelsberger,
Xuejian Shen,
Francis-Yan Cyr-Racine,
Christoph Pfrommer,
Kris Sigurdson,
Michael Boylan-Kolchin,
Annalisa Pillepich
Abstract:
We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cel…
▽ More
We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $τ$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.
△ Less
Submitted 9 January, 2019; v1 submitted 28 November, 2017;
originally announced November 2017.
-
ETHOS - An Effective Theory of Structure Formation: From dark particle physics to the matter distribution of the Universe
Authors:
Francis-Yan Cyr-Racine,
Kris Sigurdson,
Jesus Zavala,
Torsten Bringmann,
Mark Vogelsberger,
Christoph Pfrommer
Abstract:
We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross s…
▽ More
We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.
△ Less
Submitted 11 September, 2017; v1 submitted 16 December, 2015;
originally announced December 2015.
-
Scattering, Damping, and Acoustic Oscillations: Simulating the Structure of Dark Matter Halos with Relativistic Force Carriers
Authors:
Matthew R. Buckley,
Jesús Zavala,
Francis-Yan Cyr-Racine,
Kris Sigurdson,
Mark Vogelsberger
Abstract:
We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter scenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, t…
▽ More
We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter scenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, the self-interacting models have a much wider range of possible phenomenology. A long-range force in the dark matter can introduce multiple scales to the initial power spectrum, in the form of dark acoustic oscillations and an exponential cut-off in the power spectrum. Using simulations we show that the impact of these scales can remain observationally relevant up to the present day. Furthermore, the self-interaction can continue to modify the small-scale structure of the dark matter halos, reducing their central densities and creating a dark matter core. The resulting phenomenology is unique to this type of models.
△ Less
Submitted 8 May, 2014;
originally announced May 2014.
-
Numerical Simulations of the Dark Universe: State of the Art and the Next Decade
Authors:
Michael Kuhlen,
Mark Vogelsberger,
Raul Angulo
Abstract:
We present a review of the current state of the art of cosmological dark matter simulations, with particular emphasis on the implications for dark matter detection efforts and studies of dark energy. This review is intended both for particle physicists, who may find the cosmological simulation literature opaque or confusing, and for astro-physicists, who may not be familiar with the role of simula…
▽ More
We present a review of the current state of the art of cosmological dark matter simulations, with particular emphasis on the implications for dark matter detection efforts and studies of dark energy. This review is intended both for particle physicists, who may find the cosmological simulation literature opaque or confusing, and for astro-physicists, who may not be familiar with the role of simulations for observational and experimental probes of dark matter and dark energy. Our work is complementary to the contribution by M. Baldi in this issue, which focuses on the treatment of dark energy and cosmic acceleration in dedicated N-body simulations. Truly massive dark matter-only simulations are being conducted on national supercomputing centers, employing from several billion to over half a trillion particles to simulate the formation and evolution of cosmologically representative volumes (cosmic scale) or to zoom in on individual halos (cluster and galactic scale). These simulations cost millions of core-hours, require tens to hundreds of terabytes of memory, and use up to petabytes of disk storage. The field is quite internationally diverse, with top simulations having been run in China, France, Germany, Korea, Spain, and the USA. Predictions from such simulations touch on almost every aspect of dark matter and dark energy studies, and we give a comprehensive overview of this connection. We also discuss the limitations of the cold and collisionless DM-only approach, and describe in some detail efforts to include different particle physics as well as baryonic physics in cosmological galaxy formation simulations, including a discussion of recent results highlighting how the distribution of dark matter in halos may be altered. We end with an outlook for the next decade, presenting our view of how the field can be expected to progress. (abridged)
△ Less
Submitted 26 October, 2012; v1 submitted 25 September, 2012;
originally announced September 2012.
-
Characterization of Dark-Matter-induced anisotropies in the diffuse gamma-ray background
Authors:
Mattia Fornasa,
Jesus Zavala,
Miguel A. Sanchez-Conde,
Jennifer M. Siegal-Gaskins,
Timur Delahaye,
Francisco Prada,
Mark Vogelsberger,
Fabio Zandanel,
Carlos S. Frenk
Abstract:
The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying a…
▽ More
The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing template all-sky gamma-ray maps for the emission produced in the galactic halo and its substructures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art N-body simulations of cosmic structure formation, namely Millennium-II for extragalactic (sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of extrapolation to account for (sub)structures that are below the resolution limit of the simulations, allowing us to estimate the total emission all the way down to the minimal self-bound halo mass. We describe in detail the features appearing in the APS of our template maps and we estimate the effect of various uncertainties such as the value of the minimal halo mass, the fraction of substructures hosted in a halo and the shape of the DM density profile. Our results indicate that the fluctuation APS of the DM-induced emission is of the same order as the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the comparison with the measured anisotropy. We also quantify the uncertainties affecting our results, finding "theoretical error bands" spanning more than two orders of magnitude and dominated (for a given particle physics model) by our lack of knowledge of the abundance of low-mass (sub)halos.
△ Less
Submitted 2 July, 2012;
originally announced July 2012.
-
Dark Matter implications of the Fermi-LAT measurement of anisotropies in the diffuse gamma-ray background: status report
Authors:
Mattia Fornasa,
Jesus Zavala,
Miguel A. Sanchez-Conde,
Francisco Prada,
Mark Vogelsberger
Abstract:
For the first time, the Fermi-LAT measured the angular power spectrum (APS) of anisotropies in the diffuse gamma-ray background. The data is found to be broadly compatible with a model with contributions from the point sources in the 1-year catalog, the galactic diffuse background, and the extragalactic isotropic emission; however deviations are present at both large and small angular scales. In t…
▽ More
For the first time, the Fermi-LAT measured the angular power spectrum (APS) of anisotropies in the diffuse gamma-ray background. The data is found to be broadly compatible with a model with contributions from the point sources in the 1-year catalog, the galactic diffuse background, and the extragalactic isotropic emission; however deviations are present at both large and small angular scales. In this study, we complement the model with a contribution from Dark Matter (DM) whose distribution is modeled exploiting the results of the most recent N-body simulations, considering the contribution of extragalactic halos and subhalos (from Millennium-II) and of galactic substructures (from Aquarius). With the use of the Fermi Science Tools, these simulations serve as templates to produce mock gamma-ray count maps for DM gamma-ray emission, both in the case of an annihilating and a decaying DM candidate. The APS will then be computed and compared with the Fermi-LAT results to derive constraints on the DM particle physics properties. The possible systematic due to an imperfect model of the galactic foreground is also studied and taken into account properly. The present paper reports on the status of the project.
△ Less
Submitted 3 October, 2011;
originally announced October 2011.
-
The cosmic X-ray and gamma-ray background from dark matter annihilation
Authors:
Jesus Zavala,
Mark Vogelsberger,
Tracy R. Slatyer,
Abraham Loeb,
Volker Springel
Abstract:
(Abridged) The extragalactic background light (EBL) observed at multiple wavelengths is a promising tool to probe the nature of dark matter since it might contain a significant contribution from gamma-rays produced promptly by dark matter annihilation. Additionally, the electrons and positrons produced in the annihilation give energy to the CMB photons to populate the EBL with X-rays and gamma-ray…
▽ More
(Abridged) The extragalactic background light (EBL) observed at multiple wavelengths is a promising tool to probe the nature of dark matter since it might contain a significant contribution from gamma-rays produced promptly by dark matter annihilation. Additionally, the electrons and positrons produced in the annihilation give energy to the CMB photons to populate the EBL with X-rays and gamma-rays. We here create full-sky maps of the radiation from both of these contributions using the high-resolution Millennium-II simulation. We use upper limits on the contributions of unknown sources to the EBL to constrain the intrinsic properties of dark matter using a model-independent approach that can be employed as a template to test different particle physics models (including those with a Sommerfeld enhancement). These upper limits are based on observations spanning eight orders of magnitude in energy (from soft X-rays measured by CHANDRA to gamma-rays measured by Fermi), and on expectations for the contributions from blazars and star forming galaxies. To exemplify this approach, we analyze a set of benchmark Sommerfeld-enhanced models that give the correct dark matter abundance, satisfy CMB constraints, and fit the cosmic ray spectra measured by PAMELA and Fermi without any contribution from local subhalos. We find that these models are in conflict with the EBL constraints unless the contribution of unresolved subhalos is small and the annihilation signal dominates the EBL. We conclude that provided the collisionless cold dark matter paradigm is accurate, even for conservative estimates of the contribution from unresolved subhalos and astrophysical backgrounds, the EBL is at least as sensitive a probe of these types of scenarios as the CMB. Our results disfavor an explanation of the positron excess measured by PAMELA based only on dark matter annihilation in the smooth Galactic halo.
△ Less
Submitted 2 June, 2011; v1 submitted 3 March, 2011;
originally announced March 2011.
-
Consistent Scenarios for Cosmic-Ray Excesses from Sommerfeld-Enhanced Dark Matter Annihilation
Authors:
Douglas P. Finkbeiner,
Lisa Goodenough,
Tracy R. Slatyer,
Mark Vogelsberger,
Neal Weiner
Abstract:
Anomalies in direct and indirect detection have motivated models of dark matter consisting of a multiplet of nearly-degenerate states, coupled by a new GeV-scale interaction. We perform a careful analysis of the thermal freezeout of dark matter annihilation in such a scenario. We compute the range of "boost factors" arising from Sommerfeld enhancement in the local halo for models which produce the…
▽ More
Anomalies in direct and indirect detection have motivated models of dark matter consisting of a multiplet of nearly-degenerate states, coupled by a new GeV-scale interaction. We perform a careful analysis of the thermal freezeout of dark matter annihilation in such a scenario. We compute the range of "boost factors" arising from Sommerfeld enhancement in the local halo for models which produce the correct relic density, and show the effect of including constraints on the saturated enhancement from the cosmic microwave background (CMB). We find that boost factors from Sommerfeld enhancement of up to ~800 are possible in the local halo. When the CMB bounds on the saturated enhancement are applied, the maximal boost factor is reduced to ~400 for 1-2 TeV dark matter and sub-GeV force carriers, but remains large enough to explain the observed Fermi and PAMELA electronic signals. We describe regions in the DM mass-boost factor plane where the cosmic ray data is well fit for a range of final states, and show that Sommerfeld enhancement alone is enough to provide the large annihilation cross sections required to fit the data, although for light mediator masses (less than ~200 MeV) there is tension with the CMB constraints in the absence of astrophysical boost factors from substructure. Additionally, we consider the circumstances under which WIMPonium formation is relevant and find for heavy WIMPs (greater than ~2 TeV) and soft-spectrum annihilation channels it can be an important consideration; we find regions with dark matter mass greater than 2.8 TeV that are consistent with the CMB bounds and have ~600-700 present-day boost factors.
△ Less
Submitted 12 April, 2011; v1 submitted 12 November, 2010;
originally announced November 2010.
-
Relic density and CMB constraints on dark matter annihilation with Sommerfeld enhancement
Authors:
Jesus Zavala,
Mark Vogelsberger,
Simon D. M. White
Abstract:
We calculate how the relic density of dark matter particles is altered when their annihilation is enhanced by the Sommerfeld mechanism due to a Yukawa interaction between the annihilating particles. Maintaining a dark matter abundance consistent with current observational bounds requires the normalization of the s-wave annihilation cross section to be decreased compared to a model without enhanc…
▽ More
We calculate how the relic density of dark matter particles is altered when their annihilation is enhanced by the Sommerfeld mechanism due to a Yukawa interaction between the annihilating particles. Maintaining a dark matter abundance consistent with current observational bounds requires the normalization of the s-wave annihilation cross section to be decreased compared to a model without enhancement. The level of suppression depends on the specific parameters of the particle model, with the kinetic decoupling temperature having the most effect. We find that the cross section can be reduced by as much as an order of magnitude for extreme cases. We also compute the mu-type distortion of the CMB energy spectrum caused by energy injection from such Sommerfeld-enhanced annihilation. Our results indicate that in the vicinity of resonances, associated with bound states, distortions can be large enough to be excluded by the upper limit |mu|<9.0x10^(-5) found by the COBE/FIRAS experiment.
△ Less
Submitted 4 March, 2010; v1 submitted 28 October, 2009;
originally announced October 2009.
