Search

Observational cosmology is rapidly closing in on a measurement of the sum M_nu of neutrino masses, at least in the simplest cosmologies, while opening the door to probes of non-standard hot dark matter (HDM) models. By extending the method of effective distributions, we show that any collection of HDM species, with arbitrary masses, temperatures, and distribution functions, including massive neutrinos, may be represented as a single effective HDM species. Implementing this method in the FlowsForTheMasses non-linear perturbation theory for free-streaming particles, we study non-standard HDM models that contain thermal QCD axions or generic bosons in addition to standard neutrinos, as well as non-standard neutrino models wherein either the distribution function of the neutrinos or their temperature is changed. Along the way, we substantially improve the accuracy of this perturbation theory at low masses, bringing it into agreement with the high-resolution TianNu neutrino N-body simulation to about 2% at k = 0.1 h/Mpc and to within 21% up to k = 1 h/Mpc. We accurately reproduce the results of simulations including axions and neutrinos of multiple masses. Studying the differences between the normal, inverted, and degenerate neutrino mass orderings on their non-linear power, we quantify the error in the common approximation of degenerate masses. We release our code publicly at http://github.com/upadhye/FlowsForTheMassesII ., Comment: 39 pages, 18 figures

We generalise the SuperEasy linear response method, originally developed to describe massive neutrinos in cosmological $N$-body simulations, to any hot dark matter (HDM) species with arbitrary momentum distributions. The method uses analytical solutions of the HDM phase space perturbations in various limits and constructs from them a modification factor to the gravitational potential that tricks the cold particles into trajectories as if HDM particles were present in the simulation box. The modification factor is algebraic in the cosmological parameters and requires no fitting. Implementing the method in a Particle-Mesh simulation code and testing it on HDM cosmologies up to the equivalent effect of $\sum m_\nu = 0.315$ eV-mass neutrinos, we find that the generalised SuperEasy approach is able to predict the total matter and cold matter power spectra to $\lesssim 0.1\%$ relative to other linear response methods and to $\lesssim 0.25\%$ relative to particle HDM simulations. Applying the method to cosmologies with mixed neutrinos\Plus{}thermal QCD axions and neutrinos\Plus{}generic thermal bosons, we find that non-standard HDM cosmologies have no intrinsically different non-linear signature in the total matter power spectrum from standard neutrino cosmologies. However, because they predict different time dependencies even at the linear level and the differences are augmented by non-linear evolution, it remains a possibility that observations at multiple redshifts may help distinguish between them., Comment: 27 pages, 13 figures, code available https://github.com/cppccosmo/gadget-4-cppc

Cosmology is poised to measure the neutrino mass sum $M_\nu$ and has identified several smaller-scale observables sensitive to neutrinos, necessitating accurate predictions of neutrino clustering over a wide range of length scales. The FlowsForTheMasses non-linear perturbation theory for the massive neutrino power spectrum, $\Delta^2_\nu(k)$, agrees with its companion N-body simulation at the $10\%-15\%$ level for $k \leq 1~h/$Mpc. Building upon the Mira-Titan IV emulator for the cold matter, we use FlowsForTheMasses to construct an emulator for $\Delta^2_\nu(k)$ covering a large range of cosmological parameters and neutrino fractions $\Omega_{\nu,0} h^2 \leq 0.01$, which corresponds to $M_\nu \leq 0.93$~eV. Consistent with FlowsForTheMasses at the $3.5\%$ level, it returns a power spectrum in milliseconds. Ranking the neutrinos by initial momenta, we also emulate the power spectra of momentum deciles, providing information about their perturbed distribution function. Comparing a $M_\nu=0.15$~eV model to a wide range of N-body simulation methods, we find agreement to $3\%$ for $k \leq 3 k_\mathrm{FS} = 0.17~h/$Mpc and to $19\%$ for $k \leq 0.4~h/$Mpc. We find that the enhancement factor, the ratio of $\Delta^2_\nu(k)$ to its linear-response equivalent, is most strongly correlated with $\Omega_{\nu,0} h^2$, and also with the clustering amplitude $\sigma_8$. Furthermore, non-linearities enhance the free-streaming-limit scaling $\partial \log(\Delta^2_\nu / \Delta^2_{\rm m}) / \partial \log(M_\nu)$ beyond its linear value of 4, increasing the $M_\nu$-sensitivity of the small-scale neutrino density., Comment: 17 pages, 14 figures, 3 tables. Emulator code available at: https://github.com/upadhye/Cosmic-Enu

A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological parameters fit to the cosmic microwave background (CMB). The origin of this '$S_8$ tension' remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in the observational measurements and/or uncertainties in the role that baryons play. Here we carefully examine the latter possibility using the new FLAMINGO suite of large-volume cosmological hydrodynamical simulations. We project the simulations onto observable harmonic space and compare with observational measurements of the power and cross-power spectra of cosmic shear, CMB lensing, and the thermal Sunyaev-Zel'dovich (tSZ) effect. We explore the dependence of the predictions on box size and resolution, cosmological parameters including the neutrino mass, and the efficiency and nature of baryonic 'feedback'. Despite the wide range of astrophysical behaviours simulated, we find that baryonic effects are not sufficiently large to remove the $S_8$ tension. Consistent with recent studies, we find the CMB lensing power spectrum is in excellent agreement with the standard model, whilst the cosmic shear power spectrum, tSZ effect power spectrum, and the cross-spectra between shear, CMB lensing, and the tSZ effect are all in varying degrees of tension with the CMB-specified standard model. These results suggest that some mechanism is required to slow the growth of fluctuations at late times and/or on non-linear scales, but that it is unlikely that baryon physics is driving this modification., Comment: 26 pages, 12 figures, MNRAS, accepted with minor revisions

Weak lensing of the cosmic microwave background is rapidly emerging as a powerful probe of neutrinos, dark energy, and new physics. We present a fast computation of the non-linear CMB lensing power spectrum which combines non-linear perturbation theory at early times with power spectrum emulation using cosmological simulations at late times. Comparing our calculation with lightcones from the FLAMINGO 5.6 Gpc cube dark-matter-only simulation, we confirm its accuracy to 1% (2%) up to multipoles L = 3000 (L = 5000) for a nuLambdaCDM cosmology consistent with current data. Clustering suppression due to small-scale baryonic phenomena such as feedback from active galactic nuclei can reduce the lensing power by of order 10%. To our perturbation theory and emulator-based calculation we add SP(k), a new fitting function for this suppression, and confirm its accuracy compared to the FLAMINGO hydrodynamic simulations to 4% at L = 5000, with similar accuracy for massive neutrino models. We further demonstrate that scale-dependent suppression due to neutrinos and baryons approximately factorize, implying that a careful treatment of baryonic feedback can limit biasing neutrino mass constraints., Comment: 14 pages, 21 figures, 1 table

Upcoming large-scale structure surveys will measure the matter power spectrum to approximately percent level accuracy with the aim of searching for evidence for new physics beyond the standard model of cosmology. In order to avoid biasing our conclusions, the theoretical predictions need to be at least as accurate as the measurements for a given choice of cosmological parameters. However, recent theoretical work has shown that complex physical processes associated with galaxy formation (particularly energetic feedback processes associated with stars and especially supermassive black holes) can alter the predictions by many times larger than the required accuracy. Here we present $\texttt{SP(k)}$, a model for the effects of baryon physics on the non-linear matter power spectrum based on a new large suite of hydrodynamical simulations. Specifically, the ANTILLES suite consists of 400 simulations spanning a very wide range of the "feedback landscape" and show that the effects of baryons on the matter power spectrum can be understood at approaching the percent level in terms of the mean baryon fraction of haloes, at scales of up to $k \lesssim 10 \, h \, $Mpc$^{-1}$ and redshifts up to $z=3$. For the range of scales and redshifts that will be probed by forthcoming cosmic shear measurements, most of the effects are driven by galaxy group-mass haloes ($M \sim 10^{13-14}$ M$_\odot$). We present a simple Python implementation of our model, available at $\href{https://github.com/jemme07/pyspk}{\mathrm{https{:}//github.com/jemme07/pyspk}}$, which can be used to incorporate baryon effects in standard gravity-only predictions, allowing for marginalisation over baryon physics within cosmological pipelines., Comment: Refereed version accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS). 16 pages, 11 figures

Velocity dispersion of the massive neutrinos presents a daunting challenge for non-linear cosmological perturbation theory. We consider the neutrino population as a collection of non-linear fluids, each with uniform initial momentum, through an extension of the Time Renormalization Group perturbation theory. Employing recently-developed Fast Fourier Transform techniques, we accelerate our non-linear perturbation theory by more than two orders of magnitude, making it quick enough for practical use. After verifying that the neutrino mode-coupling integrals and power spectra converge, we show that our perturbation theory agrees with N-body neutrino simulations to within 10% for neutrino fractions $\Omega_{\nu,0} h^2 \leq 0.005$ up to wave numbers of k = 1 h/Mpc, an accuracy consistent with 2.5% errors in the neutrino mass determination. Non-linear growth represents a >10% correction to the neutrino power spectrum even for density fractions as low as $\Omega_{\nu,0} h^2 = 0.001$, demonstrating the limits of linear theory for accurate neutrino power spectrum predictions. Our code FlowsForTheMasses is avaliable online at github.com/upadhye/FlowsForTheMasses ., Comment: 51 pages, 17 figures. Matches accepted version. Code available at github.com/upadhye/FlowsForTheMasses

Simulation of the cosmic clustering of massive neutrinos is a daunting task, due both to their large velocity dispersion and to their weak clustering power becoming swamped by Poisson shot noise. We present a new approach, the multi-fluid hybrid-neutrino simulation, which partitions the neutrino population into multiple flows, each of which is characterised by its initial momentum and treated as a separate fluid. These fluid flows respond initially linearly to nonlinear perturbations in the cold matter, but slowest flows are later converted to a particle realisation should their clustering power exceed some threshold. After outlining the multi-fluid description of neutrinos, we study the conversion of the individual flows into particles, in order to quantify transient errors, as well as to determine a set of criteria for particle conversion. Assembling our results into a total neutrino power spectrum, we demonstrate that our multi-fluid hybrid-neutrino simulation is convergent to $<3\%$ if conversion happens at $z=19$ and agrees with more expensive simulations in the literature for neutrino fractions as high as $\Omega_\nu h^2 = 0.005$. Moreover, our hybrid-neutrino approach retains fine-grained information about the neutrinos' momentum distribution. However, the momentum resolution is currently limited by free-streaming transients excited by missing information in the neutrino particle initialisation procedure, which restricts the particle conversion to z $\gtrsim 19$ if percent-level resolution is desired., Comment: 43 pages, 19 figures; v2: 45 pages, added references, extended discussion on timing, matches version accepted by JCAP

Modern cosmological surveys are delivering datasets characterized by unprecedented quality and statistical completeness; this trend is expected to continue into the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the nonlinear regime of structure formation where cosmological simulations are the primary means of obtaining the required information. The computational cost of sufficiently resolved large-volume simulations makes it prohibitive to run very large ensembles. Nevertheless, precision emulators built on a tractable number of high-quality simulations can be used to build very fast prediction schemes to enable a variety of cosmological inference studies. We have recently introduced the Mira-Titan Universe simulation suite designed to construct emulators for a range of cosmological probes. The suite covers the standard six cosmological parameters $\{\omega_m,\omega_b, \sigma_8, h, n_s, w_0\}$ and, in addition, includes massive neutrinos and a dynamical dark energy equation of state, $\{\omega_{\nu}, w_a\}$. In this paper we present the final emulator for the matter power spectrum based on 111 cosmological simulations, each covering a (2.1Gpc)$^3$ volume and evolving 3200$^3$ particles. An additional set of 1776 lower-resolution simulations and TimeRG perturbation theory results for the power spectrum are used to cover scales straddling the linear to mildly nonlinear regimes. The emulator provides predictions at the two to three percent level of accuracy over a wide range of cosmological parameters and is publicly released as part of this paper.

We present in this work a novel and yet extremely simple method for incorporating the effects of massive neutrinos in cosmological $N$-body simulations. This so-called "SuperEasy linear response" approach is based upon analytical solutions to the collisionless Boltzmann equation in the clustering and free-streaming limits, which are then connected by a rational function interpolation function with cosmology-dependent coefficients given by simple algebraic expressions of the cosmological model parameters. The outcome is a {\it one-line modification} to the gravitational potential that requires only the cold matter density contrast as a real-time input, and that can be incorporated into any $N$-body code with a Particle--Mesh component with no additional implementation cost. To demonstrate its power, we implement the SuperEasy method in the publicly available \gadgetcode{} code, and show that for neutrino mass sums not exceeding $\sum m_\nu \simeq 1$ eV, the total matter and cold matter power spectra are in sub-1\% and sub-0.1\% agreement with those from state-of-the-art linear response simulations in literature. Aside from its minimal implementation cost, compared with existing massive neutrino simulation methods, the SuperEasy approach has better memory efficiency, incurs no runtime overhead relative to a standard $\Lambda$CDM simulation, and requires no post-processing. The minimal nature of the method allows limited computational resources to be diverted to modelling other physical effects of interest, e.g., baryonic physics via hydrodynamics., Comment: 31 pages, 11 figures, v2 corresponds to the published version in JCAP

A significant challenge for modelling the massive neutrino as a hot dark matter is its large velocity dispersion. In this work, we investigate and implement a multi-fluid perturbation theory that treats the cosmic neutrino population as a collection of fluids with a broad range of bulk velocities. These fluids respond linearly to the clustering of cold matter, which may be linear and described by standard linear perturbation theory, or non-linear, described using either higher-order perturbation theory or N-body simulations. We verify that such an alternative treatment of neutrino perturbations agrees closely with state-of-the-art neutrino linear response calculations in terms of power spectrum and bispectrum predictions. Combining multi-fluid neutrino linear response with a non-linear calculation for the cold matter clustering, we find for a reference nuLambdaCDM cosmology with neutrino mass sum of 0.93 eV an enhancement of the small-scale neutrino power by an order of magnitude relative to a purely linear calculation. The corresponding clustering enhancement in the cold matter, however, is a modest ~0.05%. Importantly, our multi-fluid approach uniquely enables us to identify that the slowest-moving 25% of the neutrino population clusters strongly enough to warrant a non-linear treatment. Such a precise calculation of neutrino clustering on small scales accompanied by fine-grained velocity information would be invaluable for experiments such as PTOLEMY that probe the local neutrino density and velocity in the solar neighbourhood., Comment: 34 pages, 15 figures, 2 tables. Matches accepted version. Code available at https://github.com/joechenUNSW/gadget4-nu_lr

Numerical simulations of massive neutrino cosmologies consistently find a spoon-like feature in the non-linear matter power spectrum ratios of cosmological models that differ only in the neutrino mass fraction f_N. Typically, the ratio approaches unity at low wave numbers k, decreases by ~ 10 f_N at k ~ 1 h/Mpc, and turns up again at large k. Using the halo model of large-scale structure, we show that this spoon feature originates in the transition from the two-halo power spectrum to the one-halo power spectrum. The former's sensitivity to f_N rises with k, while that of the latter decreases with k. The presence of this spoon feature is robust with respect to different choices of the halo mass function and the halo density profile, and does not require any parameter tuning within the halo model. We demonstrate that a standard halo model calculation is already able to predict the depth, width, and position of this spoon as well as its evolution with redshift z with remarkable accuracy. Predictions at z >= 1 can be further improved using non-linear perturbative inputs., Comment: 24 pages, 9 figures (matches accepted version)

The presence of massive neutrinos affects the growth of large-scale structure in the universe, leaving a potentially observable imprint on the abundance and properties of massive dark matter-dominated halos. Cosmological surveys detect large numbers of these halos in the form of rich groups and clusters, using the information as an input to constraining the properties of dark energy. We use a suite of N-body simulations that include the effects of massive neutrinos as well as of dynamical dark energy to study the properties of the mass function. As in our previous work, we follow an approach valid at low neutrino mass, where the neutrino overdensities are assumed to be too small to act as a significant nonlinear source term for gravity. We study how well a universal form for the halo mass function describes our numerical results, finding that the use of an appropriate linear power spectrum within the formalism yields a good match to the simulation results, correctly accounting for the (neutrino mass-dependent) suppression of the mass function., Comment: 11 pages, 5 figures, to be submitted to PRD

It is known that if the Peccei-Quinn symmetry breaking field is displaced from its minimum during inflation, the axion isocurvature spectrum is generically strongly blue tilted with a break transition to a flat spectrum. We fit this spectrum (incorporated into the "vanilla" $\Lambda$-CDM cosmological model) to the Planck and BOSS DR11 data and find a mild hint for the presence of axionic blue-tilted isocurvature perturbations. We find the best fit parameter region is consistent with all of the dark matter being composed of QCD axions in the context of inflationary cosmology with an expansion rate of order $10^{8}$ GeV, the axion decay constant of order $10^{13}$ GeV, and the initial misalignment angle of order unity. Intriguingly, isocurvature with a spectral break may at least partially explain the low-$\ell$ vs. high-$\ell$ anomalies seen in the CMB data., Comment: 25 pages

Cosmology in the near future promises a measurement of the sum of neutrino masses, a fundamental Standard Model parameter, as well as substantially-improved constraints on the dark energy. We use the shape of the BOSS redshift-space galaxy power spectrum, in combination with CMB and supernova data, to constrain the neutrino masses and the dark energy. Essential to this calculation are several recent advances in non-linear cosmological perturbation theory, including FFT methods, redshift space distortions, and scale-dependent growth. Our 95% confidence upper bound of 180 meV on the sum of masses degrades substantially to 540 meV when the dark energy equation of state and its first derivative are also allowed to vary, representing a significant challenge to current constraints. We also study the impact of additional galaxy bias parameters, finding that a greater allowed range of scale-dependent bias only slightly shifts the preferred neutrino mass value, weakens its upper bound by about 20%, and has a negligible effect on the other cosmological parameters., Comment: Matches accepted version. Code available at github.com/upadhye/redTime

We introduce a new cosmic emulator for the matter power spectrum covering eight cosmological parameters. Targeted at optical surveys, the emulator provides accurate predictions out to a wavenumber k~5/Mpc and redshift z<=2. Besides covering the standard set of LCDM parameters, massive neutrinos and a dynamical dark energy of state are included. The emulator is built on a sample set of 36 cosmological models, carefully chosen to provide accurate predictions over the wide and large parameter space. For each model, we have performed a high-resolution simulation, augmented with sixteen medium-resolution simulations and TimeRG perturbation theory results to provide accurate coverage of a wide k-range; the dataset generated as part of this project is more than 1.2Pbyte. With the current set of simulated models, we achieve an accuracy of approximately 4%. Because the sampling approach used here has established convergence and error-control properties, follow-on results with more than a hundred cosmological models will soon achieve ~1% accuracy. We compare our approach with other prediction schemes that are based on halo model ideas and remapping approaches. The new emulator code is publicly available., Comment: 12 pages, 12 figures, emulator code publicly available here: https://github.com/lanl/CosmicEmu

Gravity is the weakest fundamental interaction and the only one that has not been measured at the particle level. Traditional experimental methods, from astronomical observations to torsion balances, use macroscopic masses to both source and probe gravitational fields. Matter wave interferometers have used neutrons, atoms and molecular clusters as microscopic test particles, but initially probed the field sourced by the entire earth. Later, the gravitational field arising from hundreds of kilograms of artificial source masses was measured with atom interferometry. Miniaturizing the source mass and moving it into the vacuum chamber could improve positioning accuracy, allow the use of monocrystalline source masses for improved gravitational measurements, and test new physics, such as beyond-standard-model ("fifth") forces of nature and non-classical effects of gravity. In this work, we detect the gravitational force between freely falling cesium atoms and an in-vacuum, centimeter-sized source mass using atom interferometry with state-of-the-art sensitivity. The ability to sense gravitational-strength coupling is conjectured to access a natural lower bound for fundamental forces, thereby representing an important milestone in searches for physics beyond the standard model. A local, in-vacuum source mass is particularly sensitive to a wide class of interactions whose effects would otherwise be suppressed beyond detectability in regions of high matter density. For example, our measurement strengthens limits on a number of cosmologically-motivated scalar field models, such as chameleon and symmetron fields, by over two orders of magnitude and paves the way toward novel measurements of Newton's gravitational constant G and the gravitational Aharonov-Bohm effect

Blue axion isocurvature perturbations are both theoretically well-motivated and interesting from a detectability perspective. These power spectra generically have a break from the blue region to a flat region. Previous investigations of the power spectra were analytic, which left a gap in the predicted spectrum in the break region due to the non-applicability of the used analytic techniques. We therefore compute the isocurvature spectrum numerically for an explicit supersymmetric axion model. We find a bump that enhances the isocurvature signal for this class of scenarios. A fitting function of three parameters is constructed that fits the spectrum well for the particular axion model we study. This fitting function should be useful for blue isocurvature signal hunting in data and making experimental sensitivity forecasts., Comment: 28 pages and 7 figures

Large-scale structure surveys in the coming years will measure the redshift-space power spectrum to unprecedented accuracy, allowing for powerful new tests of the LambdaCDM picture as well as measurements of particle physics parameters such as the neutrino masses. We extend the Time-RG perturbative framework to redshift space, computing the power spectrum P_s(k,mu) in massive neutrino cosmologies with time-dependent dark energy equations of state w(z). Time-RG is uniquely capable of incorporating scale-dependent growth into the P_s(k,mu) computation, which is important for massive neutrinos as well as modified gravity models. Although changes to w(z) and the neutrino mass fraction both affect the late-time scale-dependence of the non-linear power spectrum, we find that the two effects depend differently on the line-of-sight angle mu. Finally, we use the HACC N-body code to quantify errors in the perturbative calculations. For a LambdaCDM model at redshift z=1, our procedure predicts the monopole~(quadrupole) to 1% accuracy up to a wave number 0.19h/Mpc (0.28h/Mpc), compared to 0.08h/Mpc (0.07h/Mpc) for the Kaiser approximation and 0.19h/Mpc (0.16h/Mpc) for the current state-of-the-art perturbation scheme. Our calculation agrees with the simulated redshift-space power spectrum even for neutrino masses above the current bound, and for rapidly-evolving dark energy equations of state, |dw/dz| ~ 1. Along with this article, we make our redshift-space Time-RG implementation publicly available as the code redTime., Comment: 18 pages, 17 figures, 4 tables. Matches version accepted by PRD. redTime code available at http://www.hep.anl.gov/cosmology/pert.html

Over the next decade, cosmological measurements of the large-scale structure of the Universe will be sensitive to the combined effects of dynamical dark energy and massive neutrinos. The matter power spectrum is a key repository of this information. We extend higher-order perturbative methods for computing the power spectrum to investigate these effects over quasi-linear scales. Through comparison with N-body simulations we establish the regime of validity of a Time-Renormalization Group (Time-RG) perturbative treatment that includes dynamical dark energy and massive neutrinos. We also quantify the accuracy of Standard (SPT), Renormalized (RPT) and Lagrangian Resummation (LPT) perturbation theories without massive neutrinos. We find that an approximation that neglects neutrino clustering as a source for nonlinear matter clustering predicts the Baryon Acoustic Oscillation (BAO) peak position to 0.25% accuracy for redshifts 1 < z < 3, justifying the use of LPT for BAO reconstruction in upcoming surveys. We release a modified version of the public Copter code which includes the additional physics discussed in the paper., Comment: Code available at http://www.hep.anl.gov/cosmology/pert.html . 16 pages, 10 figures. Matches version accepted by PRD

We present data from our investigation of the anomalous orange-colored afterglow that was seen in the GammeV Chameleon Afterglow Search (CHASE). These data includes information about the broad band color of the observed glow, the relationship between the glow and the temperature of the apparatus, and other data taken prior to and during the science operations of CHASE. While differing in several details, the generic properties of the afterglow from CHASE are similar to luminescence seen in some vacuum compounds. Contamination from this, or similar, luminescent signatures will likely impact the design of implementation of future experiments involving single photon detectors and high intensity light sources in a cryogenic environment., Comment: 6 pages, 5 figures, submitted to PRX