|
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.
|
|
|
Forconi, M., Giare, W., Mena, O., Ruchika, Di Valentino, E., Melchiorri, A., et al. (2024). A double take on early and interacting dark energy from JWST. J. Cosmol. Astropart. Phys., 05(5), 097–37pp.
Abstract: The very first light captured by the James Webb Space Telescope (JWST) revealed a population of galaxies at very high redshifts more massive than expected in the canonical Lambda CDM model of structure formation. Barring, among others, a systematic origin of the issue, in this paper, we test alternative cosmological perturbation histories. We argue that models with a larger matter component ohm m and/or a larger scalar spectral index n s can substantially improve the fit to JWST measurements. In this regard, phenomenological extensions related to the dark energy sector of the theory are appealing alternatives, with Early Dark Energy emerging as an excellent candidate to explain (at least in part) the unexpected JWST preference for larger stellar mass densities. Conversely, Interacting Dark Energy models, despite producing higher values of matter clustering parameters such as sigma 8 , are generally disfavored by JWST measurements. This is due to the energy -momentum flow from the dark matter to the dark energy sector, implying a smaller matter energy density. Upcoming observations may either strengthen the evidence or falsify some of these appealing phenomenological alternatives to the simplest Lambda CDM picture.
|
|
|
Blennow, M., Fernandez-Martinez, E., Mena, O., Redondo, J., & Serra, E. P. (2012). Asymmetric Dark Matter and Dark Radiation. J. Cosmol. Astropart. Phys., 07(7), 022–23pp.
Abstract: Asymmetric Dark Matter (ADM) models invoke a particle-antiparticle asymmetry, similar to the one observed in the Baryon sector, to account for the Dark Matter (DM) abundance. Both asymmetries are usually generated by the same mechanism and generally related, thus predicting DM masses around 5 GeV in order to obtain the correct density. The main challenge for successful models is to ensure efficient annihilation of the thermally produced symmetric component of such a light DM candidate without violating constraints from collider or direct searches. A common way to overcome this involves a light mediator, into which DM can efficiently annihilate and which subsequently decays into Standard Model particles. Here we explore the scenario where the light mediator decays instead into lighter degrees of freedom in the dark sector that act as radiation in the early Universe. While this assumption makes indirect DM searches challenging, it leads to signals of extra radiation at BBN and CMB. Under certain conditions, precise measurements of the number of relativistic species, such as those expected from the Planck satellite, can provide information on the structure of the dark sector. We also discuss the constraints of the interactions between DM and Dark Radiation from their imprint in the matter power spectrum.
|
|