|
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.
|
|
|
Figueroa, D. G., Florio, A., & Torrenti, F. (2024). Present and future of Cosmo Lattice. Rep. Prog. Phys., 87(9), 094901–20pp.
Abstract: We discuss the present state and planned updates of Cosmo Lattice, a cutting-edge code for lattice simulations of non-linear dynamics of scalar-gauge field theories in an expanding background. We first review the current capabilities of the code, including the simulation of interacting singlet scalars and of Abelian and non-Abelian scalar-gauge theories. We also comment on new features recently implemented, such as the simulation of gravitational waves from scalar and gauge fields. Secondly, we discuss new extensions of C osmo L attice that we plan to release publicly. We comment on new physics modules, which include axion-gauge interactions phi FF , non-minimal gravitational couplings phi R-2 , creation and evolution of cosmic-defect networks, and magnetohydrodynamics. We also discuss new technical features, including evolvers for non-canonical interactions, arbitrary initial conditions, simulations in 2+1 dimensions, and higher-accuracy spatial derivatives.
|
|
|
Easther, R., Price, L. C., & Rasero, J. (2014). Inflating an inhomogeneous universe. J. Cosmol. Astropart. Phys., 08(8), 041–16pp.
Abstract: While cosmological inflation can erase primordial inhomogeneities, it is possible that inflation may not begin in a significantly inhomogeneous universe. This issue is particularly pressing in multifield scenarios, where even the homogeneous dynamics may depend sensitively on the initial configuration. This paper presents an initial survey of the onset of inflation in multifield models, via qualitative lattice-based simulations that do not include local gravitational backreaction. Using hybrid inflation as a test model, our results suggest that small subhorizon inhomogeneities do play a key role in determining whether inflation begins in multifield scenarios. Interestingly, some configurations which do not inflate in the homogeneous limit “succeed” after inhomogeneity is included, while other initial configurations which inflate in the homogeneous limit “fail” when inhomogeneity is added.
|
|
|
Domcke, V., Garcia-Cely, C., Lee, S. M., & Rodd, N. L. (2024). Symmetries and selection rules: optimising axion haloscopes for Gravitational Wave searches. J. High Energy Phys., 03(3), 128–51pp.
Abstract: In the presence of electromagnetic fields, both axions and gravitational waves (GWs) induce oscillating magnetic fields: a potentially detectable fingerprint of their presence. We demonstrate that the response is largely dictated by the symmetries of the instruments used to search for it. Focussing on low mass axion haloscopes, we derive selection rules that determine the parametric sensitivity of different detector geometries to axions and GWs, and which further reveal how to optimise the experimental geometry to maximise both signals. The formalism allows us to forecast the optimal sensitivity to GWs in the range of 100 kHz to 100 MHz for instruments such as ABRACADABRA, BASE, ADMX SLIC, SHAFT, WISPLC, and DMRadio.
|
|
|
Di Bari, P., Ludl, P. O., & Palomares-Ruiz, S. (2016). Unifying leptogenesis, dark matter and high-energy neutrinos with right-handed neutrino mixing via Higgs portal. J. Cosmol. Astropart. Phys., 11(11), 044–41pp.
Abstract: We revisit a model in which neutrino masses and mixing are described by a two right-handed (RH) neutrino seesaw scenario, implying a strictly hierarchical light neutrino spectrum. A third decoupled RH neutrino, N-DM with mass M-DM, plays the role of cold dark matter (DM) and is produced by the mixing with a source RH neutrino, Ns with mass M-S, induced by Higgs portal interactions. The same interactions are also responsible for N-DM decays. We discuss in detail the constraints coming from DM abundance and stability conditions showing that in the hierarchical case, for M-DM >> M-S, there is an allowed window on M-DM values necessarily implying a contribution, from DM decays, to the high-energy neutrino flux recently detected by IceCube. We also show how the model can explain the matter-antimatter asymmetry of the Universe via leptogenesis in the quasi-degenerate limit. In this case, the DM mass should be within the range 300 GeV less than or similar to M-S < M-DM < 10PeV. We discuss the specific properties of this high-energy neutrino flux and show the predicted event spectrum for two exemplary cases. Although DM decays, with a relatively hard spectrum, cannot account for all the IceCube high-energy data, we illustrate how this extra source of high-energy neutrinos could reasonably explain some potential features in the observed spectrum. In this way, this represents a unified scenario for leptogenesis and DM that could be tested during the next years with more high-energy neutrino events.
|
|