Halo Distribution Research Articles

A high-resolution cosmological simulation of a strong gravitational lens

ABSTRACT We present a cosmological hydrodynamical simulation of a 1013 M⊙ galaxy group and its environment (out to 10 times the virial radius) carried out using the eagle model of galaxy formation. Exploiting a novel technique to increase the resolution of the dark matter calculation independently of that of the gas, the simulation resolves dark matter haloes and subhaloes of mass 5 × 106 M⊙. It is therefore useful for studying the abundance and properties of the haloes and subhaloes targeted in strong lensing tests of the cold dark matter model. We estimate the halo and subhalo mass functions and discuss how they are affected both by the inclusion of baryons in the simulation and by the environment. We find that the halo and subhalo mass functions have lower amplitude in the hydrodynamical simulation than in its dark-matter-only counterpart. This reflects the reduced growth of haloes in the hydrodynamical simulation due to the early loss of gas by reionization and galactic winds and, additionally, in the case of subhaloes, disruption by enhanced tidal effects within the host halo due to the presence of a massive central galaxy. The distribution of haloes is highly anisotropic reflecting the filamentary character of mass accretion on to the cluster. As a result, there is significant variation in the number of structures with viewing direction. The median number of structures near the centre of the halo, when viewed in projection, is reduced by a factor of 2 when baryons are included.

Dark matter halo sparsity of modified gravity scenarios

Modified Gravity (MG) scenarios have been advocated to account for the dark energy phenomenon in the universe. These models predict departures from General Relativity on large cosmic scales that can be tested through a variety of probes such as observations of galaxy clusters among others. Here, we investigate the imprint of MG models on the internal mass distribution of cluster-like halos as probed by the dark matter halo sparsity. To this purpose we perform a comparative analysis of the properties of the halo sparsity using N-body simulation halo catalogs of a standard flat $\Lambda$CDM model and MG scenarios from the DUSTGRAIN-pathfinder simulation suite. We find that the onset of the screening mechanism leaves a distinct signature in the redshift evolution of the ensemble average halos sparsity. Measurements of the sparsity of galaxy clusters from currently available mass estimates are unable to test MG models due to the large uncertainties on the cluster masses. We show that this should be possible in the future provided large cluster samples with cluster masses determined to better than $30\%$ accuracy level.

Faint dark matter annihilation signals and the Milky Way’s supermassive black hole

A wide range of mechanisms predict present-day s-wave dark matter (DM) annihilation cross-sections that are orders of magnitude below current experimental sensitivity. We explore the capability of DM density spikes around the Milky Way's supermassive black hole to probe such faint signals of DM annihilations, considering a range of possible spike and halo distributions. As an exemplar of a theory with a suppressed s-wave annihilation cross-section, we consider a hidden sector axion portal model of DM. In this model, the leading contribution to the annihilation cross-section in the early universe is p-wave, while s-wave annihilations occur at higher order in the coupling constant. We provide a unified treatment of DM freezeout in this model including both s- and p-wave annihilations and analytically determine the photon spectrum for the dominant DM annihilation process in the universe today. We find that Fermi and H.E.S.S. observations of the Galactic Center offer excellent sensitivity to this model over a wide range of parameter space, with prospects depending sensitively on the properties of the DM spike as well as the central halo.

The past history of galaxy clusters told by their present neighbours

ABSTRACT Galaxy clusters can play a key role in modern cosmology, provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass–environment relations. This paper presents a mass-free correlation between the present neighbour distribution of cluster-size haloes and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random haloes with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size haloes to 200 Virgo-like haloes extracted from simulations that resemble the local Universe. The number of neighbours at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like haloes, clusters with numerous neighbours within a distance of about two times their virial radius experience a transition at z ≈ 1 between an active period of mass accretion, relative to the mean, and a quiet history. In contrary, clusters with few neighbours share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about four times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbours provides invaluable insights into the past history of these objects.

Statistics of solar wind electron breakpoint energies using machine learning techniques

Solar wind electron velocity distributions at 1 au consist of a thermal “core” population and two suprathermal populations: “halo” and “strahl”. The core and halo are quasi-isotropic, whereas the strahl typically travels radially outwards along the parallel or anti-parallel direction with respect to the interplanetary magnetic field. Using Cluster-PEACE data, we analyse energy and pitch angle distributions and use machine learning techniques to provide robust classifications of these solar wind populations. Initially, we used unsupervised algorithms to classify halo and strahl differential energy flux distributions to allow us to calculate relative number densities, which are of the same order as previous results. Subsequently, we applied unsupervised algorithms to phase space density distributions over ten years to study the variation of halo and strahl breakpoint energies with solar wind parameters. In our statistical study, we find both halo and strahl suprathermal breakpoint energies display a significant increase with core temperature, with the halo exhibiting a more positive correlation than the strahl. We conclude low energy strahl electrons are scattering into the core at perpendicular pitch angles. This increases the number of Coulomb collisions and extends the perpendicular core population to higher energies, resulting in a larger difference between halo and strahl breakpoint energies at higher core temperatures. Statistically, the locations of both suprathermal breakpoint energies decrease with increasing solar wind speed. In the case of halo breakpoint energy, we observe two distinct profiles above and below 500 km s−1. We relate this to the difference in origin of fast and slow solar wind.

The dark matter component of the Gaia radially anisotropic substructure

We study the properties of the dark matter component of the radially anisotropic stellar population recently identified in the Gaia data, using magneto-hydrodynamical simulations of Milky Way-like halos from the Auriga project. We identify 10 simulated galaxies that approximately match the rotation curve and stellar mass of the Milky Way. Four of these have an anisotropic stellar population reminiscent of the Gaia structure. We find an anti-correlation between the dark matter mass fraction of this population in the Solar neighbourhood and its orbital anisotropy. We estimate the local dark matter density and velocity distribution for halos with and without the anisotropic stellar population, and use them to simulate the signals expected in future xenon and germanium direct detection experiments. We find that a generalized Maxwellian distribution fits the dark matter halo integrals of the Milky Way-like halos containing the radially anisotropic stellar population. For dark matter particle masses below approximately 10 GeV, direct detection exclusion limits for the simulated halos with the anisotropic stellar population show a mild shift towards smaller masses compared to the commonly adopted Standard Halo Model.

Oral vesiculobullous lesions associated with SARS-CoV-2 infection.

Oral vesiculobullous lesions associated with SARS-CoV-2 infection.

T-ReX: a graph-based filament detection method

Numerical simulations and observations show that galaxies are not uniformly distributed in the universe but, rather, they are spread across a filamentary structure. In this large-scale pattern, highly dense regions are linked together by bridges and walls, all of them surrounded by vast, nearly-empty areas. While nodes of the network are widely studied in the literature, simulations indicate that half of the mass budget comes from a more diffuse part of the network, which is made up of filaments. In the context of recent and upcoming large galaxy surveys, it becomes essential that we identify and classify features of the Cosmic Web in an automatic way in order to study their physical properties and the impact of the cosmic environment on galaxies and their evolution. In this work, we propose a new approach for the automatic retrieval of the underlying filamentary structure from a 2D or 3D galaxy distribution using graph theory and the assumption that paths that link galaxies together with the minimum total length highlight the underlying distribution. To obtain a smoothed version of this topological prior, we embedded it in a Gaussian mixtures framework. In addition to a geometrical description of the pattern, a bootstrap-like estimate of these regularised minimum spanning trees allowed us to obtain a map characterising the frequency at which an area of the domain is crossed. Using the distribution of halos derived from numerical simulations, we show that the proposed method is able to recover the filamentary pattern in a 2D or 3D distribution of points with noise and outliers robustness with a few comprehensible parameters.

Neural physical engines for inferring the halo mass distribution function

ABSTRACT An ambitious goal in cosmology is to forward model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ∼Mpc scales and the distribution of dark matter haloes (used as a proxy for galaxies) in the same region. Here, we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field while inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded borg (Bayesian Origin Reconstruction from Galaxies) inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. In preliminary results, we notice there could potentially be orders of magnitude improvement in modelling compared to classical biasing techniques.

Corrections to halo model due to primordial magnetic field in a universe with non-zero cosmological constant

We study the corrections to the halo model of structure formation which is traditionally described by the Press Schecter formalism due to the presence of a primordial magnetic field in a universe which has a non-zero cosmological constant $\varLambda $ and positive curvature $k=+1$. We have used the result of Lokas and Hoffman (Astrophys. J. Lett. 542:139, 2000) and our earlier work (Nigam and Varalakshmi in Astrophys. Space Sci. 362:10, 2017) to incorporate the influence of both B field and $\varLambda $. We tune the value of the primordial magnetic field strength from 0.1 nG to 10 nG while varying the parameter $\lambda $ which is defined in terms of the cosmological constant $\varLambda $ and the Hubble constant $H_{o}$ from 0.5 to 0.75. We study the change in the Press Schecter number distribution of the dark matter halos as a function of these parameters. We find that while in general the presence of $\lambda $ disrupts the halo formation in presence of background primordial magnetic field, but as the value of $\lambda $ is increased above 0.65, there is a sudden fall in the halo count. Further we find that the model is far more sensitive to the variation in $\lambda $ than to the primordial magnetic field. This essentially means that the halo count has a strong correlation with $\lambda $ for fixed field strength than with the field strength for a fixed $\lambda $.

The bias of dark matter tracers: assessing the accuracy of mapping techniques

ABSTRACT We present a comparison between approximated methods for the construction of mock catalogues based on the halo-bias mapping technique. To this end, we use as reference a high-resolution N-body simulation of 38403 dark matter particles on a 400 h−1 Mpc cube box from the Multidark suite. In particular, we explore parametric versus non-parametric bias mapping approaches and compare them at reproducing the halo distribution in terms of the two- and three-point statistics down to $\sim 10^8\, {\rm M}_{\odot }\, h^{-1}$ halo masses. Our findings demonstrate that the parametric approach remains inaccurate even including complex deterministic and stochastic components. On the contrary, the non-parametric one is indistinguishable from the reference N-body calculation in the power spectrum beyond $k=1\, h\, {\rm Mpc}^{-1}$, and in the bispectrum for typical configurations relevant to baryon acoustic oscillation analysis. We conclude that approaches which extract the full bias information from N-body simulations in a non-parametric fashion are ready for the analysis of the new generation of large-scale structure surveys.

Resolving the Metallicity Distribution of the Stellar Halo with the H3 Survey

Abstract The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, , and there is no discernible metallicity gradient over the range 6 < R gal < 100 kpc. However, the halo metallicity distribution is highly structured, including distinct metal-rich and metal-poor components at R gal < 10 kpc and R gal > 30 kpc, respectively. The Sagittarius stream dominates the metallicity distribution at 20–40 kpc for stars on prograde orbits. The Gaia–Enceladus merger remnant dominates the metallicity distribution for radial orbits to ≈30 kpc. Metal-poor stars with [Fe/H] < −2 are a small population of the halo at all distances and orbital categories. We associate the “in situ” stellar halo with stars displaying thick disk chemistry on halo-like orbits; such stars are confined to . The majority of the stellar halo is resolved into discrete features in chemical–orbital space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. The relatively high metallicity of the halo derived in this work is a consequence of the unbiased selection function of halo stars and, in combination with the recent upward revision of the total stellar halo mass, implies a Galactic halo metallicity that is typical for its mass.

The bias of cosmic voids in the presence of massive neutrinos

Cosmic voids offer an extraordinary opportunity to study the effects of massive neutrinos on cosmological scales. Because they are freely streaming, neutrinos can penetrate the interior of voids more easily than cold dark matter or baryons, which makes their relative contribution to the mass budget in voids much higher than elsewhere in the Universe. In simulations it has recently been shown how various characteristics of voids in the matter distribution are affected by neutrinos, such as their abundance, density profiles, dynamics, and clustering properties. However, the tracers used to identify voids in observations (e.g., galaxies or halos) are affected by neutrinos as well, and isolating the unique neutrino signatures inherent to voids becomes more difficult. In this paper we make use of the DEMNUni suite of simulations to investigate the clustering bias of voids in Fourier space as a function of their core density and compensation. We find a clear dependence on the sum of neutrino masses that remains significant even for void statistics extracted from halos. In particular, we observe that the amplitude of the linear void bias increases with neutrino mass for voids defined in dark matter, whereas this trend gets reversed and slightly attenuated when measuring the relative void-halo bias using voids identified in the halo distribution. Finally, we argue how the original behaviour can be restored when considering observations of the total matter distribution (e.g. via weak lensing), and comment on scale-dependent effects in the void bias that may provide additional information on neutrinos in the future.

Impact of substructure on local dark matter searches

Dark matter substructure can contribute significantly to local dark matter searches and may provide a large uncertainty in the interpretation of those experiments. For direct detection experiments, sub-halos give rise to an additional dark matter component on top of the smooth dark matter distribution of the host halo. In the case of dark matter capture in the Sun, sub-halo encounters temporarily increase the number of captured particles. Even if the encounter happened in the past, the number of dark matter particles captured by the Sun can still be enhanced today compared to expectations from the host halo as those enhancements decay over time. Using results from an analytical model of the sub-halo population of a Milky Way-like galaxy, valid for sub-halo masses between 10−5 M⊙ and 1011 M⊙, we assess the impact of sub-halos on direct dark matter searches in a probabilistic way. We find that the impact on direct detection can be sizable, with a probability of ∼ 10−3 to find an \U0001d4aa(1) enhancement of the recoil rate. In the case of the capture rate in the Sun, we find that \U0001d4aa(1) enhancements are very unlikely, with probability ≲ 10−5, and are even impossible for some dark matter masses.

Dark Matter Haloes and Subhaloes

The development of methods and algorithms to solve the N-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark matter structures over most of the universe’s history. In the best-studied case—the cold dark matter or CDM model—the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure, and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from N-body simulations and on the challenges ahead.

Constraining cluster masses from the stacked phase space distribution at large radii

Abstract Velocity dispersions have been employed as a method to measure masses of clusters. To complement this conventional method, we explore the possibility of constraining cluster masses from the stacked phase space distribution of galaxies at larger radii, where infall velocities are expected to have a sensitivity to cluster masses. First, we construct a two-component model of the three-dimensional phase space distribution of haloes surrounding clusters up to 50 $\, h^{-1}$ Mpc from cluster centres based on N-body simulations. We confirm that the three-dimensional phase space distribution shows a clear cluster mass dependence up to the largest scale examined. We then calculate the probability distribution function of pairwise line-of-sight velocities between clusters and haloes by projecting the three-dimensional phase space distribution along the line of sight with the effect of the Hubble flow. We find that this projected phase space distribution, which can directly be compared with observations, shows a complex mass dependence due to the interplay between infall velocities and the Hubble flow. Using this model, we estimate the accuracy of dynamical mass measurements from the projected phase space distribution at the transverse distance from cluster centres larger than $2\, h^{-1}$ Mpc. We estimate that, by using 1.5 × 105 spectroscopic galaxies, we can constrain the mean cluster masses with an accuracy of 14.5 per cent if we fully take account of the systematic error coming from the inaccuracy of our model. This can be improved down to 5.7 per cent by improving the accuracy of the model.

Painting halos from cosmic density fields of dark matter with physically motivated neural networks

We present a novel halo painting network that learns to map approximate 3D dark matter fields to realistic halo distributions. This map is provided via a physically motivated network with which we can learn the non-trivial local relation between dark matter density field and halo distributions without relying on a physical model. Unlike other generative or regressive models, a well motivated prior and simple physical principles allow us to train the mapping network quickly and with relatively little data. In learning to paint halo distributions from computationally cheap, analytical and non-linear density fields, we bypass the need for full particle mesh simulations and halo finding algorithms. Furthermore, by design, our halo painting network needs only local patches of dark matter density to predict the halos, and as such, it can predict the 3D halo distribution for any arbitrary simulation box size. Our neural network can be trained using small simulations and used to predict large halo distributions, as long as the resolutions are equivalent. We evaluate our model's ability to generate 3D halo count distributions which reproduce, to a high degree, summary statistics such as the power spectrum and bispectrum, of the input or reference realizations.

The Hunt for Primordial Interactions in the Large-Scale Structures of the Universe

The understanding of the primordial mechanism that seeded the cosmic structures we observe today in the sky is one of the major goals in cosmology. The leading paradigm for such a mechanism is provided by the inflationary scenario, a period of violent accelerated expansion in the very early stages of evolution of the universe. While our current knowledge of the physics of inflation is limited to phenomenological models which fit observations, an exquisite understanding of the particle content and interactions taking place during inflation would provide breakthroughs in our understanding of fundamental physics at high energies. In this review, we summarize recent theoretical progress in the modeling of the imprint of primordial interactions in the large-scale structures of the universe. We focus specifically on the effects of such interactions on the statistical distribution of dark-matter halos, providing a consistent treatment of the steps required to connect the correlations generated among fields during inflation all the way to the late-time correlations of halos.

Halo concentration, galaxy red fraction, and gas properties of optically defined merging clusters

Abstract We present multi-wavelength studies of optically defined merging clusters, based on the Hyper Suprime-Cam Subaru Strategic Program. Luminous red galaxies, tracing cluster mass distributions, enable us to identify cluster subhalos at various merging stages, and thus make a homogeneous sample of cluster mergers that is unbiased with respect to the merger boost of the intracluster medium (ICM). We define, using a peak-finding method, merging clusters with multiple peaks and single clusters with single peaks from the CAMIRA cluster catalog. Stacked weak-lensing analysis indicates that our sample of merging clusters is categorized into major mergers. The average halo concentration for the merging clusters is ∼70% smaller than that of the single-peak clusters, which agrees well with predictions of numerical simulations. The spatial distribution of subhalos is less centrally concentrated than the mass distribution of the main halo. The fractions of red galaxies in the merging clusters are not higher than those of the single-peak clusters. We find a signature of the merger boost of the ICM from the stacked Planck Sunyaev–Zel’dovich effect and ROSAT X-ray luminosity, but not in optical richness. The stacked X-ray surface brightness distribution, aligned with the main subhalo pairs of low-redshift and massive clusters, shows that the central gas core is elongated along the merger axis, and overall gas distribution is misaligned by ∼60°. The homogeneous, unbiased sample of cluster mergers and multi-wavelength follow-up studies provide a unique opportunity to make a complete picture of merger physics over the whole process.

A comparison of the R_{mathrm{h}}=ct and varLambda CDM cosmologies based on the observed halo mass function

The growth of structure may be traced via the redshift-dependent halo mass function. This quantity probes the re-ionization history and quasar abundance in the Universe, constituting an important probe of the cosmological predictions. Halos are not directly observable, however, so their mass and evolution must be inferred indirectly. The most common approach is to presume a relationship with galaxies and halos. Studies based on the assumption of a constant halo to stellar mass ratio M_h/M_* (extrapolated from z lesssim 4) reveal significant tension with varLambda CDM – a failure known as “The Impossibly Early Galaxy Problem”. But whether this ratio evolves or remains constant through redshift 4 lesssim z lesssim 10 is still being debated. To eliminate the tension with varLambda CDM, it would have to change by about 0.8 dex over this range, an issue that may be settled by upcoming observations with the James Webb Space Telescope. In this paper, we explore the possibility that this major inconsistency may instead be an indication that the cosmological model is not completely correct. We study this problem in the context of another Friedmann–Lemaître–Robertson–Walker (FLRW) model known as the R_{mathrm{h}}=ct universe, and use our previous measurement of sigma _8 from the cosmological growth rate, together with new solutions to the Einstein–Boltzmann equations, to interpret these recent halo measurements. We demonstrate that the predicted mass and redshift dependence of the halo distribution in R_{mathrm{h}}=ct is consistent with the data, even assuming a constant M_h/M_* throughout the observed redshift range (4lesssim zlesssim 10), contrasting sharply with the tension in varLambda CDM. We conclude that – if M_h/M_* turns out to be constant – the massive galaxies and their halos must have formed earlier than is possible in varLambda CDM.