Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory

The linear matter power spectrum is an essential ingredient in all theoretical models for interpreting large-scale-structure observables. Although Boltzmann codes such as CLASS or CAMB are very efficient at computing the linear spectrum, the analysis of data usually requires 104-106 evaluations, which means this task can be the most computationally expensive aspect of data analysis. Here, we address this problem by building a neural network emulator that provides the linear theory (total and cold) matter power spectrum in about one millisecond with ≈0.2%(0.5%) accuracy over redshifts z ≤ 3 (z ≤ 9), and scales10-4 ≤ k [h Mpc-1] < 50. We train this emulator with more than 200,000 measurements, spanning a broad cosmological parameter space that includes massive neutrinos and dynamical dark energy. We show that the parameter range and accuracy of our emulator is enough to get unbiased cosmological constraints in the analysis of a Euclid-like weak lensing survey. Complementing this emulator, we train 15 other emulators for the cross-spectra of various linear fields in Eulerian space, as predicted by 2nd-order Lagrangian Perturbation theory, which can be used to accelerate perturbative bias descriptions of galaxy clustering. Our emulators are specially designed to be used in combination with emulators for the nonlinear matter power spectrum and for baryonic effects, all of which are publicly available at http://www.dipc.org/bacco.

Seesaw, coherence and neutrino oscillations

We present a prescription for consistently constructing non-Fock coherent flavour neutrino states within the framework of the seesaw mechanism, and establish that the physical vacuum of massive neutrinos is a condensate of Standard Model massless neutrino states. The coherent states, involving a finite number of massive states, are derived by constructing their creation operator. This construction fulfills automatically the key requirement of coherence for the oscillations of particles to occur. We comment on the inherent non-unitarity of the oscillation probability induced by the requirement of coherence.

Clustering in massive neutrino cosmologies via Eulerian Perturbation Theory

We introduce an Eulerian Perturbation Theory to study the clustering of tracers for cosmologies in the presence of massive neutrinos. Our approach is based on mapping recently-obtained Lagrangian Perturbation Theory results to the Eulerian framework. We add Effective Field Theory counterterms, IR-resummations and a biasing scheme to compute the one-loop redshift-space power spectrum. To assess our predictions, we compare the power spectrum multipoles against synthetic halo catalogues from the QUIJOTE simulations, finding excellent agreement on scales k ≲ 0.25 h Mpc-1. One can obtain the same fitting accuracy using higher wave-numbers, but then the theory fails to give a correct estimation of the linear bias parameter. We further discuss the implications for the tree-level bispectrum. Finally, calculating loop corrections is computationally costly, hence we derive an accurate approximation wherein we retain only the main features of the kernels, as produced by changes to the growth rate. As a result, we show how FFTLog methods can be used to further accelerate the loop computations with these reduced kernels.

On the road to percent accuracy V: The non-linear power spectrum beyond ΛCDM with massive neutrinos and baryonic feedback

In the context of forthcoming galaxy surveys, to ensure unbiased constraints on cosmology and gravity when using non-linear structure information, percent-level accuracy is required when modelling the power spectrum. This calls for frameworks that can accurately capture the relevant physical effects, while allowing for deviations from ΛCDM. Massive neutrino and baryonic physics are two of the most relevant such effects. We present an integration of the halo model reaction frameworks for massive neutrinos and beyond-ΛCDM cosmologies. The integrated halo model reaction, combined with a pseudo power spectrum modelled by HMCode2020 is then compared against N-body simulations that include both massive neutrinos and an f(R) modification to gravity. We find that the framework is 4 per cent accurate down to at least k ≈ 3 h Mpc−1 for a modification to gravity of |fR0| ≤ 10−5 and for the total neutrino mass Mν ≡ ∑mν ≤ 0.15 eV. We also find that the framework is 4 per cent consistent with EuclidEmulator2 as well as the Bacco emulator for most of the considered νwCDM cosmologies down to at least k ≈ 3 h Mpc−1. Finally, we compare against hydrodynamical simulations employing HMCode2020’s baryonic feedback modelling on top of the halo model reaction. For νΛCDM cosmologies we find 2 per cent accuracy for Mν ≤ 0.48 eV down to at least k ≈ 5h Mpc−1. Similar accuracy is found when comparing to νwCDM hydrodynamical simulations with Mν = 0.06 eV. This offers the first non-linear, theoretically general means of accurately including massive neutrinos for beyond-ΛCDM cosmologies, and further suggests that baryonic, massive neutrino and dark energy physics can be reliably modelled independently.

Probing neutrino magnetic moment at the Jinping neutrino experiment

Neutrino magnetic moment (νMM) is an important property of massive neutrinos. The recent anomalous excess at few keV electronic recoils observed by the XENON1T collaboration might indicate a ∼ 2.2 × 10−11μB effective neutrino magnetic moment ($$ {\mu}_{\nu}^{\mathrm{eff}} $$ μ ν eff ) from solar neutrinos. Therefore, it is essential to carry out the νMM searches at a different experiment to confirm or exclude such a hypothesis. We study the feasibility of doing νMM measurement with 4 kton fiducial mass at Jinping neutrino experiment (Jinping) using electron recoil data from both natural and artificial neutrino sources. The sensitivity of $$ {\mu}_{\nu}^{\mathrm{eff}} $$ μ ν eff can reach < 1.2 × 10−11μB at 90% C.L. with 10-year data taking of solar neutrinos. Besides the abundance of the intrinsic low energy background 14C and 85Kr in the liquid scintillator, we find the sensitivity to νMM is highly correlated with the systematic uncertainties of pp and 85Kr. Reducing systematic uncertainties (pp and 85Kr) and the intrinsic background (14C and 85Kr) can help to improve sensitivities below these levels and reach the region of astrophysical interest. With a 3 mega-Curie (MCi) artificial neutrino source 51Cr installed at Jinping neutrino detector for 55 days, it could give us a sensitivity to the electron neutrino magnetic moment ($$ {\mu}_{\nu_e} $$ μ ν e ) with < 1.1 × 10−11μB at 90% C.L. . With the combination of those two measurements, the flavor structure of the neutrino magnetic moment can be also probed at Jinping.

Tentative sensitivity of future 0νββ-decay experiments to neutrino masses and Majorana CP phases

In the near future, the neutrinoless double-beta (0νββ) decay experiments will hopefully reach the sensitivity of a few meV to the effective neutrino mass |mββ|. In this paper, we tentatively examine the sensitivity of future 0νββ-decay experiments to neutrino masses and Majorana CP phases by following the Bayesian statistical approach. Provided experimental setups corresponding to the experimental sensitivity of |mββ| ≃ 1 meV, the null observation of 0νββ decays in the case of normal neutrino mass ordering leads to a very competitive bound on the lightest neutrino mass m1. Namely, the 95% credible interval in the Bayesian approach turns out to be 1.6 meV ≲ m1 ≲ 7.3 meV or 0.3 meV ≲ m1 ≲ 5.6 meV when the uniform prior on m1/eV or on log10(m1/eV) is adopted. Moreover, one of two Majorana CP phases is strictly constrained, i.e., 140° ≲ ρ ≲ 220° for both scenarios of prior distributions of m1. In contrast, if a relatively worse experimental sensitivity of |mββ| ≃ 10 meV is assumed, the constraint on the lightest neutrino mass becomes accordingly 0.6 meV ≲ m1 ≲ 26 meV or 0 ≲ m1 ≲ 6.1 meV, while two Majorana CP phases will be essentially unconstrained. In the same statistical framework, the prospects for the determination of neutrino mass ordering and the discrimination between Majorana and Dirac nature of massive neutrinos in the 0νββ-decay experiments are also discussed. Given the experimental sensitivity of |mββ| ≃ 10 meV (or 1 meV), the strength of evidence to exclude the Majorana nature under the null observation of 0νββ decays is found to be inconclusive (or strong), no matter which of two priors on m1 is taken.