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

The current and upcoming generation of Very Large Volume Neutrino Telescopes
 collecting unprecedented quantities of neutrino events  can be used to
explore subtle effects in oscillation physics, such as (but not restricted to)
the neutrino mass ordering. The sensitivity of an experiment to these effects
can be estimated from Monte Carlo simulations. With the very high number of
events that will be collected, there is a tradeoff between the computational
expense of running such simulations and the inherent statistical uncertainty in
the determined values. In such a scenario, it becomes impractical to produce
and use adequatelysized sets of simulated events to use with traditional
methods, such as Monte Carlo weighting. In this work we present a staged
approach to the generation of binned event distributions in order to overcome
these challenges. By combining multiple integration and smoothing techniques
which address limited statistics from simulation it arrives at reliable
analysis results using modest computational resources.

Papers on research & development towards IceCubeGen2, the next generation
neutrino observatory at South Pole, submitted to the 35th International Cosmic
Ray Conference (ICRC 2017, Busan, South Korea) by the IceCubeGen2
Collaboration.

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.

This document presents the Conceptual Design Report (CDR) put forward by an
international neutrino community to pursue the Deep Underground Neutrino
Experiment at the LongBaseline Neutrino Facility (LBNF/DUNE), a groundbreaking
science experiment for longbaseline neutrino oscillation studies and for
neutrino astrophysics and nucleon decay searches. The DUNE far detector will be
a very large modular liquid argon timeprojection chamber (LArTPC) located deep
underground, coupled to the LBNF multimegawatt wideband neutrino beam. DUNE
will also have a highresolution and highprecision near detector.

Papers submitted to the 34th International Cosmic Ray Conference (ICRC 2015,
The Hague) by the IceCubeGen2 Collaboration.

We investigate the impact of light ($\sim$ eV) sterile neutrinos in the
longbaseline experiment T2K. We show that, within the 3+1 scheme, for
massmixing parameters suggested by the shortbaseline anomalies, the
interference among the sterile and the atmospheric oscillation frequencies
induces a new term in the $\nu_\mu \to \nu_e$ transition probability, which has
the same order of magnitude of the standard 3flavor solaratmospheric
interference term. We show, for the first time, that current T2K data, taken
together with the results of the $\theta_{13}$dedicated reactor experiments,
are sensitive to two of the three CPviolating phases involved in the 3+1
scheme. Both the standard CPphase and the new one ($\delta_{13} \equiv \delta$
and $\delta_{14}$ in our parameterization choice) tend to have a common best
fit value $\delta_{13} \simeq \delta_{14} \simeq \pi/2$. Quite intriguingly,
the inclusion of sterile neutrino effects leads to better agreement between the
two estimates of $\theta_{13}$ obtained, respectively, from reactors and T2K,
which in the 3flavor framework are slightly different.

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.

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

The KamLAND and Borexino experiments have detected electron antineutrinos
produced in the decay chains of natural thorium and uranium (Th and U
geoneutrinos). We analyze the energy spectra of current geoneutrino data in
combination with solar and longbaseline reactor neutrino data, with
marginalized threeneutrino oscillation parameters. We consider the case with
unconstrained Th and U event rates in KamLAND and Borexino, as well as cases
with fewer degrees of freedom, as obtained by successively assuming for both
experiments a common Th/U ratio, a common scaling of Th+U event rates, and a
chondritic Th/U value. In combination, KamLAND and Borexino can reject the null
hypothesis (no geoneutrino signal) at 5 sigma. Interesting bounds or
indications emerge on the Th+U geoneutrino rates and on the Th/U ratio, in
broad agreement with typical Earth model expectations. Conversely, the results
disfavor the hypothesis of a georeactor in the Earth's core, if its power
exceeds a few TW. The interplay of KamLAND and Borexino geoneutrino data is
highlighted.

Solar and KamLAND data are in slight tension when interpreted in the standard
twoflavor oscillations framework and this may be alleviated allowing for a
nonzero value of the mixing angle theta_13. Here we show that, likewise,
nonstandard flavorchanging interactions (FCI), possibly intervening in the
propagation of solar neutrinos, are equally able to alleviate this tension and
therefore constitute a potential source of confusion in the determination of
theta_13. By performing a full threeflavor analysis of solar and KamLAND data
in presence of FCI we provide a quantitative description of the degeneracy
existing between theta_13 and the vectorial coupling eps_e\tau^dV
characterizing the nonstandard transitions between nu_e and nu_tau in the
forward scattering process with dtype quarks. We find that couplings with
magnitude eps_e\tau^dV ~ 10%, compatible with the existing bounds, can mimic
the nonzero values of theta_13 indicated by the latest analyses.

Most neutrino mass extensions of the standard electroweak model entail
nonstandard interactions which, in the low energy limit, can be parametrized
in term of effective fourfermion operators $\nu_\alpha \nu_\beta \bar f f $.
Typically of subweak strength, $\epsilon_{\alpha \beta} G_F$, these are
characterized by dimensionless coupling parameters, $\epsilon_{\alpha \beta}$,
which may be relatively sizeable in a wide class of schemes. Here we focus on
nonuniversal (NU) flavor conserving couplings ($\alpha = \beta$) with
electrons ($f = e$) and analyse their impact on the phenomenology of solar
neutrinos. We consistently take into account their effect both at the level of
propagation where they modify the standard MSW behavior, and at the level of
detection, where they affect the cross section of neutrino elastic scattering
on electrons. We find limits which are comparable to other existing
modelindependent constraints.

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 light of recent findings which seem to disfavor a scenario with (warm)
dark matter entirely constituted of sterile neutrinos produced via the
DodelsonWidrow (DW) mechanism, we investigate the constraints attainable for
this mechanism by relaxing the usual hypothesis that the relic neutrino
abundance must necessarily account for all of the dark matter. We first study
how to reinterpret the limits attainable from Xray nondetection and
Lymanalpha forest measurements in the case that sterile neutrinos constitute
only a fraction fs of the total amount of dark matter. Then, assuming that
sterile neutrinos are generated in the early universe solely through the DW
mechanism, we show how the Xray and Lymanalpha results jointly constrain the
massmixing parameters governing their production. Furthermore, we show how the
same data allow us to set a robust upper limit fs < 0.7 at the 2 sigma level,
rejecting the case of dominant dark matter (fs = 1) at the ~ 3 sigma level.

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.

Geoneutrinos emitted by heatproducing elements (U, Th and K) represent a
unique probe of the Earth interior. The characterization of their fluxes is
subject, however, to rather large and highly correlated uncertainties. The
geochemical covariance of the U, Th and K abundances in various Earth
reservoirs induces positive correlations among the associated geoneutrino
fluxes, and between these and the radiogenic heat. Massbalance constraints in
the Bulk Silicate Earth (BSE) tend instead to anticorrelate the radiogenic
element abundances in complementary reservoirs. Experimental geoneutrino
observables may be further (anti)correlated by instrumental effects. In this
context, we propose a systematic approach to covariance matrices, based on the
fact that all the relevant geoneutrino observables and constraints can be
expressed as linear functions of the U, Th and K abundances in the Earth's
reservoirs (with relatively wellknown coefficients). We briefly discuss here
the construction of a tentative "geoneutrino source model" (GNSM) for the U,
Th, and K abundances in the main Earth reservoirs, based on selected
geophysical and geochemical data and models (when available), on plausible
hypotheses (when possible), and admittedly on arbitrary assumptions (when
unavoidable). We use then the GNSM to make predictions about several
experiments ("forward approach"), and to show how future data can constrain  a
posteriori  the error matrix of the model itself ("backward approach"). The
method may provide a useful statistical framework for evaluating the impact and
the global consistency of prospective geoneutrino measurements and Earth
models.

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.

The Kamioka Liquid scintillator AntiNeutrino Detector (KamLAND) is sensitive
to the neutrino event spectrum from (mainly Japanese) nuclear reactors in both
the energy domain and the time domain. While the energy spectrum of KamLAND
events allows the determination of the neutrino oscillation parameters, the
time spectrum can be used to monitor known and unknown neutrino sources. By
using available monthlybinned data on eventbyevent energies in KamLAND and
on reactor powers in Japan, we perform a likelihood analysis of the neutrino
event spectra in energy and time, and find significant indications in favor of
time variations of the known reactor sources, as compared with the hypothetical
case of constant reactor neutrino flux. We also find that the KamLAND data
place interesting upper limits on the power of a speculative nuclear reactor
operating in the Earth's core (the socalled georeactor); such limits are
strengthened by including solar neutrino constraints on the neutrino mass and
mixing parameters. Our results corroborate the standard interpretation of the
KamLAND signal as due to oscillating neutrinos from known reactor sources.

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.

The Kamioka Liquid scintillator AntiNeutrino Detector (KamLAND) is
collecting antineutrino events generated by nuclear reactors and by Thorium and
Uranium decay in the Earth interior. We comment on a systematic approach to the
evaluation of the geoneutrino contribution and of its uncertainties in
KamLAND, taking into account geophysical and geochemical indications,
estimates, and data. The results can help to improve both the neutrino
oscillation analysis and the knowledge of the Earth interior. Input and
desiderata for future geoneutrino analyses are identified.

Realtime neutrino oscillation experiments such as SuperKamiokande (SK), the
Sudbury Neutrino Observatory (SNO), the Kamioka Liquid scintillator
AntiNeutrino Detector (KamLAND), and Borexino, can detect time variations of
the neutrino signal, provided that the statistics is sufficiently high. We
quantify this statement by means of a simple unbinned test, whose sensitivity
depends on the variance of the signal in the time domain, as well as on the
total number of signal and background events. The test allows a unified
discussion of the statistical uncertainties affecting current or future
measurements of eccentricityinduced variations and of daynight asymmetries
(in SK, SNO, and Borexino), as well as of reactor power variations (in
KamLAND).

We review the status of the neutrino oscillation physics (as of June 2003),
with a particular emphasis on the present knowledge of the neutrino massmixing
parameters in a three generation approach. We consider first the nu_mu>nu_tau
flavor transitions of atmospheric neutrinos. It is found that standard
oscillations provide the best description of the SK+K2K data, and that the
associated massmixing parameters are determined at 1 sigma (and dof=1) as:
Delta m^2=(2.6 +0.4) x 10^3 eV^2 and sin^2(2theta)=1.00+0.000.05. Such
indications, presently dominated by SK, could be strengthened by further K2K
data. Then we analyze the energy spectrum of reactor neutrino events recently
observed at KamLAND and combine them with solar and terrestrial neutrino data.
We find that the solution to the solar neutrino problem at large mixing angle
(LMA) is basically split into two subregions, that we denote as LMAI and
LMAII. The LMAI solution, characterized by lower values of the squared
neutrino mass gap, is favored by the global data fit. Finally, we briefly
illustrate how prospective data from the SNO and KamLAND can increase our
confidence in the occurrence of standard matter effects in the Sun, which are
starting to emerge from current data.