-
PRIMER: JWST/MIRI reveals the evolution of star-forming structures in galaxies at z<2.5
Authors:
Yipeng Lyu,
Benjamin Magnelli,
David Elbaz,
Pablo G. Pérez-González,
Camila Correa,
Emanuele Daddi,
Carlos Gómez-Guijarro,
James S. Dunlop,
Norman A. Grogin,
Anton M. Koekemoer,
Derek J. McLeod,
Shiying Lu
Abstract:
The stellar structures of star-forming galaxies (SFGs) undergo significant size growth during their mass assembly and must pass through a compaction phase as they evolve into quiescent galaxies (QGs). To shed light on the mechanisms behind this structural evolution, we study the morphology of the star-forming components of 665 SFGs at 0<z<2.5 measured using JWST/MIRI observation and compare them w…
▽ More
The stellar structures of star-forming galaxies (SFGs) undergo significant size growth during their mass assembly and must pass through a compaction phase as they evolve into quiescent galaxies (QGs). To shed light on the mechanisms behind this structural evolution, we study the morphology of the star-forming components of 665 SFGs at 0<z<2.5 measured using JWST/MIRI observation and compare them with the morphology of their stellar components taken from the literature. The stellar and star-forming components of most SFGs (66%) have extended disk-like structures that are aligned with each other and are of the same size. The star-forming components of these galaxies follow a mass-size relation, similar to that followed by their stellar components. At the highest mass, the optical Sérsic index of these SFGs increases to 2.5, suggesting the presence of a dominant stellar bulge. Because their star-forming components remain disk-like, these bulges cannot have formed by secular in-situ growth. We identify a second population of galaxies lying below the MIR mass-size relation, with compact star-forming components embedded in extended stellar components (EC galaxy). These galaxies are overall rare (15%) but become more dominant (30%) at high mass ($>10^{10.5}M_\odot$). The compact star-forming components of these galaxies are also concentrated and slightly spheroidal, suggesting that this compaction phase can build dense bulge in-situ. Finally, we identify a third population of SFGs (19%), with both compact stellar and star-forming components. The density of their stellar cores resemble those of QGs and are compatible with being the descendants of EC galaxy. Overall, the structural evolution of SFGs is mainly dominated by a secular inside-out growth, which can, however, be interrupted by violent compaction phase(s) that can build dominant stellar bulges like those in massive SFGs or QGs.
△ Less
Submitted 17 June, 2024;
originally announced June 2024.
-
TangoSIDM Project: Is the Stellar Mass Tully-Fisher relation consistent with SIDM?
Authors:
Camila Correa,
Matthieu Schaller,
Joop Schaye,
Sylvia Ploeckinger,
Josh Borrow,
Yannick Bahe
Abstract:
Self-interacting dark matter (SIDM) has the potential to significantly influence galaxy formation in comparison to the cold, collisionless dark matter paradigm (CDM), resulting in observable effects. This study aims to elucidate this influence and to demonstrate that the stellar mass Tully-Fisher relation imposes robust constraints on the parameter space of velocity-dependent SIDM models. We prese…
▽ More
Self-interacting dark matter (SIDM) has the potential to significantly influence galaxy formation in comparison to the cold, collisionless dark matter paradigm (CDM), resulting in observable effects. This study aims to elucidate this influence and to demonstrate that the stellar mass Tully-Fisher relation imposes robust constraints on the parameter space of velocity-dependent SIDM models. We present a new set of cosmological hydrodynamical simulations that include the SIDM scheme from the TangoSIDM project and the SWIFT-EAGLE galaxy formation model. Two cosmological simulations suites were generated: one (Reference model) which yields good agreement with the observed $z=0$ galaxy stellar mass function, galaxy mass-size relation, and stellar-to-halo mass relation; and another (WeakStellarFB model) in which the stellar feedback is less efficient, particularly for Milky Way-like systems. Both galaxy formation models were simulated under four dark matter cosmologies: CDM, SIDM with two different velocity-dependent cross sections, and SIDM with a constant cross section. While SIDM does not modify global galaxy properties such as stellar masses and star formation rates, it does make the galaxies more extended. In Milky Way-like galaxies, where baryons dominate the central gravitational potential, SIDM thermalises, causing dark matter to accumulate in the central regions. This accumulation results in density profiles that are steeper than those produced in CDM from adiabatic contraction. The enhanced dark matter density in the central regions of galaxies causes a deviation in the slope of the Tully-Fisher relation, which significantly diverges from the observational data. In contrast, the Tully-Fisher relation derived from CDM models aligns well with observations.
△ Less
Submitted 14 March, 2024;
originally announced March 2024.
-
Influence of local structure on relic neutrino abundances and anisotropies
Authors:
Fabian Zimmer,
Camila A. Correa,
Shin'ichiro Ando
Abstract:
Gravitational potentials of the Milky Way and extragalactic structures can influence the propagation of the cosmic neutrino background (CNB). Of particular interest to future CNB observatories, such as PTOLEMY, is the CNB number density on Earth. In this study, we have developed a simulation framework that maps the trajectories of relic neutrinos as they move through the local gravitational enviro…
▽ More
Gravitational potentials of the Milky Way and extragalactic structures can influence the propagation of the cosmic neutrino background (CNB). Of particular interest to future CNB observatories, such as PTOLEMY, is the CNB number density on Earth. In this study, we have developed a simulation framework that maps the trajectories of relic neutrinos as they move through the local gravitational environment. The potentials are based on the dark matter halos found in state-of-the-art cosmological N-body simulations, resulting in a more nuanced and realistic input than the previously employed analytical models. We find that the complex dark matter distributions, along with their dynamic evolution, influence the abundance and anisotropies of the CNB in ways unaccounted for by earlier analytical methods. Importantly, these cosmological simulations contain multiple instances of Milky Way-like halos that we employ to model a variety of gravitational landscapes. Consequently, we notice a variation in the CNB number densities that can be primarily attributed to the differences in the masses of these individual halos. For neutrino masses between $0.01$ and $0.3$ eV, we note clustering factors within the range of $1+\mathcal{O}(10^{-3})$ to $1+\mathcal{O}(1)$. Furthermore, the asymmetric nature of the underlying dark matter distributions within the halos results in not only overdense, but intriguingly, underdense regions within the full-sky anisotropy maps. Gravitational clustering appears to have a significant impact on the angular power spectra of these maps, leading to orders of magnitude more power on smaller scales beyond multipoles of $\ell = 3$ when juxtaposed against predictions by primordial fluctuations. We discuss how our results reshape our understanding of relic neutrino clustering and how this might affect observability of future CNB observatories such as PTOLEMY.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
Authors:
Matthieu Schaller,
Josh Borrow,
Peter W. Draper,
Mladen Ivkovic,
Stuart McAlpine,
Bert Vandenbroucke,
Yannick Bahé,
Evgenii Chaikin,
Aidan B. G. Chalk,
Tsang Keung Chan,
Camila Correa,
Marcel van Daalen,
Willem Elbers,
Pedro Gonnet,
Loïc Hausammann,
John Helly,
Filip Huško,
Jacob A. Kegerreis,
Folkert S. J. Nobels,
Sylvia Ploeckinger,
Yves Revaz,
William J. Roper,
Sergio Ruiz-Bonilla,
Thomas D. Sandnes,
Yolan Uyttenhove
, et al. (2 additional authors not shown)
Abstract:
Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gra…
▽ More
Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code SWIFT. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. SWIFT also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarise the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with SWIFT.
△ Less
Submitted 29 March, 2024; v1 submitted 22 May, 2023;
originally announced May 2023.
-
Guess the cheese flavour by the size of its holes: A cosmological test using the abundance of Popcorn voids
Authors:
Dante J. Paz,
Carlos M. Correa,
Sebastián R. Gualpa,
Andres N. Ruiz,
Carlos S. Bederián,
R. Dario Graña,
Nelson D. Padilla
Abstract:
We present a new definition of cosmic void and a publicly available code with the algorithm that implements it. Underdense regions are defined as free-form objects, called popcorn voids, made from the union of spheres of maximum volume with a given joint integrated underdensity contrast.The method is inspired by the excursion-set theory and consequently no rescaling processing is needed, the remov…
▽ More
We present a new definition of cosmic void and a publicly available code with the algorithm that implements it. Underdense regions are defined as free-form objects, called popcorn voids, made from the union of spheres of maximum volume with a given joint integrated underdensity contrast.The method is inspired by the excursion-set theory and consequently no rescaling processing is needed, the removal of overlapping voids and objects with sizes below the shot noise threshold is inherent in the algorithm. The abundance of popcorn voids in the matter field can be fitted using the excursion-set theory provided the relationship between the linear density contrast of the barrier and the threshold used in void identification is modified relative to the spherical evolution model. We also analysed the abundance of voids in biased tracer samples in redshift space. We show how the void abundance can be used to measure the geometric distortions due to the assumed fiducial cosmology, in a test similar to an Alcock-Paczyński test. Using the formalism derived from previous works, we show how to correct the abundance of popcorn voids for redshift-space distortion effects. Using this treatment, in combination with the excursion-set theory, we demonstrate the feasibility of void abundance measurements as cosmological probes. We obtain unbiased estimates of the target parameters, albeit with large degeneracies in the parameter space. Therefore, we conclude that the proposed test in combination with other cosmological probes has potential to improve current cosmological parameter constraints.
△ Less
Submitted 3 April, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
Cosmic voids as cosmological laboratories
Authors:
Carlos Mauricio Correa
Abstract:
Cosmic voids are promising cosmological laboratories for studying the dark energy phenomenon and alternative gravity theories. They are receiving special attention nowadays in view of the new generation of galaxy spectroscopic surveys, which are covering an unprecedented volume and redshift range. There are two primary statistics in void studies: (i) the void size function, which characterises the…
▽ More
Cosmic voids are promising cosmological laboratories for studying the dark energy phenomenon and alternative gravity theories. They are receiving special attention nowadays in view of the new generation of galaxy spectroscopic surveys, which are covering an unprecedented volume and redshift range. There are two primary statistics in void studies: (i) the void size function, which characterises the abundance of voids, and (ii) the void-galaxy cross-correlation function, which contains information about the density and velocity fields in these regions. However, it is necessary a complete description of the effects of geometrical (Alcock-Paczynski effect, AP) and dynamical (Kaiser effect, RSD) distortions around voids in order to design reliable cosmological tests based on these statistics. Observational measurements show prominent anisotropic patterns that lead to biased cosmological constraints if they are not properly modelled. This thesis addresses this problematic by presenting a theoretical and statistical framework based on dynamical and cosmological foundations capable of describing all the underlying effects involved: the expansion effect (t-RSD), the off-centring effect (v-RSD), the AP-volume effect and the ellipticity effect (e-RSD). These effects can be understood by studying the mapping of voids between real and redshift space. In this way, we lay the foundations for a proper modelling of the aforementioned statistics. In addition, we present a new cosmological test based on two perpendicular projections of the correlation function. The method is fiducial-cosmology free, which allows us to effectively break any possible degeneracy between the cosmological parameters involved. Moreover, it allows us to significantly reduce the number of mock catalogues needed to estimate covariances.
△ Less
Submitted 31 October, 2022;
originally announced October 2022.
-
Towards reconstructing the halo clustering and halo mass function of N-body simulations using neural ratio estimation
Authors:
Androniki Dimitriou,
Christoph Weniger,
Camila A. Correa
Abstract:
High-resolution cosmological N-body simulations are excellent tools for modelling the formation and clustering of dark matter haloes. These simulations suggest complex physical theories of halo formation governed by a set of effective physical parameters. Our goal is to extract these parameters and their uncertainties in a Bayesian context. We make a step towards automatising this process by direc…
▽ More
High-resolution cosmological N-body simulations are excellent tools for modelling the formation and clustering of dark matter haloes. These simulations suggest complex physical theories of halo formation governed by a set of effective physical parameters. Our goal is to extract these parameters and their uncertainties in a Bayesian context. We make a step towards automatising this process by directly comparing dark matter density projection maps extracted from cosmological simulations, with density projections generated from an analytical halo model. The model is based on a toy implementation of two body correlation functions. To accomplish this we use marginal neural ratio estimation, an algorithm for simulation-based inference that allows marginal posteriors to be estimated by approximating marginal likelihood-to-evidence ratios with a neural network. In this case, we train a neural network with mock images to identify the correct values of the physical parameters that produced a given image. Using the trained neural network on cosmological N-body simulation images we are able to reconstruct the halo mass function, to generate mock images similar to the N-body simulation images and to identify the lowest mass of the haloes of those images, provided that they have the same clustering with our training data. Our results indicate that this is a promising approach in the path towards developing cosmological simulations assisted by neural networks.
△ Less
Submitted 22 June, 2022;
originally announced June 2022.
-
TangoSIDM: Tantalizing models of Self-Interacting Dark Matter
Authors:
Camila A. Correa,
Matthieu Schaller,
Sylvia Ploeckinger,
Noemi Anau Montel,
Christoph Weniger,
Shinichiro Ando
Abstract:
We introduce the TangoSIDM project, a suite of cosmological simulations of structure formation in a $Λ$-Self-Interacting Dark Matter (SIDM) universe. TangoSIDM explores the impact of large dark matter (DM) scattering cross sections over dwarf galaxy scales. Motivated by DM interactions that follow a Yukawa potential, the cross section per unit mass, $σ/m_χ$, assumes a velocity dependent form that…
▽ More
We introduce the TangoSIDM project, a suite of cosmological simulations of structure formation in a $Λ$-Self-Interacting Dark Matter (SIDM) universe. TangoSIDM explores the impact of large dark matter (DM) scattering cross sections over dwarf galaxy scales. Motivated by DM interactions that follow a Yukawa potential, the cross section per unit mass, $σ/m_χ$, assumes a velocity dependent form that avoids violations of current constraints on large scales. We demonstrate that our implementation accurately models not only core formation in haloes, but also gravothermal core collapse. For central haloes in cosmological volumes, frequent DM particle collisions isotropise the particles orbit, making them largely spherical. We show that the velocity-dependent $σ/m_χ$ models produce a large diversity in the circular velocities of satellites haloes, with the spread in velocities increasing as the cross sections reach 20, 60 and 100 cm$^2$/g in $10^9~\rm{M}_{\odot}$ haloes. The large variation in the haloes internal structure is driven by DM particles interactions, causing in some haloes the formation of extended cores, whereas in others gravothermal core collapse. We conclude that the SIDM models from the Tango project offer a promising explanation for the diversity in the density and velocity profiles of observed dwarf galaxies.
△ Less
Submitted 22 June, 2022;
originally announced June 2022.
-
Cosmology with cosmic voids
Authors:
Carlos M. Correa,
Dante J. Paz
Abstract:
Cosmic voids constitute promising cosmological laboratories. However, a full description of all the redshift-space effects that affect observational measurements is mandatory in order to obtain unbiased cosmological constraints. We make a description in a nutshell of these effects and lay the theoretical foundations for designing reliable cosmological tests based on the void size function and the…
▽ More
Cosmic voids constitute promising cosmological laboratories. However, a full description of all the redshift-space effects that affect observational measurements is mandatory in order to obtain unbiased cosmological constraints. We make a description in a nutshell of these effects and lay the theoretical foundations for designing reliable cosmological tests based on the void size function and the void-galaxy cross-correlation function. We show that modern spectroscopic surveys offer a high signal-to-noise ratio to detect and study them.
△ Less
Submitted 26 May, 2022;
originally announced May 2022.
-
Estimating the warm dark matter mass from strong lensing images with truncated marginal neural ratio estimation
Authors:
Noemi Anau Montel,
Adam Coogan,
Camila Correa,
Konstantin Karchev,
Christoph Weniger
Abstract:
Precision analysis of galaxy-galaxy strong gravitational lensing images provides a unique way of characterizing small-scale dark matter halos, and could allow us to uncover the fundamental properties of dark matter's constituents. Recently, gravitational imaging techniques made it possible to detect a few heavy subhalos. However, gravitational lenses contain numerous subhalos and line-of-sight hal…
▽ More
Precision analysis of galaxy-galaxy strong gravitational lensing images provides a unique way of characterizing small-scale dark matter halos, and could allow us to uncover the fundamental properties of dark matter's constituents. Recently, gravitational imaging techniques made it possible to detect a few heavy subhalos. However, gravitational lenses contain numerous subhalos and line-of-sight halos, whose subtle imprint is extremely difficult to detect individually. Existing methods for marginalizing over this large population of sub-threshold perturbers to infer population-level parameters are typically computationally expensive, or require compressing observations into hand-crafted summary statistics, such as a power spectrum of residuals. Here, we present the first analysis pipeline to combine parametric lensing models and a recently-developed neural simulation-based inference technique called truncated marginal neural ratio estimation (TMNRE) to constrain the warm dark matter halo mass function cutoff scale directly from multiple lensing images. Through a proof-of-concept application to simulated data, we show that our approach enables empirically testable inference of the dark matter cutoff mass through marginalization over a large population of realistic perturbers that would be undetectable on their own, and over lens and source parameters uncertainties. To obtain our results, we combine the signal contained in a set of images with Hubble Space Telescope resolution. Our results suggest that TMNRE can be a powerful approach to put tight constraints on the mass of warm dark matter in the multi-keV regime, which will be relevant both for existing lensing data and in the large sample of lenses that will be delivered by near-future telescopes.
△ Less
Submitted 14 November, 2022; v1 submitted 18 May, 2022;
originally announced May 2022.
-
Warm Dark Matter Constraints Using Milky-Way Satellite Observations and Subhalo Evolution Modeling
Authors:
Ariane Dekker,
Shin'ichiro Ando,
Camila A. Correa,
Kenny C. Y. Ng
Abstract:
Warm dark matter (WDM) can potentially explain small-scale observations that currently challenge the cold dark matter (CDM) model, as warm particles suppress structure formation due to free-streaming effects. Observing small-scale matter distribution provides a valuable way to distinguish between CDM and WDM. In this work, we use observations from the Dark Energy Survey and PanSTARRS1, which obser…
▽ More
Warm dark matter (WDM) can potentially explain small-scale observations that currently challenge the cold dark matter (CDM) model, as warm particles suppress structure formation due to free-streaming effects. Observing small-scale matter distribution provides a valuable way to distinguish between CDM and WDM. In this work, we use observations from the Dark Energy Survey and PanSTARRS1, which observe 270 Milky-Way satellites after completeness corrections. We test WDM models by comparing the number of satellites in the Milky Way with predictions derived from the Semi-Analytical SubHalo Inference ModelIng (SASHIMI) code, which we develop based on the extended Press-Schechter formalism and subhalos' tidal evolution prescription. We robustly rule out WDM with masses lighter than 4.4 keV at 95% confidence level for the Milky-Way halo mass of $10^{12} M_\odot$. The limits are a weak function of the (yet uncertain) Milky-Way halo mass, and vary as $m_{\rm WDM}>3.6$-$5.1$ keV for $(0.6$-$2.0) \times 10^{12} M_\odot$. For the sterile neutrinos that form a subclass of WDM, we obtain the constraints of $m_{ν_s}>11.6$ keV for the Milky-Way halo mass of $10^{12} M_{\odot}$. These results based on SASHIMI do not rely on any assumptions of galaxy formation physics or are not limited by numerical resolution. The models, therefore, offer a robust and fast way to constrain the WDM models. By applying a satellite forming condition, however, we can rule out the WDM mass lighter than 9.0 keV for the Milky-Way halo mass of $10^{12} M_\odot$.
△ Less
Submitted 19 January, 2023; v1 submitted 25 November, 2021;
originally announced November 2021.
-
Simulations of $^{60}$Fe entrained in ejecta from a near-Earth supernova: Effects of observer motion
Authors:
Evgenii Chaikin,
Alexander A. Kaurov,
Brian D. Fields,
Camila A. Correa
Abstract:
Recent studies have shown that live (not decayed) radioactive $^{60}$Fe is present in deep-ocean samples, Antarctic snow, lunar regolith, and cosmic rays. $^{60}$Fe represents supernova (SN) ejecta deposited in the Solar System around $3 \, \rm Myr$ ago, and recently an earlier pulse $\approx 7 \ \rm Myr $ ago has been found. These data point to one or multiple near-Earth SN explosions that presum…
▽ More
Recent studies have shown that live (not decayed) radioactive $^{60}$Fe is present in deep-ocean samples, Antarctic snow, lunar regolith, and cosmic rays. $^{60}$Fe represents supernova (SN) ejecta deposited in the Solar System around $3 \, \rm Myr$ ago, and recently an earlier pulse $\approx 7 \ \rm Myr $ ago has been found. These data point to one or multiple near-Earth SN explosions that presumably participated in the formation of the Local Bubble. We explore this theory using 3D high-resolution smooth-particle hydrodynamical simulations of isolated supernovae with ejecta tracers in a uniform interstellar medium (ISM). The simulation allows us to trace the supernova ejecta in gas form and those eject in dust grains that are entrained with the gas. We consider two cases of diffused ejecta: when the ejecta are well-mixed in the shock and when they are not. In the latter case, we find that these ejecta remain far behind the forward shock, limiting the distance to which entrained ejecta can be delivered to $\approx 100$ pc in an ISM with $n_\mathrm{H}=0.1\; \rm cm^{-3}$ mean hydrogen density. We show that the intensity and the duration of $^{60}$Fe accretion depend on the ISM density and the trajectory of the Solar System. Furthermore, we show the possibility of reproducing the two observed peaks in $^{60}$Fe concentration with this model by assuming two linear trajectories for the Solar System with $30$-km s$^{-1}$ velocity. The fact that we can reproduce the two observed peaks further supports the theory that the $^{60}$Fe signal was originated from near-Earth SNe.
△ Less
Submitted 6 February, 2022; v1 submitted 23 September, 2021;
originally announced September 2021.
-
Redshift-space effects in voids and their impact on cosmological tests. Part II: the void-galaxy cross-correlation function
Authors:
Carlos M. Correa,
Dante J. Paz,
Nelson D. Padilla,
Ariel G. Sánchez,
Andrés N. Ruiz,
Raúl E. Angulo
Abstract:
This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above…
▽ More
This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above the shot-noise level identified with a spherical void finder. Under this mapping, void properties are affected by three effects: (i) a systematic expansion as a consequence of the distortions induced by galaxy dynamics, (ii) the Alcock-Paczynski volume effect, which manifests as an overall expansion or contraction depending on the fiducial cosmology, and (iii) a systematic off-centring along the line of sight as a consequence of the distortions induced by void dynamics. We found that correlations are also affected by an additional source of distortions: the ellipticity of voids. This is the first time that distortions due to the off-centring and ellipticity effects are detected and quantified. With a simplified test, we verified that the Gaussian streaming model is still robust provided all these effects are taken into account, laying the foundations for improvements in current models in order to obtain unbiased cosmological constraints from spectroscopic surveys. Besides this practical importance, this analysis also encodes key information about the structure and dynamics of the Universe at the largest scales. Furthermore, some of the effects constitute cosmological probes by themselves, as is the case of the void ellipticity.
△ Less
Submitted 19 January, 2022; v1 submitted 2 July, 2021;
originally announced July 2021.
-
Revealing the physical properties of gas accreting to haloes in the EAGLE simulations
Authors:
Ruby J. Wright,
Claudia del P. Lagos,
Chris Power,
Camila A. Correa
Abstract:
The inflow of cosmological gas onto haloes, while challenging to directly observe and quantify, plays a fundamental role in the baryon cycle of galaxies. Using the EAGLE suite of hydrodynamical simulations, we present a thorough exploration of the physical properties of gas accreting onto haloes -- namely, its spatial characteristics, density, temperature, and metallicity. Classifying accretion as…
▽ More
The inflow of cosmological gas onto haloes, while challenging to directly observe and quantify, plays a fundamental role in the baryon cycle of galaxies. Using the EAGLE suite of hydrodynamical simulations, we present a thorough exploration of the physical properties of gas accreting onto haloes -- namely, its spatial characteristics, density, temperature, and metallicity. Classifying accretion as ``hot'' or `` cold'' based on a temperature cut of $10^{5.5}{\rm K}$, we find that the covering fraction ($f_{\rm cov}$) of cold-mode accreting gas is significantly lower than the hot-mode, with $z=0$ $f_{\rm cov}$ values of $\approx 50\%$ and $\approx 80\%$ respectively. Active Galactic Nuclei (AGN) feedback in EAGLE reduces inflow $f_{\rm cov}$ values by $\approx 10\%$, with outflows decreasing the solid angle available for accretion flows. Classifying inflow by particle history, we find that gas on first-infall onto a halo is metal-depleted by $\approx 2$~dex compared to pre-processed gas, which we find to mimic the circum-galactic medium (CGM) in terms of metal content. We also show that high (low) halo-scale gas accretion rates are associated with metal-poor (rich) CGM in haloes below $10^{12}M_{\odot}$, and that variation in halo-scale gas accretion rates may offer a physical explanation for the enhanced scatter in the star-forming main sequence at low ($\lesssim10^{9}M_{\odot}$) and high ($\gtrsim10^{10}M_{\odot}$) stellar masses. Our results highlight how gas inflow influences several halo- and galaxy-scale properties, and the need to combine kinematic and chemical data in order to confidently break the degeneracy between accreting and outgoing gas in CGM observations.
△ Less
Submitted 13 April, 2021; v1 submitted 22 February, 2021;
originally announced February 2021.
-
The dependence of the galaxy stellar-to-halo mass relation on galaxy morphology
Authors:
Camila A. Correa,
Joop Schaye
Abstract:
We investigate the dependence of the local galaxy stellar-to-halo mass relation (SHMR) on galaxy morphology. We use data from the Sloan Digital Sky Survey DR7 with morphological classifications from Galaxy Zoo, and compare with the EAGLE cosmological simulation. At fixed halo mass in the mass range $10^{11.7}-10^{12.9}M_{\odot}$, the median stellar masses of SDSS disc galaxies are up to a factor o…
▽ More
We investigate the dependence of the local galaxy stellar-to-halo mass relation (SHMR) on galaxy morphology. We use data from the Sloan Digital Sky Survey DR7 with morphological classifications from Galaxy Zoo, and compare with the EAGLE cosmological simulation. At fixed halo mass in the mass range $10^{11.7}-10^{12.9}M_{\odot}$, the median stellar masses of SDSS disc galaxies are up to a factor of 1.4 higher than the median masses of their elliptical counterparts. However, when we switch from the stellar masses from Kauffmann et al. to those calculated by Chang et al. or Brinchmann et al., the median SHMR from discs and ellipticals coincide in this mass range. For halo masses larger than $10^{13}M_{\odot}$, discs are less massive than ellipticals in same-mass haloes, regardless of whose stellar mass estimates we use. However, we find that for these high halo masses the results for discs may be affected by central/satellite misclassifications. The EAGLE simulation predicts that discs are up to a factor of 1.5 more massive than elliptical galaxies residing in same-mass haloes less massive than $10^{13}M_{\odot}$, in agreement with the Kauffmann et al. data. Haloes with masses between $10^{11.5}$ and $10^{12}M_{\odot}$, that host disc galaxies, were assembled earlier than those hosting ellipticals. This suggests that the discs are more massive because they had more time for gas accretion and star formation. In $10^{12}-10^{12.5}M_{\odot}$ haloes, the central black holes in elliptical galaxies grew faster and became more massive than their counterparts in disc galaxies. This suggests that in this halo mass range the ellipticals are less massive because AGN feedback ejected more of the halo's gas reservoir, reducing star formation, and suppressing the (re)growth of stellar discs.
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
Redshift-space effects in voids and their impact on cosmological tests. Part I: the void size function
Authors:
Carlos M. Correa,
Dante J. Paz,
Ariel G. Sánchez,
Andrés N. Ruiz,
Nelson D. Padilla,
Raúl E. Angulo
Abstract:
Voids are promising cosmological probes. Nevertheless, every cosmological test based on voids must necessarily employ methods to identify them in redshift space. Therefore, redshift-space distortions (RSD) and the Alcock-Paczynski effect (AP) have an impact on the void identification process itself generating distortion patterns in observations. Using a spherical void finder, we developed a statis…
▽ More
Voids are promising cosmological probes. Nevertheless, every cosmological test based on voids must necessarily employ methods to identify them in redshift space. Therefore, redshift-space distortions (RSD) and the Alcock-Paczynski effect (AP) have an impact on the void identification process itself generating distortion patterns in observations. Using a spherical void finder, we developed a statistical and theoretical framework to describe physically the connection between the identification in real and redshift space. We found that redshift-space voids above the shot noise level have a unique real-space counterpart spanning the same region of space, they are systematically bigger and their centres are preferentially shifted along the line of sight. The expansion effect is a by-product of RSD induced by tracer dynamics at scales around the void radius, whereas the off-centring effect constitutes a different class of RSD induced at larger scales by the global dynamics of the whole region containing the void. The volume of voids is also altered by the fiducial cosmology assumed to measure distances, this is the AP change of volume. These three systematics have an impact on cosmological statistics. In this work, we focus on the void size function. We developed a theoretical framework to model these effects and tested it with a numerical simulation, recovering the statistical properties of the abundance of voids in real space. This description depends strongly on cosmology. Hence, we lay the foundations for improvements in current models of the abundance of voids in order to obtain unbiased cosmological constraints from redshift surveys.
△ Less
Submitted 15 October, 2020; v1 submitted 23 July, 2020;
originally announced July 2020.
-
Constraining Velocity-dependent Self-Interacting Dark Matter with the Milky Way's dwarf spheroidal galaxies
Authors:
Camila A. Correa
Abstract:
The observed anti-correlation between the central dark matter (DM) densities of the bright Milky Way (MW) dwarf spheroidal galaxies (dSphs) and their orbital pericenter distances poses a potential signature of self-interacting dark matter (SIDM). In this work we investigate this possibility by analysing the range of SIDM scattering cross section per unit mass, $σ/m_χ$, able to explain such anti-co…
▽ More
The observed anti-correlation between the central dark matter (DM) densities of the bright Milky Way (MW) dwarf spheroidal galaxies (dSphs) and their orbital pericenter distances poses a potential signature of self-interacting dark matter (SIDM). In this work we investigate this possibility by analysing the range of SIDM scattering cross section per unit mass, $σ/m_χ$, able to explain such anti-correlation. We simulate the orbital evolution of dSphs subhaloes around the MW assuming an analytical form for the gravitational potential, adopting the proper motions from the Gaia mission and including a consistent characterization of gravitational tidal stripping. The evolution of the subhaloes density profile is modelled using the gravothermal fluid formalism, where DM particle collisions induce thermal conduction that depends on $σ/m_χ$. We find that models of dSphs, such as Carina and Fornax, reproduce the observed central DM densities with fixed $σ/m_χ$ ranging between $30$ and $50$ cm$^{2}$g$^{-1}$, whereas other dSphs prefer larger values ranging between $70$ and $100$ cm$^{2}$g$^{-1}$. These cross sections correlate with the average collision velocity of DM particles within each subhalo's core, so that systems modelled with large cross sections have lower collision velocities. We fit the cross section-velocity correlation with a SIDM particle model, where a DM particle of mass $m_χ=0.648\pm 0.154$ GeV interacts under the exchange of a light mediator of mass $m_φ=0.636\pm 0.055$ MeV, with the self-interactions being described by a Yukawa potential. The outcome is a cross section-velocity relation that explains the diverse DM profiles of MW dSph satellites and is consistent with observational constraints on larger scales.
△ Less
Submitted 6 July, 2020;
originally announced July 2020.
-
The effect of gas accretion on the radial gas metallicity profile of simulated galaxies
Authors:
Florencia Collacchioni,
Claudia D. P. Lagos,
Peter D. Mitchell,
Joop Schaye,
Emily Wisnioski,
Sofía A. Cora,
Camila A. Correa
Abstract:
We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where…
▽ More
We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where higher $\dot M_{\rm accr}$ are associated with more negative slopes. The slope of the RMPs depends more strongly on $\dot M_{\rm accr}$ than on stellar mass, star formation rate or gas fraction, suggesting $\dot M_{\rm accr}$ to be a more fundamental driver of the RMP slope of galaxies. We find that eliminating the dependence on stellar mass is essential for pinning down the properties that shape the slope of the RMP. Although $\dot M_{\rm accr}$ is the main property modulating the slope of the RMP, we find that it causes other correlations that are more easily testable observationally: at fixed stellar mass, galaxies with more negative RMP slopes tend to have higher gas fractions and SFRs, while galaxies with lower gas fractions and SFRs tend to have flatter metallicity profiles within an effective radius.
△ Less
Submitted 28 May, 2020; v1 submitted 11 October, 2019;
originally announced October 2019.
-
Non-fiducial cosmological test from geometrical and dynamical distortions around voids
Authors:
Carlos M. Correa,
Dante J. Paz,
Nelson D. Padilla,
Andrés N. Ruiz,
Raúl E. Angulo,
Ariel G. Sánchez
Abstract:
We present a new cosmological test using the distribution of galaxies around cosmic voids without assuming a fiducial cosmology. The test is based on a physical model for the void-galaxy cross-correlation function projected along and perpendicular to the line of sight. We treat correlations in terms of void-centric angular distances and redshift differences between void-galaxy pairs, hence it is n…
▽ More
We present a new cosmological test using the distribution of galaxies around cosmic voids without assuming a fiducial cosmology. The test is based on a physical model for the void-galaxy cross-correlation function projected along and perpendicular to the line of sight. We treat correlations in terms of void-centric angular distances and redshift differences between void-galaxy pairs, hence it is not necessary to assume a fiducial cosmology. This model reproduces the coupled dynamical (Kaiser effect, RSD) and geometrical (Alcock-Paczynski effect, GD) distortions that affect the correlation measurements. It also takes into account the scale mixing due to the projection ranges in both directions. The model is general, so it can be applied to an arbitrary cylindrical binning scheme, not only in the case of the projected correlations. It primarily depends on two cosmological parameters: $Ω_m$, the matter fraction of the Universe today (sensitive to GD), and $β$, the ratio between the growth rate factor of density perturbations and the tracer bias (sensitive to RSD). In the context of the new generation of galaxy spectroscopic surveys, we calibrated the test using the Millennium XXL simulation for different redshifts. The method successfully recovers the cosmological parameters. We studied the effect of measuring with different projection ranges, finding robust results up to wide ranges. The resulting data covariance matrices are relatively small, which reduces the noise in the Gaussian likelihood analysis and will allow the usage of a smaller number of mock catalogues. The performance evaluated in this work indicates that the developed method is a promising test to be applied on real data.
△ Less
Submitted 21 December, 2018; v1 submitted 29 November, 2018;
originally announced November 2018.
-
The origin of the red sequence galaxy population in the EAGLE simulation
Authors:
Camila A. Correa,
Joop Schaye,
James W. Trayford
Abstract:
We investigate the evolution in colour and morphology of the progenitors of red-sequence galaxies in the EAGLE cosmological hydrodynamical simulation. We quantify colours with $u^{*}-r^{*}$ intrinsic magnitudes and morphologies with a measure of the stellar kinematics. The time when galaxies moved onto the red sequence depends on their morphology. Disc-type galaxies tend to have become red during…
▽ More
We investigate the evolution in colour and morphology of the progenitors of red-sequence galaxies in the EAGLE cosmological hydrodynamical simulation. We quantify colours with $u^{*}-r^{*}$ intrinsic magnitudes and morphologies with a measure of the stellar kinematics. The time when galaxies moved onto the red sequence depends on their morphology. Disc-type galaxies tend to have become red during the last 3 Gyr, while elliptical-type galaxies joined the red sequence earlier, with half the sample already being red 5 Gyr ago. The time-scale, $τ_{\rm{Green}}$, of colour transition through the `green valley' depends weakly on the galaxy's morphological type. Elliptical-type galaxies cross the green valley slightly faster ($τ_{\rm{Green}}{\approx}1$ Gyr) than disc-type galaxies ($τ_{\rm{Green}}{\approx}$1.5 Gyr). While $τ_{\rm{Green}}$ is similar for central and satellite galaxies, for satellites $τ_{\rm{Green}}$ decreases with increasing stellar mass to host-halo mass ratio. Coupled with our finding that galaxies tend to become green after becoming satellites, this indicates that satellite-specific processes are important for quenching red-sequence galaxies. The last time central, elliptical-type red-sequence galaxies left the blue cloud is strongly correlated with the time the luminosity of the central black hole peaked, but this is not the case for discs. This suggests that AGN feedback is important for quenching ellipticals, particularly centrals, but not for discs. We find only a weak connection between transformations in colour and morphology.
△ Less
Submitted 16 October, 2018;
originally announced October 2018.
-
The star formation rate and stellar content contributions of morphological components in the EAGLE simulations
Authors:
James W. Trayford,
Carlos S. Frenk,
Tom Theuns,
Joop Schaye,
Camila Correa
Abstract:
The Hubble sequence provides a useful classification of galaxy morphology at low redshift. However, morphologies are not static, but rather evolve as the growth of structure proceeds through mergers, accretion and secular processes. We investigate how kinematically defined disc and spheroidal structures form and evolve in the EAGLE hydrodynamic simulation of galaxy formation. At high redshift most…
▽ More
The Hubble sequence provides a useful classification of galaxy morphology at low redshift. However, morphologies are not static, but rather evolve as the growth of structure proceeds through mergers, accretion and secular processes. We investigate how kinematically defined disc and spheroidal structures form and evolve in the EAGLE hydrodynamic simulation of galaxy formation. At high redshift most galaxies of all masses are asymmetric. By redshift $z\simeq 1.5$ the Hubble sequence is established and after this time most of the stellar mass is in spheroids whose contribution to the stellar mass budget continues to rise to the present day. The stellar mass fraction in discs peaks at $z\simeq 0.5$ but overall remains subdominant at all times although discs contribute most of the stellar mass in systems of mass $M^*\sim 10^{10.5}{\rm M_\odot}$ at $z \leq 1.5$. Star formation occurs predominantly in disc structures throughout most of cosmic time but morphological transformations rearrange stars, thus establishing the low-redshift morphological mix. Morphological transformations are common and we quantify the rates at which they occur. The rate of growth of spheroids decreases at $z < 2$ while the rate of decay of discs remains roughly constant at $z < 1$. Finally, we find that the prograde component of galaxies becomes increasingly dynamically cold with time.
△ Less
Submitted 28 May, 2018; v1 submitted 8 May, 2018;
originally announced May 2018.
-
The impact of feedback and the hot halo on the rates of gas accretion onto galaxies
Authors:
Camila A. Correa,
Joop Schaye,
Freeke van de Voort,
Alan R. Duffy,
J. Stuart B. Wyithe
Abstract:
We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and in…
▽ More
We investigate the physics that drives the gas accretion rates onto galaxies at the centers of dark matter haloes using the EAGLE suite of hydrodynamical cosmological simulations. We find that at redshifts $z{\le}2$ the accretion rate onto the galaxy increases with halo mass in the halo mass range $10^{10}-10^{11.7}M_{\odot}$, flattens between the halo masses $10^{11.7}-10^{12.7}M_{\odot}$, and increases again for higher-mass haloes. However, the galaxy gas accretion does not flatten at intermediate halo masses when AGN feedback is switched off. To better understand these trends, we develop a physically motivated semi-analytic model of galaxy gas accretion. We show that the flattening is produced by the rate of gas cooling from the hot halo. The ratio of the cooling radius and the virial radius does not decrease continuously with increasing halo mass as generally thought. While it decreases up to ${\sim}10^{13}M_{\odot}$ haloes, it increases for higher halo masses, causing an upturn in the galaxy gas accretion rate. This may indicate that in high-mass haloes AGN feedback is not sufficiently efficient. When there is no AGN feedback, the density of the hot halo is higher, the ratio of the cooling and virial radii does not decrease as much and the cooling rate is higher. Changes in the efficiency of stellar feedback can also increase or decrease the accretion rates onto galaxies. The trends can plausibly be explained by the re-accretion of gas ejected by progenitor galaxies and by the suppression of black hole growth, and hence AGN feedback, by stellar feedback.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
The formation of hot gaseous haloes around galaxies
Authors:
Camila A. Correa,
Joop Schaye,
J. Stuart B. Wyithe,
Alan R. Duffy,
Tom Theuns,
Robert A. Crain,
Richard Bower
Abstract:
We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo…
▽ More
We use a suite of hydrodynamical cosmological simulations from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) project to investigate the formation of hot hydrostatic haloes and their dependence on feedback mechanisms. We find that the appearance of a strong bimodality in the probability density function (PDF) of the ratio of the radiative cooling and dynamical times for halo gas provides a clear signature of the formation of a hot corona. Haloes of total mass $10^{11.5}-10^{12}{\rm{M}}_{\odot}$ develop a hot corona independent of redshift, at least in the interval $z=0-4$ where the simulation has sufficiently good statistics. We analyse the build up of the hot gas mass in the halo, $M_{\rm{hot}}$, as a function of halo mass and redshift and find that while more energetic galactic winds powered by SNe increases $M_{\rm{hot}}$, AGN feedback reduces it by ejecting gas from the halo. We also study the thermal properties of gas accreting onto haloes and measure the fraction of shock-heated gas as a function of redshift and halo mass. We develop analytic and semianalytic approaches to estimate a `critical halo mass', $M_{\rm{crit}}$, for hot halo formation. We find that the mass for which the heating rate produced by accretion shocks equals the radiative cooling rate, reproduces the mass above which haloes develop a significant hot atmosphere. This yields a mass estimate of $M_{\rm{crit}} \approx 10^{11.7}{\rm{M}}_{\odot}$ at $z=0$, which agrees with the simulation results. The value of $M_{\rm{crit}}$ depends more strongly on the cooling rate than on any of the feedback parameters.
△ Less
Submitted 6 September, 2017;
originally announced September 2017.
-
The relation between galaxy morphology and colour in the EAGLE simulation
Authors:
Camila A. Correa,
Joop Schaye,
Bart Clauwens,
Richard G. Bower,
Robert A. Crain,
Matthieu Schaller,
Tom Theuns,
Adrien C. R. Thob
Abstract:
We investigate the relation between kinematic morphology, intrinsic colour and stellar mass of galaxies in the EAGLE cosmological hydrodynamical simulation. We calculate the intrinsic u-r colours and measure the fraction of kinetic energy invested in ordered corotation of 3562 galaxies at z=0 with stellar masses larger than $10^{10}M_{\odot}$. We perform a visual inspection of gri-composite images…
▽ More
We investigate the relation between kinematic morphology, intrinsic colour and stellar mass of galaxies in the EAGLE cosmological hydrodynamical simulation. We calculate the intrinsic u-r colours and measure the fraction of kinetic energy invested in ordered corotation of 3562 galaxies at z=0 with stellar masses larger than $10^{10}M_{\odot}$. We perform a visual inspection of gri-composite images and find that our kinematic morphology correlates strongly with visual morphology. EAGLE produces a galaxy population for which morphology is tightly correlated with the location in the colour- mass diagram, with the red sequence mostly populated by elliptical galaxies and the blue cloud by disc galaxies. Satellite galaxies are more likely to be on the red sequence than centrals, and for satellites the red sequence is morphologically more diverse. These results show that the connection between mass, intrinsic colour and morphology arises from galaxy formation models that reproduce the observed galaxy mass function and sizes.
△ Less
Submitted 20 April, 2017;
originally announced April 2017.
-
The environmental dependence of gas accretion onto galaxies: quenching satellites through starvation
Authors:
Freeke van de Voort,
Yannick M. Bahé,
Richard G. Bower,
Camila A. Correa,
Robert A. Crain,
Joop Schaye,
Tom Theuns
Abstract:
Galaxies that have fallen into massive haloes may no longer be able to accrete gas from their surroundings, a process referred to as 'starvation' or 'strangulation' of satellites. We study the environmental dependence of gas accretion onto galaxies using the cosmological, hydrodynamical EAGLE simulation. We quantify the dependence of gas accretion on stellar mass, redshift, and environment, using…
▽ More
Galaxies that have fallen into massive haloes may no longer be able to accrete gas from their surroundings, a process referred to as 'starvation' or 'strangulation' of satellites. We study the environmental dependence of gas accretion onto galaxies using the cosmological, hydrodynamical EAGLE simulation. We quantify the dependence of gas accretion on stellar mass, redshift, and environment, using halo mass and galaxy overdensity as environmental indicators. We find a strong suppression, by many orders of magnitude, of the gas accretion rate in dense environments, primarily for satellite galaxies. This suppression becomes stronger at lower redshift. However, the scatter in accretion rates is very large for satellites. This is (at least in part) due to the variation in halocentric radius, since gas accretion is more suppressed at smaller radii. Central galaxies are influenced less strongly by their environment and exhibit less scatter in their gas accretion rates. The star formation rates of both centrals and satellites show similar behaviour to their gas accretion rates. The relatively small differences between gas accretion and star formation rates demonstrate that galaxies generally exhaust their gas reservoir somewhat faster at higher stellar mass, lower redshift, and in denser environments. We conclude that the environmental suppression of gas accretion could directly result in the quenching of star formation.
△ Less
Submitted 21 December, 2016; v1 submitted 11 November, 2016;
originally announced November 2016.
-
The accretion history of dark matter halos III: A physical model for the concentration-mass relation
Authors:
Camila A. Correa,
J. Stuart B. Wyithe,
Joop Schaye,
Alan R. Duffy
Abstract:
We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model combines an analytic model for the halo mass accretion history (MAH), based on extended Press Schechter (EPS) theory, with an empirical relation between concentration and formation time obtained through fits to the results of numerical simulatio…
▽ More
We present a semi-analytic, physically motivated model for dark matter halo concentration as a function of halo mass and redshift. The semi-analytic model combines an analytic model for the halo mass accretion history (MAH), based on extended Press Schechter (EPS) theory, with an empirical relation between concentration and formation time obtained through fits to the results of numerical simulations. Because the semi-analytic model is based on EPS theory, it can be applied to wide ranges in mass, redshift and cosmology. The resulting concentration-mass (c-M) relations are found to agree with the simulations, and because the model applies only to relaxed halos, they do not exhibit the upturn at high masses or high redshifts found by some recent works. We predict a change of slope in the z=0 c-M relation at a mass scale of $10^{11}\rm{M}_{\odot}$. We find that this is due to the change in the functional form of the halo MAH, which goes from being dominated by an exponential (for high-mass halos) to a power-law (for low-mass halos). During the latter phase, the core radius remains approximately constant, and the concentration grows due to the drop of the background density. We also analyse how the c-M relation predicted by this work affects the power produced by dark matter annihilation, finding that at z = 0 the power is two orders of magnitude lower than that obtained from extrapolating best-fitting c-M relations. We provide fits to the c-M relations as well as numerical routines to compute concentrations and MAHs.
△ Less
Submitted 3 August, 2015; v1 submitted 2 February, 2015;
originally announced February 2015.
-
The accretion history of dark matter halos II: The connections with the mass power spectrum and the density profile
Authors:
Camila A. Correa,
J. Stuart B. Wyithe,
Joop Schaye,
Alan R. Duffy
Abstract:
We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines an…
▽ More
We explore the relation between the structure and mass accretion histories of dark matter halos using a suite of cosmological simulations. We confirm that the formation time, defined as the time when the virial mass of the main progenitor equals the mass enclosed within the scale radius, correlates strongly with concentration. We provide a semi-analytic model for halo mass history that combines analytic relations with fits to simulations. This model has the functional form, $M(z) = M_{0}(1+z)^αe^{βz}$, where the parameters $α$ and $β$ are directly correlated with concentration. We then combine this model for the halo mass history with the analytic relations between $α$, $β$ and the linear power spectrum derived by Correa et al. (2014) to establish the physical link between halo concentration and the initial density perturbation field. Finally, we provide fitting formulas for the halo mass history as well as numerical routines, we derive the accretion rate as a function of halo mass, and we demonstrate how the halo mass history depends on cosmology and the adopted definition of halo mass.
△ Less
Submitted 3 August, 2015; v1 submitted 18 January, 2015;
originally announced January 2015.
-
The accretion history of dark matter halos I: The physical origin of the universal function
Authors:
Camila A. Correa,
J. Stuart B. Wyithe,
Joop Schaye,
Alan R. Duffy
Abstract:
Understanding the universal accretion history of dark matter halos is the first step towards determining the origin of their structure. We use the extended Press-Schechter formalism to derive the halo mass accretion history from the growth rate of initial density perturbations. We show that the halo mass history is well described by an exponential function of redshift in the high-redshift regime.…
▽ More
Understanding the universal accretion history of dark matter halos is the first step towards determining the origin of their structure. We use the extended Press-Schechter formalism to derive the halo mass accretion history from the growth rate of initial density perturbations. We show that the halo mass history is well described by an exponential function of redshift in the high-redshift regime. However, in the low-redshift regime the mass history follows a power law because the growth of density perturbations is halted in the dark energy dominated era due to the accelerated expansion of the Universe. We provide an analytic model that follows the expression $M(z) = M_{0}(1+z)^{af(M_{0})}e^{-f(M_{0})z}$, where $M_{0} = M(z=0)$, $a$ depends on cosmology and $f(M_{0})$ depends only on the linear matter power spectrum. The analytic model does not rely on calibration against numerical simulations and is suitable for any cosmology. We compare our model with the latest empirical models for the mass accretion history in the literature and find very good agreement. We provide numerical routines for the model online.
△ Less
Submitted 2 August, 2015; v1 submitted 18 September, 2014;
originally announced September 2014.