Search

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., 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

Odious sovereign contracts inhibit developing country growth. They cause money and resources to be improperly transferred from one country to an undeserving one. Think tanks and civil society actors have long suggested that Indonesia is plagued by a sizeable odious sovereign contract account. Policy makers and scholars continue to grapple with how the odious debts doctrine may operate in law to curb this important problem. Taking into account the lack of transparency in an odious contract setting, this paper proposes a new approach premised on principles of transparency, accountability and citizen participation in public contracts. In design, it proposes the following: 1) the creation of ex-ante obligations and a public website on which financiers can disclose the key terms of their contractual arrangements with a sovereign government. This website enables a financier to signal the nature of their engagement with a sovereign counterpart; 2) the creation of an ex-post tribunal, in which private citizens have standing, to adjudicate the odiousness of a disputed sovereign contract. In the tribunal’s deliberations, disclosure by a financier and compliance with ex-ante obligations weighs in favour of a presumption of legitimate contracting, whereas non-disclosure lends itself to a presumption of odiousness. This new approach has the right incentives for participation. It will revolutionize the currency of international law and international institutions, by giving the public a mechanism to eke out odiousness in transnational sovereign dealings. It can have important implications for Indonesia and more generally, the future of transnational trade and finance.

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_\chi =\frac{\sigma_0}{\sigma_{Th}}\left(\frac{m_\chi}{100 \text{GeV}}\right)^{-1}$, is $u_\chi\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 $\sigma_8$ that is inferred in the context of $\Lambda$CDM from the Planck data, leading to agreement within 1-2$\sigma$ with weak lensing estimates of $\sigma_8$, as those from KiDS-1000. However, the presence of these interactions barely affects the value of the Hubble constant $H_0$., Comment: 31 pages, 8 figures, 3 tables

We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, Neff 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 O(e3), 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 Neff 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 Neff SM 3.0440 ±0.0002, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure (|δ Neff| ∼10-4), augmented by measurement errors in the solar mixing angle sin2θ 12 (|δ Neff| ∼10-4).

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., 24 pages, 9 figures (matches accepted version)

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

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 $\delta 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 $\delta 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 $\delta 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$., Comment: 25 pages, plus appendices and 8 figures; v3: references added, version accepted by JCAP, minor terminology changes

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

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 $\nu_H$ couples to a lighter one $\nu_l$ and a massless scalar particle $\phi$ 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 $\nu_H$, $\nu_l$, and $\phi$ 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 $\nu_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 $\Gamma_{\rm T} =(1/\tau_0) (m_{\nu H}/E_{\nu H})^3$, but at a rate $\sim (1/\tau_0) (m_{\nu H}/E_{\nu 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 $\tau_0$ from $\tau_0^{\rm old} \gtrsim 1.2 \times 10^9\, {\rm s}\, (m_{\nu H}/50\, {\rm meV})^3$ to $\tau_0 \gtrsim (4 \times 10^5 \to 4 \times 10^6)\, {\rm s}\, (m_{\nu H}/50 \, {\rm meV})^5$., Comment: 53 pages, 10 figures, v2 matches the accepted JCAP version

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 $\Lambda$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}$., Comment: 25 pages, 11 figures; accepted for publication in JCAP

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 $\Delta N_{\rm eff}$ for dark matter-photon and dark matter-neutrino elastic scattering as a function of the dark matter mass $m_\psi$ and scattering cross section $\sigma_{\psi-X}$. For $(m_\psi,\sigma_{\psi-X})$-parameter regions already explored by nonlinear probes such as the Lyman-$\alpha$ forest through collisional and/or free-streaming damping, we find shifts of $|\Delta 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 $|\Delta 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 $|\Delta N_{\rm eff}| \simeq O(10^{-6})$, in which case dark matter masses up to $m_{\psi} \sim 10 \, \textrm{MeV}$ may be probed via their effect on $N_{\rm eff}$., Comment: V2: added references, minor edits, matches published version

I give an overview of the effects of neutrinos on cosmology, focussing in particular on the role played by massive neutrinos in the evolution of cosmological perturbations. I discuss how recent observations of the cosmic microwave background anisotropies and the large-scale structure of galaxies 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, and also some recent advances in the computation of the nonlinear matter power spectrum and related observables.

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., Comment: 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

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., 30 pages, 7 figures

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 $\alpha$ and $\beta$, which determine respectively the slope of the scalar field potential and the dark matter-dark energy coupling strength, can be constrained to $\alpha < 0.17$ and $\beta < 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 $\alpha$ 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 $\Lambda$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., Comment: 25 pages, 6 figures

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., Comment: 17 pages, 5 figures

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., 44 pages, 14 figures; included scalar Boltzmann hierarchy in the massless case and plots of integral kernels; accepted by JCAP

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., v1: 42 pages, 6 figures v2: 44 pages, 8 figures, minor modifications, version accepted by JCAP

If the small and large coupling behavior of a physical system can be computed perturbatively and expressed respectively as power series in a coupling parameter $g$ and $1/g$, a Pad\'{e} approximant embracing the two series can interpolate between these two limits and provide an accurate estimate of the system's behavior in the generally intractable intermediate coupling regime. The methodology and validity of this approach are illustrated by considering several stationary state problems in quantum mechanics., Comment: RevTeX4, 7 pages (including 7 tables); v4 typos corrected

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 $\mu$-distortion in both cases. We show that the $\mu$-parameter is generally enhanced relative to $\Lambda$CDM for interactions with neutrinos, and may be either enhanced or suppressed in the case of interactions with photons. The deviations from the $\Lambda$CDM signal are potentially within the sensitivity reach of a PRISM-like experiment if $\sigma_{\textrm{DM}-\gamma} \gtrsim 1.1\times10^{-30} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $\sigma_{\textrm{DM}-\nu} \gtrsim 4.8\times 10^{-32} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for time-independent cross sections, and $\sigma^{0}_{\textrm{DM}-\gamma} \gtrsim 1.8 \times 10^{-40} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $\sigma^{0}_{\textrm{DM}-\nu} \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 $\mu$-distortion signals differ from those of warm dark matter (no deviation from $\Lambda$CDM) and a suppressed primordial power spectrum (strongly suppressed or a negative $\mu$-parameter), demonstrating that CMB spectral distortion can potentially be used to distinguish between solutions to the small-scale crisis., Comment: V2: added references, minor edits, matches published version

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., Comment: 4 pages, 1 figure, matches version published in PRL

We comment on the recent paper by A. D. Dolgov, S. H. Hansen, S. Pastor and D. V. Semikoz (DHPS) [Astropart. Phys. {\bf 14}, 79 (2000)] on the generation of neutrino asymmetries from active-sterile neutrino oscillations. We demonstrate that the approximate asymmetry evolution equation obtained therein is an expansion, up to a minor discrepancy, of the well-established static approximation equation, valid only when the supposedly new higher order correction term is small. In the regime where this so-called ``back-reaction'' term is large and artificially terminates the asymmetry growth, their evolution equation ceases to be a faithful approximation to the Quantum Kinetic Equations (QKEs) simply because pure Mikheyev-Smirnov-Wolfenstein (MSW) transitions have been neglected. At low temperatures the MSW effect is the dominant asymmetry amplifier. Neither the static nor the DHPS approach contains this important physics. Therefore we conclude that the DHPS results have sufficient veracity at the onset of explosive asymmetry generation, but are invalid in the ensuing low temperature epoch where MSW conversions are able to enhance the asymmetry to values of order $0.2 - 0.37$. DHPS do claim to find a significant final asymmetry for very large $\delta m^2$ values. However, for this regime the effective potential they employed is not valid.

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., 20 pages, 8 figures, matches version published in JCAP 1310 (2013) 020

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., Comment: 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

I review how precision data from observations of the cosmic microwave background anisotropies and the large-scale structure distribution can be used to probe particle physics. Some examples are the absolute neutrino mass scale, dark radiation, light sterile neutrinos, QCD axions, WIMP annihilation, and dark sector interactions.

We consider the implications of a violation of the equivalence principle or of Lorentz invariance in the neutrino sector for the T-asymmetry $\Delta P_T \equiv P(\nu_{\alpha} \to \nu_{\beta}) - P(\nu_{\beta} \to \nu_{\alpha})$ in a three-flavour framework. We find that additional mixing due to these mechanisms, while obeying all present bounds, can lead to a substantial enhancement, suppression, and/or sign change in $\Delta P_T$ for the preferred energies and baselines of a neutrino factory. This in turn allows for the possibility of improving existing constraints by several orders of magnitude., Comment: 2 pages; Talk given at the 4th NuFact '02 Workshop (Neutrino Factories Based On Muon Storage Rings), 1-6 Jul 2002, London, England; To appear in proceedings

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)., 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

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., 33 pages, 14 of them appendix, 4 figures, matches published version in JCAP

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., 27 pages, 8 figures and 7 tables; corrected numerical error, match the submitted version

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., 37 pages, 14 figures and 7 tables. Replacement for matching the version accepted for publication

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, Comment: 14 pages, 6 figures, v2: minor changes, matches version accepted for publication in JCAP

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., 11 pages

Active-sterile neutrino oscillations provide a mechanism by which large differences in the neutrino and antineutrino number densities can be created in the early universe. The quantum kinetic equations are employed in the study of these neutrino asymmetries, which, when solved analytically in the adiabatic limit, generate physically transparent evolution equations that are very useful for the understanding of the nature of the asymmetry growth.

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 cosmic microwave background 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{sub s} and effective number N{sub 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{sub s} is in the sub-eV range.

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., 25 pages, 6 figures, matches version published in JCAP

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_\nu 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_\nu. 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_\nu < 1.16 eV, N_eff = 4.8 \pm 2.0., Comment: 18 pages, 4 figures; v2: references added, matches published version

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., Comment: 20 pages, 8 figures

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_��, 9 pages, 3 figures, uses iopart.cls; v2: one additional figure, references added, version accepted by JCAP

We study the evolution of a cosmological baryon asymmetry produced via leptogenesis by means of the full classical Boltzmann equations, without the assumption of kinetic equilibrium and including all quantum statistical factors. Beginning with the full mode equations we derive the usual equations of motion for the right-handed neutrino number density and integrated lepton asymmetry, and show explicitly the impact of each assumption on these quantities. For the first time, we investigate also the effects of scattering of the right-handed neutrino with the top quark to leading order in the Yukawa couplings by means of the full Boltzmann equations. We find that in our full Boltzmann treatment the final lepton asymmetry can be suppressed by as much as a factor of 1.5 in the weak wash-out regime (K1), the full Boltzmann treatment and the integrated approach give nearly identical final lepton asymmetries (within 10 % of each other at K>3). Finally, we show that the opposing effects of quantum statistics on decays/inverse decays and the scattering processes tend to reduce the net importance of scattering on leptogenesis in the full treatment compared to the integrated approach., 39 pages, 8 figures, typos corrected, replaced to match published version

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., 14 pages, 4 figures; v2: matches published version

We update our previous constraints on two-component hot dark matter (axions and neutrinos), including the recent WMAP 5-year data release. Marginalising over sum m_nu provides m_a < 1.02 eV (95% C.L.) for the axion mass. In the absence of axions we find sum m_nu < 0.63 eV (95% C.L.)., 4 pages, 1 figure, uses iopart.cls; v2 matches published version

We explore and compare the performances of two nonlinear correction and scale-dependent biasing models for the extraction of cosmological information from galaxy power spectrum data, especially in the context of beyond-LCDM cosmologies. The first model is the well known Q model, first applied in the analysis of 2dFGRS data. The second, the P model, is inspired by the halo model, in which nonlinear evolution and scale-dependent biasing are encapsulated in a single non-Poisson shot noise term. We find that while both models perform equally well in providing adequate correction for a range of galaxy clustering data in standard LCDM cosmology and in extensions with massive neutrinos, the Q model can give unphysical results in cosmologies containing a subdominant free-streaming dark matter whose temperature depends on the particle mass, e.g., relic thermal axions, unless a suitable prior is imposed on the correction parameter. This last case also exposes the danger of analytic marginalisation, a technique sometimes used in the marginalisation of nuisance parameters. In contrast, the P model suffers no undesirable effects, and is the recommended nonlinear correction model also because of its physical transparency., Comment: 21 pages, 8 figures, uses iopart.cls; v2: 22 pages, matches published version

I give an overview of the effects of neutrinos on cosmology, focussing in particular on the role played by neutrinos in the evolution of cosmological perturbations. I discuss how recent observations of the cosmic microwave background and the large‐scale structure of galaxies 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.

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., Comment: 22 pages, 5 figures, uses iopart.cls; v2: 23 pages, added references, to appear in JCAP

We revisit the issue of ripples in the primordial power spectra caused by trans-Planckian physics, and the potential for their detection by future cosmological probes. We find that for reasonably large values of the first slow-roll parameter epsilon (> 0.001), a positive detection of trans-Planckian ripples can be made even if the amplitude is as low as 10^-4. Data from the Large Synoptic Survey Telescope (LSST) and the proposed future 21 cm survey with the Fast Fourier Transform Telescope (FFTT) will be particularly useful in this regard. If the scale of inflation is close to its present upper bound, a scale of new physics as high as 0.2 M_Planck could lead to observable signatures., 20 pages, 3 figures, uses iopart.cls; v2: 21 pages, added references, to appear in JCAP

The existence of an unparticle sector, weakly coupled to the standard model, would have a profound impact on supernova (SN) physics. Emission of energy into the unparticle sector from the core of SN 1987A would have significantly shortened the observed neutrino burst. The unparticle interaction with nucleons, neutrinos, electrons and muons is constrained to be so weak that it is unlikely to provide any missing-energy signature at colliders. One important exception are models where scale invariance in the hidden sector is broken by the Higgs vacuum expectation value. In this case the SN emission is suppressed by threshold effects.