|
Barenboim, G., Blinov, N., & Stebbins, A. (2021). Smallest remnants of early matter domination. J. Cosmol. Astropart. Phys., 12(12), 026–50pp.
Abstract: The evolution of the universe prior to Big Bang Nucleosynthesis could have gone through a phase of early matter domination which enhanced the growth of small-scale dark matter structure. If this period was long enough, self-gravitating objects formed prior to reheating. We study the evolution of these dense early halos through reheating. At the end of early matter domination, the early halos undergo rapid expansion and eventually eject their matter. We find that this process washes out structure on scales much larger than naively expected from the size of the original halos. We compute the density profiles of the early halo remnants and use them to construct late-time power spectra that include these non-linear effects. We evolve the resulting power spectrum to estimate the properties of microhalos that would form after matter-radiation equality. Surprisingly, cosmologies with a short period of early matter domination lead to an earlier onset of microhalo formation compared to those with a long period. In either case, dark matter structure formation begins much earlier than in the standard cosmology, with most dark matter bound in microhalos in the late universe.
|
|
|
Figueroa, D. G., Florio, A., Torrenti, F., & Valkenburg, W. (2023). CosmoLattice: A modern code for lattice simulations of scalar and gauge field dynamics in an expanding universe. Comput. Phys. Commun., 283, 108586–13pp.
Abstract: This paper describes CosmoGattice, a modern package for lattice simulations of the dynamics of interacting scalar and gauge fields in an expanding universe. CosmoGattice incorporates a series of features that makes it very versatile and powerful: i) it is written in C++ fully exploiting the object oriented programming paradigm, with a modular structure and a clear separation between the physics and the technical details, ii) it is MPI-based and uses a discrete Fourier transform parallelized in multiple spatial dimensions, which makes it specially appropriate for probing scenarios with well -separated scales, running very high resolution simulations, or simply very long ones, iii) it introduces its own symbolic language, defining field variables and operations over them, so that one can introduce differential equations and operators in a manner as close as possible to the continuum, iv) it includes a library of numerical algorithms, ranging from O(delta t(2)) to O(delta t(10)) methods, suitable for simulating global and gauge theories in an expanding grid, including the case of 'self-consistent' expansion sourced by the fields themselves. Relevant observables are provided for each algorithm (e.g. energy densities, field spectra, lattice snapshots) and we note that, remarkably, all our algorithms for gauge theories (Abelian or non-Abelian) always respect the Gauss constraint to machine precision. Program summary Program Title:: CosmoGattice CPC Library link to program files: https://doi .org /10 .17632 /44vr5xssc6 .1 Developer's repository link: http://github .com /cosmolattice /cosmolattice Licensing provisions: MIT Programming language: C++, MPI Nature of problem: The phenomenology of high energy physics in the early universe is typically characterized by non-linear dynamics, which cannot be captured accurately with analytical techniques. In order to fully understand the non-linearities developed in a given scenario, one needs to carry out lattice simulations. A number of public packages for lattice simulations have appeared over the years, but most of them are only capable of simulating scalar fields. However, realistic models of particle physics do contain other kind of field species, such as (Abelian or non-Abelian) gauge fields, whose non-linear dynamics can also play a relevant role in the early universe. Tensor modes representing gravitational waves are also naturally expected in many scenarios. Solution method: CosmoGattice represents a modern code for lattice simulations of scalar-gauge field theories in an expanding universe. It allows for the simulation of the evolution of interacting (singlet) scalar fields, charged scalar fields under U(1) and/or SU(2) gauge groups, and the corresponding associated Abelian and/or non-Abelian gauge fields. From version 1.1 onward, CosmoGattice also allows to simulate the production of gravitational waves. Simulations can be done either in a flat space-time background, or in a homogeneous and isotropic (spatially flat) expanding FLRW background. CosmoGattice provides symplectic integrators, with accuracy ranging from O (delta t(2)) up to O(delta t(10)), to simuate the non-linear dynamics of the appropriate fields in comoving three-dimensional lattices. The code is parallelized with MPI, and uses a discrete Fourier Transform parallelized in multiple spatial dimensions, which makes it a very powerful code for probing physical problems with well-separated scales. Moreover, the code has been designed as a `platform' to implement any system of dynamical equations suitable for discretization on a lattice.
|
|
|
Bennett, J. J., Buldgen, G., de Salas, P. F., Drewes, M., Gariazzo, S., Pastor, S., et al. (2021). Towards a precision calculation of the effective number of neutrinos N-eff in the Standard Model. Part II. Neutrino decoupling in the presence of flavour oscillations and finite-temperature QED. J. Cosmol. Astropart. Phys., 04(4), 073–33pp.
Abstract: We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, N-eff(SM), that quantifies the cosmological neutrinoto-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to O(e(3)), where e is the elementary electric charge, and a full evaluation of the neutrino-neutrino collision integral. We provide furthermore a detailed assessment of the uncertainties in the benchmark N(eff)(SM )value, through testing the value's dependence on (i) optional approximate modelling of the weak collision integrals, (ii) measurement errors in the physical parameters of the weak sector, and (iii) numerical convergence, particularly in relation to momentum discretisation. Our new, recommended standard-model benchmark is N-eff(SM) 3.0440 +/- 0.0002, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure (vertical bar delta N-eff vertical bar similar to 10(-4)), augmented by measurement errors in the solar mixing angle sin(2) theta(12) (vertical bar delta N-eff vertical bar similar to 10(-4)).
|
|
|
Giare, W., Di Valentino, E., Melchiorri, A., & Mena, O. (2021). New cosmological bounds on hot relics: axions and neutrinos. Mon. Not. Roy. Astron. Soc., 505(2), 2703–2711.
Abstract: Axions, if realized in nature, can be copiously produced in the early universe via thermal processes, contributing to the mass-energy density of thermal hot relics. In light of the most recent cosmological observations, we analyse two different thermal processes within a realistic mixed hot dark matter scenario which includes also massive neutrinos. Considering the axion-gluon thermalization channel, we derive our most constraining bounds on the hot relic masses m(a) < 7.46 eV and Sigma m(nu) < 0.114 eV both at 95 percent CL; while studying the axion-pion scattering, without assuming any specific model for the axion-pion interactions, and remaining in the range of validity of the chiral perturbation theory, our most constraining bounds are improved to m(a) < 0.91 eV and Sigma m(nu) < 0.105 eV, both at 95 percent CL. Interestingly, in both cases, the total neutrino mass lies very close to the inverted neutrino mass ordering prediction. If future terrestrial double beta decay and/or long-baseline neutrino experiments find that the nature mass ordering is the inverted one, this could rule out a wide region in the currently allowed thermal axion window. Our results therefore, strongly support multi messenger searches of axions and neutrino properties, together with joint analyses of their expected sensitivities.
|
|
|
Figueroa, D. G., Florio, A., Torrenti, F., & Valkenburg, W. (2021). The art of simulating the early universe. Part I. Integration techniques and canonical cases. J. Cosmol. Astropart. Phys., 04(4), 035–108pp.
Abstract: We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce the basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of O(delta t(2)) accuracy like staggered leapfrog and Verlet integration, to Runge-Kutta methods up to O(delta t(4)) accuracy, and the Yoshida and Gauss-Legendre higher-order integrators, accurate up to O(delta t(10)) We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in 3+1 dimensions, including the case of 'self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: i) Interacting scalar fields, ii) Abelian U(1) gauge theories, and iii) Non-Abelian SU(2) gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from O(delta t(2)) up to O(delta t(10)) For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, including when 'self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under SU(2) x U(1). We note that the present manuscript is meant to be part of the theoretical basis for the code CosmoLattice, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, publicly available at http://www.cosrnolattice.net.
|
|