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.

Abstract: We discuss the interpretation of the detected signal by pulsar timing array (PTA) observations as a gravitational wave background of cosmological origin. We combine NANOGrav 15-years and EPTADR2new datasets and confront them against backgrounds from supermassive black hole binaries (SMBHBs), and cosmological signals from inflation, cosmic (super)strings, first-order phase transitions, Gaussian and non-Gaussian large scalar fluctuations, and audible axions. We find that scalar-induced, and to a lesser extent audible axion and cosmic superstring signals, provide a better fit than SMBHBs. These results depend, however, on modeling assumptions, so further data and analysis are needed to reach robust conclusions. Independently of the signal origin, the data strongly constrain the parameter space of cosmological signals, for example, setting an upper bound on primordial non-Gaussianity at PTA scales as jfnlj less than or similar to 2.34 at 95% C.L.

Abstract: We consider the phenomenological implications of the violation of the Pauli exclusion principle for neutrinos, focusing on cosmological observables such as the spectrum of Cosmic Microwave Background anisotropies, Baryon Acoustic Oscillations and the primordial abundances of light elements. Neutrinos that behave (at least partly) as bosonic particles have a modified equilibrium distribution function that implies a different influence on the evolution of the Universe that, in the case of massive neutrinos, can not be simply parametrized by a change in the effective number of neutrinos. Our results show that, despite the precision of the available cosmological data, only very weak bounds can be obtained on neutrino statistics, disfavouring a more bosonic behaviour at less than 2 sigma.

Abstract: One of the targets of future cosmic microwave background (CMB) and baryon acoustic oscillation measurements is to improve the current accuracy in the neutrino sector and reach a much better sensitivity on extra dark radiation in the early Universe. In this paper, we study how these improvements can be translated into constraining power for well-motivated extensions of the standard model of elementary particles that involve axions thermalized before the quantum chromodynamics (QCD) phase transition by scatterings with gluons. Assuming a fiducial Lambda cold dark matter cosmological model, we simulate future data for Stage-IV CMB-like and Dark Energy Spectroscopic Instrument (DESI)-like surveys and analyse a mixed scenario of axion and neutrino hot dark matter. We further account also for the effects of these QCD axions on the light element abundances predicted by big bang nucleosynthesis. The most constraining forecasted limits on the hot relic masses are m(a) less than or similar to 0.92 eV and n-ary sumation m(nu) less than or similar to 0.12 eV at 95 per cent Confidence Level, showing that future cosmic observations can substantially improve the current bounds, supporting multimessenger analyses of axion, neutrino, and primordial light element properties.

Abstract: The combination charge conjugation-parity-time reversal (CPT) is a fundamental symmetry in our current understanding of nature. As such, testing CPT violation is a strongly motivated path to explore new physics. In this paper we study CPT violation in the neutrino sector, giving for the first time a bound, for a fundamental particle, in the CPT violating particle-antiparticle gravitational mass difference. We argue that cosmology is nowadays the only data sensitive to CPT violation for the neutrino-antineutrino mass splitting and we use the latest data release from Planck combined with the current baryonic-acoustic-oscillation measurement to perform a full cosmological analysis. To show the potential of the future experiments we also show the results for Euclid, a next generation large scale structure experiment.