
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 nonperturbatively, 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 Bparameter 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 modeldependent estimate $B_{CMO} \sim 1  4$,
currently used in phenomenological analyses, and improve the uncertainty on
this important phenomenological quantity.

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 highenergy 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
nonobservation 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.

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 oneloop, the relevant Green's functions with two
(quarkquark) and three (quarkquarkgluon) 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
twistedmass 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 nonperturbatively the mixings of the chromomagnetic operator with
lower dimensional operators, through proper renormalization conditions. For the
first time the 1/a**2divergent mixing of the chromomagnetic operator with the
scalar density has been determined nonperturbatively with high precision.
Moreover, the 1/adivergent mixing with the pseudoscalar density, due to the
breaking of parity within the twistedmass regularization of QCD, has been
calculated nonperturbatively and found to be smaller than its oneloop
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.

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 nonstandard 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.

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 dimensionfive, 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 quarkantiquark (2pt) and the
quarkantiquarkgluon (3pt) 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, nonperturbative
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.

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 oneloop in perturbation theory. The power divergent
coefficients controlling the mixing with operators of lower dimension have been
computed nonperturbatively, 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 Bparameter of the chromomagnetic operator is
B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~14 currently
used in phenomenological analyses.

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 nonstandard 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.

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 nonstandard 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}$.

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, twozero 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 betadecay are also discussed, for both normal and inverted orderings.

We study matrix elements of the "chromomagnetic" operator on the lattice.
This operator is contained in the strangenesschanging 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 2point
(quarkantiquark) and 3point (gluonquarkantiquark) 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.

In this paper we analyze the compatibility of twozero neutrino Majorana
textures with the recent experimental data. Differently from previous works, we
use the experimental data to fix the values of the nonvanishing mass matrix
entries and study in detail the correlations and degree of finetuning 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.

We present a threeflavour fit to the recent \nu\mu > \nu e and \nu\mu >
\nu\mu T2K oscillation data with different models for the neutrinonucleus
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.

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.

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 nexttoleading 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.

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.

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 IDSNF) was
established by the community at the ninth "International Workshop on Neutrino
Factories, superbeams, and beta beams" which was held in Okayama in August
2007. The IDSNF 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 IDSNF 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 IDSNF is, in essence, a siteindependent study.
Example sites, CERN, FNAL, and RAL, have been identified to allow sitespecific
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.

A class of discrete flavorsymmetrybased models predicts constrained
neutrino mass matrix schemes that lead to specific neutrino mass sumrules
(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.

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 nonabelian 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.

We study the contribution of $\nu_e \to \nu_\tau \to \tau \to \mu$
transitions to the wrongsign 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.

We present a model based on the A4 nonabelian discrete symmetry leading to a
predictive fiveparameter 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.

We present an extension of the Standard Model (SM) based on the discrete
flavor symmetry S3 which gives a neutrino mass matrix with twozero texture of
Fritzschtype 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.

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 Zmoment 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.