• 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 N-one-body 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 contrast-mass 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 right-handed neutrinos are added to the Standard Model of particle physics. Assuming that both the left- and right-handed 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 right-handed 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 three-flavor neutrino oscillation probabilities in matter of a constant density. A useful definition of the $\eta$-gauge neutrino mass-squared 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 mass-squared 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 by-product, 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 right-handed 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 tri-bimaximal (TBM), golden-ratio (GR) or bimaximal (BM) mixing pattern at a superhigh-energy scale, where such a mixing pattern could be realized by a flavor symmetry, and become compatible with experimental data at the low-energy 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 renomalization-group (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 CP-violating 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 bi-maximal mixing (BM), $\theta^\nu_{12} = \tan^{-1}(1/\sqrt{2}) \approx 35.3^\circ$ for the tri-bimaximal mixing (TBM) or $\theta^\nu_{12} = \tan^{-1}\left[2/(\sqrt{5} + 1)\right] \approx 31.7^\circ$ for the golden-ratio mixing (GR), and investigate the renormalization-group (RG) running effects on lepton flavor mixing parameters when this sum rule is assumed at a superhigh-energy 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 low-energy 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 Schechter-Valle theorem states that a positive observation of neutrinoless double-beta ($0\nu \beta \beta$) decays implies a finite Majorana mass term for neutrinos when any unlikely fine-tuning or cancellation is absent. In this note, we reexamine the quantitative impact of the Schechter-Valle theorem, and find that current experimental lower limits on the half-lives 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 four-loop level. Furthermore, we generalize this quantitative analysis of $0\nu \beta \beta$ decays to that of the lepton-number-violating (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 loop-induced 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}$Ge-type 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 on-shell 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 vector-like 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 vector-like 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 left-handed neutrinos $\nu^{}_{\rm L}$ and right-handed antineutrinos $\overline{\nu}^{}_{\rm R}$ in future experiments of neutrino capture on beta-decaying 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 right-handed neutrinos on the capture rate, assuming that massive neutrinos are Dirac particles and both right-handed neutrinos $\nu^{}_{\rm R}$ and left-handed 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 right-handed neutrinos in the early Universe.
  • Motivated by recent intensive experimental efforts on searching for neutrinoless double-beta 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 double-beta 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 Type-I seesaw scenario with four texture zeros in the neutrino Yukawa matrix, in the weak basis where both the charged-lepton Yukawa matrix and the Majorana mass matrix for right-handed 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 up-type 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 next-to-lightest right-handed neutrino $N_2^{}$, can greatly reduce the number of allowed textures, and disfavors the scenario that three light neutrinos are quasi-degenerate. Such a quasi-degenerate 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 type-I seesaw model, in which only two heavy right-handed 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 Frampton-Glashow-Yanagida ansatz by taking into account the renormalization-group 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 nearly-degenerate 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 two-thirds 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 Gatto-Sartori-Tonin like relation, predicts a normal hierarchy for Tri-bimaximal 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 up-quark matrix, natural in $\bs{\m Z_7 \rtimes \m Z_3}$, Bimaximal or Tri-bimaximal seesaw mixing, and $SU(5)$ unification to relate the down-quark and charged lepton Dirac Yukawa matrices using Georgi-Jarlskog mechanisms. The measured value of $\theta_{13}$ requires an asymmetric down-quark 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 non-Abelian 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 multi-agent systems.
  • In resent years, Vicsek model has attracted more and more attention and been well developed. However, the in-depth 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 noise-free 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.