Neutrino Masses in a Two Higgs Doublet Model with a U(1) Gauge Symmetry

General Two Higgs Doublet Models (2HDM) are popular Standard Model extensions but feature flavor changing interactions and lack neutrino masses. We discuss a 2HDM where neutrino masses are generated via type I seesaw and propose an extension where neutrino masses are generated via a type II seesaw mechanism with flavor changing interactions being absent via the presence of a U(1) gauge symmetry. After considering a variety of bounds such as those rising from collider and electroweak precision we show that our proposal stands as a UV complete 2HDM with a dark photon where neutrino masses and flavor changing interactions are addressed. A possible dark matter realization is also discussed.


Introduction
The Standard Model (SM) is the most accurate description of nature to the electroweak and strong interactions [1][2][3]. The discovery of a 125 GeV spin-0 state at CERN was the last piece of the puzzle in the SM [4,5] and established the existence of, as far as we know, an elementary scalar particle in nature. However, elementary scalar particles are common figures in many beyond the SM adventures, and among those Two-Higgs-Doublet stand out [6]. The ρ parameter, ρ = m 2 W /(m 2 Z cos 2 θ W ) is a powerful probe to those models featuring multiple scalar particles because they may contribute to the gauge boson masses [7] and, therefore, alter the SM prediction. The current value from global fits point to ρ = 1.00039 ± 0.00019 [3]. The gauge boson masses arise from the kinetic terms of the scalars, thus the ρ parameter can be parametrized at tree level as, where I i , Y i and v i are the isospins and hypercharges and vacuum expectation values of the scalars. From eq. (1.1) we can see that scalar doublets (I = 1/2) with Y = ±1 and scalar singlets (I = 0) with Y = 0 do not contribute to the ρ parameter, and for this reason are desired extensions of the SM.
In this work, we propose a new model, still within the scope of 2HDM but different from previous studies. We explain the absence of FCNI via an abelian gauge symmetry U (1) X and neutrino masses via type II seesaw mechanism [50][51][52]. The presence of a new abelian gauge symmetry gives rise to gauge anomalies which, in order to be canceled out, impose restrictions over the SM fermion charges under the new symmetry. Moreover, this abelian group induces the presence of a Z gauge boson. The scalar doublets develop vacuum expectation value (VEV) at the electroweak scale and the VEV of the scalar triplet cannot be large due to bounds stemming from the ρ parameter, hence the Z is necessarily light, with masses below the weak scale. This fact has important experimental implications which will be investigated. We will also briefly discuss how one could potentially accommodate a dark matter candidate in our model.
Our work is structured as follows: In section II, we introduce the 2HDM-U(1) model and revisit how neutrino masses could be easily generated via a type I seesaw mechanism; In section III we show how to implement a type II seesaw mechanism and obtain the mass spectrum of the model; In section IV we discuss some phenomenological constraints; In section V we draw our conclusions. At the end we left three sections in the appendix where details of the anomaly cancellation and spontaneous symmetry breaking mechanism are shown.
2 Type I seesaw in the 2HDM-U (1) As aforementioned, general 2HDM suffer from severe bounds rising from flavor physics. The problem of FCNI at tree level can be elegantly handled by the introduction of an extra abelian gauge symmetry. This new gauge symmetry is certainly more theoretically appealing than the usually ad-hoc Z 2 discrete symmetry which must be explicitly broken in order to not generate domain walls [53]. We will now revisit how the addition of an abelian gauge symmetry benefits 2HDM and generates neutrino masses via type I seesaw mechanism [54][55][56][57]. In this model, we have two scalar doublets Φ i ∼ (1, 2, 1, Q X i ) with the same hypercharge Y = 1 where, The fact that Q X 1 = Q X 2 leads to the scalar potential, Since the scalar doublets have different charges under U (1) X only one of them will couple to SM fermions, and we arbitrarily choose Φ 2 . In this way, we get the Yukawa Lagrangian, Thus far our model is nearly identical to the usual type I 2HDM model proposed previously in the literature. The key difference lies in the introduction of a gauge symmetry which naturally explain the origin of the Z 2 symmetry. Since we also want to accommodate neutrino masses via type I seesaw mechanism we need to add three right-handed neutrinos as follows, where we now included a singlet scalar Φ s , charged under U (1) X , to build a majorana mass term, which features a scalar potential, With these Lagrangians we can simultaneously explain neutrino masses and the absence of FCNI. However, there is still one point that needs to be addressed which is the presence of gauge anomalies arising from the introduction of a new gauge symmetry. Since the SM fermions can be charged under the new gauge group and we have introduced new chiral fermions (right-handed neutrinos), the anomaly cancellation procedure becomes non-trivial. Generally the Φ 2 charge under U (1) X is different from zero, and therefore the SM fermions should also be charged under U (1) X . The anomaly cancellation procedure is described in Appendix A. We highlight that we can in principle keep the model anomaly free without the addition of right-handed neutrinos. However, without them we would not be able to explain neutrino masses the way we wished for. Anyways, with their presence we can derive the anomaly cancellation requirements that preserve our Lagrangians as follows, where u and d are the U (1) X charges of the up and down quarks respectively, q (l) the charge of the quark (lepton) doublet, e (n) the charge of the right-handed charged leptons (neutrinos), and lastly Q X i (q X ) the U (1) X charge of the scalar doublet (singlet). The scalar singlet Φ s is responsible for generating a majorana mass term for the right-handed neutrinos and breaking the U (1) X gauge symmetry that yields a massive Z gauge boson. A natural question that rises to this conclusion is: is there a way to explain neutrino masses without adding a singlet scalar and right-handed neutrinos within the scope of 2HDM? Motivated by this question we will describe hereafter how one could accomplish that.
3 Type II seesaw in the 2HDM-U(1) A popular mechanism to explain the active neutrino masses without the presence of right-handed neutrinos is the so called type II seesaw mechanism [58,59]. In order to implement this mechanism within the scope of 2HDM the fermion charges under the gauge U (1) X symmetry need to be tied to one another to cancel out the triangle anomalies. We have seen above that one of the anomaly cancellation conditions is n = −(u + 2d), which comes from the U (1) 3 triangle anomaly as shown in the Appendix A, where n is the right-handed neutrino charge under U (1) X . Therefore, if there are no right-handed neutrinos we must set u = −2d to be free from gauge anomalies. Compared to the type I seesaw scenario, instead of having two independent charges (u and d), we now have only one, say d. That implies into, (3.1) The charge of the first doublet is free, as long as Q X 1 = Q X 2 , in order to recover the Yukawa Lagrangian (2.3) and keep the model free from FCNI. As shown in the Table 1, there is essentially only two different possibilities. One where the SM fermions are neutral and the other where they are charged under U (1) X . If a particular nonzero value is chosen for d, any other multiple of this value would produce a physically equivalent model, because a change in d can be balanced by a rescaling on the gauge coupling constant g X , so that the U (1) X interaction remains the same. In particular, taking d = −2/3 we notice that the charges of the SM fermions under U (1) X are similar to the SM weak hypercharge. In this way, it is clear that a type II seesaw realization in the 2HDM-U(1) gives rise either to a fermiophobic or a sequential Z boson.
The implementation of type II seesaw mechanism requires an SU (2) L scalar triplet ∆ ∼ (1, 3, 2, q Xt ), where the quantum numbers refers to the transformation Charges in Type II seesaw 2HDMs free from FCNI Fields properties under the symmetry group SU The field ∆ can be parameterized as, with, The SU (2) L symmetry allows the introduction of an interaction between ∆ and the leptons via, which requires ∆ to have hypercharge Y t = 2 and lepton number L t = 2, automatically forbidding interactions to quarks. The inclusion of eq. (3.4) implies, 2l + q Xt = 0, (3.5) and using eq. (3.1) we get, explaining the ∆ charge shown in Table 1.
Since ∆ carries lepton number, when the neutral scalar ∆ 0 develops a VEV, v t , lepton number is violated, and from eq. (3.4) we can easily see that it generates a majorana mass term for the neutrinos with, Thus, v t has to be very small in order to accommodate neutrino masses in the sub-eV range. In summary, with the presence of a U (1) X gauge symmetry, we can explain the absence of FCNI and accommodate neutrino masses via a type II seesaw without extra fermions, which is the main idea of this work. However, we need also to study the phenomenological implications of such proposal before concluding whether we have a feasible theoretical model. We start studying the mass spectrum of the model.

Mass Spectrum -Scalars
Our goal in this section is to study the spontaneous symmetry breaking pattern to find the physical scalars, gauge bosons, and neutrino masses. The SM charged lepton masses are the same as in the SM. That said, we begin our reasoning with the scalar sector.
The scalar sector is described by the Lagrangian, where the covariant derivatives of the scalar doublets and the triplet read, where τ a are the generators of the SU (2) L group. The scalar potential in eq. (3.8), invariant under all gauge symmetries is given by (3.11) Observe that the terms Φ T 1 iσ 2 ∆ † Φ 2 and Φ T 1 iσ 2 ∆ † Φ 1 are forbidden by the U (1) X symmetry, as we require Q X1 = Q X2 . There is only one non-hermitian term Φ T 2 iσ 2 ∆ † Φ 2 , which breaks lepton number in two units. Such lepton number violation is a common feature in seesaw type II models. It is important to note that neutrino masses are generated when ∆ 0 develops a vacuum expectation value as shown in eq.(3.7) and that would be related to lepton number violation since the scalar triplet carries lepton number. However, notice that the non-hermitian term in eq.(3.11) already explicitly violates lepton number, thus lepton number had been violated even before ∆ 0 develops a non-trivial vacuum expect value. We checked that without this nonhermitian term in the scalar potential the pseudoscalar from the scalar triplet field would remain massless, i.e. a majoron field [60].
Anyway, substituting the VEVs, in order to break spontaneously the gauge symmetries, we have the following constraint equations for a minimal point of the potential, In the Standard Model, the symmetry is spontaneously broken when the mass term flips sign. In the case of a type II seesaw, the mass term of the scalar triplet in the scalar potential does not need to flip sign to break the symmetry. We emphasize that the gauge symmetry is broken when the scalar doublets acquire a non-zero vacuum expectation value. See [61] for a detailed discussion about the type II seesaw vacuum. Moreover, we will see later on that the mass term of the scalar triplet should be positive in order to generate a pseudoscalar with positive mass.
Assuming that 2m 2 t is the dominant term between the brackets in the constraint equation (3.15), we have a seesaw relation, which leads to a naturally dwindled v t for |m 2 t | |µ t2 v 2 |. In this way, a small v t can be understood as a simply consequence of having the coefficient of the bilinear term in ∆ to be comparatively large with respect to the other energy scales of the scalar potential. Note that from eq. (3.16) we conclude that m 2 t and µ t2 should have the same sign.
In the scalar sector, Φ i and ∆ render the existence of seven physical fields: 3 CP-even scalars, h, H and H t ; one CP-odd, A; two singly charged H + , H + t and one doubly charged H ++ . The other scalar degrees of freedom are absorbed as longitudinal components by the gauge bosons, W ± , Z and Z , making them massive.
In the basis (ρ 1 , ρ 2 , ρ t ) the neutral scalars mix according to the following mass matrix, (3.17) From the diagonalization procedure of this mass matrix we will get three physical scalars, h, H and H t . We can parametrize this diagonalization in terms of three mixing angles α, α 1 and α 2 , where s α,α 1 ,α 2 and c α,α 1 ,α 2 are sine and cosine functions. We choose h to denote the 125 GeV SM-like Higgs found in the LHC [4,5]. The angles are determined by the parameters of the potential and the scalar VEVs. Fully analytic expressions for the masses and eigenvectors are complicated but we can obtain approximate results. As shown in the Appendix B.1, in the limit v t v i , the masses of the CP-even scalars are approximately, This limit v t v i that will be assumed throughout this work yields a higgs boson, h, with the correct mass is shown in Fig. 1. It is straightforward to see that we can easily find a higgs with the correct mass for couplings of order one and µ t2 either at the weak of multi-TeV scale.
As for the pseudoscalars, in the basis (η 1 , η 2 , η t ) the mass matrix is given by, Note that η 1 is decoupled and massless. Thus it can be immediately recognized as a Goldstone boson, G 1 . After the diagonalization procedure we found another Goldstone boson, G 2 . These two massless pseudoscalars represent the degrees of freedom needed to generate the Z and Z masses. In the diagonalization procedure we find the rotation matrix, where, which gives rise to two massless fields as aforementioned and a massive pseudoscalar, A, with mass, Observe that v t and µ t2 must have the same sign in order to have m 2 A > 0. We had concluded previously from eq.(3.16) that µ t2 and m 2 t should have the same sign to keep v t positive definite, thus from eq.(3.25) µ t2 must be positive to generate a positive squared mass for the pseudoscalar A. Hence, both µ t2 and m 2 t are strictly positive.
It is important to stress that even with the introduction of a new gauge symmetry, the pseudoscalar A, which is a common figure in 2HDM, remains in the spectrum. Under the assumption that v t is smaller than µ t2 the pseudoscalar can have a mass sufficiently large to evade existing bounds, as we shall discuss further.
The charged scalars mass matrix in the basis (φ + 1 , φ + 2 , ∆ + ) is, (3.26) The physical fields are given by performing the following rotation, (3.27) The Goldstone boson G + is absorbed by W + , and the physical states H + and H + t have masses, where, The doubly charged scalar ∆ ++ does not mix any other field. This mass eigenstate, which we will denote henceforth by H ++ , has a mass given by, In summary, the scalar mass spectrum is largely controlled by the relative sizes of v t , µ t2 and v i . As v i is fixed to be ∼ 100 GeV and, as we will see later on, v t is constrained to be O(1) GeV, we will always take v t v i . In this limit, the masses of h, H and H + are rather insensitive to v t and µ t2 , while the masses of H t , A, H + • µ t2 ∼ v t v i : In this case, µ t2 has little influence on the masses of h and H. H remains always heavier than h, with a mass around the 100 − 300 GeV range, for λ's of order ∼ 1. The masses of H t , A, H + t and H ++ are controlled by the ratio µ t2 /v t , with H t and A nearly mass degenerate. In particular, for µ t2 = v t , the masses are around 200 GeV. Such low masses can be dangerous in light of existing bounds [62][63][64][65][66][67][68][69][70].
• µ t2 ∼ v i : In this scenario the spectrum is shifted up and the scalar masses can be significantly larger than 100 GeV. A, H t , H + t and H ++ are mass degenerate and may reach masses in the TeV domain. For example, taking µ t2 = v 2 = 100 GeV and v t = 100 MeV, we obtain m h = 125 GeV, m H = 404 GeV, We have adopted λ 1 = 1.6, As stressed above, the masses of H t , A, H + t and H ++ are always close to each other because their masses follow m 2 Thus, for small v t , the m 2 t term is the dominant one, so that the masses are all approximately given by m 2 The mass splittings are controlled by the scalar VEVs. At leading order, (3.31) These mass splittings are noticeable only when the masses are small, i.e., for small µ t2 . For µ t2 v i , they are basically mass degenerate.
As aforementioned, the masses of h, H and H + are less sensitive to v t and µ t2 , and depend mostly on the VEVs v i and the λ i 's. Therefore, they naturally lie at the weak scale. As for H + , we find m 2 H + − 1 2 λ 4 v 2 (see eq. (B.36) in the Appendix B.2), which requires λ 4 to be negative. If we took |λ 4 | > 1, we would have charged scalar masses above 500 GeV, which can easily evade existing limits [71][72][73][74][75][76][77][78]. In order to have scalar masses above 500GeV we need couplings larger than the unit. This typically true in seesaw type II models.
Returning to the mass expressions of the CP-even scalars, eqs. (3.19) and (3.20), we see that ifthere were just the two scalar doublets, the neutral scalar masses would be given by these expressions with v t , µ t2 set to zero. However, the scalar triplet generate negative correction terms proportional to µ t2 v t , so that these scalars become lighter than they would be if there was not the triplet. However, as shown in Figure 1, the parameter space allows to fit a mass m h = 125 GeV for the Higgs boson h. Furthermore, the eqs. (3.19) and (3.20) imply an upper bound in the combination µ t2 v t : a large value of µ t2 must be balanced by a small value of v t in order to avoid negative squared masses for h and H. Looking at these equations, we conclude that in order to preserve the masses positive, v t must satisfy, where we assume λv ∼ 100 GeV. For instance, µ t2 ∼ 10 14 GeV implies v t 10 −1 eV. Thus, one may naturally generate small v t taking µ t2 at a Grand Unification scale [79]. Notice that eq. (3.32) is another kind of seesaw relation between v t and µ t2 valid in the limit v t µ t2 , which is independent of the relation in eq. (3.16).
Now that we have finished with the scalar sector, we will derive the masses of the gauge bosons.

Mass Spectrum -Gauge Bosons
The Lagrangian for the kinetic terms of the gauge fields associated to the hypercharge U (1) Y and the U (1) X symmetry is given by, where is the kinetic mixing parameter.
A canonical gauge kinetic Lagrangian is obtained through a GL(2, R) rotation on the fieldsB µ andX µ ,X so that the covariant derivatives (3.9) and (3.10) become, where G Xi = 2g Y i cos θ W + g X Q Xi and G Xt = 2g Yt cos θ W + g X q Xt .
After spontaneous symmetry breaking and performing the electroweak rotation, the spectrum of vector bosons turns out to be comprised of: the charged W ± µ ; the photon, A µ ; two neutral states, Z 0 µ and X µ , mixing to each other. They have the following mass Lagrangian where, with We see that the W ± µ bosons are already the mass-eigenstates with mass m W . The Z and Z gauge bosons on the other hand mix and lead to the following mass matrix, The diagonalization leads to, where, with ξ given by, This mixing angle is constrained to be very small by the LEP electroweak precision measurements on the Z boson pole [80]. Thus, Also, as we are interested in a light Z , we will assume the limit m 2 Z 0 X as well). In this limit, and we can write approximate expressions for the masses of Z and Z as , and, at leading order, m 2 In terms of β, defined by tan β = v 2 /v 1 (see Appendix B.2), Note that the presence of the triplet induces only a tiny correction proportional to (v t /v) 2 , so that the addition of a triplet scalar cannot generate a heavy Z , as opposed to the singlet case [49]. Thus, the Z mass lies below the electroweak scale, being controlled by the value of g X . For instance, taking tan β = 10 and Q X1 − Q X2 = 1, m Z varies from 1MeV − 1GeV, for g X in the range of 10 −3 − 10 −1 , regardless of the value of v t , as long as v t < 2 GeV.
In summary, we have proposed a type II seesaw mechanism for neutrino masses within the scope of 2HDM which prevent FCNI via gauge symmetries. Having discussed the mass spectrum of the model, we now will pay attention to some phenomenological constraints.

Phenomenological constraints 4.1 Electroweak Precision
The ρ parameter, which measures the relative intensity between the neutral and charged currents, is very accurately determined experimentally, ρ = 1.00039 ± 0.00019 [3] at 1σ level.
In the SM, the ρ parameter is equal to 1 at tree level, and its good agreement with the experimental value poses tight constraints on new physics models with extended scalar sector. In our model, the ρ parameter places an upper bound on the VEV of the triplet scalar because it contributes to the masses of Z and W ± bosons according to eqs. (3.38) and (3.48), translating into Hence at 3σ we obtain v t ≤ 2.3 GeV, where we used v 2 + 2v 2 t = 246 2 GeV 2 . As we are interested in a small v t for the generation of tiny neutrino masses, this constraint can be easily satisfied in our model. Notice that as v t becomes very small the scalar masses increase as can be seen, for instance in eq.(3.25) and eq.(3.30).

LHC
The U (1) X symmetry is spontaneously broken by the VEV of the doublets and the triplet, which also contributed to the mass generation of the Z vector boson. As v t is small and v is at the electroweak scale, the Z mass will be at the electroweak scale or below, depending on the value of g X and the other parameter such as tan β. Such a light Z is subject to a variety of experimental constraints. Notice that we have two possible Z models (U (1) N or U (1) Y ), one which resembles the sequential Z model, and other the dark photon model. Concerning the latter, LHC bounds are weakened because the Z − Z mixing is necessarily small, and that would suppress its production cross section at the LHC [81][82][83].
As for the U (1) Y model, we do not have much freedom since the SM fermions are charged under U (1) Y the production cross section is much larger. In this scenario the LHC bounds are rather restrictive. Assuming g X = 1 the LHC severely rules Z masses below 3 TeV [84]. In our model the Z mass is set by g X . In order to have Z masses around 100 GeV, g X should be around 0.1, which is not sufficiently small to evade LHC limits [85]. If we adopt g X = 0.01 we will get m Z = 1 GeV, and for such small coupling we can easily evade LHC limits [85]. We have used eq. (3.49) and assumed Q X1 of the same order of Q X2 to find the corresponding Z mass. We point out that the kinetic mixing parameter while not relevant for the Z mass it is important to determine the Z interactions with SM fermions. The conclusions drawn above are valid for sufficiently small kinetic mixing.
Regarding the scalar spectrum of our model, the most relevant ones come from LHC searches for heavy Higgs and triplet scalars. The cleanest signature signal is the doubly charged Higgs. We have then implemented the model in Madgraph [86,87] and followed the recipe described in [88]. Assuming no hierarchy in the Yukawa couplings the doubly charged scalar decays essentially, with equal branching ratios, into charged leptons. That said, we found the current LHC bound with L = 36f b −1 of integrated luminosity and performed future projects for the High Luminosity and High Energy LHC setups as summarized in the Table 2. In the light of current bounds our model is agreement with existing bounds if we take µ t2 ≤ v i which predicts masses at the TeV scale as discussed previously. We highlight that we need couplings larger than one to find charged scalar masses above the TeV scale. Therefore, LHC and its planned upgrade will be important because it will probe a large portion of the model. The presence of a doubly charged scalar is the key signature of the type II seesaw mechanism. If we find ourselves in a situation where both doubly charged scalar and Z fields are observed, our model stands as a potential well motivated environment for them.

Belle-II and KLOE2
Belle and KLOE collaborations represent e + e − colliders searching for light gauge bosons with the /2F µν F µν . The two models proposed here feature a similar term. In the U (1) N model SM fermions are uncharged under U (1) N , thus the Z will couple to SM fermions only via its mixing with the Z boson generated by the presence of the Table 2: Summary of collider bounds on the doubly charged scalar in our model using current and planned configurations. We used 13TeV of center-of-mass energy for the LHC configurations, whereas 27TeV for the high-energy upgrade. We can see that LHC and its upgrade will be paramount to probe the model up to the TeV scale. kinetic mixing. In this case, our model would a UV complete version of the simplified dark photon model [89,90]. This scenario for heavy Z masses was investigated in [91]. For the U (1) Y model, where the SM fermions are charged under the gauge symmetry, if we take g X 1 and g X < , again the model falls back to the dark photon model because the Z interactions to SM via the kinetic mixing would more pronounced, the experimental limits on dark photon become applicable to our study.
In summary, the experimental limits derived for dark photon models apply here, except in the case where g X and m Z 1 GeV. Experimental collaborations usually display their bounds in terms of 2 . In Figure 2 we display a summary of the existing (gray) and planned (color) constraints. For m Z ∼ 10 − 30 MeV, current bound impose < 10 −4 , limiting the region of parameter of our model. Anyhow, we emphasize that we can still obey such bounds by taking g X and to be very small as it is usually assumed in dark photon models.

Accelerators
There are several accelerators using electron/positron or hadronic beams which search for bremsstrahlung of dark photons or its appearance in meson decays. These bounds are inside the gray region in Figure 2. It is important to point out the future sensitivity on the flavor violating decay µ → 3e [92], in case of no signal, will give rise to the upper limit in cyan. Moreover, The Heavy Photon Search Experiment (HPS) which was already installed at SLAC collides highly energetic electrons into a tungsten target, and in the process electrons may radiate dark photons. The experimental sensitivity of HPS is shown in blue.

Low Energy Probes
The muon anomalous magnetic moment (g-2) [93], neutrino-electron scattering [94] and atomic parity violation [49] provide complementary but subdominant limits to our model. For instance, neutrino-electron scattering rules out > 10 −5 [94]. One cannot accommodate g-2 with the U (1) Y model because the electrons are charged under the gauge symmetry, and these couplings are subject to tight constraints [95][96][97]. In the U (1) N model, where our model resembles the dark photon one, the favored region to explain g-2 has already been excluded [93].

Dark Matter Possibility
In the two models described here, there are no dark matter candidates. One could simply add a vector-like fermion charged under the gauge symmetry while preserving the gauge anomaly cancellation. The dark matter relic abundance, direct detection and indirect detection signals would be governed the kinetic mixing term and the g X gauge coupling. For concreteness, if we take the U (1) N model, where SM fermions are not charged under the new gauge symmetry, the dark matter phenomenology would be similar to the dark photon portal investigated recently in the literature [98]. It has been shown that in this setup if the dark matter mass is smaller than the Z one, only s-channel processes would be relevant, and this case is nearly ruled out by current data for dark matter masses above 10 MeV. There is tiny region for ∼ 10 −5 and m Z ∼ 10 MeV which obeys on existing limits and reproduce the correct relic density. If the dark matter particle is heavier than the Z then the secluded dark matter setup arises [99], scenario which is much less restricted by data. In the near future we plan to carry out a detailed dark matter phenomenology in both models taking into account the particularities such as the presence of both mass mixing and kinetic mixing terms as well as the existence of many other scalars which might alter the overall predictions.

Discussion
There are important things to be stressed about the two models we proposed to explain neutrino masses and the absence of flavor changing interactions: • The gauge symmetry imposed to distinguish Φ 1 from Φ 2 and then allow just one scalar doublet to couple to SM fermions gives rise to two very different type of models. In the U (1) N model, the SM fermions are uncharged under the gauge symmetry, and the corresponding massive Z only couples to SM fermions via its mixing with the Z boson. In the U (1) Y setup, the Z gauge boson will have a neutral current with SM fermions determined by the U (1) Y SM fermion charges, leading to a sequential Z model [100].
• Since we have added a triplet scalar to explain neutrino masses via a type II seesaw mechanism, nothing prohibits from one to consider off-diagonal yukawa couplings, involving for instance the scalar triplet, to be non-vanishing. This will lead again to flavor changing interactions and give rise to µ → eγ, µ → 3e, µ − e conversion processes which are rather restricted by data. In particular, the product of the yukawa terms are limited to be smaller than 10 −7 [93]. Anyways, this feature is common in the models which extend the SM scalar sector and, thus, we can set to zero off-diagonal couplings involving the extra scalars without prejudice.
• The additional gauge symmetry allow us to easily introduce a dark matter candidate, a vector-like fermion, without spoiling the anomaly cancellation requirements. In summary, we argue that the addition of a gauge symmetry and a triplet scalar is well-motivated since it adds nice features to the original 2HDM proposal.
• The bounds discussed previously can be safely satisfied by taking g X and to be sufficiently small, below 10−3 similarly to dark photon models.

Conclusions
Two Higgs Doublet Models represent interesting extensions of the Standard Model. However they lack neutrino masses and a first principle explanation for the absence of flavor changing neutral interactions. Typically, a Z 2 symmetry is invoked where one scalar doublet is charged under, to prevent the other to couple to Standard Model fermions.
In this work, we proposed two models where neutrino masses are explained within the type II seesaw mechanism via the addition of a scalar triplet and a gauge symmetry that allows only one scalar doublet to couple to fermions. In this way can simultaneously explain neutrino masses and avoid flavor changing neutral interactions.
We have investigated several constraints coming from low energy probes, electroweak precision and collider. In particular, we derived collider bounds on the mass of the doubly charged scalar using current and planned LHC reach with highluminosity and high-energy configurations. We discussed which regions of parameter space are consistent with current data to conclude that both models stand as viable alternatives to the original Two Higgs Doublet Model proposal.

Acknowledgments
The authors are grateful to Carlos Yaguna, Carlos Pires and Diego Cogollo for their comments. D.A.C. acknowledges support from UFRN and MEC while T.B.M. thanks CAPES for financial support. A.G.D. is partially supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) under the grant 306636/2016-6. F.S.Q. is partially financed by ICTP-SAIFR grant 2016, UFRN and MEC. A.G.D. is grateful to UFRN for the hospitality during the early stages of this project.

A Anomaly Cancellation
In this section we will describe how the gauge anomalies can be cancelled in our model which simply adds an abelian gauge group to the SM. Taking the SM fermion charges under the new weak hypercharge, Y , to be l (left-handed leptons), q (left-handed quarks), e (right-handed leptons), u (right-handed quarks with positive isospin), d (right-handed quarks with positive isospin). Let us consider each relevant gauge anomaly individually, We conclude that, Consequently we get, These relations between the charges under U(1) are general, since we have not added any condition concerning the Lagrangians of the model. If we want to accomodate neutrino masses either via type I or type II seesaw mechanisms and explain the absence of flavor changing interactions, further relations between hypercharges under the new gauge symmetry U(1) arise as summarized in eq. (2.6) and eq.(3.1).

B Scalar masses and mixings
In this appendix we show in details how to obtain approximate expressions for the masses and mixing angles for the physical scalars in the relevant regime adopted in this paper, v i ∼ 100 GeV, v t v i . In the first section we treat the CP-even scalars, and in the next we treat the charged scalars.

B.1 CP-even scalars
In general, it is not possible to obtain analytic expressions for the diagonalization of the mass matrix M 2 CPeven , given in eq. (3.17). The masses of the neutral scalars (eigenvalues of that matrix) are given implicitly as the solutions of the polynomial equation, with a, b, c and d given by, This equation can be solved numerically once the set of parameters is fixed. For the mixing angles, it is very difficult even to furnish an equation that determine them in terms of the parameters of the potential, because of the difficulty in computing the eigenvectors of M 2 CPeven . There are some limits, however, in which these expressions are calculable. The idea is to take advantage of the different energy scales involved and decompose the original matrix into matrices whose entries belong to the same scale.
Let us decompose M 2 CPeven , in the following way, The matrix M 2 1 would be the mixing matrix of the neutral scalars in case there were just the two doublets, while M 2 2 account for the effects of the presence of the triplet. Let us first consider the limit v t , µ t2 v i . In this case, the decomposition (B.6) makes it clear that the triplet decouples from the doublets. The matrix M 2 1 remains the same and M 2 2 reduces to M 2 2 = diag(0, 0, µ t2 v 2 2 / √ 2v t ), from which we obtain immediately the mass of H t . To diagonalize M 2 1 we need only one mixing angle, so that we can make α 1 , α 2 → 0 in eq.(3.18), leading to the following physical fields, with α given by, From the eigenvalues of M 2 1 and M 2 2 , we have the masses, Now, relaxing the condition on µ t2 and allowing it to increase to the same order of v i or higher, this comparatively large value of µ t2 produces a sizable perturbation on the spectrum obtained above, but as we will see, the mass expressions are somewhat similar to the ones obtained in eqs. (B.8) and (B.9). In this case, is still possible to diagonalize the matrices M 2 1 and M 2 2 almost independently. First, as v t is always taken to be small but now µ t2 can be large, M 2 2 can be approximated by, M 2 2 is diagonalized by moving to an intermediate basis ( so that, The effect of this rotation on M 2 1 is (B.14) As sin α 2 1, at leading order we have R α 2 M 2 1 R T such that, where m 2 h and m 2 H are as given in eq. (B.8). This second rotation R α , albeit diagonalize M 2 1 , tends to disturb the previous diagonalization of M 2 2 , by generating corrections to the diagonal elements and also off-diagonal elements in M 2 2diag , This matrix is diagonal up to O(v t /v 2 ) terms, which can be discarded for sufficiently small v t . As we are interested in extracting the leading order correction to the masses of h and H, we shall keep the diagonal elements.
In summary what we have done is, From, we can read the scalar masses, Note that the expressions obtained agree with those given in eq. (B.8)-(B.9) if we take µ t2 v i , as expected. The main difference in the expressions in these two limits is the presence of the correction terms −2 √ 2 sin 2 αµ t2 v t and −2 √ 2 cos 2 αµ t2 v t in m 2 h and m 2 H , respectively, which pushes down their values, making h and H lighter than it would be in the absence of the triplet. Notice also that in this approximation, we managed to perform the diagonalization using only two mixing angles, α and α 2 , instead of the three angles needed in the general case.

B.2 Charged scalars
In this section we will apply to the charged scalars mass matrix the same method used in the previous section for the neutral scalars. First note that the matrix (3.26), can be decomposed as,

(B.22)
In the matrix M 2 2+ are contained all the mixing effects among the doublets and the triplet. If M 2 2+ vanished, there would be mixing only between the two doublets, as described by M 2 1+ , leading to a charged Goldstone boson and a charged physical scalar, as in the usual 2HDM. Let's again consider the limit v t , µ t2 v i , in which M 2 1+ remains the same and M 2 2+ reduces to M 2 2+ = diag(0, 0, √ 2µ t2 v 2 2 /2v t −λ tt1 v 2 1 /2− λ tt2 v 2 2 /2), so that the triplet completely decouples from the doublets, which still mix with themselves. In this case, it is necessary only one angle in the diagonalization and we can make β 1 , β 2 → 0 in eq.(3.27), leading to the following physical fields, which can be put in the form tan 2β = 2 tan β/(1 − tan 2 β) by dividing numerator and denominator by v 2 1 , so that The masses of H + and H + t in this approximation are, As mentioned early, for µ t2 v t the mass of H + t is small and may be in tension with existing bounds, so that in this limit µ t2 /v t > 1 is favored. Now, allowing µ t2 be large but keeping v t small, M 2 2+ reduces to, The effect of this rotation on M 2 1+ is As sin β 2 1, at leading order we have R β 2 M 2 1+ R T β 2 M 2 1+ , and the rotation R β 2 does not change M 2 1+ , as we wanted. Then, rotating M 2 1+ by a matrix R β , with β given by eq. (B.24), we move to the physical basis, from which we obtain the masses for H + and H + t , As a consistency check, notice that taking the limit v t v i directly in the eqs. (3.28) and (3.29), we obtain as result the eqs. (B.36) and (B.37). Finally, note that for the diagonalization in this limit we need to use only two mixing angles, β and β 2 , instead of the three angles needed in the general case.