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: 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)).

Abstract: Leptoquarks are ubiquitous in several extensions of the Standard Model and seem to be able to accommodate the universality-violation-driven B-meson-decay anomalies and the (g-2)(mu) discrepancy interpreted as deviations from the Standard Model predictions. In addition, the search for lepton-flavour violation in the charged sector is, at present, a major research program that could also be facilitated by the dynamics generated by leptoquarks. In this article, we consider a rather wide framework of both scalar and vector leptoquarks as the generators of lepton-flavour violation in processes involving the tau lepton. We single out its couplings to leptoquarks, thus breaking universality in the lepton sector, and we integrate out leptoquarks at tree level, generating the corresponding dimension-6 operators of the Standard Model Effective Field Theory. In ref. [1] we obtained model-independent bounds on the Wilson coefficients of those operators contributing to lepton-flavour-violating hadron tau decays and l-tau conversion in nuclei, with l = e, mu. Hence, we use those results to translate the bounds into the couplings of leptoquarks to the Standard Model fermions.

Abstract: Dark-matter particles that annihilate or decay can undergo complex sequences of processes, including strong and electromagnetic radiation, hadronisation, and hadron de-cays, before particles that are stable on astrophysical time scales are produced. Antiprotons produced in this way may leave footprints in experiments such as AMS-02. Several groups have reported an excess of events in the antiproton flux in the rigidity range of 10-20 GV. However, the theoretical modeling of baryon production is not straightforward and relies in part on phenomenological models in Monte Carlo event generators. In this work, we assess the impact of QCD uncertainties on the spectra of antiprotons from dark-matter annihila-tion. As a proof-of-principle, we show that for a two-parameter model that depends only on the thermally-averaged annihilation cross section ((o -v)) and the dark-matter mass (Mx), QCD uncertainties can affect the best-fit mass by up to ti 14% (with large uncertainties for large DM masses), depending on the choice of Mx and the annihilation channel (bb over bar or W+W-), and (o -v) by up to ti 10%. For comparison, changes to the underlying diffusion parameters are found to be within 1%-5%, and the results are also quite resilient to the choice of cosmic-ray propagation model. These findings indicate that QCD uncertainties need to be included in future DM analyses. To facilitate full-fledged analyses, we provide the spectra in tabulated form including QCD uncertainties and code snippets to perform mass interpolations and quick DM fits. The code can be found in this GitHub [1] repository.

Abstract: NEXT is a new experiment to search for neutrinoless double beta decay using a 100 kg radio-pure high-pressure gaseous xenon TPC. The detector requires excellent energy resolution, which can be achieved in a Xe TPC with electroluminescence readout. Hamamatsu R8520-06SEL photomultipliers are good candidates for the scintillation readout. The performance of this photomultiplier, used as VUV photosensor in a gas proportional scintillation counter, was investigated. Initial results for the detection of primary and secondary scintillation produced as a result of the interaction of 5.9 keV X-rays in gaseous xenon, at room temperature and at pressures up to 3 bar, are presented. An energy resolution of 8.0% was obtained for secondary scintillation produced by 5.9 keV X-rays. No significant variation of the primary scintillation was observed for different pressures (1, 2 and 3 bar) and for electric fields up to 0.8 V cm(-1) torr(-1) in the drift region, demonstrating negligible recombination luminescence. A primary scintillation yield of 81 +/- 7 photons was obtained for 5.9 keV X-rays, corresponding to a mean energy of 72 +/- 6 eV to produce a primary scintillation photon in xenon.