Anomaly-free Abelian gauge symmetries with Dirac seesaws

We perform a systematic analysis of Standard Model extensions with an additional anomaly-free gauge $U(1)$ symmetry, to generate tree-level Dirac neutrino masses. An anomaly-free symmetry demands nontrivial conditions on the charges of the unavoidable new states. An intensive scan was performed, looking for solutions generating neutrino masses by the type-I and type-II tree-level Dirac seesaw mechanism, via operators with dimension 5 and 6, that correspond to active or dark symmetries. Special attention was paid to the cases featuring no extra massless chiral fermions or multicomponent dark matter with unconditional stability.


Introduction
Despite its enormous success, the standard model (SM) has to be extended in order to account for neutrino masses and dark matter (DM). Regarding neutrinos, the oscillation data is compatible with both Majorana or Dirac neutrino masses [1], with no strong preference for either of the two possibilities. Even if most of the literature assumes that neutrinos are Majorana in nature (see, e.g., Ref. [2] for a review), the mass generation mechanism for Dirac neutrinos has recently received increased attention.
To explain Dirac neutrino masses, right-handed neutrinos (RHNs) have to be introduced. Additionally, an extra local symmetry is also required to guarantee proper total lepton number conservation [3]. Even so, the required Yukawa couplings are typically very suppressed, of the order O (10 −10 ), if Dirac neutrino masses are induced directly from the SM Higgs mechanism [4,5]. Nevertheless, if the symmetry forbids the tree-level contribution driven by the SM Higgs, a Dirac-seesaw mechanism can be implemented [6,7]. For example, the type-I Dirac-seesaw [8][9][10] could appear in the context of anomaly-free gauge U (1) B−L [3] or U (1) R [11] symmetries. Specific models have been studied for producing neutrino masses via the type-II Dirac Seesaw [12,13] with both anomaly-free gauge U (1) B−L symmetries [14,15] or with anomaly-free Abelian dark symmetries [16]. The extra chiral singlet fermions required to cancel out the anomalies can be part of a hidden sector with DM candidates, as in the case type-I [17,18] or type-II [19] Majorana seesaw. At loop level, the heavy particles in the radiative seesaw can be fully associated to an Abelian gauge dark symmetry U (1) D with the lightest of them as DM candidate [20][21][22][23], as well to an active Abelian gauge symmetry U (1) X , like an U (1) B−L [4,[23][24][25]. The studies of one-loop Dirac neutrino masses have typically focused on finding specific anomaly-free solutions of this two kinds of symmetries, see, e.g., Refs. [11,[26][27][28][29]. Only few studies have performed systematic analysis of SM extensions with an additional anomaly-free gauge U (1) symmetry to generate Dirac neutrino masses [30,31].
Here we continue the effort of performing a systematic analysis of SM extensions, presenting a complete set of relevant anomaly-free solutions to the general problem of the generation of neutrino masses by the type-I and type-II tree-level Dirac seesaw mechanism, via dimension 5 and 6 operators. A full set of relevant solutions is obtained. Our method can be easily applied to find the full set of anomaly-free solutions to well defined phenomenological problems.
In this work, we look for anomaly-free solutions to SM extensions with an additional U (1) gauge symmetry, giving rise to tree-level Dirac neutrino masses. For that purpose, in section 2 we briefly revise the conditions to have a non-anomalous U (1) gauge symmetry. In section 3 we discuss solutions giving rise to neutrino masses, by the type-I and type-II tree-level Dirac seesaw mechanism, via operators with dimensions 5 and 6, corresponding to active or dark symmetries. Finally, in section 4 our conclusions are presented.

Anomaly conditions
We consider an extension of the SM with an additional U (1) X gauge symmetry, and N right-handed chiral fields ψ ρ singlets under the SM SU (3) c ⊗SU (2) L ⊗U (1) Y group, with charges n ρ under the U (1) X , where ρ = 1, · · · , N . Additionally, we assume that the SM right-handed chiral fermions transform under the active U (1) X symmetry, with charges denoted with the same name of the field. 1 To avoid having an anomalous U (1) X , the three linear anomaly conditions to be fulfilled are As they only depend on the SM fermions, three of their X-charges can be expressed in terms of the other two [23,[32][33][34], chosen to be e and L, as We note that the quadratic anomaly condition in U (1) X is trivially satisfied. However, the mixed gaugegravitational [Grav] 2 U (1) X and the cubic [U (1) X ] 3 anomalies do depend on the extra fermion charges n ρ , and therefore two additional conditions have to be imposed in order to avoid an anomalous U (1) X [23]: where m ≡ e + 2L. Equation (4) can be conveniently rewritten as Finally, we note that the SM Higgs must have an X-charge to guarantee that SM quarks and charged leptons acquire masses through the standard Higgs mechanism. 2 Along these lines, we also assume that the singlet chiral fermions ψ ρ only acquire mass through the spontaneous symmetry breaking (SSB) of the extra U (1) X symmetry. This excludes solutions with vector-like states. We note that the existence of fields charged under both hypercharge and U (1) X induce at loop level the kinetic mixing operator L ⊃ 2 B µν X µν , where B µν and X µν are the field strengths related to the U (1) Y and the extra U (1) X , respectively. The dimensionless parameter depends on the masses of the particles in the loop, as well as their specific charge assignment and the gauge couplings under the two U (1) symmetries [35][36][37].
It is interesting to note that the conditions in Eq. (5) are completely equivalent to the ones coming from a scenario where the SM is extended with a dark U (1) D gauge symmetry with N = N + 3 righthanded singlet chiral fermions, N of them with the charges n ρ and three with charge m, and where the SM is invariant (hence a dark symmetry) [31]. Even if comparable, there is a major technical advantage of the latter approach: If the SM is extended with and additional dark U (1) D gauge symmetry (under which it is uncharged), and N right-handed chiral fields singlets under the SM group, the U (1) D is not anomalous if the Diophantine equations

Dirac seesaw models
In this section we look for anomaly-free U (1) X gauge extensions of the SM, with N singlet chiral fermions, realizing the effective Dirac neutrino mass operators [12,38] at tree-level. In the two-component spinor Figure 1: Diagram leading to tree-level type-I Dirac neutrino masses, via the dimension 6 operator described in Eq. (9). notation, they can be written as and δ = 1 or 2 for dimension 5 (D-5) or 6 (D-6) operators, respectively. Here h αi ν correspond to dimensionless induced couplings, ν Rα are at least two RHNs (α = 1, 2, . . .) with the same X-charge ν, L i are the lepton doublets with X-charge −L, H is the SM Higgs doublet with X-charge h = L − m, S is the complex singlet scalar responsible for the SSB of the anomaly-free gauge symmetry with X-charge s = −(ν + m)/δ, respectively, and Λ is a scale of new physics, which is parametrically the typical mass scale of the new (heavy) states. In general, after the SSB, a remnant Z |s| discrete symmetry is left, which guarantees the stability of a potential DM candidate [21].

Type-I Dirac Seesaw
The realizations of type-I Dirac seesaw with D-5 operators automatically imply a vector-like pair of singlet fermions, and is therefore not possible with only chiral fermions. Hence, we will not consider that case any longer in this study. However, the situation is different for D-6 operators. Figure 1 presents the diagram realizing the D-6 effective Lagrangian in Eq. (9) for type-I seesaw, if one only allows SSB masses for the singlet chiral fermions. It is important to note that this process can only take place for active symmetries (a dark symmetry would imply a state χ R with charge r = 0). In that figure, the flux of charges in each vertex imply that The charge of S can be expressed as a function of l and m as and the chiral fermion charges obey It is important to note that the first condition in Eq. (10) does not imply vector-like fermions, even if it contains particles with opposite charges, since the charge m is associated to SM doublets. In this way, for any chiral solution (i.e., without opposite charges) we have two minimal ways of building a solution satisfying the condition in Eq. (12): • If m is already present in the solution as a non-repeated X-charge (m, . . .), we can add two sets of opposite sign charges m, such that where −m is the X-charge of two new chiral fields χ R1 , χ R2 . This is the minimal requirement for a rank-2 neutrino mass matrix for light neutrinos. We note that if the original solution satisfies the Diophantine equations, the second will trivially do it as well.
• If m is not present in the solution we can add three sets of opposite sign charges m, and identify the X-charges of three new chiral fields χ R1 , χ R2 and χ R3 as −m.
As mentioned previously, in order to limit the total number of solutions that cancel the anomaly induced by the additional U (1) X , the following restrictions are taken into account: 1. By construction, all new chiral fermions have to be charged under U (1) X , i.e., solutions with vanishing charges are disregarded.
2. For the chiral fields, the maximal charge allowed (in absolute value) is 30.
3. Solutions with vector-like fermions are disregarded. We emphasize that the first condition in Eq. (10) together with the requirement of a rank-2 neutrino mass matrix demand at least 2 opposite pairs of charges. However, this does not imply vector-like fermions, since the charges m are associated to SM doublets.
4. At least two charges have to be equal. Their corresponding fields are identified with the RHNs.

5.
A second set of at least two equal charges have to exist. Their corresponding fields are identified with the heavy left-handed chiral fermions. This is due to the need of a neutrino mass matrix for the light neutrinos of at least rank 2.
6. A third set of three equal charges is required. Their corresponding fields are identified with the SM doublets charged under the new U (1) X .
7. We restrict ourselves to N ≤ 9 fields, with charges satisfying the two Diophantine conditions in Eq. (8), and take the minimal charge (in absolute value) to be positive. We note that there are no solutions for N ≤ 5 with at least two sets of equal charges [39,40].
8. The numbers of left-handed chiral fermions χ R and χ L has to be the same, for their masses to be generated by the SSB of the S.
9. The charge assignment may not allow all chiral fields to acquire masses via the SSB. We only consider solutions which have extra massless chiral fields all of them acquiring masses through an extra singlet scalar S with X-charge s . This implies that these charges cannot be repeated. 10. We want RHN masses to be generated by tree-level Dirac seesaw. That implies that all vertices between S, the RHN and the other chiral fields should be forbidden by the symmetries.
11. In the case of new set of chiral fermions which get masses through S the lightest fermion must be a viable DM candidate .
12. In order to have a viable type-I Dirac seesaw, new heavy chiral fermions have to satisfy the condition in Eq. (12). For type-II Dirac seesaw, the conditions to be imposed appear in Eqs. (15) and (17) for active or dark symmetries, respectively, as will be seen in the following.
Additionally, we note that for a fixed number of chiral fields, different solutions could share the same qualitative behavior. For example, the second solution in Table 1 is (1, Table 1. The solutions of the Diophantine equations satisfying all the previously enumerated conditions are shown in Table 1. The solutions for N extra chiral fermions are parametrized as a function of two sets of integers and k (first three columns). The fourth column shows the charge assignments, whereas the fifth the general common denominator (GCD) of the original solution. Technically speaking, the solutions were found using the package anomalies 3 [31]. We note that even if most of the solutions contain massless chiral fermion, there are two solutions without (N highlighted in bold). Regarding these solutions without massless fermions, a few comments are in order: • The first solution corresponds to N = 9 and has the charge assignment (6, 6, 6, −6, −6, 2, −5,   All other solutions presented in Table 1 have a number of massless chiral fermions. They can be either extra relativistic degrees of freedom, or additional DM candidates if they acquire mass from another mechanism. Concerning the solutions with multicomponent DM, we explore the cases which feature at least two DM candidates with unconditional stability [41]. This happens when there are two remnant symmetries such that Z |s| ∼ = Z p ⊗ Z q with Z p ⊗ Z q coprimes, which guaranteed the stability of each lightest state under Z p and Z q respectively, without imposing any kinematical restriction. For the two DM candidates associated to the set of chiral fields ψ i and χ j , we consider the first two possibilities for |s| [41] • Z 6 ∼ = Z 2 ⊗ Z 3 : solutions with at least a set of chiral fields with ψ i ∼ [ω 2 6 ∨ ω 4 6 ] under Z 6 , and at least a set of chiral fields with χ i ∼ ω 3 6 under Z 6 , • Z 10 ∼ = Z 2 ⊗ Z 5 : solutions with at least a set of chiral fields with ψ i ∼ [ω 2 10 ∨ ω 6 10 ∨ ω 8 10 ] under Z 10 and at least a set of chiral fields with χ i ∼ ω 5 10 under Z 10 , where ω |s| = e i 2π/|s| . The solutions with unconditional stability are highlighted with a bold font in the column s of Table 1. In this case, we have only a Z 10 solution in which we have a first set of singlet chiral fermions ψ i with charges (ω 2 10 , ω 8 10 ) and a second singlet chiral fermion χ with charge ω 5 10 . This give to arise a Dirac fermion DM candidate ψ 1 (ψ 2 ) † T protected by Z 5 and a Majorana fermion DM Finally, the type-I Dirac seesaw realizations of the effective operator of D-6 in Eq. (9) have to have a sufficiently rich h αi ν structure to explain the full neutrino oscillation data. That can be guaranteed by having a rank 2 or 3 Dirac neutrino mass matrix, via the inclusion of a proper set of chiral fermions for each solution. For example, consider the first solution (1, −4, −4, 9, 9, −11), which is promoted to (1, 1, 1, −1, −1, −4, −4, 9, 9, −11) by adding two pairs of extra states (1, −1). We assign m = 1, r = −1, l = −4 and ν = 9, such that s = −5 gives masses to two Dirac neutrinos. For avoiding the chiral fermion with charge −11 to be massless, an extra scalar S with X-charge 22 can be introduced to give a Majorana mass.

Type-II Dirac Seesaw
Contrary to type-I, type-II Dirac seesaw with only chiral fermions can be realized via D-5 operators. Figure 2 presents the diagrams realizing the D-5 (left panel) and D-6 (right panel) effective Lagrangian in Eq. (9) for type-II seesaw, if one only allows SSB masses for the singlet chiral fermions.
For the case D-5, the flux of the X-charges in each vertex for an active symmetry satisfy which implies that ν + m + s = 0 .
For D-6, the type-II seesaw simply implies a change of charge of s to s/2, as in the upper vertex there are two ongoing S instead of a single one.  Table 2: Type-II Dirac neutrino masses for an active symmetry: Set of charges satisfying the Diophantine equations together with the conditions enumerated in the text, for N extra singlet chiral fermions, featuring Dirac neutrino masses generated by D-5 operators. The solutions without massless chiral fermions are highlighted with a bold font.

Active Symmetry
The solutions of the Diophantine equations satisfying all the previously enumerated conditions for a type-II Dirac neutrino masses and an active symmetry, are shown in Table 2 for the D-5 operator. We note that even if most of the solutions contain massless chiral fermions, there is a single solution without (highlighted in bold). It corresponds to a case with N = 6 new chiral fermions, with the charge assignment (1, 2, −6, −6, − 6, 8, 9, 9, −11). It contains two RHNs with charge ν = 9, and three states with charge m = −6. The other four chiral fermions form two Dirac states (1, 2) and (9, −11) that obtain mass via the scalar S with charge s = −ν − m = −3, and that could be viable DM candidates. N Table 3: Type-II Dirac neutrino masses for a dark symmetry: Set of charges satisfying the Diophantine equations together with the conditions enumerated in the text, for N extra singlet chiral fermions, featuring Dirac neutrino masses generated by D-5 operators. None of the present solutions contain massless chiral fermions.

Dark symmetry
Contrary to the type-I seesaw, the type-II seesaw can accommodate a dark symmetry [16]. 4 For the case D-5, the flux of charges in each vertex is In general, for realizing the Type-II Dirac Seesaw, one requires with m for a gauge U (1) D symmetry, δ = 1 or 2 for D-5 or D-6 realizations, as in the upper vertex there are two ongoing S instead of a single one. The solutions of the Diophantine equations satisfying all the previously enumerated conditions for a type-II Dirac seesaw neutrino mass mechanism and an dark symmetry, are shown in Table 3, for the D-5 operator. A total of 19 solutions where found, 3 of them (highlighted in bold) without massless fermions. 5 Regarding these solutions without massless fermions, a few comments are in order: • The first solution corresponds to N = 6 and has the charge assignment (1, −2, −3, 5, 5, −6). It contains two RHNs with charge ν = 5. There are also two extra Dirac fermions (1, −6) and (−2, 4 We note that this case reduces to the one with an active symmetry in the limit L = 0 = m. 5 The first two solutions were recently presented in Ref. [16]. −3) that get mass via the SSB by the scalar S with charge s = −ν = −5, each of them being independent DM candidates protected by a residual Z 5 symmetry.
We also have four solutions with unconditional stability. Concerning the second with Z 6 ∼ = Z 2 ⊗ Z 3 , some phenomenological considerations are in order: the particle content of the model, along with the charges of the remnant symmetries, are presented in Table 4. The interaction Lagrangian includes where and After the spontaneous symmetry breaking, the neutrino mass matrix reads and therefore, v η must be small to allow for sizeable Yukawa couplings. This condition can be easily is expected to be much more smaller than v for large M η .
Contrary to the scotogenic model, the dark sector is completely independent of the heavy particles associated to the type-II Dirac seesaw. Moreover, because the Z 6 symmetry has two subgroups, the subsequent two dark sectors are completely independent between them, and therefore, we have at least two independent DM candidates. The first one, protected by the subgroup Z 3 , is the Dirac fermion ψ = ψ 1 (ψ 2 ) † T , while the second one, protected by the subgroup Z 2 , is the Majorana fermion Ξ = χ (χ) † T . For simplicity, we assume that other potential DM candidates are heavier than ψ and Ξ and have very small densities in the early universe, primarily by DM conversion into ψ and Ξ mediated by Z and S [30]. Besides the DM conversion processes, the proper relic density for the Majorana DM candidate through the annihilation ΞΞ → SS (m S < |m Ξ |), and compatible with direct detection constraints from PandaX-4T [42], was analyzed for the first model in Table 3 in Ref. [16]. The authors did the same analysis for the second model in Table 3 with one of Dirac DM candidates through the annihilation ψψ → Z D Z D (m D < |m ψ |).
We can also introduce an extra singlet scalar in the dark sector with a D-charge different from ±6, φ, to allow the Yukawa coupling between one of the heavy chiral fields with the right-handed neutrinos [43], as for example y φα ψ 1 ν Rα φ for a φ with D-charge 4. This coupling leads to a Dirac neutrino portal scenario which can explain fermion DM in the context of the type-II Dirac seesaw [44].

Conclusions
Studies on tree-level Dirac neutrino masses have typically focused on finding specific anomaly-free solutions for a given kind of symmetry, either for an active or dark symmetry. Alternatively, in the present work a complete set of relevant anomaly-free solutions to the general problem of the generation of Dirac neutrino masses at tree level with chiral singlet fermions has been presented. In particular, we restricted the analysis to solutions satisfying a set of general conditions enumerated in the text. An intensive scan was performed, looking for solutions generating neutrino masses by the type-I and type-II treelevel Dirac seesaw mechanism, via operators with dimension 5 and 6, that correspond to active or dark symmetries. Each of the presented solutions leads to a unique model with specific phenomenological implications.
It is interestingly to note that type-I Dirac seesaw can only take place for active symmetries, if one demands all extra fermions to be charged under the new symmetry. Additionally, type-I Dirac seesaw with dimension-5 operators automatically implies a vector-like pair of singlet fermions, and is therefore not possible with only chiral fermions. However, for dimension-6 operators we found a set of 36 solutions of the Diophantine equations (i.e. anomaly-free solutions) satisfying general conditions enumerated in the text (see Table 1). Among them, only 2 solutions with all extra fermions getting mass via the spontaneous symmetry breaking of the new Higgs field. The massless fermions of the other solutions can either contribute to the relativistic degrees of freedom ∆N eff in the early universe [23], or acquire masses after the introduction of an extra singlet scalar, becoming independent DM candidates [17].
Contrary to type-I, type-II Dirac seesaw with only chiral fermions can be realized via dimension-5 operators. For the case of an active symmetry and dimension-5 operators, 9 solutions were found, only one of them without massless chiral fermions (see Table 2). Alternatively, for the case of a dark symmetry, 19 solutions were found, 3 of them without massless chiral fermions (see Table 3).
Most of the solutions found could feature multicomponent DM. Special attention was brought to those with at least two DM candidates with unconditional stability, which guarantees the viability of DM without imposing any kinematical restriction.