Abstract: Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library.

Abstract: We compute the leading-logarithmic correction to the neutrino mass matrix in the Standard Model effective field theory (SMEFT) to dimension seven. In the limit of negligible lepton and down-type quark Yukawa couplings, it receives contributions from the Weinberg dimension-five operator as well as from 11 dimension-six and five dimension-seven independent interactions. Two of the main implications we derive from this result are the following. First, we find dimension-seven operators which, despite violating lepton number, do not renormalize neutrino masses at one loop. And second, we demonstrate that the presence of dimension-six operators around the TeV scale can modify the Standard Model prediction by up to O(50%). Our result comprises also one step forward towards the renormalization of the SMEFT to order v(3)/Lambda(3).

Abstract: We study the capabilities of the DUNE near detector to probe deviations from unitarity of the leptonic mixing matrix, the 3+1 sterile formalism and Non-Standard Interactions affecting neutrino production and detection. We clarify the relation and possible mappings among the three formalisms at short-baseline experiments, and we add to current analyses in the literature the study of the nu(mu)-> nu(tau) appearance channel. We study in detail the impact of spectral uncertainties on the sensitivity to new physics using the DUNE near detector, which has been widely overlooked in the literature. Our analysis shows that this plays an important role on the results and, in particular, that it can lead to a strong reduction in the sensitivity to sterile neutrinos from nu(mu)-> nu(e) transitions, by more than two orders of magnitude. This stresses the importance of a joint experimental and theoretical effort to improve our understanding of neutrino nucleus cross sections, as well as hadron production uncertainties and beam focusing effects. Nevertheless, even with our conservative and more realistic implementation of systematic uncertainties, we find that an improvement over current bounds in the new physics frameworks considered is generally expected if spectral uncertainties are below the 5% level.

Abstract: In addition to spectacular signatures such as black hole superradiance and the rotation of CMB polarization, the plenitude of axions appearing in the string axiverse may have potentially dangerous implications. An example is the cosmological overproduction of relic axions and moduli by the misalignment mechanism, more pronounced in regions where the signals mentioned above may be observable, that is for large axion decay constant. In this work, we study the minimal requirements to soften this problem and show that the fundamental requirement is a long period of low-scale inflation. However, in this case, if the inflationary Hubble scale is lower than around O(100) eV, no relic DM axion is produced in the early Universe. Cosmological production of some axions may be activated, via the misalignment mechanism, if their potential minimum changes between inflation and today. As a particular example, we study in detail how the maximal-misalignment mechanism dilutes the effect of dangerous axions and allows the production of axion DM in a controlled way. In this case, the potential of the axion that realises the mechanism shifts by a factor increment theta = pi between the inflationary epoch and today, and the axion starts to oscillate from the top of its potential. We also show that axions with masses m(a) similar to O(1 – 100) H-0 realising the maximal-misalignment mechanism generically behave as dark energy with a decay constant that can take values well below the Planck scale, avoiding problems associated to super-Planckian scales. Finally, we briefly study the basic phenomenological implications of the mechanism and comment on the compatibility of this type of maximally-misaligned quintessence with the swampland criteria.

Abstract: In the context of lepton flavor universality violation (LFUV) studies, we fully derive a general tensor formalism to investigate the role that left- and right-handed neutrino new-physics (NP) terms may have in b -> c tau(nu) over bar (tau) transitions. We present, for several extensions of the Standard Model (SM), numerical results for the Lambda(b) -> Lambda(c)tau(nu) over bar (tau) semileptonic decay, which is expected to be measured with precision at the LHCb. This reaction can be a new source of experimental information that can help to confirm, or maybe rule out, LFUV presently seen in (B) over bar meson decays. The present study analyzes observables that can help in distinguishing between different NP scenarios that otherwise provide very similar results for the branching ratios, which are our currently best hints for LFUV. Since the tau lepton is very short-lived, we consider three subsequent tau-decay modes, two hadronic pi nu(tau) and rho nu(tau) and one leptonic mu(nu) over bar (mu)nu(tau), which have been previously studied for (B) over bar -> D(*) decays. Within the tensor formalism that we have developed in previous works, we re-obtain the expressions for the differential decay width written in terms of visible (experimentally accessible) variables of the massive particle created in the tau decay. There are seven different tau angular and spin asymmetries that are defined in this way and that can be extracted from experiment. Those asymmetries provide observables that can help in constraining possible SM extensions.