Predictions for the neutrino parameters in the minimal model extended by linear combination of U(1)$_{L_e-L_\mu}$, U(1)$_{L_\mu-L_\tau}$ and U(1)$_{B-L}$ gauge symmetries

TL;DR

This paper explores minimal extensions of the Standard Model with specific U(1) gauge symmetries, deriving precise predictions for neutrino parameters, mass hierarchies, and implications for leptogenesis based on constrained mass matrix structures.

Contribution

It introduces a novel framework combining multiple U(1) symmetries to restrict neutrino mass matrices, leading to unique predictions for neutrino parameters and compatibility with recent data.

Findings

Predicted all neutrino parameters as functions of mixing angles and mass differences.

Identified seven minimal models consistent with current neutrino oscillation data.

Provided predictions for the sum of neutrino masses and effective Majorana mass.

Abstract

We study the minimal extensions of the Standard Model by a linear combination of U(1)$_{L_e-L_\mu}$, U(1)$_{L_\mu-L_\tau}$ and U(1)$_{B-L}$ gauge symmetries, where three right-handed neutrinos and one U(1)-breaking SU(2)$_L$ singlet or doublet scalar are introduced. Because of the dependence on the lepton flavor, the structures of both Dirac and Majorana mass matrices of neutrinos are restricted. In particular, the two-zero minor and texture structures in the mass matrix for the active neutrinos are interesting. Analyzing these structures, we obtain uniquely all the neutrino parameters, namely the Dirac CP phase $\delta$, the Majorana CP phases $\alpha_{2,3}$ and the mass eigenvalues of the light neutrinos $m_i$ as functions of the neutrino mixing angles $\theta_{12}$, $\theta_{23}$, and $\theta_{13}$, and the squared mass differences $\Delta m^2_{21}$ and $\Delta m^2_{31}$. In 7 minimal models which are consistent with the recent neutrino oscillation data, we also obtain the predictions for the sum of the neutrino masses $\Sigma_i m_i$ and the effective Majorana neutrino mass $\langle m_{\beta \beta} \rangle$ and compare them with the current experimental limits. In addition, we also discuss the implication of our results for leptogenesis.

Tables

Table 1. Table 1 : U(1) Y ′ superscript 𝑌 ′ {}_{Y^{\prime}} charges of SM fermions and right-handed neutrinos

field	quarks	leptons			Higgs
		$e,{\nu }_e,N_e$	$\mu ,{\nu }_{\mu },N_{\mu }$	$\tau ,{\nu }_{\tau },N_{\tau }$
charge	$x_q$	$x_e$	$x_{\mu }$	$x_{\tau }$	0

Table 2. Table 2 : Gauged U(1) flavor symmetries which realize the TZM or TZT structures and their structural patterns. The notations for the structural patterns of the TZM and TZT structures  [ 75 , 76 ] are written in Tab.  3 . The index 𝐑 𝐑 \mathbf{R} ( 𝝂 𝝂 \bm{\nu} ) means that the corresponding symmetry realizes TZM (TZT) structure of the light neutrino mass matrix. Numbers in parentheses represent the U(1) Y ′ superscript 𝑌 ′ {}_{Y^{\prime}} charges of the U(1) Y ′ superscript 𝑌 ′ {}_{Y^{\prime}} -breaking scalar singlets and doublets and “ − - ” means that the corresponding U(1) Y ′ superscript 𝑌 ′ {}_{Y^{\prime}} cannot realize the TZM or TZT structure.

	$Y^{\prime }$	Structural pattern
		Singlet	Doublet
$e\mu \tau$	$L_e-L_{\mu }$	${𝐄}_{\mathrm{𝟏}}^{𝐑}\left(+1\right)$	${𝐀}_{\mathrm{𝟐}}^{𝝂}\left(+1\right)$	${𝐃}_{\mathrm{𝟏}}^{𝝂}\left(-1\right)$
	$L_{\mu }-L_{\tau }$	${𝐂}^{𝐑}\left(+1\right)$	${𝐁}_{\mathrm{𝟑}}^{𝝂}\left(+1\right)$	${𝐁}_{\mathrm{𝟒}}^{𝝂}\left(-1\right)$
	$L_e-L_{\tau }$	${𝐄}_{\mathrm{𝟐}}^{𝐑}\left(+1\right)$	${𝐀}_{\mathrm{𝟏}}^{𝝂}\left(+1\right)$	${𝐃}_{\mathrm{𝟐}}^{𝝂}\left(-1\right)$
$Be\mu \tau$	$B-3L_e-L_{\mu }+L_{\tau }$	${𝐀}_{\mathrm{𝟏}}^{𝐑}\left(+2\right)$	$-$	$-$
	$B-3L_e+L_{\mu }-L_{\tau }$	${𝐀}_{\mathrm{𝟐}}^{𝐑}\left(+2\right)$	$-$	$-$
	$B-L_e-3L_{\mu }+L_{\tau }$	${𝐁}_{\mathrm{𝟑}}^{𝐑}\left(+2\right)$	$-$	$-$
	$B-L_e+L_{\mu }-3L_{\tau }$	${𝐁}_{\mathrm{𝟒}}^{𝐑}\left(+2\right)$	$-$	$-$
	$B+L_e-3L_{\mu }-L_{\tau }$	${𝐃}_{\mathrm{𝟏}}^{𝐑}\left(+2\right)$	$-$	$-$
	$B+L_e-L_{\mu }-3L_{\tau }$	${𝐃}_{\mathrm{𝟐}}^{𝐑}\left(+2\right)$	$-$	$-$
	$B-3L_e$	${𝐅}_{\mathrm{𝟏}}^{𝐑}\left(+6\right)$	$-$	$-$
	$B-\frac{3}{2}{L}_{\mu }-\frac{3}{2}{L}_{\tau }$	${𝐅}_{\mathrm{𝟏}}^{𝐑}\left(+3\right)$	$-$	$-$
	$B-3L_{\mu }$	${𝐅}_{\mathrm{𝟐}}^{𝐑}\left(+6\right)$	$-$	$-$
	$B-\frac{3}{2}{L}_e-\frac{3}{2}{L}_{\tau }$	${𝐅}_{\mathrm{𝟐}}^{𝐑}\left(+3\right)$	$-$	$-$
	$B-3L_{\tau }$	${𝐅}_{\mathrm{𝟑}}^{𝐑}\left(+6\right)$	$-$	$-$
	$B-\frac{3}{2}{L}_e-\frac{3}{2}{L}_{\mu }$	${𝐅}_{\mathrm{𝟑}}^{𝐑}\left(+3\right)$	$-$	$-$

Table 3. Table 3 : The list of the structural patterns of the two-zero minor and texture structure. We follow the notation in  [ 75 , 76 ] .

structure index	${𝐀}_{\mathrm{𝟏}}$	${𝐀}_{\mathrm{𝟐}}$	${𝐁}_{\mathrm{𝟑}}$	${𝐁}_{\mathrm{𝟒}}$	$𝐂$	${𝐃}_{\mathrm{𝟏}}$
pattern
	${𝐃}_{\mathrm{𝟐}}$	${𝐄}_{\mathrm{𝟏}}$	${𝐄}_{\mathrm{𝟐}}$	${𝐅}_{\mathrm{𝟏}}$	${𝐅}_{\mathrm{𝟐}}$	${𝐅}_{\mathrm{𝟑}}$

Table 4. Table 4 : Values for the neutrino oscillation parameters we use in this paper. We take them from the NuFIT v4.0 result with the Super-Kamiokande atmospheric data   [ 104 , 105 ] .

	Normal Ordering		Inverted Ordering
Parameter	Best fit $\pm 1\sigma$	3$\sigma$ range	Best fit $\pm 1\sigma$	3$\sigma$ range
${\sin }^2{\theta }_{12}$	${0.310}_{-0.012}^{+0.013}$	0.275–0.350	${0.310}_{-0.012}^{+0.013}$	0.275–0.350
${\sin }^2{\theta }_{23}$	${0.582}_{-0.019}^{+0.015}$	0.428–0.624	${0.582}_{-0.018}^{+0.015}$	0.433–0.623
${\sin }^2{\theta }_{13}$	${0.02240}_{-0.00066}^{+0.00065}$	0.02044–0.02437	${0.02263}_{-0.00066}^{+0.00065}$	0.02067–0.02461
$\Delta m_{21}^2/{10}^{-5}{\text{eV}}^2$	${7.39}_{-0.20}^{+0.21}$	6.79–8.01	${7.39}_{-0.20}^{+0.21}$	6.79–8.01
$\Delta m_{3\ell }^2/{10}^{-3}{\text{eV}}^2$	${2.525}_{-0.031}^{+0.033}$	$2.431$–$2.622$	$-{2.512}_{-0.031}^{+0.034}$	$-$(2.606–$2.413$)
$\delta [{}_{}{}^{\circ }]$	${217}_{-28}^{+40}$	135–366	${280}_{-28}^{+25}$	196–351

Table 5. Table 5 : The consistency between the neutrino oscillation parameters and the mass ordering. “ ✓ ✓ \checkmark ” represents the compatibility of the TZM structure with NO or IO when the best fit value of θ 23 subscript 𝜃 23 \theta_{23} is used, and “ × \times ” does the opposite. “ △ △ \triangle ” in the 𝐁 𝟑 𝐑 superscript subscript 𝐁 3 𝐑 \mathbf{B_{3}^{R}} and 𝐁 𝟒 𝐑 superscript subscript 𝐁 4 𝐑 \mathbf{B_{4}^{R}} cases mean that the NO or IO cannot be realized when we use the best fit value of θ 23 subscript 𝜃 23 \theta_{23} , but they are realized when we use the smallest value in the 3 ​ σ 3 𝜎 3\sigma region of θ 23 subscript 𝜃 23 \theta_{23} . The subscript 𝐑 𝐑 \mathbf{R} means that the neutrino mass matrix has the two-zero minor structure.

structure index	${𝐀}_{\mathrm{𝟏}}^{𝐑}$	${𝐀}_{\mathrm{𝟐}}^{𝐑}$	${𝐁}_{\mathrm{𝟑}}^{𝐑}$	${𝐁}_{\mathrm{𝟒}}^{𝐑}$	${𝐂}^{𝐑}$	${𝐃}_{\mathrm{𝟏}}^{𝐑}$	${𝐃}_{\mathrm{𝟐}}^{𝐑}$	${𝐄}_{\mathrm{𝟏}}^{𝐑}$	${𝐄}_{\mathrm{𝟐}}^{𝐑}$	${𝐅}_{\mathrm{𝟏}}^{𝐑}$	${𝐅}_{\mathrm{𝟐}}^{𝐑}$	${𝐅}_{\mathrm{𝟑}}^{𝐑}$
NO	$\times$	$\times$	$✓$	$△$	$✓$	$✓$	$✓$	$\times$	$\times$	$\times$	$\times$	$\times$
IO	$\times$	$\times$	$△$	$✓$	$\times$	$\times$	$\times$	$\times$	$\times$	$\times$	$\times$	$\times$

Table 6. Table 6 : The consistency between the models with the TZM structure and the experimental constraints, namely the neutrino mixing angles, sum of the neutrino masses and effective Majorana neutrino mass. We also show the possibility to generate the correct sign of the baryon asymmetry. “N” and “I” represent the NO and IO, respectively. “ ✓ ✓ \checkmark ” and “ × \times ” represent whether models are consistent with the experimental constraints or not. “ △ △ \triangle ” means that, in the column of mixing, the models predict the NO or the IO if θ 23 subscript 𝜃 23 \theta_{23} is taken to be in the range of 3 ​ σ 3 𝜎 3\sigma , and in that of mass, there is a strong tension between the prediction and the Planck limit. Moreover, in that of leptogenesis, it means that the limited parameter regions realize the correct sign and enough size of the baryon asymmetry. “ □ □ \square ” means that the models predict that the 0 ​ ν ​ β ​ β 0 𝜈 𝛽 𝛽 0\nu\beta\beta process never occurs.

	Structural index
	${𝐀}_{\mathrm{𝟏}}^{𝐑}$	${𝐀}_{\mathrm{𝟐}}^{𝐑}$	${𝐁}_{\mathrm{𝟑}}^{𝐑}$		${𝐁}_{\mathrm{𝟒}}^{𝐑}$		${𝐂}^{𝐑}$		${𝐃}_{\mathrm{𝟏}}^{𝐑}$		${𝐃}_{\mathrm{𝟐}}^{𝐑}$		${𝐄}_{\mathrm{𝟏}}^{𝐑}$	${𝐄}_{\mathrm{𝟐}}^{𝐑}$	${𝐅}_{\mathrm{𝟏}}^{𝐑}$	${𝐅}_{\mathrm{𝟐}}^{𝐑}$	${𝐅}_{\mathrm{𝟑}}^{𝐑}$
	N$\cdot$I	N$\cdot$I	N	I	N	I	N	I	N	I	N	I	N$\cdot$I	N$\cdot$I	N$\cdot$I	N$\cdot$I	N$\cdot$I
Mixing	$\times$	$\times$	$✓$	$△$	$△$	$✓$	$✓$	$\times$	$✓$	$\times$	$✓$	$\times$	$\times$	$\times$	$\times$	$\times$	$\times$
Mass	$-$	$-$	$\times$	$\times$	$\times$	$\times$	$△$	$-$	$✓$	$-$	$✓$	$-$	$-$	$-$	$-$	$-$	$-$
$0\nu \beta \beta$	$-$	$-$	$-$	$-$	$-$	$-$	$✓$	$-$	$\square$	$-$	$\square$	$-$	$-$	$-$	$-$	$-$	$-$
Leptogenesis	$-$	$-$	$-$	$-$	$-$	$-$	$✓$	$-$	$✓$	$-$	$△$	$-$	$-$	$-$	$-$	$-$	$-$

Table 7. Table 7 : The list of the functions R 2 , 3 subscript 𝑅 2 3 R_{2,3} .

structure index	$R_i$
${𝐂}^{𝐑}$	$R_2$	$-\frac{2{\sin }^2{\theta }_{12}\cos 2{\theta }_{23}+\sin 2{\theta }_{12}\sin 2{\theta }_{23}\sin {\theta }_{13}{e}^{i\delta }}{2{\cos }^2{\theta }_{12}\cos 2{\theta }_{23}-\sin 2{\theta }_{12}\sin 2{\theta }_{23}\sin {\theta }_{13}{e}^{i\delta }}$
	$R_3$	$-\frac{\sin {\theta }_{13}{e}^{2i\delta }\left[2\cos 2{\theta }_{12}\cos 2{\theta }_{23}\sin {\theta }_{13}-\sin 2{\theta }_{12}\sin 2{\theta }_{23}\left(e^{-i\delta }+{\sin }^2{\theta }_{13}{e}^{i\delta }\right)\right]}{{\cos }^2{\theta }_{13}\left[2{\cos }^2{\theta }_{12}\cos 2{\theta }_{23}-\sin 2{\theta }_{12}\sin 2{\theta }_{23}\sin {\theta }_{13}{e}^{i\delta }\right]}$
${𝐃}_{\mathrm{𝟏}}^{𝐑}$	$R_2$	$\frac{\sin {\theta }_{12}\left(e^{i\delta }\cos {\theta }_{12}\sin {\theta }_{23}\sin {\theta }_{13}+\sin {\theta }_{12}\cos {\theta }_{23}\right)}{\cos {\theta }_{12}\left(e^{i\delta }\sin {\theta }_{12}\sin {\theta }_{23}\sin {\theta }_{13}-\cos {\theta }_{12}\cos {\theta }_{23}\right)}$
	$R_3$	$-\frac{e^{i\delta }\sin {\theta }_{13}\left(e^{i\delta }\cos {\theta }_{12}\sin {\theta }_{23}\sin {\theta }_{13}+\sin {\theta }_{12}\cos {\theta }_{23}\right)}{\cos {\theta }_{12}\sin {\theta }_{23}{\cos }^2{\theta }_{13}}$
${𝐃}_{\mathrm{𝟐}}^{𝐑}$	$R_2$	$\frac{\sin {\theta }_{12}\left(e^{i\delta }\cos {\theta }_{12}\cos {\theta }_{23}\sin {\theta }_{13}-\sin {\theta }_{12}\sin {\theta }_{23}\right)}{\cos {\theta }_{12}\left(e^{i\delta }\sin {\theta }_{12}\cos {\theta }_{23}\sin {\theta }_{13}+\cos {\theta }_{12}\sin {\theta }_{23}\right)}$
	$R_3$	$\frac{e^{i\delta }\sin {\theta }_{13}\left(-e^{i\delta }\cos {\theta }_{12}\cos {\theta }_{23}\sin {\theta }_{13}+\sin {\theta }_{12}\sin {\theta }_{23}\right)}{\cos {\theta }_{12}\cos {\theta }_{23}{\cos }^2{\theta }_{13}}$

Table 8. Table 8 : The list of the formulae of the neutrino mass m 1 subscript 𝑚 1 m_{1} .

structure index	$m_1$
${𝐂}^{𝐑}$	$\delta m{\left[\frac{4s_{12}^2{\cos }^{2}2{\theta }_{23}+4s_{12}^3{c}_{12}{s}_{13}\sin 4{\theta }_{23}\cos \delta +s_{13}^2{\sin }^{2}2{\theta }_{12}{\sin }^{2}2{\theta }_{23}}{2\left(2\cos 2{\theta }_{12}{\cos }^{2}2{\theta }_{23}-s_{13}\sin 2{\theta }_{12}\sin 4{\theta }_{23}\cos \delta \right)}\right]}^{\frac{1}{2}}$
${𝐃}_{\mathrm{𝟏}}^{𝐑}$	$\delta m{\left[\frac{s_{12}^2\left(2c_{12}^2{s}_{23}^2{s}_{13}^2+2s_{12}^2{c}_{23}^2+\sin 2{\theta }_{12}\sin 2{\theta }_{23}{s}_{13}\cos \delta \right)}{2c_{23}\left(\cos 2{\theta }_{12}{c}_{23}-\sin 2{\theta }_{12}{s}_{23}{s}_{13}\cos \delta \right)}\right]}^{\frac{1}{2}}$
${𝐃}_{\mathrm{𝟐}}^{𝐑}$	$\delta m{\left[\frac{s_{12}^2\left(2c_{12}^2{c}_{23}^2{s}_{13}^2+2s_{12}^2{s}_{23}^2-\sin 2{\theta }_{12}\sin 2{\theta }_{23}{s}_{13}\cos \delta \right)}{2s_{23}\left(\cos 2{\theta }_{12}{s}_{23}+\sin 2{\theta }_{12}{c}_{23}{s}_{13}\cos \delta \right)}\right]}^{\frac{1}{2}}$

Equations

$$9x_q+x_e+x_{\mu }+x_{\tau }=0.$$

$$\displaystyle Y^{\prime}=\left\{\begin{array}[]{ll}x_{e}L_{e}+x_{\mu}L_{\mu}-(x_{e}+x_{\mu})L_{\tau}&\ \ \ (x_{q}=0)\\ B+x_{e}L_{e}+x_{\mu}L_{\mu}-(3+x_{e}+x_{\mu})L_{\tau}&\ \ \ (x_{q}=\frac{1}{3})\end{array}\right.~{}.$$

$$\displaystyle Q_{Y^{\prime}}(\text{Dirac \& charged lepton}):~{}\left\{\begin{array}[]{ll}\begin{pmatrix}0&-x_{e}+x_{\mu}&-2x_{e}-x_{\mu}\\ x_{e}-x_{\mu}&0&-x_{e}-2x_{\mu}\\ 2x_{e}+x_{\mu}&x_{e}+2x_{\mu}&0\end{pmatrix}&\text{for}~{}Y^{\prime}=e\mu\tau\\ \begin{pmatrix}0&-x_{e}+x_{\mu}&-3-2x_{e}-x_{\mu}\\ x_{e}-x_{\mu}&0&-3-x_{e}-2x_{\mu}\\ 3+2x_{e}+x_{\mu}&3+x_{e}+2x_{\mu}&0\end{pmatrix}&\text{for}~{}Y^{\prime}=Be\mu\tau\end{array},\right.$$

$$\displaystyle Q_{Y^{\prime}}(\text{Majorana}):~{}\left\{\begin{array}[]{ll}\begin{pmatrix}-2x_{e}&-x_{e}-x_{\mu}&x_{\mu}\\ -x_{e}-x_{\mu}&-2x_{\mu}&x_{e}\\ x_{\mu}&x_{e}&2(x_{e}+x_{\mu})\end{pmatrix}&\text{for}~{}Y^{\prime}=e\mu\tau\\ \begin{pmatrix}-2x_{e}&-x_{e}-x_{\mu}&3+x_{\mu}\\ -x_{e}-x_{\mu}&-2x_{\mu}&3+x_{e}\\ 3+x_{\mu}&3+x_{e}&2(3+x_{e}+x_{\mu})\end{pmatrix}&\text{for}~{}Y^{\prime}=Be\mu\tau\end{array},\right.$$

$$\Delta \mathcal{L}=$$

$$-{\lambda }_e{N}_e^c(L_e\cdot H)-{\lambda }_{\mu }{N}_{\mu }^c(L_{\mu }\cdot H)-{\lambda }_{\tau }{N}_{\tau }^c(L_{\tau }\cdot H)$$

$$-M_{e\tau }{N}_e^c{N}_{\tau }^c-\frac{1}{2}{\lambda }_{ee}\sigma N_e^c{N}_e^c-{\lambda }_{e\mu }{\sigma }^{\ast }{N}_e^c{N}_{\mu }^c-\frac{1}{2}{\lambda }_{\tau \tau }{\sigma }^{\ast }{N}_{\tau }^c{N}_{\tau }^c+\text{h.c.},$$

$${\mathcal{M}}_D=\frac{v}{\sqrt{2}}\left(\begin{array}{ccc}{\lambda }_e& 0& 0\\ 0& {\lambda }_{\mu }& 0\\ 0& 0& {\lambda }_{\tau }\end{array}\right),{\mathcal{M}}_R=\left(\begin{array}{ccc}{\lambda }_{ee}⟨\sigma ⟩& {\lambda }_{e\mu }⟨\sigma ⟩& M_{e\tau }\\ {\lambda }_{e\mu }⟨\sigma ⟩& 0& 0\\ M_{e\tau }& 0& {\lambda }_{\tau \tau }⟨\sigma ⟩\end{array}\right),{\mathcal{M}}_{\ell }=\frac{v}{\sqrt{2}}\left(\begin{array}{ccc}{y}_e& 0& 0\\ 0& y_{\mu }& 0\\ 0& 0& y_{\tau }\end{array}\right).$$

$${\mathcal{M}}_{{\nu }_L}\simeq -{\mathcal{M}}_D{\mathcal{M}}_R^{-1}{\mathcal{M}}_D^T.$$

$$U_{\text{PMNS}}^T{\mathcal{M}}_{{\nu }_L}{U}_{\text{PMNS}}=\text{diag}(m_1,m_2,m_3),$$

$$U_{\text{PMNS}}=\left(\begin{array}{ccc}{c}_{12}{c}_{13}& s_{12}{c}_{13}& s_{13}{e}^{-i\delta }\\ -s_{12}{c}_{23}-c_{12}{s}_{23}{s}_{13}{e}^{i\delta }& c_{12}{c}_{23}-s_{12}{s}_{23}{s}_{13}{e}^{i\delta }& s_{23}{c}_{13}\\ s_{12}{s}_{23}-c_{12}{c}_{23}{s}_{13}{e}^{i\delta }& -c_{12}{s}_{23}-s_{12}{c}_{23}{s}_{13}{e}^{i\delta }& c_{23}{c}_{13}\end{array}\right)\left(\begin{array}{ccc}1& & \\ & e^{i\frac{{\alpha }_2}{2}}& \\ & & e^{i\frac{{\alpha }_3}{2}}\end{array}\right),$$

$${\mathcal{M}}_{{\nu }_L}^{-1}=-({\mathcal{M}}_D^{-1}{)}^T{\mathcal{M}}_R{\mathcal{M}}_D^{-1}=U_{\text{PMNS}}\text{diag}(m_1^{-1},m_2^{-1},m_3^{-1})U_{\text{PMNS}}^T.$$

$$\frac{1}{m_1}{V}_{\mu 1}^2+\frac{1}{m_2}{V}_{\mu 2}^2{e}^{i{\alpha }_2}+\frac{1}{m_3}{V}_{\mu 3}^2{e}^{i{\alpha }_3}$$

$$\frac{1}{m_1}{V}_{\mu 1}{V}_{\tau 1}+\frac{1}{m_2}{V}_{\mu 2}{V}_{\tau 2}{e}^{i{\alpha }_2}+\frac{1}{m_3}{V}_{\mu 3}{V}_{\tau 3}{e}^{i{\alpha }_3}$$

$${\mathcal{M}}_{{\nu }_L}=-{\mathcal{M}}_D{\mathcal{M}}_R^{-1}{\mathcal{M}}_D^T=U_{\text{PMNS}}^{\ast }\text{diag}(m_1,m_2,m_3)U_{\text{PMNS}}^{†}.$$

$$e^{i{\alpha }_2}$$

$$R_2\equiv \frac{V_{\mu 1}{V}_{e2}^{\ast }}{V_{\mu 2}{V}_{e1}^{\ast }},R_3\equiv \frac{V_{\mu 1}{V}_{e3}^{\ast }}{V_{\mu 3}{V}_{e1}^{\ast }},$$

$$\frac{m_2}{m_1}=\frac{1}{|R_2(\delta )|},\frac{m_3}{m_1}=\frac{1}{|R_3(\delta )|}.$$

$$ϵ|R_2(\delta ){|}^2(1-|R_3(\delta ){|}^2)=|R_3(\delta ){|}^2(1-|R_2(\delta ){|}^2)$$

$${\alpha }_2=\text{Arg}\left[\frac{m_2}{m_1}{R}_2(\delta )\right],{\alpha }_3=\text{Arg}\left[\frac{m_3}{m_1}{R}_3(\delta )\right].$$

$$\displaystyle\begin{array}[]{ccc}\mathcal{M}_{\nu}&&\mathcal{M}_{\nu}^{-1}\\ \hline\cr\mathbf{B_{3}^{\nu}}:\left(\begin{array}[]{ccc}*&0&*\\ 0&0&*\\ &*&*\end{array}\right)&\Longleftrightarrow&\mathbf{B_{4}^{R}}:\left(\begin{array}[]{ccc}*&*&0\\ &*&*\\ 0&*&0\end{array}\right)\\ \mathbf{B_{4}^{\nu}}:\left(\begin{array}[]{ccc}*&*&0\\ &*&*\\ 0&*&0\end{array}\right)&\Longleftrightarrow&\mathbf{B_{3}^{R}}:\left(\begin{array}[]{ccc}*&0&*\\ 0&0&*\\ &*&*\end{array}\right)\\ \hline\cr\end{array}$$

$$⟨m_{\beta \beta }⟩\equiv \left|\sum _i{\left(U_{\text{PMNS}}\right)}_{ei}^2{m}_i\right|=\left|c_{12}^2{c}_{13}^2{m}_1+s_{12}^2{c}_{13}^2{e}^{i{\alpha }_2}{m}_2+s_{13}^2{e}^{i({\alpha }_3-2\delta )}{m}_3\right|.$$

$$⟨m_{\beta \beta }⟩$$

$$=\left|{\left(U_{\text{PMNS}}^{\ast }\text{diag}(m_1,m_2,m_3)U_{\text{PMNS}}^{†}\right)}_{ee}\right|$$

$$=\left|{\left({\mathcal{M}}_{{\nu }_L}\right)}_{ee}\right|$$

$$=0.$$

$$\Delta \mathcal{L}=-\sum _{i=1}^3\sum _{\alpha =e,\mu ,\tau }{\hat{\lambda }}_{i\alpha }{\hat{N}}_i^c(L_{\alpha }\cdot H)-\frac{1}{2}\sum _{i=1}^3{M}_i{\hat{N}}_i^c{\hat{N}}_i^c+h.c.,$$

$${\mathcal{M}}_R$$

$${\hat{N}}_i^c$$

$${\hat{\lambda }}_{i\alpha }$$

$${ϵ}_1$$

$$f(x)$$

$$({\lambda }_e,{\lambda }_{\mu },{\lambda }_{\tau }{)}^T=\lambda (\cos \theta ,\sin \theta \cos \varphi ,\sin \theta \sin \varphi {)}^T\equiv \lambda {\mathbf{n}}^T.$$

$${ϵ}_1$$

$${\hat{n}}_{i\alpha }$$

Full text

UT-19-17

**Predictions for the neutrino parameters in the minimal model extended by linear combination of U(1)${}_{L_{e}-L_{\mu}}$, U(1)${}_{L_{\mu}-L_{\tau}}$ and U(1)B-L gauge symmetries

**

Kento Asai

Department of Physics, University of Tokyo, Bunkyo-ku, Tokyo 133–0033, Japan

We study the minimal extensions of the Standard Model by a linear combination of U(1)${}_{L_{e}-L_{\mu}}$, U(1)${}_{L_{\mu}-L_{\tau}}$ and U(1)B-L gauge symmetries, where three right-handed neutrinos and one U(1)-breaking SU(2)L singlet or doublet scalar are introduced. Because of the dependence on the lepton flavor, the structures of both Dirac and Majorana mass matrices of neutrinos are restricted. In particular, the two-zero minor and texture structures in the mass matrix for the active neutrinos are interesting. Analyzing these structures, we obtain uniquely all the neutrino parameters, namely the Dirac CP phase $\delta$, the Majorana CP phases $\alpha_{2,3}$ and the mass eigenvalues of the light neutrinos $m_{i}$ as functions of the neutrino mixing angles $\theta_{12}$, $\theta_{23}$, and $\theta_{13}$, and the squared mass differences $\Delta m^{2}_{21}$ and $\Delta m^{2}_{31}$. In 7 minimal models which are consistent with the recent neutrino oscillation data, we also obtain the predictions for the sum of the neutrino masses $\Sigma_{i}m_{i}$ and the effective Majorana neutrino mass $\braket{m_{\beta\beta}}$ and compare them with the current experimental limits. In addition, we also discuss the implication of our results for leptogenesis.

1 Introduction

Two decades have passed since the Super-Kamiokande observed the neutrino oscillation in 1998. From that time on, a lot of experiments, such as T2K [1, 2], KamLAND [3], and so on, are conducted, and then the neutrino mixing angles and the squared mass differences are vigorously measured. The observations of neutrino oscillations claim that neutrinos have non-zero masses and we have to modify the neutrino sector in the Standard Model (SM). During these two decades, a lot of modifications were proposed to solve this puzzle. The simplest way to explain the neutrino masses is to introduce right-handed neutrinos into the SM. If they couple to the SM leptons and the Higgs field, the active neutrinos acquire the Dirac masses as the other fermions do. Moreover, the smallness of the neutrino masses can be naturally understood by the seesaw mechanism [4, 5, 6, 7]. For instance, if neutrinos are Majorana particles, and the Majorana masses are much heavier than the Dirac masses, then the tiny neutrino masses are realized by the type-I seesaw mechanism.

On the other hand, it is often considered to extend the SM gauge sector by U(1) gauge symmetries [8, 9, 10, 11]. For instance, it is well known that U(1)B-L can be introduced to the SM without anomalies when three right-handed neutrinos are introduced to the SM. In addition, it is well known that the U(1)${}_{L_{\alpha}-L_{\beta}}$ gauge symmetries, where $L_{\alpha}$ represents the lepton number of generation associated with $\alpha\ (=e,\mu,\tau)$, can be also introduced to the SM without regard to whether three right-handed neutrinos are introduced or not. Above all, the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry is often discussed in the contexts of the muon $g-2$ anomaly [12, 13, 14, 15, 16, 17], flavor physics [18, 19], dark matter [20, 21, 22, 23, 24, 25, 26], and so on. Other recent related studies on the gauged U(1)${}_{L_{\alpha}-L_{\beta}}$ models are found in Refs. [27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55].

In models extended by a lepton flavor-dependent U(1) gauge symmetry (hereafter written as U(1)${}_{Y^{\prime}}$), the structure of the neutrino mass matrix is tightly restricted. Then if we extend only the gauge sector and the extra U(1)${}_{Y^{\prime}}$ gauge symmetry is preserved, the neutrino mass matrix is too sparse to explain the neutrino oscillation data, and therefore we need to break the U(1)${}_{Y^{\prime}}$. To do so, one usually introduces a SU(2)L singlet or doublet scalar with a non-zero U(1)${}_{Y^{\prime}}$ charge such that its VEV spontaneously breaks the U(1)${}_{Y^{\prime}}$ gauge symmetry and gives a mass to the U(1)${}_{Y^{\prime}}$ gauge boson. We, however, note that even though the U(1)${}_{Y^{\prime}}$ gauge symmetry is spontaneously broken, the structure of the neutrino mass matrix is still highly constrained if we introduce only one U(1)${}_{Y^{\prime}}$-breaking scalar field and consider only renormalizable interactions. Particularly, in some U(1)${}_{Y^{\prime}}$ cases, including U(1)${}_{L_{\mu}-L_{\tau}}$, (the inverse of) the neutrino mass matrix has zeros in its two components, and such a structure is called two-zero texture (minor). For previous studies on the neutrino sector of the gauged U(1)${}_{L_{\alpha}-L_{\beta}}$ models, see Refs [56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66]. In Refs. [58, 59], the relations between lepton flavor-dependent U(1) gauge symmetries and structures of the neutrino mass matrix are comprehensively discussed and, previously, in Refs. [67, 68, 69], the relations between two-zero minor and texture structures and constrains on the neutrino parameters are discussed. In Refs. [64, 65], we have discussed them in detail, and have given the predictions for unknown parameters, such as the CP phases, the neutrino masses, and the effective neutrino mass, in the case of the minimal extended model by a U(1)${}_{L_{\alpha}-L_{\beta}}$ gauge symmetry, that we call “the minimal gauged U(1)${}_{L_{\alpha}-L_{\beta}}$ model”, and we have found that these models are already ruled out or are driven into a corner.

In this paper, we have extended the U(1)${}_{L_{\alpha}-L_{\beta}}$ gauge symmetries to a linear combination of the U(1)${}_{L_{e}-L_{\mu}}$, U(1)${}_{L_{\mu}-L_{\tau}}$, and U(1)B-L, and have found viable minimal models. Such a linear combination is the most general lepton flavor-dependent U(1), which can be introduced to the SM gauge sector without anomalies and is consistent with renormalizable Yukawa terms when three right-handed neutrinos are added to the SM. This kind of U(1)s are also discussed, for instance, in Refs. [70, 71, 72, 73]. We analyze the minimal models with two-zero minor or texture structure of the neutrino mass matrix in the same way as we did in Refs. [64, 65], and then we find that only three minimal models are consistent with the recent experiments. We obtain all the CP phases in the PMNS matrix and the mass eigenstates of the light neutrinos as functions of the neutrino oscillation parameters. We also discuss the implication of our results for leptogenesis and give the predictions of the baryon asymmetry of the Universe.

This paper is organized as follows. In Sec. 2, we review the anomaly cancellation conditions and then introduce the minimal gauged U(1) models and explain only fifteen U(1) gauge symmetries realize two-zero minor (TZM) or texture (TZT) structure. In Sec. 3, we show the analyses of the TZM structure following Ref. [64, 65]. We then show in Sec. 4 the predictions for the Dirac CP phase $\delta$, the sum of the neutrino masses $\Sigma_{i}m_{i}$, and the effective mass for the neutrinoless double beta decay $\braket{m_{\beta\beta}}$. In Sec. 5, we discuss the implications for the leptogenesis. Finally, our conclusions are summarized in Sec. 6. In Appendix, we show the explicit expressions for the values and equations which appear in the analyses of the neutrino mass matrices.

2 Minimal gauged U(1) models and neutrino mass matrices

2.1 Anomaly cancellation conditions

For a start, we discuss the conditions that lepton flavor-dependent U(1)${}_{Y^{\prime}}$ gauge symmetries do not cause anomalies, where $Y^{\prime}$ represents the charge of the U(1)${}_{Y^{\prime}}$. We assume as follows:

• •

All quarks have the same U(1)${}_{Y^{\prime}}$ charges.

• •

Left-handed and right-handed fermions have the same U(1)${}_{Y^{\prime}}$ charges.

• •

Charged leptons and neutrinos with the same flavor have the same U(1)${}_{Y^{\prime}}$ charge.

• •

U(1)${}_{Y^{\prime}}$ charge of the Higgs field is zero.

Tab. 1 summarizes the above assumptions.

Under these assumptions, the cubic, gravitational, U(1)${}_{Y}\times$U(1)${}_{Y^{\prime}}\times$U(1)${}_{Y^{\prime}}$ and SU(3)C $\times$SU(3)${}_{C}\times$U(1)${}_{Y^{\prime}}$ anomalies vanish. This is because left-handed and right-handed fermions have the same U(1)${}_{Y^{\prime}}$ charges and the opposite contributions to these anomalies. Therefore, as shown in Refs. [58, 59, 74], the SU(2)${}_{L}\times$SU(2)${}_{L}\times$U(1)${}_{Y^{\prime}}$ and U(1)${}_{Y}\times$U(1)${}_{Y}\times$U(1)${}_{Y^{\prime}}$ anomalies lead the same nontrivial condition of the anomaly cancellation as follows:

[TABLE]

Thus, the charge $Y^{\prime}$ can be written as follows:

[TABLE]

Therefore, only two types of lepton flavor-dependent U(1) gauge symmetries, i.e. U(1)${}_{x_{e}L_{e}+x_{\mu}L_{\mu}-(x_{e}+x_{\mu})L_{\tau}}$ and U(1)${}_{B+x_{e}L_{e}+x_{\mu}L_{\mu}-(3+x_{e}+x_{\mu})L_{\tau}}$ (hereinafter abbreviated to U(1)eμτ and U(1)Beμτ, respectively), can be introduced to the models without anomalies.

Introducing these gauge symmetries, the mass terms of the neutrino Dirac, Majorana and charged lepton are strictly restricted. In the $(e,\mu,\tau)$ basis, the U(1)${}_{Y^{\prime}}$ charges of the neutrino Dirac Yukawa and the charged lepton terms are

[TABLE]

where the $(\alpha,\beta)$ entry in the above matrices represents the U(1)${}_{Y^{\prime}}$ charge of the fermion bilinear term $N_{\alpha}^{c}L_{\beta}$ or $e_{\alpha}^{c}L_{\beta}$, with $\alpha,\beta$ the flavor indices. On the other hand, those of the Majorana mass term of the right-handed neutrinos are

[TABLE]

where the $(\alpha,\beta)$ entry in the above matrices represents the U(1)${}_{Y^{\prime}}$ charge of the fermion bilinear term $N_{\alpha}^{c}N_{\beta}^{c}$. From Eqs. (7) and (10), we find that the Dirac, charged lepton, and Majorana mass matrices are sparse, at least block diagonal, even in any $Y^{\prime}$ charge assignments. Then required values of the neutrino mixing angles cannot be obtained from the simple model, as long as the extra U(1)${}_{Y^{\prime}}$ is preserved and only renormalizable interactions are considered. Therefore, we introduce one SU(2)L scalar singlet or doublet with non-zero U(1)${}_{Y^{\prime}}$ charge 111We can choose a SU(2)L triplet as the U(1)${}_{Y^{\prime}}$ breaking scalar. However, we find, in that case, the resultant neutrino mass matrix cannot have the TZM or the TZT structure, and then we do not consider the triplet case in this paper.. Its VEV leads the U(1)${}_{Y^{\prime}}$ symmetry breaking, and then the Majorana or Dirac mass terms are generated. These mass matrices have various structures depending on U(1)${}_{Y^{\prime}}$ gauge symmetries and additional scalars, and the TZM and TZT structures are especially interesting. This is because the TZM or TZT structure enables us to predict the low-energy neutrino parameters. When (the inverse of) the neutrino mass matrix has more than two-zeros, the vanishing conditions lead too many constraints and there is no solution. On the other hand, when (the inverse of) the neutrino mass matrix has one or no zero, the model has less predictive power and this case is less interesting. That is why we focus on the TZM and TZT cases in this paper. In Tab. 2, we summarize the U(1)${}_{Y^{\prime}}$s which realize the TZM or TZT structure in the minimal models.

The structural patterns in Tab. 2 are summarized in Tab. 3 and we follow the notation in Refs. [75, 76].

As we see in Tab. 2, there are only 15 patterns, and only the minimal gauged U(1)${}_{L_{\alpha}-L_{\beta}}$ models with a doublet scalar can realize the TZT structures. These doublet cases have been discussed in detail in Ref. [65], and so we explain them briefly and show only the results for completeness in this paper. Moreover, we can use the same method as Ref. [64, 65] to analyze the neutrino mass matrices in the cases of general U(1)${}_{Y^{\prime}}$s because the difference between the U(1)${}_{L_{\mu}-L_{\tau}}$ and the other U(1)${}_{Y^{\prime}}$ cases are only the positions of the zero components.

2.2 Minimal gauged U(1) models

Here, using U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$ case, we show the realization of the TZM structure in the minimal gauged U(1)${}_{Y^{\prime}}$ model with a U(1)${}_{Y^{\prime}}$-breaking singlet scalar. The interaction terms in the leptonic sector are then given by

[TABLE]

where the dots indicate the contraction of the SU(2)L indices. After the Higgs field $H$ and the singlet scalar $\sigma$ acquire VEVs222We can always take the VEV of $\sigma$ to be real by using U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$ transformations. $\braket{H}=v/\sqrt{2}$ and $\braket{\sigma}$, the Dirac, Majorana and charged lepton mass matrices are obtained as follows:

[TABLE]

Throughout this work, we assume that the non-zero components in the Majorana mass matrix $\mathcal{M}_{R}$ are much larger than those in the Dirac matrix $\mathcal{M}_{D}$ so that the mass matrix of the active neutrinos is given by the seesaw formula [4, 5, 6, 7]

[TABLE]

We note that any size of the VEV of the singlet scalar field $\sigma$ can be allowed as long as the Majorana masses are much larger than the Dirac masses.

The minimal gauged U(1)${}_{Y^{\prime}}$ models are constrained by many kinds of the experiments relevant to the U(1)${}_{Y^{\prime}}$ gauge boson $Z^{\prime}$ [44, 47, 54]. For example, $Z^{\prime}$ with relatively large gauge coupling like $\gtrsim\mathcal{O}(10^{-3})$ is constrained by the $e^{+}e^{-}$ colliders, such as BABAR [77], Belle [78, 79], KLOE [80, 81, 82], etc. When $Z^{\prime}$ couples to electrons, these experiments give strict constraints. Otherwise, when $Z^{\prime}$ couples to muons, for example, the BABAR experiment search for $e\bar{e}\to\mu\bar{\mu}Z^{\prime}$, $Z^{\prime}\to\mu\bar{\mu}$ constrains $Z^{\prime}$ in around GeV mass range [83]. Other than that, there is the Borexino experiment [84, 85], which gives a bound on the $\nu\mathchar 45\relax e$ interactions and the tests of neutrino trident production [18, 86, 87, 88]. For relatively light mass ($m_{Z^{\prime}}\lesssim 100$ MeV) and small coupling constant ($g_{Z^{\prime}}\lesssim 10^{-3}$), the minimal models are constrained by the electron beam-dump experiments [89, 90], such as SLAC E137 [91], SLAC E141 [92], Fermilab E774 [93], Orsay [94], etc, and the proton beam-dump experiments, such as CHARM [95, 96], U70/Nu-Cal [97, 98], etc.

The doublet case realizes the TZT structure. As we have discussed in Ref. [65], the difference between the singlet and doublet cases is that a doublet scalar with U(1)${}_{Y^{\prime}}$ charge $+1$ realizes different TZT structure from one with U(1)${}_{Y^{\prime}}$ charge $-1$ 333In the singlet case, the exchange of U(1)${}_{Y^{\prime}}$ charge $+1$ for $-1$ just corresponds to that of $\sigma$ for $\sigma^{*}$ and the realized neutrino mass structure does not change.. Furthermore, the U(1)${}_{Y^{\prime}}$-breaking scale should be lower than the electroweak scale, and only the U(1)${}_{L_{\mu}-L_{\tau}}$ avoids the experimental constraints relevant to $Z^{\prime}$.

3 Vanishing conditions and neutrino parameters

In this section, following the analyses of Refs. [64, 65], we analyze the conditions from two-zero minor (texture) structure in the active neutrino mass matrix and see the two conditional equations which the low-energy neutrino parameters should satisfy. By solving these equations, we show that four unknown parameters can be written as functions of the rest of the parameters. As we mentioned before, we explain how to analyze neutrino mass matrix, using the U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$ case, and this analyses can be applied to the other U(1)${}_{Y^{\prime}}$s.

3.1 Vanishing conditions in the neutrino mass matrix

As we said in the previous section, we assume that the Majorana masses are much larger than the Dirac masses so that the mass matrix of the active neutrinos is given by the seesaw formula (Eq. (13)). On the other hand, the charged mass matrix $\mathcal{M}_{\ell}$ is diagonal in this model, and thus we can obtain the mass eigenvalues of the light neutrinos, diagonalizing this matrix by a unitary matrix $U_{\text{PMNS}}$:

[TABLE]

where $U_{\text{PMNS}}$ is the so-called Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [99, 100, 101, 102]. The PMNS matrix is parametrized as

[TABLE]

where $c_{ij}\equiv\cos\theta_{ij}$ and $s_{ij}\equiv\sin\theta_{ij}$ for $\theta_{ij}=[0,\pi/2]$, $\delta,\alpha_{2},\alpha_{3}=[0,2\pi]$, and we have ordered $m_{1}<m_{2}$ without loss of generality. We follow the convention of the Particle Data Group [103], where $m_{2}^{2}-m_{1}^{2}\ll|m_{3}^{2}-m_{1}^{2}|$ and $m_{1}<m_{2}<m_{3}$ (Normal Ordering, NO) or $m_{3}<m_{1}<m_{2}$ (Inverted Ordering, IO). In Tab. 4, we show the best fit values for the neutrino oscillation parameters and $1\sigma\ (3\sigma)$ ranges based on the global analysis of neutrino data from the NuFIT v4.0 result with the Super-Kamiokande atmospheric data  [104, 105]. Here, we take $\ell=1$ for NO and $\ell=2$ for IO in $\Delta m_{3\ell}^{2}$ [106].

In this paper, we consider only the $m_{i}\neq 0$ cases like Refs. [64, 65]. This is because if $m_{i}=0$ ($i=1$ or $3$), the mass matrix for the light neutrinos $\mathcal{M}_{\nu_{L}}$ is block-diagonal, and we cannot have desired mixing angles. From Eqs. (13) and (14), the following relation is given:

[TABLE]

$\mathcal{M}_{D}$ is diagonal and $(\mu,\mu)$ and $(\mu,\tau)$ components in $\mathcal{M}_{R}$ vanish, and so these components of the right-hand side also should vanish. Then the minimal gauged U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$ model is one of the concrete realizations of the two-zero-minor structure. From these vanishing conditions, the following two equations are given:

[TABLE]

where the unitary matrix $V$ is defined by $U_{\text{PMNS}}=V\cdot\text{diag}(1,e^{i\alpha_{2}/2},e^{i\alpha_{3}/2})$. We notice that neither the VEV of the U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$-symmetry breaking singlet scalar $\langle\sigma\rangle$ nor Majorana masses $M_{e\tau}$ appear in these conditions explicitly. This is the most important point in this analysis. This is because the following discussions and results based on the above conditions are independent of the U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$-symmetry breaking scale and the Majorana mass scale. Needless to say, to explain the lightness of the active neutrinos, we use the seesaw mechanism in Eq. (13), and so these scales lie around much greater than the scale of the light neutrinos. Moreover, it is shown in Ref. [64] that the two-zero minor structure is preserved under the renormalization group flow when the charged lepton mass matrix is diagonal.

In the doublet cases, the TZT structure of the neutrino mass matrix gives the following vanishing conditions:

[TABLE]

In the case of the minimal gauged U(1)${}_{L_{\mu}-L_{\tau}}$ model with a doublet scalar with U(1)${}_{L_{\mu}-L_{\tau}}$ charge $+1$, both sides of Eq. (19) have zeros in the $(e,\mu)$ and $(\mu,\mu)$ components. Then, we obtain the vanishing conditions from Eq. (19) and can analyze these in the same way as the singlet case.

3.2 Analyses of vanishing conditions

Eqs. (17) and (18) are two complex equations. Therefore, by solving these conditional equations, we can obtain four real parameters, namely the Dirac CP phase $\delta$, Majorana CP phases $\alpha_{2,3}$, and mass eigenvalue of the active neutrino $m_{1}$, as functions of the neutrino oscillation parameters $\theta_{12}$, $\theta_{23}$, $\theta_{13}$, $\Delta m^{2}_{21}$, and $\Delta m^{2}_{3\ell}$. We analyze these conditions, following Refs. [64, 65]. From Eqs. (17) and (18), we obtain

[TABLE]

with

[TABLE]

where we have used $\widetilde{V}^{T}=V^{-1}$ and $\textrm{det}V=1$ with $\widetilde{V}$ being the cofactor matrix of $V$.444From $\tilde{V}^{T}=V^{-1}$, we obtain $V^{*}=\tilde{V}$ , and for instance, the following relation holds: $V^{*}_{e1}=V_{\mu 2}V_{\tau 3}-V_{\mu 3}V_{\tau 2}$. In Appendix A.1, we give the explicit expressions for $R_{2}$ and $R_{3}$ in terms of neutrino oscillation parameters. By taking absolute values of Eqs. (20), we obtain the mass ratios $m_{2}/m_{1}$, $m_{3}/m_{1}$ as functions of the Dirac CP phase:

[TABLE]

In Fig. 1, we show the mass ratios $m_{2}/m_{1}$ and $m_{3}/m_{1}$ against the Dirac CP phase $\delta$.

The bands show uncertainty coming from the 1$\sigma$ error in $\theta_{12}$, which gives dominant contributions. For input parameters of the neutrino mixing angles, we use the NO values in Tab. 4. We find that the NO: $m_{1}<m_{2}<m_{3}$ is realized in $0\leq\delta\leq 2\pi$ from this figure, while the IO is not.

We have calculated the mass ratios in the other U(1)${}_{Y^{\prime}}$ extended models in the same way as the above example, and then have found that only five U(1)${}_{Y^{\prime}}$ symmetries realize the NO and two do the IO555This result that five U(1)${}_{Y^{\prime}}$ symmetries realize the NO or the IO is consistent with Refs. [58, 59].. We show the summary of the consistency between the neutrino oscillation parameters and mass ordering in Tab. 5.

“$\checkmark$” represents the compatibility of the TZM structure with the NO or IO when the best fit value of $\theta_{23}$ is used, and “$\times$” does the opposite. “$\triangle$” in the $\mathbf{B_{3}^{R}}$ and $\mathbf{B_{4}^{R}}$ cases mean that the NO or IO cannot be realized when we use the best fit value of $\theta_{23}$, but they are realized when we use the smallest value in the $3\sigma$ region of $\theta_{23}$. From this result, we discuss only the U(1)${}_{L_{\mu}-L_{\tau}}[\mathbf{C^{R}}]$, U(1)${}_{B-L_{e}-3L_{\mu}+L_{\tau}}[\mathbf{B_{3}^{R}}]$, U(1)${}_{B-L_{e}+L_{\mu}-3L_{\tau}}[\mathbf{B_{4}^{R}}]$, U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}[\mathbf{D_{1}^{R}}]$ and U(1)${}_{B+L_{e}-L_{\mu}-3L_{\tau}}[\mathbf{D_{2}^{R}}]$ cases, henceforth.

Returning to the U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}$ [$\mathbf{D_{1}^{R}}$] case, from the above discussion, we fix the mass ordering to the NO and determine $m_{1}$ and $\delta$ as functions of neutrino oscillation parameters. From Eqs. (22), an equation for $\cos\delta$ is obtained by

[TABLE]

where $\epsilon\equiv\Delta m_{21}^{2}/\Delta m_{31}^{2}$, and it leads $\delta$ as functions of neutrino oscillation parameters. We give an explicit expression of Eq. (23) in Eq. (A.2) in Appendix A.2. Substituting the best fit values of the mixing angles $\theta_{ij}$ in Tab. 4 into Eq. (A.2) and solving it numerically, we obtain $\cos\delta\simeq 0.102$ for the best fit values of the neutrino oscillation parameters, which means $\delta=0.468\pi$ or $1.532\pi$.

Substituting this result into Eqs. (22), we also obtain $m_{1}$ as a function of neutrino oscillation parameters. An explicit expression for $m_{1}$ is given in Appendix A.3. Moreover, comparing the phases of $R_{2,3}$ in Eqs. (20) respectively, we can obtain the Majorana CP phases:

[TABLE]

It is shown in Ref. [64] that $|R_{2}(\delta)|$ and $|R_{3}(\delta)|$ do not change under $\delta\to-\delta$, and so $m_{1}$ also does not change whether we choose $\delta$ or $-\delta$. Furthermore, we have $\alpha_{2,3}(-\delta)=-\alpha_{2,3}(\delta)$. These discussions are applicable to the other U(1)${}_{Y^{\prime}}$ symmetries and we obtain the same conclusions.

4 Predictions for the neutrino parameters

In the previous section, we have shown that, by analyzing the two-zero minor and texture structures, we can obtain the neutrino mass, Dirac CP phase and Majorana CP phases as functions of the neutrino oscillation parameters, namely the three mixing angles and two squared mass differences. Here, we show the results of calculations with the errors in the neutrino oscillation parameters. For input values, we use the values given in Tab. 4 [104, 105]. Moreover, we show the result of the effective neutrino mass relevant to the neutrinoless double beta decay. The results in the U(1)${}_{L_{\mu}-L_{\tau}}$ case are not new ones, but ones of our previous work [65]. However, we show them in this paper for completeness.

4.1 Neutrino masses

In the singlet cases, the U(1)${}_{Y^{\prime}}$ gauge symmetries which realize the NO or the IO are only five, as shown in Tab. 5. In Fig. 2, we plot the sum of the neutrino masses as functions of $\theta_{23}$ in the red line. $\theta_{23}$ is varied in the $3\sigma$ region. Fig. 2 (2(a)) is the same as Fig. 1 in Ref. [65]. The vertical gray dashed line represents the best fit value of $\theta_{23}$, and the vertical gray dotted lines (the plot range) indicate the $1\sigma$ ($3\sigma$) region. The dark (light) red bands show the uncertainty coming from the $1\sigma$ ($3\sigma$) errors of $\theta_{13}$ and $\Delta m^{2}_{31}$ in the $\mathbf{C^{R}}$ and the other structures, respectively. The contribution to $\sum m_{i}$ from the errors of the other parameters are subdominant, and so we take them to be the best fit values. The horizontal dashed line shows the present limit on the sum of the neutrino masses by the Planck experiment: $\sum_{i}m_{i}<0.12$ eV (Planck TT+lowP+lensing+ext) [107].

As we see in Fig. 2 (2(a)), in the $\mathbf{C^{R}}$ case, there is a strong tension between the prediction and Planck limit, and the allowed region of $\theta_{23}$ is only around $\theta_{23}\simeq 52^{\circ}$. For such a $\theta_{23}$, the sum of the neutrino masses $\sum_{i}m_{i}$ is so large that it implies the quasi-degenerate mass spectrum. In both the $\mathbf{D_{1}^{R}}$ and $\mathbf{D_{2}^{R}}$ cases, $\sum_{i}m_{i}$ is about $0.065$ eV and lighter than the limit. But Eq. (23) has solutions for $\theta_{23}\gtrsim 44^{\circ}$ in the $\mathbf{D_{1}^{R}}$ case, and for $\theta_{23}\lesssim 46^{\circ}$ in the $\mathbf{D_{2}^{R}}$ case, as we see in Figs. 2 (2(b)) and (2(c)). On the other hand, as we see in Figs. 2 (2(d)), (2(e)), (2(f)) and (2(g)), in these cases, the predicted values of $\sum_{i}m_{i}$ conflict with the limit when we allow the parameters to be varied in $3\sigma$. Therefore these four cases are excluded.

From the above results, in the singlet cases, only three lepton flavor-dependent U(1) gauge symmetries, U(1)${}_{L_{\mu}-L_{\tau}}~{}[\mathbf{C^{R}}]$, U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}~{}[\mathbf{D_{1}^{R}}]$ and U(1)${}_{B+L_{e}-L_{\mu}-3L_{\tau}}~{}[\mathbf{D_{2}^{R}}]$, realize the NO without conflicting with the existing experiments and no symmetry does the IO. Furthermore, as also shown in Ref. [65], the $\mathbf{C^{R}}$ case has strong tension with the Planck 2018 limit and will soon be tested in the future neutrino experiments.

On the other hand, in the doublet cases, only the U(1)${}_{L_{\mu}-L_{\tau}}$ gauge symmetry realizes the two-zero texture structure and the possible mass orderings without conflicting with the various experiments relevant to $Z^{\prime}$. This case has been discussed in the previous work [65] and the results are shown there, but for completeness, we show them here.

In Fig. 3, which is the same as Fig. 2 in Ref. [65], we plot the sum of the neutrino masses in the doublet models as functions of $\theta_{23}$ in the red line. The dark (light) red bands show the uncertainty coming from the $1\sigma$ ($3\sigma$) errors of $\Delta m_{31}^{2}$. The contribution to $\sum m_{i}$ from the errors of the other parameters are subdominant, and so we take them to be the best fit values.

Comparing Fig. 3 with Fig. 2, we find that the predicted values in the $\mathbf{B_{3}^{\nu}}$ ($\mathbf{B_{4}^{\nu}}$) cases are the same as those in the $\mathbf{B_{4}^{R}}$ ($\mathbf{B_{3}^{R}}$). This is because the inverse of the matrix which has the $\mathbf{B_{3}}$ structure has the $\mathbf{B_{4}}$ structure and vice versa.

[TABLE]

In any case, as we have concluded in the previous work [65], all the doublet cases result in too heavy neutrino masses and they are excluded by the Planck 2018 limit. Therefore, we show the prediction of only the singlet cases hereafter.

4.2 Dirac CP phase

In Fig. 4, we plot the Dirac CP phase as functions of $\theta_{23}$ in the red line. $\theta_{23}$ is varied in the $3\sigma$ region. Fig. 4 (4(a)) is the same as Fig. 3 (a) in Ref. [65]. The dark (light) red bands show the uncertainty coming from the $1\sigma$ ($3\sigma$) errors of $\theta_{12}$. We notice that the dominant contribution to the uncertainty comes from the error in $\theta_{12}$. On the other hand, the other parameters give subdominant contribution, and so we take them to be the best fit values. We also show the $1\sigma$ ($3\sigma$) favored region of $\delta$ in the dark (light) horizontal green bands.

As we see in Fig. 4 (4(a)), in the $\mathbf{C^{R}}$ case, the upper band of the predicted values is in the experimentally favored region. For $\theta_{23}\simeq 52^{\circ}$, the prediction is $\delta\simeq 1.3\pi$ and it is in the $1\sigma$ range of the favored region of $\delta$. In the $\mathbf{D_{1}^{R}}$ and $\mathbf{D_{2}^{R}}$ cases, the predictions depend on $\theta_{12}$ as much as $\theta_{23}$ and the other parameters have considerable contributions. Moreover, the quadratic equations of $\cos\delta$ have no solution in some regions of $\theta_{23}$ and $\mathbf{D_{2}^{R}}$ cannot realize the NO when we take all the neutrino oscillation parameters to be the best fit values.

4.3 Effective Majorana neutrino mass

Lastly, we discuss neutrinoless double-beta decay ($0\nu\beta\beta$) and the possibility of testing the singlet models666See also Ref. [61].. The rate of neutrinoless double-beta decay is proportional to the square of the effective Majorana neutrino mass $\braket{m_{\beta\beta}}$, which is defined by

[TABLE]

As we note in the previous section, not only the neutrino mass $m_{1}$ and Dirac CP phase $\delta$, but also the Majorana CP phases $\alpha_{2,3}$ are uniquely determined as functions of the neutrino oscillation parameters $\theta_{12}$, $\theta_{23}$, $\theta_{13}$, $\Delta m_{21}^{2}$ and $\Delta m_{31}^{2}$. Thus we can calculate the effective Majorana neutrino mass $\braket{m_{\beta\beta}}$ without ambiguity. In Fig. 5, we show the predicted value of $\braket{m_{\beta\beta}}$, in the $\mathbf{C^{R}}$, as a function of $\theta_{23}$ in the red line. This figure is also shown in Ref. [65]. $\theta_{23}$ is varied in the $3\sigma$ region. The dark (light) red bands show the uncertainty coming from the $1\sigma$ ($3\sigma$) errors of the other parameters. We also show in the light blue band the current bound on $\braket{m_{\beta\beta}}$ given by the KamLAND-Zen experiment, $\braket{m_{\beta\beta}}<0.061\mathchar 45\relax 0.165$ eV [3], where the uncertainty stems from the estimation of the nuclear matrix element for 136Xe.

As we see in Fig. 5, $\braket{m_{\beta\beta}}$ is predicted to be $\simeq 0.016$ eV for $\theta_{23}\simeq 52^{\circ}$, which is well below the present KamLAND-Zen limit. Future experiments are expected to have sensitivities as low as $\mathcal{O}(0.01)$ eV [108], and thus it is expected that the $\mathbf{C^{R}}$ case will be tested in the future.

On the other hand, the $\mathbf{D_{1}^{R}}$ and $\mathbf{D_{2}^{R}}$ cases predict that $\braket{m_{\beta\beta}}$ vanishes and it means that neutrinoless double-beta decay process never occurs. This is because when the mass matrix of the neutrinos has the $\mathbf{D_{1}^{R}}$ ($\mathbf{D_{2}^{R}}$) structure, that matrix also has the $\mathbf{A_{2}^{\nu}}$ ($\mathbf{A_{1}^{\nu}}$) structure, respectively, and then ($e,e$) component of the neutrino mass matrix vanishes. Thus the following relation is satisfied:

[TABLE]

Therefore, the existing experimental limits on $\braket{m_{\beta\beta}}$ give no constraint to the $\mathbf{D_{1}^{R}}$ and $\mathbf{D_{2}^{R}}$ cases. On the other hand, if neutrinoless double-beta decay processes are detected in future experiments, these cases are excluded.

5 Implications for Leptogenesis

The minimal gauged U(1)${}_{Y^{\prime}}$ models contain 3 right-handed neutrinos which couple to the SM leptons. Therefore, if they are heavy enough, the observed baryon asymmetry of the Universe can be explained by the leptogenesis scenario [109], which is one of the most promising mechanisms.

In the previous section, we have shown that the neutrino masses, the Dirac CP phase, and Majorana CP phases are determined as functions of the neutrino oscillation parameters in the minimal gauged U(1)${}_{Y^{\prime}}$ models. This smallness of the degree of freedom affects the parameters relevant to the baryon asymmetry generated through leptogenesis.

In this section, following the analyses of Ref. [64], we firstly show the predicted signs of the asymmetry parameter $\epsilon$, which is one of the most important parameters in leptogenesis scenario, and then discuss whether the correct sign of the baryon asymmetry can be generated in the minimal gauged U(1)${}_{Y^{\prime}}$ models. Next, we briefly discuss the possibility of generating enough amount of the baryon asymmetry in the thermal leptogenesis scenario. For simplicity, we do not take into account flavor effects [110, 111, 112].

In the minimal gauged U(1)${}_{Y^{\prime}}$ model, the absolute neutrino mass $m_{1}$, Dirac CP phase $\delta$ and Majorana CP phases $\alpha_{2,3}$ are not free parameters and determined as functions of the neutrino oscillation parameters $\theta_{12},\theta_{23},\theta_{13},\Delta m_{21}^{2},\Delta m_{31}^{2}$ and sign of $\sin\delta$. Then, from the seesaw formula: $\mathcal{M}_{\nu_{L}}\simeq-\mathcal{M}_{D}\mathcal{M}_{R}^{-1}\mathcal{M}_{D}^{T}$, the Majorana mass matrix $\mathcal{M}_{R}$ is tightly constrained, and the free parameters in $\mathcal{M}_{R}$ are only the three Dirac Yukawa couplings of the neutrinos. By diagonalizing the masses of the right-handed neutrinos, we can rewrite the Lagrangian (11) as follows:

[TABLE]

where $\hat{N}_{i}^{c}$ are the right-handed neutrinos in the basis where the Majorana mass matrix $\mathcal{M}_{R}$ is diagonal and $\hat{\lambda}_{i\alpha}$ are the Dirac Yukawa couplings of the neutrinos in that basis, and $M_{i}$ are the masses which are taken to be real and positive. These values are written by

[TABLE]

Here, for simplicity, we assume that the heavy neutrinos have mass hierarchy and the lightest right-handed neutrino $N_{1}$ is much lighter than the others. In this case, the lepton asymmetry is generated by the decay of the lightest right-handed neutrino, and then the asymmetry parameter $\epsilon_{1}$ is given, at the leading order, by [113, 114, 115]

[TABLE]

As shown in Ref. [64], this asymmetry parameter have a correlation with the sign of $\sin\delta$. This is because that the Majorana CP phases flip: $\alpha_{2,3}\to-\alpha_{2,3}$ under sign flipping of the Dirac CP phase: $\delta\to-\delta$, and then the PMNS matrix transforms as $U_{\text{PMNS}}\to U_{\text{PMNS}}^{*}$. On the other hand, the absolute neutrino mass $m_{1}$ depends on only $\cos\delta$, not $\sin\delta$, and so $m_{1}$ does not change under sign flipping of $\delta$. Then $\mathcal{M}_{\nu_{L}}$, $\mathcal{M}_{R}$, $\Omega$, and $\hat{\lambda}$ also transform as $\mathcal{M}_{\nu_{L}}\to\mathcal{M}_{\nu_{L}}^{*}$, $\mathcal{M}_{R}\to\mathcal{M}_{R}^{*}$, $\Omega\to\Omega^{*}$, and $\hat{\lambda}\to\hat{\lambda}^{*}$. Therefore, from Eq. (47), the sign of the asymmetry parameter $\epsilon_{1}$ is flipped: $\epsilon_{1}\to-\epsilon_{1}$, and it means that there is a one-to-one correspondence between the sign of the Dirac CP phase and that of the baryon asymmetry in the Universe. In the following discussion, we take the Dirac CP phase to be $\delta>\pi$, which is favored by the experiments.

When we focus on only the sign of the baryon asymmetry, only the ratio between the Yukawa couplings $\lambda_{e}$, $\lambda_{\mu}$, $\lambda_{\tau}$ are important. Now we parametrize the Yukawa couplings as follows:

[TABLE]

For this parametrization, the asymmetry parameter in Eq. (47) is written by

[TABLE]

From Eq. (50), the sign of $\epsilon_{1}$ is independent of the scale of the Yukawa couplings $\lambda$, and it is determined by the ratios $\theta$, $\phi$. Note that the sphaleron process predicts $n_{B}/n_{L}<0$ [116], and then $\epsilon_{1}<0$ leads to the observed baryon asymmetry of the Universe.

In Fig. 6, we show the parameter region of the $(\theta,\phi)$ plane where the asymmetry parameter $\epsilon_{1}$ is negative. In Figs. 6 (6(a)), (6(b)) and (6(c)), we use $\theta_{23}=52^{\circ},49.7^{\circ}\text{(best fit value)},43^{\circ}$, respectively, and take the other parameters to be the best fit values. The red shaded areas show the parameter region of the $(\theta,\phi)$ plane where the asymmetry parameter $\epsilon_{1}$ is negative: $\epsilon_{1}<0$. Moreover, we also show the contours of the right-handed neutrino mass ratio $M_{2}/M_{1}$.

These figures show that the $\mathbf{C^{R}}$ and $\mathbf{D_{1}^{R}}$ cases predict $\epsilon_{1}<0$ in the large regions of the $(\theta,\phi)$ plane. On the other hand, $\mathbf{D_{2}^{R}}$ case realizes $\epsilon_{1}<0$ only on the thin red line: $\phi\simeq\pi/4$ & $\theta\lesssim\pi/4$ ($\lambda_{\mu}\simeq\lambda_{\tau}<\lambda_{e}$) and in the red triangle: $\theta>\pi/4$.

For the estimation of the final baryon asymmetry generated by the leptogenesis scenario, the production mechanism of the right-handed neutrinos and the efficiency factor $\kappa_{f}$ are important. In the case of thermal leptogenesis, there are weak and strong washout regimes. The $\mathbf{C^{R}}$ case predicts $m_{1}\gtrsim 0.03$ eV, and the $\mathbf{D_{1}^{R}}$ and $\mathbf{D_{2}^{R}}$ cases do $m_{1}\gtrsim 0.004$ eV, and then they are classified in the strong washout regime777See , e.g., [117]. In that regime, the efficiency factor $\kappa_{f}$ is given by [117]

[TABLE]

where $\tilde{m}_{1}\equiv(\mathcal{M}^{\dagger}\mathcal{M}_{D})_{11}/M_{1}$ is the effective neutrino mass. In general, the relation $\tilde{m}_{1}\geq m_{1}$ holds [118], and then $\kappa_{f}$ are roughly given as follows;

[TABLE]

The generated baryon asymmetry is given by [117]

[TABLE]

where $g_{*}$ is the degree of freedom of relativistic fields. Therefore, the asymmetry parameter needs to be $\epsilon_{1}\sim\mathcal{O}(10^{-5}\mathchar 45\relax 10^{-6})$ to explain the observed baryon asymmetry $Y_{\Delta B}\sim 8.7\times 10^{-11}$ [107].

In Figs. 7, we show the sizes of the asymmetry parameters $-\epsilon_{1}/\lambda^{2}$ as functions of $\phi$. As you see in Eq. (50), for fixed $\theta$, the asymmetry parameter $\epsilon_{1}$ scales as $\epsilon_{1}\propto\lambda^{2}$ in this parametrization, and so we normalize $\epsilon_{1}$ by $\lambda^{2}$. The neutrino oscillation parameters except for $\theta_{23}$ are taken to be the best fit values, and $\theta_{23}$ are taken to be $\theta_{23}=52^{\circ},49.7^{\circ}\text{(best fit value)},43^{\circ}$, respectively. The sign of the Dirac CP phase is taken to be negative ($\delta>\pi$).

In this setup, these figures show that the $\mathbf{C^{R}}$ and $\mathbf{D_{1}^{R}}$ cases predict sizable values $|\epsilon_{1}|>10^{-5}$. On the other hand, it is difficult for the $\mathbf{D_{2}^{R}}$ case to give a sizable negative value, and only for $2\lesssim\tan\theta\lesssim 200$, that is, in the red triangle area in Fig. 6 (6(c)), the predicted value of $-\epsilon_{1}/\lambda^{2}$ is greater than $10^{-6}$. In a part of this area, as we see in Fig. 6 (6(c)), the right-handed neutrino masses $M_{1},M_{2}$ are degenerate, and so we have to consider the resonant leptogenesis [119, 120, 121, 122] for more correct discussion.

6 Conclusions

In this work, we have extended U(1)${}_{L_{\alpha}-L_{\beta}}$ gauge symmetries in Refs. [64, 65] to more general lepton flavor-dependent U(1)${}_{Y^{\prime}}$, namely linear combinations of U(1)${}_{L_{e}-L_{\mu}}$, U(1)${}_{L_{\mu}-L_{\tau}}$ and U(1)B-L, and have studied the two-zero minor and texture structures in the minimal gauged U(1)${}_{Y^{\prime}}$ models. In these models, because of the gauge symmetries, the structures of the Dirac and Majorana mass matrices are tightly restricted. These restrictions connect the low-energy neutrino parameters each other, and four of them are given as functions of the rest of the parameters, namely the neutrino mixing angles and squared mass differences. Using these relations, we obtained the prediction for the Dirac CP phase, sum of the neutrino masses, and effective Majorana neutrino mass. Furthermore, among the minimal U(1)${}_{Y^{\prime}}$ gauged models which realize the TZM or TZT, we found that only three models with one U(1)${}_{Y^{\prime}}$-breaking singlet scalar for NO are not in conflict with the experimental constraints. We summarize the results in Tab. 6. “N” and “I” represent the NO and IO, respectively. “$\checkmark$” and “$\times$” represent whether models are consistent with the experimental constraints or not. “$\triangle$” means that, in the column of mixing, the models predict the NO or the IO if $\theta_{23}$ is taken to be in the range of $3\sigma$, and in that of mass, there is a strong tension between the prediction and the Planck limit. Moreover, in that of leptogenesis, it means that the limited parameter regions realize the correct sign and enough size of the baryon asymmetry. “$\square$” means that the models predict that the $0\nu\beta\beta$ process never occurs.

Among them, U(1)${}_{L_{\mu}-L_{\tau}}[\mathbf{C^{R}}]$ case predicts heavy neutrino masses and has a strong tension with the Planck 2018 limit. Therefore, this case will soon be tested by the future experiments of $\sum_{i}m_{i}$ and $\theta_{23}$. On the other hand, U(1)${}_{B+L_{e}-3L_{\mu}-L_{\tau}}[\mathbf{D_{1}^{R}}]$ and U(1)${}_{B+L_{e}-L_{\mu}-3L_{\tau}}[\mathbf{D_{2}^{R}}]$ cases predict lighter $\sum_{i}m_{i}$ than the Planck limit and no neutrinoless double beta decay. Then if neutrinoless double beta decay is observed in the future, these cases will be excluded completely.

We have also discussed the implications of the minimal gauged U(1)${}_{Y^{\prime}}$ models for the leptogenesis scenario. Because of the two-zero minor conditions, the neutrino Majorana mass matrix in the minimal models is severely restricted and have only three free parameters. We found that the $\mathbf{C^{R}}$ and $\mathbf{D_{1}^{R}}$ cases generate the correct sign of the baryon asymmetry through the leptogenesis scenario in the large parameter spaces, and predict enough size of the asymmetry parameters in the thermal leptogenesis scenario. On the other hand, the $\mathbf{D_{2}^{R}}$ case generates the correct sign of the baryon asymmetry only in the limited regions of the parameter space, and the observed value of the baryon asymmetry can be generated in this case when the right-handed neutrino masses are degenerate.

Acknowledgments

We thank Koichi Hamaguchi and Natsumi Nagata for valuable discussions and suggestions. This work was supported by JSPS KAKENHI Grant Number 19J13812.

Data Availability Statement

This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The values of the neutrino oscillation parameters used during this analysis are shown in Tab. 4, and they can be obtained from Ref. [104, 105]]

Appendix

Appendix A Miscellaneous formulae

Here we give formulae that are useful for the study of the neutrino mass structures in three minimal gauged U(1)${}_{Y^{\prime}}$ models which do not conflict with the experimental results, such as the neutrino oscillation experiments and the Planck 2018. The formulae in the $\mathbf{C^{R}}$ case are also shown in Ref. [64].

A.1 $R_{2}$ and $R_{3}$

Here, we show the list of the functions relevant to the ratios $m_{2,3}/m_{1}$, namely $R_{2,3}(\delta)$ in Tab. 7.

A.2 Equation for $\cos\delta$

In this subsection, we show cubic and quadratic equations whose real solution in terms of $x$ gives $\cos\delta$.

In the $\mathbf{C^{R}}$ case, the Dirac CP phase can be obtained from the following cubic equation:

[TABLE]

where

[TABLE]

In the $\mathbf{D_{1}^{R}}$ case, the Dirac CP phase can be obtained from the following quadratic equation:

[TABLE]

This equation can be solved easily and one of the solutions is certainly greater than one. The other solution can give the Dirac CP phase $\delta$ if, and only if it satisfies $-1\leq x\leq 1$:

[TABLE]

In the $\mathbf{D_{2}^{R}}$ case, the Dirac CP phase can be obtained from the following quadratic equation:

[TABLE]

This equation can be solved easily and one of the solutions is certainly greater than one. The other solution can give the Dirac CP phase $\delta$ if, and only if it satisfies $-1\leq x\leq 1$:

[TABLE]

A.3 Neutrino mass $m_{1}$

Here, we show the list of the formulae of the neutrino mass $m_{1}$ in Tab. 8.