Search

We analyze how neutrino oscillation and coherent elastic neutrino-nucleus scattering data impact the global SMEFT fit. We first review the mapping between the SMEFT parameters and the so-called NSI framework, commonly considered in the neutrino literature. We also present a detailed discussion of how the measurements for the normalization of neutrino fluxes and cross sections, that will also be affected by the new physics, indirectly impact the measured oscillation probabilities. We then analyze two well-motivated simplified scenarios. Firstly, we study a lepton flavour conserving case, usually assumed in global SMEFT analyses, showing the complementarity of neutrino oscillation and CE$\nu$NS experiments with other low-energy observables. We find that the inclusion of neutrino data allows to constrain previously unbounded SMEFT operators involving the tau flavour and confirm the improvement of the constraint on a combination of Wilson coefficients previously identified. Moreover, we find that neutrino oscillation constraints on NSI are improved when embedded in the global SMEFT framework. Secondly, we study a lepton flavour violating scenario and find that neutrino data also improves over previously derived global constraints thanks to its sensitivity to new combinations of Wilson coefficients., Comment: 25 pages + Appendix

Cosmological neutrino mass bounds are becoming increasingly stringent. The latest limit within $\Lambda$CDM from Planck 2018+ACT lensing+DESI is $\sum m_\nu < 0.072\,{\rm eV}$ at 95\% CL, very close to the minimum possible sum of neutrino masses ($\sum m_\nu > 0.06\,{\rm eV}$), hinting at vanishing or even ``negative'' cosmological neutrino masses. In this context, it is urgent to carefully evaluate the origin of these cosmological constraints. In this paper, we investigate the robustness of these results in three ways: i) we check the role of potential anomalies in Planck CMB and DESI BAO data; ii) we compare the results for frequentist and Bayesian techniques, as very close to physical boundaries subtleties in the derivation and interpretation of constraints can arise; iii) we investigate how deviations from $\Lambda$CDM, potentially alleviating these anomalies, can alter the constraints. From a profile likelihood analysis, we derive constraints in agreement at the $\sim 10\%$ level with Bayesian posteriors. We find that the weak preference for negative neutrino masses is mostly present for Planck 18 data, affected by the well-known `lensing anomaly'. It disappears when the new Planck 2020 HiLLiPoP is used, leading to significantly weaker constraints. Additionally, the pull towards negative masses in DESI data stems from the $z=0.7$ bin, which contains a BAO measurement in $\sim 3\sigma$ tension with Planck expectations. Without this bin, and in combination with HiLLiPoP, the bound relaxes to $\sum m_\nu < 0.11\,{\rm eV}$ at 95\% CL. The recent preference for dynamical dark energy alleviates this tension and further weakens the bound. As we are at the dawn of a neutrino mass discovery from cosmology, it will be very exciting to see if this trend is confirmed by future data., Comment: 17 pages + Appendices. Matches version accepted for publication in PRD

Grain boundaries (GBs) in polycrystalline materials are powerful sinks for irradiation defects. While standard theories assume that the sink efficiency of a grain boundary is defined solely by its character before irradiation, recent evidence conclusively shows that the irradiation sink efficiency is a highly dynamic property controlled by the intrinsic metastability of GBs under far-from-equilibrium irradiation conditions. In this paper, we reveal that the denuded (i.e., defect-free) zone, typically the signature of a strong sink, can collapse as irradiation damage accumulates. We propose a radiation damage evolution model that captures this behavior based on the emergence of a series of irradiation defect-enabled metastable GB microstate changes that dynamically alter the ability of the GB to absorb further damage. We show that these microstate changes control further defect absorption and give rise to the formation of a defect network that manifests itself as a net Nye-tensor signal detectable via lattice curvature experiments.

Dislocation mobility, which dictates the response of dislocations to an applied stress, is a fundamental property of crystalline materials that governs the evolution of plastic deformation. Traditional approaches for deriving mobility laws rely on phenomenological models of the underlying physics, whose free parameters are in turn fitted to a small number of intuition-driven atomic scale simulations under varying conditions of temperature and stress. This tedious and time-consuming approach becomes particularly cumbersome for materials with complex dependencies on stress, temperature, and local environment, such as body-centered cubic crystals (BCC) metals and alloys. In this paper, we present a novel, uncertainty quantification-driven active learning paradigm for learning dislocation mobility laws from automated high-throughput large-scale molecular dynamics simulations, using Graph Neural Networks (GNN) with a physics-informed architecture. We demonstrate that this Physics-informed Graph Neural Network (PI-GNN) framework captures the underlying physics more accurately compared to existing phenomenological mobility laws in BCC metals.

We perform a global analysis of the bounds from charged lepton flavour violating observables to new physics. We parametrize generic new physics through the Effective Field Theory formalism and perform global fits beyond the common one-operator-at-a-time analyses to investigate how much present data is able to constrain the full parameter space. We particularly focus on leptonic and semileptonic operators with light quarks, identifying unbounded flat directions, detailing how many are present and which operators are involved. The analysis is performed in the general LEFT formalism, which contains all possible low-energy effective operators relevant for lepton flavour violation, as well as in more restricted scenarios, when operators come from a SMEFT completion. We find that flat directions play no role in the fully leptonic four-fermion operators. Conversely, they significantly hinder the ability to derive global bounds on semileptonic operators, with several flat or at least very poorly constrained directions preventing to fully constrain the parameter space. These results are particularly affected by the proper inclusion of uncertainties in the parameters describing $\mu-e$ conversion, which decrease the number of well-constrained directions in operator space when treated as nuisance parameters in the fit. While present data is able to provide global constraints on all operators only in the more restricted scenarios we investigated, very strong correlations among the parameters must exist to avoid conflict with the different observables. We provide correlation matrices approximating our full results as a useful tool to compare present data with particular UV completions., Comment: 28 pages + Appendices, 8 figures. Matches published version in EPJC. Global constraints and correlations available at https://github.com/dnaredo/cLFV_GlobalBounds.git

In this work we develop a stochastic algorithm to integrate the Cahn-Hilliard equations. The algorithm is based on Gillespie's stochastic simulation algorithm, also known as kinetic Monte Carlo. The deterministic integration of the phase field equations leads to the closest minimum of free energy and does not overcome free energy barriers. However, in the nucleation and growth regime of the phase diagram, free energy barriers need to be thermally overcome for the system to phase separate to reach equilibrium. We show in this work that the proposed stochastic integration algorithm indeed allows the system to overcome free energy barriers. We discuss the results in terms of fluctuation distributions, grid sizes and efficiency.

High-entropy alloys have drawn recent interest for their promising mechanical properties and irradiation resistance. Various properties, namely transport properties, are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for impurity concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy but clear systematic patterns for FeAl., Comment: 13 pages, 6 figures

Dislocations are the primary carriers of plasticity in metallic material. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate for an edge dislocation to overcome the elastic interaction with a prismatic loop in tungsten. Beyond the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius behavior. We have computed the activation entropy analytically along the minimum energy path following Schoeck's methods [1], which capture the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI) [2], another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck's formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck's and PAFI results compared to molecular dynamics.

We present an updated and improved global fit analysis of current flavor and electroweak precision observables to derive bounds on unitarity deviations of the leptonic mixing matrix and on the mixing of heavy neutrinos with the active flavours. This new analysis is motivated by new and updated experimental results on key observables such as $V_{ud}$, the invisible decay width of the $Z$ boson and the $W$ boson mass. It also improves upon previous studies by considering the full correlations among the different observables and explicitly calibrating the test statistic, which may present significant deviations from a $\chi^2$ distribution. The results are provided for three different Type-I seesaw scenarios: the minimal scenario with only two additional right-handed neutrinos, the next to minimal one with three extra neutrinos, and the most general one with an arbitrary number of heavy neutrinos that we parametrize via a generic deviation from a unitary leptonic mixing matrix. Additionally, we also analyze the case of generic deviations from unitarity of the leptonic mixing matrix, not necessarily induced by the presence of additional neutrinos. This last case relaxes some correlations among the parameters and is able to provide a better fit to the data. Nevertheless, inducing only leptonic unitarity deviations avoiding both the correlations implied by the right-handed neutrino extension as well as more strongly constrained operators is challenging and would imply significantly more complex UV completions., Comment: 27 pages + appendices, 7 figures, 7 tables

We construct a contact interaction model for the $\eta$ and $\eta'$ mesons in the SDE-BSE approach to QCD and compute several static properties of these mesons and their transition form factors. We find that this model gives an excellent description of the $\eta$ and $\eta'$ static properties, namely their masses, decay width and decay constants. However, a contact interaction disagrees with experimental data for $Q^2$ greater than 2 GeV$^2$, and produces transition form factors in conflict with perturbative QCD prediction. This is not surpring and the reasons for this are explained, Comment: 10 pages, 4 figures, submitted to PRD

We investigate the relation between essential divisors and F-blowups, in particular, address the problem whether all essential divisors appear on the $e$-th F-blowup for large enough $e$. Focusing on the case of normal affine toric varieties, we establish a simple sufficient condition for a divisor over the given toric variety to appear on the normalized limit F-blowup as a prime divisor. As a corollary, we show that if a normal toric variety has a crepant resolution, then the above problem has a positive answer, provided that we use the notion of essential divisors in the sense of Bouvier and Gonzalez-Sprinberg. We also provide an example of toric threefold singularities for which a non-essential divisor appears on an F-blowup., Comment: 24 pages. v2: minor revision, v3: fixed errors that was pointed out by a referee and based on confusion about several notions of essential divisors, v4: modified the definition of ciritical arrows, added some descriptions to help understanding of the discussion

The existence of right-handed neutrinos, or heavy neutral leptons (HNLs), is strongly motivated by the observation of neutrino masses and mixing. The mass of these new particles could lie below the electroweak scale, making them accessible to low-energy laboratory experiments. Additional new physics at high energies can mediate new interactions between the Standard Model particles and HNLs, and is most conveniently parametrized by the neutrino Standard Model Effective Field Theory, or $\nu$SMEFT for short. In this work, we consider the dimension six $\nu$SMEFT operators involving one HNL field in the mass range of $\mathcal{O}(1)$ MeV $

We propose a systematic method to construct crystal-based molecular structures often needed as input for computational chemistry studies. These structures include crystal ``slabs" with periodic boundary conditions (PBCs) and non-periodic solids such as Wulff structures. We also introduce a method to build crystal slabs with orthogonal PBC vectors. These methods are integrated into our code, Los Alamos Crystal Cut (LCC), which is open source and thus fully available to the community. Examples showing the use of these methods are given throughout the manuscript.

Low-scale seesaw variants protected by lepton number symmetry provide a natural explanation of the smallness of neutrino masses but, unlike their higher-scale counterparts, with potentially testable phenomenology. The approximate lepton number symmetry arranges the heavy neutrinos in pseudo-Dirac pairs, which might be accessible at collider or even beam dump experiments if their mass is low enough and their mixing with the active neutrinos sufficiently large. Despite their pseudo-Dirac nature, their small mass splittings may lead to oscillations that prevent the cancellation of their potential lepton-number-violating signals. Interestingly, these small splittings may also resonantly enhance the production of a lepton number asymmetry for low-scale leptogenesis scenarios or, for extremely degenerate states, lead to an asymmetry large enough to resonantly produce a keV sterile neutrino dark matter candidate with the correct relic abundance via the Shi-Fuller mechanism. In this work we explore the parameter space of the different low-scale seesaw mechanisms and study the size of these splittings, given their important and interesting phenomenological consequences. While all low-scale seesaw variants share the same dimension 5 and 6 operators when integrating out the heavy states, we point out that the mass splitting of the pseudo-Dirac pairs are very different in different realizations such as the inverse or linear seesaw. This different phenomenology could offer a way to discriminate between low-scale seesaw realizations., Comment: 27 pages, 6 figures. Matches published version in JHEP

The Inverse Seesaw mechanism is dynamically realised within the Minimal Lepton Flavour Violation context. Lepton number, whose breaking is spontaneously realised, is generalised to a global Abelian factor of the whole flavour symmetry, that also plays the role of the Peccei-Quinn symmetry. The associated Goldstone boson is a Majoraxion that solves the Strong CP problem and represents a Dark Matter candidate. Three distinct scenarios are identified in terms of flavour symmetry and transformation properties of the exotic neutral leptons that enrich the Standard Model spectrum. The associated phenomenology is studied, focusing on the deviations from unitarity of the PMNS mixing matrix. The strongest constraints arise from the determination of the number of active neutrinos through the invisible width of the $Z$, the comparison of the measured $W$ boson mass with its prediction in terms of the Fermi constant from muon decay, and the null searches for the radiative rare muon decay and $\mu\to e$ conversion in nuclei. The heavy neutral leptons may have masses of a few TeV, leaving open the possibility for a direct detection at future colliders. The impact of the recent measurement of the $W$ mass at the CDF II detector has also been considered, which, in one of the scenarios, points to a sharp prediction for the masses of the heavy neutral leptons at about $2-3$ TeV., Comment: 36 pages, 9 figures, 6 tables

We point out that right-handed neutrinos can resolve the tension between the latest CDF II measurement of the $W$-boson mass, $M_W$, and the standard model. Integrating out the new states yields a single $d=6$ operator, which induces a deviation from unitarity in the PMNS matrix. This alters the extraction of the Fermi constant from muon decay and increases the prediction for $M_W$, in line with the CDF II result. Non-unitarity of the PMNS matrix would also affect beta, meson, and tau decays. We find that the CDF II value for $M_W$ can be explained without conflicting with lepton flavour universality constraints or the invisible decay width of the $Z$. However, the so-called Cabibbo angle anomaly is worsened if right-handed neutrinos are the origin of the $d=6$ operator. The situation improves if the operator coefficient is left unconstrained, implying additional sources of new physics, but a common explanation of both anomalies is in tension with universality bounds., Comment: 5 pages, 1 figure, REVTeX4

We derive a simple analytical line integral expression for the relaxation volume tensor of an arbitrary interconnected dislocation network. This quantity determines the magnitude of dislocation contribution to the dimensional changes and volumetric swelling of a material, and highlights the fundamental dual role of dislocations as sources of internal strain as well as carriers of plastic deformation. To illustrate applications of the method, we compute the relaxation volume of a stacking fault tetrahedron, a defect commonly occurring in fcc metals; the volume of an unusual tetrahedral configuration formed by the (a/2)<111> and a<001> dislocations in a bcc metal; and estimate the relative contribution of extended dislocations to the volume relaxation of heavily irradiated tungsten., Comment: 15 pages, 6 figures

Most of the studies related to human mobility are focused on intra-country mobility. However, there are many scenarios (e.g., spreading diseases, migration) in which timely data on international commuters are vital. Mobile phones represent a unique opportunity to monitor international mobility flows in a timely manner and with proper spatial aggregation. This work proposes using roaming data generated by mobile phones to model incoming and outgoing international mobility. We use the gravity and radiation models to capture mobility flows before and during the introduction of non-pharmaceutical interventions. However, traditional models have some limitations: for instance, mobility restrictions are not explicitly captured and may play a crucial role. To overtake such limitations, we propose the COVID Gravity Model (CGM), namely an extension of the traditional gravity model that is tailored for the pandemic scenario. This proposed approach overtakes, in terms of accuracy, the traditional models by 126.9% for incoming mobility and by 63.9% when modeling outgoing mobility flows.

We show that the normalization of the Nash blow-up of order n of the toric surface singularity An can be factorized by the minimal resolution of An. The result is obtained using the combinatorial description of these objects.

We consider the inverse Seesaw scenario for neutrino masses with the approximate Lepton number symmetry broken dynamically by a scalar with Lepton number two. We show that the Majoron associated to the spontaneous symmetry breaking can alleviate the Hubble tension through its contribution to $\Delta N_\text{eff}$ and late decays to neutrinos. Among the additional fermionic states required for realizing the inverse Seesaw mechanism, sterile neutrinos at the keV-MeV scale can account for all the dark matter component of the Universe if produced via freeze-in from the decays of heavier degrees of freedom., Comment: 22 pages, 2 figures

In this work we examine the driving force for edge dislocations to climb in $\alpha$-Fe from atomistic and mesoscale viewpoints. We study the bias for the climb process depending on the dislocation orientation and the applied stress as coming from both the gradient of the chemical potential and the transport coefficients. Both terms are modified by the applied stress and therefore contribute to climb. Surprisingly, even though the vacancy migration barrier distribution is modified by the external stress as obtained by nudged-elastic band calculations, the mobilities resulting from a kinetic Monte Carlo model applied on the obtained energy landscape are indistinguishable, independently of the stress. Moreover, an object kinetic Monte Carlo (OKMC) model including the effect of the dislocation strain field to first order shows indeed a slight anisotropic component in the diffusion in more complex dislocation configurations. However, the OKMC results highlight the fact that the thermodynamic component is the dominant driving force. We conclude that in $\alpha$-Fe under thermal conditions, the main source of bias is given by the difference in vacancy chemical potentials, which is small enough to hinder the process for dynamic atomistic simulations.

During early 2020, the SARS-CoV-2 virus rapidly spread worldwide, forcing many governments to impose strict lockdown measures to tackle the pandemic. This significantly changed people's mobility and habits, subsequently impacting how they use telecommunication networks. In this paper, we investigate the effects of the COVID-19 emergency on a UK Mobile Network Operator (MNO). We quantify the changes in users' mobility and investigate how this impacted the cellular network usage and performance. Our analysis spans from the entire country to specific regions, and geodemographic area clusters. We also provide a detailed analysis for London. Our findings bring insights at different geo-temporal granularity on the status of the cellular network, from the decrease in data traffic volume in the cellular network and lower load on the radio network, counterposed to a surge in the conversational voice traffic volume., Comment: This paper will appear in ACM Internet Measurement Conference (IMC'20). Andra Lutu, Diego Perino, Marcelo Bagnulo, Enrique Frias-Martinez, and Javad Khangosstar. 2020. A Characterization of the COVID-19 Pandemic Impact on a Mobile Network Operator Traffic. In ACM Internet Measurement Conference (IMC'20), October 27-29, 2020, Virtual Event, USA. ACM, New York, NY, USA, 15 pages

The recent tension between local and early measurements of the Hubble constant can be explained in a particle physics context. A mechanism is presented where this tension is alleviated due to the presence of a Majoron, arising from the spontaneous breaking of Lepton Number. The lightness of the active neutrinos is consistently explained. Moreover, this mechanism is shown to be embeddable in the Minimal (Lepton) Flavour Violating context, providing a correct description of fermion masses and mixings, and protecting the flavour sector from large deviations from the Standard Model predictions. A QCD axion is also present to solve the Strong CP problem. The Lepton Number and the Peccei-Quinn symmetries naturally arise in the Minimal (Lepton) Flavour Violating setup and their spontaneous breaking is due to the presence of two extra scalar singlets. The Majoron phenomenology is also studied in detail. Decays of the heavy neutrinos and the invisible Higgs decay provide the strongest constraints in the model parameter space., Comment: 18 pages, 1 figure, 4 tables. v3: improved discussions, typos fixed

The simplest extension of the SM to account for the observed neutrino masses and mixings is the addition of at least two singlet fermions (or right-handed neutrinos). If their masses lie at or below the GeV scale, such new fermions would be produced in meson decays. Similarly, provided they are sufficiently heavy, their decay channels may involve mesons in the final state. Although the couplings between mesons and heavy neutrinos have been computed previously, significant discrepancies can be found in the literature. The aim of this paper is to clarify such discrepancies and provide consistent expressions for all relevant effective operators involving mesons with masses up to 2 GeV. Moreover, the effective Lagrangians obtained for both the Dirac and Majorana scenarios are made publicly available as FeynRules models so that fully differential event distributions can be easily simulated. As an application of our setup, we numerically compute the expected sensitivity of the DUNE near detector to these heavy neutral leptons., Comment: v5: Minor modifications. Published in EPJC. FeynRules models in ancillary files

The evolution of point defect concentrations under irradiation is controlled by their diffusion properties, and by their formation and elimination mechanisms. The latter include the mutual recombination of vacancies and interstitials, and the elimination of point defects at sinks. We show here that the modelling of this evolution by means of atomistic kinetic Monte Carlo (AKMC) simulations, necessarily using small system sizes, introduces strong space and time correlations between the vacancies and interstitials, which may strongly affect the recombination rate and the point defect concentrations. In such situations, standard rate theory models fail to predict the actual point defect concentrations. The effect is especially strong when the elimination of point defects occurs only by recombination, but can still be significant in the presence of sinks. We propose a new Correlated Pair Theory that fully takes into account the correlations between vacancy and interstitial pairs and predicts point defect concentrations in good agreement with AKMC simulations, even in very small systems. The Correlated Pair Theory can be used to modify the elimination rates in AKMC simulations to yield point defect concentrations as predicted by the standard rate theory, i.e. representative of large systems, even when using small simulation boxes., Comment: 13 pages, 7 figures

Mobile phone data has enabled the timely and fine-grained study human mobility. Call Detail Records, generated at call events, allow building descriptions of mobility at different resolutions and with different spatial, temporal and social granularity. Individual trajectories are the basis for long-term observation of mobility patterns and identify factors of human dynamics. Here we propose a systematic analysis to characterize mobility network flows and topology and assess their impact into individual traces. Discrete flow-based descriptors are used to classify and understand human mobility patterns at multiple scales. This framework is suitable to assess urban planning, optimize transportation, measure the impact of external events and conditions, monitor internal dynamics and profile users according to their movement patterns.