
We discuss the minimal seesaw model for the Dirac CP violating phase of the
lepton mixing matrix. We introduce two righthanded Majorana neutrinos and
obtain several textures of the trimaximal lepton mixing matrices. Moreover, we
discuss the observed baryon asymmetry of the universe through the leptogenesis
mechanism. As the result, we obtain the specific model which predicts the
negative sign of maximal Dirac CP violating phase and normal hierarchy of
neutrino masses.

We revisit the $A_4$ model for leptons in light of new result of NuFIT 3.2.
We introduce a new flavon $\eta$ transforming as $A_4$ singlet $1'$ or $1"$
which couples to both charged leptons and neutrinos in nextleading order
operators. The model consists of the five parameters: the lightest neutrino
mass $m_1$, the vacuum expectation value of $\eta$ and three CP violating
phases after inputting the experimental values of $\Delta m_{\rm atm}^2$ and
$\Delta m_{\rm sol}^2$. The model with the $1"$ singlet flavon gives the
prediction of $\sin^2 \theta_{12}$ around the best fit of NuFIT 3.2 while
keeping near the maximal mixing of $\theta_{23}$. Inputting the experimental
mixing angles with the $1\,\sigma$ errorbar, the Dirac CP violating phase is
clearly predicted to be $\delta_\text{CP}=50^\circ 120^\circ$, which will be
tested by the precise observed value in the future. In order to get the best
fit value $\sin^2\theta_{23}=0.538$, the sum of three neutrino masses is
predicted to be larger than $90\,$meV. The cosmological observation for the sum
of neutrino masses will also provide a crucial test of our predictions. It is
remarked that the model is consistent with the experimental data only for the
normal hierarchy of neutrino masses.

We study a model with a downtype SU(2) singlet vectorlike quark (VLQ) as a
minimal extension of the standard model (SM). In this model, flavor changing
neutral currents (FCNCs) arise at tree level and the unitarity of the $3\times
3$ CabibboKobayashiMaskawa (CKM) matrix does not hold. In this paper, we
constrain the FCNC coupling from $b\rightarrow s$ transitions, especially
$B_s\rightarrow \mu^+\mu^$ and $\bar{B}\rightarrow X_s\gamma$ processes. In
order to analyze these processes, we derive an effective Lagrangian which is
valid below the electroweak symmetry breaking scale. For this purpose, we first
integrate out the VLQ field and derive an effective theory by matching Wilson
coefficients up to oneloop level. Using the effective theory, we construct the
effective Lagrangian for $b\rightarrow s\gamma^{(*)}$. It includes the effects
of the SM quarks and the violation of the CKM unitarity. We show the
constraints on the magnitude of the FCNC coupling and its phase by taking
account of the current experimental data on $\Delta M_{B_s}$,
$\mathrm{Br}[B_s\rightarrow\mu^+\mu^]$, $\mathrm{Br}[\bar{B}\rightarrow
X_s\gamma]$ and CKM matrix elements as well as theoretical uncertainties. We
find that the constraint from the $\mathrm{Br}[B_s\rightarrow\mu^+\mu^]$ is
more stringent than that from the $\mathrm{Br}[\bar{B}\rightarrow X_s\gamma$].
We also obtain the bound for the mass of the VLQ and the strength of the Yukawa
couplings related to the FCNC coupling of $b\rightarrow s$ transition. Using
the CKM elements which satisfy above constraints, we show how the unitarity is
violated on the complex plane.

We discuss a supersymmetric model with discrete flavor symmetry $A_4\times
Z_3$. The additional scalar fields which contribute masses of leptons in the
Yukawa terms are introduced in this model. We analyze their scalar potential
and find that they have various vacuum structures. We show the relations among
24 different vacua and classify them into two types. We derive expressions of
the lepton mixing angles, Dirac CP violating phase and Majorana phases for the
two types. The model parameters which are allowed by the experimental data of
the lepton mixing angles are different for each type. We also study the
constraints on the model parameters which are related to Majorana phases. The
different allowed regions of the model parameters for the two types are shown
numerically for a given region of two combinations of the CP violating phases.

We study the minimal seesaw model, where two righthanded Majorana neutrinos
are introduced, focusing on the CP violating phase. In addition, we take the
trimaximal mixing pattern for the neutrino flavor where the charged lepton mass
matrix is diagonal. Thanks to this symmetric framework, the $3\times 2$ Dirac
neutrino mass matrix is given in terms of a few parameters. Numerical studies
reveal that the observation of the CP violating phase can determine the flavor
structure of the Dirac neutrino mass matrix in the minimal seesaw model. In
particular, new minimal Dirac neutrino mass matrices are proposed in the case
of $\rm TM_1$, which is derived by the additional 23 family mixing to the
tribimaximal mixing basis in the normal hierarchy of neutrino masses. Our
analyses include the Littlest seesaw model by King {\it et al.}, which is one
of the specific one in our results. Furthermore, it is remarked that our
$3\times 2$ Dirac neutrino mass matrix is reproduced by introducing gauge
singlet flavons with the specific alignments of the VEV's. These alignments
suggest the residual symmetry of $S_4$ group.

We discuss the correlation between the CP violating Dirac phase of the lepton
mixing matrix and the cosmological baryon asymmetry based on the leptogenesis
in the minimal seesaw model with two righthanded Majorana neutrinos and the
trimaximal mixing for neutrino flavors. The sign of the CP violating Dirac
phase at low energy is fixed by the observed cosmological baryon asymmetry
since there is only one phase parameter in the model. According to the recent
T2K and NO$\nu$A data of the CP violation, the Dirac neutrino mass matrix of
our model is fixed only for the normal hierarchy of neutrino masses.

Flavor symmetries are useful to realize fermion flavor structures in the
standard model. In particular, discrete $A_4$ symmetry is used to realize
lepton flavor structures, and some scalars which are called flavon are
introduced to break this symmetry. In many models, flavons are assumed to be
much heavier than the electroweak scale. However, our previous work showed that
flavon mass around 100 GeV is allowed by experimental constraints in the $A_4$
symmetric model with residual $Z_3$ symmetry. In this paper, we discuss
collider search of such a light flavon $\varphi_T$. We find that an electron 
photon collision, as a considerable option at the international linear
collider, has advantages to search for the signals. At the electron  photon
collider flavons are produced as $e^\gamma \to l^ \varphi_T$ and decay into
two charged leptons. Then we analyze signals of flavorconserving finalstate
$\tau^+ \tau^ e^$, and flavorviolating finalstates $\tau^+ \mu^ \mu^$ and
$\mu^+ \tau^ \tau^$ by carrying out numerical simulation. For the former
finalstate, SM background can be strongly suppressed by imposing cuts on the
invariant masses of finalstate leptons. For the later finalstates, SM
background is extremely small, because in the SM there are no such
flavorviolating finalstates. We then find that sufficient discovery
significance can be obtained, even if flavons are heavier than the lower limits
from flavor physics.

In $M$estimation under standard asymptotics, the weak convergence combined
with the polynomial type large deviation estimate of the associated statistical
random field Yoshida (2011) provides us with not only the asymptotic
distribution of the associated $M$estimator but also the convergence of its
moments, the latter playing an important role in theoretical statistics. In
this paper, we study the above program for statistical random fields of
multiple and also possibly mixedrates type in the sense of Radchenko (2008)
where the associated statistical random fields may be nondifferentiable and
may fail to be locally asymptotically quadratic. Consequently, a very strong
mode of convergence of a wide range of regularized $M$estimators is ensured.
The results are applied to regularized estimation of an ergodic diffusion
observed at high frequency.

The model of neutrino mass matrix with minimal texture is now tightly
constrained by experiment so that it can yield a prediction for the phase of CP
violation. This phase is predicted to lie in the range $\delta_{CP}=0.77\pi 
1.24\pi$. If neutrino oscillation experiment would find the CP violation phase
outside this range, this means that the minimaltexture neutrino mass matrix,
the element of which is all real, fails and the neutrino mass matrix must be
complex, i.e., the phase must be present that is responsible for leptogenesis.

We present a model with $A_4 \times U(1)_{PQ}$ lepton flavor symmetry which
explains the origin of the lepton flavor structure and also solves the strong
$CP$ problem. Standard model gauge singlet fields, socalled "flavons", charged
under the $A_4 \times U(1)_{PQ}$ symmetry are introduced and are coupled with
the lepton and the Higgs sectors. The flavon vacuum expectation values (VEVs)
trigger spontaneous breaking of the $A_4 \times U(1)_{PQ}$ symmetry. The
breaking pattern of the $A_4$ accounts for the tribimaximal neutrino mixing
and the deviation from it due to the nonzero $\theta_{13}$ angle, and the
breaking of the $U(1)_{PQ}$ gives rise to a pseudoNambuGoldstone boson,
axion, whose VEV cancels the QCD $\theta$ term. We investigate the breaking of
the $A_4 \times U(1)_{PQ}$ symmetry through an analysis on the scalar potential
and further discuss the properties of the axion in the model, including its
decay constant, mass and coupling with photons. It is shown that the axion
decay constant is related with the righthanded neutrino mass through the
flavon VEVs. Experimental constraints on the axion and their implications are
also studied.

We discuss the neutrino flavor structures in the Occam's razor approach for
the Dirac neutrino mass matrices. We assume that the charged lepton mass matrix
takes a diagonal base, while the righthanded Majorana neutrino mass matrix is
also diagonal and we consider nine patterns of the four zero textures for the
Dirac neutrinos mass matrices. We numerically analyze the lefthanded Majorana
neutrino mass matrices for nine patterns of the Dirac neutrino mass matrices
and we found two interesting patterns where both normal and inverted neutrino
mass hierarchies can be realized. We also found if the neutrino mass is normal
hierarchy, this scenario will be likely to be excluded by the measurements of
the Dirac CP violating phase in the T2K and NO$\nu $A neutrino experiments, for
instance. On the other hand, if the neutrino mass is inverted hierarchy, this
scenario will be also likely to be excluded by the measurements of neutrinoless
double beta decay experiments, e.g., KamLANDZen experiment.

We present a doubly parametric extension of the standard FroggattNielsen
(FN) mechanism. As is well known, mass matrices of the up and downtype quark
sectors and the charged lepton sector in the standard model can be parametrized
well by a parameter $\lambda$ which is usually taken to be the sine of the
Cabibbo angle ($\lambda = \sin\theta_\text{C} \simeq 0.225$). However, in the
neutrino sector, there is still room to realize the two neutrino mass squared
differences $\Delta m_\text{sol}^2$ and $\Delta m_\text{atm}^2$, two large
mixing angles $\theta _{12}$ and $\theta _{23}$, and nonzero $\theta _{13}$.
Then we consider an extension with an additional parameter $\rho$ in addition
to $\lambda$. Taking the relevant FN charges for a power of $\lambda~(=0.225)$
and additional FN charges for a power of $\rho$, which we assume to be less
than one, we can reproduce the ratio of the two neutrino mass squared
differences and three mixing angles. In the normal neutrino mass hierarchy, we
show several patterns for taking relevant FN charges and the magnitude of
$\rho$. We find that if $\sin \theta_{23}$ is measured more precisely, we can
distinguish each pattern. This is testable in the near future, for example in
neutrino oscillation experiments. In addition, we predict the Dirac
CPviolating phase for each pattern.

We present a modified Altarelli and Feruglio $A_4$ model where an additional
$A_4$ singletprime flavon is introduced. In this model, nonzero $\theta
_{13}$ is given by this additional $A_4$ singletprime flavon which breaks
tribimaximal mixing. In the framework of the supersymmetry with $U(1)_R$
symmetry, we obtain vacuum expectation values (VEVs) and VEV alignments of
flavons through driving fields. It is considered that flavon induces
distinctive flavor violating process if flavon mass is light. Assuming mass of
SUSY particles are sufficiently heavy so that the SUSY contributions can be
negligible, we discuss the flavor violating Yukawa interaction through flavon
exchange in the charged lepton sector. According to the potential analysis, the
VEV of flavon breaks $A_4$ down to $Z_3$ in the charged lepton sector and
relation among flavon masses is determined. Thanks for the residual $Z_3$
symmetry, many lepton flavor violating decay modes are forbidden except for
$\tau \rightarrow \mu \mu \bar{e}$ and $\tau \rightarrow e e \bar{\mu}$. A mass
limit of the flavon from these threebody decay modes is $60$~GeV taking into
account the current experimental lower bounds at the Belle experiment. In our
model, we predict a ratio of the branching ratios $\tau \rightarrow \mu \mu
\bar{e}$ and $\tau \rightarrow e e \bar{\mu}$ by using known charged lepton
masses. We also find that the production cross section for the flavon can be
$\mathcal{O}(1)$ fb. Thus the flavon would be found at the LHC run 2 by
searching for 4tau lepton signal.

We discuss the neutrino mass matrix based on the Occam's Razor approach in
the framework of the seesaw mechanism. We impose four zeros in the Dirac
neutrino mass matrix, which give the minimum number of parameters needed for
the observed neutrino masses and lepton mixing angles, while the charged lepton
mass matrix and the righthanded Majorana neutrino mass matrix are taken to be
real diagonal ones. The lowenergy neutrino mass matrix has only seven physical
parameters. We show successful predictions for the mixing angle $\theta_{13}$
and the CP violating phase $\delta_{CP}$ with the normal mass hierarchy of
neutrinos by using the experimental data on the neutrino mass squared
differences, the mixing angles $\theta_{12}$ and $\theta_{23}$. The most
favored region of $\sin\theta_{13}$ is around $0.13\sim 0.15$, which is
completely consistent with the observed value. The CP violating phase
$\delta_{CP}$ is favored to be close to $\pm \pi/2$. We also discuss the
Majorana phases as well as the effective neutrino mass for the neutrinoless
doublebeta decay $m_{ee}$, which is around $7\sim 8$ meV. It is extremely
remarkable that we can perform a "complete experiment" to determine the
lowenergy neutrino mass matrix, since we have only seven physical parameters
in the neutrino mass matrix. In particular, two CP violating phases in the
neutrino mass matrix are directly given by two CP violating phases at high
energy. Thus, assuming the leptogenesis we can determine the sign of the cosmic
baryon in the universe from the lowenergy experiments for the neutrino mass
matrix.

Nonconcave penalized maximum likelihood methods, such as the Bridge, the
SCAD, and the MCP, are widely used because they not only do parameter
estimation and variable selection simultaneously but also have a high
efficiency as compared to the Lasso. They include a tuning parameter which
controls a penalty level, and several information criteria have been developed
for selecting it. While these criteria assure the model selection consistency
and so have a high value, it is a severe problem that there are no appropriate
rules to choose the one from a class of information criteria satisfying such a
preferred asymptotic property. In this paper, we derive an information
criterion based on the original definition of the AIC by considering the
minimization of the prediction error rather than the model selection
consistency. Concretely speaking, we derive a function of the score statistic
which is asymptotically equivalent to the nonconcave penalized maximum
likelihood estimator, and then we provide an asymptotically unbiased estimator
of the KullbackLeibler divergence between the true distribution and the
estimated distribution based on the function. Furthermore, through simulation
studies, we check that the performance of the proposed information criterion
gives almost the same as or better than that of the crossvalidation.

We propose two new simple lepton flavor models in the framework of the $S_4$
flavor symmetry. The neutrino mass matrices, which are given by two complex
parameters, lead to the inverted mass hierarchy. The charged lepton mass matrix
has the 12 lepton flavor mixing, which gives the nonvanishing reactor angle
$\theta_{13}$. These models predict the Dirac phase and the Majorana phases,
which are testable in the future experiments. The predicted magnitudes of the
effective neutrino mass for the neutrinoless double beta decay are in the
regions as $32~\text{meV}\lesssim m_{ee}\lesssim 49~\text{meV}$ and
$34~\text{meV}\lesssim m_{ee}\lesssim 59~\text{meV}$, respectively. These
values are close to the expected reaches of the coming experiments. The total
sum of the neutrino masses are predicted in both models as
$0.0952~\text{eV}\lesssim \sum m_i\lesssim 0.101~\text{eV}$ and
$0.150~\text{eV}\lesssim \sum m_i\lesssim 0.160~\text{eV}$, respectively.

We propose an update estimation method for a diffusion parameter from
highfrequency dependent data under a nuisance drift element. We ensure the
asymptotic equivalence of the estimator to the corresponding quasiMLE, which
has the asymptotic normality and the asymptotic efficiency. We give a
simulation example to illustrate the theory.

We study lepton flavor violating Higgs decays in two models, with the
recently found hint for Higgs $\to \mu \tau$ at CMS as a benchmark value for
the branching ratio. The first model uses the discrete flavor symmetry group
$A_4$, broken at the electroweak scale, while the second is renormalizable and
based on the Abelian gauge group $L_\mu  L_\tau$. Within the models we find
characteristic predictions for other nonstandard Higgs decay modes, charged
lepton flavor violating decays and correlations of the branching ratios with
neutrino oscillation parameters.

We study the CP violation in the deviation from the tribimaximal mixing
(TBM) of neutrinos. We examine nontrivial relations among the mixing angles
and the CP violating Dirac phase in the typical four cases of the deviation
from the TBM. The first two cases are derived by the additional rotation of the
23 or 13 generations of neutrinos in the TBM basis. The other two cases are
given by the additional rotation of the 13 or 12 generations of charged
leptons with the TBM neutrinos. These four cases predict different relations
among three mixing angles and the CP violating Dirac phase. The rotation of the
23 generations of neutrinos in the TBM basis predicts $\sin ^2\theta
_{12}<1/3$, and the CP violating Dirac phase to be $\pm (0.09\pi\sim 0.76\pi)$
for NH ($\pm (0.15\pi\sim 0.73\pi) \ \text{for IH}$) depending on $\sin
^2\theta _{23}$. The rotation of the 13 generations of neutrinos in the TBM
basis gives $\sin ^2\theta _{12}>1/3$. The CP violating Dirac phase is not
constrained by the input of the present experimental data. For the case of the
13 and 12 rotations of charged leptons in the TBM basis, the CP violating
Dirac phase is predicted in $\pm(0.35\pi\sim 0.60\pi)$ depending on $\sin
^2\theta _{12}$ for both NH and IH cases. We also discuss the specific case
that $\theta_{13}$ is related with the Cabibbo angle $\lambda$ such as
$\sin\theta_{13}=\lambda/\sqrt{2}$, in which the maximal CP violation is
preferred. The CP violating Dirac phase can distinguish the lepton flavor
mixing patterns at T2K and NO$\nu$A experiments in the future.

The origin of fermion mass hierarchies and mixings is one of the unresolved
and most difficult problem in highenergy physics. One possibility to address
the flavour problem is by extending the Standard Model to include a family
symmetry. In the recent years it has become very popular to use nonAbelian
discrete flavour symmetries because of their power in the prediction of the
large leptonic mixing angles relevant for neutrino oscillation experiments.
Here we give an introduction to the flavour problem and to discrete groups
which have been used to attempt a solution for it. We review the current status
of models in the light of the recent measurement of the reactor angle and we
consider different model building directions taken. The use of the flavons or
multi Higgs scalars in model building is discussed as well as the direct vs.
indirect approaches. We also focus on the possibility to distinguish
experimentally flavour symmetry models by means of mixing sum rules and mass
sum rules. In fact, we illustrate in this review the complete path from
mathematics, via model building, to experiments, so that any reader interested
to start working in the field could use this text as a starting point in order
to get a broad overview of the different subject areas.

We study CP violations in the B_s> K+K and Bs>K0K0 decays in order to find
the contribution of the supersymmetry, which comes from the gluinosquark
mediated flavor changing current. We obtain the allowed region of the squark
flavor mixing parameters by putting the experimental data, the mass difference
Delta M_Bs, the CP violating phase phi_s in Bs to J/psi phi decay and the b to
s gamma branching ratio. In addition to these data, we take into account the
constraint from the asymmetry of B0>K+pi because the Bs>K+K decay is related
with the B0>K+pi decay by replacing the spectator s with d. Under these
constraints, we predict the magnitudes of the CP violation in the Bs>K+K and
Bs>K0K0 decays. The predicted region of the CP violation C_{K+K} is strongly
cut from the direct CP violation of barB0 to Kpi+, therefore, the deviation
from the SM prediction of C_{K+K} is not found. On the other hand, the CP
violation S_{K+K} is possibly deviated from the SM prediction considerably, in
the region of 0.1 0.5. Since the standard model predictions of C_{K0bar K0}
and S_{K0bar K0} are very small, the squark contribution can be detectable in
C_{K0bar K0} and S_{K0bar K0}. These magnitudes are expected in the region
C_{K0bar K0}=0.060.06 and S_{K0bar K0}=0.50.3. More precise data of these
CP violations provide us a crucial test for the gluinosquark mediated flavor
changing current.

We study the sensitivity of the squark flavor mixing to the CP violating
phenomena of $K$, $B^0$ and $B_s$ mesons in the framework of the splitfamily
scenario, where the first and second family squarks are very heavy, ${\cal
O}(10)$TeV, on the other hand, the third family squark masses are at ${\cal
O}(1)$TeV. In order to constrain the gluinosbottomquark mixing parameters, we
input the experimental data of the CP violations of $K$, $B^0$, and $B_s$
mesons, that is $\epsilon_K$, $\phi_d$, and $\phi_s$. The experimental upper
bound of the chromoEDM of the strange quark is also input. In addition, we
take account of the observed values $\Delta M_{B^0}$, $\Delta M_{B_s}$, the CKM
mixing $V_{ub}$, and the branching ratio of $b\to s\gamma$. The allowed
region of the mixing parameters are obtained as $\delta_{13}^{dL(dR)}=0\sim
0.01$ and $\delta_{23}^{dL(dR)}=0\sim 0.04$. By using these values, the
deviations from the SM are estimated in the CP violations of the $B^0$ and
$B_s$ decays. The deviation from the SM one is tiny in the CP asymmetries of
$B^0\to \phi K_S$ and $B^0\to \eta 'K^0$ due to the chromoEDM of the strange
quark. On the other hand, the CP asymmetries $B_s \to \phi \phi$ and $B_s \to
\phi \eta '$ could be largely deviated from the SM predictions. We also predict
the time dependent CP asymmetry of $B^0\to K^0\bar K^0$ and the semileptonic
CP asymmetries of $B^0 \to \mu ^X$ and $B_s \to \mu ^X$. We expect those
precise measurements at Belle II, which will provide us interesting tests for
the squark flavor mixing.

We study the contribution of the gluinosquark mediated flavor changing
process for the CP violation in $b\to s$ and $b\to d$ transitions facing on
recent experimental data. The mass insertion parameters of squarks are
constrained by the branching ratios of $b\to s\gamma $ and $b\to d\gamma $
decays. In addition, the time dependent CP asymmetries of $B^0\to \phi K_S$ and
$B^0\to \eta ' K^0$ decays severely restrict the allowed region of the mass
insertion parameter for the $b\to s$ transition. By using these constraints
with squark and gluino masses of 1.5 TeV, we predict the CP asymmetries of
$B_s\to \phi \phi $, $B_s\to \eta '\phi $, and $B^0\to K^0\bar K^0$ decays, as
well as the CP asymmetries in $b\to s\gamma $ and $b\to d\gamma $ decays. The
CP violation in the $B_s\to \phi \phi$ decay is expected to be large owing to
the squark flavor mixing, which will be tested at LHCb soon.

We present neutrino mass matrix textures in a minimal framework of the typeI
seesaw mechanism where two righthanded Majorana neutrinos are introduced in
order to reproduce experimental results of neutrino oscillations. The textures
can lead to experimentally favored leptonic mixing angles described by the
tribimaximal mixing with one additional rotation. We present minimal and next
to minimal textures for the normal mass hierarchy case in a context of the
texture zero. A minimal texture in the inverted hierarchy case is also
constructed, which does not have any vanishing entries in a Dirac neutrino mass
matrix. We also discuss some cases that model parameters in the textures are
supposed to be a neutrino mass ratio and/or the Cabibbo angle. Predicted
regions of mixing angles, a leptonic CPviolation parameter, and an effective
mass for the neutrinoless double beta decay are presented in all textures.

We show that the neutrino mass matrix with the minimal texture leads to the
neutrino mass and mixing that are empirically determined, including
$\theta_{13}$ that is measured recently at a high precision. This mass matrix,
where the Majorana nature of the neutrino is essential, only allows the normal
hierarchy of the neutrino mass, excluding either inverse mass hierarchy or
degenerate mass neutrinos, and hence predicts the effective mass of double beta
decay to lie within the range $m_{ee}=3.75.6$ eV.