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: Neutrinos are one of the most promising messengers for signals of new physics Beyond the Standard Model (BSM). On the theoretical side, their elusive nature, combined with their unknown mass mechanism, seems to indicate that the neutrino sector is indeed opening a window to new physics. On the experimental side, several long-standing anomalies have been reported in the past decades, providing a strong motivation to thoroughly test the standard three-neutrino oscillation paradigm. In this Snowmass21 white paper, we explore the potential of current and future neutrino experiments to explore BSM effects on neutrino flavor during the next decade.

Abstract: One of the basic propositions of quantum field theory is Lorentz invariance. The spontaneous breaking of Lorentz symmetry at a high energy scale can be studied at low energy extensions like the Standard model in a model-independent way through effective field theory (EFT). The present and future Long-baseline neutrino experiments can give a scope to observe such a Planck-suppressed physics of Lorentz invariance violation (LIV). A proposed long baseline experiment, Protvino to ORCA (dubbed “P2O”) with a baseline of 2595 km, is expected to provide good sensitivities to unresolved issues, especially neutrino mass ordering. P2O can offer good statistics even with a moderate beam power and runtime, owing to the very large (similar to 6 Mt) detector volume at KM3NeT/ ORCA. Here we discuss in detail, how the individual LIV parameters affect neutrino oscillations at P2O and DUNE baselines at the level of probability and derive analytical expressions to understand interesting degeneracies and other features. We estimate increment Delta chi(2) sensitivities to the LIV parameters, analyzing their correlations among each other, and also with the standard oscillation parameters. We calculate these results for P2O alone and also carry out a combined analysis of P2O with DUNE. We point out crucial features in the sensitivity contours and explain them qualitatively with the help of the relevant probability expressions derived here. Finally we estimate constraints on the individual LIV parameters at 95% confidence level (C.L.) intervals stemming from the combined analysis of P2O and DUNE datasets, and highlight the improvement over the existing constraints. We also find out that the additional degeneracy induced by the LIV parameter a(ee) around -22 x 10(-23) GeV is lifted by the combined analysis at 95% C.L.