Minimal Dirac Neutrino Mass Models from $U(1)_R$ Gauge Symmetry and Left-Right Asymmetry at Colliders

TL;DR

This paper introduces minimal models for Dirac neutrino masses using a $U(1)_R$ gauge symmetry, predicting distinctive collider signatures and linking neutrino mass with dark matter, while analyzing experimental bounds and phenomenology.

Contribution

It presents novel minimal $U(1)_R$-based models for Dirac neutrino masses, including a radiative scotogenic model linking neutrinos and dark matter, with distinctive collider signatures.

Findings

Predicted sizable Left-Right asymmetry at colliders.

Analyzed discovery potential of $Z'$ at LHC and ILC.

Explored dark matter stability and lepton flavor violation.

Abstract

In this work, we propose minimal realizations for generating Dirac neutrino masses in the context of a right-handed abelian gauge extension of the Standard Model. Utilizing only $U(1)_R$ symmetry, we address and analyze the possibilities of Dirac neutrino mass generation via (a) \textit{tree-level seesaw} and (b) \textit{radiative correction at the one-loop level}. One of the presented radiative models implements the attractive \textit{scotogenic} model that links neutrino mass with Dark Matter (DM), where the stability of the DM is guaranteed from a residual discrete symmetry emerging from $U(1)_R$. Since only the right-handed fermions carry non-zero charges under the $U(1)_R$, this framework leads to sizable and distinctive Left-Right asymmetry as well as Forward-Backward asymmetry discriminating from $U(1)_{B-L}$ models and can be tested at the colliders. We analyze the current experimental bounds and present the discovery reach limits for the new heavy gauge boson $Z^{\prime}$ at the LHC and ILC. Furthermore, we also study the associated charged lepton flavor violating processes, dark matter phenomenology and cosmological constraints of these models.

Tables

Table 1. Table I : Quantum numbers of the fermions and the SM Higgs doublet.

Multiplets	$SU{\left(3\right)}_C\times SU{\left(2\right)}_L\times U{\left(1\right)}_Y\times U{\left(1\right)}_R$
Quarks	\pbox10cm $Q_L^{}{}_i{}^{}(3,2,\frac{1}{6},0)$
$u_R^{}{}_i{}^{}(3,1,\frac{2}{3},R_H)$
$d_R^{}{}_i{}^{}(3,1,-\frac{1}{3},\left\{-R_H\right\})$
Leptons	\pbox10cm $L_L^{}{}_i{}^{}(1,2,-\frac{1}{2},0)$
$\ell _R^{}{}_i{}^{}(1,1,-1,-R_H)$
$\nu _R^{}{}_i{}^{}(1,1,0,\{R_{{\nu }_1},R_{{\nu }_2},R_{{\nu }_3}\})$
Higgs	\pbox10cm $H(1,2,\frac{1}{2},R_H)$

Table 2. Table II : Couplings of the fermions with the new gauge boson. Here we use the notation: c 2 ​ w = cos ⁡ ( 2 ​ θ w ) subscript 𝑐 2 𝑤 2 subscript 𝜃 𝑤 c_{2w}=\cos(2\theta_{w}) . 𝒩 L , R subscript 𝒩 𝐿 𝑅 \mathcal{N}_{L,R} is any vector-like fermion singlet under the SM and carries R 𝒩 subscript 𝑅 𝒩 R_{\mathcal{N}} charge under U ​ ( 1 ) R 𝑈 subscript 1 𝑅 U(1)_{R} . If a model does not contain vector-like fermions, we set R 𝒩 = 0 subscript 𝑅 𝒩 0 R_{\mathcal{N}}=0 .

Fermion, $\psi$	Coupling, $g_{\psi }$
Quarks	\pbox10cm $g_{u_L}=-\frac{1}{6}\frac{g}{c_w}\left(1+2c_{2w}\right)s_X$
$g_{d_L}=\frac{1}{6}\frac{g}{c_w}\left(2+c_{2w}\right)s_X$
$g_{u_R}=\frac{2}{3}\frac{g}{c_w}{s}_w^2{s}_X+g_R{c}_X{R}_H$
$g_{d_R}=-\frac{1}{3}\frac{g}{c_w}{s}_w^2{s}_X-g_R{c}_X{R}_H$
Leptons	\pbox10cm $g_{{\nu }_L}=-\frac{1}{2}\frac{g}{c_w}{s}_X$
$g_{{\ell }_L}=\frac{1}{2}\frac{g}{c_w}{c}_{2w}{s}_X$
$g_{{\ell }_R}=-\frac{g}{c_w}{s}_w^2{s}_X-g_R{c}_X{R}_H$
$g_{{\nu }_{R_i}}=g_R{c}_X{R}_{{\nu }_i}$
Vector-like fermions	\pbox10cm $g_{𝒩}=g_R{c}_X{R}_{𝒩}$

Table 3. Table III : Quantum numbers of the fermions and the scalars in Dirac seesaw model.

Multiplets	$SU{\left(3\right)}_C\times SU{\left(2\right)}_L\times U{\left(1\right)}_Y\times U{\left(1\right)}_R$
Leptons	\pbox10cm $L_L^{}{}_i{}^{}(1,2,-\frac{1}{2},0)$
$\ell _R^{}{}_i{}^{}(1,1,-1,-1)$
$\nu _R^{}{}_i{}^{}(1,1,0,\{-5,4,4\})$
Scalars	\pbox10cm $H(1,2,\frac{1}{2},1)$
$\chi (1,1,0,3)$
Vector-like fermion	\pbox10cm ${𝒩}_{L,R}(1,1,0,1)$

Table 4. Table IV : Quantum numbers of the fermions and the scalars in radiative Dirac model.

Multiplets	$SU{\left(3\right)}_C\times SU{\left(2\right)}_L\times U{\left(1\right)}_Y\times U{\left(1\right)}_R$
Leptons	\pbox10cm $L_L^{}{}_i{}^{}(1,2,-\frac{1}{2},0)$
$\ell _R^{}{}_i{}^{}(1,1,-1,-1)$
$\nu _R^{}{}_i{}^{}(1,1,0,\{-5,4,4\})$
Scalars	\pbox10cm $H(1,2,\frac{1}{2},1)$
$\chi (1,1,0,3)$
$S_1^{+}(1,1,1,0)$
$S_2^{+}(1,1,1,-3)$

Table 5. Table V : Quantum numbers of the fermions and the scalars in scotogenic Dirac neutrino mass model.

Multiplets	$SU{\left(3\right)}_C\times SU{\left(2\right)}_L\times U{\left(1\right)}_Y\times U{\left(1\right)}_R$
Leptons	\pbox10cm $L_L^{}{}_i{}^{}(1,2,-\frac{1}{2},0)$
$\ell _R^{}{}_i{}^{}(1,1,-1,-1)$
$\nu _R^{}{}_i{}^{}(1,1,0,\{-5,4,4\})$
Scalars	\pbox10cm $H(1,2,\frac{1}{2},1)$
$\chi (1,1,0,3)$
$S(1,1,0,-\frac{7}{2})$
$\eta (1,2,\frac{1}{2},\frac{1}{2})$
Vector-like fermion	\pbox10cm ${𝒩}_{L,R}(1,1,0,\frac{1}{2})$

Equations

$${\mathcal{L}}_Y\supset y_u{\overline{Q}}_L\tilde{H}{u}_R+y_d{\overline{Q}}_{L}H{d}_R+y_e{\overline{L}}_{L}H{\ell }_R+h.c.$$

$$R_u=-R_d=-R_{\ell }=R_H.$$

$$[U(1{)}_R]:R_{{\nu }_1}+R_{{\nu }_2}+R_{{\nu }_3}=3R_H,$$

$$[U(1{)}_R{]}^3:R_{{\nu }_1}^3+R_{{\nu }_2}^3+R_{{\nu }_3}^3=3R_H^3.$$

$$D_{\mu }H=\left({\partial }_{\mu }-ig{W}_{\mu }-i{g}^{\prime }{Y}_H{B}_{\mu }-i{g}_R{R}_H{X}_{\mu }\right)H,$$

$$D_{\mu }\chi =\left({\partial }_{\mu }-i{g}_R{R}_{\chi }{X}_{\mu }\right)\chi .$$

$$⟨H⟩=\left(\begin{array}{c}0\\ \frac{v_H}{\sqrt{2}}\end{array}\right),⟨\chi ⟩=\frac{v_{\chi }}{\sqrt{2}},$$

$$\left(\begin{array}{ccc}B& W_3& X\end{array}\right)\left(\frac{v_H^2}{4}\right)\left(\begin{array}{ccc}{g}^{\prime 2}& -g^{\prime }g& 2g^{\prime }{g}_R{R}_H\\ -g^{\prime }g& g^2& -2g{g}_R{R}_{\chi }\\ 2g^{\prime }{g}_R{R}_H& -2g{g}_R{R}_{\chi }& 4g_R^2{R}_H^2(1+r_v^2)\end{array}\right)\left(\begin{array}{c}B\\ W_3\\ X\end{array}\right).$$

$$M_{Z,Z^{\prime }}=\frac{g{v}_H}{2c_w}{\left(\frac{1}{2}\left[1+r_X^2{c}_w^2(1+r_v^2)\right]\mp \left[\frac{r_X{c}_w}{\sin (2{\theta }_X)}\right]\right)}^{\frac{1}{2}},$$

$$r_X=(2g_R{R}_H)/g,$$

$$\sin (2{\theta }_X)=\frac{2r_X{c}_w}{{\left({\left[2r_X{c}_w\right]}^2+{\left[(1+r_v^2)r_X^2{c}_w^2-1\right]}^2\right)}^{\frac{1}{2}}}.$$

$$\left(\begin{array}{c}B\\ W_3\\ X\end{array}\right)=U\left(\begin{array}{c}A\\ Z\\ Z^{\prime }\end{array}\right),U=\left(\begin{array}{ccc}{c}_w& -s_w{c}_X& s_w{s}_X\\ s_w& c_w{c}_X& -c_w{s}_X\\ 0& s_X& c_X\end{array}\right).$$

$$\left|\Delta M_Z\right|=\left|M_Z^{SM}\left(1-\sqrt{\frac{r_v^2}{1+r_v^2}}\right)\right|\le 2.1\mathrm{MeV}.$$

$$\mathcal{L}\supset g_{\psi }\overline{\psi }{\gamma }^{\mu }{Z}_{\mu }^{\prime }\psi .$$

$${\mathcal{L}}_Y\supset y^H{\overline{L}}_L\tilde{H}{\mathcal{N}}_R+M_{\mathcal{N}}{\overline{\mathcal{N}}}_L{\mathcal{N}}_R+y^{\chi }{\overline{\mathcal{N}}}_L{\nu }_R{\chi }^{\ast }+h.c.$$

$$V=-{\mu }_H^2{H}^{†}H+\lambda (H^{†}H{)}^2-{\mu }_{\chi }^2{\chi }^{\ast }\chi +{\lambda }_1({\chi }^{\ast }\chi {)}^2+{\lambda }_2{H}^{†}H{\chi }^{\ast }\chi .$$

$${\mathcal{L}}_Y\supset \left(\begin{array}{cc}{\overline{\nu }}_L& {\overline{\mathcal{N}}}_L\end{array}\right)M_{\nu ,\mathcal{N}}\left(\begin{array}{c}{\nu }_R\\ {\mathcal{N}}_R\end{array}\right),M_{\nu ,\mathcal{N}}=\left(\begin{array}{cc}0& \frac{v_H}{\sqrt{2}}{y}^H\\ \frac{v_{\chi }}{\sqrt{2}}{y}^{\chi }& M_{\mathcal{N}}\end{array}\right).$$

$$m_{\nu }\sim \frac{v_H{v}_{\chi }}{2}\frac{y^H{y}^{\chi }}{M_{\mathcal{N}}}.$$

$${\mathcal{L}}_Y\supset y^H{\overline{L}}_L{\ell }_{R}H+y^{S_1}\overline{L_L^c}ϵL_L{S}_1^{+}+y^{S_2}\overline{{\nu }_R^c}{\ell }_R{S}_2^{+}+h.c.$$

$$V$$

$$+{\lambda }_{\chi }({\chi }^{\ast }\chi {)}^2+{\lambda }_3|S_1^{+}{|}^2|S_2^{+}{|}^2+{\lambda }_4|S_1^{+}{|}^2{H}^{†}H+{\lambda }_5|S_2^{+}{|}^2{H}^{†}H+{\lambda }_6{H}^{†}H{\chi }^{\ast }\chi .$$

$${m_{\nu }}_{ab}=\frac{\sin (2\theta )}{16{\pi }^2}\text{ln}\left(\frac{m_{H_2}^2}{m_{H_1}^2}\right)y_{ai}^{S_1}{m_E}_i{y}_{ib}^{S_2}.$$

$$y^{\eta }{\overline{L}}_L{\mathcal{N}}_R\tilde{\eta }+M_{\mathcal{N}}{\overline{\mathcal{N}}}_L{\mathcal{N}}_R+y^S{\overline{\mathcal{N}}}_L{\nu }_{R}S+h.c.$$

$$V$$

$$+{\lambda }_1{H}^{†}H{\eta }^{†}\eta +{\lambda }_2{H}^{†}H{S}^{\ast }S+{\lambda }_3{H}^{†}H{\chi }^{\ast }\chi +{\lambda }_4{\eta }^{†}\eta S^{\ast }S+{\lambda }_5{\eta }^{†}\eta {\chi }^{\ast }\chi +{\lambda }_6{\chi }^{\ast }\chi S^{\ast }S$$

$$+({\lambda }_7{H}^{†}\eta {\eta }^{†}H+h.c.)+({\lambda }_D{\eta }^{†}H\chi S+h.c.).$$

$${m_{\nu }}_{ab}$$

$$-\frac{1}{16{\pi }^2}\frac{\sin {\theta }^{\prime }\cos {\theta }^{\prime }}{2}{y^{\eta }}_{ai}{M_{\mathcal{N}}}_i{y^S}_{ib}\left(F\left[\frac{m_{A_2^0}^2}{{M_{\mathcal{N}}^2}_i}\right]-F\left[\frac{m_{A_1^0}^2}{{M_{\mathcal{N}}^2}_i}\right]\right).$$

$$\theta =\frac{1}{2}{\sin }^{-1}\left(\frac{{\lambda }_D{v}_H{v}_{\chi }}{m_{H_2^0}^2-m_{H_1^0}^2}\right),{\theta }^{\prime }=\frac{1}{2}{\sin }^{-1}\left(\frac{{\lambda }_D{v}_H{v}_{\chi }}{m_{A_2^0}^2-m_{A_1^0}^2}\right).$$

$${\beta }_R=\frac{1}{16{\pi }^2}{b}_R{g}_R^3.$$

$$b_R=\sum _{f_i}\frac{4}{3}\kappa N_g{S}_2(f_i)+\sum _{s_i}\frac{1}{6}\eta S_2(s_i).$$

$${\Lambda }_{Landau}={\mu }_0{e}^{\frac{16{\pi }^2}{2b_R{\left(g_R({\mu }_0)\right)}^2}}.$$

$$\Gamma \left({\ell }_{\alpha }\to {\ell }_{\beta }+\gamma \right)=\frac{\alpha }{4{\left(16{\pi }^2\right)}^2}\frac{m_{\alpha }^5}{144}\left[{\left(\frac{{\cos }^2\theta }{m_{H_1}^2}+\frac{{\sin }^2\theta }{m_{H_2}^2}\right)}^2{\left|y_{i\alpha }^{S_1}{y}_{i\beta }^{S_1\ast }\right|}^2+{\left(\frac{{\sin }^2\theta }{m_{H_1}^2}+\frac{{\cos }^2\theta }{m_{H_2}^2}\right)}^2{\left|y_{i\alpha }^{S_2}{y}_{i\beta }^{S_2\ast }\right|}^2\right].$$

$$\Gamma \left(l_{\alpha }\to l_{\beta }+\gamma \right)=\frac{\alpha }{4}\frac{{\left|y_{\alpha i}^{\eta }{y}_{\beta i}^{\eta \ast }\right|}^2}{{\left(16{\pi }^2\right)}^2}\frac{{\left(m_{\alpha }^2-m_{\beta }^2\right)}^3\left(m_{\alpha }^2+m_{\beta }^2\right)}{m_{\alpha }^3{m}_{H^{+}}^4}{\left[f_B(t)\right]}^2.$$

Full text

1. 1 Introduction

2. 2 Framework

3. 3 Dirac Neutrino Mass Models

4. 3.1 Tree-level Dirac Seesaw

5. 3.2 Simplest one-loop implementation

6. 3.3 Scotogenic Dirac neutrino mass

7. 4 Running of the $U(1)_{R}$ Gauge Coupling

8. 5 Lepton Flavor Violation

9. 6 Dark Matter Phenomenology

10. 7 Collider Implications

11. 7.1 Constraint on Heavy Gauge Boson $Z^{\prime}$ from LEP

12. 7.2 Heavy Gauge Boson $Z^{\prime}$ at the LHC

13. 7.3 Heavy Gauge Boson $Z^{\prime}$ at the ILC

14. 7.3.1 Forward-Backward Asymmetry

15. 7.3.2 Left-Right Asymmetry

16. 8 Constraint from Cosmology

17. 9 Conclusions

1 Introduction

Neutrino oscillation data [1] indicates that at-least two neutrinos have tiny masses. The origin of the neutrino mass is one of the unsolved mysteries in Particle Physics. The minimal way to obtain the non-zero neutrino masses is to introduce three right-handed neutrinos that are singlets under the Standard Model (SM). Consequently, Dirac neutrino mass term at the tree-level is allowed and has the form: $\mathcal{L}_{Y}\supset y_{\nu}\overline{L}_{L}\widetilde{H}\nu_{R}$. However, this leads to unnaturally small Yukawa couplings for neutrinos ($y_{\nu}\leq 10^{-11}$). There have been many proposals to naturally induce neutrino mass mostly by using the seesaw mechanism [2, 3, 4, 5, 6] or via radiative mechanism [7]. Most of the models of neutrino mass generation assume that the neutrinos are Majorana 111For a recent review on models based on Majorana neutrinos see Ref. [8]. For Majorana neutrino mass models within the context of simple grand unified theories see Ref. [9]. type in nature. Whether neutrinos are Dirac or Majorana type particles is still an open question. This issue can be resolved by neutrinoless double beta decay experiments [10]. However, up-to-now there is no concluding evidence from these experiments.

Recently, there has been a growing interest in models where neutrinos are assumed to be Dirac particles. Many of these models use ad hoc discrete symmetries [11, 12, 13, 14, 16, 17, 18, 15, 19, 20, 21] to forbid the aforementioned unnaturally small tree-level Yukawa term as well as Majorana mass terms. However, it is more appealing to forbid all these unwanted terms utilizing simple gauge extension of the SM instead of imposing discrete or continuous global symmetries. This choice is motivated by the fact that contrary to gauge symmetries, global symmetries are known not to be respected by the gravitational interactions [22, 23, 24, 25, 26].

In this work, we extend the SM with $U(1)_{R}$ gauge symmetry, under which only the SM right-handed fermions are charged and the left-handed fermions transform trivially. This realization is very simple in nature and has several compelling features to be discussed in great details. Introducing only the three right-handed neutrinos all the gauge anomalies can be canceled and $U(1)_{R}$ symmetry can be utilized to forbid all the unwanted terms to build desired models of Dirac neutrino mass. Within this framework, by employing the $U(1)_{R}$ symmetry we construct a tree-level Dirac seesaw model [27] and two models where neutrino mass appears at the one-loop level. One of these loop models presented in this work is the most minimal model of radiative Dirac neutrino mass [28] and the second model uses the scotogenic mechanism [29] that links two seemingly uncorrelated phenomena: neutrino mass with Dark Matter (DM). As we will discuss, the stability of the DM in the latter scenario is a consequence of a residual $\mathcal{Z}_{2}$ discrete symmetry that emerges from the spontaneous breaking of the $U(1)_{R}$ gauge symmetry.

Among other simple possibilities, one can also extend the SM with $U(1)_{B-L}$ gauge symmetry [30] for generating the Dirac neutrino mass [31, 32, 33, 34, 28, 35]. Both of the two possibilities are attractive and can be regarded as the minimal gauge extensions of the SM. However, the phenomenology of $U(1)_{R}$ model is very distinctive compared to the $U(1)_{B-L}$ case. In the literature, gauged $U(1)_{B-L}$ symmetry has been extensively studied whereas gauged $U(1)_{R}$ extension has received very little attention.

Unlike the $U(1)_{B-L}$ case, in our set-up, the SM Higgs doublet is charged under this $U(1)_{R}$ symmetry to allow the desired Yukawa interactions to generate mass for the charged fermions, this leads to interactions with the new gauge boson that is absent in $U(1)_{B-L}$ model. The running of the Higgs quartic coupling gets modified due to having such interactions with the new gauge boson $Z^{\prime}$ that can make the Higgs vacuum stable [36]. Due to the same reason, the SM Higgs phenomenology also gets altered [37].

We show by detail analysis that despite their abelian nature, $U(1)_{R}$ and $U(1)_{B-L}$ have distinguishable phenomenology. The primary reason that leads to different features is: $U(1)_{R}$ gauge boson couples only to the right-handed chiral fermions, whereas $U(1)_{B-L}$ is chirality-universal. As a consequence, $U(1)_{R}$ model leads to large left-right (LR) asymmetry and also forward-backward (FB) asymmetry that can be tested in the current and future colliders that make use of the polarized initial states, such as in ILC. We also comment on the differences of our $U(1)_{R}$ scenario with the other $U(1)_{R}$ models existing in the literature. Slightly different features emerge as a result of different charge assignment of the right-handed neutrinos in our set-up for the realization of Dirac neutrino mass. In the existing $U(1)_{R}$ models, flavor universal charge assignment for the right-handed neutrinos are considered and neutrinos are assumed to be Majorana particles. Whereas, in our set-up, neutrinos are Dirac particles that demands non-universal charge assignment for the right-handed neutrinos under $U(1)_{R}$. Neutrinos being Dirac in nature also leads to null neutrinoless double beta decay signal.

The originality of this work is, by employing only the gauged $U(1)_{R}$ symmetry, we construct Dirac neutrino masses at the tree-level and one-loop level (with or without DM) which has not been done before and, by a detailed study of the phenomenology associated to the new heavy gauge boson, we show that $U(1)_{R}$ model is very promising to be discovered in the future colliders. Due to the presence of the TeV or sub-TeV scale BSM particles, these models can give rise to sizable rate for the charged lepton flavor violating processes which we also analyze. On top of that, we bring both the dark matter and the neutrino mass generation issues under one umbrella without imposing any additional symmetry and, work out the associated dark matter phenomenology. We also discuss the cosmological consequences due to the presence of the light right-handed neutrinos in our framework.

The paper is organized as follows. In Section 2, we discuss the framework where SM is extended by an abelian gauge symmetry $U(1)_{R}$. In Section 3, we present the minimal Dirac neutrino mass models in details, along with the particle spectrum and charge assignments. In Section 4, we discuss the running of the $U(1)_{R}$ coupling. Charged lepton flavor violating processes are analyzed in Section 5. We have also done the associated dark matter phenomenology in Section 6 for the scotogenic model. Furthermore, we analyze the collider implications in Section 7. In Section 8, we study the constraints from cosmological measurement and finally, we conclude in Section 9.

2 Framework

Our framework is a very simple extension of the SM: an abelian gauge extension under which only the right-handed fermions are charged. Such a charge assignment is anomalous, however, all the gauge anomalies can be canceled by the minimal extension of the SM with just three right-handed neutrinos. Within this framework the minimal choice to generate the charged fermion masses is to utilize the already existing SM Higgs doublet, hence the associated Yukawa couplings have the form:

[TABLE]

As a result, the choice of the $U(1)_{R}$ charges of the right-handed fermions of the SM must be universal and obey the following relationship:

[TABLE]

Here $R_{k}$ represents the $U(1)_{R}$ charge of the particle $k$. Hence, all the charges are determined once $R_{H}$ is fixed, which can take any value. The anomaly is canceled by the presence of the right-handed neutrinos that in general can carry non-universal charge under $U(1)_{R}$. Under the symmetry of the theory, the quantum numbers of all the particles are shown in Table I.

In our set-up, all the anomalies automatically cancel except for the following two:

[TABLE]

This system has two different types of solutions. The simplest solution corresponds to the case of flavor universal charge assignment that demands: $R_{\nu{1,2,3}}=R_{H}$ which has been studied in the literature [38, 39, 40, 41, 42]. In this work, we adopt the alternative choice of flavor non-universal solution and show that the predictions and phenomenology of this set-up can be very different from the flavor universal scenario. We compare our model with the other $U(1)_{R}$ extensions, as well as $U(1)_{B-L}$ extensions of the SM. As already pointed out, a different charge assignment leads to distinct phenomenology in our model and can be distinguished in the neutrino and collider experiments.

Since SM is a good symmetry at the low energies, $U(1)_{R}$ symmetry needs to be broken around $O(10)$ TeV scale or above. We assume that $U(1)_{R}$ gets broken spontaneously by the VEV of a SM singlet $\chi(1,1,0,R_{\chi})$ that must carry non-zero charge ($R_{\chi}\neq 0$) under $U(1)_{R}$. As a result of this symmetry breaking, the imaginary part of $\chi$ will be eaten up by the corresponding gauge boson $X_{\mu}$ to become massive. Since EW symmetry also needs to break down around the $O(100)$ GeV scale, one can compute the masses of the gauge bosons from the covariant derivatives associated with the SM Higgs $H$ and the SM singlet scalar $\chi$:

[TABLE]

As a consequence of the symmetry breaking, the neutral components of the gauge bosons will all mix with each other. Inserting the following VEVs:

[TABLE]

one can compute the neutral gauge boson masses as:

[TABLE]

Where, $r_{v}=\frac{R_{\chi}v_{\chi}}{R_{H}v_{H}}$ and the well-known relation $\tan\theta_{w}=g^{\prime}/g$ and furthermore $v_{H}=246$ GeV. In the above mass matrix denoted by $M^{2}$, one of the gauge bosons remains massless, which must be identified as the photon field, $A_{\mu}$. Moreover, two massive states appear which are the SM $Z$-boson and a heavy $Z^{\prime}$-boson ($M_{Z}<M_{Z^{\prime}}$). The corresponding masses are given by:

[TABLE]

here we define:

[TABLE]

Which clearly shows that for $g_{R}=0$, the mass of the SM gauge boson is reproduced: $M^{SM}_{Z}=\frac{1}{2}v_{H}(g^{2}+g^{\prime 2})^{1/2}=\frac{1}{2}gv_{H}/c_{w}$. To find the corresponding eigenstates, we diagonalize the mass matrix as: $M^{2}=U^{\dagger}M^{2}_{diag}U^{\ast}$, with:

[TABLE]

From Eq. (2.9) one can see that the mass of the SM $Z$-boson gets modified as a consequence of $U(1)_{R}$ gauge extension. Precision measurement of the SM $Z$-boson puts bound on the scale of the new physics. From the experimental measurements, the bound on the lower limit of the new physics scale can be found by imposing the constraint $\Delta M_{Z}\leq 2.1$ MeV [43]. For our case, this bound can be translated into:

[TABLE]

With $M^{SM}_{Z}=91.1876$ GeV [43], we find $v_{\chi}\geq\left(\frac{v_{H}R_{H}}{R_{\chi}}\right)21708.8$. Which corresponds to $v_{\chi}\geq 12.08$ TeV for $R_{H}=1$ and $R_{\chi}=3$ (this charge assignment for the SM Higgs doublet $H$ and the SM singlet scalar $\chi$ that breaks $U(1)_{R}$ will be used in Secs. 3 and 7).

Furthermore, the coupling of all the fermions with the new gauge boson can be computed from the following relevant part of the Lagrangian:

[TABLE]

The couplings $g_{\psi}$ of all the fermions in our theory are collected in Table II and will be useful for our phenomenological study performed later in the text. Note that the couplings of the left-handed SM fermions are largely suppressed compared to the right-handed ones, since they are always proportional to $\sin\theta_{X}$ and $\theta_{X}$ must be small and is highly constrained by the experimental data.

Based on the framework introduced in this section, we construct various minimal models of Dirac neutrino masses in Sec 3 and study various phenomenology in the subsequent sections.

3 Dirac Neutrino Mass Models

By adopting the set-up as discussed above in this section, we construct models of Dirac neutrino masses. Within this set-up, if the solution $R_{\nu_{i}}=R_{H}$ is chosen which is allowed by the anomaly cancellation conditions, then tree-level Dirac mass term $y_{\nu}v_{H}\overline{\nu}_{L}\nu_{R}$ is allowed and observed oscillation data requires tiny Yukawa couplings of order $y_{\nu}\sim 10^{-11}$. This is expected not to be a natural scenario, hence due to aesthetic reason we generate naturally small Dirac neutrino mass by exploiting the already existing symmetries in the theory. This requires the implementation of the flavor non-universal solution of the anomaly cancellation conditions, in such a scenario $U(1)_{R}$ symmetry plays the vital role in forbidding the direct Dirac mass term and also all Majorana mass terms for the neutrinos.

In this section, we explore three different models within our framework where neutrinos receive naturally small Dirac mass either at the tree-level or at the one-loop level. Furthermore, we also show that the stability of DM can be assured by a residual discrete symmetry resulting from the spontaneous symmetry breaking of $U(1)_{R}$. In the literature, utilizing $U(1)_{R}$ symmetry, two-loop Majorana neutrino mass is constructed with the imposition of an additional $\mathcal{Z}_{2}$ symmetry in [38, 39] and three types of seesaw cases are discussed, standard type-I seesaw in [40], type-II seesaw in [41] and inverse seesaw model in [42]. In constructing the inverse seesaw model, in addition to $U(1)_{R}$, additional flavor dependent U(1) symmetries are also imposed in [42]. In all these models, neutrinos are assumed to be Majorana particles which is not the case in our scenario.

3.1 Tree-level Dirac Seesaw

In this sub-section, we focus on the tree-level neutrino mass generation via Dirac seesaw mechanism [27]222For correlating Dirac seesaw with leptogenesis, see for example [99, 100].. For the realization of this scenario, we introduce three generations of vector-like fermions that are singlets under the SM: $\mathcal{N}_{L,R}(1,1,0,R_{\mathcal{N}})$. In this model, the quantum numbers of the multiplets are shown in Table III and the corresponding Feynman diagram for neutrino mass generation is shown in Fig. 1. This choice of the particle content allows one to write the following Yukawa coupling terms relevant for neutrino mass generation:

[TABLE]

Here, we have suppressed the generation and the group indices. And the Higgs potential is given by:

[TABLE]

When both the $U(1)_{R}$ and EW symmetries are broken, the part of the above Lagrangian responsible for neutrino mass generation can be written as:

[TABLE]

Where $M_{\nu,\mathcal{N}}$ is a $6\times 6$ matrix and, since $\nu_{R_{1}}$ carries a different charge we have $y^{\chi}_{i1}=0$. The bare mass term $M_{\mathcal{N}}$ of the vector-like fermions can in principle be large compared to the two VEVs, $M_{\mathcal{N}}\gg v_{H,\chi}$, assuming this scenario the light neutrino masses are given by:

[TABLE]

Assuming $v_{\chi}=10$ TeV, $y^{H}=y^{\chi}\sim 10^{-3}$, to get $m_{\nu}=0.1$ eV one requires $M_{\mathcal{N}}\sim 10^{10}$ GeV. Dirac neutrino mass generation of this type from a generic point of view without specifying the underline symmetry is discussed in [17].

In this scenario two chiral massless states appear, one of them is $\nu_{R_{1}}$, which is a consequence of its charge being different from the other two generations. In principle, all three generations of neutrinos can be given Dirac mass if the model is extended by a second SM singlet $\chi^{\prime}(1,1,0,-6)$. When this field acquires an induced VEV all neutrinos become massive. This new SM singlet scalar, if introduced, gets an induced VEV from a cubic coupling of the form: $\mu\chi^{2}\chi^{\prime}+h.c.$. Alternatively, without specifying the ultraviolet completion of the model, a small Dirac neutrino mass for the massless chiral states can be generated via the dimension-5 operator $\overline{\mathcal{N}}_{L}\nu_{R}\langle\chi\rangle\langle\chi\rangle/\Lambda$ once $U(1)_{R}$ is broken spontaneously.

3.2 Simplest one-loop implementation

In this sub-section, we consider the most minimal [28] model of radiative Dirac neutrino mass in the context of $U(1)_{R}$ symmetry. Unlike the previous sub-section, we do not introduce any vector-like fermions, hence neutrino mass does not appear at the tree-level. All tree-level Dirac and Majorana neutrino mass terms are automatically forbidden due to $U(1)_{R}$ symmetry reasons. This model consists of two singly charged scalars $S^{+}_{i}$ to complete the loop-diagram and a neutral scalar $\chi$ to break the $U(1)_{R}$ symmetry, the particle content with their quantum numbers is presented in Table IV.

With this particle content, the gauge invariant terms in the Yukawa sector responsible for generating neutrino mass are given by:

[TABLE]

And the complete Higgs potential is given by:

[TABLE]

By making use of the existing cubic term $V\supset\mu S_{2}^{+}S_{1}^{-}\chi+h.c.$ one can draw the desired one-loop Feynman diagram that is presented in Fig. 2. The neutrino mass matrix in this model is given by:

[TABLE]

Here $\theta$ represents the mixing between the singly charged scalars and $m_{H_{i}}$ represents the mass of the physical state $H^{+}_{i}$. Here we make a crude estimation of the neutrino masses: for $\theta=0.1$ radian, $m_{H_{2}}/m_{H_{1}}=1.1$ and $y^{S_{i}}\sim 10^{-3}$ one gets the correct order of neutrino mass $m_{\nu}=0.1$ eV.

This is the most minimal radiative Dirac neutrino mass mechanism which was constructed by employing a $\mathcal{Z}_{2}$ symmetry in [44] and just recently in [33, 28] by utilizing $U(1)_{B-L}$ symmetry. As a result of the anti-symmetric property of the Yukawa couplings $y^{S_{1}}$, one pair of chiral states remains massless to all orders, higher dimensional operators cannot induce mass to all the neutrinos. As already pointed out, neutrino oscillation data is not in conflict with one massless state.

3.3 Scotogenic Dirac neutrino mass

The third possibility of Dirac neutrino mass generation that we discuss in this sub-section contains a DM candidate. The model we present here belongs to the radiative scotogenic [29] class of models and contains a second Higgs doublet in addition to two SM singlets. Furthermore, a vector-like fermion singlet under the SM is required to complete the one-loop diagram. The particle content of this model is listed in Table V and the associated loop-diagram is presented in Fig. 3.

The relevant Yukawa interactions are given as follows:

[TABLE]

And the complete Higgs potential is given by:

[TABLE]

The SM singlet $S$ and the second Higgs doublet $\eta$ do not acquire any VEV and the loop-diagram is completed by making use of the quartic coupling $V\supset\lambda_{D}\eta^{\dagger}H\chi S+h.c.$. Here for simplicity, we assume that the SM Higgs does not mix with the other CP-even states, consequently, the mixing between $S^{0}$ and $\eta^{0}$ originates from the quartic coupling $\lambda_{D}$ (and similarly for the CP-odd states). Then the neutrino mass matrix is given by:

[TABLE]

Where the mixing angle $\theta$ ( $\theta^{\prime}$) between the CP-even (CP-odd) states are given by:

[TABLE]

For a rough estimation we assume no cancellation among different terms occurs. Then by setting $m_{H}=1$ TeV, $M_{\mathcal{N}}=10^{3}$ TeV, $\lambda_{D}=0.1$, $v_{\chi}=10$ TeV, $y^{\eta,S}\sim 10^{-3}$ one can get the correct order of neutrino mass $m_{\nu}\sim 0.1$ eV.

Since $\nu_{R_{1}}$ carries a charge of $-5$, a pair of chiral states associated with this state remains massless. However, in this scotogenic version, unlike the simplest one-loop model presented in the previous sub-section, all the neutrinos can be given mass by extending the model further. Here just for completeness, we discuss a straightforward extension, even though this is not required since one massless neutrino is not in conflict with the experimental data. If the model defined by Table V is extended by two SM singlets $\chi^{\prime}(1,1,0,-6)$ and a $S^{\prime}(1,1,0,\frac{11}{2})$, all the neutrinos will get non-zero mass. The VEV of the field $\chi^{\prime}$ can be induced by the allowed cubic term of the form $\mu\chi^{2}\chi^{\prime}+h.c.$ whereas, $S^{\prime}$ does not get any induced VEV.

Here we comment on the DM candidate present in this model. As aforementioned, we do not introduce new symmetries by hand to stabilize the DM. In search of finding the unbroken symmetry, first, we rescale all the $U(1)_{R}$ charges of the particles in the theory given in Table V including the quark fields in such a way that the magnitude of the minimum charge is unity. From this rescaling, it is obvious that when the $U(1)_{R}$ symmetry is broken spontaneously by the VEV of the $\chi$ field that carries six units of rescaled charge leads to: $U(1)_{R}\to\mathcal{Z}_{6}$. However, since the SM Higgs doublet carries a charge of two units under this surviving $\mathcal{Z}_{6}$ symmetry, its VEV further breaks this symmetry down to: $\mathcal{Z}_{6}\to\mathcal{Z}_{2}$. This unbroken discrete $\mathcal{Z}_{2}$ symmetry can stabilize the DM particle in our theory. Under this residual symmetry, all the SM particles are even, whereas only the scalars $S,\eta$ and vector-like fermions $\mathcal{N}_{L,R}$ are odd and can be the DM candidate. Phenomenology associated with the DM matter in this scotogenic model will be discussed in Sec. 6.

4 Running of the $U(1)_{R}$ Gauge Coupling

In this section, we briefly discuss the running of the $U(1)_{R}$ gauge coupling $g_{R}$, at the one-loop level in our framework. The associated $\beta$-function can be written as:

[TABLE]

Where the coefficient $b_{R}$ can be calculated from [45]:

[TABLE]

The first (second) sum is over the fermions (scalars), $f_{i}$ ($s_{i}$). Here, $\kappa=1/2$ for Weyl fermions, $N_{g}$ is the number of fermion generations, $\eta=2$ for complex scalars and $S_{2}$ are the Dynkin indices of the representations with the appropriate multiplicity factors. By solving Eq. (4.27), the Landau pole can be found straightforwardly:

[TABLE]

The scale of the Landau pole depends on the value of the coupling $g_{R}$, at the input scale $\mu_{0}$. Depending on the choice, both the $\Lambda_{Landau}<M_{Planck}$ and $\Lambda_{Landau}>M_{Planck}$ scenarios can emerge.

Utilizing the basic set-up defined in Sec. 2, we have constructed three different models in Sec. 3, which correspond to three different coefficients $b_{R}=\{179/3,56,731/12\}$ for the Dirac seesaw, simplest one-loop, and Scotogenic models respectively. For demonstration purpose, we choose $\mu_{0}=10$ TeV and show the scale $\Lambda_{Landau}$ as a function of gauge coupling in Fig. 4 for the three different models discussed in this work. As expected, the higher the value of $g_{R}$, smaller the $\Lambda_{Landau}$ gets.

5 Lepton Flavor Violation

In this section, we pay special attention to the charged lepton flavor violation (cLFV) which is an integral feature of these Dirac neutrino mass models. These lepton flavor violating processes provide stringent constraints on TeV-scale extensions of the standard model and, as a consequence put restrictions on the free parameters of our theories. For the first model we discussed, where neutrino masses are generated via Dirac seesaw mechanism, the cLFV decay rates induced by the neutrino mixings (cf. Fig. 5) are highly suppressed by the requirement that the scale of new physics (vector-like fermions $\mathcal{N}_{L,R}$) is at $10^{15}$ GeV to satisfy the neutrino oscillation data, with Yukawa couplings being order one, and hence, are well below the current experimental bounds. Here, we can safely ignore cLFV processes associated with Dirac seesaw model. On the other hand, in the simplest one-loop Dirac neutrino mass model and in the scotogenic model, several new contributions appear due to the additional contributions from charged scalars (cf. Fig. 5), which could lead to sizable cLFV rates.

The cLFV decay processes $\ell_{\alpha}\to\ell_{\beta}+\gamma$ arise from one-loop diagrams are shown in Fig. 5. Let us now focus on the major cLFV processes $\ell_{\alpha}\to\ell_{\beta}+\gamma$ in the simplest one-loop Dirac neutrino mass model. Processes of these types are most dominantly mediated by both the $SU(2)_{L}$ singlet charged scalars ($H_{1,2}^{\pm}$). However, the charged scalar $S_{1}^{\pm}$ determines the chirality of the initial and final-state charged leptons to be left-handed, whereas $S_{2}^{\pm}$ mediated process fixes the chirality to be right-handed and hence there will be no interference between these two contributions. The Yukawa term $y^{S_{1}}$ is anti-symmetric in nature, whereas $y^{S_{2}}$ has completely arbitrary elements in the second and third rows (recall the restriction $y^{S_{2}}_{i1}=0$). We can always make such a judicious choice that no more than one entry in a given row of $y^{S_{2}}$ can be large and thus we can suppress the contribution from the charged scalar $H_{2}^{\pm}$ for the cLFV processes. The expression for $\ell_{\alpha}\to\ell_{\beta}+\gamma$ decay rates can be expressed as333The general expression for this decay rate can be found in Ref. [46, 47]. :

[TABLE]

In Fig. 6, we have shown the contour plots for branching ratio predictions for the cLFV processes: $\mu\to e+\gamma$ (top left), $\tau\to e+\gamma$ (top right) and $\tau\to\mu+\gamma$ (bottom) as a function of mass $(m_{H_{1}})$ and Yukawa $\left|y^{S_{1}}_{i\alpha}y_{i\beta}^{S_{1}*}\right|$ plane in simplest one-loop Dirac neutrino mass model. Red solid lines indicate the current bounds on branching ratios: 4.2 $\times 10^{-13}$ [48] for the $\mu\to e+\gamma$ (top left) process, 3.3 $\times 10^{-8}$ [49] for the $\tau\to e+\gamma$ (top right) process and 4.4 $\times 10^{-8}$ [49] for the $\tau\to\mu+\gamma$ (top right) process. Red dashed lines indicate the future projected bounds on the branching ratios: 6 $\times 10^{-14}$ [50] for the $\mu\to e+\gamma$ (top left), 3 $\times 10^{-9}$ [51] for the $\tau\to e+\gamma$ (top right) and 3 $\times 10^{-9}$ [51] for the $\tau\to\mu+\gamma$ (top right) processes respectively. For simplicity, we choose $m_{H_{2}}=m_{H_{1}}+100$ GeV. As we can see from the Fig. 6, $\mu\to e+\gamma$ is the most constraining cLFV process in this model. Since this could lead to sizable rates, it can be tested in the upcoming experiments.

Similarly, we analyze the major cLFV processes in scotogenic Dirac neutrino mass model. The representative Feynman diagram for the cLFV process $\ell_{\alpha}\to\ell_{\beta}+\gamma$ is shown in Fig. 5 (right diagram). Here also, charged Higgs $H^{\pm}$, which is the part of the $SU(2)_{L}$ doublet $\eta$, mainly contributes to the cLFV process $\ell_{\alpha}\to\ell_{\beta}+\gamma$ (cf. Fig. 5). The decay rate for $\ell_{\alpha}\to\ell_{\beta}+\gamma$ solely depends on the two mass terms $m_{H^{+}},m_{\mathcal{N}}$ and Yukawa term $y^{\eta}$. The decay width expression for this process can be written as:

[TABLE]

Here, $t=m_{F}^{2}/m_{B}^{2}$, and the function $f_{B}(t)$ is expressed as [46, 47]

[TABLE]

In Fig. 7, we have shown the branching ratio predictions for the different cLFV processes: $\mu\to e+\gamma$ (top left), $\tau\to e+\gamma$ (top right) and $\tau\to\mu+\gamma$ (bottom) as a function of mass $(m_{H^{+}})$ in scotogenic one-loop Dirac neutrino mass model for three benchmark values of Yukawas: $\sqrt{\left|y^{\eta}_{\alpha i}y_{\beta i}^{\eta*}\right|}=10^{-1},10^{-2}$ and $10^{-3}$. For our analysis, we set the vector-like fermion mass $m_{\mathcal{N}}$ to be 5 TeV. The $\mu\to e\gamma$ process imposes the most stringent bounds. In this set-up, for the Yukawas: $\sqrt{\left|y^{\eta}_{\alpha i}y_{\beta i}^{\eta*}\right|}=10^{-1},10^{-2}$ and $10^{-3}$, we get charged Higgs mass bounds to be $m_{H^{+}}=$ 3.1 TeV, 4.6 TeV and 5 TeV respectively. As we can see from Fig. 7, most of the parameter space in this model is well-consistent with these cLFV processes and which can be testable at the future experiments. We have shown the future projection reach for these cLFV processes by red dashed lines in Fig. 7.

6 Dark Matter Phenomenology

In this section, we briefly discuss the Dark Matter phenomenology in the scotogenic Dirac neutrino mass model. As aforementioned, in this model, a $\mathcal{Z}_{2}$ subgroup of the original $U(1)_{R}$ symmetry remains unbroken that can stabilize the DM particle. Under this residual symmetry, all the SM particles are even, whereas only the scalars $S,\eta$ and vector-like Dirac fermion $\mathcal{N_{L,R}}$ are odd and the lightest among these can be the DM candidate. DM phenomenology associated with the neutral component of inert scalar doublet, $\eta$ is extensively studied in Ref. [52] in a different set-up and corresponding study has been done for the neutral singlet scalar, $S$ in Ref. [53, 54]. In the following analysis, we consider $\mathcal{N}_{1}$ to be the lightest among all of these particles, hence serves as a good candidate for DM (for simplicity we will drop the subscript from $\mathcal{N}_{1}$ in the following). We aim to study the DM phenomenology associated with the vector-like Dirac fermion $\mathcal{N_{L,R}}$ here. Due to Dirac nature of the dark matter, the phenomenology associated with it is very different from the Majorana fermionic dark matter scenario [55].

In our case, $\mathcal{N}$ pairs can annihilate through s-channel $Z^{\prime}$ exchange process to a pair of SM fermions and right-handed neutrinos. Furthermore, if $m_{DM}>m_{Z^{\prime}}$, then $\mathcal{N}$ may also annihilate directly into pairs of on-shell $Z^{\prime}$ bosons, which subsequently decay to SM fermions. It can also annihilate to SM fermions and right-handed neutrinos via $t-$ channel scalar ($S,\eta_{0},\eta^{+}$) exchanges. The representative Feynman diagrams for the annihilation of DM particle are shown in Fig. 8. It is important to mention that for the Majorana fermionic dark matter case, the annihilation rate is $p-$ wave ($\sim v^{2}$) suppressed since the vector coupling to a self-conjugate particle vanishes, on the contrary, the annihilation rate is not suppressed for the Dirac scenario ($s$-wave). The non-relativistic form for this annihilation cross-section can be found here [58]. In Fig. 9, we analyze the dark matter relic abundance as a function of dark matter mass $m_{DM}$ for various gauge couplings $g_{R}$ (left) and $Z^{\prime}$ boson masses (right). Horizontal red and blue lines represent WMAP [56] relic density constraint $0.094\leq\Omega_{\rm DM}h^{2}\leq 0.128$ and the PLANCK constraint $0.112\leq\Omega_{\rm DM}h^{2}\leq 0.128$ [57] respectively. For simplicity, we set $m_{Z^{\prime}}=10$ TeV (left) and provide the relic abundance prediction for two different values of gauge coupling ($g_{R}$= 0.1 and 0.277). For the right plot in Fig. 9, DM relic abundance is analyzed for two different values of the $Z^{\prime}$ masses $m_{Z^{\prime}}=10$ and 20 TeV setting $g_{R}$= 0.1. As expected, we can satisfy the WMAP [56] relic density constraint $0.094\leq\Omega_{\rm DM}h^{2}\leq 0.128$ and the PLANCK constraint $0.112\leq\Omega_{\rm DM}h^{2}\leq 0.128$ [57] for most of the parameter space in our model as long as $m_{DM}$ is not too far away from $m_{Z^{\prime}}/2$ mass. Throughout our DM analysis, we make sure that we are consistent with the SM $Z-$ boson mass correction constraint while choosing specific $g_{R}$ and $m_{Z^{\prime}}$ values.

In addition to the relic density, we also take into account the constraints from DM direct detection experiments. In case of Majorana fermionic dark matter, at the tree-level, the spin-independent DM-nucleon scattering cross-section vanishes. However, at the loop-level, the spin-independent operators can be generated and hence it is considerably suppressed. The dominant direct detection signal remains the spin-dependent DM-nucleon scattering cross-section which for the Majorana fermionic dark matter is four times that for the Dirac-fermionic dark matter case. In general, the $Z^{\prime}$ interactions induce both spin-independent (SI) and spin-dependent (SD) scattering with nuclei. The representative Feynman diagram for the DM-nucleon scattering is shown in Fig. 10. Particularly, in the scotogenic Dirac neutrino mass model, DM can interact with nucleon through $t-$ channel $Z^{\prime}$ exchange. Hence, large coherent spin-independent scattering may occur since both dark matter and the valence quarks of nucleons possess vector interactions with $Z^{\prime}$ and this process is severely constrained by present direct detection experiment bounds. The DM-nucleon scattering cross-section is estimated in Ref. [58]. In Fig. 11, we analyze the spin-independent dark matter-nucleon scattering cross-section, $\sigma$ (in pb) as a function of the dark matter mass $m_{DM}$ with different gauge coupling $g_{R}=0.2,0.277$. For this plot, we set $m_{Z^{\prime}}=10$ TeV. Yellow, blue and green color solid lines represent current direct detection cross-section limits from LUX-2017 [59], XENON1T [60] and PandaX-II (2017) [61] experiments respectively. As can be seen from Fig. 11, we can satisfy all the present direct detection experiment bounds as long as we are consistent with the other severe bounds on mass $m_{Z^{\prime}}$ and $g_{R}$ arising from colliders to be discussed in the next section.

7 Collider Implications

Models with extra $U(1)_{R}$ implies a new $Z^{\prime}$ neutral boson, which contains a plethora of phenomenological implications at colliders. Here we mainly focus on the phenomenology of the heavy gauge boson $Z^{\prime}$ emerging from $U(1)_{R}$.

7.1 Constraint on Heavy Gauge Boson $Z^{\prime}$ from LEP

There are two kinds of $Z^{\prime}$ searches: indirect and direct. In case of indirect searches, one can look for deviations from the SM which might be associated with the existence of a new gauge boson $Z^{\prime}$. This generally involves precision EW measurements below and above the Z-pole. $e^{+}e^{-}$ collision at LEP experiment [62] above the Z boson mass provides significant constraints on contact interactions involving $e^{+}e^{-}$ and fermion pairs. One can integrate out the new physics and express its influence via higher-dimensional (generally dim-6) operators. For the process $e^{+}e^{-}\to f\bar{f}$, contact interactions can be parameterized by an effective Lagrangian, $\mathcal{L}_{eff}$, which is added to the SM Lagrangian and has the form:

[TABLE]

Where $\Lambda$ is the new physics scale, $\delta_{ef}$ is the Kronecker delta function, $f$ indicates all the fermions in the model and $\eta$ takes care of the chirality structure coefficients. The exchange of the new $Z^{\prime}$ boson state emerging from $U(1)_{R}$ can be stated in a similar way:

[TABLE]

Due to the nature of $U(1)_{R}$ gauge symmetry, the above interaction favors only the right-handed chirality structure. Thus, the constraint on the scale of the contact interaction for the process $e^{+}e^{-}\to l^{+}l^{-}$ from LEP measurements [62] will indirectly impose bound on $Z^{\prime}$ mass and the gauge coupling $(g_{R})$ that can be translated into:

[TABLE]

Other processes such as $e^{+}e^{-}\to c\bar{c}$ and $e^{+}e^{-}\to b\bar{b}$ impose somewhat weaker bounds than the ones quoted in Eq. 7.35.

7.2 Heavy Gauge Boson $Z^{\prime}$ at the LHC

Now we analyze the physics of the heavy neutral gauge boson $Z^{\prime}$ at the Large Hadron Collider (LHC). At the LHC, $Z^{\prime}$ can be resonantly produced via the quark fusion process $q\bar{q}\to Z^{\prime}$ since the coupling of $Z^{\prime}$ with right-handed quarks ($u_{R},d_{R}$) are not suppressed. After resonantly produced at the LHC, $Z^{\prime}$ will decay into SM fermions and also to the exotic scalars ($S_{2}^{+}S_{2}^{-},\chi\chi$) or fermions ($\mathcal{N}\mathcal{N}$) depending on the model if kinematically allowed444Even if we include $Z^{\prime}\to\mathcal{N}\mathcal{N},\rm~{}S_{2}^{+}S_{2}^{-},\rm~{}\chi\chi$ decay modes, the branching fraction ($\sim 4\%$) for $Z^{\prime}\to e^{+}e^{-}/{\mu}^{+}{\mu}^{-}$ mode does not change much.. The present lack of any signal for di-lepton resonances at the LHC dictates the stringent bound on the $Z^{\prime}$ mass and $U(1)_{R}$ coupling constant $g_{R}$ in our model as the production cross-section solely depends on these two free parameters. Throughout our analysis, we consider that the mixing $Z-Z^{\prime}$ angle is not very sensitive ($s_{X}=0$). In order to obtain the constraints on these parameter space, we use the dedicated search for new resonant high-mass phenomena in di-electron and di-muon final states using $36.1$ fb*-1* of proton-proton collision data, collected at $\sqrt{s}=13$ TeV by the ATLAS collaboration [63]. The searches for high mass phenomena in di-jet final states [64] will also impose bound on the model parameter space, but it is somewhat weaker than the di-lepton searches due to large QCD background. For our analysis, we implement our models in FeynRules$\char 95\relax$v2.0 package [65] and simulate the events for the process $pp\to Z^{\prime}\to e^{+}e^{-}({\mu}^{+}{\mu}^{-})$ with MadGraph5$\char 95\relax$aMC@NLO$\char 95\relax$v3$\char 95\relax$0$\char 95\relax$1 code [66]. Then, using parton distribution function (PDF) NNPDF23$\char 95\relax$lo$\char 95\relax$as$\char 95\relax$0130 [67], the cross-section and cut efficiencies are estimated. Since no significant deviation from the SM prediction is observed in experimental searches [63] for high-mass phenomena in di-lepton final states, the upper limit on the cross-section is derived from the experimental analyses [63] using $\sigma\times$ BR $=N_{rec}/(A\times\epsilon\times\int Ldt)$, where $N_{rec}$ is the number of reconstructed heavy $Z^{\prime}$ candidate, $\sigma$ is the resonant production cross-section of the heavy $Z^{\prime}$, BR is the branching ratio of $Z^{\prime}$ decaying into di-lepton final states , $A\times\epsilon$ is the acceptance times efficiency of the cuts for the analysis. In Fig. 12, we have shown the upper limits on the cross-section at 95$\%$ C.L. for the process $pp\to Z^{\prime}\to l^{+}l^{-}$ as a function of the di-lepton invariant mass using ATLAS results [63] at $\sqrt{s}$ = 13 TeV with 36.1 fb*-1* integrated luminosity. Red solid, dashed and dotted lines in Fig. 12 indicate the model predicted cross-section for three different values of $U(1)_{R}$ gauge coupling constant $g_{R}=0.5,~{}0.3,~{}0.1$ respectively. We find that $Z^{\prime}$ mass should be heavier555For related works see also [68, 69]. than $4.4,3.9$ and $2.9$ TeV for three different values of $U(1)_{R}$ gauge coupling constant $g_{R}=0.5,~{}0.3$ and $0.1$.

In Fig. 13, we have shown all the current experimental bounds in $M_{Z^{\prime}}-g_{R}$ plane. Red meshed zone is excluded from the current experimental di-lepton searches [63]. The cyan meshed zone is forbidden from the LEP constraint[62] and the blue meshed zone is excluded from the limit on SM Z boson mass correction: $\frac{1}{3}M_{Z^{\prime}}/g_{R}>12.082$ TeV as aforementioned. We can see from Fig. 13 that the most stringent bound in $M_{Z^{\prime}}-g_{R}$ plane is coming from direct $Z^{\prime}$ searches at the LHC. After imposing all the current experimental bounds, we analyze the future discovery prospect of this heavy gauge boson $Z^{\prime}$ within the allowed parameter space in $M_{Z^{\prime}}-g_{R}$ plane looking at the prompt di-lepton resonance signature at the LHC. We find that a wider region of parameter space in $M_{Z^{\prime}}-g_{R}$ plane can be tested at the future collider experiment. Black, green, purple and brown dashed lines represent the projected discovery reach at $5\sigma$ significance at 13 TeV LHC for 100 fb*-1*, 300 fb*-1*, 500 fb*-1* and 1 ab*-1* luminosities. On the top of that, the right-handed chirality structure of $U(1)_{R}$ can be investigated at the LHC by measuring Forward-Backward (FB) and top polarization asymmetries in $Z^{\prime}\to t\bar{t}$ mode [70] and which can discriminate our $U(1)_{R}$ $Z^{\prime}$ interaction from the other $Z^{\prime}$ interactions in $U(1)_{B-L}$ model. The investigation of other exotic decay modes $(\mathcal{N}\mathcal{N},\rm~{}\chi\chi,\rm~{}S_{2}^{+}S_{2}^{-})$ of heavy $Z^{\prime}$ is beyond the scope of this article and shall be presented in a future work since these will lead to remarkable multi-lepton or displaced vertex signature [71, 72, 73, 74, 75, 76, 77] at the colliders.

7.3 Heavy Gauge Boson $Z^{\prime}$ at the ILC

Due to the point-like structure of leptons and polarized initial and final state fermions, lepton colliders like ILC will provide much better precision of measurements. The purpose of the $Z^{\prime}$ search at the ILC would be either to help identifying any $Z^{\prime}$ discovered at the LHC or to extend the $Z^{\prime}$ discovery reach (in an indirect fashion) following effective interaction. Even if the mass of the heavy gauge boson $Z^{\prime}$ is too heavy to directly probe at the LHC, we will show that by measuring the process $e^{+}e^{-}\to f^{+}f^{-}$, the effective interaction dictated by Eq. 7.34 can be tested at the ILC. Furthermore, analysis with the polarized initial states at ILC can shed light on the chirality structure of the effective interaction and thus it can distinguish between the heavy gauge boson $Z^{\prime}$ emerging from $U(1)_{R}$ extended model and the $Z^{\prime}$ from other $U(1)$ extended model such as $U(1)_{B-L}$. The process $e^{+}e^{-}\to f^{+}f^{-}$ typically exhibits asymmetries in the distributions of the final-state particles isolated by the angular- or polarization-dependence of the differential cross-section. These asymmetries can thus be utilized as a sensitive measurement of differences in interaction strength and to distinguish a small asymmetric signal at the lepton colliders. In the following, the asymmetries (Forward-Backward asymmetry, Left-Right asymmetry) related to this work will be described in great detail.

7.3.1 Forward-Backward Asymmetry

The differential cross-section in Eq. 7.3.1 is asymmetric in polar angle, leading to a difference of cross-sections for $Z^{\prime}$ decays between the forward and backward hemispheres. Earlier, LEP experiment [62] used Forward-backward asymmetries to measure the difference in the interaction strength of the $Z$-boson between left-handed and right-handed fermions, which gives a precision measurement of the weak mixing angle. Here we will show that our framework leads to sizable and distinctive Forward-Backward (FB) asymmetry discriminating from other models and which can be tested at the ILC, since only the right-handed fermions carry non-zero charges under the $U(1)_{R}$. For earlier analysis of FB asymmetry in the context of other models as well as model-independent analysis see for example Refs. [78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 40, 42].

At the ILC, $Z^{\prime}$ effects have been studied for the following processes:

[TABLE]

where $\sigma_{i}=\pm 1$ are the helicities of initial (final)-state leptons and $k_{i}$’s are the momenta. Since the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process is the most sensitive one at the ILC, we will focus on this process only for the rest of our analysis. One can write down the corresponding helicity amplitudes as:

[TABLE]

where $s=(k_{1}+k_{2})^{2}=(k_{3}+k_{4})^{2}$, $s_{Z}=s-m_{Z}^{2}+im_{Z}\Gamma_{Z}$, and $\cos\theta$ indicates the scattering polar angle. $e^{2}=4\pi\alpha$ with $\alpha=$ QED coupling constant, $c_{R}=\tan\theta_{W}$ and $c_{L}=-\cot 2\theta_{W}$ and $\theta_{W}$ is the weak mixing angle.

For a purely polarized initial state, the differential cross-section is expressed as:

[TABLE]

Then the differential cross-section for the partially polarized initial state with a degree of polarization $P_{e^{-}}$ for the electron beam and $P_{e^{+}}$ for the positron beam can be written as [78, 40]:

[TABLE]

One can now define polarized cross-section $\sigma_{L,R}$ (for the realistic values at the ILC [89]) as:

[TABLE]

Using this one can study the initial state polarization-dependent forward-backward asymmetry as:

[TABLE]

where

[TABLE]

where $\mathcal{L}$ represents the integrated luminosity, $\epsilon$ indicates the efficiency of observing the events, and $c_{\rm max}$ is a kinematical cut chosen to maximize the sensitivity. For our analysis we consider $\epsilon=1$, and $c_{\rm max}=0.95$. Then we estimate the sensitivity to $Z^{\prime}$ contribution by:

[TABLE]

where $A_{FB}^{SM+Z^{\prime}}$ and $A_{FB}^{SM}$ are FB asymmetry originated from both the SM and $Z^{\prime}$ contribution and from the SM case only. Next, it is compared with the statistical error of the asymmetry (in only SM case) $\delta A_{FB}$ [78, 40]:

[TABLE]

In Fig. 14, we analyze the strength of FB asymmetry $\Delta A_{FB}$ as a function of VEV $v_{\chi}(=M_{Z^{\prime}}/{3g_{R}})$ for both left and right-handed polarized cross-sections of the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process. In order to compare, we have done the analysis for both the cases: $Z^{\prime}$ from both $U(1)_{R}$ and $U(1)_{B-L}$ cases. We have considered the center of mass energy for the ILC at $\sqrt{s}=500$ GeV and the integrated luminosity $\mathcal{L}$ is set to be 1 ab*-1* (5 ab*-1*) for the left (right) panel of Fig. 14. The grey shaded region corresponds to excluded region from the SM $Z$ boson mass correction. Red dashed (solid) line represents $\Delta A_{FB}$ for $U(1)_{R}$ case for left (right) handed polarized cross-sections of the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process, whereas blue dotted (dashed) line indicates $\Delta A_{FB}$ for $U(1)_{B-L}$ case for left (right) handed polarized cross-sections. From Fig. 14, we find that in case of $U(1)_{R}$ model, it provides significant difference of $\Delta A_{FB}$ for $\sigma_{R}$ and $\sigma_{L}$ due to the right-handed chirality structure of $Z^{\prime}$ interaction from $U(1)_{R}$, while in the case of $U(1)_{B-L}$ model, it provides small difference. Hence by comparing the difference of $\Delta A_{FB}$ for differently polarized cross-section $\sigma_{R}$ and $\sigma_{L}$ at the ILC, we can easily discriminate the $Z^{\prime}$ interaction from $U(1)_{R}$ and $U(1)_{B-L}$ model. As we can see from Fig. 14 that there are significant region for $M_{Z^{\prime}}/3g_{R}>12.082$ TeV which can give more than $2\sigma$ sensitivity for FB asymmetry by looking at $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process at the ILC. We can also expect much higher sensitivity while combining different final fermionic states such as other leptonic modes ($e^{+}e^{-},{\tau}^{+}{\tau}^{−}$) as well as hadronic modes $jj$. Moreover, the sensitivity to $Z^{\prime}$ interactions can be enhanced by analyzing the scattering angular distribution in details, although it is beyond the scope of our paper.

7.3.2 Left-Right Asymmetry

The simplest example of the EW asymmetry for an experiment with a polarized electron beam is the left-right asymmetry $A_{LR}$, which measures the asymmetry at the initial vertex. Since there is no dependence on the final state fermion couplings, one can get an advantage by looking at LR asymmetry at lepton collider. Another advantage of this LR asymmetry measurement is that it is barely sensitive to the details of the detector. As long as at each value of $\cos{\theta}$, its detection efficiency of fermions is the same as that for anti-fermions, the efficiency effects should be canceled within the ratio because the $Z^{\prime}$ decays into a back-to-back fermion-antifermion pair and about the midplane perpendicular to the beam axis, the detector was designed to be symmetric. For earlier studies on LR asymmetry in different contexts, one can see for example Refs. [78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 90]. LR asymmetry is defined as:

[TABLE]

where $N_{L}$ is the number of events in which initial-state particle is left-polarized, while $N_{R}$ is the corresponding number of right-polarized events.

[TABLE]

Similarly, one can estimate the sensitivity to $Z^{\prime}$ contribution in LR asymmetry by [79, 90, 82]:

[TABLE]

with a statistical error of the asymmetry $\delta A_{{LR}}$ , given [79, 90, 82] as

[TABLE]

In Fig. 15, we analyze the strength of LR asymmetry $\Delta A_{LR}$ for the $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process as a function of VEV $v_{\chi}(=M_{Z^{\prime}}/{3g_{R}})$. In order to distinguish $Z^{\prime}$ interaction, we have analysed both the cases: $Z^{\prime}$ emerging from both $U(1)_{R}$ and $U(1)_{B-L}$ cases. We have considered the center of mass energy for the ILC at $\sqrt{s}=500$ GeV and the integrated luminosity $\mathcal{L}$ is set to be 1 ab*-1* (5 ab*-1*) for the left (right) panel of Fig. 15. The grey shaded region corresponds to excluded region from the SM $Z$ boson mass correction. Red (blue) solid line represents $\Delta A_{LR}$ for $U(1)_{R}$ ($U(1)_{B-L}$) case. From Fig. 15, we find that in case of $U(1)_{R}$ model, it provides remarkably large LR asymmetry $\Delta A_{LR}$ due to the right-handed chirality structure of $Z^{\prime}$ interaction from $U(1)_{R}$, while in case of $U(1)_{B-L}$ model, it gives a smaller contribution. Hence by comparing the difference of $\Delta A_{LR}$ at the ILC, we can easily discriminate the $Z^{\prime}$ interaction from $U(1)_{R}$ and $U(1)_{B-L}$ model. As we can see from Fig. 15 that there is a significant region for $M_{Z^{\prime}}/3g_{R}>12.082$ TeV which can give more than $3\sigma$ sensitivity for LR asymmetry by looking at $e^{+}e^{-}\to\mu^{+}\mu^{-}$ process at the ILC. Even if, we can achieve $5\sigma$ sensitivity for a larger parameter space in our framework if integrated luminosity of ILC is upgraded to $5$ ab*-1*. Although, measurement of both the FB and LR asymmetries at the ILC can discriminate $Z^{\prime}$ interaction for $U(1)_{R}$ model from other $U(1)$ extended models such as $U(1)_{B-L}$ model, it is needless to mention that the LR asymmetry provides much better sensitivity than the FB asymmetry in our case. In Fig. 16, we have shown the survived parameter space in $M_{Z^{\prime}}-g_{R}$ plane satisfying all existing bounds and which can be probed at the ILC in future by looking at LR asymmetry strength. Green and yellow shaded zones correspond to sensitivity confidence levels 1$\sigma$ and 2$\sigma$ by measuring LR asymmetry for $U(1)_{R}$ extended model at the ILC. For higher $Z^{\prime}$ mass (above $\sim$ 10 TeV), it is too heavy to directly produce and probe at the LHC looking at prompt di-lepton signature. On the other hand, ILC can probe the heavy $Z^{\prime}$ effective interaction and LR asymmetry can pin down/distinguish our $U(1)_{R}$ model from other existing $U(1)$ extended model for a large region of the parameter space. Thus, $Z^{\prime}$ search at the ILC would help to identify the origin of $Z^{\prime}$ boson as well as to extend the $Z^{\prime}$ discovery reach following effective interaction.

8 Constraint from Cosmology

In the previous section, we have extensively analyzed the collider implications of the new gauge boson $Z^{\prime}$. In this section, we aim to study the constraints on the mass of the new gauge boson from cosmological measurements and compare with the collider bounds. Since the right-handed neutrinos carry non-zero $U(1)_{R}$ charge in our set-up, they couple to the SM sector via the $Z^{\prime}$ boson interactions. Furthermore, since they are either massless or very light, they contribute to the relativistic degrees of freedom $N_{eff}$, hence in principle can increase the expansion rate of the Universe. Their contribution to this process is parametrized by $\Delta N_{eff}$ and to compute it we follow the procedure discussed in Ref. [91]. After $\nu_{R}$ states decouple, specifically for $T<T_{dec}^{\nu_{L}}<T^{\nu_{R}}_{dec}$ ($T_{dec}^{\nu_{L/R}}$ represents the decoupling temperature of the $\nu_{L/R}$ neutrinos) their total contribution is given by:

[TABLE]

here $N_{\nu_{R}}$ is the number of massless or light right-handed neutrinos, $g(T)$ is the relativistic degrees of freedom at temperature T, with the well-known quantities $g(T_{dec}^{\nu_{L}})=43/4$ and $T_{dec}^{\nu_{L}}=2.3$ MeV [92]. For the following computation, we take the temperature-dependent degrees of freedom from the data listed in Table S2 of Ref. [93], and by utilizing the cubic spline interpolation method, we present $g$ as a function of $T$ in Fig. 17 (left plot).

The current cosmological measurement of this quantity is $N_{eff}=2.99^{+0.34}_{-0.33}$ [94], which is completely consistent with the SM prediction $N^{SM}_{eff}=3.045$ [95]. These data limit the contribution of the right-handed neutrinos to be $\Delta N_{eff}<0.285$. However, future measurements [96] can put even tighter constraints on this deviation $\Delta N_{eff}<0.06$. The right-handed neutrinos decouple from the thermal bath when the interaction rate drops below the expansion rate of the Universe:

[TABLE]

Here the Hubble expansion parameter is defined as:

[TABLE]

where $M_{Pl}$ is the Planck mass and $g_{\nu_{R}}=2$ is the spin degrees of freedom of the right-handed neutrinos. And the interaction rate that keeps the right-handed neutrinos at the thermal bath is given by:

[TABLE]

Here, the Fermi-Dirac distribution is $f_{\nu_{R}}(p)=1/(e^{p/T}+1)$, the number density is $n_{\nu_{R}}=\left(3/(2\pi^{2})\right)\zeta(3)T^{3}$, $s=2pq(1-\cos\theta)$ and $v=1-\cos\theta$. Furthermore, the annihilation cross-section $\sigma(\overline{\nu_{R}}\nu_{R}\to\overline{f}_{i}f_{i})$ is as follows:

[TABLE]

Where $N^{f}_{C}$ and $Q_{f}$ represent the color degrees of freedom and the charge under the $U(1)_{R}$ for a fermion $f$ respectively.

By plugging Eqs. (8.57)-(8.59) in Eq. (8.56) and then solving numerically, we present our result of $\Delta N_{eff}$ as a function of $M_{Z^{\prime}}/g_{R}$ in Fig. 17 (right plot). From this figure, one sees that cosmology provides strong bound on the mass of the new gauge boson based on the associated decoupling temperature of the right-handed neutrinos. The blue curve corresponds to the contribution of all the three right-handed neutrinos and the red dashed line represents the current experimental upper bound on the deviation of $\Delta N_{eff}$. This bound puts the restriction $M_{Z^{\prime}}/g_{R}\gtrsim 26.5$ TeV, which is quite stronger than the LEP bound $M_{Z^{\prime}}/g_{R}\gtrsim 3.59$ TeV, however, lies within the constraint provided by the SM $Z$-boson mass correction $M_{Z^{\prime}}/g_{R}\gtrsim 36.2$ TeV. The framework presented in this work puts larger bound on the mass of the new gauge boson from cosmology due to large charge assignment of the right-handed neutrinos compared to the conventional $U(1)_{B-L}$ models with universal charge, $M_{Z^{\prime}}/g_{B-L}\gtrsim 14$ TeV [97, 98].

9 Conclusions

We believe that the scale of new physics is not far from the EW scale and a simple extension of the SM should be able to address a few of the unsolved problems of the SM. Adopting this belief, in this work, we have explored the possibility of one of the most minimal gauge extensions of the SM which is $U(1)_{R}$ that is responsible for generating Dirac neutrino mass and may also stabilize the DM particle. Cancellations of the gauge anomalies are guaranteed by the presence of the right-handed neutrinos that pair up with the left-handed partners to form Dirac neutrinos. Furthermore, this $U(1)_{R}$ symmetry is sufficient to forbid all the unwanted terms for constructing naturally light Dirac neutrino mass models without imposing any additional symmetries by hand. The chiral non-universal structure of our framework induces asymmetries, such as forward-backward asymmetry and especially left-right asymmetry that are very distinct compared to any other $U(1)$ models. By performing detailed phenomenological studies of the associated gauge boson, we have derived the constraints on the $U(1)_{R}$ model parameter space and analyzed the prospect of its testability at the collider such as at LHC and ILC. We have shown that a heavy $Z^{\prime}$ (emerging from $U(1)_{R}$), even if its mass is substantially higher than the center of mass energy available at the ILC, would manifest itself at tree-level by its propagator effects producing sizable contributions to the LR asymmetry or FB asymmetry. This can be taken as an initial guide to explore the $U(1)_{R}$ model at colliders. These models can lead to large lepton flavor violating observables which we have studied and they could give a complementary test for these models. In this work, we have also analyzed the possibility of having a viable Dirac fermionic DM candidate stabilized by the residual discrete symmetry originating from $U(1)_{R}$, which connects to SM via $Z^{\prime}$ portal coupling in a framework that also cater for neutrino mass generation. The DM phenomenology is shown to be crucially dictated by the interaction of $\mathcal{N}$ with $Z^{\prime}$. Furthermore, we have inspected the constraints coming from the cosmological measurements and compared this result with the different collider bounds. For a comparison, here we provide a benchmark point by fixing the gauge coupling $g_{R}=0.056$. With this, the current upper bound on the $Z^{\prime}$ mass is $M_{Z^{\prime}}>4.25$ TeV from 13 TeV LHC data with $36.1fb^{-1}$ luminosity, and the future projection reach limit translates into $M_{Z^{\prime}}>4.67$ TeV with $100fb^{-1}$ luminosity. Whereas for the same value of the gauge coupling, the ILC has the discovery reach of $4.63$ TeV at the $2\sigma$ confidence level looking at the left-right asymmetry. The corresponding bounds from LEP, $Z-$boson mass correction and from cosmology are $M_{Z^{\prime}}>0.2,2,1.49$ TeV respectively, which are somewhat weaker compared to LHC and ILC bounds. To summarize, the presented Dirac neutrino mass models are well motivated and have rich phenomenology.

Acknowledgments

We thank K. S. Babu, Bhupal Dev and S. Nandi for useful discussions. The work of SJ and VPK was in part supported by US Department of Energy Grant Number DE-SC 0016013. The work of SJ was also supported in part by the Neutrino Theory Network Program. SJ thanks the Theoretical Physics Department at Washington University in St. Louis for warm hospitality during the completion of this work.