Abstract: We present a detailed study of the dynamics associated with the ghost sector of quenched QCD in the Landau gauge, where the relevant dynamical equations are supplemented with key inputs originating from large-volume lattice simulations. In particular, we solve the coupled system of Schwinger-Dyson equations that governs the evolution of the ghost dressing function and the ghost-gluon vertex, using as input for the gluon propagator lattice data that have been cured from volume and discretization artifacts. In addition, we explore the soft gluon limit of the same system, employing recent lattice data for the three-gluon vertex that enters in one of the diagrams defining the Schwinger-Dyson equation of the ghost-gluon vertex. The results obtained from the numerical treatment of these equations are in excellent agreement with lattice data for the ghost dressing function, once the latter have undergone the appropriate scale-setting and artifact elimination refinements. Moreover, the coincidence observed between the ghost-gluon vertex in general kinematics and in the soft gluon limit reveals an outstanding consistency of physical concepts and computational schemes.

Abstract: We present a detailed description of the electromagnetic filter for the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Starting with an initial estimate for the orbital magnetic moment, the higher-order drift process of E x B is configured to balance the gradient-B drift motion of the electron in such a way as to guide the trajectory into the standing voltage potential along the mid-plane of the filter. As a function of drift distance along the length of the filter, the filter zooms in with exponentially increasing precision on the transverse velocity component of the electron kinetic energy. This yields a linear dimension for the total filter length that is exceptionally compact compared to previous techniques for electromagnetic filtering. The parallel velocity component of the electron kinetic energy oscillates in an electrostatic harmonic trap as the electron drifts along the length of the filter. An analysis of the phase-space volume conservation validates the expected behavior of the filter from the adiabatic invariance of the orbital magnetic moment and energy conservation following Liouville's theorem for Hamiltonian systems. (C) 2019 Elsevier B.V. All rights reserved.

Abstract: We present a detailed analysis of the kinetic and mass terms associated with the Landau gauge gluon propagator in the presence of dynamical quarks, and a comprehensive dynamical study of certain special kinematic limits of the three-gluon vertex. Our approach capitalizes on results from recent lattice simulations with (2+1) domain wall fermions, a novel nonlinear treatment of the gluon mass equation, and the nonperturbative reconstruction of the longitudinal three-gluon vertex from its fundamental Slavnov-Taylor identities. Particular emphasis is placed on the persistence of the suppression displayed by certain combinations of the vertex form factors at intermediate and low momenta, already known from numerous pure Yang-Mills studies. One of our central findings is that the inclusion of dynamical quarks moderates the intensity of this phenomenon only mildly, leaving the asymptotic low-momentum behavior unaltered, but displaces the characteristic “zero crossing” deeper into the infrared region. In addition, the effect of the three-gluon vertex is explored at the level of the effective gauge coupling, whose size is considerably reduced with respect to its counterpart obtained from the ghost-gluon vertex. The main upshot of the above considerations is the further confirmation of the tightly interwoven dynamics between the two- and three-point sectors of QCD.

Abstract: We present a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, beta-decay, and neutrinoless double-beta-decay (0 nu beta beta) experiments. We directly constrain the full 3 + 1 active-sterile mixing matrix elements vertical bar U-alpha 4 vertical bar(2) , with alpha is an element of (e,mu,tau), and the mass-squared splitting Delta m(41)(2) (math) m(4)(2) – m(1)(2). We find that results for a 3 + 1 case differ from previously studied 1 + 1 scenarios where the sterile is coupled to only one of the neutrinos, which is largely explained by parameter space volume effects. Limits on the mass splitting and the mixing matrix elements are currently dominated by the cosmological datasets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below vertical bar U-alpha 4 vertical bar(2) less than or similar to 10(-3) . Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos arc in serious tension with these bounds, irrespective of prior assumptions. We also translate the bounds from the cosmological analysis into constraints on the parameters probed by laboratory searches, such as m(beta) or m(beta)(beta), the effective mass parameters probed by beta-decay and 0 nu beta beta searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of m(beta) < 0.09 eV and m(beta)(beta )< 0.07 eV at 95% CL, more stringent than the limits from current laboratory experiments.

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.