
We present an uptodate global analysis of data coming from neutrino
oscillation and nonoscillation experiments, as available in April 2018, within
the standard framework including three massive and mixed neutrinos. We discuss
in detail the status of the threeneutrino (3nu) massmixing parameters, both
known and unknown. Concerning the latter, we find that: normal ordering (NO) is
favored over inverted ordering (IO) at 3sigma level; the Dirac CP phase is
constrained within ~15% (~9%) uncertainty in NO (IO) around nearlymaximal
CPviolating values; the octant of the largest mixing angle and the absolute
neutrino masses remain undetermined. We briefly comment on other unknowns
related to theoretical and experimental uncertainties (within 3nu) or possible
new states and interactions (beyond 3nu).

In the Ratra scenario of inflationary magnetogenesis, the kinematic coupling
between the photon and the inflaton undergoes a nonanalytical jump at the end
of inflation. Using smooth interpolating analytical forms of the coupling
function, we show that such unphysical jump does not invalidate the main
prediction of the model, which still represents a viable mechanism for
explaining cosmic magnetization. Nevertheless, there is a spurious result
associated with the nonanaliticity of the coupling, to wit, the prediction that
the spectrum of created photons has a powerlaw decay in the ultraviolet
regime. This issue is discussed using both semiclassical approximation and
smooth coupling functions.

Within the standard 3nu massmixing framework, we present an uptodate
global analysis of neutrino oscillation data (as of January 2016), including
the latest available results from experiments with atmospheric neutrinos
(SuperKamiokande and IceCube DeepCore), at accelerators (first T2K antinu and
NOvA nu runs in both appearance and disappearance mode), and at shortbaseline
reactors (Daya Bay and RENO far/near spectral ratios), as well as a reanalysis
of older KamLAND data in the light of the "bump" feature recently observed in
reactor spectra. We discuss improved constraints on the five known oscillation
parameters (delta m^2, Delta m^2, sin^2theta_12, sin^2theta_13,
sin^2theta_23), and the status of the three remaining unknown parameters: the
mass hierarchy, the theta_23 octant, and the possible CPviolating phase delta.
With respect to previous global fits, we find that the reanalysis of KamLAND
data induces a slight decrease of both delta m^2 and sin^2theta_12, while the
latest accelerator and atmospheric data induce a slight increase of Delta
m^2. Concerning the unknown parameters, we confirm the previous intriguing
preference for negative values of sin(delta) [with bestfit values around
sin(delta) ~ 0.9], but we find no statistically significant indication about
the theta_23 octant or the mass hierarchy (normal or inverted). Assuming an
alternative (socalled LEM) analysis of NOvA data, some delta ranges can be
excluded at >3 sigma, and the normal mass hierarchy appears to be slightly
favored at 90% C.L. We also describe in detail the covariances of selected
pairs of oscillation parameters. Finally, we briefly discuss the implications
of the above results on the three nonoscillation observables sensitive to the
(unknown) absolute nu mass scale: the sum of nu masses, the effective nu_e
mass, and the effective Majorana mass.

Nuclear reactors provide intense sources of electron antineutrinos,
characterized by fewMeV energy E and unoscillated spectral shape Phi(E).
Highstatistics observations of reactor neutrino oscillations over
mediumbaseline distances L ~ O(50) km would provide unprecedented
opportunities to probe both the longwavelength massmixing parameters (delta
m^2 and theta_12) and the shortwavelength ones (Delta m^2 and theta_13),
together with the subtle interference effects associated with the neutrino mass
hierarchy (either normal or inverted). In a given experimental setting  here
taken as in the JUNO project for definiteness  the achievable hierarchy
sensitivity and parameter accuracy depend not only on the accumulated
statistics but also on systematic uncertainties, which include (but are not
limited to) the massmixing priors and the normalizations of signals and
backgrounds. We examine, in addition, the effect of introducing smooth
deformations of the detector energy scale, E > E'(E), and of the reactor flux
shape, Phi(E) > Phi'(E), within reasonable error bands inspired by
stateoftheart estimates. It turns out that energyscale and fluxshape
systematics can noticeably affect the performance of a JUNOlike experiment,
both on the hierarchy discrimination and on precision oscillation physics. It
is shown that a significant reduction of the assumed energyscale and
fluxshape uncertainties (by, say, a factor of 2) would be highly beneficial to
the physics program of mediumbaseline reactor projects. Our results also shed
some light on the role of the inversebeta decay threshold, of geoneutrino
backgrounds, and of matter effects in the analysis of future reactor
oscillation data.

A new general purpose fixed target facility is proposed at the CERN SPS
accelerator which is aimed at exploring the domain of hidden particles and make
measurements with tau neutrinos. Hidden particles are predicted by a large
number of models beyond the Standard Model. The high intensity of the SPS
400~GeV beam allows probing a wide variety of models containing light
longlived exotic particles with masses below ${\cal O}$(10)~GeV/c$^2$,
including very weakly interacting lowenergy SUSY states. The experimental
programme of the proposed facility is capable of being extended in the future,
e.g. to include direct searches for Dark Matter and Lepton Flavour Violation.

The proposed PINGU project (Precision IceCube Next Generation Upgrade) is
expected to collect O(10^5) atmospheric muon and electron neutrino in a few
years of exposure, and to probe the neutrino mass hierarchy through its imprint
on the event spectra in energy and direction. In the presence of nonnegligible
and partly unknown shape systematics, the analysis of highstatistics spectral
variations will face subtle challenges that are largely unprecedented in
neutrino physics. We discuss these issues both on general grounds and in the
currently envisaged PINGU configuration, where we find that possible shape
uncertainties at the (few) percent level can noticeably affect the sensitivity
to the hierarchy. We also discuss the interplay between the mixing angle
theta_23 and the PINGU sensitivity to the hierarchy. Our results suggest that
more refined estimates of spectral uncertainties are needed in nextgeneration,
largevolume atmospheric neutrino experiments.

The standard threeneutrino (3nu) oscillation framework is being increasingly
refined by results coming from different sets of experiments, using neutrinos
from solar, atmospheric, accelerator and reactor sources. At present, each of
the known oscillation parameters [the two squared mass gaps (delta m^2, Delta
m^2) and the three mixing angles (theta_12}, theta_13, theta_23)] is dominantly
determined by a single class of experiments. Conversely, the unknown parameters
[the mass hierarchy, the theta_23 octant and the CPviolating phase delta] can
be currently constrained only through a combined analysis of various
(eventually all) classes of experiments. In the light of recent new results
coming from reactor and accelerator experiments, and of their interplay with
solar and atmospheric data, we update the estimated Nsigma ranges of the known
3nu parameters, and revisit the status of the unknown ones. Concerning the
hierarchy, no significant difference emerges between normal and inverted mass
ordering. A slight overall preference is found for theta_23 in the first octant
and for nonzero CP violation with sin delta < 0; however, for both parameters,
such preference exceeds 1 sigma only for normal hierarchy. We also discuss the
correlations and stability of the oscillation parameters within different
combinations of data sets.

Proposed mediumbaseline reactor neutrino experiments offer unprecedented
opportunities to probe, at the same time, the massmixing parameters which
govern $\nu_e$ oscillations both at short wavelength (delta m^2 and theta_{12})
and at long wavelength (Delta m^2 and theta_{13}), as well as their tiny
interference effects related to the mass hierarchy (i.e., the relative sign of
Delta m^2 and delta m^2). In order to take full advantage of these
opportunities, precision calculations and refined statistical analyses of event
spectra are required. In such a context, we revisit several input ingredients,
including: nucleon recoil in inverse beta decay and its impact on energy
reconstruction and resolution, hierarchy and matter effects in the oscillation
probability, spread of reactor distances, irreducible backgrounds from
geoneutrinos and from far reactors, and degeneracies between energy scale and
spectrum shape uncertainties. We also introduce a continuous parameter alpha,
which interpolates smoothly between normal hierarchy (alpha=+1) and inverted
hierarchy (alpha=1). The determination of the hierarchy is then transformed
from a test of hypothesis to a parameter estimation, with a sensitivity given
by the distance of the true case (either alpha=+1 or alpha=1) from the
undecidable case (alpha=0). Numerical experiments are performed for the
specific set up envisaged for the JUNO project, assuming a realistic sample of
O(10^5) reactor events. We find a typical sensitivity of ~2 sigma to the
hierarchy in JUNO, which, however, can be challenged by energy scale and
spectrum shape systematics, whose possible conspiracy effects are investigated.
The prospective accuracy reachable for the other massmixing parameters is also
discussed.

[Abridged] We use data on massive galaxy clusters ($M_{\rm cluster} > 8
\times 10^{14} h^{1} M_\odot$ within a comoving radius of $R_{\rm cluster} =
1.5 h^{1}\Mpc$) in the redshift range $0.05 \lesssim z \lesssim 0.83$ to place
constraints, simultaneously, on the nonrelativistic matter density parameter
$\Omega_m$, on the amplitude of mass fluctuations $\sigma_8$, on the index $n$
of the powerlaw spectrum of the density perturbations, and on the Hubble
constant $H_0$, as well as on the equationofstate parameters $(w_0,w_a)$ of a
smooth dark energy component.
For the first time, we properly take into account the dependence on redshift
and cosmology of the quantities related to cluster physics: the critical
density contrast, the growth factor, the mass conversion factor, the virial
overdensity, the virial radius and, most importantly, the cluster number count
derived from the observational temperature data.
We show that, contrary to previous analyses, cluster data alone prefer low
values of the amplitude of mass fluctuations, $\sigma_8 \leq 0.69 (1\sigma
C.L.)$, and large amounts of nonrelativistic matter, $\Omega_m \geq 0.38
(1\sigma C.L.)$, in slight tension with the $\Lambda$CDM concordance
cosmological model, though the results are compatible with $\Lambda$CDM at
$2\sigma$. In addition, we derive a $\sigma_8$ normalization relation,
$\sigma_8 \Omega_m^{1/3} = 0.49 \pm 0.06 (2\sigma C.L.)$.

We perform a global analysis of neutrino oscillation data, including
highprecision measurements of the neutrino mixing angle theta_13 at reactor
experiments, which have confirmed previous indications in favor of theta_13>0.
Recent data presented at the Neutrino 2012 Conference are also included. We
focus on the correlations between theta_13 and the mixing angle theta_23, as
well as between theta_13 and the neutrino CPviolation phase delta. We find
interesting indications for theta_23< pi/4 and possible hints for delta ~ pi,
with no significant difference between normal and inverted mass hierarchy.

The neutrino mixing angle theta(13) is at the focus of current neutrino
research. From a global analysis of the available oscillation data in a
3neutrino framework, we previously reported [Phys. Rev. Lett. 101, 141801
(2008)] hints in favor of theta(13)>0 at the 90 % C.L. Such hints are
consistent with the recent indications of nu(mu)>nu(e) appearance in the T2K
and MINOS longbaseline accelerator experiments. Our global analysis of all the
available data currently provides >3 sigma evidence for nonzero theta(13), with
1sigma ranges sin^2 theta(13) = 0.021+0.007 or 0.025+0.007, depending on
reactor neutrino flux systematics. Updated ranges are also reported for the
other 3neutrino oscillation parameters (delta m^2, sin^2 theta(12)) and (Delta
m^2, sin^2 theta(23)).

We analyze the magnituderedshift data of type Ia supernovae included in the
Union and Union2 compilations in the framework of an anisotropic Bianchi type I
cosmological model and in the presence of a dark energy fluid with anisotropic
equation of state. We find that the amount of deviation from isotropy of the
equation of state of dark energy, the skewness \delta, and the present level of
anisotropy of the largescale geometry of the Universe, the actual shear
\Sigma_0, are constrained in the ranges 0.16 < \delta < 0.12 and 0.012 <
\Sigma_0 < 0.012 (1\sigma C.L.) by Union2 data. Supernova data are then
compatible with a standard isotropic universe (\delta = \Sigma_0 = 0), but a
large level of anisotropy, both in the geometry of the Universe and in the
equation of state of dark energy, is allowed.

At the previous Venice meeting NOVE 2008, we discussed possible hints in
favor of a nonzero value for the unknown neutrino mixing angle theta(13),
emerging from the combination of solar and longbaseline reactor data, as well
as from the combination of atmospheric, CHOOZ and longbaseline accelerator
nu_mu>nu_mu data. Recent MINOS 2009 results in the nu_mu>nu_e appearance
channel also seem to support such hints. A combination of all current
oscillation data provides, as preferred range, sin^2 theta(13) = 0.02 + 0.01
(1\sigma). We review several issues raised by such hints in the last year, and
comment on their possible nearfuture improvements and tests.

Nailing down the unknown neutrino mixing angle theta_13 is one of the most
important goals in current lepton physics. In this context, we perform a global
analysis of neutrino oscillation data, focusing on theta_13, and including
recent results [Neutrino 2008, Proceedings of the XXIII International
Conference on Neutrino Physics and Astrophysics, Christchurch, New Zealand,
2008 (unpublished)]. We discuss two converging hints of theta_13>0, each at the
level of ~1sigma: an older one coming from atmospheric neutrino data, and a
newer one coming from the combination of solar and longbaseline reactor
neutrino data. Their combination provides the global estimate sin^2(theta_13) =
0.016 + 0.010 (1sigma), implying a preference for \theta_13>0 with
nonnegligible statistical significance (~90% C.L.). We discuss possible
refinements of the experimental data analyses, which might sharpen such
intriguing indication.

We present updated values for the massmixing parameters relevant to neutrino
oscillations, with particular attention to emerging hints in favor of
theta_13>0. We also discuss the status of absolute neutrino mass observables,
and a possible approach to constrain theoretical uncertainties in neutrinoless
double beta decay. Desiderata for all these issues are also briefly mentioned.

In corecollapse supernovae, neutrinos and antineutrinos are initially
subject to significant selfinteractions induced by weak neutral currents,
which may induce strongcoupling effects on the flavor evolution (collective
transitions). The interpretation of the effects is simplified when selfinduced
collective transitions are decoupled from ordinary matter oscillations, as for
the matter density profile that we discuss. In this case, approximate
analytical tools can be used (pendulum analogy, swap of energy spectra). For
inverted neutrino mass hierarchy, the sequence of effects involves:
synchronization, bipolar oscillations, and spectral split. Our simulations
shows that the main features of these regimes are not altered when passing from
simplified (angleaveraged) treatments to full, multiangle numerical
experiments.

In this followup to Phys. Rev. D 75, 053001 (2007) [arXiv:hepph/0608060] we
report updated constraints on neutrino massmixing parameters, in light of
recent neutrino oscillation data (KamLAND, SNO, and MINOS) and cosmological
observations (WMAP 5year and other data). We discuss their interplay with the
final 0nu2beta decay results in 76Ge claimed by part of the HeidelbergMoscow
Collaboration, using recent evaluations of the corresponding nuclear matrix
elements, and their uncertainties. We also comment on the 0nu2beta limits in
130Te recently set by Cuoricino, and on prospective limits or signals from the
KATRIN experiment.

In the dense supernova core, selfinteractions may align the flavor
polarization vectors of neutrinos and antineutrinos, and induce collective
flavor transformations. Different alignment ansatzes are known to describe
approximately the phenomena of synchronized or bipolar oscillations, and the
split of neutrino energy spectra. We discuss another phenomenon observed in
some numerical experiments in inverted hierarchy, showing features akin to a
lowenergy split of antineutrino spectra. The phenomenon appears to be
approximately described by another alignment ansatz which, in the considered
scenario, reduces the (nonadiabatic) dynamics of all energy modes to only two
neutrino plus two antineutrino modes. The associated spectral features,
however, appear to be fragile when passing from single to multiangle
simulations.

It has been speculated that quantum gravity might induce a "foamy" spacetime
structure at small scales, randomly perturbing the propagation phases of
freestreaming particles (such as kaons, neutrons, or neutrinos). Particle
interferometry might then reveal nonstandard decoherence effects, in addition
to standard ones (due to, e.g., finite source size and detector resolution.) In
this work we discuss the phenomenology of such nonstandard effects in the
propagation of electron neutrinos in the Sun and in the longbaseline reactor
experiment KamLAND, which jointly provide us with the best available probes of
decoherence at neutrino energies E ~ few MeV. In the solar neutrino case, by
means of a perturbative approach, decoherence is shown to modify the standard
(adiabatic) propagation in matter through a calculable damping factor. By
assuming a powerlaw dependence of decoherence effects in the energy domain
(E^n with n = 0,+/1,+/2), theoretical predictions for twofamily neutrino
mixing are compared with the data and discussed. We find that neither solar nor
KamLAND data show evidence in favor of nonstandard decoherence effects, whose
characteristic parameter gamma_0 can thus be significantly constrained. In the
"Lorentzinvariant" case n=1, we obtain the upper limit gamma_0<0.78 x 10^26
GeV at 95% C.L. In the specific case n=2, the constraints can also be
interpreted as bounds on possible matter density fluctuations in the Sun, which
we improve by a factor of ~ 2 with respect to previous analyses.

In the light of recent neutrino oscillation and nonoscillation data, we
revisit the phenomenological constraints applicable to three observables
sensitive to absolute neutrino masses: The effective neutrino mass in single
beta decay (m_beta); the effective Majorana neutrino mass in neutrinoless
double beta decay (m_2beta); and the sum of neutrino masses in cosmology
(Sigma). In particular, we include the constraints coming from the first Main
Injector Neutrino Oscillation Search (MINOS) data and from the Wilkinson
Microwave Anisotropy Probe (WMAP) threeyear (3y) data, as well as other
relevant cosmological data and priors. We find that the largest neutrino
squared mass difference is determined with a 15% accuracy (at 2sigma) after
adding MINOS to world data. We also find upper bounds on the sum of neutrino
masses Sigma ranging from ~2 eV (WMAP3y data only) to ~0.2 eV (all
cosmological data) at 2sigma, in agreement with previous studies. In addition,
we discuss the connection of such bounds with those placed on the matter power
spectrum normalization parameter sigma_8. We show how the partial degeneracy
between Sigma and sigma_8 in WMAP3y data is broken by adding further
cosmological data, and how the overall preference of such data for relatively
high values of sigma_8 pushes the upper bound of Sigma in the subeV range.
Finally, for various combination of data sets, we revisit the (in)compatibility
between current Sigma and m_2beta constraints (and claims), and derive
quantitative predictions for future single and double beta decay experiments.

John N. Bahcall championed solar neutrino physics for many years. Thanks to
his pioneering and longlasting contributions, this field of research has not
only reached maturity, but has also opened a new window on physics beyond the
standard electroweak model through the phenomenon of neutrino flavor
oscillations. We briefly outline some recent accomplishments in the field, and
also discuss a couple of issues that do not seem to fit in the ``standard
picture,'' namely, the chemical controversy at the solar surface, and possible
implications of recent gallium radioactive source experiments.

We present a comprehensive phenomenological analysis of a vast amount of data
from neutrino flavor oscillation and nonoscillation searches, performed within
the standard scenario with three massive and mixed neutrinos, and with
particular attention to subleading effects. The detailed results discussed in
this review represent a stateoftheart, accurate and uptodate (as of August
2005) estimate of the threeneutrino massmixing parameters.

We present a brief review of the current status of neutrino mass and mixing
parameters, based on a comprehensive phenomenological analysis of neutrino
oscillation and nonoscillation searches, within the standard threeneutrino
mixing framework.

We derive upper limits on the sum of neutrino masses from an updated
combination of data from Cosmic Microwave Background experiments and Galaxy
Redshifts Surveys. The results are discussed in the context of threeflavor
neutrino mixing and compared with neutrino oscillation data, with upper limits
on the effective neutrino mass in Tritium beta decay from the Mainz and Troitsk
experiments and with the claimed lower bound on the effective Majorana neutrino
mass in neutrinoless double beta decay from the HeidelbergMoscow experiment.

In the context of threeflavor neutrino mixing, we present a thorough study
of the phenomenological constraints applicable to three observables sensitive
to absolute neutrino masses: The effective neutrino mass in Tritium beta decay
(m_beta); the effective Majorana neutrino mass in neutrinoless double beta
decay (m_2beta); and the sum of neutrino masses in cosmology (Sigma). We
discuss the correlations among these variables which arise from the combination
of all the available neutrino oscillation data, in both normal and inverse
neutrino mass hierarchy. We set upper limits on m_beta by combining updated
results from the Mainz and Troitsk experiments. We also consider the latest
results on m_2beta from the HeidelbergMoscow experiment, both with and without
the lower bound claimed by such experiment. We derive upper limits on Sigma
from an updated combination of data from the Wilkinson Microwave Anisotropy
Probe (WMAP) satellite and the 2 degrees Fields (2dF) Galaxy Redshifts Survey,
with and without Lymanalpha forest data from the Sloan Digital Sky Survey
(SDSS), in models with a nonzero running of the spectral index of primordial
inflationary perturbations. The results are discussed in terms of
twodimensional projections of the globally allowed region in the
(m_beta,m_2beta,Sigma) parameter space, which neatly show the relative impact
of each data set. In particular, the (in)compatibility between Sigma and
m_2beta constraints is highlighted for various combinations of data. We also
briefly discuss how future neutrino data (both oscillatory and nonoscillatory)
can further probe the currently allowed regions.