Showing papers on "Structure formation published in 2021"

Progress on cosmological magnetic fields.

Abstract: A variety of observations impose upper limits at the nano Gauss level on magnetic fields that are coherent on inter-galactic scales while blazar observations indicate a lower bound ∼10-16G. Such magnetic fields can play an important astrophysical role, for example at cosmic recombination and during structure formation, and also provide crucial information for particle physics in the early Universe. Magnetic fields with significant energy density could have been produced at the electroweak phase transition. The evolution and survival of magnetic fields produced on sub-horizon scales in the early Universe, however, depends on the magnetic helicity which is related to violation of symmetries in fundamental particle interactions. The generation of magnetic helicity requires new CP violating interactions that can be tested by accelerator experiments via decay channels of the Higgs particle.

Structure formation in large-volume cosmological simulations of fuzzy dark matter: impact of the non-linear dynamics

Abstract: An ultra-light bosonic particle of mass around $10^{-22}\,\mathrm{eV}/c^2$ is of special interest as a dark matter candidate, as it both has particle physics motivations, and may give rise to notable differences in the structures on highly non-linear scales due to the manifestation of quantum-physical wave effects on macroscopic scales, which could address a number of contentious small-scale tensions in the standard cosmological model, $\Lambda$CDM. Using a spectral technique, we here discuss simulations of such fuzzy dark matter (FDM), including the full non-linear wave dynamics, with a comparatively large dynamic range and for larger box sizes than considered previously. While the impact of suppressed small-scale power in the initial conditions associated with FDM has been studied before, the characteristic FDM dynamics are often neglected; in our simulations, we instead show the impact of the full non-linear dynamics on physical observables. We focus on the evolution of the matter power spectrum, give first results for the FDM halo mass function directly based on full FDM simulations, and discuss the computational challenges associated with the FDM equations. FDM shows a pronounced suppression of power on small scales relative to cold dark matter (CDM), which can be understood as a damping effect due to 'quantum pressure'. In certain regimes, however, the FDM power can exceed that of CDM, which may be interpreted as a reflection of order-unity density fluctuations occurring in FDM. In the halo mass function, FDM shows a significant abundance reduction below a characteristic mass scale only. This could in principle alleviate the need to invoke very strong feedback processes in small galaxies to reconcile $\Lambda$CDM with the observed galaxy luminosity function, but detailed studies that also include baryons will be needed to ultimately judge the viability of FDM.

Bounds on warm dark matter from Schwarzschild primordial black holes

Abstract: We consider light dark matter candidates originated from the evaporation of Schwarzschild primordial black holes, with masses in the range $$10^{-5}$$ – $$10^9$$ g. These candidates are beyond standard model particles with negligible couplings to the other particles, so that they interact only gravitationally. Belonging to the category of warm dark matter, they nevertheless spoil structure formation, with a softer impact for increasing values of the candidate spin. Requiring such candidates to fully account for the observed dark matter, we find that the scenario of black hole domination is ruled out for all spin values up to 2. For the scenario of radiation domination, we derive upper limits on the parameter $$\beta $$ (the primordial black hole energy density at formation over the radiation one), which are less stringent the higher the candidate spin is.

The Abell 3391/95 galaxy cluster system - A 15 Mpc intergalactic medium emission filament, a warm gas bridge, infalling matter clumps, and (re-) accelerated plasma discovered by combining SRG/eROSITA data with ASKAP/EMU and DECam data

Abstract: Context. Inferences about dark matter, dark energy, and the missing baryons all depend on the accuracy of our model of large-scale structure evolution. In particular, with cosmological simulations in our model of the Universe, we trace the growth of structure, and visualize the build-up of bigger structures from smaller ones and of gaseous filaments connecting galaxy clusters.Aims. Here we aim to reveal the complexity of the large-scale structure assembly process in great detail and on scales from tens of kiloparsecs up to more than 10 Mpc with new sensitive large-scale observations from the latest generation of instruments. We also aim to compare our findings with expectations from our cosmological model.Methods. We used dedicated SRG/eROSITA performance verification (PV) X-ray, ASKAP/EMU Early Science radio, and DECam optical observations of a ~15 deg2 region around the nearby interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process, the morphological complexity in the inner parts of the clusters, and the (re-)acceleration of plasma. We also used complementary Sunyaev-Zeldovich (SZ) effect data from the Planck survey and custom-made Galactic total (neutral plus molecular) hydrogen column density maps based on the HI4PI and IRAS surveys. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite.Results. We trace the irregular morphology of warm and hot gas of the main clusters from their centers out to well beyond their characteristic radii, r 200 . Between the two main cluster systems, we observe an emission bridge on large scale and with good spatial resolution. This bridge includes a known galaxy group but this can only partially explain the emission. Most gas in the bridge appears hot, but thanks to eROSITA’s unique soft response and large field of view, we discover some tantalizing hints for warm, truly primordial filamentary gas connecting the clusters. Several matter clumps physically surrounding the system are detected. For the “Northern Clump,” we provide evidence that it is falling towards A3391 from the X-ray hot gas morphology and radio lobe structure of its central AGN. Moreover, the shapes of these X-ray and radio structures appear to be formed by gas well beyond the virial radius, r 100 , of A3391, thereby providing an indirect way of probing the gas in this elusive environment. Many of the extended sources in the field detected by eROSITA are also known clusters or new clusters in the background, including a known SZ cluster at redshift z = 1. We find roughly an order of magnitude more cluster candidates than the SPT and ACT surveys together in the same area. We discover an emission filament north of the virial radius of A3391 connecting to the Northern Clump. Furthermore, the absorption-corrected eROSITA surface brightness map shows that this emission filament extends south of A3395 and beyond an extended X-ray-emitting object (the “Little Southern Clump”) towards another galaxy cluster, all at the same redshift. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA PV observation field. The Northern and Southern Filament are each detected at >4σ . The Planck SZ map additionally appears to support the presence of both new filaments. Furthermore, the DECam galaxy density map shows galaxy overdensities in the same regions. Overall, the new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution; the similarities of features found in a similar system in the Magneticum simulation are striking. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure.

Scaling Properties of Galaxy Groups

Abstract: Galaxy groups and poor clusters are more common than rich clusters, and host the largest fraction of matter content in the Universe. Hence, their studies are key to understand the gravitational and thermal evolution of the bulk of the cosmic matter. Moreover, because of their shallower gravitational potential, galaxy groups are systems where non-gravitational processes (e.g., cooling, AGN feedback, star formation) are expected to have a higher impact on the distribution of baryons, and on the general physical properties, than in more massive objects, inducing systematic departures from the expected scaling relations. Despite their paramount importance from the astrophysical and cosmological point of view, the challenges in their detection have limited the studies of galaxy groups. Upcoming large surveys will change this picture, reassigning to galaxy groups their central role in studying the structure formation and evolution in the Universe, and in measuring the cosmic baryonic content. Here, we review the recent literature on various scaling relations between X-ray and optical properties of these systems, focusing on the observational measurements, and the progress in our understanding of the deviations from the self-similar expectations on groups' scales. We discuss some of the sources of these deviations, and how feedback from supernovae and/or AGNs impacts the general properties and the reconstructed scaling laws. Finally, we discuss future prospects in the study of galaxy groups.

COSMIC BIRTH: Efficient Bayesian Inference of the Evolving Cosmic Web from Galaxy Surveys

Abstract: We present COSMIC BIRTH: COSMological Initial Conditions from Bayesian Inference Reconstructions with THeoretical models: an algorithm to reconstruct the primordial and evolved cosmic density fields from galaxy surveys on the light-cone. The displacement and peculiar velocity fields are obtained from forward modelling at different redshift snapshots given some initial cosmic density field within a Gibbs-sampling scheme. This allows us to map galaxies, observed in a light-cone, to a single high redshift and hereby provide tracers and the corresponding survey completeness in Lagrangian space including phase-space mapping. These Lagrangian tracers in turn permit us to efficiently obtain the primordial density field, making the COSMIC BIRTH code general to any structure formation model. Our tests are restricted for the time being to Augmented Lagrangian Perturbation theory. We show how to robustly compute the non-linear Lagrangian bias from clustering measurements in a numerical way, enabling us to get unbiased dark matter field reconstructions at initial cosmic times. We also show that we can greatly recover the information of the dark matter field from the galaxy distribution based on a detailed simulation. Novel key ingredients to this approach are a higher-order Hamiltonian sampling technique and a non-diagonal Hamiltonian mass-matrix. This technique could be used to study the Eulerian galaxy bias from galaxy surveys and could become an ideal baryon acoustic reconstruction technique. In summary, this method represents a general reconstruction technique, including in a self-consistent way a survey mask, non-linear and non-local bias and redshift space distortions, with an efficiency about 10 times superior to previous comparable methods.

ETHOS ‐ an effective theory of structure formation: Impact of dark acoustic oscillations on cosmic dawn

Abstract: Upcoming data of the 21-cm hydrogen line during cosmic dawn ($z\ensuremath{\sim}10--30$) will revolutionize our understanding of the astrophysics of the first galaxies. Here we present a case study on how to exploit those same measurements to learn about the nature of dark matter (DM) at small scales. Focusing on the effective theory of structure formation (ETHOS) paradigm, we run a suite of simulations covering a broad range of DM microphysics, connecting the output of $N$-body simulations to dedicated 21-cm simulations to predict the evolution of the 21-cm signal across the entire cosmic dawn. We find that observatories targeting both the global signal and the 21-cm power spectrum are sensitive to all ETHOS models we study, and can distinguish them from CDM if the suppression wave number is smaller than $k\ensuremath{\approx}300\text{ }\text{ }h/\mathrm{Mpc}$, even when accounting for feedback with a phenomenological model. This is an order of magnitude smaller comoving scales than currently constrained by other datasets, including the Lyman-$\ensuremath{\alpha}$ forest. Moreover, if a prospective 21-cm detection confirmed a deficiency of power at small scales, we show that ETHOS models with strong dark acoustic oscillations can be discriminated from the pure suppression of warm dark matter, showing the power of 21-cm data to understand the behavior of DM at the smallest physical scales.

Cannibal domination and the matter power spectrum

Abstract: Decoupled hidden sectors can easily and generically result in a period of cannibal domination, during which the dominant component of the Universe has an equation of state intermediate between radiation and matter due to self-heating by number-changing interactions. We present for the first time the consequences of a cannibal-dominated era prior to big bang nucleosynthesis for structure formation on small scales. We find that an early cannibal-dominated era imprints a characteristic peak on the dark matter power spectrum, with scale and amplitude directly determined by the mass, lifetime, and number-changing interaction strength of the cannibal field. This enhancement to the small-scale matter power spectrum will generate early-forming dark matter microhalos, and we provide a detailed and transparent map between the properties of the cannibal species and the characteristic mass and formation time of these structures. These relations demonstrate how the internal workings of a hidden sector leave a potentially observable imprint on the matter power spectrum even if dark matter has no direct couplings to the Standard Model.

Self-interacting neutrino portal dark matter

Abstract: Author(s): Lamprea, JM; Peinado, E; Smolenski, S; Wudka, J | Abstract: We present a realistic, simple and natural model of strongly-interacting dark matter based on the neutrino-portal paradigm. The strong interactions at small velocities are generated by the exchange of dark photons, and produce the observed core-like DM distribution in galactic centers; this effect could be spoiled by the formation of DM bound states (also due to dark-photon effects), which we avoid by requiring the DM candidates to be light, with masses below O(10 GeV). The mixing of the dark photon with the Z and ordinary photon is strongly suppressed by introducing a softly-broken discrete symmetry similar to charge conjugation, which also ensures that the dark photon life-time is short enough to avoid restrictions derived form big-bang nucleosynthesis and large-scale structure formation. Other constraints are accommodated without the need of fine tuning, in particular nucleon scattering occurs only at one loop, so direct detection cross sections are naturally suppressed. Neutrino masses are generated through the inverse see saw.

Euclid: Impact of non-linear and baryonic feedback prescriptions on cosmological parameter estimation from weak lensing cosmic shear

Abstract: Upcoming surveys will map the growth of large-scale structure with unprecented precision, improving our understanding of the dark sector of the Universe. Unfortunately, much of the cosmological information is encoded on small scales, where the clustering of dark matter and the effects of astrophysical feedback processes are not fully understood. This can bias the estimates of cosmological parameters, which we study here for a joint analysis of mock Euclid cosmic shear and Planck cosmic microwave background data. We use different implementations for the modelling of the signal on small scales and find that they result in significantly different predictions. Moreover, the different non-linear corrections lead to biased parameter estimates, especially when the analysis is extended into the highly non-linear regime, with the Hubble constant, H0, and the clustering amplitude, σ8, affected the most. Improvements in the modelling of non-linear scales will therefore be needed if we are to resolve the current tension with more and better data. For a given prescription for the non-linear power spectrum, using different corrections for baryon physics does not significantly impact the precision of Euclid, but neglecting these correction does lead to large biases in the cosmological parameters. In order to extract precise and unbiased constraints on cosmological parameters from Euclid cosmic shear data, it is therefore essential to improve the accuracy of the recipes that account for non-linear structure formation, as well as the modelling of the impact of astrophysical processes that redistribute the baryons.Key words: gravitational lensing: weak / large-scale structure of Universe / cosmological parameters⋆ This paper is published on behalf of the Euclid Consortium.

Scalar dark matter candidates revisited

Abstract: We revisit the possibility of light-scalar dark matter, in the MeV to GeV mass bracket and coupled to electrons through fermion or vector mediators, in light of significant experimental and observational advances that probe new physics below the GeV scale. We establish new limits from electron colliders and fixed-target beams and derive the strength of loop-induced processes that are probed by precision physics, among other laboratory probes. In addition, we compute the cooling bound from SN1987A, consider self-scattering, structure formation, and cosmological constraints as well as the limits from dark matter--electron scattering in direct detection experiments. We then show that the combination of constraints largely excludes the possibility that the galactic annihilation of these particles may explain the long-standing International Gamma-Ray Astrophysics Laboratory excess of 511 keV photons as observed in the Galactic bulge. As a caveat to these conclusions, we identify the resonant annihilation regime where the vector mediator goes nearly on shell.

Bayesian forward modelling of cosmic shear data

Abstract: We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. We test and validate our framework with realistic simulated shear data, demonstrating that the method recovers the unbiased matter distribution and the correct lensing and matter power spectrum. While the cosmology is fixed in this test, and the method employs a prior power spectrum, we demonstrate that the lensing results are sensitive to the true power spectrum when this differs from the prior. In this case, the density field samples are generated with a power spectrum that deviates from the prior, and the method recovers the true lensing power spectrum. The method also recovers the matter power spectrum across the sky, but as currently implemented, it cannot determine the radial power since isotropy is not imposed. In summary, our method provides physically plausible inference of the dark matter distribution from cosmic shear data, allowing us to extract information beyond the two-point statistics and exploiting the full information content of the cosmological fields.

Gravitational Wave and CMB Probes of Axion Kination

Abstract: Rotations of an axion field in field space provide a natural origin for an era of kination domination, where the energy density is dominated by the kinetic term of the axion field, preceded by an early era of matter domination. Remarkably, no entropy is produced at the end of matter domination and hence these eras of matter and kination domination may occur even after Big Bang Nucleosynthesis. We derive constraints on these eras from both the cosmic microwave background and Big Bang Nucleosynthesis. We investigate how this cosmological scenario affects the spectrum of possible primordial gravitational waves and find that the spectrum features a triangular peak. We discuss how future observations of gravitational waves can probe the viable parameter space, including regions that produce axion dark matter by the kinetic misalignment mechanism or the baryon asymmetry by axiogenesis. For QCD axion dark matter produced by the kinetic misalignment mechanism, a modification to the inflationary gravitational wave spectrum occurs above 0.01 Hz and, for high values of the energy scale of inflation, the prospects for discovery are good. We briefly comment on implications for structure formation of the universe.

Formation of inflaton halos after inflation

Abstract: The early Universe may have passed through an extended period of matter-dominated expansion following inflation and prior to the onset of radiation domination. Subhorizon density perturbations grow gravitationally during such an epoch, collapsing into bound structures if it lasts long enough. The strong analogy between this phase and structure formation in the present-day Universe allows the use of $N$-body simulations and approximate methods for halo formation to model the fragmentation of the inflaton condensate into inflaton halos. For a simple model we find that these halos have masses of up to 20 kg and radii of the order of ${10}^{\ensuremath{-}20}\text{ }\text{ }\mathrm{m}$, roughly ${10}^{\ensuremath{-}24}$ seconds after the big bang. We find that the $N$-body halo mass function matches predictions of the mass-peak patch method and the Press-Schechter formalism within the expected range of scales. A long matter-dominated phase would imply that reheating and thermalization occurs in a universe with large variations in density, potentially modifying the dynamics of this process. In addition, large overdensities can source gravitational waves and may lead to the formation of primordial black holes.

On the formation and stability of fermionic dark matter haloes in a cosmological framework

Abstract: The formation and stability of collisionless self-gravitating systems is a long standing problem, which dates back to the work of D. Lynden-Bell on violent relaxation, and extends to the issue of virialization of dark matter (DM) halos. An important prediction of such a relaxation process is that spherical equilibrium states can be described by a Fermi-Dirac phase-space distribution, when the extremization of a coarse-grained entropy is reached. In the case of DM fermions, the most general solution develops a degenerate compact core surrounded by a diluted halo. As shown recently, the latter is able to explain the galaxy rotation curves while the DM core can mimic the central black hole. A yet open problem is whether this kind of astrophysical core-halo configurations can form at all, and if they remain stable within cosmological timescales. We assess these issues by performing a thermodynamic stability analysis in the microcanonical ensemble for solutions with given particle number at halo virialization in a cosmological framework. For the first time we demonstrate that the above core-halo DM profiles are stable (i.e. maxima of entropy) and extremely long lived. We find the existence of a critical point at the onset of instability of the core-halo solutions, where the fermion-core collapses towards a supermassive black hole. For particle masses in the keV range, the core-collapse can only occur for $M_{\rm vir} \gtrsim 10^9 M_\odot$ starting at $z_{\rm vir}\approx 10$ in the given cosmological framework. Our results prove that DM halos with a core-halo morphology are a very plausible outcome within nonlinear stages of structure formation.

Observational constraints on dark matter scattering with electrons

Abstract: We present new observational constraints on the elastic scattering of dark matter with electrons for dark matter masses between 10 keV and 1 TeV. We consider scenarios in which the momentum-transfer cross section has a power-law dependence on the relative particle velocity, with a power-law index $n \in \{-4,-2,0,2,4,6\}$. We search for evidence of dark matter scattering through its suppression of structure formation. Measurements of the cosmic microwave background temperature, polarization, and lensing anisotropy from \textit{Planck} 2018 data and of the Milky Way satellite abundance measurements from the Dark Energy Survey and Pan-STARRS1 show no evidence of interactions. We use these data sets to obtain upper limits on the scattering cross section, comparing them with exclusion bounds from electronic recoil data in direct detection experiments. Our results provide the strongest bounds available for dark matter--electron scattering derived from the distribution of matter in the Universe, extending down to sub-MeV dark matter masses, where current direct detection experiments lose sensitivity.

The origin of metal-poor stars on prograde disc orbits in FIRE simulations of Milky Way-mass galaxies

Abstract: Author(s): Santistevan, IB; Wetzel, A; Sanderson, RE; El-Badry, K; Samuel, J; Faucher-Giguere, CA | Abstract: In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disc. However, recent observations of stars in the MW show that metal-poor ([Fe/H] ≲ -2) stars are preferentially on prograde orbits with respect to the disc. Using the Feedback In Realistic Environments 2 (FIRE-2) suite of cosmological zoom-in simulations of MW/M31-mass galaxies, we investigate the prevalence and origin of prograde metal-poor stars. Almost all (11 of 12) of our simulations have metal-poor stars on preferentially prograde orbits today and throughout most of their history: we thus predict that this is a generic feature of MW/M31-mass galaxies. The typical prograde-to-retrograde ratio is ~2:1, which depends weakly on stellar metallicity at [Fe/H] ~ -1. These trends predicted by our simulations agree well with MW observations. Prograde metal-poor stars originate largely from a single Large/Small Magellanic Cloud (LMC/SMC)-mass gas-rich merger 7-12.5 Gyr ago, which deposited existing metal-poor stars and significant gas on an orbital vector that sparked the formation of and/or shaped the orientation of a long-lived stellar disc, giving rise to a prograde bias for all low-metallicity stars. We find subdominant contributions from in situ stars formed in the host galaxy before this merger, and in some cases, additional massive mergers. We find few clear correlations between any properties of our MW/M31-mass galaxies at z = 0 and the degree of this prograde bias as a result of diverse merger scenarios.

Primordial black hole formation with full numerical relativity

Abstract: We study the formation of black holes from subhorizon and superhorizon perturbations in a matter dominated universe with 3+1D numerical relativity simulations. We find that there are two primary mechanisms of formation depending on the initial perturbation's mass and geometry -- via $\textit{direct collapse}$ of the initial overdensity and via $\textit{post-collapse accretion}$ of the ambient dark matter. In particular, for the latter case, the initial perturbation does not have to satisfy the hoop conjecture for a black hole to form. In both cases, the duration of the formation the process is around a Hubble time, and the initial mass of the black hole is $M_\mathrm{BH} \sim 10^{-2} H^{-1} M_\mathrm{Pl}^2$. Post formation, we find that the PBH undergoes rapid mass growth beyond the self-similar limit $M_\mathrm{BH}\propto H^{-1}$, at least initially. We argue that this implies that most of the final mass of the PBH is accreted from its ambient surroundings post formation.

Maximum Absorption of the Global 21 cm Spectrum in the Standard Cosmological Model

Abstract: The absorption feature in the global spectrum is likely the first 21cm observable from the cosmic dawn, which provides valuable insights into the earliest history of structure formation. We run a set of high-resolution hydrodynamic simulations of early structure formation to assess the effect of non-linear structure formation on the maximum absorption level (i.e. assuming the spin temperature coupling is saturated) of the global 21 cm spectrum in the standard cosmological framework. We ignore the star formation and feedbacks, which also tends to reduce the absorption signal, but take into account the inevitable non-linear density fluctuations in the intergalactic medium (IGM), shock heating and Compton heating which can reduce the absorption level. We found that the combination of these reduced the maximum absorption signal by $\sim 15\%$ at redshift 17, as compared with the homogeneous or linearly fluctuating IGM. These effects have to be carefully accounted for when interpreting the observational results, especially when considering the necessity of introducing new physics.

Gravitational wave induced baryon acoustic oscillations

Abstract: We study the impact of gravitational waves originating from a first order phase transition on structure formation. To do so, we perform a second order perturbation analysis in the $1+3$ covariant framework and derive a wave equation in which second order, adiabatic density perturbations of the photon-baryon fluid are sourced by the gravitational wave energy density during radiation domination and on sub-horizon scales. The scale on which such waves affect the energy density perturbation spectrum is found to be proportional to the horizon size at the time of the phase transition times its inverse duration. Consequently, structure of the size of galaxies and bigger can only be affected in this way by relatively late phase transitions at $\ge 10^{6}\,\text{s}$. Using cosmic variance as a bound we derive limits on the strength $\alpha$ and the relative duration $(\beta/H_*)^{-1}$ of phase transitions as functions of the time of their occurrence which results in a new exclusion region for the energy density in gravitational waves today. We find that the cosmic variance bound forbids only relative long lasting phase transitions, e.g. $\beta/H_*\lesssim 6.8$ for $t_*\approx 5\times10^{11}\,\text{s}$, which exhibit a substantial amount of supercooling $\alpha>20$ to affect the matter power spectrum.

Setting the stage: structures from Gaussian random fields

Abstract: We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Λ cold dark matter, ΛCDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z = 0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrization of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics, and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localized in real space, our method is a powerful tool to study structure formation in cosmological simulations.

Kinetic field theory: Non-linear cosmic power spectra in the mean-field approximation

Abstract: We use the recently developed Kinetic Field Theory (KFT) for cosmic structure formation to show how non-linear power spectra for cosmic density fluctuations can be calculated in a mean-field approximation to the particle interactions. Our main result is a simple, closed and analytic, approximate expression for this power spectrum. This expression has two parameters characterising non-linear structure growth which can be calibrated within KFT itself. Using this self-calibration, the non-linear power spectrum agrees with results obtained from numerical simulations to within typically $\lesssim10\,\%$ up to wave numbers $k\lesssim10\,h\,\mathrm{Mpc}^{-1}$ at redshift $z = 0$. Adjusting the two parameters to optimise agreement with numerical simulations, the relative difference to numerical results shrinks to typically $\lesssim 5\,\%$. As part of the derivation of our mean-field approximation, we show that the effective interaction potential between dark-matter particles relative to Zel'dovich trajectories is sourced by non-linear cosmic density fluctuations only, and is approximately of Yukawa rather than Newtonian shape.

Vector modes in ΛCDM: the gravitomagnetic potential in dark matter haloes from relativistic N-body simulations

Abstract: We investigate the transverse modes of the gravitational and velocity fields in Λ cold dark matter, based on a high-resolution simulation performed using the adaptive-mesh refinement general-relativistic N-body code GRAMSES. We study the generation of vorticity in the dark matter velocity field at low redshift, providing fits to the shape and evolution of its power spectrum over a range of scales. By analysing the gravitomagnetic vector potential, that is absent in Newtonian simulations, in dark matter haloes with masses ranging from ∼1012.5 to ∼1015 h−1 M⊙, we find that its magnitude correlates with the halo mass, peaking in the inner regions. Nevertheless, on average, its ratio against the scalar gravitational potential remains fairly constant, below percent level, decreasing roughly linearly with redshift and showing a weak dependence on halo mass. Furthermore, we show that the gravitomagnetic acceleration in haloes peaks towards the core and reaches almost 10−10h cm s−2 in the most massive halo of the simulation. However, regardless of the halo mass, the ratio between the gravitomagnetic force and the standard gravitational force is typically at around the 10−5 level inside the haloes, again without significant radius dependence. This result confirms that the gravitomagnetic effects have negligible impact on structure formation, even for the most massive structures, although its behaviour in low-density regions remains to be explored. Likewise, the impact on observations remains to be understood in the future.

Gravitational Lensing in Cosmology

Abstract: Gravitational lensing is a unique and direct probe of mass in the universe. It depends only on the law of gravity and does not depend on the dynamical state nor the composition of matter. Thus, it is used to study the distribution of the dark matter in the lensing object. Combined with the traditional observations such as optical and X-ray, it gives us useful informations of the structure formation in the universe. The lensing observables depend also on the global geometry as well as large scale structure of the universe. Therefore it is possible to withdraw useful constraints on the cosmological parameters once the distribution of lensing mass is accurately known. Since the first discovery of the lensing event by a galaxy in 1979, various kinds of lensing phenomena caused by star, galaxy, cluster of galaxies and large scale structure have been observed and are used to study mass distribution in various scales and cosmology. Thus, the gravitational lensing is now regarded as an indispensable research field in the observational cosmology. In this paper, we give an instructive introduction to gravitational lensing and its applications to cosmology.

Persistent homology of the cosmic web – I. Hierarchical topology in ΛCDM cosmologies

Abstract: Using a set of $\Lambda$CDM simulations of cosmic structure formation, we study the evolving connectivity and changing topological structure of the cosmic web using state-of-the-art tools of multiscale topological data analysis (TDA). We follow the development of the cosmic web topology in terms of the evolution of Betti number curves and feature persistence diagrams of the three (topological) classes of structural features: matter concentrations, filaments and tunnels, and voids. The Betti curves specify the prominence of features as a function of density level, and their evolution with cosmic epoch reflects the changing network connections between these structural features. The persistence diagrams quantify the longevity and stability of topological features. In this study we establish, for the first time, the link between persistence diagrams, the features they show, and the gravitationally driven cosmic structure formation process. By following the diagrams' development over cosmic time, the link between the multiscale topology of the cosmic web and the hierarchical buildup of cosmic structure is established. The sharp apexes in the diagrams are intimately related to key transitions in the structure formation process. The apex in the matter concentration diagrams coincides with the density level at which, typically, they detach from the Hubble expansion and begin to collapse. At that level many individual islands merge to form the network of the cosmic web and a large number of filaments and tunnels emerge to establish its connecting bridges. The location trends of the apex possess a self-similar character that can be related to the cosmic web's hierarchical buildup. We find that persistence diagrams provide a significantly higher and more profound level of information on the structure formation process than more global summary statistics like Euler characteristic or Betti numbers.

Sommerfeld-enhanced dark matter searches with dwarf spheroidal galaxies

Abstract: We study observable signals from dark matter that self-annihilates via Sommerfeld effect in dwarf spheroidal galaxies (dSphs). Since the effect of the Sommerfeld enhancement depends on the velocity of dark matter, it is crucial to determine the profile of dSphs to compute the J-factor, i.e., the line-of-sight integral of density squared. In our study we use the prior distributions of the parameters for satellite density profiles in order to determine the J-factor, making most out of the recent developments in the N-body simulations and semi-analytical modeling for the structure formation. As concrete model, we analyze fermionic dark matter that annihilates via a light scalar and Wino dark matter in supersymmetric models. We find that, with the more realistic prior distributions that we adopt in this study, the J-factor of the most promising dwarf galaxies is decreased by a factor of a few, compared with earlier estimates based on non-informative priors. Nevertheless, the Cherenkov Telescope Array should be able to detect thermal Wino dark matter by pointing it toward best classical or ultrafaint dwarf galaxies for 500 hours.

A new mechanism for primordial black hole formation during reheating.

Abstract: The Reheating process at the end of inflation is often modeled by an oscillating scalar field which shows a background dust-like behaviour, prompting the analysis of gravitational collapse and black hole formation in this era to be approached by the spherical collapse of standard structure formation. In the scalar field dark matter structure formation process virialized halos halt the direct collapse, resulting in halos with condensed central cores at the de Broglie scale of the dominant scalar field. We show that a similar process can take place during reheating, leading to the formation of primordial black holes (PBHs). We study the formation of PBHs through the gravitational further collapse of structures virialized during reheating, looking at the collapse of either the whole structure, or that of the central core within these configurations. We compute the threshold amplitude for the density contrast to undergo this process, for both free and self-interacting scalar fields. We discuss the relevance of our results for the abundance of PBHs at the lower end of the mass spectrum.

The anisotropy of the power spectrum in periodic cosmological simulations

Abstract: The classical gravitational force on a torus is anisotropic and always lower than Newton's $1/r^2$ law We demonstrate the effects of periodicity in dark matter only $N$-body simulations of spherical collapse and standard $\Lambda$CDM initial conditions Periodic boundary conditions cause an overall negative and anisotropic bias in cosmological simulations of cosmic structure formation The effect is most significant when the largest mildly non-linear scales are comparable to the linear size of the simulation box, as often is the case for high-resolution hydrodynamical simulations Spherical collapse morphs into a shape similar to an octahedron The anisotropic growth distorts the large-scale $\Lambda$CDM dark matter structures We introduce the direction-dependent power spectrum invariant under the octahedral group of the simulation volume and show that the results break spherical symmetry

Neutrino as the dark force

Abstract: We point out a novel role for the Standard Model neutrino in dark matter phenomenology where the exchange of neutrinos generates a long-range potential between dark matter particles. The resulting dark matter self interaction could be sufficiently strong to impact small-scale structure formation, without the need of any light dark force carrier. This is a generic feature of theories where dark matter couples to the visible sector through the neutrino portal. It is highly testable with improved decay rate measurements at future $Z$, Higgs, and $\ensuremath{\tau}$ factories, as well as precision cosmology.

Frozen-in fermionic singlet dark matter in non-standard cosmology with a decaying fluid

Abstract: We perform a detailed study of dark matter production via freeze-in under the assumption that some fluid dominates the early Universe before depositing its energy to the plasma causing entropy injection. As a dark matter candidate we consider a fermionic singlet that is produced through its interactions with a scalar particle in the thermal plasma. The fluid alters the expansion rate of the Universe, as well as the scaling of the temperature, which significantly affects the evolution of both the number density and the mean momentum of the dark matter particle. We identify and discuss in detail the effects of the evolution of these quantities by considering several examples representing dark matter production at different stages of expansion and entropy injection. We find that, since the dark matter density is reduced when the entropy injection to the plasma continues after freeze-in, in order to reproduce its observational value an enhanced rate of dark matter production is required relative to standard cosmology. Furthermore, the impact of the assumed non-standard cosmological history on the dark matter mean momentum can result in either a relaxed or a tightened bound on the dark matter mass from large structure formation data.