Search

Secure confirmation that a gravitational wave (GW) has been gravitationally lensed would bring together these two pillars of General Relativity for the first time. This breakthrough is challenging for many reasons, including: GW sky localization uncertainties dwarf the angular scale of gravitational lensing, the mass and structure of gravitational lenses is diverse, the mass function of stellar remnant compact objects is not yet well constrained, and GW detectors do not operate continuously. We introduce a new approach that is agnostic to the mass and structure of the lenses, compare the efficiency of different methods for lensed GW discovery, and explore detection of lensed kilonova counterparts as a direct method for localising candidates. Our main conclusions are: (1) lensed neutron star mergers (NS-NS) are magnified into the "mass gap" between NS and black holes, therefore selecting candidates from public GW alerts with high mass gap probability is efficient, (2) the rate of detectable lensed NS-NS will approach one per year in the mid-2020s, (3) the arrival time difference between lensed NS-NS images is $1\,\rm sec\lesssim\Delta t\lesssim1\,year$, and thus well-matched to the operations of GW detectors and optical telescopes, (4) lensed kilonova counterparts are faint at peak (e.g.\ $r_{\rm AB}\simeq24-26$ in the mid-2020s), fade quickly ($d, Comment: Accepted by MNRAS. 20 pages, 10 figures

Journal of cosmology and astroparticle physics 01(01), 009 (2023). doi:10.1088/1475-7516/2023/01/009, We show that very compact axion mini-clusters can form in models where axion-like-particle (ALP) dark matter is produced via the kinetic misalignment mechanism, which is well-motivated in pre-inflationary $U(1)$ symmetry breaking scenarios. This is due to ALP fragmentation. We predict denser halos than what has been obtained so far in the literature from standard misalignment in post-inflationary $U(1)$ breaking scenarios or from large misalignment. The main reason is that adiabatic fluctuations are significant at early times; therefore, even if amplification from parametric resonance effects is moderate, the final size of ALP fluctuations is larger in kinetic misalignment. We compare halo mass functions and halo spectra obtained in kinetic misalignment, large misalignment, and standard misalignment, respectively. Our analysis does not depend on the specific model realization of the kinetic misalignment mechanism. We present our results generally as a function of the ALP mass and the ALP decay constant only. We show that a sizable region of this ALP parameter space can be tested by future experiments that probe small-scale structures., Published by IOP, London

It is generally assumed within the standard cosmological model that initial density perturbations are Gaussian at all scales. However, primordial quantum diffusion unavoidably generates non-Gaussian, exponential tails in the distribution of inflationary perturbations. These exponential tails have direct consequences for the formation of collapsed structures in the universe, as has been studied in the context of primordial black holes. We show that these tails also affect the very-large-scale structures, making heavy clusters like "El Gordo", or large voids like the one associated with the cosmic microwave background cold spot, more probable. We compute the halo mass function and cluster abundance as a function of redshift in the presence of exponential tails. We find that quantum diffusion generically enlarges the number of heavy clusters and depletes subhalos, an effect that cannot be captured by the famed $f_{\mathrm{NL}}$ corrections. These late-universe signatures could thus be fingerprints of quantum dynamics during inflation that should be incorporated in $N$-body simulations and checked against astrophysical data., 5 pages without appendices (total 8 pages), 3 figures, matches published version in Phys. Rev. Letters

The unusually large "dwarf" galaxy Crater II, with its small velocity dispersion, $\simeq 3$ km/s, defies expectations that low mass galaxies should be small and dense. We combine the latest stellar and velocity dispersion profiles finding Crater II has a prominent dark core of radius $\simeq 0.71^{+0.09}_{-0.08}$ kpc, surrounded by a low density halo, with a transition visible between the core and the halo. We show that this profile matches the distinctive core-halo profile predicted by "Wave Dark Matter" as a Bose-Einstein condensate, $\psi$DM, where the ground state soliton core is surrounded by a tenuous halo of interfering waves, with a marked density transition predicted between the core and halo. Similar core-halo structure is seen in most dwarf spheroidal galaxies (dSph), but with smaller cores, $\simeq 0.25$ kpc and higher velocity dispersions, $\simeq 9$km/s, and we argue here that Crater II may have been a typical dSph that has lost most of its halo mass to tidal stripping, so its velocity dispersion is lower by a factor of 3 and the soliton is wider by a factor of 3, following the inverse scaling required by the Uncertainty Principle. This tidal solution for Crater II in the context of $\psi$DM, is supported by its small pericenter of $\simeq 20$ kpc established by Gaia, implying significant tidal stripping of Crater II by the Milky Way is expected., Comment: Accepted version by MNRAS. 10 pages, 6 figures and 1 table

In this work, we investigate the structural properties, distribution and abundance of ΛCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify – for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass – dependencies of subhalo properties on mass, maximum circular velocity, Vmax, host halo mass, and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. Subhalo structural properties are codified via a concentration parameter, cV, that does not depend on any pre-defined density profile and relies only on Vmax. We derive the cV−Vmax relation and find an important dependence on distance of the subhalo to the host halo centre. Interestingly, we also find subhaloes of the same mass to be significantly more concentrated when they reside inside more massive hosts. Finally, we investigate the redshift evolution of cV, and provide accurate fits. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Thus, we expect our work to be particularly useful for any future research involving dark matter halo substructure. © 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., The work of AM, MASC, and AAS was supported by the Spanish Agencia Estatal de Investigación through the grants PGC2018-095161-B-I00 and IFT Centro de Excelencia Severo Ochoa SEV-2016-0597, the Atracción de Talento contracts no. 2016-T1/TIC-1542 and 2020-5A/TIC-19725 granted by the Comunidad de Madrid in Spain, and the MultiDark Consolider Network FPA2017-90566-REDC. The work of AM was also supported by ‘Comunidad de Madrid S2018/NMT-4291 TEC2SPACE-CM’. The work of AAS was also supported by the Spanish Ministry of Science and Innovation through the grant FPI-UAM 2018. FP thanks the support of the Spanish Ministry of Science and Innovation funding grant PGC2018- 101931-B-I00. TI was supported by MEXT as ‘Program for Promoting Researches on the Supercomputer Fugaku’ (JPMXP1020200109), MEXT/JSPS KAKENHI Grant Number JP19KK0344, JP20H05245, JP21H01122 and IAAR Research Support Program, Chiba University, Japan. SAC acknowledges funding from Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET, PIP-0387), Agencia Nacional de Promoción de la Investigación, el Desarrollo Tecnológico y la Innovación (Agencia I+D + i, PICT-2018-3743), and Universidad Nacional de La Plata (G11-150), Argentina. EJ acknowledges financial support from CNRS. TO is supported by JSPS KAKENHI Grant Number 20K22360 and 21H05449. This work used the skun6@IAA facility (http://www.skiesanduniverses.org) managed by the Instituto de Astrofísica de Andalucía (CSIC). The equipment was funded by the Spanish Ministry of Science EU-FEDER infrastructure grants EQC2018-004366-P and EQC2019-006089-P. The Uchuu, Shin-Uchuu, and Phi-4096 simulations were carried out on Aterui II supercomputer at Center for Computational Astrophysics, CfCA, of National Astronomical Observatory of Japan. The numerical analysis were partially carried out on XC40 at the Yukawa Institute Computer Facility in Kyoto University., With funding from the Spanish government through the "Severo Ochoa Centre of Excellence" accreditation (CEX2021-001131-S).

We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass Mν = 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of Mν. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

The spatial distribution of Milky Way (MW) subhaloes provides an important set of observables for testing cosmological models. These include the radial distribution of luminous satellites, planar configurations, and the abundance of dark subhaloes whose existence or absence is key to distinguishing amongst dark matter models. We use the COCO $N$-body simulations of cold dark matter (CDM) and 3.3keV thermal relic warm dark matter (WDM) to predict the satellite spatial distribution in the limit that the impact of baryonic physics is minimal. We demonstrate that the radial distributions of CDM and 3.3keV-WDM luminous satellites are identical if the minimum pre-infall halo mass to form a galaxy is $>10^{8.5}\mathrm{M}_{\odot}$. The distribution of dark subhaloes is significantly more concentrated in WDM due to the absence of low mass, recently accreted substructures that typically inhabit the outer parts of a MW halo in CDM. We show that subhaloes of mass $[10^{7},10^{8}]\mathrm{M}_{\odot}$ and within 30kpc of the centre are the stripped remnants of larger haloes in both models. Therefore their abundance in WDM is $3\times$ higher than one would anticipate from the overall WDM subhalo population. We estimate that differences between CDM and WDM concentration-mass relations can be probed for subhalo-stream impact parameters $, 15 pages, 4 figures. Accepted version of paper. Contact: lovell@hi.is

We quantify the information content of the non-linear matter power spectrum, the halo mass function, and the void size function, using the Quijote $N$-body simulations. We find that these three statistics exhibit very different degeneracies amongst the cosmological parameters, and thus the combination of all three probes enables the breaking of degeneracies, in turn yielding remarkably tight constraints. We perform a Fisher analysis using the full covariance matrix, including all auto- and cross-correlations, finding that this increases the information content for neutrino mass compared to a correlation-free analysis. The multiplicative improvement of the constraints on the cosmological parameters obtained by combining all three probes compared to using the power spectrum alone are: 137, 5, 8, 20, 10, and 43, for $\Omega_m$, $\Omega_b$, $h$, $n_s$, $\sigma_8$, and $M_\nu$, respectively. The marginalized error on the sum of the neutrino masses is $\sigma(M_\nu)=0.018\,{\rm eV}$ for a cosmological volume of $1\,(h^{-1}{\rm Gpc})^3$, using $k_{\max}=0.5\,h{\rm Mpc}^{-1}$, and without CMB priors. We note that this error is an underestimate insomuch as we do not consider super-sample covariance, baryonic effects, and realistic survey noises and systematics. On the other hand, it is an overestimate insomuch as our cuts and binning are suboptimal due to restrictions imposed by the simulation resolution. Given upcoming galaxy surveys will observe volumes spanning $\sim 100\,(h^{-1}{\rm Gpc})^3$, this presents a promising new avenue to measure neutrino mass without being restricted by the need for accurate knowledge of the optical depth, which is required for CMB-based measurements. Furthermore, the improved constraints on other cosmological parameters, notably $\Omega_m$, may also be competitive with CMB-based measurements., Comment: 11+5 pages, 10 figures, 2 tables

We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integrals that otherwise involve small-scale, non-perturbative modes. We then present results on the inferred value of the linear power spectrum normalization $\sigma_{8}$ from rest-frame halo catalogs using both second- and third-order bias expansions, imposing uniform priors on all bias parameters. Due to the perfect bias-$\sigma_{8}$ degeneracy at linear order, constraints on $\sigma_{8}$ rely entirely on nonlinear information. The results show the expected convergence behavior when lowering the cutoff in wavenumber, $\Lambda$. When including modes up to $k \leq \Lambda = 0.1\,h\,{\rm Mpc}^{-1}$ in the second-order case, $\sigma_{8}$ is recovered to within $\lesssim 6\,\%$ for a range of halo masses and redshifts. The systematic bias shrinks to $4\,\%$ or less for the third-order bias expansion on the same range of scales. Together with additional evidence we provide, this shows that the residual mismatch in $\sigma_{8}$ can be attributed to higher-order bias contributions. We conclude that the EFT likelihood is able to infer unbiased cosmological constraints, within expected theoretical systematic errors, from physical biased tracers on quasilinear scales, Comment: 36 pages, 10 figures, 2 tables; main results in Figs. 2, 4, 6; v2: Sec. 3 clarified, Sec. 5.4 expanded; matches published JCAP version

We study the core mass-halo mass relation of bosonic dark matter halos, in the form of self-gravitating Bose-Einstein condensates, harboring a supermassive black hole. We use the ``velocity dispersion tracing'' relation according to which the velocity dispersion in the core ${v}_{c}^{2}\ensuremath{\sim}G{M}_{c}/{R}_{c}$ is of the same order as the velocity dispersion in the halo ${v}_{h}^{2}\ensuremath{\sim}G{M}_{h}/{r}_{h}$ (this relation can be justified from thermodynamical arguments) and the approximate analytical mass-radius relation of the quantum core in the presence of a central black hole obtained in our previous paper [P.H. Chavanis, Eur. Phys. J. Plus 134, 352 (2019)]. For a given minimum halo mass $({M}_{h}{)}_{\mathrm{min}}\ensuremath{\sim}{10}^{8}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ determined by the observations, the only free parameter of our model is the scattering length ${a}_{s}$ of the bosons (their mass $m$ is then determined in order to match the characteristics of the minimum halo). For noninteracting bosons and for bosons with a repulsive self-interaction, we find that the core mass ${M}_{c}$ increases with the halo mass ${M}_{h}$ and achieves a maximum value $({M}_{c}{)}_{\mathrm{max}}$ at some halo mass $({M}_{h}{)}_{*}$ before decreasing. The whole series of equilibria is stable. For bosons with an attractive self-interaction, we find that the core mass achieves a maximum value $({M}_{c}{)}_{\mathrm{max}}$ at some halo mass $({M}_{h}{)}_{*}$ before decreasing. The series of equilibria becomes unstable above a maximum halo mass $({M}_{h}{)}_{\mathrm{max}}\ensuremath{\ge}({M}_{h}{)}_{*}$. In the absence of black hole $({M}_{h}{)}_{\mathrm{max}}=({M}_{h}{)}_{*}$. At that point, the quantum core (similar to a dilute axion star) collapses. We perform a similar study for fermionic dark matter halos. We find that they behave similarly to bosonic dark matter halos with a repulsive self-interaction, the Pauli exclusion principle for fermions playing the role of the repulsive self-interaction for bosons.

Cosmic voids, the underdense regions in the universe, are particularly sensitive to diffuse density components such as cosmic neutrinos. This sensitivity is enhanced by the match between void sizes and the free-streaming scale of massive neutrinos. Using the massive neutrino simulations \texttt{MassiveNuS}, we investigate the effect of neutrino mass on dark matter halos as a function of environment. We find that the halo mass function depends strongly on neutrino mass and that this dependence is more pronounced in voids than in high-density environments. An observational program that measured the characteristic mass of the most massive halos in voids should be able to place novel constraints on the sum of the masses of neutrinos $\sum m_\nu$. The neutrino mass effect in the simulations is quite strong: In a 512$^3$ $h^{-3}$ Mpc$^3$ survey, the mean mass of the 1000 most massive halos in the void interiors is $(4.82 \pm 0.11) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.6$ eV and $(8.21 \pm 0.13) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.1$ eV. Subaru (SuMIRe), Euclid and WFIRST will have both spectroscopic and weak lensing surveys. Covering volumes at least 50 times larger than our simulations, they should be sensitive probes of neutrino mass through void substructure., Comment: 5 pages, 2 figures

We consider models of scalar dark matter with a generic interaction potential and non-canonical kinetic terms of the K-essence type that are subleading with respect to the canonical term. We analyze the low-energy regime and derive, in the nonrelativistic limit, the effective equations of motions. In the fluid approximation they reduce to the conservation of matter and to the Euler equation for the velocity field. We focus on the case where the scalar field mass $10^{-21} \ll m \lesssim 10^{-4} \, {\rm eV}$ is much larger than for fuzzy dark matter, so that the quantum pressure is negligible on cosmological and galactic scales, while the self-interaction potential and non-canonical kinetic terms generate a significant repulsive pressure. At the level of cosmological perturbations, this provides a dark-matter density-dependent speed of sound. At the nonlinear level, the hydrostatic equilibrium obtained by balancing the gravitational and scalar interactions imply that virialized structures have a solitonic core of finite size depending on the speed of sound of the dark matter fluid. For the most relevant potential in $\lambda_4 \phi^4/4$ or K-essence with a $(\partial \phi)^4$ interaction, the size of such stable cores cannot exceed 60 kpc. Structures with a density contrast larger than $10^6$ can be accommodated with a speed of sound $c_s\lesssim 10^{-6}$. We also consider the case of a cosine self-interaction, as an example of bounded nonpolynomial self-interaction. This gives similar results in low-mass and low-density halos whereas solitonic cores are shown to be absent in massive halos., Comment: 25 pages

International audience; Lensing of the cosmic microwave background (CMB) is now a well-developed probe of the clustering of the large-scale mass distribution over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission will allow production of a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that will follow from its power spectrum and the cross-correlation with other clustering data. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous measurements of the baryon acoustic oscillation feature in the clustering of galaxies, three times smaller than the minimum total mass allowed by neutrino oscillation measurements. Lensing has applications across many other science goals of CORE, including the search for B-mode polarization from primordial gravitational waves. Here, lens-induced B-modes will dominate over instrument noise, limiting constraints on the power spectrum amplitude of primordial gravitational waves. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60 %. This can be improved to 70 % by combining lensing and measurements of the cosmic infrared background from CORE, leading to an improvement of a factor of 2.5 in the error on the amplitude of primordial gravitational waves compared to no delensing (in the null hypothesis of no primordial B-modes). Lensing measurements from CORE will allow calibration of the halo masses of the tens of thousands of galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. The 19 frequency channels proposed for CORE will allow accurate removal of Galactic emission from CMB maps. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE.

We examine the efficacy of the halo mass function as a probe of $f(R)$, symmetron and DGP gravity. We develop an excursion-set method to generalise a range of popular HMF fitting functions from GR to screened MG, considering the HMF dependence on the critical density parameter $\delta_{c}$. In particular, we propose a variety of new methods to account for the environmental dependence of chameleon screening and compare their accuracy to existing ones via $N$-body data. Finally, using the nested sampling routine MultiNest, we examine two possible interpretations of the MG $N$-body results: whether they can be interpreted as $\Lambda$CDM HMFs and whether the MG HMFs display any universality. We find that the effects of MG can be interpreted as a change in best-fit parameters in the $\Lambda$CDM HMF---i.e. can be mistaken for GR---for all of the fitting functions. Alternatively, the relation can be inverted to judge the universality of the HMF---not in terms of its purported invariance to $z$, $H_{0}$, $\Omega_{m}$, $\Omega_{\Lambda}$ as in GR, but its independence on the underlying theory of gravity. Although we found no completely universal HMF, the parameter values did cluster according to the type of screening mechanism. The results suggest that a single, best-fit HMF might be used for each type of screening, independent of the parameters in the MG model. This demonstrates that the additional complexity of the gravitational collapse in screened MG theories cannot always be accounted for using the techniques developed in GR. However, the variations are small enough that it is unnecessary to develop new fitting functions and calibrate them on a case-by-case basis., Comment: 61 pages, 21 figures Replaced introduction. Added background detail on excursion set theory. Shifted technical aspects to appendices. Results unchanged but explained with clearer references to figures

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged], 44 pages, 12 figures

We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $\nu\Lambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($\sim5-10\%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $\Omega_\nu-\sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{\rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $\Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories., Comment: 25 pages, 7 figures

Some of the dark matter in the Universe is made up of massive neutrinos. Their impact on the formation of large scale structure can be used to determine their absolute mass scale from cosmology, but to this end accurate numerical simulations have to be developed. Due to their relativistic nature, neutrinos pose additional challenges when one tries to include them in N-body simulations that are traditionally based on Newtonian physics. Here we present the first numerical study of massive neutrinos that uses a fully relativistic approach. Our N-body code, gevolution, is based on a weak-field formulation of general relativity that naturally provides a self-consistent framework for relativistic particle species. This allows us to model neutrinos from first principles, without invoking any ad-hoc recipes. Our simulation suite comprises some of the largest neutrino simulations performed to date. We study the effect of massive neutrinos on the nonlinear power spectra and the halo mass function, focusing on the interesting mass range between 0.06 eV and 0.3 eV and including a case for an inverted mass hierarchy., Comment: 26 pages + appendix, 17 figures; v2: minor revision, published version; code available at https://github.com/gevolution-code/gevolution-1.1

The influence of the large scale structure on host halos may be studied by examining the angular infall pattern of subhalos. In particular, since warm and cold dark matter cosmologies predict different abundances and internal properties for halos at the low mass end of the mass function, it is interesting to examine if there are differences in how these low mass halos are accreted. The accretion events are defined as the moment a halo becomes a substructure, namely when it crosses its host's virial radius. We quantify the cosmic web at each point by the shear tensor and examine where, with respect to its eigenvectors, such accretion events occur in cold ($\Lambda$CDM) and warm (1keV sterile neutrino WDM) dark matter cosmological models. We find that the CDM and WDM subhalos are preferentially accreted along the principal axis of the shear tensor corresponding to the direction of weakest collapse. The beaming strength is modulated by the host and subhalo masses and by the redshift at which the accretion event occurs. Although strongest for the most massive hosts and subhalos at high redshift, the preferential infall is found to be always aligned with the axis of weakest collapse, thus we say that it has universal nature. We compare the strength of beaming in the WDM cosmology with the one found in the $\Lambda$CDM scenario. While the main findings remain the same, the accretion in the WDM model for the most massive host halos appears more beamed than in $\Lambda$CDM cosmology across all the redshifts., Comment: 12 pages, 9 figures

We use state-of-art measurements of the galaxy luminosity function (LF) at z=6, 7 and 8 to derive constraints on warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ULADM) models alternative to the cold dark matter (CDM) paradigm. To this purpose we have run a suite of high-resolution N-body simulations to accurately characterise the low mass-end of the halo mass function and derive DM model predictions of the high-z luminosity function. In order to convert halo masses into UV-magnitudes we introduce an empirical approach based on halo abundance matching which allows us to model the LF in terms of the amplitude and scatter of the ensemble average star formation rate halo mass relation of each DM model, $\langle {\rm SFR}({\rm M_{ h}},z)\rangle$. We find that independent of the DM scenario the average SFR at fixed halo mass increases from z=6 to 8, while the scatter remains constant. At halo mass ${\rm M_{h}}\gtrsim 10^{12}\,{\rm M}_\odot$ h$^{-1}$ the average SFR as function of halo mass follows a double power law trend that is common to all models, while differences occur at smaller masses. In particular, we find that models with a suppressed low-mass halo abundance exhibit higher SFR compared to the CDM results. Using deviance statistics we obtain a lower limit on the WDM thermal relic particle mass, $m_{\rm WDM}\gtrsim 1.5$ keV at $2\sigma$. In the case of LFDM models, the phase transition redshift parameter is bounded to $z_t\gtrsim 8\cdot 10^5$ at $2\sigma$. We find ULADM best-fit models with axion mass $m_a\gtrsim 1.6\cdot 10^{-22}$ eV to be well within $2\sigma$ of the deviance statistics. We remark that measurements at $z=6$ slightly favour a flattening of the LF at faint UV-magnitudes. This tends to prefer some of the non-CDM models in our simulation suite, although not at a statistically significant level to distinguish them from CDM., Comment: 22 pages, 15 figures; added separate analysis using LF data at z=6 from Livermore et al. (2016) strongly constraining Axion DM models. PRD accepted version

We study the imprint of non-standard dark energy (DE) and dark matter (DM) models on the 21cm intensity map power spectra from high-redshift neutral hydrogen (HI) gas. To this purpose we use halo catalogs from N-body simulations of dynamical DE models and DM scenarios which are statistically indistinguishable from the standard Cold Dark Matter model with Cosmological Constant (LCDM) using currently available cosmological observations. We limit our analysis to halo catalogs at redshift z = 1 and 2.3 which are common to all simulations. For each catalog we model the HI distribution by using a simple prescription to associate the HI gas mass to N-body halos. We find that the DE models leave a distinct signature on the HI spectra across a wide range of scales, which correlates with differences in the halo mass function and the onset of the non-linear regime of clustering. In the case of the non-standard DM model significant differences of the HI spectra with respect to the LCDM model only arise from the suppressed abundance of low mass halos. These cosmological model dependent features also appear in the 21cm spectra. In particular, we find that future SKA measurements can distinguish the imprints of DE and DM models at high statistical significance., to be submitted, comments welcome

We compute the mass function of galactic dark matter halos for different values of the Warm Dark Matter (WDM) particle mass m_X and compare it with the abundance of ultra-faint galaxies derived from the deepest UV luminosity function available so far at redshift z~2. The magnitude limit M_UV=-13 reached by such observations allows us to probe the WDM mass functions down to scales close to or smaller than the half-mass mode mass scale ~10^9 M_sun. This allowed for an efficient discrimination among predictions for different m_X which turn out to be independent of the star formation efficiency adopted to associate the observed UV luminosities of galaxies to the corresponding dark matter masses. Adopting a conservative approach to take into account the existing theoretical uncertainties in the galaxy halo mass function, we derive a robust limit m_X>1.8 keV for the mass of thermal relic WDM particles when comparing with the measured abundance of the faintest galaxies, while m_X>1.5 keV is obtained when we compare with the Schechter fit to the observed luminosity function. The corresponding lower limit for sterile neutrinos depends on the modeling of the production mechanism; for instance m_sterile > 4 keV holds for the Shi-Fuller mechanism. We discuss the impact of observational uncertainties on the above bound on m_X. As a baseline for comparison with forthcoming observations from the HST Frontier Field, we provide predictions for the abundance of faint galaxies with M_UV=-13 for different values of m_X and of the star formation efficiency, valid up to z~4., 14 pages, 3 figures. Accepted for publication in The Astrophysical Journal