Impact Of Baryonic Physics Research Articles

The spatial distribution of Milky Way satellites, gaps in streams, and the nature of dark matter

ABSTRACT The spatial distribution of Milky Way (MW) subhaloes provides an important set of observables for testing cosmological models. These include the radial distribution of luminous satellites, planar configurations, and the abundance of dark subhaloes whose existence or absence is key to distinguishing among dark matter models. We use the cocoN-body simulations of cold dark matter (CDM) and 3.3 keV thermal relic warm dark matter (WDM) to predict the satellite spatial distribution in the limit that the impact of baryonic physics is minimal. We demonstrate that the radial distributions of CDM and 3.3 keV-WDM luminous satellites are identical if the minimum pre-infall halo mass to form a galaxy is >108.5 ${\, \rm M_\odot }$. The distribution of dark subhaloes is significantly more concentrated in WDM due to the absence of low mass, recently accreted substructures that typically inhabit the outer parts of a MW halo in CDM. We show that subhaloes of mass [107, 108] ${\, \rm M_\odot }$ and within 30 kpc of the centre are the stripped remnants of larger haloes in both models. Therefore, their abundance in WDM is 3× higher than one would anticipate from the overall WDM subhalo population. We estimate that differences between CDM and WDM concentration–mass relations can be probed for subhalo–stream impact parameters <2 kpc. Finally, we find that the impact of WDM on planes of satellites is likely negligible. Comprehensive comparisons with observations will require further work with high resolution, self-consistent hydrodynamical simulations.

Galaxy formation in the brane world I: overview and first results

ABSTRACT We carry out ‘full-physics’ hydrodynamical simulations of galaxy formation in the normal-branch Dvali–Gabadadze–Porrati (nDGP) braneworld model using a new modified version of the arepo code and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) that represent, respectively, weak and moderate departures from general relativity (GR), in boxes of sizes $62$ and $25\, h^{-1}{\rm Mpc}$ using 2 × 5123 dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to ∼25 per cent on large scales. Similar to their impact in the standard cosmological model (Λ cold dark matter), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to ∼15 per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction, and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above ∼0.1 times the halo radius.

Cosmological simulations of the same spiral galaxy: the impact of baryonic physics

ABSTRACT The interplay of star formation (SF) and supernova (SN) feedback in galaxy formation is a key element for understanding galaxy evolution. Since these processes occur at small scales, it is necessary to have sub-grid models that recover their evolution and environmental effects at the scales reached by cosmological simulations. In this work, we present the results of the Mochima simulation, where we simulate the same spiral galaxy inhabiting a Milky Way (MW) size halo in a cosmological environment changing the sub-grid models for SN feedback and SF. We test combinations of the Schmidt law and a multifreefall based SF with delayed cooling feedback or mechanical feedback. We reach a resolution of 35 pc in a zoom-in box of 36 Mpc. For this, we use the code $\rm{\small RAMSES}$ with the implementation of gas turbulence in time and trace the local hydrodynamical features of the star-forming gas. Finally, we compare the galaxies at redshift 0 with global and interstellar medium observations in the MW and local spiral galaxies. The simulations show successful comparisons with observations. Nevertheless, diverse galactic morphologies are obtained from different numerical implementations. We highlight the importance of detailed modelling of the SF and feedback processes, especially for simulations with a resolution that start to reach scales relevant for molecular cloud physics. Future improvements could alleviate the degeneracies exhibited in our simulated galaxies under different sub-grid models.

Exploring extensions to the standard cosmological model and the impact of baryons on small scales

ABSTRACT It has been claimed that the standard model of cosmology (ΛCDM) cannot easily account for a number of observations on relatively small scales, motivating extensions to the standard model. Here, we introduce a new suite of cosmological simulations that systematically explores three plausible extensions: warm dark matter, self-interacting dark matter, and a running of the scalar spectral index of density fluctuations. Current observational constraints are used to specify the additional parameters that come with these extensions. We examine a large range of observable metrics on small scales, including the halo mass function, density, and circular velocity profiles, the abundance of satellite subhaloes, and halo concentrations. For any given metric, significant degeneracies can be present between the extensions. In detail, however, the different extensions have quantitatively distinct mass and radial dependencies, suggesting that a multiprobe approach over a range of scales can be used to break the degeneracies. We also demonstrate that the relative effects on the radial density profiles in the different extensions (compared to the standard model) are converged down to significantly smaller radii than are the absolute profiles. We compare the derived cosmological trends with the impact of baryonic physics using the EAGLE and ARTEMIS simulations. Significant degeneracies are also present between baryonic physics and cosmological variations (with both having similar magnitude effects on some observables). Given the inherent uncertainties both in the modelling of galaxy formation physics and extensions to ΛCDM, a systematic and simultaneous exploration of both is strongly warranted.

The impact of baryonic physics and massive neutrinos on weak lensing peak statistics

ABSTRACT We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass Mν = 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of Mν. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

The Impact of Baryonic Physics on the Kinetic Sunyaev–Zel’dovich Effect

Abstract Poorly understood “baryonic physics” impacts our ability to predict the power spectrum of the kinetic Sunyaev–Zel’dovich (kSZ) effect. We study this in a sample high-resolution simulation of galaxy formation and feedback, Illustris. The high resolution of Illustris allows us to probe the kSZ power spectrum on multipoles . Strong AGN feedback in Illustris nearly wipes out gas fluctuations at and at late times, likely somewhat underpredicting the kSZ power generated at . The post-reionization kSZ power spectrum for Illustris is well-fit by over 10,000, somewhat lower than most other reported values but consistent with the analysis of Shaw et al. Our analysis of the bias of free electrons reveals subtle effects associated with the multi-phase gas physics and stellar fractions that affect even linear scales. In particular, there are fewer electrons in biased galaxies, due to gas-cooling and star formation, and this leads to an electron bias of less than one, even at low wavenumbers. The combination of bias and electron fraction that determines the overall suppression is relatively constant, , but more simulations are needed to see if this is Illustris-specific. By separating the kSZ power into different terms, we find that at least 6% (10%) of the signal at ℓ = 3000 (10,000) comes from non-Gaussian connected four-point density and velocity correlations, , even without correcting for the Illustris simulation box-size. A challenge going forward will be accurately modeling long-wave velocity modes simultaneously with Illustris-like high resolution to capture the complexities of galaxy formation and its correlations with large-scale flows.

The structure and assembly history of cluster-sized haloes in self-interacting dark matter

We perform dark-matter-only simulations of 28 relaxed massive cluster-sized haloes for Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, to study structural differences between the models at large radii, where the impact of baryonic physics is expected to be very limited. We find that the distributions for the radial profiles of the density, ellipsoidal axis ratios, and velocity anisotropies ($\beta$) of the haloes differ considerably between the models (at the $\sim1\sigma$ level), even at $\gtrsim10\%$ of the virial radius, if the self-scattering cross section is $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$. Direct comparison with observationally inferred density profiles disfavours SIDM for $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, but in an intermediate radial range ($\sim3\%$ of the virial radius), where the impact of baryonic physics is uncertain. At this level of the cross section, we find a narrower $\beta$ distribution in SIDM, clearly skewed towards isotropic orbits, with no SIDM (90\% of CDM) haloes having $\beta>0.12$ at $7\%$ of the virial radius. We estimate that with an observational sample of $\sim30$ ($\sim10^{15}$ M$_\odot$) relaxed clusters, $\beta$ can potentially be used to put competitive constraints on SIDM, once observational uncertainties improve by a factor of a few. We study the suppression of the memory of halo assembly history in SIDM clusters. For $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, we find that this happens only in the central halo regions ($\sim1/4$ of the scale radius of the halo), and only for haloes that assembled their mass within this region earlier than a formation redshift $z_f\sim2$. Otherwise, the memory of assembly remains and is reflected in ways similar to CDM, albeit with weaker trends.

Implications of hydrodynamical simulations for the interpretation of direct dark matter searches

In recent years, realistic hydrodynamical simulations of galaxies like the Milky Way have become available, enabling a reliable estimate of the dark matter density and velocity distribution in the Solar neighborhood. We review here the status of hydrodynamical simulations and their implications for the interpretation of direct dark matter searches. We focus in particular on: the criteria to identify Milky Way-like galaxies; the impact of baryonic physics on the dark matter velocity distribution; the possible presence of substructures like clumps, streams, or dark disks; and on the implications for the direct detection of dark matter with standard and nonstandard interactions.

Projected axis ratios of galaxy clusters in the Horizon-AGN simulation: Impact of baryon physics and comparison with observations

Abstract We characterize the non-sphericity of galaxy clusters by the projected axis ratio of spatial distribution of star, dark matter, and X-ray surface brightness (XSB). We select 40 simulated groups and clusters of galaxies with mass larger than 5 × 1013 M⊙ from the Horizon simulation that fully incorporates the relevant baryon physics, in particular, the active galactic nucleus feedback. We find that the baryonic physics around the central region of galaxy clusters significantly affects the non-sphericity of dark matter distribution even beyond the central region, approximately up to half of the virial radius. Therefore it is very difficult to predict the probability density function (PDF) of the projected axis ratio of XSB from dark-matter-only N-body simulations as attempted in previous studies. Indeed, we find that the PDF derived from our simulated clusters exhibits much better agreement with that from the observed X-ray clusters. This indicates that our present methodology to estimate the non-sphericity directly from the Horizon simulation is useful and promising. Further improvements in both numerical modeling and observational data will establish the non-sphericity of clusters as a cosmological test complementary to more conventional statistics based on spherically averaged quantities.

Cosmic shear as a probe of galaxy formation physics

We evaluate the potential for current and future cosmic shear measurements from large galaxy surveys to constrain the impact of baryonic physics on the matter power spectrum. We do so using a model-independent parameterization that describes deviations of the matter power spectrum from the dark-matter-only case as a set of principal components that are localized in wavenumber and redshift. We perform forecasts for a variety of current and future datasets, and find that at least ~90% of the constraining power of these datasets is contained in no more than nine principal components. The constraining power of different surveys can be quantified using a figure of merit defined relative to currently available surveys. With this metric, we find that the final Dark Energy Survey dataset (DES Y5) and the Hyper Suprime Cam Survey will be roughly an order of magnitude more powerful than existing data in constraining baryonic effects. Upcoming Stage IV surveys (LSST, Euclid, and WFIRST) will improve upon this by a further factor of a few. We show that this conclusion is robust to marginalization over several key systematics. The ultimate power of cosmic shear to constrain galaxy formation is dependent on understanding systematics in the shear measurements at small (sub-arcminute) scales. If these systematics can be sufficiently controlled, cosmic shear measurements from DES Y5 and other future surveys have the potential to provide a very clean probe of galaxy formation and to strongly constrain a wide range of predictions from modern hydrodynamical simulations.

The impact of baryonic physics on the structure of dark matter haloes: the view from the FIRE cosmological simulations

We study the distribution of cold dark matter (CDM) in cosmological simulations from the FIRE (Feedback In Realistic Environments) project, for $M_{\ast}\sim10^{4-11}\,M_{\odot}$ galaxies in $M_{\rm h}\sim10^{9-12}\,M_{\odot}$ halos. FIRE incorporates explicit stellar feedback in the multi-phase ISM, with energetics from stellar population models. We find that stellar feedback, without "fine-tuned" parameters, greatly alleviates small-scale problems in CDM. Feedback causes bursts of star formation and outflows, altering the DM distribution. As a result, the inner slope of the DM halo profile ($\alpha$) shows a strong mass dependence: profiles are shallow at $M_{\rm h}\sim10^{10}-10^{11}\,M_{\odot}$ and steepen at higher/lower masses. The resulting core sizes and slopes are consistent with observations. This is broadly consistent with previous work using simpler feedback schemes, but we find steeper mass dependence of $\alpha$, and relatively late growth of cores. Because the star formation efficiency $M_{\ast}/M_{\rm h}$ is strongly halo mass dependent, a rapid change in $\alpha$ occurs around $M_{\rm h}\sim 10^{10}\,M_{\odot}$ ($M_{\ast}\sim10^{6}-10^{7}\,M_{\odot}$), as sufficient feedback energy becomes available to perturb the DM. Large cores are not established during the period of rapid growth of halos because of ongoing DM mass accumulation. Instead, cores require several bursts of star formation after the rapid buildup has completed. Stellar feedback dramatically reduces circular velocities in the inner kpc of massive dwarfs; this could be sufficient to explain the "Too Big To Fail" problem without invoking non-standard DM. Finally, feedback and baryonic contraction in Milky Way-mass halos produce DM profiles slightly shallower than the Navarro-Frenk-White profile, consistent with the normalization of the observed Tully-Fisher relation.

Counts of galaxy clusters as cosmological probes: the impact of baryonic physics

The halo mass function from N-body simulations of collisionless matter is generally used to retrieve cosmological parameters from observed counts of galaxy clusters. This neglects the observational fact that the baryonic mass fraction in clusters is a random variable that, on average, increases with the total mass (within an overdensity of 500). Considering a mock catalog that includes tens of thousands of galaxy clusters, as expected from the forthcoming generation of surveys, we show that the effect of a varying baryonic mass fraction will be observable with high statistical significance. The net effect is a change in the overall normalization of the cluster mass function and a milder modification of its shape. Our results indicate the necessity of taking into account baryonic corrections to the mass function if one wants to obtain unbiased estimates of the cosmological parameters from data of this quality. We introduce the formalism necessary to accomplish this goal. Our discussion is based on the conditional probability of finding a given value of the baryonic mass fraction for clusters of fixed total mass. Finally, we show that combining information from the cluster counts with measurements of the baryonic mass fraction in a small subsample of clusters (including only a few tens of objects) will nearly optimally constrain the cosmological parameters.

Constraining primordial non-Gaussianity with cosmological weak lensing: shear and flexion

We examine the cosmological constraining power of future large-scale weak lensing surveys on the model of the ESA planned mission Euclid, with particular reference to primordial non-Gaussianity. Our analysis considers several different estimators of the projected matter power spectrum, based on both shear and flexion. We review the covariance and Fisher matrix for cosmic shear and evaluate those for cosmic flexion and for the cross-correlation between the two. The bounds provided by cosmic shear alone are looser than previously estimated, mainly due to the reduced sky coverage and background number density of sources for the latest Euclid specifications. New constraints for the local bispectrum shape, marginalized over σ8, are at the level of ΔfNL ∼ 100, with the precise value depending on the exact multipole range that is considered in the analysis. We consider three additional bispectrum shapes, for which the cosmic shear constraints range from ΔfNL ∼ 340 (equilateral shape) up to ΔfNL ∼ 500 (orthogonal shape). Also, constraints on the level of non-Gaussianity and on the amplitude of the matter power spectrum σ8 are almost perfectly anti-correlated, except for the orthogonal bispectrum shape for which they are correlated. The competitiveness of cosmic flexion constraints against cosmic shear ones depends by and large on the galaxy intrinsic flexion noise, that is still virtually unconstrained. Adopting the very high value that has been occasionally used in the literature results in the flexion contribution being basically negligible with respect to the shear one, and for realistic configurations the former does not improve significantly the constraining power of the latter. Since the shear shot noise is white, while the flexion one decreases with decreasing scale, by considering high enough multipoles the two contributions have to become comparable. Extending the analysis up to ℓmax = 20,000 cosmic flexion, while being still subdominant, improves the shear constraints by ∼ 10% when added. However on such small scales the highly non-linear clustering of matter, the impact of baryonic physics, and the non-Gaussian part of the covariance matrix make any error estimation uncertain. By considering lower, and possibly more realistic, values of the flexion intrinsic shape noise results in flexion constraining power being a factor of ∼ 2 better than that of shear, and the bounds on σ8 and fNL being improved by a factor of ∼ 3 upon their combination.

The impact of baryonic physics on the shape and radial alignment of substructures in cosmological dark matter haloes

We use two simulations performed within the Constrained Local UniversE Simulation (CLUES) project to study both the shape and radial alignment of (the dark matter component of) subhaloes; one of the simulations is a dark matter only model while the other run includes all the relevant gas physics and star formation recipes. We find that the involvement of gas physics does not have a statistically significant effect on either property -- at least not for the most massive subhaloes considered in this study. However, we observe in both simulations including and excluding gasdynamics a (pronounced) evolution of the dark matter shapes of subhaloes as well as of the radial alignment signal since infall time. Further, this evolution is different when positioned in the central and outer regions of the host halo today; while subhaloes tend to become more aspherical in the central 50% of their host's virial radius, the radial alignment weakens in the central regime while strengthening in the outer parts. We confirm that this is due to tidal torquing and the fact that subhaloes at pericentre move too fast for the alignment signal to respond.

Impact of baryon physics on dark matter structures: a detailed simulation study of halo density profiles

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

Substructures in hydrodynamical cluster simulations

The abundance and structure of dark matter subhalos has been analyzed extensively in recent studies of dark matter-only simulations, but comparatively little is known about the impact of baryonic physics on halo substructures. We here extend the SUBFIND algorithm for substructure identification such that it can be reliably applied to dissipative hydrodynamical simulations that include star formation. This allows, in particular, the identification of galaxies as substructures in simulations of clusters of galaxies, and a determination of their content of gravitationally bound stars, dark matter, and hot and cold gas. Using a large set of cosmological cluster simulations, we present a detailed analysis of halo substructures in hydrodynamical simulations of galaxy clusters, focusing in particular on the influence both of radiative and non-radiative gas physics, and of non-standard physics such as thermal conduction and feedback by galactic outflows. We also examine the impact of numerical nuisance parameters such as artificial viscosity parameterizations. We find that diffuse hot gas is efficiently stripped from subhalos when they enter the highly pressurized cluster atmosphere. This has the effect of decreasing the subhalo mass function relative to a corresponding dark matter-only simulation. These effects are mitigated in radiative runs, where baryons condense in the central subhalo regions and form compact stellar cores. However, in all cases, only a very small fraction, of the order of one percent, of subhalos within the cluster virial radii preserve a gravitationally bound hot gaseous atmosphere. (abridged)