Dark Matter Overdensities Research Articles

Dark matter, black holes, and gravitational waves

The formation and growth of black holes can strongly influence the distribution of dark matter around them. I discuss here the different types of dark matter overdensities around black holes, including dark matter cusps, spikes, mounds, crests, and gravitational atoms. I then review recent results on the evolution of a black holes binary in presence of dark matter, focusing on the energy transfer between binary and dark matter induced by dynamical friction. Finally, I present the prospects for studying dark matter with gravitational wave observations, and argue that future interferometers might be able to detect and characterize dark matter overdensities around black holes.

The Supersonic Project: Lighting Up the Faint End of the JWST UV Luminosity Function

The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe, when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation (z ≳ 10). We study low-mass galaxy formation, including this “stream velocity,” using high-resolution AREPO hydrodynamics simulations and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function down to extremely faint and low-mass galaxies (M UV ≳ −15, 104 M ⊙ ≤ M * ≤ 108 M ⊙). We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint end of the UVLF at z = 12. We demonstrate that JWST observations are close to this enhanced regime and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories.

Measuring dark matter spikes around primordial black holes with Einstein Telescope and Cosmic Explorer

Future ground-based gravitational wave observatories will be ideal probes of the environments surrounding black holes with masses $1\ensuremath{-}10{M}_{\ensuremath{\bigodot}}$. Binary black hole mergers with mass ratios of order $q={m}_{2}/{m}_{1}\ensuremath{\lesssim}{10}^{\ensuremath{-}3}$ can remain in the frequency band of such detectors for months or years, enabling precision searches for modifications of their gravitational waveforms with respect to vacuum inspirals. As a concrete example of an environmental effect, we consider here a population of binary primordial black holes which are expected to be embedded in dense cold dark matter spikes. We provide a viable formation scenario for these systems compatible with all observational constraints and predict upper and lower limits on the merger rates of small-mass-ratio pairs. Given a detected signal of one such system by either Einstein Telescope or Cosmic Explorer, we show that the properties of the binary and of the dark matter spike can be measured to excellent precision with one week's worth of data, if the effect of the dark matter spike on the waveform is taken into account. However, we show that there is a risk of biased parameter inference or missing the events entirely if the effect of the predicted dark matter overdensity around these objects is not properly accounted for.

Comparing accretion disks and dark matter spikes in intermediate mass ratio inspirals

Intermediate mass ratio inspirals (IMRIs) will be observable with space-based gravitational wave detectors such as the Llaser interferometer space antenna. To this end, the environmental effects in such systems have to be modeled and understood. These effects can include (baryonic) accretion disks and dark matter (DM) overdensities, so-called spikes. For the first time, we model an IMRI system with both an accretion disk and a DM spike present and compare their effects on the inspiral and the emitted gravitational wave signal. We study the eccentricity evolution, employ the braking index and derive the dephasing index, which turn out to be complementary observational signatures. They allow us to disentangle the accretion disk and DM spike effects in the IMRI system and can be utilized to study environmental effects in general.

Iterative mean-field approach to the spherical collapse of dark matter haloes

ABSTRACT Gravitational collapse of dark matter overdensities leads to the formation of dark matter haloes which embed galaxies and galaxy clusters. An intriguing feature of dark matter haloes is that their density profiles closely follow a universal form irrespective of the initial condition or the corresponding growth history. This represents a class of dynamical systems with emergent universalities. We propose an ‘iterative mean-field approach’ to compute the solutions of the gravitational collapse dynamics. This approach iteratively searches for the evolution of the interaction field ϕ(t) – in this case the enclosed mass profile M(r, t) – that is consistent with the dynamics, thus that ϕ(t) is the fix-point of the iterative mapping, $\mathcal {H}(\phi) = \phi$. The formalism replaces the N-body interactions with one-body interactions with the coarse-grained interaction field, and thus shares the spirit of the mean-field theory in statistical physics. This ‘iterative mean-field approach’ combines the versatility of numerical simulations and the comprehensiveness of analytical solutions, and is particularly powerful in searching for and understanding intermediate asymptotic states in a wide range of dynamical systems where the solutions can not be obtained through the traditional self-similar analysis.

Structure formation in dark matter particle production cosmology

We investigate a cosmological scenario in which the dark matter particles can be created during the evolution of the Universe. By regarding the Universe as an open thermodynamic system and using non-equilibrium thermodynamics, we examine the mechanism of gravitational particle production. In this setup, we study the large-scale structure (LSS) formation of the Universe in the Newtonian regime of perturbations and derive the equations governing the evolution of the dark matter overdensities. Then, we implement the cosmological data from Planck 2018 CMB measurements, SNe Ia and BAO observations, as well as the Riess et al. (2019) local measurement for H0 to provide some cosmological constraints for the parameters of our model. We see that the best case of our scenario (χtot2=3834.40) fits the observational data better than the baseline ΛCDM model (χtot2=3838.00) at the background level. We moreover estimate the growth factor of linear perturbations and show that the best case of our model (χfσ82=39.85) fits the LSS data significantly better than the ΛCDM model (χfσ82=45.29). Consequently, our model also makes a better performance at the level of the linear perturbations compared to the standard cosmological model. Although the improvement in χtot for Model 2 compared to ΛCDM at the background level may arise somewhat from the additional degree of freedom, the better consistency of this model at the level of linear perturbations results mainly from the impact of particle production on the sound speed at low redshifts.

A universal profile for stacked filaments from cold dark matter simulations

ABSTRACT We study the stacked filaments connecting group-mass halo pairs, using dark-matter-only N-body simulations. We calculate the dark matter overdensity profile of these stacked filaments at different redshifts as a function of the distance perpendicular to the filament axis. A four-parameter universal functional form, including three comoving scale radii and one amplitude parameter (core density), provides a good fit out to a radius of 20h−1cMpc for stacked filaments over a range of redshifts, lengths, and masses. The scale radii are approximately independent of redshift but increase as power laws with the comoving filament length. Lastly, we compare the scaling of the filament mass measured directly from the simulations to the predicted scaling from the halo–halo–matter three-point correlation function as a function of redshift and of the mass of the halo pairs. We find that both measured scalings are similar to, but somewhat shallower than the predictions, by 10 per cent and 30 per cent, respectively. These results provide a template to interpret present and upcoming observational results based on stacking, for example, weak lensing, thermal and kinetic Sunyaev-Zel’dovich, or X-ray observations.

LyMAS reloaded: improving the predictions of the large-scale Lyman-α forest statistics from dark matter density and velocity fields

ABSTRACT We present LyMAS2, an improved version of the ‘Lyman-α Mass Association Scheme’ aiming at predicting the large-scale 3D clustering statistics of the Lyman-α forest (Ly α) from moderate-resolution simulations of the dark matter (DM) distribution, with prior calibrations from high-resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN suite simulations, (100 Mpc h)−3 comoving volume, using Wiener filtering, combining information from DM density and velocity fields (i.e. velocity dispersion, vorticity, line-of-sight 1D-divergence and 3D-divergence). All new predictions have been done at z = 2.5 in redshift space, while considering the spectral resolution of the SDSS-III BOSS Survey and different DM smoothing (0.3, 0.5, and 1.0 Mpc h−1 comoving). We have tried different combinations of DM fields and found that LyMAS2, applied to the Horizon-noAGN DM fields, significantly improves the predictions of the Ly α 3D clustering statistics, especially when the DM overdensity is associated with the velocity dispersion or the vorticity fields. Compared to the hydrodynamical simulation trends, the two-point correlation functions of pseudo-spectra generated with LyMAS2 can be recovered with relative differences of ∼5 per cent even for high angles, the flux 1D power spectrum (along the light of sight) with ∼2 per cent and the flux 1D probability distribution function exactly. Finally, we have produced several large mock BOSS spectra (1.0 and 1.5 Gpc h−1) expected to lead to much more reliable and accurate theoretical predictions.

Circularization versus eccentrification in intermediate mass ratio inspirals inside dark matter spikes

Inspirals of an intermediate mass black hole (IMBH) and a solar mass type object will be observable by space based gravitational wave detectors such as the Laser Interferometer Space Antenna. A dark matter overdensity around an IMBH---a dark matter spike---can affect the orbital evolution of the system. We consider here such intermediate mass ratio inspirals on eccentric orbits, experiencing dynamical friction of the dark matter spike. We find that by including the relative velocities of the dark matter particles, the dynamical friction tends to circularize the orbit, in contrast to previous inquiries. We derive a general condition for circularization or eccentrification for any given dissipative force. In addition to the dephasing, we suggest using the circularization rate as another probe of the dark matter spike. Observing these effects would be an indicator for the particle nature of dark matter.

Measuring the dark matter environments of black hole binaries with gravitational waves

Large dark matter overdensities can form around black holes of astrophysical and primordial origin as they form and grow. This "dark dress" inevitably affects the dynamical evolution of binary systems, and induces a dephasing in the gravitational waveform that can be probed with future interferometers. In this paper, we introduce a new analytical model to rapidly compute gravitational waveforms in presence of an evolving dark matter distribution. We then present a Bayesian analysis determining when dressed black hole binaries can be distinguished from GR-in-vacuum ones and how well their parameters can be measured, along with how close they must be to be detectable by the planned Laser Interferometer Space Antenna (LISA). We show that LISA can definitively distinguish dark dresses from standard binaries and characterize the dark matter environments around astrophysical and primordial black holes for a wide range of model parameters. Our approach can be generalized to assess the prospects for detecting, classifying, and characterizing other environmental effects in gravitational wave physics.

Investigating cosmological GAN emulators using latent space interpolation

ABSTRACT Generative adversarial networks (GANs) have been recently applied as a novel emulation technique for large-scale structure simulations. Recent results show that GANs can be used as a fast and efficient emulator for producing novel weak lensing convergence maps as well as cosmic web data in 2D and 3D. However, like any algorithm, the GAN approach comes with a set of limitations, such as an unstable training procedure, inherent randomness of the produced outputs, and difficulties when training the algorithm on multiple data sets. In this work, we employ a number of techniques commonly used in the machine learning literature to address the mentioned limitations. Specifically, we train a GAN to produce weak lensing convergence maps and dark matter overdensity field data for multiple redshifts, cosmological parameters, and modified gravity models. In addition, we train a GAN using the newest Illustris data to emulate dark matter, gas, and internal energy distribution data simultaneously. Finally, we apply the technique of latent space interpolation as a tool for understanding the feature space of the GAN algorithm. We show that the latent space interpolation procedure allows the generation of outputs with intermediate cosmological parameters that were not included in the training data. Our results indicate a 1–20 per cent difference between the power spectra of the GAN-produced and the test data samples depending on the data set used and whether Gaussian smoothing was applied. Similarly, the Minkowski functional analysis indicates a good agreement between the emulated and the real images for most of the studied data sets.

The Coupling of Galactic Dark Matter Halos with Stellar Bars

Abstract Resonant torques couple stellar bars to dark matter halos. Here we use high-resolution numerical simulations to demonstrate long-term angular momentum transfer between stellar bars and dark matter orbits of varying orientation. We show that bar-driven reversals of dark matter orbit orientations can play a surprisingly large role in the evolution of the bar pattern speed. In predominantly prograde (corotating) halos, dark matter orbits become trapped in the stellar bar, forming a parallel dark matter bar. This dark matter bar reaches more than double the vertical height of the stellar bar. In halos dominated by retrograde orbits, a dark matter wake forms oriented perpendicular to the stellar bar. These dark matter overdensities provide a novel space to look for dark matter annihilation or decay signals. We predict that the Milky Way hosts a dark matter bar aligned with the stellar bar as well as a dark matter wake, the near side of which should extend from the Galactic center to a galactic longitude of l ≈ 323°.

Detecting dark matter around black holes with gravitational waves: Effects of dark-matter dynamics on the gravitational waveform

A dark matter overdensity around a black hole may significantly alter the dynamics of the black hole's merger with another compact object. We consider here intermediate mass-ratio inspirals of stellar-mass compact objects with intermediate-mass black holes "dressed" with dark matter. We first demonstrate that previous estimates based on a fixed dark-matter dress are unphysical for a range of binaries and dark-matter distributions by showing that the total energy dissipated by the compact object through dynamical friction, as it inspirals through the dense dark matter environment towards the black hole, is larger than the gravitational binding energy of the dark-matter dress itself. We then introduce a new formalism that allows us to self-consistently follow the evolution of the dark-matter dress due to its gravitational interaction with the binary. We show that the dephasing of the gravitational waveform induced by dark matter is smaller than previously thought, but is still potentially detectable with the LISA space interferometer. The gravitational waves from such binaries could provide powerful diagnostics of the particle nature of dark matter.

Clustering of Lyα Emitters around Quasars at z ∼ 4*

Abstract The strong observed clustering of z > 3.5 quasars indicates that they are hosted by massive ( ) dark matter halos. Assuming that quasars and galaxies trace the same large-scale structures, this should manifest as strong clustering of galaxies around quasars. Previous works on high-redshift quasar environments have failed to find convincing evidence for these overdensities. Here we conduct a survey for Lyα emitters (LAEs) in the environs of 17 quasars at z ∼ 4 probing scales of . We measure an average LAE overdensity of , which we quantify by fitting the quasar–LAE cross-correlation function. We find consistency with a power-law shape with correlation length for a fixed slope of γ = 1.8 and rule out a zero clustering hypothesis at the 95% confidence level. We also measure the LAE autocorrelation length and find cMpc (γ = 1.8), which is times higher than the value measured in blank fields. Taken together, our results clearly indicate that LAEs are significantly clustered around z ∼ 4 quasars. We compare the observed clustering with the expectation from a deterministic bias model, whereby LAEs and quasars probe the same underlying dark matter overdensities, and find that our measurements fall short of the predicted overdensities by a factor of . We discuss possible explanations for this discrepancy, including large-scale quenching or the presence of excess dust in galaxies near quasars. Finally, the large cosmic variance from field to field observed in our sample (10/17 fields are actually underdense) cautions one from overinterpreting studies of z ∼ 6 quasar environments based on a single or handful of quasar fields.

Comparing approximate methods for mock catalogues and covariance matrices II: power spectrum multipoles

ABSTRACT We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constraints related to growth rate of structure, Alcock–Paczynski distortions, and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE–COLA, pinocchio, and peakpatch), methods that rely on halo assignment schemes on to dark matter overdensities calibrated with a target N-body run (halogen, patchy), and two standard assumptions about the full density probability distribution function (Gaussian and lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5{{\ \rm per\ cent}}$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10{{\ \rm per\ cent}}$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10{{\ \rm per\ cent}}$ of those from the N-body-based covariance.

Formation of First Galaxies inside Density Peaks and Voids under the Influence of Dark Matter–Baryon Streaming Velocity. I. Initial Condition and Simulation Scheme

Abstract We present a systematic study of the cosmic variance that existed in the formation of the first stars and galaxies. We focus on the cosmic variance induced by the large-scale density and velocity environment engraved at the epoch of recombination. The density environment is predominantly determined by the dark-matter overdensity, and the velocity environment by the dark matter–baryon streaming velocity. Toward this end, we introduce a new cosmological initial condition generator, BCCOMICS, which solves the quasi-linear evolution of small-scale perturbations in the large-scale density and streaming-velocity environment and generates the initial condition for dark matter and baryons, as either particles or grid data at a specific redshift. We also describe a scheme to simulate the formation of the first galaxies inside density peaks and voids, where a local environment is treated as a separate universe. The resulting cosmic variance in the number density of minihalos and the amount of cooling mass are presented as an application. Density peaks become a site for enhanced formation of the first galaxies, which compete with the negative effect from the dark matter–baryon streaming velocity on structure formation.

No Evidence for Millimeter Continuum Source Overdensities in the Environments of z ≳ 6 Quasars

Abstract Bright high-redshift quasars (z > 6) hosting supermassive black holes (M BH > 108 M ⊙) are expected to reside in massive host galaxies embedded within some of the earliest and most massive galaxy overdensities. We analyze 1.2 mm ALMA dust continuum maps of 35 bright quasars at 6 < z < 7 and search the primary beam for excess dust continuum emission from sources with L IR ≳ 1012 L ⊙ as evidence for early protoclusters. We compare the detection rates of continuum sources at ≥5σ significance in the fields surrounding the quasars (A eff = 4.3 arcmin2) with millimeter number counts in blank field surveys. We discover 15 mm sources in the fields excluding the quasars themselves, corresponding to an overdensity of δ gal ≡ (N gal − N exp)/N exp = −0.07 ± 0.56, consistent with no detected overdensity of dusty galaxies within 140 physical kpc of the quasars. However, the apparent lack of continuum overdensity does not negate the hypothesis that quasars live in overdense environments, as evidenced by strong [C ii] overdensities found on the same scales as similarly selected quasars. The small field of view of ALMA could miss a true overdensity if it exists on scales larger than 1 cMpc, if the quasar is not centered in the overdensity, or if quasar feedback plays a role close to the quasar, but it is most likely that the large line-of-sight volume probed by a continuum survey will wash out a true overdensity signal. We discuss the necessary factors in determining the bias with which dusty star-forming galaxies trace true dark matter overdensities in order to improve upon overdensity searches in the dust continuum.

A massive core for a cluster of galaxies at a redshift of 4.3.

Massive galaxy clusters have been found that date to times as early as three billion years after the Big Bang, containing stars that formed at even earlier epochs1-3. The high-redshift progenitors of these galaxy clusters-termed 'protoclusters'-can be identified in cosmological simulations that have the highest overdensities (greater-than-average densities) of dark matter4-6. Protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts 7 . However, recent detections of possible protoclusters hosting such starbursts8-11 do not support the kind of rapid cluster-core formation expected from simulations 12 : the structures observed contain only a handful of starbursting galaxies spread throughout a broad region, with poor evidence for eventual collapse into a protocluster. Here we report observations of carbon monoxide and ionized carbon emission from the source SPT2349-56. We find that this source consists of at least 14 gas-rich galaxies, all lying at redshifts of 4.31. We demonstrate that each of these galaxies is forming stars between 50 and 1,000 times more quickly than our own Milky Way, and that all are located within a projected region that is only around 130 kiloparsecs in diameter. This galaxy surface density is more than ten times the average blank-field value (integrated over all redshifts), and more than 1,000 times the average field volume density. The velocity dispersion (approximately 410 kilometres per second) of these galaxies and the enormous gas and star-formation densities suggest that this system represents the core of a cluster of galaxies that was already at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 could be building one of the most massive structures in the Universe today.

Galaxy bias from galaxy–galaxy lensing in the DES science verification data

We present a measurement of galaxy-galaxy lensing around a magnitude-limited ($i_{AB} < 22.5$) sample of galaxies from the Dark Energy Survey Science Verification (DES-SV) data. We split these lenses into three photometric-redshift bins from 0.2 to 0.8, and determine the product of the galaxy bias $b$ and cross-correlation coefficient between the galaxy and dark matter overdensity fields $r$ in each bin, using scales above 4 Mpc/$h$ comoving, where we find the linear bias model to be valid given our current uncertainties. We compare our galaxy bias results from galaxy-galaxy lensing with those obtained from galaxy clustering (Crocce et al. 2016) and CMB lensing (Giannantonio et al. 2016) for the same sample of galaxies, and find our measurements to be in good agreement with those in Crocce et al. (2016), while, in the lowest redshift bin ($z\sim0.3$), they show some tension with the findings in Giannantonio et al. (2016). We measure $b\cdot r$ to be $0.87\pm 0.11$, $1.12 \pm 0.16$ and $1.24\pm 0.23$, respectively for the three redshift bins of width $\Delta z = 0.2$ in the range $0.2<z <0.8$, defined with the photometric-redshift algorithm BPZ. Using a different code to split the lens sample, TPZ, leads to changes in the measured biases at the 10-20\% level, but it does not alter the main conclusion of this work: when comparing with Crocce et al. (2016) we do not find strong evidence for a cross-correlation parameter significantly below one in this galaxy sample, except possibly at the lowest redshift bin ($z\sim 0.3$), where we find $r = 0.71 \pm 0.11$ when using TPZ, and $0.83 \pm 0.12$ with BPZ.

NO OVERDENSITY OF LYMAN-ALPHA EMITTING GALAXIES AROUND A QUASAR AT z ∼ 5.7

ABSTRACT Bright quasars, observed when the universe was less than one billion years old (z &gt; 5.5), are known to host massive black holes (∼109 M ⊙) and are thought to reside in the center of massive dark matter overdensities. In this picture, overdensities of galaxies are expected around high-redshift quasars. However, observations based on the detection of Lyman-break galaxies (LBGs) around these quasars do not offer a clear picture: this may be due to the uncertain redshift constraints of LBGs, which are solely selected through broadband filters. To circumvent such uncertainties, we here perform a search for Lyman-alpha emitting galaxies (LAEs) in the field of the quasar PSO J215.1512–16.0417 at z ∼ 5.73, through narrowband deep imaging with FORS2 at the Very Large Telescope. We study an area of 37 arcmin2, i.e., ∼206 comoving Mpc2 at the redshift of the quasar. We find no evidence of an overdensity of LAEs in the quasar field with respect to blank-field studies. Possible explanations for these findings may be that our survey volume is too small, or that the strong ionizing radiation from the quasar hinders galaxy formation in its immediate proximity. Another possibility is that these quasars are not situated in the dense environments predicted by some simulations.