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.

Abstract: Light primordial black holes (PBHs) with masses smaller than 10(9) g (10(-24) M-circle dot) evaporate before the onset of Big-Bang nucleosynthesis, rendering their detection rather challenging. If efficiently produced, they may have dominated the universe energy density. We study how such an early matter-dominated era can be probed successfully using gravitational waves (GW) emitted by local and global cosmic strings. While previous studies showed that a matter era generates a single-step suppression of the GW spectrum, we instead find a double-step suppression for local-string GW whose spectral shape provides information on the duration of the matter era. The presence of the two steps in the GW spectrum originates from GW being produced through two events separated in time: loop formation and loop decay, taking place either before or after the matter era. The second step – called the knee – is a novel feature which is universal to any early matter-dominated era and is not only specific to PBHs. Detecting GWs from cosmic strings with LISA, ET, or BBO would set constraints on PBHs with masses between 10(6) and 10(9) g for local strings with tension G μ= 10(-11), and PBHs masses between 10(4) and 10(9) g for global strings with symmetry-breaking scale eta = 10(15) GeV. Effects from the spin of PBHs are discussed.

Abstract: We propose a minimal gauged U(1)(B-L) extension of the minimal supersymmetric Standard Model (MSSM) which resolves the cosmological moduli problem via thermal inflation, and realizes late-time Affleck-Dine leptogensis so as to generate the right amount of baryon asymmetry at the end of thermal inflation. The present relic density of dark matter can be explained by sneutrinos, MSSM neutralinos, axinos, or axions. Cosmic strings from U(1)(B-L) breaking are very thick, and so the expected stochastic gravitational wave background from cosmic string loops has a spectrum different from the one in the conventional Abelian-Higgs model, as would be distinguishable at least at LISA and DECIGO. The characteristic spectrum is due to a flat potential, and may be regarded as a hint of supersymmetry. Combined with the resolution of moduli problem, the expected signal of gravitational waves constrains the U(1)(B-L) breaking scale to be O(10(12-13)) GeV. Interestingly, our model provides a natural possibility for explaining the observed ultra-high-energy cosmic rays thanks to the fact that the core width of strings in our scenario is very large, allowing a large enhancement of particle emissions from the cusps of string loops. Condensation of LHu flat-direction inside of string cores arises inevitably and can also be the main source of the ultra-high-energy cosmic rays accompanied by ultra-high-energy lightest supersymmetric particles.

Abstract: A brief review of different quantum Monte Carlo simulations of small (p-H-2)(N) clusters is presented. The clusters are viewed as a set of N structureless p-H-2 molecules, interacting via an isotropic pairwise potential. Properties as superfluidity, magic numbers, radial structure, excitation spectra, and abundance production of (p-H-2)(N) clusters are discussed and, whenever possible, a comparison with He-4(N) droplets is presented. All together, the simulations indicate that temperature has a paradoxical effect of the properties of (p-H-2)(N) clusters, as they are solid-like at high T and liquid-like at low T, due to quantum delocalization at the lowest temperature.