
In the presence of both direct and indirect unitarity violation in the lepton
mixing matrix, we derive a complete set of series expansion formulas for
neutrino oscillation probabilities in matter of constant density. Expansions in
the mass hierarchy parameter $\alpha \equiv \Delta m_{21}^{2} / \Delta
m_{31}^{2}$ and those unitarity violation parameters $s^{2}_{ij}$ (for i = 1,
2, 3 and j = 4, 5, 6) up to the first order are studied in this paper. We
analyse the accuracy of the analytical series expansion formulas in different
regions of L / E. A detailed numerical analysis is also performed, of which the
different effects of the direct and the indirect unitarity violation are
particularly emphasized. We also study in this paper the summed
$\nu^{}_{\alpha} \rightarrow \nu^{}_{e, \nu, \tau}$ probabilities, whose
deviation from the unity provides a definite signal of the unitarity violation.

The Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton multipurpose
underground liquid scintillator detector, was proposed with the determination
of the neutrino mass hierarchy as a primary physics goal. It is also capable of
observing neutrinos from terrestrial and extraterrestrial sources, including
supernova burst neutrinos, diffuse supernova neutrino background, geoneutrinos,
atmospheric neutrinos, solar neutrinos, as well as exotic searches such as
nucleon decays, dark matter, sterile neutrinos, etc. We present the physics
motivations and the anticipated performance of the JUNO detector for various
proposed measurements. By detecting reactor antineutrinos from two power plants
at 53km distance, JUNO will determine the neutrino mass hierarchy at a 34
sigma significance with six years of running. The measurement of antineutrino
spectrum will also lead to the precise determination of three out of the six
oscillation parameters to an accuracy of better than 1\%. Neutrino burst from a
typical corecollapse supernova at 10 kpc would lead to ~5000
inversebetadecay events and ~2000 allflavor neutrinoproton elastic
scattering events in JUNO. Detection of DSNB would provide valuable information
on the cosmic starformation rate and the average corecollapsed neutrino
energy spectrum. Geoneutrinos can be detected in JUNO with a rate of ~400
events per year, significantly improving the statistics of existing geoneutrino
samples. The JUNO detector is sensitive to several exotic searches, e.g. proton
decay via the $p\to K^++\bar\nu$ decay channel. The JUNO detector will provide
a unique facility to address many outstanding crucial questions in particle and
astrophysics. It holds the great potential for further advancing our quest to
understanding the fundamental properties of neutrinos, one of the building
blocks of our Universe.

According to different effects on neutrino oscillations, the unitarity
violation in the MNSP matrix can be classified into the direct unitarity
violation and the indirect unitarity violation which are induced by the
existence of the light and the heavy sterile neutrinos respectively. Of which
subeV sterile neutrinos are of most interesting. We study in this paper the
possibility of searching for subeV sterile neutrinos in the precision reactor
antineutrino oscillation experiments with three different baselines at around
500 m, 2 km and 60 km. We find that the antineutrino survival probabilities
obtained in the reactor experiments are sensitive only to the direct unitarity
violation and offer very concentrated sensitivity to the two parameters
$\theta^{}_{14}$ and $\Delta m^{2}_{41}$. If such light sterile neutrinos do
exist, the activesterile mixing angle $\theta^{}_{14}$ could be acquired by
the combined rate analysis at all the three baselines and the masssquared
difference $\Delta m^{2}_{41}$ could be obtained by taking the Fourier
transformation to the L / E spectrum. Of course, for such measurements to
succeed, both high energy resolution and large statistics are essentially
important.

We point out that the observed neutrino mixing pattern at low energies is
very likely to originate from the 3 times 3 lepton flavor mixing matrix U which
possesses the exact mutau permutation symmetry U_{mu i} = U_{tau i} (for
i=1,2,3) at a superhigh energy scale Lambda_{mu tau} \sim 10^{14} GeV. The
deviation of theta_{23} from 45^\circ and that of delta from 270^\circ in the
standard parametrization of U are therefore a natural consequence of small
mutau symmetry breaking via the renormalizationgroup equations (RGEs) running
from Lambda_{mu tau} down to the electroweak scale Lambda_{EW} \sim 10^2 GeV.
In fitting current experimental data we find that the RGEcorrected value of
theta_{23} is uniquely correlated with the neutrino mass ordering: theta_{23}
\simeq 42.4^\circ reported by Capozzi et al (or theta_{23} \simeq 48.9^\circ
reported by Forero et al) at Lambda_{EW} can arise from theta_{23} = 45^\circ
at Lambda_{mu tau} in the minimal supersymmetric standard model if the neutrino
mass ordering is inverted (or normal). Accordingly, the preliminary bestfit
results of delta at Lambda_{EW} can also evolve from delta = 270^\circ at
Lambda_{mu tau} no matter whether the massive neutrinos are Dirac or Majorana
particles.

The neutrinodominated accretion flow (NDAF) with accretion rates \dot{M} =
0.01  10 M_{\sun} s^{1} is a plausible candidate for the central engine of
gammaray bursts (GRBs). This hyperaccretion disk is optically thin to
neutrinos in the radial direction, therefore the neutrinos produced at one
radius can travel for a long distance in the disk. Those neutrinos can thus be
absorbed with certain probability by the disk matter at the other radius and
heat the disk there. The effect of this "global neutrino heating" has been
ignored in previous works and is the focus of this paper. We find that around
the "ignition" radius r_{ign}, the global neutrino heating rate could be
comparable to or even larger than the local viscous heating rate thus must be
an important process. Two possible consequences are in order if the "global
neutrino heating" is taken into account: i) the temperature of the disk is
slightly raised and the "ignition" radius r_{ign} slightly shifts to a larger
radius, both lead to the increasing of the total neutrino flux; ii) what is
more interesting is that, the temperature of the ADAF just beyond r_{ign} may
be raised above the virial temperature thus the accretion will be suppressed.
In this case, the activity of the black hole is expected to oscillate between
an active and inactive phases. The timescale of the active phases is estimated
to be \sim 1 second. If the timescale of the inactive phase is comparable to or
less than this value, this intermittent activity may explain the slow
variability component of the GRBs. Selfconsistent global calculations of NDAFs
with the "global neutrino heating" included are required in the future to more
precisely evaluate this effect.

We present an overview on some basic properties of massive neutrinos and
focus on their flavor issues, including the mass spectrum, flavor mixing
pattern and CP violation. The lepton flavor structures are explored by taking
account of the observed value of the smallest neutrino mixing angle
\theta_{13}. The impact of \theta_{13} on the running behaviors of other flavor
mixing parameters is discussed in some detail. The seesawinduced enhancement
of the electromagnetic dipole moments for three Majorana neutrinos is also
discussed in a TeV seesaw scenario.

The recent observation of the smallest neutrino mixing angle $\theta_{13}$ in
the Daya Bay and RENO experiments motivates us to examine whether $\theta_{13}
\simeq 9^\circ$ at the electroweak scale can be generated from $\theta_{13} =
0^\circ$ at a superhighenergy scale via the radiative corrections. We find
that it is difficult but not impossible in the minimal supersymmetric standard
model (MSSM), and a relatively large $\theta_{13}$ may have some nontrivial
impacts on the running behaviors of the other two mixing angles and
CPviolating phases. In particular, we demonstrate that the CPviolating phases
play a crucial role in the evolution of the mixing angles by using the oneloop
renormalizationgroup equations of the Dirac or Majorana neutrinos in the MSSM.
We also take the "correlative" neutrino mixing pattern with $\theta_{12} \simeq
35.3^\circ$, $\theta_{23} = 45^\circ$ and $\theta_{13} \simeq 9.7^\circ$ at a
presumable flavor symmetry scale as an example to illustrate that the three
mixing angles can receive comparably small radiative corrections and thus
evolve to their bestfit values at the electroweak scale if the CPviolating
phases are properly adjusted.

Enlightened by the idea of the 3 times 3 CKM angle matrix proposed recently
by Harrison et al., we introduce the Dirac angle matrix Phi and the Majorana
angle matrix Psi in the lepton sector for Dirac and Majorana neutrinos
respectively. We show that in presence of the CP violation, the angle matrix
Phi or Psi is entirely equivalent to the complex MNS matrix V itself, but has
the advantage of being real, phase rephasing invariant, directly associated to
the leptonic unitarity triangles (UTs) and do not depend on any particular
parametrization of V. In this paper, we further analyzed how the angle matrices
evolve with the energy scale. The oneloop Renormalization Group Equations
(RGEs) of Phi, Psi and some other rephasing invariant parameters are derived
and the numerical analysis is performed to compare between the case of Dirac
and Majorana neutrinos. Different neutrino mass spectra are taken into account
in our calculation. We find that apparently different from the case of Dirac
neutrinos, for Majorana neutrinos the RGevolutions of Phi, Psi and the
Jarlskog strongly depend on the Majoranatype CPviolating parameters and are
quite sensitive to the sign of Delta m^{2}_{31}. They may receive significant
radiative corrections in the MSSM if three neutrino masses are nearly
degenerate.

Current cosmological data drop an interesting hint about the existence of
subeV sterile neutrinos, which should be a part of the cosmic neutrino
background (C$\nu$B). We point out that such light sterile neutrinos may leave
a distinct imprint on the electron energy spectrum in the capture of relic
electron neutrinos by means of radioactive betadecaying nuclei. We examine
possible signals of sterile neutrinos relative to active neutrinos,
characterized by their masses and sensitive to their number densities, in the
reaction $\nu^{}_e + ~^3{\rm H} \to ~^3{\rm He} + e^$ against the
corresponding tritium beta decay. We stress that this kind of direct laboratory
detection of the C$\nu$B and its sterile component might not be hopeless in the
long term.

The 3 X 3 CKM matrix defines six unitarity triangles in the complex plane,
which will be carefully explored in the LHCb experiment and at the SuperB
factory. We calculate the running effects of nine different inner angles and
eighteen different sides of the six triangles from the electroweak scale to a
superhighenergy scale by using the oneloop renormalizationgroup equations,
and demonstrate that all the nine angles are stable against radiative
corrections. In particular, we find that the CPviolating angle \alpha is most
insensitive to the changes of energy scales.

We write out the generic Dirac neutrino mass operator which possesses the
FriedbergLee (FL) symmetry and find that its corresponding neutrino mass
matrix is asymmetric. Following a simple way to break the FL symmetry, we
calculate the neutrino mass eigenvalues and show that the resultant neutrino
mixing pattern is nearly tribimaximal. Imposing the Hermitian condition on the
neutrino mass matrix, we also show that the simplified ansatz is consistent
with current experimental data and favors the normal neutrino mass hierarchy.

We propose a new pattern of the neutrino mixing matrix which can be
parametrized as the product of an arbitrary Hermitian matrix and the wellknown
tribimaximal mixing matrix. In this scenario, nontrivial values of the
smallest neutrino mixing angle \theta_13 and the CPviolating phases entirely
arise from the nonunitary corrections. We present a complete set of series
expansion formulas for neutrino oscillation probabilities both in vacuum and in
matter of constant density. We do a numerical analysis to show the nonunitary
effects on neutrino oscillations. The possibility of determining small
nonunitary perturbations and CPviolating phases is discussed by measuring
neutrino oscillation probabilities and constructing "deformed unitarity
triangles". Some brief comments on the nonunitary neutrino mixing matrix in
the typeII seesaw models are also given.

Current experimental data allow the zero value for one neutrino mass, either
m_1 =0 or m_3 =0. This observation implies that a realistic neutrino mass
texture can be established by starting from the limit (a) m_1 = m_2 =0 and m_3
\neq 0 or (b) m_1 = m_2 \neq 0 and m_3 =0. In both cases, we may introduce a
particular perturbation which ensures the resultant neutrino mixing matrix to
be the tribimaximal mixing pattern or its viable variations. We find that it
is natural to incorporate this kind of neutrino mass matrix in the minimal
TypeII seesaw model with only one heavy righthanded Majorana neutrino N. We
show that it is possible to account for the cosmological baryon number
asymmetry in the m_3 =0 case via thermal leptogenesis, in which the
CPviolating asymmetry of N decays is attributed to the electron flavor.

A natural extension of the standard model to accommodate massive neutrinos is
to introduce one Higgs triplet and three righthanded Majorana neutrinos,
leading to a 6 \times6 neutrino mass matrix. We show that three light Majorana
neutrinos (i.e., the mass eigenstates of \nu_e, \nu_\mu and \nu_\tau) are
exactly massless, if and only if M_L = M_D M^{1}_R M^T_D exactly holds in this
seesaw model. We propose three simple TypeII seesaw scenarios with broken A_4
\times U(1)_X flavor symmetry to interpret the observed neutrino mass spectrum
and neutrino mixing pattern. Such a TeVscale neutrino model can be tested in
two complementary ways: (1) searching for possible collider signatures of
lepton number violation induced by the righthanded Majorana neutrinos and
doublycharged Higgs particles; and (2) searching for possible consequences of
unitarity violation of the 3\times 3 neutrino mixing matrix in the future
longbaseline neutrino oscillation experiments.

We propose two TypeII seesaw scenarios for the neutrino mass matrix in the
leftright symmetric model, in which the Higgs triplet Yukawa coupling matrix
takes the appealing FriedbergLee texture. We show that the nearly
tribimaximal neutrino mixing pattern, which is especially favored by current
neutrino oscillation data, can be obtained from both scenarios. We also show
that the cosmological baryon number asymmetry can naturally be interpreted in
these two scenarios via the flavorindependent leptogenesis mechanism.

A natural extension of the standard $SU(2)_{\rm L} \times U(1)_{\rm Y}$ gauge
model to accommodate massive neutrinos is to introduce one Higgs triplet and
three righthanded Majorana neutrinos, leading to a $6\times 6$ neutrino mass
matrix which contains three $3\times 3$ submatrices $M_{\rm L}$, $M_{\rm D}$
and $M_{\rm R}$. We show that three light Majorana neutrinos (i.e., the mass
eigenstates of $\nu_e$, $\nu_\mu$ and $\nu_\tau$) are exactly massless in this
model, if and only if $M_{\rm L} = M_{\rm D} M_{\rm R}^{1} M_{\rm D}^T$
exactly holds. This nogo theorem implies that small but nonvanishing neutrino
masses may result from a significant but incomplete cancellation between
$M_{\rm L}$ and $M_{\rm D} M_{\rm R}^{1} M_{\rm D}^T$ terms in the TypeII
seesaw formula, provided three righthanded Majorana neutrinos are of ${\cal
O}(1)$ TeV and experimentally detectable at the LHC. We propose three simple
TypeII seesaw scenarios with the $A_4 \times U(1)_{\rm X}$ flavor symmetry to
interpret the observed neutrino mass spectrum and neutrino mixing pattern. Such
a TeVscale neutrino model can be tested in two complementary ways: (1)
searching for possible collider signatures of lepton number violation induced
by the righthanded Majorana neutrinos and doublycharged Higgs particles; and
(2) searching for possible consequences of unitarity violation of the $3\times
3$ neutrino mixing matrix in the future longbaseline neutrino oscillation
experiments.

Current experimental data allow the zero value for one neutrino mass, either
m_1 = 0 or m_3 = 0. This observation implies that a realistic neutrino mass
texture can be established by starting from the limit (a) m_1 = m_2 = 0 and m_3
\neq 0 or (b) m_1 = m_2 \neq 0 and m_3 = 0. In both cases, we may introduce a
particular perturbation which ensures the resultant neutrino mixing matrix to
be the tribimaximal mixing pattern or its viable variations with all entries
being formed from small integers and their square roots. We find that it is
natural to incorporate this kind of neutrino mass matrix in the minimal TypeII
seesaw model with only one heavy righthanded Majorana neutrino N in addition
to the SU(2)_L Higgs triplet \Delta_L. We show that it is possible to account
for the cosmological baryon number asymmetry in the m_3 =0 case via thermal
leptogenesis, in which the oneloop vertex correction to N decays is mediated
by \Delta_L and the CPviolating asymmetry of N decays is attributed to the
electron flavor.

We consider an effective Majorana neutrino mass operator with the
FriedbergLee symmetry; i.e., it is invariant under the transformation
\nu_\alpha > \nu_\alpha + z (for \alpha = e, \mu, \tau) with z being a
spacetime independent constant element of the Grassmann algebra. We show that
this new flavor symmetry can be broken in such a nontrivial way that the
lightest neutrino remains massless but an experimentallyfavored neutrino
mixing pattern is achievable. In particular, we get a novel prediction for the
unknown neutrino mixing angle \theta_13 in terms of two known angles:
\sin\theta_13 = \tan\theta_12 (1 \tan\theta_23)/ (1+\tan\theta_23). The
model can simply be generalized to accommodate CP violation and be combined
with the seesaw mechanism.

Taking the standard parametrization of threeflavor neutrino mixing, we
carefully examine the evolution of three CPviolating phases $(\delta,
\alpha^{}_1, \alpha^{}_2)$ with energy scales in the realistic limit
$\theta^{}_{13} \to 0$. If $m^{}_3$ vanishes, we find that the oneloop
renormalizationgroup equation (RGE) of $\delta$ does not diverge and its
running has no quasifixed point. When $m^{}_3 \neq 0$ holds, we show that the
continuity condition derived by Antusch {\it et al} is always valid, no matter
whether the $\tau$dominance approximation is taken or not. The RGE running of
$\delta$ undergoes a quasifixed point determined by a nontrivial input of
$\alpha^{}_2$ in the limit $m^{}_1 \to 0$. If three neutrino masses are nearly
degenerate, it is also possible to arrive at a quasifixed point in the RGE
evolution of $\delta$ from the electroweak scale to the seesaw scale or vice
versa. Furthermore, the continuity condition and the quasifixed point of
CPviolating phases in another useful parametrization are briefly discussed.

We argue that the tribimaximal neutrino mixing pattern $V_0$ or its
generalized form $V'_0$, which includes two arbitrary Majorana phases of CP
violation, may result from an underlying flavor symmetry at a superhigh energy
scale close to the seesaw scale ($\sim 10^{14}$ GeV). Taking the working
assumption that three neutrino masses are nearly degenerate, we calculate
radiative corrections to $V_0$ and $V'_0$ in their evolution down to the
electroweak scale ($\sim 10^2$ GeV). Three mixing angles of $V_0$ or $V'_0$ are
essentially stable against radiative corrections in the standard model (SM). In
the minimal supersymmetric standard model (MSSM), however, $V_0$ is in general
disfavored and $V'_0$ can be compatible with current neutrino oscillation data
if its two Majorana phases $\alpha^{}_1$ and $\alpha^{}_2$ are properly
finetuned. We also find that it is possible to radiatively generate the
CPviolating phase $\delta$ from $\alpha^{}_1$ and $\alpha^{}_2$, and \delta
may keep on staying at its fixed point in either the SM or the MSSM.

Three CPviolating phases of the $3\times 3$ lepton flavor mixing matrix $V$
are entangled with one another in the renormalizationgroup evolution from the
seesaw scale ($\Lambda_{\rm SS} \sim 10^{14}$ GeV) to the electroweak scale
($\Lambda_{\rm EW} \sim 10^2$ GeV). Concerning the Dirac phase $\delta$, we
show that $\delta =90^\circ$ at $\Lambda_{\rm EW}$ can be radiatively generated
from $\delta =0^\circ$ at $\Lambda_{\rm SS}$ in the minimal supersymmetric
standard model, if three neutrino masses are nearly degenerate. As for the
Majorana phases $\rho$ and $\sigma$, it is also possible to radiatively
generate $\rho =90^\circ$ or $\sigma = 90^\circ$ at $\Lambda_{\rm EW}$ from
$\rho =0^\circ$ or $\sigma = 0^\circ$ at $\Lambda_{\rm SS}$. The oneloop
renormalizationgroup equations for the Jarlskog invariant and two offdiagonal
asymmetries of $V$ are derived, and their running behaviors from $\Lambda_{\rm
SS}$ to $\Lambda_{\rm EW}$ are numerically illustrated.