Source author record

Thomas Tram

Thomas Tram appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

$ν\scriptstyle{\rm CO}N{\rm CEPT}$: Cosmological neutrino simulations from the non-linear Boltzmann hierarchy

In this paper the non-linear effect of massive neutrinos on cosmological structures is studied in a conceptually new way. We have solved the non-linear continuity and Euler equations for the neutrinos on a grid in real space in $N$-body simulations, and closed the Boltzmann hierarchy at the non-linear Euler equation using the stress and pressure perturbations from linear theory. By comparing with state-of-the art cosmological neutrino simulations, we are able to simulate the non-linear neutrino power spectrum very accurately. This translates into a negligible error in the matter power spectrum, and so our $\scriptstyle{\rm CO}N{\rm CEPT}$ code is ideally suited for extracting the neutrino mass from future high precision non-linear observational probes such as Euclid.

preprint2022arXiv

Decaying warm dark matter revisited

Decaying dark matter models provide a physically motivated way of channeling energy between the matter and radiation sectors. In principle, this could affect the predicted value of the Hubble constant in such a way as to accommodate the discrepancies between CMB inferences and local measurements of the same. Here, we revisit the model of warm dark matter decaying non-relativistically to invisible radiation. In particular, we rederive the background and perturbation equations starting from a decaying neutrino model and describe a new, computationally efficient method of computing the decay product perturbations up to large multipoles. We conduct MCMC analyses to constrain all three model parameters, for the first time including the mass of the decaying species, and assess the ability of the model to alleviate the Hubble and $σ_8$ tensions, the latter being the discrepancy between the CMB and weak gravitational lensing constraints on the amplitude of matter fluctuations on an $8 h^{-1}$ Mpc$^{-1}$ scale. We find that the model reduces the $H_0$ tension from $\sim 4 σ$ to $\sim 3 σ$ and neither alleviates nor worsens the $S_8 \equiv σ_8 (Ω_m/0.3)^{0.5}$ tension, ultimately showing only mild improvements with respect to $Λ$CDM. However, the values of the model-specific parameters favoured by data is found to be well within the regime of relativistic decays where inverse processes are important, rendering a conclusive evaluation of the decaying warm dark matter model open to future work.

preprint2020arXiv

MeV-scale reheating temperature and cosmological production of light sterile neutrinos

We investigate how sterile neutrinos with a range of masses influence cosmology in MeV-scale reheating temperature scenarios. By computing the production of sterile neutrinos through the combination of mixing and scattering in the early Universe, we find that light sterile neutrinos, with masses and mixings as inferred from short-baseline neutrino oscillation experiments, are consistent with big-bang nucleosynthesis (BBN) and cosmic microwave background (CMB) radiation for the reheating temperature of ${\cal O}(1)$ MeV if the parent particle responsible for reheating decays into electromagnetic components (radiative decay). In contrast, if the parent particle mainly decays into hadrons (hadronic decay), the bound from BBN becomes more stringent. In this case, the existence of the light sterile neutrinos can be cosmologically excluded, depending on the mass and the hadronic branching ratio of the parent particle.

preprint2016arXiv

Cosmological $N$-body simulations including radiation perturbations

Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian approximation, thus discarding general relativistic effects such as the coupling between matter and radiation ($\equiv$ photons and neutrinos). In this paper we investigate novel hybrid simulations which incorporate interactions between radiation and matter to the leading order in General Relativity, whilst evolving the matter dynamics in full non-linearity according to Newtonian theory. Our hybrid simulations come with a relativistic space-time and make it possible to investigate structure formation in a unified framework. In the present work we focus on simulations initialized at $z=99$, and show that the extracted matter power spectrum receives up to $3\%$ corrections on very large scales through radiation. Our numerical findings compare favourably with linear analytical results from Fidler et al. (2016), from which we deduce that there cannot be any significant non-linear mode-coupling induced through linear radiation corrections.

preprint2016arXiv

Pseudoscalar - sterile neutrino interactions: reconciling the cosmos with neutrino oscillations

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

preprint2016arXiv

Reconciling CMB and structure growth measurements with dark energy interactions

We study a coupled quintessence model with pure momentum exchange and present the effects of such an interaction on the Cosmic Microwave Background (CMB) and matter power spectrum. For a wide range of negative values of the coupling parameter $β$ structure growth is suppressed and the model can reconcile the tension between Cosmic Microwave Background observations and structure growth inferred from cluster counts. We find that this model is as good as $Λ$CDM for CMB and baryon acoustic oscillation (BAO) data, while the addition of cluster data makes the model strongly preferred, improving the best-fit $χ^2$-value by more than $16$.

preprint2016arXiv

Relativistic Interpretation of Newtonian Simulations for Cosmic Structure Formation

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

preprint2016arXiv

The Intrinsic Matter Bispectrum in $Λ$CDM

We present a fully relativistic calculation of the matter bispectrum at second order in cosmological perturbation theory assuming a Gaussian primordial curvature perturbation. For the first time we perform a full numerical integration of the bispectrum for both baryons and cold dark matter using the second-order Einstein-Boltzmann code, SONG. We review previous analytical results and provide an improved analytic approximation for the second-order kernel in Poisson gauge which incorporates Newtonian nonlinear evolution, relativistic initial conditions, the effect of radiation at early times and the cosmological constant at late times. Our improved kernel provides a percent level fit to the full numerical result at late times for most configurations, including both equilateral shapes and the squeezed limit. We show that baryon acoustic oscillations leave an imprint in the matter bispectrum, making a significant impact on squeezed shapes.

preprint2015arXiv

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

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.

preprint2015arXiv

Cosmology with self-interacting sterile neutrinos and dark matter - A pseudoscalar model

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. Such sterile neutrinos are incompatible with cosmology because they suppress structure formation unless they can be prevented from thermalising in the early Universe or removed by subsequent decay or annihilation. Here we present a novel scenario in which both sterile neutrinos and dark matter are coupled to a new, light pseudoscalar. This can prevent thermalisation of sterile neutrinos and make dark matter sufficiently self-interacting to have an impact on galactic dynamics and possibly resolve some of the known problems with the standard cold dark matter scenario. Even more importantly it leads to a strongly self-interacting plasma of sterile neutrinos and pseudoscalars at late times and provides an excellent fit to CMB data. The usual cosmological neutrino mass problem is avoided by sterile neutrino annihilation to pseudoscalars. The preferred value of $H_0$ is substantially higher than in standard $Λ$CDM and in much better agreement with local measurements.

preprint2015arXiv

General relativistic corrections to $N$-body simulations and the Zel'dovich approximation

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

preprint2015arXiv

Robustness of cosmic neutrino background detection in the cosmic microwave background

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

preprint2015arXiv

Sterile neutrinos with pseudoscalar self-interactions and cosmology

Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave Background and large scale structure. This paper discusses how charging sterile neutrinos under a new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available data which is way better than the standard $Λ$CDM model with one additional fully thermalized sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the locally measured value.

preprint2015arXiv

The Cosmic Linear Anisotropy Solving System (CLASS) IV: Efficient implementation of non-cold relics

We present a new flexible, fast and accurate way to implement massive neutrinos, warm dark matter and any other non-cold dark matter relics in Boltzmann codes. For whatever analytical or numerical form of the phase-space distribution function, the optimal sampling in momentum space compatible with a given level of accuracy is automatically found by comparing quadrature methods. The perturbation integration is made even faster by switching to an approximate viscous fluid description inside the Hubble radius, which differs from previous approximations discussed in the literature. When adding one massive neutrino to the minimal cosmological model, CLASS becomes just 1.5 times slower, instead of about 5 times in other codes (for fixed accuracy requirements). We illustrate the flexibility of our approach by considering a few examples of standard or non-standard neutrinos, as well as warm dark matter models.

preprint2014arXiv

A note of clarification: BICEP2 and Planck are not in tension

The apparent discrepancy between the value of the tensor-to-scalar ratio reported by the BICEP2 collaboration, $r = 0.20^{+0.07}_{-0.05}$ at 68% CL, and the Planck upper limit, $r < 0.11$ at 95% CL, has attracted a great deal of attention. In this short note, we show that this discrepancy is mainly due to an `apples to oranges' comparison. The result reported by BICEP2 was measured at a pivot scale $k_* = 0.05$ Mpc$^{-1}$, assuming $n_t = 0$, whereas the Planck limit was provided at $k_* = 0.002$ Mpc$^{-1}$, assuming the slow-roll consistency relation $n_t = -r/8$. One should obviously compare the BICEP2 and Planck results under the same circumstances. By imposing $n_t = 0$, the Planck constraint at $k_* = 0.05$ Mpc$^{-1}$ becomes $r < 0.135$ at $95%$ CL, which can be compared directly with the BICEP2 result. Once a plausible dust contribution to the BICEP2 signal is taken into account (DDM2 model), $r$ is reduced to $r = 0.16^{+0.06}_{-0.05}$ and the discrepancy becomes of order $1.3σ$ only.

preprint2014arXiv

Fast and accurate CMB computations in non-flat FLRW universes

We present a new method for calculating CMB anisotropies in a non-flat Friedmann universe, relying on a very stable algorithm for the calculation of hyperspherical Bessel functions, that can be pushed to arbitrary precision levels. We also introduce a new approximation scheme which gradually takes over in the flat space limit, and significant speeds-up calculations. Our method is implemented in the Boltzmann code CLASS. It can be used to benchmark the accuracy of the CAMB code in curved space, which is found to match expectations. For default precision settings, corresponding to 0.1% for scalar temperature spectra and 0.2% for scalar polarisation spectra, our code is two to three times faster, depending on curvature. We also simplify the temperature and polarisation source terms significantly, so the different contributions to the $C_\ell$'s are easy to identify inside the code.

preprint2014arXiv

Strongest model-independent bound on the lifetime of Dark Matter

Dark Matter is essential for structure formation in the late Universe so it must be stable on cosmological time scales. But how stable exactly? Only assuming decays into relativistic particles, we report an otherwise model independent bound on the lifetime of Dark Matter using current cosmological data. Since these decays affect only the low-$\ell$ multipoles of the CMB, the Dark Matter lifetime is expected to correlate with the tensor-to-scalar ratio $r$ as well as curvature $Ω_k$. We consider two models, including $r$ and $r+Ω_k$ respectively, versus data from Planck, WMAP, WiggleZ and Baryon Acoustic Oscillations, with or without the BICEP2 data (if interpreted in terms of primordial gravitational waves). This results in a lower bound on the lifetime of CDM given by 160Gyr (without BICEP2) or 200Gyr (with BICEP2) at 95% confidence level.

preprint2013arXiv

Can active-sterile neutrino oscillations lead to chaotic behavior of the cosmological lepton asymmetry?

While the cosmic baryon asymmetry has been measured at high accuracy to be 6.1*10^-10, a corresponding lepton asymmetry could be as large as 10^-2 if it hides in the neutrino sector. It has been known for some time that a large asymmetry can be generated from a small initial asymmetry given the existence of a sterile neutrino with a mass less than the mass of the active neutrino. While the magnitude of the final lepton asymmetry is deterministic, its sign has been conjectured to be chaotic in nature. This has been proven in the single momentum approximation, also known as the quantum rate equations, but has up to now not been established using the full momentum dependent quantum kinetic equations. Here we investigate this problem by solving the quantum kinetic equations for a system of 1 active and 1 sterile neutrino on an adaptive grid. We show that by increasing the resolution oscillations in the lepton asymmetry are eliminated and that the sign of the lepton asymmetry is in fact deterministic. This paper also serves as a launch paper for the adaptive solver LASAGNA which is available at http://users-phys.au.dk/steen.

preprint2013arXiv

How secret interactions can reconcile sterile neutrinos with cosmology

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. However, such neutrinos seem incompatible with cosmology because they have too large an impact on cosmic structure formation. Here we show that new interactions in the sterile neutrino sector can prevent their production in the early Universe and reconcile short baseline oscillation experiments with cosmology.

preprint2013arXiv

The effective gravitational decoupling between dark matter and the CMB

We present a detailed and self-contained analytical derivation of the evolution of sub-horizon cosmological perturbations before decoupling, based on previous work by S. Weinberg. These solutions are valid in the minimal LCDM scenario, to first order in perturbation theory, in the tight-coupling limit and neglecting neutrino shear stress. We compare them to exact numerical solutions computed by a Boltzmann code, and we find the two to be in very good agreement. The analytic solutions show explicitly that CDM and the baryon-photon fluid effectively behave as separate self-gravitating fluids until the epoch of baryon drag. This in turn leads to the surprising conclusion that the CMB is much less sensitive to the clustering properties of minimally coupled Dark Matter models than what would be naively expected.

preprint2012arXiv

Thermalisation of light sterile neutrinos in the early universe

Recent cosmological data favour additional relativistic degrees of freedom beyond the three active neutrinos and photons, often referred to as 'dark' radiation. Light sterile neutrinos is one of the prime candidates for such additional radiation. However, constraints on sterile neutrinos based on the current cosmological data have been derived using simplified assumptions about thermalisation of the sterile neutrino at the Big Bang Nucleosynthesis (BBN) epoch. These assumptions are not necessarily justified and here we solve the full quantum kinetic equations in the (1 active + 1 sterile) scenario and derive the number of thermalised species just before BBN begins (T~1MeV) for null (L=0) and large (L=0.01) initial lepton asymmetry and for a range of possible mass-mixing parameters. We find that the full thermalisation assumption during the BBN epoch is justified for initial small lepton asymmetry only. Partial or null thermalisation occurs when the initial lepton asymmetry is large.

preprint2011arXiv

Sommerfeld Enhancement of DM Annihilation: Resonance Structure, Freeze-Out and CMB Spectral Bound

In the last few years there has been some interest in WIMP Dark Matter models featuring a velocity dependent cross section through the Sommerfeld enhancement mechanism, which is a nonrelativistic effect due to massive bosons in the dark sector. In the first part of this article, we find analytic expressions for the boost factor for three different model potentials, the Coulomb potential, the spherical well and the spherical cone well and compare with the numerical solution of the Yukawa potential. We find that the resonance pattern of all the potentials can be cast into the same universal form. In the second part of the article we perform a detailed computation of the Dark Matter relic density for models having Sommerfeld enhancement by solving the Boltzmann equation numerically. We calculate the expected distortions of the CMB blackbody spectrum from WIMP annihilations and compare these to the bounds set by FIRAS. We conclude that only a small part of the parameter space can be ruled out by the FIRAS observations.

preprint2011arXiv

The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes

Boltzmann codes are used extensively by several groups for constraining cosmological parameters with Cosmic Microwave Background and Large Scale Structure data. This activity is computationally expensive, since a typical project requires from 10'000 to 100'000 Boltzmann code executions. The newly released code CLASS (Cosmic Linear Anisotropy Solving System) incorporates improved approximation schemes leading to a simultaneous gain in speed and precision. We describe here the three approximations used by CLASS for basic LambdaCDM models, namely: a baryon-photon tight-coupling approximation which can be set to first order, second order or to a compromise between the two; an ultra-relativistic fluid approximation which had not been implemented in public distributions before; and finally a radiation streaming approximation taking reionisation into account.