
The standard model of cosmology predicts the existence of cosmic neutrino
background in the present Universe. To detect cosmic relic neutrinos in the
vicinity of the Earth, it is necessary to evaluate the gravitational clustering
effects on relic neutrinos in the Milky Way. Here we introduce a reweighting
technique in the Nonebody simulation method, so that a single simulation can
yield neutrino density profiles for different neutrino masses and phase space
distributions. In light of current experimental results that favor small
neutrino masses, the neutrino number density contrast around the Earth is found
to be almost proportional to the square of neutrino mass. The density
contrastmass relation and the reweighting technique are useful for studying
the phenomenology associated with the future detection of the cosmic neutrino
background.

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.

Motivated by tremendous progress in neutrino oscillation experiments, we
derive a new set of simple and compact formulas for threeflavor neutrino
oscillation probabilities in matter of a constant density. A useful definition
of the $\eta$gauge neutrino masssquared difference $\Delta^{}_* \equiv \eta
\Delta^{}_{31} + (1\eta) \Delta^{}_{32}$ is introduced, where $\Delta^{}_{ji}
\equiv m^2_j  m^2_i$ for $ji = 21, 31, 32$ are the ordinary neutrino
masssquared differences and $0 \leq \eta \leq 1$ is a real and positive
parameter. Expanding neutrino oscillation probabilities in terms of $\alpha
\equiv \Delta^{}_{21}/\Delta^{}_*$, we demonstrate that the analytical formulas
can be remarkably simplified for $\eta = \cos^2 \theta^{}_{12}$, with
$\theta_{12}^{}$ being the solar mixing angle. As a byproduct, the mapping
from neutrino oscillation parameters in vacuum to their counterparts in matter
is obtained at the order of ${\cal O}(\alpha^2)$. Finally, we show that our
approximate formulas are not only valid for an arbitrary neutrino energy and
any baseline length, but also still maintaining a high level of accuracy.

The Littlest Seesaw (LS) model involves two righthanded neutrinos and a very
constrained Dirac neutrino mass matrix, involving one texture zero and two
independent Dirac masses, leading to a highly predictive scheme in which all
neutrino masses and the entire PMNS matrix is successfully predicted in terms
of just two real parameters. We calculate the renormalisation group (RG)
corrections to the LS predictions, with and without supersymmetry, including
also the threshold effects induced by the decoupling of heavy Majorana
neutrinos both analytically and numerically. We find that the predictions for
neutrino mixing angles and mass ratios are rather stable under RG corrections.
For example we find that the LS model with RG corrections predicts close to
maximal atmospheric mixing, $\theta_{23}=45^\circ \pm 1^\circ$, in most
considered cases, in tension with the latest NOvA results. The techniques used
here apply to other seesaw models with a strong normal mass hierarchy.

In light of the latest neutrino oscillation data, we examine whether the
leptonic flavor mixing matrix can take on an exact form of tribimaximal (TBM),
goldenratio (GR) or bimaximal (BM) mixing pattern at a superhighenergy scale,
where such a mixing pattern could be realized by a flavor symmetry, and become
compatible with experimental data at the lowenergy scale. Within the framework
of the Minimal Supersymmetric Standard Model (MSSM), the only hope for
realizing such a possibility is to count on the corrections from the
renomalizationgroup (RG) running. In this work we focus on these radiative
corrections, and fully explore the allowed parameter space for each of these
mixing patterns. We find that when the upper bound on the sum of neutrino
masses $\Sigma^{}_\nu \equiv m^{}_1 + m^{}_2 + m^{}_3 < 0.23~\text{eV}$ at the
$95\%$ confidence level from Planck 2015 is taken into account, none of these
mixing patterns can be identified as the leptonic mixing matrix below the
seesaw threshold. If this cosmological upper bound on the sum of neutrino
masses were relaxed, the TBM and GR mixing patterns would still be compatible
with the latest neutrino oscillation data at the $3\sigma$ level, but not at
the $1\sigma$ level. Even in this case, no such a possibility exists for the BM
mixing.

Parametric instability is an intrinsic risk in high power laser
interferometer gravitational wave detectors, in which the optical cavity modes
interact with the acoustic modes of the mirrors leading to exponential growth
of the acoustic vibration. In this paper, we investigate the potential
parametric instability for a proposed next generation gravitational wave
detector based on cooled silicon test masses. It is shown that there would be
about 2 unstable modes per test mass, with the highest parametric gain of ~76.
The importance of developing suitable instability suppression schemes is
emphasized.

The simple correlation among three lepton flavor mixing angles
$(\theta^{}_{12}, \theta^{}_{13}, \theta^{}_{23})$ and the leptonic Dirac
CPviolating phase $\delta$ is conventionally called a sum rule of lepton
flavor mixing, which may be derived from a class of neutrino mass models with
flavor symmetries. In this paper, we consider the solar lepton mixing sum rule
$\theta^{}_{12} \approx \theta^{\nu}_{12} + \theta^{}_{13} \cos \delta$, where
$\theta^\nu_{12}$ stems from a constant mixing pattern in the neutrino sector
and takes the value of $\theta^\nu_{12} = 45^\circ$ for the bimaximal mixing
(BM), $\theta^\nu_{12} = \tan^{1}(1/\sqrt{2}) \approx 35.3^\circ$ for the
tribimaximal mixing (TBM) or $\theta^\nu_{12} = \tan^{1}\left[2/(\sqrt{5} +
1)\right] \approx 31.7^\circ$ for the goldenratio mixing (GR), and investigate
the renormalizationgroup (RG) running effects on lepton flavor mixing
parameters when this sum rule is assumed at a superhighenergy scale. For
illustration, we work within the framework of the minimal supersymmetric
standard model (MSSM), and implement the Bayesian approach to explore the
posterior distribution of $\delta$ at the lowenergy scale, which becomes quite
broad when the RG running effects are significant. Moreover, we also discuss
the compatibility of the above three mixing scenarios with current neutrino
oscillation data, and observe that radiative corrections can increase such a
compatibility for the BM scenario, resulting in a weaker preference for the TBM
and GR ones.

The SchechterValle theorem states that a positive observation of
neutrinoless doublebeta ($0\nu \beta \beta$) decays implies a finite Majorana
mass term for neutrinos when any unlikely finetuning or cancellation is
absent. In this note, we reexamine the quantitative impact of the
SchechterValle theorem, and find that current experimental lower limits on the
halflives of $0\nu \beta \beta$decaying nuclei have placed a restrictive
upper bound on the Majorana neutrino mass $\delta m^{ee}_\nu < 7.43 \times
10^{29}~{\rm eV}$ radiatively generated at the fourloop level. Furthermore,
we generalize this quantitative analysis of $0\nu \beta \beta$ decays to that
of the leptonnumberviolating (LNV) meson decays $M^ \to {M^\prime}^+ +
\ell^_\alpha + \ell^_\beta$ (for $\alpha$, $\beta$ = $e$ or $\mu$). Given the
present upper limits on these rare LNV decays, we have derived the loopinduced
Majorana neutrino masses $\delta m^{ee}_\nu < 9.7 \times 10^{18}~{\rm eV}$,
$\delta m^{e\mu}_\nu < 1.6 \times 10^{15}~{\rm eV}$ and $\delta m^{\mu
\mu}_\nu < 1.0 \times 10^{12}~{\rm eV}$ from $K^ \to \pi^+ + e^ + e^$,
$K^ \to \pi^+ + e^ + \mu^$ and $K^ \to \pi^+ + \mu^ + \mu^$,
respectively. A partial list of radiative neutrino masses from the LNV decays
of $D$, $D_s^{}$ and $B$ mesons is also given.

The atomic transition from an excited state ${\rm e}\rangle$ to the ground
state ${\rm g}\rangle$ by emitting a neutrino pair and a photon, i.e., ${\rm
e}\rangle \to {\rm g}\rangle + \gamma\rangle + \nu^{}_i\rangle +
\bar{\nu}^{}_j\rangle$ with $i, j = 1, 2, 3$, has been proposed by Yoshimura
and his collaborators as an alternative way to determine the absolute scale
$m^{}_0$ of neutrino masses. More recently, a statistical analysis of the fine
structure of the photon spectrum from this atomic process has been performed
[N. Song {\it et al.}, Phys.\ Rev.\ D {\bf 93}, 013020 (2016)] to
quantitatively examine the experimental requirements for a realistic
determination of absolute neutrino masses. In this paper, we show how to
improve the statistical analysis and demonstrate that the previously required
detection time can be reduced by one order of magnitude for the case of a
$3\sigma$ determination of $m_0^{} \sim 0.01~\text{eV}$ with an accuracy better
than $10\%$. Such an improvement is very encouraging for further investigations
on measuring absolute neutrino masses through atomic processes.

We study the purely leptonic decays of $W^\pm \to e^\pm e^\pm \mu^\mp \nu$
and $\mu^\pm \mu^\pm e^\mp \nu$ produced at the LHC, induced by sterile
neutrinos with mass $m_N$ below $M_W$ in the intermediate state. Since the
final state neutrino escapes detection, one cannot tell whether this process
violates lepton number, what would indicate a Majorana character for the
intermediate sterile neutrino. Our study shows that when the sterile neutrino
mixings with electrons and muons are different enough, one can still
discriminate between the Dirac and Majorana character of this intermediate
neutrino by simply counting and comparing the above decay rates. After
performing collider simulations and statistical analysis, we find that at the
$14~\text{TeV}$ LHC with an integrated luminosity of $3000~\text{fb}^{1}$, for
two benchmark scenarios $m_N$ = 20 GeV and 50 GeV, at least a $3\sigma$ level
of exclusion on the Dirac case can be achieved for disparities as mild as e.g.
$U_{Ne}^2 < 0.7~ U_{N\mu}^2$ or $U_{N\mu}^2 < 0.7~ U_{N e}^2$, provided
that $U_{Ne}^2$, $U_{N\mu}^2$ are both above $\sim 2\times 10^{6}$.

Motivated by recent intensive experimental efforts on searching for
neutrinoless double beta decays, we present a detailed quantitative analysis on
the prospect of resolving neutrino mass ordering in the next generation
$^{76}$Getype experiments.

Recently, both ATLAS and CMS collaborations at the CERN Large Hadron Collider
(LHC) have announced their observations of an excess of diphoton events around
the invariant mass of $750~{\rm GeV}$ with a local significance of $3.6\sigma$
and $2.6\sigma$, respectively. In this paper, we interpret the diphoton excess
as the onshell production of a real singlet scalar in the $pp \to S \to \gamma
\gamma$ channel. To accommodate the observed production rate, we further
introduce a vectorlike fermion $F$, which is carrying both color and electric
charges. The viable regions of model parameters are explored for this simple
extension of the Standard Model (SM). Moreover, we revisit the problem of
electroweak vacuum stability in the same scenario, and find that the
requirement for the electroweak vacuum stability up to high energy scales
imposes serious constraints on the Yukawa coupling of the vectorlike fermion
and the quartic couplings of the SM Higgs boson and the new singlet scalar.
Consequently, a successful explanation for the diphoton excess and the absolute
stability of electroweak vacuum cannot be achieved simultaneously in this
economical setup.

It remains to be determined experimentally if massive neutrinos are Majorana
or Dirac particles. In this connection, it has been recently suggested that the
detection of cosmic neutrino background of lefthanded neutrinos $\nu^{}_{\rm
L}$ and righthanded antineutrinos $\overline{\nu}^{}_{\rm R}$ in future
experiments of neutrino capture on betadecaying nuclei (e.g., $\nu^{}_e +
{^3{\rm H}} \to {^3}{\rm He} + e^$ for the PTOLEMY experiment) is likely to
distinguish between Majorana and Dirac neutrinos, since the capture rate is
twice larger in the former case. In this paper, we investigate the possible
impact of righthanded neutrinos on the capture rate, assuming that massive
neutrinos are Dirac particles and both righthanded neutrinos $\nu^{}_{\rm R}$
and lefthanded antineutrinos $\overline{\nu}^{}_{\rm L}$ can be efficiently
produced in the early Universe. It turns out that the capture rate can be
enhanced at most by $28\%$ due to the presence of relic $\nu^{}_{\rm R}$ and
$\overline{\nu}^{}_{\rm L}$ with a total number density of $95~{\rm cm}^{3}$,
which should be compared to the number density $336~{\rm cm}^{3}$ of cosmic
neutrino background. The enhancement has actually been limited by the latest
cosmological and astrophysical bounds on the effective number of neutrino
generations $N^{}_{\rm eff} = 3.14^{+0.44}_{0.43}$ at the $95\%$ confidence
level. For illustration, two possible scenarios have been proposed for thermal
production of righthanded neutrinos in the early Universe.

Motivated by recent intensive experimental efforts on searching for
neutrinoless doublebeta decays, we perform a detailed analysis of the physics
potential of the experiments based on $^{76}\mathrm{Ge}$. Assuming no signals,
current and future experiments could place a $90\%$ lower limit on the half
life $T^{0\nu}_{1/2} \gtrsim 4\times 10^{26}~{\rm yr}$ and $T^{0\nu}_{1/2}
\gtrsim 7\times 10^{27}~{\rm yr}$, respectively. Then, how to report an
evidence for neutrinoless doublebeta decays is addressed by following the
Bayesian statistical approach. For the first time, we present a quantitative
description of experimental power to distinguish between normal and inverted
neutrino mass orderings. Taking an exposure of $10^{4}~{\rm kg}\cdot{\rm yr}$
and a background rate of $10^{4}~{\rm counts}/({\rm keV}\cdot{\rm kg}\cdot{\rm
yr})$, we find that a moderate evidence for normal neutrino mass ordering
(i.e., with a Bayes factor ${\cal B}$ given by $\ln({\cal B}) \simeq 2.5$ or a
probability about $92.3\%$ according to the Jeffreys scale) can be achieved if
the true value of effective neutrino mass $m^{}_{\beta\beta}$ turns out to be
below $0.01~{\rm eV}$.

We consider a restricted TypeI seesaw scenario with four texture zeros in
the neutrino Yukawa matrix, in the weak basis where both the chargedlepton
Yukawa matrix and the Majorana mass matrix for righthanded neutrinos are
diagonal and real. Inspired by grand unified theories, we further require the
neutrino Yukawa matrix to exhibit a similar hierarchical pattern to that in the
uptype quark Yukawa matrix. With such a hierarchy requirement, we find that
leptogenesis, which would operate in a $N_2^{}$dominated scenario with the
asymmetry generated by the nexttolightest righthanded neutrino $N_2^{}$, can
greatly reduce the number of allowed textures, and disfavors the scenario that
three light neutrinos are quasidegenerate. Such a quasidegenerate scenario of
light neutrinos may soon be tested in upcoming neutrino experiments.

In light of the latest neutrino oscillation data, we revisit the minimal
scenario of typeI seesaw model, in which only two heavy righthanded Majorana
neutrinos are introduced to account for both tiny neutrino masses and the
baryon number asymmetry in our Universe. In this framework, we carry out a
systematic study of the FramptonGlashowYanagida ansatz by taking into account
the renormalizationgroup running of neutrino mixing parameters and the flavor
effects in leptogenesis. We demonstrate that the normal neutrino mass ordering
is disfavored even in the minimal supersymmetric standard model with a large
value of $\tan \beta$, for which the running effects could be significant.
Furthermore, it is pointed out that the original scenario with a hierarchical
mass spectrum of heavy Majorana neutrinos contradicts with the upper bound
derived from a naturalness criterion, and the resonant mechanism with
nearlydegenerate heavy Majorana neutrinos can be a possible way out.

We present a model in which the Supersymmetric Standard Model is augmented by
the family symmetry $\bs{\m Z_7 \rtimes \m Z_3}$. Motivated by $SO(10)$, where
the charge twothirds and neutral Dirac Yukawa matrices are related, we
propose, using family symmetry, a special form for the seesaw Majorana matrix;
it contains a squared correlated hierarchy, allowing it to mitigate the severe
hierarchy of the quark sector. It is reproduced naturally by the invariant
operators of $\bs{\m Z_7 \rtimes \m Z_3}$, with the hierarchy carried by
familon fields. In addition to relating the hierarchy of the $\Delta I_{\rm
w}=1/2$ to the $\Delta I_{\rm w}=0$ sector, it contains a GattoSartoriTonin
like relation, predicts a normal hierarchy for Tribimaximal and Golden Ratio
mixings, and gives specific values for the light neutrino masses.

We present the results of a numerical search for the Dirac Yukawa matrices of
the Standard Model, consistent with the quark and lepton masses and their
mixing angles. We assume a diagonal upquark matrix, natural in $\bs{\m Z_7
\rtimes \m Z_3}$, Bimaximal or Tribimaximal seesaw mixing, and $SU(5)$
unification to relate the downquark and charged lepton Dirac Yukawa matrices
using GeorgiJarlskog mechanisms. The measured value of $\theta_{13}$ requires
an asymmetric downquark Yukawa matrix. Satisfying the measured values of both
$\theta_{13}$ and the electron mass restricts the number of solutions,
underlying the importance of the recent measurement of the reactor angle.

We present a toy model in which the Higgs sector fields transform as
nonAbelian representations of a family symmetry group, and consider the
possibility that the extra family partners of the Higgs particles act as
messengers for both supersymmetry and family symmetry breakings. Although such
mediation schemes generically produce family dependent soft supersymmetry
breaking terms at the messenger scale, we demonstrate the existence of a
focusing mechanism which may erase such hierarchies through renormalization
group running.

In recent years, many efforts have been addressed on collision avoidance of
collectively moving agents. In this paper, we propose a modified version of the
Vicsek model with adaptive speed, which can guarantee the absence of
collisions. However, this strategy leads to an aggregated state with slowly
moving agents. We therefore further introduce a certain repulsion, which
results in both faster consensus and longer safe distance among agents, and
thus provides a powerful mechanism for collective motions in biological and
technological multiagent systems.

In resent years, Vicsek model has attracted more and more attention and been
well developed. However, the indepth analysis on the convergence time are
scarce thus far. In this paper, we study some certain factors that mainly
govern the convergence time of Vicsek model. By extensively numerical
simulations, we find the convergence time scales in a power law with $r^2\ln N$
in the noisefree case, where $r$ and $N$ are horizon radius and the number of
particles. Furthermore, to accelerate the convergence, we propose a new model
in which the speed of each particle is variable. The convergence time can be
remarkably shortened compared with the standard Vicsek model.