Showing 1–50 of 69 results for author: Wong, Y Y Y

Towards a precision calculation of $N_{\rm eff}$ in the Standard Model III: Improved estimate of NLO contributions to the collision integral

Authors: Marco Drewes, Yannis Georis, Michael Klasen, Luca Paolo Wiggering, Yvonne Y. Y. Wong

Abstract: We compute the dominant QED correction to the neutrino-electron interaction rate in the vicinity of neutrino decoupling in the early universe, and estimate its impact on the effective number of neutrino species $N_{\rm eff}$ in cosmic microwave background anisotropy observations. We find that the correction to the interaction rate is at the sub-percent level, consistent with a recent estimate by J… ▽ More We compute the dominant QED correction to the neutrino-electron interaction rate in the vicinity of neutrino decoupling in the early universe, and estimate its impact on the effective number of neutrino species $N_{\rm eff}$ in cosmic microwave background anisotropy observations. We find that the correction to the interaction rate is at the sub-percent level, consistent with a recent estimate by Jackson and Laine. Relative to that work we include the electron mass in our computations, but restrict our analysis to the enhanced $t$-channel contributions. The fractional change in $N_{\rm eff}^{\rm SM}$ due to the rate correction is of order $10^{-5}$ or below, i.e., about a factor of 30 smaller than that recently claimed by Cielo {\it et al.}, and below the nominal computational uncertainties of the current benchmark value of $N_{\rm eff}^{\rm SM} = 3.0440 \pm 0.0002$. We therefore conclude that aforementioned number remains to be the state-of-the-art benchmark for $N_{\rm eff}^{\rm SM}$ in the standard model of particle physics. △ Less

Submitted 14 June, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 34 pages, 8 figures, 2 tables; v2: version accepted by JCAP

Report number: CPPC-2024-01, MS-TP-24-06

Journal ref: JCAP 06 (2024) 032

Cosmic-Enu: An emulator for the non-linear neutrino power spectrum

Authors: Amol Upadhye, Juliana Kwan, Ian G. McCarthy, Jaime Salcido, Kelly R. Moran, Earl Lawrence, Yvonne Y. Y. Wong

Abstract: Cosmology is poised to measure the neutrino mass sum $M_ν$ 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, $Δ^2_ν(k)$, agrees with its companion N-body simulation at the $10\%-15\%$ l… ▽ More Cosmology is poised to measure the neutrino mass sum $M_ν$ 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, $Δ^2_ν(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 $Δ^2_ν(k)$ covering a large range of cosmological parameters and neutrino fractions $Ω_{ν,0} h^2 \leq 0.01$, which corresponds to $M_ν\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_ν=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 $Δ^2_ν(k)$ to its linear-response equivalent, is most strongly correlated with $Ω_{ν,0} h^2$, and also with the clustering amplitude $σ_8$. Furthermore, non-linearities enhance the free-streaming-limit scaling $\partial \log(Δ^2_ν/ Δ^2_{\rm m}) / \partial \log(M_ν)$ beyond its linear value of 4, increasing the $M_ν$-sensitivity of the small-scale neutrino density. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 17 pages, 14 figures, 3 tables. Emulator code available at: https://github.com/upadhye/Cosmic-Enu

Revisiting Cosmological Constraints on Supersymmetric SuperWIMPs

Authors: Meera Deshpande, Jan Hamann, Dipan Sengupta, Martin White, Anthony G. Williams, Yvonne Y. Y. Wong

Abstract: SuperWIMPs are extremely weakly interacting massive particles that inherit their relic abundance from late decays of frozen-out parent particles. Within supersymmetric models, gravitinos and axinos represent two of the most well-motivated superWIMPs. In this paper we revisit constraints on these scenarios from a variety of cosmological observations that probe their production mechanisms as well as… ▽ More SuperWIMPs are extremely weakly interacting massive particles that inherit their relic abundance from late decays of frozen-out parent particles. Within supersymmetric models, gravitinos and axinos represent two of the most well-motivated superWIMPs. In this paper we revisit constraints on these scenarios from a variety of cosmological observations that probe their production mechanisms as well as the superWIMP kinematic properties in the early Universe. We consider in particular observables of Big Bang Nucleosynthesis and the Cosmic Microwave Background (spectral distortion and anisotropies), which limit the fractional energy injection from the late decays, as well as warm and mixed dark matter constraints derived from the Lyman-$α$ forest and other small-scale structure observables. We discuss complementary constraints from collider experiments, and argue that cosmological considerations rule out a significant part of the gravitino and the axino superWIMP parameter space. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 17 Pages, 7 figures

Axion minivoids and implications for direct detection

Authors: Benedikt Eggemeier, Ciaran A. J. O'Hare, Giovanni Pierobon, Javier Redondo, Yvonne Y. Y. Wong

Abstract: In the scenario in which QCD axion dark matter is produced after inflation, the Universe is populated by large inhomogeneities on very small scales. Eventually, these fluctuations will collapse gravitationally to form dense axion miniclusters that trap up to $\sim$75% of the dark matter within asteroid-mass clumps. Axion miniclusters are physically tiny however, so haloscope experiments searching… ▽ More In the scenario in which QCD axion dark matter is produced after inflation, the Universe is populated by large inhomogeneities on very small scales. Eventually, these fluctuations will collapse gravitationally to form dense axion miniclusters that trap up to $\sim$75% of the dark matter within asteroid-mass clumps. Axion miniclusters are physically tiny however, so haloscope experiments searching for axions directly on Earth are much more likely to be probing ``minivoids'' -- the space in between miniclusters. This scenario seems like it ought to spell doom for haloscopes, but while these minivoids might be underdense, they are not totally devoid of axions. Using Schrödinger-Poisson and N-body simulations to evolve from realistic initial field configurations, we quantify the extent to which the local ambient dark matter density is suppressed in the post-inflationary scenario. We find that a typical experimental measurement will sample an axion density that is only around 10% of the expected galactic dark matter density. Our results also have implications for experimental campaigns lasting longer than a few years, as well as broadband haloscopes that have sensitivity to transient signatures. We show that for a $\mathcal{O}$(year)-long integration times, the measured dark matter density should be expected to vary by 20--30%. △ Less

Submitted 9 May, 2023; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 11 pages + appendix, 10 figures. Version accepted by PRD

Snowmass Neutrino Frontier Report

Authors: Patrick Huber, Kate Scholberg, Elizabeth Worcester, Jonathan Asaadi, A. Baha Balantekin, Nathaniel Bowden, Pilar Coloma, Peter B. Denton, André de Gouvêa, Laura Fields, Megan Friend, Steven Gardiner, Carlo Giunti, Julieta Gruszko, Benjamin J. P. Jones, Georgia Karagiorgi, Lisa Kaufman, Joshua R. Klein, Lisa W. Koerner, Yusuke Koshio, Jonathan M. Link, Bryce R. Littlejohn, Ana A. Machado, Pedro A. N. Machado, Kendall Mahn , et al. (34 additional authors not shown)

Abstract: This report summarizes the current status of neutrino physics and the broad and exciting future prospects identified for the Neutrino Frontier as part of the 2021 Snowmass Process. This report summarizes the current status of neutrino physics and the broad and exciting future prospects identified for the Neutrino Frontier as part of the 2021 Snowmass Process. △ Less

Submitted 8 December, 2022; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: 49 pages, contribution to: 2021 Snowmass Summer Study. Minor updates

Flows For The Masses: A multi-fluid non-linear perturbation theory for massive neutrinos

Authors: Joe Zhiyu Chen, Amol Upadhye, Yvonne Y. Y. Wong

Abstract: 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-line… ▽ More 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 $Ω_{ν,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 $Ω_{ν,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 . △ Less

Submitted 16 May, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 51 pages, 17 figures. Matches accepted version. Code available at github.com/upadhye/FlowsForTheMasses

Hybrid multi-fluid-particle simulations of the cosmic neutrino background

Authors: Joe Zhiyu Chen, Markus R. Mosbech, Amol Upadhye, Yvonne Y. Y. Wong

Abstract: 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… ▽ More 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 $Ω_ν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. △ Less

Submitted 14 February, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 43 pages, 19 figures; v2: 45 pages, added references, extended discussion on timing, matches version accepted by JCAP

Report number: CPPC-2022-10

Gravitational-wave event rates as a new probe for dark matter microphysics

Authors: Markus R. Mosbech, Alexander C. Jenkins, Sownak Bose, Celine Boehm, Mairi Sakellariadou, Yvonne Y. Y. Wong

Abstract: We show that gravitational waves have the potential to unravel the microphysical properties of dark matter due to the dependence of the binary black hole merger rate on cosmic structure formation, which is itself highly dependent on the dark matter scenario. In particular, we demonstrate that suppression of small-scale structure -- such as that caused by interacting, warm, or fuzzy dark matter --… ▽ More We show that gravitational waves have the potential to unravel the microphysical properties of dark matter due to the dependence of the binary black hole merger rate on cosmic structure formation, which is itself highly dependent on the dark matter scenario. In particular, we demonstrate that suppression of small-scale structure -- such as that caused by interacting, warm, or fuzzy dark matter -- leads to a significant reduction in the rate of binary black hole mergers at redshifts $z\gtrsim5$. This shows that future gravitational-wave observations will provide a new probe of physics beyond the $Λ$CDM model. △ Less

Submitted 18 July, 2023; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 15 pages, 12 figures, matches version published in PRD

Report number: KCL-PH-TH-2022-42

Probing dark matter interactions with 21cm observations

Authors: Markus R. Mosbech, Celine Boehm, Yvonne Y. Y. Wong

Abstract: Similarly to warm dark matter which features a cut-off in the matter power spectrum due to free-streaming, many interacting dark matter models predict a suppression of the matter power spectrum on small length scales through collisional damping. Forecasts for 21cm line intensity mapping have shown that an instrument like the SKA will be able to probe a suppression of power in warm dark matter scen… ▽ More Similarly to warm dark matter which features a cut-off in the matter power spectrum due to free-streaming, many interacting dark matter models predict a suppression of the matter power spectrum on small length scales through collisional damping. Forecasts for 21cm line intensity mapping have shown that an instrument like the SKA will be able to probe a suppression of power in warm dark matter scenarios in a statistically significant way. Here we investigate the implications of these findings on interacting dark matter scenarios, particularly dark matter-neutrino interactions, which we use as an example. Using a suite of cosmological $N$-body simulations, we demonstrate that interacting scenarios show a suppression of the non-linear power spectrum similar to warm dark matter models. This implies that 21cm line intensity mapping will be able to set the strongest limits yet on dark matter-neutrino scattering, improving the constraints by two orders of magnitude over current Lyman-$α$ bounds, and by four orders of magnitude over cosmic microwave background and baryon acoustic oscillations limits. However, to distinguish between warm dark matter and interacting scenarios, our simulations show that percent-level precision measurements of the matter power spectrum at redshifts $z\gtrsim15$ are necessary, as the key features of interacting scenarios are washed out by non-linear evolution at later times. △ Less

Submitted 26 February, 2023; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 21 pages, 7 figures. Version accepted for publication

Weaker yet again: mass spectrum-consistent cosmological constraints on the neutrino lifetime

Authors: Joe Zhiyu Chen, Isabel M. Oldengott, Giovanni Pierobon, Yvonne Y. Y. Wong

Abstract: We consider invisible neutrino decay $ν_H \to ν_l + φ$ in the ultra-relativistic limit and compute the neutrino anisotropy loss rate relevant for the cosmic microwave background (CMB) anisotropies. Improving on our previous work which assumed massless $ν_l$ and $φ$, we reinstate in this work the daughter neutrino mass $m_{νl}$ in a manner consistent with the experimentally determined neutrino mass… ▽ More We consider invisible neutrino decay $ν_H \to ν_l + φ$ in the ultra-relativistic limit and compute the neutrino anisotropy loss rate relevant for the cosmic microwave background (CMB) anisotropies. Improving on our previous work which assumed massless $ν_l$ and $φ$, we reinstate in this work the daughter neutrino mass $m_{νl}$ in a manner consistent with the experimentally determined neutrino mass splittings. We find that a nonzero $m_{νl}$ introduces a new phase space factor in the loss rate $Γ_{\rm T}$ proportional to $(Δm_ν^2/m_{ν_H}^2)^2$ in the limit of a small squared mass gap between the parent and daughter neutrinos, i.e., $Γ_{\rm T} \sim (Δm_ν^2/m_{νH}^2)^2 (m_{νH}/E_ν)^5 (1/τ_0)$, where $τ_0$ is the $ν_H$ rest-frame lifetime. Using a general form of this result, we update the limit on $τ_0$ using the Planck 2018 CMB data. We find that for a parent neutrino of mass $m_{νH} \lesssim 0.1 {\rm eV}$, the new phase space factor weakens the constraint on its lifetime by up to a factor of 50 if $Δm_ν^2$ corresponds to the atmospheric mass gap and up to $10^{5}$ if the solar mass gap, in comparison with naive estimates that assume $m_{νl}=0$. The revised constraints are (i) $τ^0 \gtrsim (6 \to 10) \times 10^5~{\rm s}$ and $τ^0 \gtrsim (400 \to 500)~{\rm s}$ if only one neutrino decays to a daughter neutrino separated by, respectively, the atmospheric and the solar mass gap, and (ii) $τ^0 \gtrsim (2 \to 3) \times 10^7~{\rm s}$ in the case of two decay channels with one near-common atmospheric mass gap. In contrast to previous, naive limits which scale as $m_{νH}^5$, these mass spectrum-consistent $τ_0$ constraints are remarkably independent of the parent mass and open up a swath of parameter space within the projected reach of IceCube and other neutrino telescopes in the next two decades. △ Less

Submitted 16 June, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 23 pages, 6 figures. v3 accepted for publication by EPJC. Updated the paragraph explaining the separable ansatz on page 5 and added some extra text and equations in appendix B for the Legendre polynomials expansion. Minor additions and small correction to figure 3 bounds

White Paper on Light Sterile Neutrino Searches and Related Phenomenology

Authors: M. A. Acero, C. A. Argüelles, M. Hostert, D. Kalra, G. Karagiorgi, K. J. Kelly, B. Littlejohn, P. Machado, W. Pettus, M. Toups, M. Ross-Lonergan, A. Sousa, P. T. Surukuchi, Y. Y. Y. Wong, W. Abdallah, A. M. Abdullahi, R. Akutsu, L. Alvarez-Ruso, D. S. M. Alves, A. Aurisano, A. B. Balantekin, J. M. Berryman, T. Bertólez-Martínez, J. Brunner, M. Blennow , et al. (147 additional authors not shown)

Abstract: This white paper provides a comprehensive review of our present understanding of experimental neutrino anomalies that remain unresolved, charting the progress achieved over the last decade at the experimental and phenomenological level, and sets the stage for future programmatic prospects in addressing those anomalies. It is purposed to serve as a guiding and motivational "encyclopedic" reference,… ▽ More This white paper provides a comprehensive review of our present understanding of experimental neutrino anomalies that remain unresolved, charting the progress achieved over the last decade at the experimental and phenomenological level, and sets the stage for future programmatic prospects in addressing those anomalies. It is purposed to serve as a guiding and motivational "encyclopedic" reference, with emphasis on needs and options for future exploration that may lead to the ultimate resolution of the anomalies. We see the main experimental, analysis, and theory-driven thrusts that will be essential to achieving this goal being: 1) Cover all anomaly sectors -- given the unresolved nature of all four canonical anomalies, it is imperative to support all pillars of a diverse experimental portfolio, source, reactor, decay-at-rest, decay-in-flight, and other methods/sources, to provide complementary probes of and increased precision for new physics explanations; 2) Pursue diverse signatures -- it is imperative that experiments make design and analysis choices that maximize sensitivity to as broad an array of these potential new physics signatures as possible; 3) Deepen theoretical engagement -- priority in the theory community should be placed on development of standard and beyond standard models relevant to all four short-baseline anomalies and the development of tools for efficient tests of these models with existing and future experimental datasets; 4) Openly share data -- Fluid communication between the experimental and theory communities will be required, which implies that both experimental data releases and theoretical calculations should be publicly available; and 5) Apply robust analysis techniques -- Appropriate statistical treatment is crucial to assess the compatibility of data sets within the context of any given model. △ Less

Submitted 17 May, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021 by the NF02 Topical Group (Understanding Experimental Neutrino Anomalies). Submitted to J. Phys. G as a Major Report

Simulations of axion-like particles in the post-inflationary scenario

Authors: Ciaran A. J. O'Hare, Giovanni Pierobon, Javier Redondo, Yvonne Y. Y. Wong

Abstract: Axions and axion-like particles (ALPs) are some of the most popular candidates for dark matter, with several viable production scenarios that make different predictions. In the scenario in which the axion is born after inflation, its field develops significant inhomogeneity and evolves in a highly nonlinear fashion. Understanding the eventual abundance and distribution of axionic dark matter in th… ▽ More Axions and axion-like particles (ALPs) are some of the most popular candidates for dark matter, with several viable production scenarios that make different predictions. In the scenario in which the axion is born after inflation, its field develops significant inhomogeneity and evolves in a highly nonlinear fashion. Understanding the eventual abundance and distribution of axionic dark matter in this scenario therefore requires dedicated numerical simulations. So far the community has focused its efforts on simulations of the QCD axion, a model that predicts a specific temperature dependence for the axion mass. Here, we go beyond the QCD axion, and perform a suite of simulations over a range of possible temperature dependencies labelled by a power-law index. We study the complex dynamics of the axion field, including the scaling of cosmic strings and domain walls, the spectrum of non-relativistic axions, the lifetime and internal structure of axitons, and the seeds of miniclusters. In particular, we quantify how much the string-wall network contributes to the dark matter abundance as a function of how quickly the axion mass grows. We find that a temperature-independent model produces 25\% more dark matter than the standard misalignment calculation. In contrast to this generic ALP, QCD axion models are almost six times less efficient at producing dark matter. Given the flourishing experimental campaign to search for ALPs, these results have potentially wide implications for direct and indirect searches. △ Less

Submitted 30 March, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: Version accepted by PRD, added small paragraph (section IIG) on the comparison with previous simulations. 19 pages, additional visualisations at https://www.youtube.com/channel/UCqosWwoOKrdmNKBbIwVUuGg

Journal ref: Phys. Rev. D 105, 055025 (2022)

Towards a precision calculation of $N_{\rm eff}$ in the Standard Model II: Neutrino decoupling in the presence of flavour oscillations and finite-temperature QED

Authors: Jack J. Bennett, Gilles Buldgen, Pablo F. de Salas, Marco Drewes, Stefano Gariazzo, Sergio Pastor, Yvonne Y. Y. Wong

Abstract: We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, $N_{\rm eff}^{\rm SM}$, that quantifies the cosmological neutrino-to-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to ${\cal O}(e^3)$, where $e$ is the elementary electric… ▽ More We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, $N_{\rm eff}^{\rm SM}$, that quantifies the cosmological neutrino-to-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to ${\cal O}(e^3)$, where $e$ is the elementary electric charge, and a full evaluation of the neutrino--neutrino collision integral. We provide furthermore a detailed assessment of the uncertainties in the benchmark $N_{\rm eff}^{\rm SM}$ value, through testing the value's dependence on (i)~optional approximate modelling of the weak collision integrals, (ii)~measurement errors in the physical parameters of the weak sector, and (iii)~numerical convergence, particularly in relation to momentum discretisation. Our new, recommended standard-model benchmark is $N_{\rm eff}^{\rm SM} = 3.0440 \pm 0.0002$, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure ($|δN_{\rm eff}| \sim10^{-4}$), augmented by measurement errors in the solar mixing angle $\sin^2θ_{12}$ ($|δN_{\rm eff}| \sim10^{-4}$). △ Less

Submitted 19 March, 2021; v1 submitted 4 December, 2020; originally announced December 2020.

Comments: 31 pages, 7 figures. Version accepted for publication in JCAP

Report number: CPPC-2020-10

Journal ref: JCAP 04 (2021) 073

One line to run them all: SuperEasy massive neutrino linear response in $N$-body simulations

Authors: Joe Zhiyu Chen, Amol Upadhye, Yvonne Y. Y. Wong

Abstract: 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… ▽ More 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_ν\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 $Λ$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. △ Less

Submitted 2 May, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: 31 pages, 11 figures, v2 corresponds to the published version in JCAP

Report number: CPPC-2020-07

Journal ref: JCAP04(2021)078

The cosmic neutrino background as a collection of fluids in large-scale structure simulations

Authors: Joe Zhiyu Chen, Amol Upadhye, Yvonne Y. Y. Wong

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: 34 pages, 15 figures, 2 tables. Matches accepted version. Code available at https://github.com/joechenUNSW/gadget4-nu_lr

Report number: CPPC-2020-08

Journal ref: JCAP03(2021)065

The full Boltzmann hierarchy for dark matter-massive neutrino interactions

Authors: Markus R. Mosbech, Celine Boehm, Steen Hannestad, Olga Mena, Julia Stadler, Yvonne Y. Y. Wong

Abstract: The impact of dark matter-neutrino interactions on the measurement of the cosmological parameters has been investigated in the past in the context of massless neutrinos exclusively. Here we revisit the role of a neutrino-dark matter coupling in light of ongoing cosmological tensions by implementing the full Boltzmann hierarchy for three massive neutrinos. Our tightest 95% CL upper limit on the str… ▽ More The impact of dark matter-neutrino interactions on the measurement of the cosmological parameters has been investigated in the past in the context of massless neutrinos exclusively. Here we revisit the role of a neutrino-dark matter coupling in light of ongoing cosmological tensions by implementing the full Boltzmann hierarchy for three massive neutrinos. Our tightest 95% CL upper limit on the strength of the interactions, parameterized via $u_χ=\frac{σ_0}{σ_{Th}}\left(\frac{m_χ}{100 \text{GeV}}\right)^{-1}$, is $u_χ\leq3.34 \cdot 10^{-4}$, arising from a combination of Planck TTTEEE data, Planck lensing data and SDSS BAO data. This upper bound is, as expected, slightly higher than previous results for interacting massless neutrinos, due to the correction factor associated with neutrino masses. We find that these interactions significantly relax the lower bounds on the value of $σ_8$ that is inferred in the context of $Λ$CDM from the Planck data, leading to agreement within 1-2$σ$ with weak lensing estimates of $σ_8$, as those from KiDS-1000. However, the presence of these interactions barely affects the value of the Hubble constant $H_0$. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: 31 pages, 8 figures, 3 tables

Invisible neutrino decay in precision cosmology

Authors: Gabriela Barenboim, Joe Zhiyu Chen, Steen Hannestad, Isabel M. Oldengott, Thomas Tram, Yvonne Y. Y. Wong

Abstract: We revisit the topic of invisible neutrino decay in the precision cosmological context, via a first-principles approach to understanding the cosmic microwave background and large-scale structure phenomenology of such a non-standard physics scenario. Assuming an effective Lagrangian in which a heavier standard-model neutrino $ν_H$ couples to a lighter one $ν_l$ and a massless scalar particle $φ$ vi… ▽ More We revisit the topic of invisible neutrino decay in the precision cosmological context, via a first-principles approach to understanding the cosmic microwave background and large-scale structure phenomenology of such a non-standard physics scenario. Assuming an effective Lagrangian in which a heavier standard-model neutrino $ν_H$ couples to a lighter one $ν_l$ and a massless scalar particle $φ$ via a Yukawa interaction, we derive from first principles the complete set of Boltzmann equations, at both the spatially homogeneous and the first-order inhomogeneous levels, for the phase space densities of $ν_H$, $ν_l$, and $φ$ in the presence of the relevant decay and inverse decay processes. With this set of equations in hand, we perform a critical survey of recent works on cosmological invisible neutrino decay in both limits of decay while $ν_H$ is ultra-relativistic and non-relativistic. Our two main findings are: (i) in the non-relativistic limit, the effective equations of motion used to describe perturbations in the neutrino--scalar system in the existing literature formally violate momentum conservation and gauge invariance, and (ii) in the ultra-relativistic limit, exponential damping of the anisotropic stress does not occur at the commonly-used rate $Γ_{\rm T} =(1/τ_0) (m_{νH}/E_{νH})^3$, but at a rate $\sim (1/τ_0) (m_{νH}/E_{νH})^5$. Both results are model-independent. The impact of the former finding on the cosmology of invisible neutrino decay is likely small. The latter, however, implies a significant revision of the cosmological limit on the neutrino lifetime $τ_0$ from $τ_0^{\rm old} \gtrsim 1.2 \times 10^9\, {\rm s}\, (m_{νH}/50\, {\rm meV})^3$ to $τ_0 \gtrsim (4 \times 10^5 \to 4 \times 10^6)\, {\rm s}\, (m_{νH}/50 \, {\rm meV})^5$. △ Less

Submitted 30 March, 2021; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 53 pages, 10 figures, v2 matches the accepted JCAP version

Spoon or slide? The non-linear matter power spectrum in the presence of massive neutrinos

Authors: Steen Hannestad, Amol Upadhye, Yvonne Y. Y. Wong

Abstract: 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… ▽ More 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. △ Less

Submitted 30 November, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 24 pages, 9 figures (matches accepted version)

Towards a precision calculation of the effective number of neutrinos $N_{\rm eff}$ in the Standard Model I: The QED equation of state

Authors: Jack J. Bennett, Gilles Buldgen, Marco Drewes, Yvonne Y. Y. Wong

Abstract: We revisit several aspects of Standard Model physics at finite temperature that drive the theoretical value of the cosmological parameter $N_{\rm eff}$, the effective number of neutrinos in the early universe, away from 3. Our chief focus is finite-temperature corrections to the equation of state of the QED plasma in the vicinity of neutrino decoupling at $T \sim 1$ MeV, where $T$ is the photon te… ▽ More We revisit several aspects of Standard Model physics at finite temperature that drive the theoretical value of the cosmological parameter $N_{\rm eff}$, the effective number of neutrinos in the early universe, away from 3. Our chief focus is finite-temperature corrections to the equation of state of the QED plasma in the vicinity of neutrino decoupling at $T \sim 1$ MeV, where $T$ is the photon temperature. Working in the instantaneous decoupling approximation, we recover at ${\cal O}(e^2)$, where $e$ is the elementary electric charge, the well-established correction of $δN_{\rm eff}^{(2)} \simeq 0.010$ across a range of plausible neutrino decoupling temperatures, in contrast to an erroneous claim in the recent literature which found twice as large an effect. At ${\cal O}(e^3)$ we find a new and significant correction of $δN_{\rm eff}^{(3)} \simeq -0.001$ that has so far not been accounted for in any precision calculation of $N_{\rm eff}$, significant because this correction is potentially larger than the change in $N_{\rm eff}$ induced between including and excluding neutrino oscillations in the transport modelling. In addition to the QED equation of state, we make a first pass at quantifying finite-temperature QED corrections to the weak interaction rates that directly affect the neutrino decoupling process, and find that the ${\cal O}(e^2)$ thermal electron mass correction induces a change of $δN_{\rm eff}^{m_{\rm th}} \lesssim 10^{-4}$. A complete assessment of the various effects considered in this work on the final value of $N_{\rm eff}$ will necessitate an account of neutrino energy transport beyond the instantaneous decoupling approximation. However, relative to $N_{\rm eff} = 3.044$ obtained in the most recent such calculation, we expect the new effects found in this work to lower the number to $N_{\rm eff} = 3.043$. △ Less

Submitted 30 November, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 25 pages, plus appendices and 8 figures; v3: references added, version accepted by JCAP, minor terminology changes

Fitting functions on the cheap: the relative nonlinear matter power spectrum

Authors: Steen Hannestad, Yvonne Y. Y. Wong

Abstract: We propose an alternative approach to the construction of fitting functions to the nonlinear matter power spectrum extracted from $N$-body simulations based on the relative matter power spectrum $δ(k,a)$, defined as the fractional deviation in the absolute matter power spectrum produced by a target cosmology away from a reference $Λ$CDM prediction. From the computational perspective, $δ(k,a)$ is f… ▽ More We propose an alternative approach to the construction of fitting functions to the nonlinear matter power spectrum extracted from $N$-body simulations based on the relative matter power spectrum $δ(k,a)$, defined as the fractional deviation in the absolute matter power spectrum produced by a target cosmology away from a reference $Λ$CDM prediction. From the computational perspective, $δ(k,a)$ is fairly insensitive to the specifics of the simulation settings, and numerical convergence at the 1%-level can be readily achieved without the need for huge computing capacity. Furthermore, $δ(k,a)$ exhibits several interesting properties that enable a piece-wise construction of the full fitting function, whereby component fitting functions are sought for single-parameter variations and then multiplied together to form the final product. Then, to obtain 1%-accurate absolute power spectrum predictions for any target cosmology only requires that the community as a whole invests in producing one single ultra-precise reference $Λ$CDM absolute power spectrum, to be combined with the fitting function to produce the desired result. To illustrate the power of this approach, we have constructed the fitting function RelFit using only five relatively inexpensive $w$CDM simulations (box length $L=256 h^{-1}$Mpc, $N=1024^3$ particles, initialised at $z_i=49$). In a 6-parameter space spanning $\{ω_m,A_s,n_s,w,ω_b,h\}$, the output relative power spectra of RelFit are consistent with the predictions of the CosmicEmu emulator to 1% or better for a wide range of cosmologies up to $k\simeq 10$/Mpc. Thus, our approach could provide an inexpensive and democratically accessible route to fulfilling the 1%-level accuracy demands of the upcoming generation of large-scale structure probes, especially in the exploration of "non-standard" or "exotic" cosmologies on nonlinear scales. △ Less

Submitted 23 February, 2020; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: 37 pages, 16 figures; v2: 40 pages, two new sections on code comparison and relation to the stable clustering ansatz, two new figures, results and conclusions unchanged, version accepted by JCAP

Prior dependence of cosmological constraints on dark matter-radiation interactions

Authors: James A. D. Diacoumis, Yvonne. Y. Y. Wong

Abstract: We explore the issue of prior dependence in the context of one-sided constraints on the dark matter-photon and dark matter-neutrino elastic scattering cross-sections derived from cosmic microwave background (CMB) anisotropies measurements. Testing in particular the linear flat, Jeffreys, and logarithmic flat priors, we find that the former two yield upper limits on the cross-sections that are mutu… ▽ More We explore the issue of prior dependence in the context of one-sided constraints on the dark matter-photon and dark matter-neutrino elastic scattering cross-sections derived from cosmic microwave background (CMB) anisotropies measurements. Testing in particular the linear flat, Jeffreys, and logarithmic flat priors, we find that the former two yield upper limits on the cross-sections that are mutually consistent to within 20%. In contrast, bounds derived under the assumption of the logarithmic flat prior are strongly sensitive to the choice of the lower prior boundary. Indeed, surveying the recent literature, we find that this pathology of the logarithmic prior has resulted in published constraints that are up to an order of magnitude artificially tighter than they should objectively be. Our revised `objective' constraints from the 2015 data of the Planck CMB mission on the present-day scattering cross-sections are $σ_{\rm DM-γ} < 1.72 \times 10^{-6} \, σ_{\rm T} \, (m_{\rm DM}/{\rm GeV})$ and $σ_{\rm DM-γ} < 2.74 \times 10^{-15} \, σ_{\rm T} \, (m_{\rm DM}/{\rm GeV})$ for dark matter-photon interactions scaling as $a^0$ and $a^{-2}$ respectively, where $a$ is the scale factor, and $σ_{\rm T}$ the total Thomson scattering cross-section. Their dark matter-neutrino counterparts read $σ_{\rm DM-ν} < 2.14 \times 10^{-6} \, σ_{\rm T} \, (m_{\rm DM}/{\rm GeV})$ and $σ_{\rm DM-ν}< 2.46 \times 10^{-15} \, σ_{\rm T} \, (m_{\rm DM}/{\rm GeV})$. All have been computed assuming the Jeffreys prior. △ Less

Submitted 14 May, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Comments: 17 pages, 3 figures, 7 tables, comments welcome; v2: references added, small changes to the numerical results, version accepted by JCAP

Trading kinetic energy: How late kinetic decoupling of dark matter changes $N_{\textrm{eff}}$

Authors: James A. D. Diacoumis, Yvonne Y. Y. Wong

Abstract: Elastic scattering between dark matter particles and a relativistic species such as photons or neutrinos leads to a transfer of energy from the latter due to their intrinsically different temperature scaling relations. In this work, we point out that this siphoning of energy from the radiation bath manifests as a change in the effective number of neutrinos $N_{\rm eff}$, and compute the expected s… ▽ More Elastic scattering between dark matter particles and a relativistic species such as photons or neutrinos leads to a transfer of energy from the latter due to their intrinsically different temperature scaling relations. In this work, we point out that this siphoning of energy from the radiation bath manifests as a change in the effective number of neutrinos $N_{\rm eff}$, and compute the expected shift $ΔN_{\rm eff}$ for dark matter-photon and dark matter-neutrino elastic scattering as a function of the dark matter mass $m_ψ$ and scattering cross section $σ_{ψ-X}$. For $(m_ψ,σ_{ψ-X})$-parameter regions already explored by nonlinear probes such as the Lyman-$α$ forest through collisional and/or free-streaming damping, we find shifts of $|ΔN_{\rm eff}| \simeq O(10^{-2})$, which may be within the reach of the proposed CMB-S4 experiment. For most of the as-yet-unexplored parameter space, however, we expect $|ΔN_{\rm eff}| \lesssim O(10^{-3})$. An ideal 21 cm tomography survey of the dark ages limited only by cosmic variance is potentially sensitive to $|ΔN_{\rm eff}| \simeq O(10^{-6})$, in which case dark matter masses up to $m_ψ \sim 10 \, \textrm{MeV}$ may be probed via their effect on $N_{\rm eff}$. △ Less

Submitted 6 January, 2019; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: V2: added references, minor edits, matches published version

Journal ref: JCAP01(2019)001

Using CMB spectral distortions to distinguish between dark matter solutions to the small-scale crisis

Authors: James A. D. Diacoumis, Yvonne Y. Y. Wong

Abstract: The dissipation of small-scale perturbations in the early universe produces a distortion in the blackbody spectrum of cosmic microwave background photons. In this work, we propose to use these distortions as a probe of the microphysics of dark matter on scales $1\,\textrm{Mpc}^{-1}\lesssim k \lesssim 10^{4}\,\textrm{Mpc}^{-1}$. We consider in particular models in which the dark matter is kinetical… ▽ More The dissipation of small-scale perturbations in the early universe produces a distortion in the blackbody spectrum of cosmic microwave background photons. In this work, we propose to use these distortions as a probe of the microphysics of dark matter on scales $1\,\textrm{Mpc}^{-1}\lesssim k \lesssim 10^{4}\,\textrm{Mpc}^{-1}$. We consider in particular models in which the dark matter is kinetically coupled to either neutrinos or photons until shortly before recombination, and compute the photon heating rate and the resultant $μ$-distortion in both cases. We show that the $μ$-parameter is generally enhanced relative to $Λ$CDM for interactions with neutrinos, and may be either enhanced or suppressed in the case of interactions with photons. The deviations from the $Λ$CDM signal are potentially within the sensitivity reach of a PRISM-like experiment if $σ_{\textrm{DM}-γ} \gtrsim 1.1\times10^{-30} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $σ_{\textrm{DM}-ν} \gtrsim 4.8\times 10^{-32} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for time-independent cross sections, and $σ^{0}_{\textrm{DM}-γ} \gtrsim 1.8 \times 10^{-40} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $σ^{0}_{\textrm{DM}-ν} \gtrsim 2.5 \times 10^{-47} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for cross sections scaling as temperature squared, coinciding with the parameter regions in which late kinetic decoupling may serve as a solution to the small-scale crisis. Furthermore, these $μ$-distortion signals differ from those of warm dark matter (no deviation from $Λ$CDM) and a suppressed primordial power spectrum (strongly suppressed or a negative $μ$-parameter), demonstrating that CMB spectral distortion can potentially be used to distinguish between solutions to the small-scale crisis. △ Less

Submitted 11 September, 2017; v1 submitted 21 July, 2017; originally announced July 2017.

Comments: V2: added references, minor edits, matches published version

Journal ref: JCAP 1709 (2017) no.09, 011

Interacting neutrinos in cosmology: exact description and constraints

Authors: Isabel M. Oldengott, Thomas Tram, Cornelius Rampf, Yvonne Y. Y. Wong

Abstract: We consider the impact of neutrino self-interactions described by an effective four-fermion coupling on cosmological observations. Implementing the exact Boltzmann hierarchy for interacting neutrinos first derived in [arxiv:1409.1577] into the Boltzmann solver CLASS, we perform a detailed numerical analysis of the effects of the interaction on the cosmic microwave background (CMB) anisotropies, an… ▽ More We consider the impact of neutrino self-interactions described by an effective four-fermion coupling on cosmological observations. Implementing the exact Boltzmann hierarchy for interacting neutrinos first derived in [arxiv:1409.1577] into the Boltzmann solver CLASS, we perform a detailed numerical analysis of the effects of the interaction on the cosmic microwave background (CMB) anisotropies, and compare our results with known approximations in the literature. While we find good agreement between our exact approach and the relaxation time approximation used in some recent studies, the popular $\left( c_{\text{eff}}^2,c_{\text{vis}}^2 \right)$-parameterisation fails to reproduce the correct scale dependence of the CMB temperature power spectrum. We then proceed to derive constraints on the effective coupling constant $G_{\text{eff}}$ using currently available cosmological data via an MCMC analysis. Interestingly, our results reveal a bimodal posterior distribution, where one mode represents the standard $Λ$CDM limit with $G_{\rm eff} \lesssim 10^8 \, G_{\rm F}$, and the other a scenario in which neutrinos self-interact with an effective coupling constant $G_{\rm eff} \simeq 3 \times 10^9 \, G_{\rm F}$. △ Less

Submitted 28 November, 2017; v1 submitted 7 June, 2017; originally announced June 2017.

Comments: 25 pages, 11 figures; accepted for publication in JCAP

Journal ref: JCAP 1711 (2017) no.11, 027

Active-sterile neutrino oscillations in the early Universe with full collision terms

Authors: Steen Hannestad, Rasmus Sloth Hansen, Thomas Tram, Yvonne Y. Y. Wong

Abstract: Sterile neutrinos are thermalised in the early Universe via oscillations with the active neutrinos for certain mixing parameters. The most detailed calculation of this thermalisation process involves the solution of the momentum-dependent quantum kinetic equations, which track the evolution of the neutrino phase space distributions. Until now the collision terms in the quantum kinetic equations ha… ▽ More Sterile neutrinos are thermalised in the early Universe via oscillations with the active neutrinos for certain mixing parameters. The most detailed calculation of this thermalisation process involves the solution of the momentum-dependent quantum kinetic equations, which track the evolution of the neutrino phase space distributions. Until now the collision terms in the quantum kinetic equations have always been approximated using equilibrium distributions, but this approximation has never been checked numerically. In this work we revisit the sterile neutrino thermalisation calculation using the full collision term, and compare the results with various existing approximations in the literature. We find a better agreement than would naively be expected, but also identify some issues with these approximations that have not been appreciated previously. These include an unphysical production of neutrinos via scattering and the importance of redistributing momentum through scattering, as well as details of Pauli blocking. Finally, we devise a new approximation scheme, which improves upon some of the shortcomings of previous schemes. △ Less

Submitted 17 June, 2015; originally announced June 2015.

Comments: 30 pages, 7 figures

Report number: INT-PUB-15-025

How CMB and large-scale structure constrain chameleon interacting dark energy

Authors: Daniel Boriero, Subinoy Das, Yvonne Y. Y. Wong

Abstract: We explore a chameleon type of interacting dark matter-dark energy scenario in which a scalar field adiabatically traces the minimum of an effective potential sourced by the dark matter density. We discuss extensively the effect of this coupling on cosmological observables, especially the parameter degeneracies expected to arise between the model parameters and other cosmological parameters, and t… ▽ More We explore a chameleon type of interacting dark matter-dark energy scenario in which a scalar field adiabatically traces the minimum of an effective potential sourced by the dark matter density. We discuss extensively the effect of this coupling on cosmological observables, especially the parameter degeneracies expected to arise between the model parameters and other cosmological parameters, and then test the model against observations of the cosmic microwave background (CMB) anisotropies and other cosmological probes. We find that the chameleon parameters $α$ and $β$, which determine respectively the slope of the scalar field potential and the dark matter-dark energy coupling strength, can be constrained to $α< 0.17$ and $β< 0.19$ using CMB data alone. The latter parameter in particular is constrained only by the late Integrated Sachs-Wolfe effect. Adding measurements of the local Hubble expansion rate $H_0$ tightens the bound on $α$ by a factor of two, although this apparent improvement is arguably an artefact of the tension between the local measurement and the $H_0$ value inferred from Planck data in the minimal $Λ$CDM model. The same argument also precludes chameleon models from mimicking a dark radiation component, despite a passing similarity between the two scenarios in that they both delay the epoch of matter-radiation equality. Based on the derived parameter constraints, we discuss possible signatures of the model for ongoing and future large-scale structure surveys. △ Less

Submitted 12 May, 2015; originally announced May 2015.

Comments: 25 pages, 6 figures

Report number: BI-TP 2015/13

Future cosmological sensitivity for hot dark matter axions

Authors: Maria Archidiacono, Tobias Basse, Jan Hamann, Steen Hannestad, Georg Raffelt, Yvonne Y. Y. Wong

Abstract: We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch,… ▽ More We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch, assumed here to occur at a temperature $T_{\rm QCD} \sim 170$ MeV; this leaves an axion population of such low density that its late-time cosmological impact is negligible. For larger axion masses, $m_a \gtrsim 0.15$ eV, where axions remain in equilibrium until after the QCD phase transition, we find that a Euclid-like survey combined with Planck CMB data can detect $m_a$ at very high significance. Our conclusions are robust against assumptions about prior knowledge of the neutrino mass. Given that the proposed IAXO solar axion search is sensitive to $m_a\lesssim 0.2$ eV, the axion mass range probed by cosmology is nicely complementary. △ Less

Submitted 11 February, 2015; originally announced February 2015.

Comments: 17 pages, 5 figures

Higher-order massive neutrino perturbations in large-scale structure

Authors: Florian Führer, Yvonne Y. Y. Wong

Abstract: We develop a higher-order perturbation theory for large-scale structure formation involving a free-streaming hot or warm dark matter species. We focus on the case of mixed cold dark matter and massive neutrinos, although our approach is applicable also to a single warm dark matter species. In order to capture the suppressed growth of neutrino density perturbations on small scales, we account for… ▽ More We develop a higher-order perturbation theory for large-scale structure formation involving a free-streaming hot or warm dark matter species. We focus on the case of mixed cold dark matter and massive neutrinos, although our approach is applicable also to a single warm dark matter species. In order to capture the suppressed growth of neutrino density perturbations on small scales, we account for the full momentum dependence of the phase space distribution using the Vlaslov equation, and derive from it a formal closed-form nonlinear equation for the neutrino density. Using a systematic perturbative expansion of this equation we compute high-order corrections to the neutrino density contrast without the explicit need to track the perturbed neutrino momentum distribution. We calculate the leading-order total matter bispectrum for several neutrino masses. Using our result as a benchmark, we test the accuracy of the fluid approximation and a linear approximation used in perturbative and N-body analyses, as well as a new hybrid approach that combines the exact linear evolution with the nonlinear structure of the fluid equations. Aiming at $\lesssim1%$ accuracy, we find that the total matter bispectrum with a low neutrino mass m=0.046 eV can be reproduced by all but the fluid approximation, while for larger neutrino masses m=0.46-0.93 eV only the hybrid approach has the desired accuracy on a large range of scales. This result serves as a cautionary note that approximate nonlinear models of neutrino clustering that reproduce the gross features of some observables may not suffice for precision calculations, nor are they guaranteed to apply to other observables. All of the approximation schemes fail to reproduce the bispectrum of the neutrino density perturbations at better than 20% accuracy across all scales, indicating that an exact treatment of nonlinear neutrino perturbations is necessary. △ Less

Submitted 3 March, 2015; v1 submitted 8 December, 2014; originally announced December 2014.

Comments: v1: 42 pages, 6 figures v2: 44 pages, 8 figures, minor modifications, version accepted by JCAP

Getting leverage on inflation with a large photometric redshift survey

Authors: Tobias Basse, Jan Hamann, Steen Hannestad, Yvonne Y. Y. Wong

Abstract: We assess the potential of a future large-volume photometric redshift survey to constrain observational inflationary parameters using three large-scale structure observables: the angular shear and galaxy power spectra, and the cluster mass function measured through weak lensing. When used in combination with Planck-like CMB measurements, we find that the spectral index n_s can be constrained to a… ▽ More We assess the potential of a future large-volume photometric redshift survey to constrain observational inflationary parameters using three large-scale structure observables: the angular shear and galaxy power spectra, and the cluster mass function measured through weak lensing. When used in combination with Planck-like CMB measurements, we find that the spectral index n_s can be constrained to a 1 sigma precision of up to 0.0025. The sensitivity to the running of the spectral index can potentially improve to 0.0017, roughly a factor of five better than the present 1 sigma~constraint from Planck and auxiliary CMB data, allowing us to test the assumptions of the slow-roll scenario with unprecedented accuracy. Interestingly, neither CMB+shear nor CMB+galaxy nor CMB+clusters alone can achieve this level of sensitivity; it is the combined power of all three probes that conspires to break the different parameter degeneracies inherent in each type of observations. We make our forecast software publicly available via download or upon request from the authors. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Comments: 22 pages, 6 figures; the forecast software can be downloaded from http://jhamann.web.cern.ch/jhamann/simdata/simdata.tar.gz

Boltzmann hierarchy for interacting neutrinos I: formalism

Authors: Isabel M. Oldengott, Cornelius Rampf, Yvonne Y. Y. Wong

Abstract: Starting from the collisional Boltzmann equation, we derive for the first time and from first principles the Boltzmann hierarchy for neutrinos including interactions with a scalar particle. Such interactions appear, for example, in majoron-like models of neutrino mass generation. We study two limits of the scalar mass: (i) An extremely massive scalar whose only role is to mediate an effective 4-fe… ▽ More Starting from the collisional Boltzmann equation, we derive for the first time and from first principles the Boltzmann hierarchy for neutrinos including interactions with a scalar particle. Such interactions appear, for example, in majoron-like models of neutrino mass generation. We study two limits of the scalar mass: (i) An extremely massive scalar whose only role is to mediate an effective 4-fermion neutrino-neutrino interaction, and (ii) a massless scalar that can be produced in abundance and thus demands its own Boltzmann hierarchy. In contrast to, e.g., the first-order Boltzmann hierarchy for Thomson-scattering photons, our interacting neutrino/scalar Boltzmann hierarchies contain additional momentum-dependent collision terms arising from a non-negligible energy transfer in the neutrino-neutrino and neutrino-scalar interactions. This necessitates that we track each momentum mode of the phase space distributions individually, even if the particles were massless. Comparing our hierarchy with the commonly used $(c_{\rm eff}^2,c_{\rm vis}^2)$-parameterisation, we find no formal correspondence between the two approaches, which raises the question of whether the latter parameterisation even has an interpretation in terms of particle scattering. Lastly, although we have invoked majoron-like models as a motivation for our study, our treatment is in fact generally applicable to all scenarios in which the neutrino and/or other ultrarelativistic fermions interact with scalar particles. △ Less

Submitted 16 March, 2015; v1 submitted 4 September, 2014; originally announced September 2014.

Comments: 44 pages, 14 figures; included scalar Boltzmann hierarchy in the massless case and plots of integral kernels; accepted by JCAP

How real-time cosmology can distinguish between different anisotropic models

Authors: Luca Amendola, Ole Eggers Bjaelde, Wessel Valkenburg, Yvonne Y. Y. Wong

Abstract: We present a new analysis on how to distinguish between isotropic and anisotropic cosmological models based on tracking the angular displacements of a large number of distant quasars over an extended period of time, and then performing a multipole-vector decomposition of the resulting displacement maps. We find that while the GAIA mission operating at its nominal specifications does not have suffi… ▽ More We present a new analysis on how to distinguish between isotropic and anisotropic cosmological models based on tracking the angular displacements of a large number of distant quasars over an extended period of time, and then performing a multipole-vector decomposition of the resulting displacement maps. We find that while the GAIA mission operating at its nominal specifications does not have sufficient angular resolution to resolve anisotropic universes from isotropic ones using this method within a reasonable timespan of ten years, a next-generation GAIA-like survey with a resolution ten times better should be equal to the task. Distinguishing between different anisotropic models is however more demanding. Keeping the observational timespan to ten years, we find that the angular resolution of the survey will need to be of order 0.1 micro-arcsecond in order for certain rotating anisotropic models to produce a detectable signature that is also unique to models of this class. However, should such a detection become possible, it would immediately allow us to rule out large local void models. △ Less

Submitted 15 December, 2013; v1 submitted 9 July, 2013; originally announced July 2013.

Comments: 25 pages, 4 figures; v2: 27 pages, 5 figures in new colour scheme, extended discussion on the uniqueness of the observer position, main conclusions unchanged, matches version accepted by JCAP

Axion hot dark matter bounds after Planck

Authors: Maria Archidiacono, Steen Hannestad, Alessandro Mirizzi, Georg Raffelt, Yvonne Y. Y. Wong

Abstract: We use cosmological observations in the post-Planck era to derive limits on thermally produced cosmological axions. In the early universe such axions contribute to the radiation density and later to the hot dark matter fraction. We find an upper limit m_a < 0.67 eV at 95% C.L. after marginalising over the unknown neutrino masses, using CMB temperature and polarisation data from Planck and WMAP res… ▽ More We use cosmological observations in the post-Planck era to derive limits on thermally produced cosmological axions. In the early universe such axions contribute to the radiation density and later to the hot dark matter fraction. We find an upper limit m_a < 0.67 eV at 95% C.L. after marginalising over the unknown neutrino masses, using CMB temperature and polarisation data from Planck and WMAP respectively, the halo matter power spectrum extracted from SDSS-DR7, and the local Hubble expansion rate H_0 released by the Carnegie Hubble Program based on a recalibration of the Hubble Space Telescope Key Project sample. Leaving out the local H_0 measurement relaxes the limit somewhat to 0.86 eV, while Planck+WMAP alone constrain the axion mass to 1.01 eV, the first time an upper limit on m_a has been obtained from CMB data alone. Our axion limit is therefore not very sensitive to the tension between the Planck-inferred H_0 and the locally measured value. This is in contrast with the upper limit on the neutrino mass sum, which we find here to range from 0.27 eV at 95% C.L. combining all of the aforementioned observations, to 0.84 eV from CMB data alone. △ Less

Submitted 21 October, 2013; v1 submitted 2 July, 2013; originally announced July 2013.

Comments: 20 pages, 8 figures, matches version published in JCAP 1310 (2013) 020

Report number: MPP-2013-113

Journal ref: JCAP 1310 (2013) 020

Dark energy properties from large future galaxy surveys

Authors: Tobias Basse, Ole Eggers Bjaelde, Jan Hamann, Steen Hannestad, Yvonne Y. Y. Wong

Abstract: We perform a detailed forecast on how well a {\sc Euclid}-like survey will be able to constrain dark energy and neutrino parameters from a combination of its cosmic shear power spectrum, galaxy power spectrum, and cluster mass function measurements. We find that the combination of these three probes vastly improves the survey's potential to measure the time evolution of dark energy. In terms of a… ▽ More We perform a detailed forecast on how well a {\sc Euclid}-like survey will be able to constrain dark energy and neutrino parameters from a combination of its cosmic shear power spectrum, galaxy power spectrum, and cluster mass function measurements. We find that the combination of these three probes vastly improves the survey's potential to measure the time evolution of dark energy. In terms of a dark energy figure-of-merit defined as $(σ(w_{\mathrm p}) σ(w_a))^{-1}$, we find a value of 690 for {\sc Euclid}-like data combined with {\sc Planck}-like measurements of the cosmic microwave background (CMB) anisotropies in a 10-dimensional cosmological parameter space, assuming a $Λ$CDM fiducial cosmology. For the more commonly used 7-parameter model, we find a figure-of-merit of 1900 for the same data combination. We consider also the survey's potential to measure dark energy perturbations in models wherein the dark energy is parameterised as a fluid with a nonstandard non-adiabatic sound speed, and find that in an \emph{optimistic} scenario in which $w_0$ deviates by as much as is currently observationally allowed from $-1$, models with $\hat{c}_\mathrm{s}^2 = 10^{-6}$ and $\hat{c}_\mathrm{s}^2 = 1$ can be distinguished at more than $2σ$ significance. We emphasise that constraints on the dark energy sound speed from cluster measurements are strongly dependent on the modelling of the cluster mass function; significantly weaker sensitivities ensue if we modify our model to include fewer features of nonlinear dark energy clustering. Finally, we find that the sum of neutrino masses can be measured with a $1 σ$ precision of 0.015~eV, (abridged) △ Less

Submitted 29 April, 2014; v1 submitted 8 April, 2013; originally announced April 2013.

Comments: 26 pages, 5 figures, matches JCAP version

Measuring neutrino masses with a future galaxy survey

Authors: Jan Hamann, Steen Hannestad, Yvonne Y. Y. Wong

Abstract: We perform a detailed forecast on how well a Euclid-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of sum m_nu ~ 0.06 eV in the normal hierarchy can be detected at 1.5 sigma to 2.5 si… ▽ More We perform a detailed forecast on how well a Euclid-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of sum m_nu ~ 0.06 eV in the normal hierarchy can be detected at 1.5 sigma to 2.5 sigma significance, depending on the model complexity, using a combination of galaxy and cosmic shear power spectrum measurements in conjunction with CMB temperature and polarisation observations from Planck. With better knowledge of the galaxy bias, the significance of the detection could potentially reach 5.4 sigma. Interestingly, neither Planck+shear nor Planck+galaxy alone can achieve this level of sensitivity; it is the combined effect of galaxy and cosmic shear power spectrum measurements that breaks the persistent degeneracies between the neutrino mass, the physical matter density, and the Hubble parameter. Notwithstanding this remarkable sensitivity to sum m_nu, Euclid-like shear and galaxy data will not be sensitive to the exact mass spectrum of the neutrino sector; no significant bias (< 1 sigma) in the parameter estimation is induced by fitting inaccurate models of the neutrino mass splittings to the mock data, nor does the goodness-of-fit of these models suffer any significant degradation relative to the true one (Delta chi_eff ^2< 1). △ Less

Submitted 7 November, 2012; v1 submitted 5 September, 2012; originally announced September 2012.

Comments: v1: 29 pages, 10 figures. v2: 33 pages, 12 figures; added sections on shape evolution and constraints in more complex models, accepted for publication in JCAP

Confronting the sound speed of dark energy with future cluster surveys

Authors: Tobias Basse, Ole Eggers Bjaelde, Steen Hannestad, Yvonne Y. Y. Wong

Abstract: Future cluster surveys will observe galaxy clusters numbering in the hundred thousands. We consider this work how these surveys can be used to constrain dark energy parameters: in particular, the equation of state parameter w and the non-adiabatic sound speed c_s^2. We demonstrate that, in combination with Cosmic Microwave Background (CMB) observations from Planck, cluster surveys such as that in… ▽ More Future cluster surveys will observe galaxy clusters numbering in the hundred thousands. We consider this work how these surveys can be used to constrain dark energy parameters: in particular, the equation of state parameter w and the non-adiabatic sound speed c_s^2. We demonstrate that, in combination with Cosmic Microwave Background (CMB) observations from Planck, cluster surveys such as that in the ESA Euclid project will be able to determine a time-independent w with subpercent precision. Likewise, if the dark energy sound horizon falls within the length scales probed by the cluster survey, then c_s^2 can be pinned down to within an order of magnitude. In the course of this work, we also investigate the process of dark energy virialisation in the presence of an arbitrary sound speed. We find that dark energy clustering and virialisation can lead to dark energy contributing to the total cluster mass at approximately the 0.1% level at maximum. △ Less

Submitted 2 May, 2012; originally announced May 2012.

Comments: 24 pages, 7 figures, prepared for submission to JCAP

Light Sterile Neutrinos: A White Paper

Authors: K. N. Abazajian, M. A. Acero, S. K. Agarwalla, A. A. Aguilar-Arevalo, C. H. Albright, S. Antusch, C. A. Arguelles, A. B. Balantekin, G. Barenboim, V. Barger, P. Bernardini, F. Bezrukov, O. E. Bjaelde, S. A. Bogacz, N. S. Bowden, A. Boyarsky, A. Bravar, D. Bravo Berguno, S. J. Brice, A. D. Bross, B. Caccianiga, F. Cavanna, E. J. Chun, B. T. Cleveland, A. P. Collin , et al. (162 additional authors not shown)

Abstract: This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data. This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data. △ Less

Submitted 18 April, 2012; originally announced April 2012.

Lagrangian perturbations and the matter bispectrum II: the resummed one-loop correction to the matter bispectrum

Authors: Cornelius Rampf, Yvonne Y. Y. Wong

Abstract: This is part two in a series of papers in which we investigate an approach based on Lagrangian perturbation theory (LPT) to study the non-linear evolution of the large-scale structure distribution in the universe. Firstly, we compute the matter bispectrum in real space using LPT up one-loop order, for both Gaussian and non-Gaussian initial conditions. In the initial position limit, we find that th… ▽ More This is part two in a series of papers in which we investigate an approach based on Lagrangian perturbation theory (LPT) to study the non-linear evolution of the large-scale structure distribution in the universe. Firstly, we compute the matter bispectrum in real space using LPT up one-loop order, for both Gaussian and non-Gaussian initial conditions. In the initial position limit, we find that the one-loop bispectrum computed in this manner is identical to its counterpart obtained from standard Eulerian perturbation theory (SPT). Furthermore, the LPT formalism allows for a simple reorganisation of the perturbative series corresponding to the resummation of an infinite series of perturbations in SPT. Applying this method, we find a resummed one-loop bispectrum that compares favourably with results from N-body simulations. We generalise the resummation method also to the computation of the redshift-space bispectrum up to one loop. △ Less

Submitted 14 June, 2012; v1 submitted 19 March, 2012; originally announced March 2012.

Comments: 33 pages, 14 of them appendix, 4 figures, matches published version in JCAP

Report number: TTK-12-10

Journal ref: JCAP06(2012)018

Evidence for dark matter modulation in CoGeNT?

Authors: Chiara Arina, Jan Hamann, Roberto Trotta, Yvonne Y Y Wong

Abstract: We investigate the question of whether the recent modulation signal claimed by CoGeNT is best explained by the dark matter (DM) hypothesis from a Bayesian model comparison perspective. We consider five phenomenological explanations for the data: no modulation signal, modulation due to DM, modulation due to DM compatible with the total CoGeNT rate, and a signal coming from other physics with a free… ▽ More We investigate the question of whether the recent modulation signal claimed by CoGeNT is best explained by the dark matter (DM) hypothesis from a Bayesian model comparison perspective. We consider five phenomenological explanations for the data: no modulation signal, modulation due to DM, modulation due to DM compatible with the total CoGeNT rate, and a signal coming from other physics with a free phase but annual period, or with a free phase and a free period. In each scenario, we assign to the free parameters physically motivated priors. We find that while the DM models are weakly preferred to the no modulation model, but when compared to models where the modulation is due to other physics, the DM hypothesis is favoured with odds ranging from 185:1 to 560:1. This result is robust even when astrophysical uncertainties are taken into account and the impact of priors assessed. Interestingly, the odds for the DM model in which the modulation signal is compatible with the total rate against a DM model in which this prior is not implemented is only 5:8, in spite of the former's prediction of a modulation amplitude in the energy range 0.9 to 3.0 keVee that is significantly smaller than the value observed by CoGeNT. Classical hypothesis testing also rules out the null hypothesis of no modulation at the 1.6 sigma to 2.3 sigma level, depending on the details of the alternative. Lastly, we investigate whether anisotropic velocity distributions can help to mitigate the tension between the CoGeNT total and modulated rates, and find encouraging results. △ Less

Submitted 30 January, 2012; v1 submitted 14 November, 2011; originally announced November 2011.

Comments: 27 pages, 8 figures and 7 tables; corrected numerical error, match the submitted version

Neutrino mass in cosmology: status and prospects

Authors: Yvonne Y. Y. Wong

Abstract: I give an overview of the effects of neutrino masses in cosmology, focussing on the role they play in the evolution of cosmological perturbations. I discuss how recent observations of the cosmic microwave background anisotropies and the large-scale matter distribution can probe neutrino masses with greater precision than current laboratory experiments. I describe several new techniques that will b… ▽ More I give an overview of the effects of neutrino masses in cosmology, focussing on the role they play in the evolution of cosmological perturbations. I discuss how recent observations of the cosmic microwave background anisotropies and the large-scale matter distribution can probe neutrino masses with greater precision than current laboratory experiments. I describe several new techniques that will be used to probe cosmology in the future, as well as recent advances in the computation of the nonlinear matter power spectrum and related observables. △ Less

Submitted 6 November, 2011; originally announced November 2011.

Comments: 41 pages, 5 figures, written for the Annual Review of Nuclear and Particle Science in January 2011. Publisher limits the number of references. Apologies if I missed yours

Report number: TTK-11-56

Journal ref: Ann.Rev.Nucl.Part.Sci.(2011)61:69-98

Sterile neutrinos with eV masses in cosmology -- how disfavoured exactly?

Authors: Jan Hamann, Steen Hannestad, Georg G. Raffelt, Yvonne Y. Y. Wong

Abstract: We study cosmological models that contain sterile neutrinos with eV-range masses as suggested by reactor and short-baseline oscillation data. We confront these models with both precision cosmological data (probing the CMB decoupling epoch) and light-element abundances (probing the BBN epoch). In the minimal LambdaCDM model, such sterile neutrinos are strongly disfavoured by current data because th… ▽ More We study cosmological models that contain sterile neutrinos with eV-range masses as suggested by reactor and short-baseline oscillation data. We confront these models with both precision cosmological data (probing the CMB decoupling epoch) and light-element abundances (probing the BBN epoch). In the minimal LambdaCDM model, such sterile neutrinos are strongly disfavoured by current data because they contribute too much hot dark matter. However, if the cosmological framework is extended to include also additional relativistic degrees of freedom -- beyond the three standard neutrinos and the putative sterile neutrinos, then the hot dark matter constraint on the sterile states is considerably relaxed. A further improvement is achieved by allowing a dark energy equation of state parameter w<-1. While BBN strongly disfavours extra radiation beyond the assumed eV-mass sterile neutrino, this constraint can be circumvented by a small nu_e degeneracy. Any model containing eV-mass sterile neutrinos implies also strong modifications of other cosmological parameters. Notably, the inferred cold dark matter density can shift up by 20 to 75% relative to the standard LambdaCDM value. △ Less

Submitted 28 September, 2011; v1 submitted 20 August, 2011; originally announced August 2011.

Comments: 14 pages, 6 figures, v2: minor changes, matches version accepted for publication in JCAP

Report number: MPP-2011-100, TTK-11-29

Journal ref: JCAP 1109 (2011) 034

A Bayesian view of the current status of dark matter direct searches

Authors: Chiara Arina, Jan Hamann, Yvonne Y. Y. Wong

Abstract: Bayesian statistical methods offer a simple and consistent framework for incorporating uncertainties into a multi-parameter inference problem. In this work we apply these methods to a selection of current direct dark matter searches. We consider the simplest scenario of spin-independent elastic WIMP scattering, and infer the WIMP mass and cross-section from the experimental data with the essential… ▽ More Bayesian statistical methods offer a simple and consistent framework for incorporating uncertainties into a multi-parameter inference problem. In this work we apply these methods to a selection of current direct dark matter searches. We consider the simplest scenario of spin-independent elastic WIMP scattering, and infer the WIMP mass and cross-section from the experimental data with the essential systematic uncertainties folded into the analysis. We find that when uncertainties in the scintillation efficiency of Xenon100 have been accounted for, the resulting exclusion limit is not sufficiently constraining to rule out the CoGeNT preferred parameter region, contrary to previous claims. In the same vein, we also investigate the impact of astrophysical uncertainties on the preferred WIMP parameters. We find that within the class of smooth and isotropic WIMP velocity distributions, it is difficult to reconcile the DAMA and the CoGeNT preferred regions by tweaking the astrophysics parameters alone. If we demand compatibility between these experiments, then the inference process naturally concludes that a high value for the sodium quenching factor for DAMA is preferred. △ Less

Submitted 13 September, 2011; v1 submitted 25 May, 2011; originally announced May 2011.

Comments: 37 pages, 14 figures and 7 tables. Replacement for matching the version accepted for publication

Cosmological and Astrophysical Neutrino Mass Measurements

Authors: K. N. Abazajian, E. Calabrese, A. Cooray, F. De Bernardis, S. Dodelson, A. Friedland, G. M. Fuller, S. Hannestad, B. G. Keating, E. V. Linder, C. Lunardini, A. Melchiorri, R. Miquel, E. Pierpaoli, J. Pritchard, P. Serra, M. Takada, Y. Y. Y. Wong

Abstract: Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach. Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach. △ Less

Submitted 25 March, 2011; originally announced March 2011.

Comments: 11 pages

Journal ref: Astropart. Phys. 35, 177 (2011)

Spherical collapse of dark energy with an arbitrary sound speed

Authors: Tobias Basse, Ole Eggers Bjaelde, Yvonne Y. Y. Wong

Abstract: We consider a generic type of dark energy fluid, characterised by a constant equation of state parameter w and sound speed c_s, and investigate the impact of dark energy clustering on cosmic structure formation using the spherical collapse model. Along the way, we also discuss in detail the evolution of dark energy perturbations in the linear regime. We find that the introduction of a finite sound… ▽ More We consider a generic type of dark energy fluid, characterised by a constant equation of state parameter w and sound speed c_s, and investigate the impact of dark energy clustering on cosmic structure formation using the spherical collapse model. Along the way, we also discuss in detail the evolution of dark energy perturbations in the linear regime. We find that the introduction of a finite sound speed into the picture necessarily induces a scale-dependence in the dark energy clustering, which in turn affects the dynamics of the spherical collapse in a scale-dependent way. As with other, more conventional fluids, we can define a Jeans scale for the dark energy clustering, and hence a Jeans mass M_J for the dark matter which feels the effect of dark energy clustering via gravitational interactions. For bound objects (halos) with masses M >> M_J, the effect of dark energy clustering is maximal. For those with M << M_J, the dark energy component is effectively homogeneous, and its role in the formation of these structures is reduced to its effects on the Hubble expansion rate. To compute quantitatively the virial density and the linearly extrapolated threshold density, we use a quasi-linear approach which is expected to be valid up to around the Jeans mass. We find an interesting dependence of these quantities on the halo mass M, given some w and c_s. The dependence is the strongest for masses lying in the vicinity of M ~ M_J. Observing this M-dependence will be a tell-tale sign that dark energy is dynamic, and a great leap towards pinning down its clustering properties. △ Less

Submitted 27 April, 2012; v1 submitted 31 August, 2010; originally announced September 2010.

Comments: 25 pages, 6 figures, matches version published in JCAP

Cosmology seeking friendship with sterile neutrinos

Authors: Jan Hamann, Steen Hannestad, Georg G. Raffelt, Irene Tamborra, Yvonne Y. Y. Wong

Abstract: Precision cosmology and big-bang nucleosynthesis mildly favor extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to derive credible regions f… ▽ More Precision cosmology and big-bang nucleosynthesis mildly favor extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to derive credible regions for the assumed common mass scale m_s and effective number N_s of thermally excited sterile neutrino states. Our results are compatible with the existence of one or perhaps two sterile neutrinos, as suggested by LSND and MiniBooNE, if m_s is in the sub-eV range. △ Less

Submitted 26 October, 2010; v1 submitted 28 June, 2010; originally announced June 2010.

Comments: 4 pages, 1 figure, matches version published in PRL

Report number: MPP-2010-73, TTK-10-36

Journal ref: Phys.Rev.Lett.105:181301,2010

Neutrinos in Non-linear Structure Formation - The Effect on Halo Properties

Authors: Jacob Brandbyge, Steen Hannestad, Troels Haugboelle, Yvonne Y. Y. Wong

Abstract: We use N-body simulations to find the effect of neutrino masses on halo properties, and investigate how the density profiles of both the neutrino and the dark matter components change as a function of the neutrino mass. We compare our neutrino density profiles with results from the N-one-body method and find good agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos is… ▽ More We use N-body simulations to find the effect of neutrino masses on halo properties, and investigate how the density profiles of both the neutrino and the dark matter components change as a function of the neutrino mass. We compare our neutrino density profiles with results from the N-one-body method and find good agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos is not saturated. Finally we study how the halo mass function changes as a function of the neutrino mass and compare our results with the Sheth-Tormen semi-analytic formulae. Our results are important for surveys which aim at probing cosmological parameters using clusters, as well as future experiments aiming at measuring the cosmic neutrino background directly. △ Less

Submitted 23 April, 2010; originally announced April 2010.

Comments: 20 pages, 8 figures

Neutrino and axion hot dark matter bounds after WMAP-7

Authors: Steen Hannestad, Alessandro Mirizzi, Georg G. Raffelt, Yvonne Y. Y. Wong

Abstract: We update cosmological hot dark matter constraints on neutrinos and hadronic axions. Our most restrictive limits use 7-year data from the Wilkinson Microwave Anisotropy Probe for the cosmic microwave background anisotropies, the halo power spectrum (HPS) from the 7th data release of the Sloan Digital Sky Survey, and the Hubble constant from Hubble Space Telescope observations. We find 95% C.L. upp… ▽ More We update cosmological hot dark matter constraints on neutrinos and hadronic axions. Our most restrictive limits use 7-year data from the Wilkinson Microwave Anisotropy Probe for the cosmic microwave background anisotropies, the halo power spectrum (HPS) from the 7th data release of the Sloan Digital Sky Survey, and the Hubble constant from Hubble Space Telescope observations. We find 95% C.L. upper limits of \sum m_ν<0.44 eV (no axions), m_a<0.91 eV (assuming \sum m_ν=0), and \sum m_ν<0.41 eV and m_a<0.72 eV for two hot dark matter components after marginalising over the respective other mass. CMB data alone yield \sum m_ν<1.19 eV (no axions), while for axions the HPS is crucial for deriving m_a constraints. This difference can be traced to the fact that for a given hot dark matter fraction axions are much more massive than neutrinos. △ Less

Submitted 12 July, 2010; v1 submitted 5 April, 2010; originally announced April 2010.

Comments: 9 pages, 3 figures, uses iopart.cls; v2: one additional figure, references added, version accepted by JCAP

Report number: MPP-2010-36, TTK-10-28

Journal ref: JCAP 1008:001,2010

Cosmological parameters from large scale structure - geometric versus shape information

Authors: Jan Hamann, Steen Hannestad, Julien Lesgourgues, Cornelius Rampf, Yvonne Y. Y. Wong

Abstract: The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape info… ▽ More The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape information currently contained in the SDSS DR7 halo power spectrum (HPS) is by far superseded by geometric information derived from the baryonic features. An immediate corollary is that contrary to popular beliefs, the upper limit on the neutrino mass m_νpresently derived from LSS combined with cosmic microwave background (CMB) data does not in fact arise from the possible small-scale power suppression due to neutrino free-streaming, if we limit the model framework to minimal LambdaCDM+m_ν. However, in more complicated models, such as those extended with extra light degrees of freedom and a dark energy equation of state parameter w differing from -1, shape information becomes crucial for the resolution of parameter degeneracies. This conclusion will remain true even when data from the Planck surveyor become available. In the course of our analysis, we introduce a new dewiggling procedure that allows us to extend consistently the use of the SDSS HPS to models with an arbitrary sound horizon at decoupling. All the cases considered here are compatible with the conservative 95%-bounds \sum m_ν< 1.16 eV, N_eff = 4.8 \pm 2.0. △ Less

Submitted 15 July, 2010; v1 submitted 21 March, 2010; originally announced March 2010.

Comments: 18 pages, 4 figures; v2: references added, matches published version

Report number: CERN-PH-TH/2010-064, LAPTH-016/10, TTK-10-27

Cosmological axion bounds

Authors: Steen Hannestad, Jan Hamann, Alessandro Mirizzi, Georg G. Raffelt, Yvonne Y. Y. Wong

Abstract: We discuss current cosmological constraints on axions, as well as future sensitivities. Bounds on axion hot dark matter are discussed first, and subsequently we discuss both current and future sensitivity to models in which axions play the role as cold dark matter, but where the Peccei-Quinn symmetry is not restored during reheating. We discuss current cosmological constraints on axions, as well as future sensitivities. Bounds on axion hot dark matter are discussed first, and subsequently we discuss both current and future sensitivity to models in which axions play the role as cold dark matter, but where the Peccei-Quinn symmetry is not restored during reheating. △ Less

Submitted 29 October, 2009; originally announced October 2009.

Comments: 4 pages, 2 figures, To appear in the proceedings of 5th Patras Workshop on Axions, WIMPs and WISPs, Durham 13-17 July 2009

Isocurvature forecast in the anthropic axion window

Authors: J. Hamann, S. Hannestad, G. G. Raffelt, Y. Y. Y. Wong

Abstract: We explore the cosmological sensitivity to the amplitude of isocurvature fluctuations that would be caused by axions in the "anthropic window" where the axion decay constant f_a >> 10^12 GeV and the initial misalignment angle Theta_i << 1. In a minimal Lambda-CDM cosmology extended with subdominant scale-invariant isocurvature fluctuations, existing data constrain the isocurvature fraction to al… ▽ More We explore the cosmological sensitivity to the amplitude of isocurvature fluctuations that would be caused by axions in the "anthropic window" where the axion decay constant f_a >> 10^12 GeV and the initial misalignment angle Theta_i << 1. In a minimal Lambda-CDM cosmology extended with subdominant scale-invariant isocurvature fluctuations, existing data constrain the isocurvature fraction to alpha < 0.09 at 95% C.L. If no signal shows up, Planck can improve this constraint to 0.042 while an ultimate CMB probe limited only by cosmic variance in both temperature and E-polarisation can reach 0.017, about a factor of five better than the current limit. In the parameter space of f_a and H_I (Hubble parameter during inflation) we identify a small region where axion detection remains within the reach of realistic cosmological probes. △ Less

Submitted 16 June, 2009; v1 submitted 3 April, 2009; originally announced April 2009.

Comments: 14 pages, 4 figures; v2: matches published version

Report number: CERN-PH-TH/2009-044, LAPTH-1323/09, MPP-2009-37

Journal ref: JCAP 0906:022,2009

Higher order corrections to the large scale matter power spectrum in the presence of massive neutrinos

Authors: Yvonne Y. Y. Wong

Abstract: We present the first systematic derivation of the one-loop correction to the large scale matter power spectrum in a mixed cold+hot dark matter cosmology with subdominant massive neutrino hot dark matter. Starting with the equations of motion for the density and velocity fields, we derive perturbative solutions to these quantities and construct recursion relations for the interaction kernels, not… ▽ More We present the first systematic derivation of the one-loop correction to the large scale matter power spectrum in a mixed cold+hot dark matter cosmology with subdominant massive neutrino hot dark matter. Starting with the equations of motion for the density and velocity fields, we derive perturbative solutions to these quantities and construct recursion relations for the interaction kernels, noting and justifying all approximations along the way. We find interaction kernels similar to those for a cold dark matter-only universe, but with additional dependences on the neutrino energy density fraction f_nu and the linear growth functions of the incoming wavevectors. Compared with the f_nu=0 case, the one-loop corrected matter power spectrum for a mixed dark matter cosmology exhibits a decrease in small scale power exceeding the canonical ~8 f_nu suppression predicted by linear theory, a feature also seen in multi-component N-body simulations. △ Less

Submitted 6 October, 2008; v1 submitted 4 September, 2008; originally announced September 2008.

Comments: 22 pages, 5 figures, uses iopart.cls; v2: 23 pages, added references, to appear in JCAP

Report number: MPP-2008-111

Journal ref: JCAP 0810:035,2008