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.
|
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)).
|
Cosme, C., Figueroa, D. G., & Loayza, N. (2023). Gravitational wave production from preheating with trilinear interactions. J. Cosmol. Astropart. Phys., 05(5), 023–30pp.
Abstract: We investigate the production of gravitational waves (GWs) during preheating with monomial/polynomial inflationary potentials, considering a trilinear coupling & phi;x2 between a singlet inflaton & phi; and a daughter scalar field x. For sufficiently large couplings, the trilinear interaction leads to an exponential production of x particles and, as a result, a large stochastic GW background (SGWB) is generated throughout the process. We study the linear and non-linear dynamics of preheating with lattice simulations, following the production of GWs through all relevant stages. We find that large couplings lead to SGWBs with amplitudes today that can reach up to h2 �(0) GW <^> 5 & BULL; 10-9. These backgrounds are however peaked at high frequencies fp > 5 & BULL; 106 Hz, which makes them undetectable by current/planned GW observatories. As the amount of GWs produced is in any case remarkable, we discuss the prospects for probing the SGWB indirectly by using constraints on the effective number of relativistic species in the universe Neff.
|
Domcke, V., Ema, Y., & Sandner, S. (2024). Perturbatively including inhomogeneities in axion inflation. J. Cosmol. Astropart. Phys., 03(3), 019–24pp.
Abstract: Axion inflation, i.e. an axion-like inflaton coupled to an Abelian gauge field through a Chern-Simons interaction, comes with a rich and testable phenomenology. This is particularly true in the strong backreaction regime, where the gauge field production heavily impacts the axion dynamics. Lattice simulations have recently demonstrated the importance of accounting for inhomogeneities of the axion field in this regime. We propose a perturbative scheme to account for these inhomogeneities while maintaining high computational efficiency. Our goal is to accurately capture deviations from the homogeneous axion field approximation within the perturbative regime as well as self -consistently determine the onset of the nonperturbative regime.
|
Hernandez, P., Lopez-Pavon, J., Rius, N., & Sandner, S. (2022). Bounds on right-handed neutrino parameters from observable leptogenesis. J. High Energy Phys., 12(12), 012–58pp.
Abstract: We revisit the generation of a matter-antimatter asymmetry in the minimal extension of the Standard Model with two singlet heavy neutral leptons (HNL) that can explain neutrino masses. We derive an accurate analytical approximation to the solution of the complete linearized set of kinetic equations, which exposes the non-trivial parameter dependencies in the form of parameterization-independent CP invariants. The identification of various washout regimes relevant in different regions of parameter space sheds light on the relevance of the mass corrections in the interaction rates and clarifies the correlations of baryogenesis with other observables. In particular, by requiring that the measured baryon asymmetry is reproduced, we derive robust upper or lower bounds on the HNL mixings depending on their masses, and constraints on their flavour structure, as well as on the CP-violating phases of the PMNS mixing matrix, and the amplitude of neutrinoless double-beta decay. We also find certain correlations between low and high scale CP phases. Especially emphasizing the testable part of the parameter space we demonstrate that our findings are in very good agreement with numerical results. The methods developed in this work can help in exploring more complex scenarios.
|
