Cold Dark Matter Density Research Articles

A dark matter profile to model diverse feedback-induced core sizes of ΛCDM haloes

ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation.

CMB distance priors revisited: effects of dark energy dynamics, spatial curvature, primordial power spectrum, and neutrino parameters

As a physical and sufficient compression of the full CMB data, the CMB distance priors, or shift parameters, have been widely used and provide a convenient way to include CMB data when obtaining cosmological constraints. In this paper, we revisit this data vector and examine its stability under different cosmological models. We find that the CMB distance priors are an accurate substitute for the full CMB data when probing dark energy dynamics. This is true when the primordial power spectrum model is directly generalized from the power spectrum of the model used in the derivation of the distance priors from the CMB data. We discover a difference when a non-flat model with the untilted primordial inflation power spectrum is used to measure the distance priors. This power spectrum is a radical change from the more conventional tilted primordial power spectrum and violates fundamental assumptions for the reliability of the CMB shift parameters. We also investigate the performance of CMB distance priors when the sum of neutrino masses ∑ mν and the effective number of relativistic species Neff are allowed to vary. Our findings are consistent with earlier results: the neutrino parameters can change the measurement of the sound horizon from CMB data, and thus the CMB distance priors. We find that when the neutrino model is allowed to vary, the cold dark matter density ωc and Neff need to be included in the set of parameters that summarize CMB data, in order to reproduce the constraints from the full CMB data. We present an updated and expanded set of CMB distance priors which can reproduce constraints from the full CMB data within 1σ, and are applicable to models with massive neutrinos, as well as non-standard cosmologies.

Model of metastable EeV dark matter

We propose a model where a long-lived pseudoscalar EeV particle can be produced with sufficient abundance so as to account for the cold dark matter density, despite having a Planck mass suppressed coupling to the thermal bath. Connecting this state to a hidden sterile neutrino sector through derivative couplings, induced by higher dimensional operators, allows one to account for light neutrino masses while having a lifetime that can be much larger than the age of the Universe. Moreover, the same derivative coupling accounts for the production of dark matter in the very first instant of the reheating. Given the sensitivity of the IceCube and ANITA collaborations, we study the possible signatures of such a model in the form of Ultra-High-Energy Cosmic Rays in the neutrino sector, and show that such signals could be detected in the near future.

Spherical collapse in coupled quintessence with a ΛCDM background

In this work we study the growth of cold dark matter density perturbations in the nonlinear regime on a conformally coupled quintessence model in which the background is designed to mimic a $\Lambda$CDM cosmology. The spherical collapse of overdense regions is analyzed. We highlight the role of the coupling on the overall dynamics, trace the evolution of the density contrast throughout the cosmic history and compute perturbative parameters such as the critical density contrast. We find that the coupling has the influence of delaying the collapse due to the slower growth of matter perturbations. We follow to compute the cluster number counts using the Press-Schechter and Sheth-Tormen mass functions. In both cases, the transfer of energy between the dark energy field and dark matter suppresses the number of objects at low redshifts and enhances the number at high redshifts. Finally, we compute the expected cluster number counts for the future eROSITA mission and the current South Pole Telescope survey.

Superstring-inspired particle cosmology: inflation, neutrino masses, leptogenesis, dark matter & the SUSY scale

We develop a string-inspired model for particle cosmology, based on a flipped SU(5)×U(1) gauge group formulated in a no-scale supergravity framework. The model realizes Starobinsky-like inflation, which we assume to be followed by strong reheating, with the GUT symmetry being broken subsequently by a light `flaton' field whose decay generates a second stage of reheating. We discuss the production of gravitinos and the non-thermal contribution made by their decays to the density of cold dark matter, which is assumed to be provided by the lightest neutralino. We also discuss the masses of light and heavy neutrinos and leptogenesis. As discussed previously [1], a key rôle is played by a superpotential coupling between the inflaton, matter and GUT Higgs fields, called λ6. We scan over possible values of λ6, exploring the correlations between the possible values of observables. We emphasize that the release of entropy during the GUT transition allows large regions of supersymmetry-breaking parameter space that would otherwise lead to severe overdensity of dark matter. Furthermore, we find that the Big Bang nucleosynthesis lower limit on the reheating temperature of ∼ 1 MeV restricts the supersymmetry-breaking scale to a range \U0001d4aa(10) TeV that is consistent with the absence of supersymmetric particles at the LHC.

Coupled vector dark energy

We provide a general framework for studying the evolution of background and cosmological perturbations in the presence of a vector field Aμ coupled to cold dark matter (CDM) . We consider an interacting Lagrangian of the form Q f(X) Tc, where Q is a coupling constant, f is an arbitrary function of X=−AμAμ/2, and Tc is a trace of the CDM energy-momentum tensor. The matter coupling affects the no-ghost condition and sound speed of linear scalar perturbations deep inside the sound horizon, while those of tensor and vector perturbations are not subject to modifications. The existence of interactions also modifies the no-ghost condition of CDM density perturbations. We propose a concrete model of coupled vector dark energy with the tensor propagation speed equivalent to that of light. In comparison to the Q=0 case, we show that the decay of CDM to the vector field leads to the phantom dark energy equation of state wDE closer to −1. This alleviates the problem of observational incompatibility of uncoupled models in which wDE significantly deviates from −1. The maximum values of wDE reached during the matter era are bounded from the CDM no-ghost condition of future de Sitter solutions.

Secondary CMB anisotropies from magnetized haloes

Magnetized plasmas within haloes of galaxies leave their footprint on the polarized anisotropies of the cosmic microwave background. The two dominant effects of astrophysical haloes are Faraday rotation, which generates rotation of the plane of linear polarization, and Faraday conversion, which induces a leakage from linear polarization to circular polarization. We revisit these sources of secondary anisotropies by computing the angular power spectra of the Faraday rotation angle and the Faraday conversion rate by the large-scale structures. To this end, we use the halo model and we pay special attention to the impact of magnetic field projections. Assuming magnetic fields of haloes to be uncorrelated, we found a vanishing two-halo term, and angular power spectra peaking at multipoles ℓ ∼ 104. The Faraday rotation angle is dominated by the contribution of thermal electrons. For the Faraday conversion rate, we found that both thermal electrons and relativistic, non-thermal electrons contribute equally in the most optimistic case for the density and Lorentz factor of relativistic electrons, while in more pessimistic cases the thermal electrons give the dominant contribution. Assuming the magnetic field to be independent of the halo mass, the angular power spectra for both effects roughly scale with the amplitude of matter perturbations as ∼σ38, and with a very mild dependence with the density of cold dark matter. Introducing a dependence of the magnetic field strength with the halo mass leads to an increase of the scaling at large angular scales (above a degree) with the amplitude of matter fluctuations up to ∼σ9.58 for Faraday rotation and ∼σ158 for Faraday conversion for a magnetic field strength scaling linearly with the halo mass. Introducing higher values of the magnetic field for galaxies, as compared to clusters, instead leads to a decrease of such a scaling at arcminute scales down to ∼σ0.98 for Faraday rotation.

Forecasts of cosmological constraints from Type Ia supernovae including the weak-lensing convergence

We investigate how the cosmological constraints from SNe Ia are improved by including the effects of weak-lensing convergence. To do so, we introduce the lognormal function as the convergence PDF modeling the lensing scatter of SN Ia magnitude, and apply a sample selection for SNe Ia to avoid strongly lensed samples. Comparing with the contribution of other uncertainties (e.g., the intrinsic magnitude scatter), we find that the lensing effect is dominant at z > 1. Then forecasting the parameter constraints for the Wide-Field InfraRed Survey Telescope survey, we show that considering the weak-lensing effect, the constraints on the density parameters Ωm or ΩΛ, and the dark energy equation of state w are improved, especially for SNe Ia samples at higher redshift z > 1. Furthermore, we see that the degeneracy between the total mass of neutrino Σ mν and the (cold) dark matter density parameter Ωc can be resolved.

Null test for interactions in the dark sector

Since there is no known symmetry in Nature that prevents a non-minimal coupling between the dark energy (DE) and cold dark matter (CDM) components, such a possibility constitutes an alternative to standard cosmology, with its theoretical and observational consequences being of great interest. In this paper we propose a new null test on the standard evolution of the dark sector based on the time dependence of the ratio between the CDM and DE energy densities which, in the standard $\Lambda$CDM scenario, scales necessarily as $a^{-3}$. We use the latest measurements of type Ia supernovae, cosmic chronometers and angular baryonic acoustic oscillations to reconstruct the expansion history using model-independent Machine Learning techniques, namely, the Linear Model formalism and Gaussian Processes. We find that while the standard evolution is consistent with the data at $3\sigma$ level, some deviations from the $\Lambda$CDM model are found at low redshifts, which may be associated with the current tension between local and global determinations of $H_0$.

Accurate fitting functions for peculiar velocity spectra in standard and massive-neutrino cosmologies

We estimate the velocity field in a large set of N-body simulations including massive neutrino particles, and measure the auto-power spectrum of the velocity divergence field as well as the cross-power spectrum between the cold dark matter density and the velocity divergence. We perform these measurements at four different redshifts and within four different cosmological scenarios, covering a wide range in neutrino masses. We find that the nonlinear correction to the velocity power spectra largely depends on the degree of nonlinear evolution with no specific dependence on the value of neutrino mass. We provide a fitting formula based on the value of the rms of the matter fluctuations in spheres of 8 h−1 Mpc, describing the nonlinear corrections with 3% accuracy on scales below k = 0.7 h Mpc−1.

Cosmological constraints on parametrized interacting dark energy

We reconsider the dynamics of the Universe in the presence of interactions in the cosmological dark sector. A class of interacting models is introduced via a real function fr of the ratio r between the energy densities of the (pressureless) cold dark matter (CDM) and dark energy (DE). The subclass of models for which the ratio r depends only on the scale factor is shown to be equivalent to unified models of the dark sector, i.e. models for which the CDM and DE components can be combined in order to form a unified dark fluid. For specific choices of the function fr we recover several models already studied in the literature. We analyze various special cases of this type of interacting models using a suitably modified version of the CLASS code combined with Monte Python in order to constrain the parameter space with the data from supernova of type SNe Ia (JLA), the Hubble constant H0, cosmic chronometers (CC), baryon acoustic oscillations (BAO) and data from the Planck satellite (Planck TT). Our analysis shows that even if data from the late Universe ( H0, SNe Ia and CC) indicate an interaction in the dark sector, the data related to the early Universe (BAO and Planck TT) constrain this interaction substantially, in particular for cases in which the background dynamics is strongly affected.

Observational Constraints on the Tilted Spatially Flat and the Untilted Nonflat ϕCDM Dynamical Dark Energy Inflation Models

Abstract We constrain spatially flat tilted and nonflat untilted scalar field (ϕ) dynamical dark energy inflation (ϕCDM) models by using Planck 2015 cosmic microwave background (CMB) anisotropy measurements and recent baryonic acoustic oscillation distance observations, Type Ia supernovae apparent magnitude data, Hubble parameter measurements, and growth rate data. We assume an inverse power-law scalar field potential energy density . We find that the combination of the CMB data with the four non-CMB data sets significantly improves parameter constraints and strengthens the evidence for nonflatness in the nonflat untilted ϕCDM case from 1.8σ for the CMB measurements only to more than 3.1σ for the combined data. In the nonflat untilted ϕCDM model, current observations favor a spatially closed universe with spatial curvature contributing about two-thirds of a percent of the present cosmological energy budget. The flat tilted ϕCDM model is a 0.4σ better fit to the data than is the standard flat tilted ΛCDM model: current data allow for the possibility that dark energy is dynamical. The nonflat tilted ϕCDM model is in better accord with the Dark Energy Survey bounds on the rms amplitude of mass fluctuations now (σ 8) as a function of the nonrelativistic matter density parameter now (Ω m ) but it does not provide as good a fit to the larger-multipole Planck 2015 CMB anisotropy data as does the standard flat tilted ΛCDM model. A few cosmological parameter value measurements differ significantly when determined using the tilted flat and the untilted nonflat ϕCDM models, including the cold dark matter density parameter and the reionization optical depth.

Dark matter casts light on the early Universe

We show how knowledge of the cold dark matter (CDM) density can be used, in conjunction with measurements of the parameters of a scenario for beyond the Standard Model (BSM) physics, to provide information about the evolution of the Universe before Big Bang Nucleosynthesis (BBN). As examples of non-standard evolution, we consider models with a scalar field that may decay into BSM particles, and quintessence models. We illustrate our calculations using various supersymmetric models as representatives of classes of BSM scenarios in which the CDM density is either larger or smaller than the observed density when the early Universe is assumed to be radiation-dominated. In the case of a decaying scalar field, we show how the CDM density can constrain the initial scalar density and the reheating temperature after it decays in BSM scenarios that would yield overdense dark matter in standard radiation-dominated cosmology, and how the decays of the scalar field into BSM particles can be constrained in scenarios that would otherwise yield underdense CDM. We also show how the early evolution of the quintessence field can be constrained in BSM scenarios.

Lensing reconstruction in post-Born cosmic microwave background weak lensing

The study of the Cosmic Microwave Background (CMB) lensing potential has established itself by now as a robust way of probing the physics of large-scale structure growth. The most common estimators of the lensing potential are derived under the assumption of Gaussianity of the matter distribution and in the Born approximation of the photon diffusion. In this paper we study the performance of quadratic estimators when applied to realistic sky maps extracted from multiple-lens ray tracing techniques in cosmological $N$-body simulations. These are expected to model accurately the effects due to the non-Gaussianity of the matter distribution induced by its nonlinearity and the deviation from the Born approximation. We show that both these effects on their own lead to reconstruction biases, but these tend to partially cancel each other when both these effects are considered together. We forecast the impact of these biases on the estimation of cosmological parameters for future high-sensitivity CMB experiments like CMB-S4. We find that the cold dark matter density, $\Omega_\textrm{cdm}$, the optical depth to reionization $\tau$, the amplitude of primordial inflationary perturbations, $A_s$ and the sum of neutrino masses, $M_\nu$, could be biased at the 1-2$\sigma$ level, if no external data set is used. We also observe a reduction of the bias if external data like baryon acoustic oscillations (BAO) is included.

Primordial black holes with an accurate QCD equation of state

Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately 0.7M⊙, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately 0.3M⊙<M<1.4M⊙ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, 10−3M⊙<MPBH<103M⊙, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs.

“Invisible” QCD axion rolling through the QCD phase transition

Visible matter in the current Universe is a consequence of the phase transition of the strong force, quantum chromodynamics (QCD). This phase transition has occurred at the Universe temperature around Tc≃165 MeV while it was expanding. Strongly interacting matter particles are quarks above Tc, while they are pions, protons and neutrons below Tc. The spin degrees of freedom 37 (u and d quarks and gluons) just above Tc are converted to 3 (pions) after the phase transition. This phase transition might have been achieved mostly at supercooled temperatures. The supercooling was provided by the expansion of the Universe. We obtain the effective bubble formation rate α(T)≈104−5 MeV and the completion temperature of the phase change (to the hadronic phase), Tf≃126 MeV. During the phase transition, the scale factor R has increased by a factor of 2.4. This provides a key knowledge on the energy density of “invisible” QCD axion at the full hadronic-phase commencement temperature Tf, and allows for us to estimate the current energy density of cold dark matter composed of “invisible” QCD axions.

Scale-dependent bias and bispectrum in neutrino separate universe simulations

Cosmic background neutrinos have a large velocity dispersion, which causes the evolution of long-wavelength density perturbations to depend on scale. This scale-dependent growth leads to the well-known suppression in the linear theory matter power spectrum that is used to probe neutrino mass. In this paper, we study the impact of long-wavelength density perturbations on small-scale structure formation. By performing separate universe simulations where the long-wavelength mode is absorbed into the local expansion, we measure the responses of the cold dark matter (CDM) power spectrum and halo mass function, which correspond to the squeezed-limit bispectrum and halo bias. We find that the scale-dependent evolution of the long-wavelength modes causes these quantities to depend on scale and provide simple expressions to model them in terms of scale and the amount of massive neutrinos. Importantly, this scale-dependent bias reduces the suppression in the linear halo power spectrum due to massive neutrinos by 13 and 26% for objects of bias $\bar{b}=2$ and $\bar{b} \gg1$, respectively. We demonstrate with high statistical significance that the scale-dependent halo bias ${\it cannot}$ be modeled by the CDM and neutrino density transfer functions at the time when the halos are identified. This reinforces the importance of the temporal nonlocality of structure formation, especially when the growth is scale dependent.

Interactive mixture of inhomogeneous dark fluids driven by dark energy: a dynamical system analysis

We examine the evolution of an inhomogeneous mixture of non-relativistic pressureless cold dark matter (CDM), coupled to dark energy (DE) characterised by the equation of state parameter w<-1/3, with the interaction term proportional to the DE density. This coupled mixture is the source of a spherically symmetric Lemaître–Tolman–Bondi (LTB) metric admitting an asymptotic Friedman–Lemaître–Robertson–Walker (FLRW) background. Einstein’s equations reduce to a 5-dimensional autonomous dynamical system involving quasi–local variables related to suitable averages of covariant scalars and their fluctuations. The phase space evolution around the critical points (past/future attractors and five saddles) is examined in detail. For all parameter values and both directions of energy flow (CDM to DE and DE to CDM) the phase space trajectories are compatible with a physically plausible early cosmic times behaviour near the past attractor. This result compares favourably with mixtures with interaction driven by the CDM density, whose past evolution is unphysical for DE to CDM energy flow. Numerical examples are provided describing the evolution of an initial profile that can be associated with idealised structure formation scenarios.

Likelihood analysis of the sub-GUT MSSM in light of LHC 13-TeV data

We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale M_mathrm{in} below the supersymmetric grand unification scale M_mathrm{GUT}, as can occur in mirage mediation and other models. In addition to M_mathrm{in}, such ‘sub-GUT’ models have the 4 parameters of the CMSSM, namely a common gaugino mass m_{1/2}, a common soft supersymmetry-breaking scalar mass m_0, a common trilinear mixing parameter A and the ratio of MSSM Higgs vevs tan beta , assuming that the Higgs mixing parameter mu > 0. We take into account constraints on strongly- and electroweakly-interacting sparticles from sim 36/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for M_mathrm{in}sim 10^5 to 10^9 ,, mathrm {GeV}, with M_mathrm{in}sim M_mathrm{GUT} disfavoured by Delta chi ^2 sim 3 due to the mathrm{BR}(B_{s, d} rightarrow mu ^+mu ^-) constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with m_{tilde{chi }^0_{1}} sim 1 ,, mathrm {TeV}, with annihilation via heavy Higgs bosons H / A and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the (g-2)_mu constraint, finding similar results whether or not it is included.

A Comparison of Cosmological Parameters Determined from CMB Temperature Power Spectra from the South Pole Telescope and the Planck Satellite

Abstract The Planck cosmic microwave background temperature data are best fit with a ΛCDM model that mildly contradicts constraints from other cosmological probes. The South Pole Telescope (SPT) 2540 SPT-SZ survey offers measurements on sub-degree angular scales (multipoles ) with sufficient precision to use as an independent check of the Planck data. Here we build on the recent joint analysis of the SPT-SZ and Planck data in Hou et al. by comparing ΛCDM parameter estimates using the temperature power spectrum from both data sets in the SPT-SZ survey region. We also restrict the multipole range used in parameter fitting to focus on modes measured well by both SPT and Planck, thereby greatly reducing sample variance as a driver of parameter differences and creating a stringent test for systematic errors. We find no evidence of systematic errors from these tests. When we expand the maximum multipole of SPT data used, we see low-significance shifts in the angular scale of the sound horizon and the physical baryon and cold dark matter densities, with a resulting trend to higher Hubble constant. When we compare SPT and Planck data on the SPT-SZ sky patch to Planck full-sky data but keep the multipole range restricted, we find differences in the parameters n s and . We perform further checks, investigating instrumental effects and modeling assumptions, and we find no evidence that the effects investigated are responsible for any of the parameter shifts. Taken together, these tests reveal no evidence for systematic errors in SPT or Planck data in the overlapping sky coverage and multipole range and at most weak evidence for a breakdown of ΛCDM or systematic errors influencing either the Planck data outside the SPT-SZ survey area or the SPT data at .