Evolution Of Dark Matter Research Articles

Novel structures and collapse of solitons in nonminimally gravitating dark matter halos

Ultralight dark matter simulations predict condensates with short-range correlation, known as solitons or boson stars, at the centers of dark matter halos. This paper investigates the formation and collapse of dark matter solitons influenced by nonminimal gravitational effects, characterized by gradient-dependent self-interactions of dark matter and an additional source in Poisson's equation for gravity. Our simulations suggest that the initial evolution of dark matter resembles that without nonminimal gravitational effects. However, regions with negative potential curvature may develop, and solitons will collapse when their densities reach certain critical values for both positive and negative coupling constants. With strong nonminimal gravitational effects, we verify that linear density perturbations could grow on both large and small scales, potentially enhancing structure formation.

Overview and public data release of the augmented Auriga Project: cosmological simulations of dwarf and Milky Way-mass galaxies

ABSTRACT We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical ‘zoom-in’ simulations of 40 Milky Way-mass haloes and 26 dwarf galaxy–mass haloes run with the moving-mesh code arepo. Auriga adopts the Lambda cold dark matter cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully–Fisher relation, the star-forming (SF) main sequence, and H i gas fraction/disc thickness. Finally, we show that SF gas discs appear to build rotation and velocity dispersion rapidly for $z\gtrsim 3$ before they ‘settle’ into ever-increasing rotation-dispersion ratios ($V/\sigma$). This evolution appears to be in rough agreement with some kinematic measurements from H$\alpha$ observations, and demonstrates an application of how to utilize the released data.

A solar investigation of multicomponent dark matter

If multiple thermal weakly interacting massive particle (WIMP) dark matter candidates exist, then their capture and annihilation dynamics inside a massive star such as Sun could change from conventional method of study. With a simple correction to time evolution of dark matter (DM) number abundance inside the Sun for multiple dark matter candidates, significant changes in DM annihilation flux depending on annihilation, direct detection cross-section, internal conversion and their contribution to relic abundance are reported in present work.

Dark matter in a bi-metric universe

We study the possibility to describe dark matter in a model of the universe with two-scale factors and a nonstandard Poisson bracket structure characterized by the deformation parameter [Formula: see text]. The dark matter evolution is analyzed in the early stages of the universe, and its relic density is obtained via the Freeze-In and Freeze-Out mechanism. We show that by fixing [Formula: see text] and the initial ratio of energy densities present in the different sectors of the universe, the space of thermal average annihilation cross-sections and dark matter masses compatible with the standard cosmology prior to Big Bang Nucleosynthesis (BBN), is enlarged. This feature of the model is compatible with nonstandard cosmology.

Field-based physical inference from peculiar velocity tracers

ABSTRACT We present a proof-of-concept Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. Using a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the late-time density and velocity fields as natural byproducts. We implement this field-based physical inference approach by adapting the Bayesian Origin Reconstruction from Galaxies ($\small {\rm BORG}$) algorithm, which explores the high-dimensional posterior through the use of Hamiltonian Monte Carlo sampling. We test the self-consistency of the method using random sets of tracers, and assess its accuracy in a more complex scenario where peculiar velocity tracers are mock haloes drawn from $\small {\rm GADGET2}$ N-body simulations. We find that our framework self-consistently infers the initial conditions, density and velocity fields, and shows some robustness to model mis-specification. Compared with the approach of constrained Gaussian random fields/Wiener filtering, the hierarchical model produces more accurate final density and velocity field reconstructions. It also allows us to constrain the initial conditions by peculiar velocity observations, complementing in this aspect other field-based approaches based on alternative cosmological observables such as galaxy clustering or weak lensing.

Qualitative study of Lyra cosmologies with spatial curvature

We investigate the spatially curved Friedmann–Robertson–Walker space–time with Lyra geometry model. The cosmological solutions are studied by using dynamical system methodology to explore the universe evolution in the model. The constructed autonomous system has been explored to extract the physically viable cosmological solutions by studying the behavior of eigenvalues at the fixed points of the system. Late-time attractor solution corresponding to the accelerated expansion exist in the model and thus, the model may explain the late-time accelerating expansion of universe. The model is also composed of deceleration to acceleration transition with Λ cold dark matter evolution at late-times. The evolution of cosmological parameters like deceleration parameter , effective equation of state parameter, total matter density and role of energy conditions with scale factor have been studied along with statefinder diagnostic to explain the universe evolution in the model. By dynamical system analysis, it has been observed that the displacement field plays a dominant role during the early times and, the universe is dominated by stiff matter like component. Implications of bouncing behavior in the model with role of energy conditions and continuity equation have been discussed in detail along with their physical consequences. • Lyra geometry model with FRW metric having spatial curvature is investigated. • Dynamical system method has been used to extract the cosmic dynamics of model. • The model explains late-time accelerating evolution of universe. • Cosmological parameters are studied by integrating the dynamical variables. • We illustrate that the fixed point corresponding to radiation phase is absent.

Exploring the High-redshift PBH-ΛCDM Universe: Early Black Hole Seeding, the First Stars and Cosmic Radiation Backgrounds

Abstract We explore the observational implications of a model in which primordial black holes (PBHs) with a broad birth mass function ranging in mass from a fraction of a solar mass to ∼106 M ⊙, consistent with current observational limits, constitute the dark matter (DM) component in the universe. The formation and evolution of dark matter and baryonic matter in this PBH-Λ cold dark matter (ΛCDM) universe are presented. In this picture, PBH-DM mini-halos collapse earlier than in standard ΛCDM, baryons cool to form stars at z ∼ 15–20, and growing PBHs at these early epochs start to accrete through Bondi capture. The volume emissivity of these sources peaks at z ∼ 20 and rapidly fades at lower redshifts. As a consequence, PBH DM could also provide a channel to make early black hole seeds and naturally account for the origin of an underlying DM halo–host galaxy and central black hole connection that manifests as the M bh–σ correlation. To estimate the luminosity function and contribution to integrated emission power spectrum from these high-redshift PBH-DM halos, we develop a halo occupation distribution model. In addition to tracing the star formation and reionization history, it permits us to evaluate the cosmic infrared and X-ray backgrounds. We find that accretion onto PBHs/active galactic nuclei successfully accounts for the detected backgrounds and their cross-correlation, with the inclusion of an additional IR stellar emission component. Detection of the deep IR source count distribution by the James Webb Space Telescope could reveal the existence of this population of high-redshift star-forming and accreting PBH DM.

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N-body mock challenge for the eBOSS emission line galaxy sample

ABSTRACT Cosmological growth can be measured in the redshift space clustering of galaxies targeted by spectroscopic surveys. Accurate prediction of clustering of galaxies will require understanding galaxy physics, which is a very hard and highly non-linear problem. Approximate models of redshift space distortion (RSD) take a perturbative approach to solve the evolution of dark matter and galaxies in the universe. In this paper, we focus on extended Baryon Oscillation Spectroscopic (eBOSS) emission line galaxies (ELGs) that live in intermediate mass haloes. We create a series of mock catalogues using haloes from the Multidark and outer rim dark matter only N-body simulations. Our mock catalogues include various effects inspired by baryonic physics such as assembly bias and the characteristics of satellite galaxies kinematics, dynamics, and statistics deviating from dark matter particles. We analyse these mocks using the TNS RSD model in Fourier space and the convolution Lagrangian perturbation theory (CLPT) in configuration space. We conclude that these two RSD models provide an unbiased measurement of RSD within the statistical error of our mocks. We obtain the conservative theoretical systematic uncertainty of $3.3{{\ \rm per\ cent}}$, $1.8{{\ \rm per\ cent}}$, and $1.5{{\ \rm per\ cent}}$ in fσ8, α∥, and α⊥, respectively, for the TNS and CLPT models. We note that the estimated theoretical systematic error is an order of magnitude smaller than the statistical error of the eBOSS ELG sample and hence are negligible for the purpose of the current eBOSS ELG analysis.

On the phase-space structure of galaxy clusters from cosmological simulations

ABSTRACT Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) haloes and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion σv with the density ρ can be used to define the pseudo-entropy $S(r)=\sigma _\mathrm{v}^2/\rho ^{\, 2/3}$, whose profile is well described by a simple power law $S\propto \, r^{\, \alpha }$. We analyse a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyse four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profiles as traced by substructures, while DM particles profiles are characterized by an almost universal behaviour, thus suggesting that the level of pseudo-entropy could represent a potential low-scatter mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light and the stars bound to the brightest cluster galaxy: the combination of these two pseudo-entropy profiles is well described by a single power law out to almost the entire cluster virial radius.

The dark matter bispectrum from effective viscosity and one-particle irreducible vertices

Dark matter evolution during the process of cosmological structure formation can be described in terms of a one-particle irreducible effective action at a characteristic scale km and a loop expansion below this scale, based on the effective propagators and vertices. We calculate the form of the effective vertices and compute the bispectrum of density perturbations within a one-loop approximation. We find that the effective vertices play a subdominant role as compared to the effective viscosity and sound velocity that modify the (inverse) propagators. For the bispectrum we reproduce the results of standard perturbation theory in the range where it is applicable, and find a slightly improved agreement with N-body simulations at larger wavenumbers.

Real scalar dark matter: relativistic treatment

A stable real scalar provides one of the simplest possibilities to account for dark matter. We consider the regime where its coupling to the Standard Model fields is negligibly small. Due to self-coupling, the scalar field can reach thermal or at least kinetic equilibrium, in which case the system is characterized by its temperature and effective chemical potential. We perform a fully relativistic analysis of dark matter evolution, thermalization conditions and different freeze-out regimes, including the chemical potential effects. To this end, we derive a relativistic Bose-Einstein analog of the Gelmini-Gondolo formula for a thermal averaged cross section. Finally, we perform a comprehensive parameter space analysis to determine regions consistent with observational constraints. Dark matter can be both warm and cold in this model.

The IllustrisTNG simulations: public data release

We present the full public release of all data from the TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive black holes from early time to the present day, z=0. Each of the flagship runs—TNG50, TNG100, and TNG300—are accompanied by halo/subhalo catalogs, merger trees, lower-resolution and dark-matter only counterparts, all available with 100 snapshots. We discuss scientific and numerical cautions and caveats relevant when using TNG.The data volume now directly accessible online is ∼750 TB, including 1200 full volume snapshots and ∼80,000 high time-resolution subbox snapshots. This will increase to ∼1.1 PB with the future release of TNG50. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform enabling user computation with local access to TNG data, alleviating the need to download large datasets.

Resolution convergence in cosmological hydrodynamical simulations using adaptive mesh refinement

We have explored the evolution of gas distributions from cosmological simulations carried out using the RAMSES adaptive mesh refinement (AMR) code, to explore the effects of resolution on cosmological hydrodynamical simulations. It is vital to understand the effect of both the resolution of initial conditions and the final resolution of the simulation. Lower initial resolution simulations tend to produce smaller numbers of low mass structures. This will strongly affect the assembly history of objects, and has the same effect of simulating different cosmologies. The resolution of initial conditions is an important factor in simulations, even with a fixed maximum spatial resolution. The power spectrum of gas in simulations using AMR diverges strongly from the fixed grid approach - with more power on small scales in the AMR simulations - even at fixed physical resolution and also produces offsets in the star formation at specific epochs. This is because before certain times the upper grid levels are held back to maintain approximately fixed physical resolution, and to mimic the natural evolution of dark matter only simulations. Although the impact of hold back falls with increasing spatial and initial-condition resolutions, the offsets in the star formation remain down to a spatial resolution of 1 kpc. These offsets are of order of 10-20%, which is below the uncertainty in the implemented physics but are expected to affect the detailed properties of galaxies. We have implemented a new grid-hold-back approach to minimize the impact of hold back on the star formation rate.

MassiveNuS: cosmological massive neutrino simulations

The non-zero mass of neutrinos suppresses the growth of cosmic structure on small scales. Since the level of suppression depends on the sum of the masses of the three active neutrino species, the evolution of large-scale structure is a promising tool to constrain the total mass of neutrinos and possibly shed light on the mass hierarchy. In this work, we investigate these effects via a large suite of N-body simulations that include massive neutrinos using an analytic linear-response approximation: the Cosmological Massive Neutrino Simulations (MassiveNuS). The simulations include the effects of radiation on the background expansion, as well as the clustering of neutrinos in response to the nonlinear dark matter evolution. We allow three cosmological parameters to vary: the neutrino mass sum Mν in the range of 0–0.6 eV, the total matter density Ωm, and the primordial power spectrum amplitude As. The rms density fluctuation in spheres of 8 comoving Mpc/h (σ8) is a derived parameter as a result. Our data products include N-body snapshots, halo catalogues, merger trees, ray-traced galaxy lensing convergence maps for four source redshift planes between zs=1–2.5, and ray-traced cosmic microwave background lensing convergence maps. We describe the simulation procedures and code validation in this paper. The data are publicly available at http://columbialensing.org.

New constraints on dark matter production during kination

Our ignorance of the period between the end of inflation and the beginning of Big Bang Nucleosynthesis limits our understanding of the origins and evolution of dark matter. One possibility is that the Universe's energy density was dominated by a fast-rolling scalar field while the radiation bath was hot enough to thermally produce dark matter. We investigate the evolution of the dark matter density and derive analytic expressions for the dark matter relic abundance generated during such a period of kination. Kination scenarios in which dark matter does not reach thermal equilibrium require $\langle \sigma v \rangle < 2.7\times 10^{-38} \,\mathrm{cm^3\,s^{-1}}$ to generate the observed dark matter density while allowing the Universe to become radiation dominated by a temperature of $3 \, \mathrm{MeV}$. Kination scenarios in which dark matter does reach thermal equilibrium require $\langle \sigma v \rangle > 3\times 10^{-26} \,\mathrm{cm^3\,s^{-1}}$ in order to generate the observed dark matter abundance. We use observations of dwarf spheroidal galaxies by the Fermi Gamma-Ray Telescope and observations of the Galactic Center by the High Energy Stereoscopic System to constrain these kination scenarios. Combining the unitarity constraint on $\langle \sigma v \rangle$ with these observational constraints sets a lower limit on the temperature at which the Universe can become radiation dominated following a period of kination if ${\langle \sigma v \rangle > 3\times 10^{-31} \,\mathrm{cm^3\,s^{-1}}}$. This lower limit is between ${0.05 \, \mathrm{GeV}}$ and ${1 \, \mathrm{GeV}}$, depending on the dark matter annihilation channel.

THE BARYON CYCLE AT HIGH REDSHIFTS: EFFECTS OF GALACTIC WINDS ON GALAXY EVOLUTION IN OVERDENSE AND AVERAGE REGIONS

ABSTRACT We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at z ∼ 6–12 in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent from the “normal” region, and to shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on the environment, the more physically motivated VW model does exhibit this effect. In addition, VW can reproduce the observed specific SFR (sSFR) and the sSFR–stellar mass relation, which CW and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the intergalactic medium (IGM). The difference lies in the volume-filling factor of hot, high-metallicity gas, which is near unity for CW, while such gas remains confined in massive filaments for VW, and locked up in galaxies for NW. Such gas is nearly absent from the normal region. Although all wind models suffer from deficiencies, the VW model seems to be promising in correlating the outflow properties with those of host galaxies. Further constraints on the state of the IGM at high z are needed to separate different wind models.

Spurious haloes and discreteness-driven relaxation in cosmological simulations

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

The (Anti)Social Life of a Digital Journal in the Era of Global Info-Capitalism

By reflecting on running the open-access online “race” journal darkmatter, the author examines the crisis of academic journal publishing. This crisis is manifested by the transformations in the economics of journal production, especially with the rise of open-access publishing and the neoliberalization of the academy; the developments in digital communications; and the global consolidation of Western corporate publishing power in “info-capitalism.” The evolution of darkmatter is poised in terms of negotiating this crisis as opportunity and risk. The digital journal is able to produce a range of multimedia outputs addressing the interdisciplinary fields of race and postcolonial study. The challenge of emerging Web projects and the proliferation of online information potentially repositions darkmatter as the privileged site for the publication of scholarly research. The essay speculates on this critical journal in a time of “information overload” as a site of “autonomous resistance” in the increasing commodification of knowledge.

THE CONVERGENCE OF PARTICLE-IN-CELL SCHEMES FOR COSMOLOGICAL DARK MATTER SIMULATIONS

ABSTRACT Particle methods are a ubiquitous tool for solving the Vlasov–Poisson equation in comoving coordinates, which is used to model the gravitational evolution of dark matter (DM) in an expanding universe. However, these methods are known to produce poor results on idealized test problems, particularly at late times, after the particle trajectories have crossed. To investigate this, we have performed a series of one- and two-dimensional “Zel’dovich pancake” calculations using the popular particle-in-cell (PIC) method. We find that PIC can indeed converge on these problems provided that the following modifications are made. The first modification is to regularize the singular initial distribution function by introducing a small but finite artificial velocity dispersion. This process is analogous to artificial viscosity in compressible gas dynamics, and, as with artificial viscosity, the amount of regularization can be tailored so that its effect outside of a well-defined region—in this case, the high-density caustics—is small. The second modification is the introduction of a particle remapping procedure that periodically reexpresses the DM distribution function using a new set of particles. We describe a remapping algorithm that is third-order accurate and adaptive in phase space. This procedure prevents the accumulation of numerical errors in integrating the particle trajectories from growing large enough to significantly degrade the solution. Once both of these changes are made, PIC converges at second order on the Zel’dovich pancake problem, even at late times, after many caustics have formed. Furthermore, the resulting scheme does not suffer from the unphysical, small-scale “clumping” phenomenon known to occur on the pancake problem when the perturbation wavevector is not aligned with one of the Cartesian coordinate axes.

The role of the dark matter haloes on the cosmic star formation rate

The cosmic star formation rate (CSFR) represents the fraction of gas that is converted into stars within a certain comoving volume and at a given time t. However the evolution of the dark matter haloes and its relationship with the CSFR is not yet clear. In this context, we have investigated the role of the dark halo mass function - DHMF - in the process of gas conversion into stars. We observed a strong dependence between the fraction of baryons in structures, fb, and the specific mass function used for describing the dark matter haloes. In some cases, we have obtained fb greater than one at redshift z=0. This result indicates that the evolution of dark matter, described by the specific DHMF, could not trace the baryonic matter without a bias parameter. We also observed that the characteristic time-scale for star formation, τ, is strongly dependent on the considered DHMF, when the model is confronted against the observational data. Also, as part of this work it was released, under GNU general public license, a Python package called ‘pycosmicstar’ to study the CSFR and its relationship with the DHMF.