• ### $K \to \pi$ matrix elements of the chromomagnetic operator on the lattice(1712.09824)

March 12, 2018 hep-ph, hep-lat
We present the results of the first lattice QCD calculation of the $K \to \pi$ matrix elements of the chromomagnetic operator $O_{CM} = g\, \bar s\, \sigma_{\mu\nu} G_{\mu\nu} d$, which appears in the effective Hamiltonian describing $\Delta S = 1$ transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been determined non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing. Our result for the B-parameter of the chromomagnetic operator at the physical pion and kaon point is $B_{CMO}^{K \pi} = 0.273 ~ (70)$, while in the SU(3) chiral limit we obtain $B_{CMO} = 0.072 ~ (22)$. Our findings are significantly smaller than the model-dependent estimate $B_{CMO} \sim 1 - 4$, currently used in phenomenological analyses, and improve the uncertainty on this important phenomenological quantity.
• ### Origin of the high energy neutrino flux at IceCube(1703.10786)

Nov. 27, 2017 hep-ph, astro-ph.HE
We discuss the spectrum of the different components in the astrophysical neutrino flux reaching the Earth and the possible contribution of each component to the high-energy IceCube data. We show that the diffuse flux from cosmic ray interactions with gas in our galaxy implies just 2 events among the 54 event sample. We argue that the neutrino flux from cosmic ray interactions in the intergalactic (intracluster) space depends critically on the transport parameter $\delta$ describing the energy dependence in the diffusion coefficient of galactic cosmic rays. Our analysis motivates a E^{-2.1} neutrino spectrum with a drop at PeV energies that fits well the data, including the non-observation of the Glashow resonance at 6.3 PeV. We also show that a cosmic ray flux described by an unbroken power law may produce a neutrino flux with interesting spectral features (bumps and breaks) related to changes in the cosmic ray composition.
• ### The chromomagnetic operator on the lattice(1506.00361)

June 1, 2015 hep-ph, hep-lat
We present our study of the renormalization of the chromomagnetic operator,O(CM), which appears in the effective Hamiltonian describing Delta S = 1 transitions in and beyond the Standard Model. We have computed, perturbatively to one-loop, the relevant Green's functions with two (quark-quark) and three (quark-quark-gluon) external fields, at nonzero quark masses, using both the lattice and dimensional regularizations. The perturbative computation on the lattice is carried out using the maximally twisted-mass action for the fermions, while for the gluons we employed the Symanzik improved gauge action for different sets of values of the Symanzik coefficients. We have identified all the operators which can possibly mix with O(CM), including lower dimensional and non gauge invariant operators, and we have calculated those elements of the mixing matrix which are relevant for the renormalization of O(CM). We have also performed numerical lattice calculations to determine non-perturbatively the mixings of the chromomagnetic operator with lower dimensional operators, through proper renormalization conditions. For the first time the 1/a**2-divergent mixing of the chromomagnetic operator with the scalar density has been determined non-perturbatively with high precision. Moreover, the 1/a-divergent mixing with the pseudoscalar density, due to the breaking of parity within the twisted-mass regularization of QCD, has been calculated non-perturbatively and found to be smaller than its one-loop perturbative estimate. The QCD simulations have been carried out using the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2 + 1 + 1 dynamical quarks, which include in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to their physical values.
• ### Probing new physics scenarios in accelerator and reactor neutrino experiments(1411.5330)

May 11, 2015 hep-ph
We perform a detailed combined fit to the $\overline \nu_e \rightarrow \overline \nu_e$ disappearence data of the Daya Bay experiment and the appearance $\nu_{\mu} \rightarrow \nu_{e}$ and disappearance $\nu_{\mu} \rightarrow \nu_{\mu}$ data of the Tokai to Kamioka (T2K) one in the presence of two models of new physics affecting neutrino oscillations, namely a model where sterile neutrinos can propagate in a large compactified extra dimension and a model where non-standard interactions (NSI) affect the neutrino production and detection. We find that the Daya Bay $\oplus$ T2K data combination constrains the largest radius of the compactified extra dimensions to be $R\lesssim 0.17$ $\mu {\rm m}$ at 2$\sigma$ C.L. (for the inverted ordering of the neutrino mass spectrum) and the relevant NSI parameters in the range ${\mathcal O}(10^{-3})-{\mathcal O}(10^{-2})$, for particular choices of the charged parity violating phases.
• ### Perturbative and non-perturbative renormalization results of the Chromomagnetic Operator on the Lattice(1412.1718)

Dec. 4, 2014 hep-lat
The Chromomagnetic operator (CMO) mixes with a large number of operators under renormalization. We identify which operators can mix with the CMO, at the quantum level. Even in dimensional regularization (DR), which has the simplest mixing pattern, the CMO mixes with a total of 9 other operators, forming a basis of dimension-five, Lorentz scalar operators with the same flavor content as the CMO. Among them, there are also gauge noninvariant operators; these are BRST invariant and vanish by the equations of motion, as required by renormalization theory. On the other hand using a lattice regularization further operators with $d \leq 5$ will mix; choosing the lattice action in a manner as to preserve certain discrete symmetries, a minimul set of 3 additional operators (all with $d<5$) will appear. In order to compute all relevant mixing coefficients, we calculate the quark-antiquark (2-pt) and the quark-antiquark-gluon (3-pt) Green's functions of the CMO at nonzero quark masses. These calculations were performed in the continuum (dimensional regularization) and on the lattice using the maximally twisted mass fermion action and the Symanzik improved gluon action. In parallel, non-perturbative measurements of the $K-\pi$ matrix element are being performed in simulations with 4 dynamical ($N_f = 2+1+1$) twisted mass fermions and the Iwasaki improved gluon action.
• ### K -> pi matrix elements of the chromagnetic operator on the lattice(1412.1351)

Dec. 3, 2014 hep-lat
We present preliminary results of the first lattice QCD calculation of the K -> pi matrix elements of the chromomagnetic operator O_{CM}=g sbar sigma_{munu} G_{munu} d, which appears in the effective Hamiltonian describing Delta S=1 transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one-loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been computed non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with N_f=2+1+1 dynamical quarks at three values of the lattice spacing. Our preliminary result for the B-parameter of the chromomagnetic operator is B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~1-4 currently used in phenomenological analyses.
• ### The Daya Bay and T2K results on $\sin^2 2 \theta_{13}$ and Non-Standard Neutrino Interactions(1405.0416)

Aug. 9, 2014 hep-ph
We show that the relatively large best fit value of $\sin^2 2 \theta_{13} = 0.14 \, (0.17)$ measured in the T2K experiment for fixed values of i) the Dirac CP violation phase $\delta = 0$, and ii) the atmospheric neutrino mixing parameters $\theta_{23} = \pi/4$, $|\Delta m^2_{32}| = 2.4 \times 10^{-3} \; {\rm eV}^2$, can be reconciled with the Daya Bay result $\sin^2 2 \theta_{13} = 0.090 \pm 0.009$ if the effects of non-standard neutrino interactions (NSI) in the relevant $\bar \nu_e \to \bar \nu_e$ and $\nu_\mu \to \nu_e$ oscillation probabilities are taken into account.
• ### Constraining new physics scenarios in neutrino oscillations from Daya Bay data(1403.5507)

March 21, 2014 hep-ph
We perform for the first time a detailed fit to the $\bar \nu_e \to \bar \nu_e$ disappearance data of the Daya Bay experiment to constrain the parameter space of models where sterile neutrinos can propagate in a large compactified extra dimension (LED) and models where non-standard interactions affect the neutrino production and detection (NSI). We find that the compactification radius $R$ in LED scenarios can be constrained at the level of $0.57 \, \mu m$ for normal ordering and of $0.19\, \mu m$ for inverted ordering, at 2$\sigma$ confidence level. For the NSI model, reactor data put a strong upper bound on the parameter $\varepsilon_{ee}$ at the level of $\sim 10^{-3}$, whereas the main effect of $\varepsilon_{e\mu}$ and $\varepsilon_{e\tau}$ is a worsening of the determination of $\theta_{13}$.
• ### Two-zero Majorana textures in the light of the Planck results(1401.3207)

Jan. 14, 2014 hep-ph
The recent results of the Planck experiment put a stringent constraint on the sum of the light neutrino masses, m1+m2+m3 < 0.23 eV (95 % CL). On the other hand, two-zero Majorana mass matrix textures predict strong correlations among the atmospheric angle and the sum of the masses. We use the Planck result to show that, for the normal hierarchy case, the texture with vanishing (2,2) and (3,3) elements is ruled out at a high confidence level; in addition, we emphasize that a future measurement of the octant of the atmospheric mixing angle (or the one sigma determination of it based on recent fit to neutrino data) will put severe constraint on the possible structure of the Majorana mass matrix. The implication of the above mentioned correlations for neutrinoless double beta-decay are also discussed, for both normal and inverted orderings.
• ### The chromomagnetic operator on the lattice(1311.5057)

Nov. 20, 2013 hep-lat
We study matrix elements of the "chromomagnetic" operator on the lattice. This operator is contained in the strangeness-changing effective Hamiltonian which describes electroweak effects in the Standard Model and beyond. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with other operators of equal and lower dimensionality, including also non gauge invariant quantities; it is thus quite a challenge to extract from lattice simulations a clear signal for the hadronic matrix elements of this operator. We compute all relevant mixing coefficients to one loop in lattice perturbation theory; this necessitates calculating both 2-point (quark-antiquark) and 3-point (gluon-quark-antiquark) Green's functions at nonzero quark masses. We use the twisted mass lattice formulation, with Symanzik improved gluon action. For a comprehensive presentation of our results, along with detailed explanations and a more complete list of references, we refer to our forthcoming publication [1].
• The EUROnu project has studied three possible options for future, high intensity neutrino oscillation facilities in Europe. The first is a Super Beam, in which the neutrinos come from the decay of pions created by bombarding targets with a 4 MW proton beam from the CERN High Power Superconducting Proton Linac. The far detector for this facility is the 500 kt MEMPHYS water Cherenkov, located in the Fr\'ejus tunnel. The second facility is the Neutrino Factory, in which the neutrinos come from the decay of {\mu}+ and {\mu}- beams in a storage ring. The far detector in this case is a 100 kt Magnetised Iron Neutrino Detector at a baseline of 2000 km. The third option is a Beta Beam, in which the neutrinos come from the decay of beta emitting isotopes, in particular 6He and 18Ne, also stored in a ring. The far detector is also the MEMPHYS detector in the Fr\'ejus tunnel. EUROnu has undertaken conceptual designs of these facilities and studied the performance of the detectors. Based on this, it has determined the physics reach of each facility, in particular for the measurement of CP violation in the lepton sector, and estimated the cost of construction. These have demonstrated that the best facility to build is the Neutrino Factory. However, if a powerful proton driver is constructed for another purpose or if the MEMPHYS detector is built for astroparticle physics, the Super Beam also becomes very attractive.
• ### Fine-tuning and naturalness issues in the two-zero neutrino mass textures(1204.2706)

Oct. 17, 2012 hep-ph
In this paper we analyze the compatibility of two-zero neutrino Majorana textures with the recent experimental data. Differently from previous works, we use the experimental data to fix the values of the non-vanishing mass matrix entries and study in detail the correlations and degree of fine-tuning among them, which is also a measure of how naturally a given texture is able to describe all neutrino data. This information is then used to expand the textures in powers of the Cabibbo angle; extracting random O(1) coefficients, we show that only in few cases such textures reproduce the mixing parameters in their 3 sigma ranges.
• ### Revisiting the T2K data using different models for the neutrino-nucleus cross sections(1203.3335)

Sept. 19, 2012 hep-ph, hep-ex, nucl-th
We present a three-flavour fit to the recent \nu\mu --> \nu e and \nu\mu --> \nu\mu T2K oscillation data with different models for the neutrino-nucleus cross section. We show that, even for a limited statistics, the allowed regions and best fit points in the (\theta_{13},\delta_{CP}) and (\theta_{23},\Delta m^2_{atm}) planes are affected if, instead of using the Fermi Gas model to describe the quasielastic cross section, we employ a model including the multinucleon emission channel.
• ### EUROnu-WP6 2010 Report(1209.2825)

Sept. 13, 2012 hep-ph
This is a summary of the work done by the Working Package 6 (Physics) of the EU project "EUROnu" during the second year of activity of the project.
• ### D14 - A Common Origin of the Cabibbo Angle and the Lepton Mixing Angle theta^l_13(1204.0715)

May 18, 2012 hep-ph
It has been shown that the Cabibbo angle can be predicted in terms of group theoretical quantities, if the dihedral group D14 plays the role of a flavor symmetry. We extend a supersymmetric D14 model to the lepton sector and show that theta^nu_13 and the deviation of theta^nu_23 from maximal mixing in the neutrino sector originate, similar to the Cabibbo angle in the quark sector, from a mismatch of different subgroups of D14 and are of the size of the Cabibbo angle. The mixing angles in the charged lepton sector are small. Thus, the lepton mixing angle theta^l_13 is naturally in its experimentally preferred range and theta^l_23 within its 3 sigma range. The solar mixing angle is of order one and the charged lepton mass hierarchy is correctly reproduced. Leading order results are only slightly perturbed, if next-to-leading order corrections are taken into account.
• This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.
• ### Predicting leptonic CP violation in the light of Daya Bay result(1203.2535)

March 12, 2012 hep-ph, hep-ex
In the light of the recent Daya Bay result the reactor angle is about 9 degrees, we reconsider the model presented in arXiv:1005.3482 showing that, when all neutrino oscillation parameters are taken at their best fit values of Schwetz et al and the reactor angle to be the central value of Daya Bay, the predicted value of the CP phase is approximately 45 degrees.
• ### Proceedings of the first workshop on Flavor Symmetries and consequences in Accelerators and Cosmology (FLASY2011)(1201.5525)

Jan. 26, 2012 hep-ph
The main goals of the first "Workshop on FLAvor SYmmetries and consequences in accelerators and cosmology" (FLASY) was to summarize the theoretical status of flavor symmetries, bringing together young researchers in the field to stimulate discussions and new collaborations, with the aim of investigating possible new physics scenarios to be tested at the LHC, as well as in future neutrino, cosmology experiments and dark matter searches.
• The International Design Study for the Neutrino Factory (the IDS-NF) was established by the community at the ninth "International Workshop on Neutrino Factories, super-beams, and beta- beams" which was held in Okayama in August 2007. The IDS-NF mandate is to deliver the Reference Design Report (RDR) for the facility on the timescale of 2012/13. In addition, the mandate for the study [3] requires an Interim Design Report to be delivered midway through the project as a step on the way to the RDR. This document, the IDR, has two functions: it marks the point in the IDS-NF at which the emphasis turns to the engineering studies required to deliver the RDR and it documents baseline concepts for the accelerator complex, the neutrino detectors, and the instrumentation systems. The IDS-NF is, in essence, a site-independent study. Example sites, CERN, FNAL, and RAL, have been identified to allow site-specific issues to be addressed in the cost analysis that will be presented in the RDR. The choice of example sites should not be interpreted as implying a preferred choice of site for the facility.
• ### Constraining neutrinoless double beta decay(1111.5614)

Nov. 23, 2011 hep-ph, hep-ex, nucl-ex
A class of discrete flavor-symmetry-based models predicts constrained neutrino mass matrix schemes that lead to specific neutrino mass sum-rules (MSR). We show how these theories may constrain the absolute scale of neutrino mass, leading in most of the cases to a lower bound on the neutrinoless double beta decay effective amplitude.
• ### Stability of dark matter from the D4xZ2 flavor group(1104.0178)

April 14, 2011 hep-ph, astro-ph.CO, astro-ph.HE
We study a model based on the dihedral group D4 in which the dark matter is stabilized by the interplay between a remnant Z2 symmetry, of the same spontaneously broken non-abelian group, and an auxiliary Z2 introduced to eliminate unwanted couplings in the scalar potential. In the lepton sector the model is compatible with normal hierarchy only and predicts a vanishing reactor mixing angle. Since m1=0, we also have a simple prediction for the effective mass in terms of the solar angle. There also exists a large portion of the model parameter space where the upper bounds on lepton flavor violating processes are not violated. We incorporate quarks in the same scheme finding that a description of the CKM mixing matrix is possible and that semileptonic K and D decays mediated by flavour changing neutral currents are under control.
• ### The $\tau$-contamination of the golden muon sample at the Neutrino Factory(1005.2275)

Feb. 10, 2011 hep-ph
We study the contribution of $\nu_e \to \nu_\tau \to \tau \to \mu$ transitions to the wrong-sign muon sample of the {\em golden} channel of the Neutrino Factory. Muons from tau decays are not really a background, since they contain information from the oscillation signal, and represent a small fraction of the sample. However, if not properly handled they introduce serious systematic error, in particular if the detector/analysis are sensitive to muons of low energy. This systematic effect is particularly troublesome for large $\theta_{13} \geq 1^\circ$ and prevents the use of the Neutrino Factory as a precision facility for large $\theta_{13}$. Such systematic error disappears if the tau contribution to the golden muon sample is taken into account. The fact that the fluxes of the Neutrino Factory are exactly calculable permits the knowledge of the tau sample due to the $\nu_e \to \nu_\tau$ oscillation. We then compute the contribution to the muon sample arising from this sample {\em in terms of the apparent muon energy}. This requires the computation of a migration matrix $M_{ij}$ which describes the contributions of the tau neutrinos of a given energy $E_i$, to the muon neutrinos of an apparent energy $E_j$. We demonstrate that applying $M_{ij}$ to the data permits the full correction of the otherwise intolerable systematic error.
• ### Neutrino phenomenology and stable dark matter with A4(1011.1371)

Feb. 7, 2011 hep-ph
We present a model based on the A4 non-abelian discrete symmetry leading to a predictive five-parameter neutrino mass matrix and providing a stable dark matter candidate. We found an interesting correlation among the atmospheric and the reactor angles which predicts theta_23 ~ pi/4 for very small reactor angle and deviation from maximal atmospheric mixing for large theta_13. Only normal neutrino mass spectrum is possible and the effective mass entering the neutrinoless double beta decay rate is constrained to be |m_ee| > 4 10^{-4} eV.
• ### Fritzsch neutrino mass matrix from S3 symmetry(1005.3482)

Nov. 3, 2010 hep-ph
We present an extension of the Standard Model (SM) based on the discrete flavor symmetry S3 which gives a neutrino mass matrix with two-zero texture of Fritzsch-type and nearly diagonal charged lepton mass matrix. The model is compatible with the normal hierarchy only and predicts the sine squared of the reactor angle to be 0.01 at the best fit values of solar and atmospheric parameters and maximal leptonic CP violation.
• ### Atmospheric lepton fluxes at very high energy(1010.5084)

Oct. 25, 2010 hep-ph, astro-ph.HE
The observation of astrophysical neutrinos requires a detailed understanding of the atmospheric neutrino background. Since neutrinos are produced in meson decays together with a charged lepton, important constraints on this background can be obtained from the measurement of the atmospheric muon flux. Muons, however, can also be produced as mu+ mu- pairs by purely electromagnetic processes. We use the Z-moment method to study and compare the contributions to the atmospheric muon and neutrino fluxes from different sources (pi/K decay, charmed and unflavored hadron decay, and photon conversion into a muon pair). We pay special attention to the contribution from unflavored mesons (eta, eta', rho0, omega and phi). These mesons are abundant in air showers, their lifetimes are much shorter than those of charged pions or kaons, and they have decay branching ratios of order 10^-4 into final states containing a muon pair. We show that they may be the dominant source of muons at E_mu >10^3 TeV.