
We borrow the general idea of renormalizationgroup equations (RGEs) to
understand how neutrino masses and flavor mixing parameters evolve when
neutrinos propagate in a medium, highlighting a meaningful possibility that the
genuine flavor quantities in vacuum can be extrapolated from their
mattercorrected counterparts to be measured in some realistic neutrino
oscillation experiments. Taking the matter parameter $a \equiv 2\sqrt{2} \
G^{}_{\rm F} N^{}_e E$ to be an arbitrary scalelike variable with $N^{}_e$
being the net electron number density and $E$ being the neutrino beam energy,
we derive a complete set of differential equations for the effective neutrino
mixing matrix $V$ and the effective neutrino masses $\widetilde{m}^{}_i$ (for
$i = 1, 2, 3$). Given the standard parametrization of $V$, the RGEs for
$\{\widetilde{\theta}^{}_{12}, \widetilde{\theta}^{}_{13},
\widetilde{\theta}^{}_{23}, \widetilde{\delta}\}$ in matter are formulated for
the first time. We demonstrate some useful differential invariants which retain
the same form from vacuum to matter, including the wellknown Naumov and Toshev
relations. The RGEs of the partial $\mu$$\tau$ asymmetries, the offdiagonal
asymmetries and the sides of unitarity triangles of $V$ are also obtained as a
byproduct.

In an analytical way of studying matter effects on neutrino oscillations, the
Naumov and Toshev relations have been derived to respectively link the Jarlskog
invariant of CP violation and the Dirac phase in the standard parametrization
of the $3\times 3$ flavor mixing matrix to their mattercorrected counterparts.
Here we show that there exist similar relations for Dirac neutrinos and charged
leptons evolving with energy scales via the oneloop renormalizationgroup (RG)
equations in the taudominance approximation, and for the running behaviors of
up and downtype quarks in the topdominance approximation, provided a
different parametrization is taken into account.

We embed $\mu\tau$ reflection symmetry into the minimal seesaw formalism,
where two righthanded neutrinos are added to the Standard Model of particle
physics. Assuming that both the left and righthanded neutrino fields
transform under $\mu\tau$ reflection symmetry, we obtain the required forms of
the neutrino Dirac mass matrix and the Majorana mass matrix for the
righthanded neutrinos. To investigate the neutrino phenomenology at low
energies, we first consider the breaking of $\mu\tau$ reflection symmetry due
to the renormalization group running, and then systematically study various
breaking schemes by introducing explicit breaking terms at high energies.

If finite but tiny masses of the three active neutrinos are generated via the
canonical seesaw mechanism with three heavy sterile neutrinos, the 3\times 3
PontecorvoMakiNakagawaSakata neutrino mixing matrix V will not be exactly
unitary. This kind of indirect unitarity violation can be probed in a precision
reactor antineutrino oscillation experiment, but it may be entangled with
terrestrial matter effects as both of them are very small. We calculate the
probability of \overline{\nu}_e \to \overline{\nu}_e oscillations in a good
analytical approximation, and find that, besides the zerodistance effect, the
effect of unitarity violation is always smaller than matter effects, and their
entanglement does not appear until the nexttoleadingorder oscillating terms
are taken into account. Given a 20kiloton JUNOlike liquid scintillator
detector, we reaffirm that terrestrial matter effects should not be neglected
but indirect unitarity violation makes no difference, and demonstrate that the
experimental sensitivities to the neutrino mass ordering and a precision
measurement of \theta_{12} and \Delta_{21} \equiv m^2_2  m^2_1 are robust.

Given the Dirac neutrino mass term, we explore the constraint conditions
which allow the corresponding mass matrix to be invariant under the \mu\tau
reflection transformation, leading us to the phenomenologically favored
predictions \theta_{23} = \pi/4 and \delta = 3\pi/2 in the standard
parametrization of the 3\times 3 lepton flavor mixing matrix. If such a flavor
symmetry is realized at a superhigh energy scale \Lambda_{\mu\tau}, we
investigate how it is spontaneously broken via the oneloop
renormalizationgroup equations (RGEs) running from \Lambda_{\mu\tau} down to
the Fermi scale \Lambda_{\rm F}. Such quantum corrections to the neutrino
masses and flavor mixing parameters are derived, and an analytical link is
established between the Jarlskog invariants of CP violation at
\Lambda_{\mu\tau} and \Lambda_{\rm F}. Some numerical examples are also
presented in both the minimal supersymmetric standard model and the typeII
twoHiggsdoublet model, to illustrate how the octant of \theta_{23}, the
quadrant of \delta and the neutrino mass ordering are correlated with one
another as a result of the RGEinduced \mu\tau reflection symmetry breaking
effects.

If the neutrino mass spectrum turns out to be m^{}_3 < m^{}_1 < m^{}_2, one
may choose to relabel it as m^{\prime}_1 < m^{\prime}_2 < m^{\prime}_3 such
that all the masses of fundamental fermions with the same electrical charges
are in order. In this case the columns of the 3\times 3 lepton flavor mixing
matrix U should be reordered accordingly, and the resulting pattern U^\prime
may involve one or two large mixing angles in the standard parametrization or
its variations. Since the Majorana neutrino mass matrix keeps unchanged in such
a mass relabeling, a possible \mu\tau reflection symmetry is respected in this
connection and its breaking effects are modelindependently constrained at the
3\sigma level by using current experimental data.

Constraining neutrino mass remains an elusive challenge in modern physics.
Precision measurements are expected from several upcoming cosmological probes
of largescale structure. Achieving this goal relies on an equal level of
precision from theoretical predictions of neutrino clustering. Numerical
simulations of the nonlinear evolution of cold dark matter and neutrinos play
a pivotal role in this process. We incorporate neutrinos into the cosmological
Nbody code CUBEP3M and discuss the challenges associated with pushing to the
extreme scales demanded by the neutrino problem. We highlight code
optimizations made to exploit modern high performance computing architectures
and present a novel method of data compression that reduces the phasespace
particle footprint from 24 bytes in single precision to roughly 9 bytes. We
scale the neutrino problem to the Tianhe2 supercomputer and provide details of
our production run, named TianNu, which uses 86% of the machine (13,824 compute
nodes). With a total of 2.97 trillion particles, TianNu is currently the
world's largest cosmological Nbody simulation and improves upon previous
neutrino simulations by two orders of magnitude in scale. We finish with a
discussion of the unanticipated computational challenges that were encountered
during the TianNu runtime.

The democracy of quark flavors is a wellmotivated flavor symmetry, but it
must be properly broken in order to explain the observed quark mass spectrum
and flavor mixing pattern. We reconstruct the texture of flavor democracy
breaking and evaluate its strength in a novel way, by assuming a parallelism
between the Q=+2/3 and Q=1/3 quark sectors and using a nontrivial
parametrization of the flavor mixing matrix. Some phenomenological implications
of such democratic quark mass matrices, including their variations in the
hierarchy basis and their evolution from the electroweak scale to a
superhighenergy scale, are also discussed.

If massive neutrinos are the Majorana particles and have a normal mass
ordering, the effective mass term $\langle m\rangle^{}_{ee}$ of a neutrinoless
doublebeta ($0\nu 2\beta$) decay may suffer significant cancellations among
its three components and thus sink into a decline, resulting in a "well" in the
threedimensional graph of $\langle m\rangle^{}_{ee}$ against the smallest
neutrino mass $m^{}_1$ and the relevant Majorana phase $\rho$. We present a new
and complete analytical understanding of the fine issues inside such a well,
and discover a novel threshold of $\langle m\rangle^{}_{ee}$ in terms of the
neutrino masses and flavor mixing angles: $\langle m\rangle^{}_{ee}^{}_* =
m^{}_3 \sin^2\theta^{}_{13}$ in connection with $\tan\theta^{}_{12} =
\sqrt{m^{}_1/m^{}_2}$ and $\rho =\pi$. This threshold point, which links the
{\it local} minimum and maximum of $\langle m\rangle^{}_{ee}$, can be used to
signify observability or sensitivity of the future $0\nu 2\beta$decay
experiments. Given current neutrino oscillation data, the possibility of
$\langle m\rangle^{}_{ee} < \langle m\rangle^{}_{ee}^{}_*$ is found to be
very small.

Measurements of neutrino mass in cosmological observations rely on two point
statistics that are hindered by significant degeneracies with the optical depth
and galaxy bias. The relative velocity effect between cold dark matter and
neutrinos induces a large scale dipole into the matter density field and may be
able to provide orthogonal constraints to standard techniques. We numerically
investigate this dipole in the TianNu Simulation, which contains cold dark
matter and 50 meV neutrinos. We first compute the dipole using a new linear
response technique where we treat the displacement caused by the relative
velocity as a phase in Fourier space and then integrate the matter power
spectrum over redshift. Then, we compute the dipole numerically in real space
using the simulation density and velocity fields. We find excellent agreement
between the linear response and Nbody methods. Utilizing the dipole as an
observational tool will require two tracers of the matter distribution that are
differently biased with respect to the neutrino density.

Astrophysical techniques have pioneered the discovery of neutrino mass
properties. Current cosmological observations give an upper bound on neutrino
masses by attempting to disentangle the small neutrino contribution from the
sum of all matter using precise theoretical models. We discover the
differential neutrino condensation effect in our TianNu Nbody simulation.
Neutrino masses can be inferred using this effect by comparing galaxy
properties in regions of the universe with different neutrino relative
abundance (i.e. the local neutrino to cold dark matter density ratio). In
"neutrinorich"' regions, more neutrinos can be captured by massive halos
compared to "neutrinopoor" regions. This effect differentially skews the halo
mass function and opens up the path to independent neutrino mass measurements
in current or future galaxy surveys.

Given an acceleratorbased neutrino experiment with the beam energy E
\lesssim 1 GeV, we expand the probabilities of \nu_\mu \to \nu_e and \overline
{\nu}_\mu \to \overline {\nu}_e oscillations in matter in terms of two small
quantities \Delta_{21}/\Delta_{31} and A/\Delta_{31}, where \Delta_{21} \equiv
m^2_2  m^2_1 and \Delta_{31} \equiv m^2_3  m^2_1 are the neutrino
masssquared differences, and A measures the strength of terrestrial matter
effects. Our analytical approximations are numerically more accurate than those
made by Freund in this energy region, and thus they are particularly applicable
for the study of leptonic CP violation in the lowenergy MOMENT, ESS\nuSM and
T2K oscillation experiments. As a byproduct, the new analytical approximations
help us to easily understand why the mattercorrected Jarlskog parameter
\widetilde{\cal J} peaks at the resonance energy E_* \simeq 0.14 GeV (or 0.12
GeV) for the normal (or inverted) neutrino mass hierarchy, and how the three
Dirac unitarity triangles are deformed due to the terrestrial matter
contamination. We also affirm that a mediumbaseline neutrino oscillation
experiment with the beam energy E lying in the E_* \lesssim E \lesssim 2 E_*
range is capable of exploring leptonic CP violation with little matterinduced
suppression.

We have carefully examined, in both analytical and numerical ways, how small
the terrestrial matter effects can be in a given mediumbaseline reactor
antineutrino oscillation experiment like JUNO or RENO50. Taking the ongoing
JUNO experiment for example, we show that the inclusion of terrestrial matter
effects may reduce the sensitivity of the neutrino mass ordering measurement by
\Delta \chi^2_{\rm MO} \simeq 0.6, and a neglect of such effects may shift the
bestfit values of the flavor mixing angle \theta_{12} and the neutrino
masssquared difference \Delta_{21} by about 1\sigma to 2\sigma in the future
data analysis. In addition, a preliminary estimate indicates that a 2\sigma
sensitivity of establishing the terrestrial matter effects can be achieved for
about 10 years of data taking at JUNO with the help of a proper near detector
implementation.

Behind the observed pattern of lepton flavor mixing is a partial or
approximate mutau flavor symmetry  a milestone on our road to the true
origin of neutrino masses and flavor structures. In this review article we
first describe the features of mutau permutation and reflection symmetries,
and then explore their various consequences on model building and neutrino
phenomenology. We pay particular attention to soft mutau symmetry breaking,
which is crucial for our deeper understanding of the fine effects of flavor
mixing and CP violation.

Given the bestfit results of six neutrino oscillation parameters, we plot
the Dirac and Majorana unitarity triangles (UTs) of the 3\times 3 lepton flavor
mixing matrix to show their real shapes for the first time. The connections of
the Majorana UTs with neutrinoantineutrino oscillations and neutrino decays
are explored, and the possibilities of right or isosceles UTs are discussed. In
the neutrino mass limit of m_1 \to 0 or m_3 \to 0, which is allowed by current
experimental data, we show how the six triangles formed by the effective
Majorana neutrino masses \langle m\rangle_{\alpha\beta} (for \alpha, \beta = e,
\mu, \tau) and their corresponding component vectors look like in the complex
plane. The relations of such triangles to the Majorana phases and to the
leptonnumberviolating decays H^{++} \to \alpha^+ \beta^+ in the typeII
seesaw mechanism are also illustrated.

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.

We carry out a new study of quark mass matrices $M^{}_{\rm u}$ (uptype) and
$M^{}_{\rm d}$ (downtype) which are Hermitian and have four zero entries, and
find a new part of the parameter space which was missed in the previous works.
We identify two more specific fourzero patterns of $M^{}_{\rm u}$ and
$M^{}_{\rm d}$ with fewer free parameters, and present two toy flavorsymmetry
models which can help realize such special and interesting quark flavor
structures. We also show that the texture zeros of $M^{}_{\rm u}$ and
$M^{}_{\rm d}$ are essentially stable against the evolution of energy scales in
an analytical way by using the oneloop renormalizationgroup equations.

The Majorana nature of massive neutrinos will be crucially probed in the
nextgeneration experiments of the neutrinoless doublebeta ($0\nu 2\beta$)
decay. The effective mass term of this process, $\langle m\rangle^{}_{ee}$, may
be contaminated by new physics. So how to interpret a discovery or null result
of the $0\nu 2\beta$ decay in the foreseeable future is highly nontrivial. In
this paper we introduce a novel threedimensional description of $\langle
m\rangle_{ee}^{}$, which allows us to see its sensitivity to the lightest
neutrino mass and two Majorana phases in a transparent way. We take a look at
to what extent the free parameters of $\langle m\rangle_{ee}^{}$ can be well
constrained provided a signal of the $0\nu 2\beta$ decay is observed someday.
To fully explore lepton number violation, all the six effective Majorana mass
terms $\langle m\rangle_{\alpha\beta}^{}$ (for $\alpha, \beta = e, \mu, \tau$)
are calculated and their lower bounds are illustrated with the twodimensional
contour figures. The effect of possible new physics on the $0\nu 2\beta$ decay
is also discussed in a modelindependent way. We find that the result of
$\langle m\rangle_{ee}^{}$ in the normal (or inverted) neutrino mass ordering
case modified by the new physics effect may somewhat mimic that in the inverted
(or normal) mass ordering case in the standard threeflavor scheme. Hence a
proper interpretation of a discovery or null result of the $0\nu 2\beta$ decay
may demand extra information from some other measurements.

This essay is intended to provide a brief description of the peculiar
properties of neutrinos within and beyond the standard theory of weak
interactions. The focus is on the flavor oscillations of massive neutrinos,
from which one has achieved some striking knowledge about their mass spectrum
and flavor mixing pattern. The experimental prospects towards probing the
absolute neutrino mass scale, possible Majorana nature and CPviolating effects
will also be addressed.

The fact that quarks of the same electric charge possess a mass hierarchy is
a big puzzle in particle physics, and it must be highly correlated with the
hierarchy of quark flavor mixing. This review article is intended to provide a
brief description of some important issues regarding quark masses, flavor
mixing and CP violation. A comparison between the salient features of quark and
lepton flavor mixing structures is also made.

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.

We find that the lepton flavor mixing matrix $U$ should possess a partial
$\mu$$\tau$ permutation symmetry $U^{}_{\mu 1} = U^{}_{\tau 1}$, and the
latter predicts a novel correlation between the Dirac CPviolating phase
$\delta$ and three flavor mixing angles $\theta^{}_{12}$, $\theta^{}_{13}$ and
$\theta^{}_{23}$ in the standard parametrization. Inputting the bestfit values
of these angles reported by Capozzi {\it et al}, we obtain the prediction
$\delta \simeq 255^\circ$ in the normal neutrino mass ordering, which is in
good agreement with the bestfit result $\delta \simeq 250^\circ$. In this
connection the inverted neutrino mass ordering is slightly disfavored. If this
partial $\mu$$\tau$ symmetry is specified to be $U^{}_{\mu 1} = U^{}_{\tau
1} =1/\sqrt{6}~$, one can reproduce the phenomenologicallyfavored relation
$\sin^2\theta^{}_{12} = \left(1  2\tan^2\theta^{}_{13}\right)/3$ and a viable
twoparameter description of $U$ which were first uncovered in 2006. Moreover,
we point out that the octant of $\theta^{}_{23}$ and the quadrant of $\delta$
can be resolved thanks to the slight violation of $U^{}_{\mu 2} = U^{}_{\tau
2}$ and $U^{}_{\mu 3} = U^{}_{\tau 3}$ either at the tree level or from
radiative corrections.

I give a theoretical overview of some basic properties of massive neutrinos
in these lectures. Particular attention is paid to the origin of neutrino
masses, the pattern of lepton flavor mixing, the feature of leptonic CP
violation and the electromagnetic properties of massive neutrinos. I highlight
the TeV seesaw mechanisms as a possible bridge between neutrino physics and
collider physics in the era characterized by the Large Hadron Collider.

The neutrinoless doublebeta ($0\nu\beta\beta$) decay is a unique process to
identify the Majorana nature of massive neutrinos, and its rate depends on the
size of the effective Majorana neutrino mass $\langle m\rangle_{ee}$. We put
forward a novel "couplingrod" diagram to describe $\langle m\rangle_{ee}$ in
the complex plane, by which the effects of the neutrino mass ordering and
CPviolating phases on $\langle m\rangle_{ee}$ are intuitively understood. We
show that this geometric language allows us to easily obtain the maximum and
minimum of $\langle m\rangle_{ee}$. It remains usable even if there is a kind
of new physics contributing to $\langle m\rangle_{ee}$, and it can also be
extended to describe the effective Majorana masses $\langle m\rangle_{e\mu}$,
$\langle m\rangle_{e\tau}$, $\langle m\rangle_{\mu\mu}$, $\langle
m\rangle_{\mu\tau}$ and $\langle m\rangle_{\tau\tau}$ which may appear in some
other leptonnumberviolating processes.

The origin of tiny neutrino masses and the baryon number asymmetry of the
Universe are naturally interpreted by the canonical seesaw and leptogenesis
mechanisms, in which there are the heavy Majorana neutrino mass matrix $M_{\rm
R}$, the Dirac neutrino mass matrix $M_{\rm D}$, the chargedlepton mass matrix
$M_\ell$ and the effective (light) neutrino mass matrix $M_\nu$. We find that
${\rm Im}\left(\det\left[ M^\dagger_{\rm D} M_{\rm D}, M^\dagger_{\rm R} M_{\rm
R} \right]\right)$, ${\rm Im}\left(\det\left[ M_\ell M^\dagger_\ell, M_\nu
M^\dagger_\nu \right]\right)$ and ${\rm Im}\left(\det\left[ M_\ell
M^\dagger_\ell, M_{\rm D} M^\dagger_{\rm D} \right]\right)$ can serve for a
basisindependent measure of CP violation associated with
leptonnumberviolating decays of heavy neutrinos, flavor oscillations of light
neutrinos and leptonflavorviolating decays of charged leptons, respectively.
We first calculate these quantities with the help of a standard parametrization
of the $6\times 6$ flavor mixing matrix, and then discuss their implications on
both leptogenesis and CP violation at low energy scales. A comparison with the
weakbasis invariants of leptogenesis as proposed by Branco {\it et al} is also
made.