No Arabic abstract
We compute non-linear corrections to the matter power spectrum taking the time- and scale-dependent free-streaming length of neutrinos into account. We adopt a hybrid scheme that matches the full Boltzmann hierarchy to an effective two-fluid description at an intermediate redshift. The non-linearities in the neutrino component are taken into account by using an extension of the time-flow framework. We point out that this remedies a spurious behaviour that occurs when neglecting non-linear terms for neutrinos. This behaviour is related to how efficiently short modes decouple from long modes and can be traced back to the violation of momentum conservation if neutrinos are treated linearly. Furthermore, we compare our results at next to leading order to various other methods and quantify the accuracy of the fluid description. Due to the correct decoupling behaviour of short modes, the two-fluid scheme is a suitable starting point to compute higher orders in perturbations or for resummation methods.
We present a new method for simulating cosmologies that contain massive particles with thermal free streaming motion, such as massive neutrinos or warm/hot dark matter. This method combines particle and fluid descriptions of the thermal species to eliminate the shot noise known to plague conventional N-body simulations. We describe this method in detail, along with results for a number of test cases to validate our method, and check its range of applicability. Using this method, we demonstrate that massive neutrinos can produce a significant scale-dependence in the large-scale biasing of deep voids in the matter field. We show that this scale-dependence may be quantitatively understood using an extremely simple spherical expansion model which reproduces the behavior of the void bias for different neutrino parameters.
We present a comprehensive derivation of linear perturbation equations for different matter species, including photons, baryons, cold dark matter, scalar fields, massless and massive neutrinos, in the presence of a generic conformal coupling. Starting from the Lagrangians, we show how the conformal transformation affects the dynamics. In particular, we discuss how to incorporate consistently the scalar coupling in the equations of the Boltzmann hierarchy for massive neutrinos and the subsequent fluid approximations. We use the recently proposed K-mouflage model as an example to demonstrate the numerical implementation of our linear perturbation equations. K-mouflage is a new mechanism to suppress the fifth force between matter particles induced by the scalar coupling, but in the linear regime the fifth force is unsuppressed and can change the clustering of different matter species in different ways. We show how the CMB, lensing potential and matter power spectra are affected by the fifth force, and find ranges of K-mouflage parameters whose effects could be seen observationally. We also find that the scalar coupling can have the nontrivial effect of shifting the amplitude of the power spectra of the lensing potential and density fluctuations in opposite directions, although both probe the overall clustering of matter. This paper can serve as a reference for those who work on generic coupled scalar field cosmology, or those who are interested in the cosmological behaviour of the K-mouflage model.
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 formers 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.
Transformers with linearised attention (linear Transformers) have demonstrated the practical scalability and effectiveness of outer product-based Fast Weight Programmers (FWPs) from the 90s. However, the original FWP formulation is more general than the one of linear Transformers: a slow neural network (NN) continually reprograms the weights of a fast NN with arbitrary NN architectures. In existing linear Transformers, both NNs are feedforward and consist of a single layer. Here we explore new variations by adding recurrence to the slow and fast nets. We evaluate our novel recurrent FWPs (RFWPs) on two synthetic algorithmic tasks (code execution and sequential ListOps), Wikitext-103 language models, and on the Atari 2600 2D game environment. Our models exhibit properties of Transformers and RNNs. In the reinforcement learning setting, we report large improvements over LSTM in several Atari games. Our code is public.
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 $Lambda$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-$Lambda$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 5% accurate down to at least $kapprox 3 , h/{rm Mpc}$ for a modification to gravity of $|f_{rm R0}|leq 10^{-5}$ and for the total neutrino mass $M_ u equiv sum m_ u leq 0.15$ eV. We also find that the framework is 4(1)% consistent with the Bacco (EuclidEmulator2) emulator for $ u w$CDM cosmologies down to at least $k approx 3 , h$/Mpc. Finally, we compare against hydrodynamical simulations employing HMCode2020s baryonic feedback modelling on top of the halo model reaction. For $ u Lambda$CDM cosmologies we find 2% accuracy for $M_ u leq 0.48$eV down to at least $kapprox 5h$/Mpc. Similar accuracy is found when comparing to $ u w$CDM hydrodynamical simulations with $M_ u = 0.06$eV. This offers the first non-linear and theoretically general means of accurately including massive neutrinos for beyond-$Lambda$CDM cosmologies, and further suggests that baryonic effects can be reliably modelled independently of massive neutrino and dark energy physics. These extensions have been integrated into the publicly available ReACT code.