Reviving sub-TeV SU ( 2 ) L lepton doublet dark matter

In this work we study the hybrid kind of dark matter (DM) production mechanism where both thermal and non-thermal contribution at two different epochs set the DM relic abundance. This hybrid set up in turn shifts the parameter space of DM in contrast to pure thermal DM scenario. We review such production mechanism in the context of the SU ( 2 ) L lepton doublet dark matter ((cid:2)) augmented with an additional singlet dark scalar ( S ). The neutral component of the dark doublet can serve as a stable DM candidate and in pure thermal scenario, it is under-abundant as well as excluded from direct detection constraints due to its strong gauge interactions in the sub-TeV mass regime. However, in addition to the thermal contribution, the late time non-thermal DM production from the decay of the long-lived dark scalar S helps to fulﬁll the deﬁcit in DM abundance. On the other hand, the strong gauge mediated direct detection constraint can be evaded with the help of a SU ( 2 ) L triplet scalar(with Y = 2 ), resulting a pseudo-Dirac DM. To realize our proposed scenario we impose a discrete Z 2 symmetry under which both (cid:2) and S are odd while rest of the ﬁelds are even. We ﬁnd the lepton doublet pseudo-Dirac DM with mass ∼ 450 − 1200 GeV, compatible with the observed relic density, direct, indirect, and existing collider search constraints.


Introduction
The existence of dark matter (DM) is strongly suggested by several astrophysical and cosmological observations at a wide range of length scales concluding that about 80-85% of total matter density is dominated by DM [1][2][3][4].In spite of this great observational evidence, the Standard Model (SM) of particle physics at present set up fails to explain the particle DM.Additionally, the SM also can not explain the exisa e-mail: spspg2655@iacs.res.in(corresponding author) b e-mail: skjeesun48@gmail.comtence of neutrino mass and mixing as suggested by different neutrino oscillation experiments [5][6][7][8][9][10].To overcome these two shortcomings of the SM, various beyond standard model (BSM) scenarios have been proposed.The issue of neutrino masses and their mixing angles can be resolved by the three seesaw mechanisms [11][12][13][14][15][16][17][18][19][20].Weakly interacting massive particles (WIMP) [21][22][23][24][25][26] is the most popular and widely studied thermal DM candidate whose interaction strength with SM particles is assumed to be of the order of electroweak interactions to explain the observed relic density.However, the null results at various direct detection experiments [27][28][29][30][31][32][33][34][35][36][37][38][39] open the possibilities of alternate paradigms to explain the DM relic density.Among the alternate theories, the feebly interacting massive particle (FIMP) is a very attractive viable candidate whose interactions with SM plasma are considered to be too small to keep them in a thermal bath [40][41][42][43][44][45].Rather FIMPs are produced non-thermally from the decay or annihilation of bath particles and number density gradually freezes in and such a scenario is known as non-thermal DM scenario.Such tiny interactions of FIMP with SM particles can be the possible reason for the non-observations of DM in different detection experiments such as Panda [27], XENON [28], LUX [29].
Despite the fact that the dedicated direct search experiments have put the thermal WIMP paradigm in a corner and non-thermal DM scenarios are being explored widely as an alternate theory, yet there are several approaches to revive thermal DM scenario itself [46][47][48][49][50][51][52].Examples include strongly interacting massive particle (SIMP) [46], assuming non-standard cosmology [47], bouncing dark matter [49], freeze out from inverse decays [50].Among the different methods to resurrect thermal DM scenarios with different phenomenological implications, we are interested in a hybrid scenario where both the thermal and non-thermal contributions set the DM relic abundance.In such cases, the DM candidate can have sizeable interaction with SM bath particles giving rise to suppressed number density but non-thermal production from another source helps to meet observed relic giving rise to different observational consequences [51,53,54].At this point, it is worth mentioning that there exist some well-motivated particle physics models which can account for single component thermal DM but suffers from the under-abundance issue in a certain mass range i.e. can not explain 100% of the observed relic density [55][56][57][58].And for such models, the earlier mentioned hybrid scenario is beneficial where apart from the thermal density additional contribution is required to satisfy the whole fraction of DM relic density.In this work to realize such a mechanism we target the vector-like SU (2) L lepton doublet dark matter model [58].The striking feature of such a model is that it has only one free parameter, the DM mass.The interaction couplings with SM particles are gauge couplings and are not free parameters.However, in such a scenario, the strong gauge-mediated interaction of DM with SM bath leads to suppressed number density making it under abundant for masses below 1.2 TeV [59]. 1 We are interested in the sub-TeV range of SU (2) L lepton doublet dark matter and revisit such setup whether non-thermal contribution from additional sources can elevate the under-abundant region.
As mentioned above, our aim in this work is to study the phenomenology of a vector like SU (2) L lepton doublet, = (ψ 0 ψ − ) T where the neutral component, ψ 0 serves the role of the dark matter.In order to realize sub-TeV lepton doublet type DM we extend the dark sector with another dark real scalar, S. The stability of the lightest dark sector particle, ψ 0 which acts as DM can be ensured with the help of additional discrete symmetry, Z 2 , under which both the fields ( , S) are odd and the rest of the particles are even.The scalar S has Higgs portal coupling with SM particles due to which S freezes out from the thermal bath.On the other hand, ψ 0 also freezes out from the bath decided by all number-changing processes including co-annihilation.Due to the chosen parameter space S freezes out at an earlier time than the time when ψ 0 freezes out.However, S has Yukawa coupling with ψ 0 and for M S > M ψ 0 , S can decay to ψ 0 increasing ψ 0 number density.To repopulate ψ 0 number density from the non-thermal late decay of S the Yukawa coupling of S with ψ 0 should be 10 −9 .And throughout our analysis we choose it to be 10 −10 for not tampering with the prediction of light abundances from Big Bang nucleosynthesis (BBN) [70].We discuss all the observational constraints like direct detection, indirect detection, and collider search for our setup.
The remaining part of this paper is organized as follows.First, we discuss the general framework of the DM production mechanism in Sect. 2. In Sect.3, we study the phenomenology of Dirac lepton doublet DM model.We then present a detailed discussion of the pseudo-Dirac lepton doublet DM, including DM abundance, direct, indirect, and collider search constraints in Sect. 4. Finally, we conclude in Sect. 5. We also include some analysis relevant for our discussion in Appendices A-D.

A general framework of hybrid set up
In this work, our goal is to investigate the under-abundant region of the SU (2) L lepton doublet DM via late-time DM generation from a heavy thermal dark sector particle.To alleviate the under-abundance issue we discuss here a hybrid type of DM production mechanism with the help of another heavy dark sector particle which was in also thermal equilibrium at the early time of the Universe and yields non-zero number density via the freeze-out mechanism.Before going to the relic density analysis for specific SU (2) L lepton doublet DM we portray here a general discussion of the hybrid DM production mechanism. 2he general framework of the DM production is shown by cartoon diagrams in Figs. 1 and 2.Here A is the heavier dark sector particle and B plays the role of the stable DM candidate.In Fig. 1 we show the interaction for such an oversimplified setup.Initially, at very high temperatures both A and B were in thermal equilibrium with the bath particles with the help of the possible interactions A A ↔ SM SM and B B ↔ SM SM respectively.First, the heavier dark sector particle A is thermally produced which we call phase-I.In the next step, DM B freezes out from the thermal bath which is named as phase-II.At this point, if there were no interaction between A and B both of them could serve as two component dark matter.But here due to the presence of the non-thermal decay of A(A → B + SM) finally the number density of A is further added up to the abundance of DM B in phase-III.Note that in this case, particle A is not dominating the energy density of the universe at the time of decay.Hence we can neglect the entropy injection in the pre-BBN era due to the late decay of A [73].In certain scenarios where the mother particle holds a significant portion of energy density, the decay products can have impact the freeze out abundances, as discussed in Ref. [74] in contrast to our case.
For this general framework, we briefly discuss the numerical approach for evaluating the abundance of DM, B (M A > M B ).The evolution of co-moving density, Y i (Y i = n i /s where n i , s are number density and the entropy density) for both the components A and B in the early universe as a function of temperature can be described by the following coupled Boltzmann equation (BEQ), where the function G(z) defined as: . z = M sc /T is a dimensionless variable and M sc is an arbitrary mass scale to scale the temperature.
H is the Hubble expansion rate and is defined as GeV [21].g * ρ and g s are the degrees of freedom associated with energy density and entropy density respectively [21].Y eq A and Y eq B are the equilibrium co-moving number density of the species A and B respectively.The thermal average cross-section of the annihilation process X X ↔ SM SM (with X = A, B) is denoted by σ v X X↔SM SM and A→B SM signifies the thermally averaged decay of A → B + SM.The annihilation process A A ↔ SM SM is responsible to keep the particle A in the thermal bath at the early universe while for B it is B B ↔ SM SM.The evolution of the co-moving number densities for such a general setup is shown in Fig. 2. The blue and red lines indicate the number densities for species A and B respectively.The dashed and solid lines correspond to Y eq i and Y i (i = A, B) respectively.As the temperature drops below around M i , respective i particle density freezes out from the thermal bath and yields non-zero density, Y i .Due to the chosen parameter space, the heavier particle freezes out early followed by the lightest dark sector particle which is named as phase-I and phase-II respectively in Fig. 2.And finally due to the chosen dark Yukawa coupling the density of heavier species A dilutes to B which is phrased as phase-III in the above Fig. 2. Therefore one can take the final abundance of B approximately as ) be the co-moving number density of ith particle at the time of freeze-out So the relic density of B will be Note that the DM observed abundance measured by the Planck collaboration to be DM h 2 = 0.120 ± 0.001 [4].

Dirac dark matter
To realize the hybrid setup of DM production discussed in the earlier section we extend the SM particle content with a vector like SU (2) L lepton doublet, = (ψ 0 ψ − ) T and an additional real scalar singlet, S ,where the neutral component ψ 0 can act as a DM.For the minimal case (without S), the lepton doublet is the only dark sector field, and the neutral component serves as a stable DM candidate, stabilized by an extended Z 2 symmetry.It has been found that the neutral component of lepton doublet alone can not accommodate the observed relic for DM mass M 1.2 TeV [59].The reason behind the under-abundance issue is the large interaction rate with SM particles due to the gauge interactions.For M 1.2 TeV DM ψ 0 becomes over-abundant.To probe the under-abundant lepton doublet DM mass (sub-TeV) region we propose the hybrid setup and augment the particle content with additional dark sector particle S. Both the dark states ( and S) are odd and the SM particles are even under the additional Z 2 symmetry to ensure the stability of DM.The relevant part of the Lagrangian involving DM along with heavier dark scalar S is described as follows, The first term of the Lagrangian, L involves the interaction between the SU (2) L lepton doublet and the SM gauge bosons.Here g 2 and g 1 are the gauge coupling corresponding to the SM SU (2) L and U (1) Y gauge groups respectively while the corresponding gauge bosons are denoted by W μ and B μ .L i (i = e, μ, τ ) represents SM lepton doublet whereas is SM Higgs doublet.The gauge interactions are responsible for the thermal freeze out of the neutral component lepton doublet DM, ψ 0 .M is the bare mass of the lepton doublet.The masses of the neutral(ψ 0 ) and charged(ψ ± ) dark leptons can be expressed as [75], Here δm is mass the splitting between ψ ± and ψ, generated from 1 loop quantum correction as [75]: and generally is of the order of pion mass (δm ∼ m π ).L S is the interaction Lagrangian for the dark real scalar singlet,S where S is interacting with SM through Higgs portal interactions as S 2 † .With M S > 0 and λ S > 0, the scalar field S does not acquire any vacuum expectation value (vev), thus the Z 2 symmetry remains intact.M S is the physical mass of S.
The last term of the Lagrangian, L −S involving the Yukawa interaction between dark sector particles( and S) and the SM leptons is the most important term for our analysis.Depending on the value of y s the Yukawa interaction L −S will give rise to late time re-population for the lightest dark sector particle (ψ 0 or S).We assume the mass of the scalar S is heavier than the masses of the fermion states i.e.M S > M ψ ± + M .With the mass hierarchy the heavier dark sector particles can decay into the lightest state, ψ 0 as S → ψ 0 ν and S → ψ − + → ψ 0 π − + , where l denotes charged leptons.As a result the lightest neutral fermion, ψ 0 acts as a stable DM candidate thanks to the Yukawa coupling, y s .We choose small values of y s such that the decays start after both S and ψ 0 freeze out from the thermal bath and also before the BBN which will be addressed shortly.Mapping the dark sector with the discussion made in Sect. 2 one notices that S and ψ 0 resemble A and B respectively.For simplicity, we consider a universal Yukawa coupling of S (y s ) with the lepton sector.The corresponding decay widths of the heavy dark scalar, S to ψ 0 + ν and ψ ± ∓ is given by the following: We shall now identify the relevant parameters which control the phenomenology of lepton doublet DM in this framework.The model has mainly four relevant free parameters and they are as follows: In the rest of our discussions, we consider that M S > M , with ψ 0 serving as a DM.In addition to the aforementioned free parameters, the interaction between and the Standard Model (SM) sector is controlled by the SM gauge couplings (g 2 and g 1 ).The interaction between S and the SM is determined by the Higgs portal coupling, λ S .The Yukawa coupling y s plays a significant role in determining the repopulation of dark matter ψ 0 from S, as well as setting the lifetime of S (τ S ).Before delving into the detailed analysis of dark matter (DM), let us first provide a brief overview of various experimental limits imposed on the model parameters.For M < M Z /2, the invisible decay width of the Z boson receives additional contribution beyond what is predicted by the SM.However, the observational data on the invisible decay width of the Z agrees remarkably well with the SM prediction, thus strongly necessitating M > M Z /2 [10].Additionally, LEP-2 already excluded masses for exotic charged fermions, M ψ ± , below ∼ 102.7 GeV [76].In this framework, both the charged and neutral leptons are degenerate, meaning M ψ ± ≈ M ψ 0 ≡ M .As a result, the contribution of the lepton doublet to the EW precision parameters (S, T, and U) is consistent with the observed bound [77].Furthermore, the Higgs invisible decay width, as measured by LHC [78], imposes constraints on the Higgs portal coupling λ S when M S < m h /2.However, in this framework we consider that M S > M , thus rendering the Higgs invisible decay constraint inapplicable.Hence, in our discussion, we establish a lower bound on the mass parameters as follows: M S > M > 102.7 GeV.We will also consider the limits imposed by the observed DM relic abundance DM h 2 = 0.120 ± 0.001 by the PLANCK Collaboration [4], as well as direct search constraints from XENON [28], PANDA [27], and recent LZ data [79], along with indirect search constraints from FERMI-LAT [34] and MAGIC [35], on the model parameters.
The neutral component of the lepton doublet, with the mass hierarchy M < M S behave as a stable DM candidate.The relevant Boltzmann equations are where x = M S /T is the dimensionless variable.Y S is the comoving abundances of S. Note that in Eq. ( 10) we write the Boltzmann equation for Y signifying the total abundances of ψ 0 and ψ ± .But due to the presence of strong co-annihilation between nearly mass degenerate charged and a neutral component of like ψ 0 ψ ± , ψ + ψ − → SM SM, Y ψ ± becomes zero at the moment of freeze out leaving only non zero Here, x F = M S /T FO and T FO denotes the freeze out temperature of respective species.Y eq and Y eq S are the equilibrium co-moving number density of and S respectively.σ v S denotes the thermal average crosssection of S to the bath particles (SM), S S → SM SM.And σ v eff represents the effective thermal average cross-section of the dark lepton doublet associated with the annihilation processes, ψ 0 ψ 0 → SM SM and the co-annihilation processes, ψ 0 ψ ± , ψ + ψ − → SM SM.The effective thermal average cross-section, σ v eff can be expressed as follows [80,81]: The internal degrees of freedom g ψ 0 and g ψ ± are associated with the dark lepton states, ψ 0 and ψ ± respectively.The thermal average of the total decay width of S, ( (S → ± ψ ∓ )+ (S → ψ 0 + ν)) is denoted by S→ .It is important to note that for numerical analysis we have adopted the σ v eff as a function of temperature for a given set of parameters using the open-code micrOmega [82].To generate the model files for micrOmega, we first implement the model in the public code FeynRule [83].Using the σ v eff (T ), the analytical expression of σ v S for the annihilation processes of the Higgs portal S, S S → SM SM [84], and the analytical expression of S→ (in Eq. ( 7)), we solve the above coupled BEQ as described in Eqs. ( 9) and (10).
The evolution of the co-moving number densities of S and ψ 0 are described by Eqs. ( 9) and ( 10) respectively.The first term on the right-hand side (R.H.S) of Eq. ( 9) decides the freeze out of S as the temperature drops T < M S .Similarly, the first term on the R.H.S of Eq. ( 10) decides the freeze out of ψ 0 as the temperature drops T < M 0 ψ with Y ψ ± = 0 at the time of freeze out.However, the presence of a tiny Yukawa coupling y s among fermion DM, heavy dark scalar, and SM leptons gives rise to the late time decay of S into fermion DM and eventually dilutes the number density of S as depicted by the second term of Eq. ( 9).Similar to this, the second term in the R.H.S of the Eq. ( 10) governs the re-population of ψ 0 from the decay of S. Note that from the late decay of S, ψ ± is also produced.But due to the small mass splitting between ψ ± and ψ 0 (δm ∼ O(m π )), ψ ± promptly decays to ψ 0 (ψ ± → ψ 0 + π ± ) [85] and eventually all the number densities get converted to ψ 0 density.As previously stated, we set the coupling y s and the mass M S in such a way that the decay is active after the freeze out of ψ 0 and after the decay gets completed Y S is totally converted to Y ψ 0 .It is worth noting that in the absence of the Yukawa interaction, both ψ 0 and S becomes stable and can act as two component DM [86].However, we are interested in sub-TeV fermion doublet DM, which can accommodate the entire observed DM abundance by PLANCK [4] in the presence of a heavy dark scalar, S, with the tiny Yukawa coupling y s .Before going to numerical solutions of Y S , Y ψ 0 we mention some important issues regarding y s in the next paragraph.
While dealing with such late time decays, S → ± ψ ∓ , ν ψ 0 (ν ψ 0 ), one should be careful so that the decay products do not alter the predictions of BBN.The charged fermion ψ ± decays further to π ± and ψ 0 with 97% branching ratio and even the pions can decay to leptons.These excess pions and leptons may interact with nucleons during BBN and tamper the abundances of light elements attracting strong constraints [87].To be on the safer side and simplify our analysis, we restrict our analysis by demanding that the decay is completed at a much higher temperature than BBN (i.e.τ S t BBN ).In Fig. 3, we show the parameter space in M S −M vs. τ S plane where τ S is the life time of S. Here, we consider three different values of y s = 10 −10 , 10 −11 , 10 −12 shown by blue, red, and green lines respectively.The solid and dashed line correspond to M S = 1000 GeV and 2000 GeV respectively.The gray shaded region corresponds to the parameter space where the decay of S will occur after BBN.According to Eq. ( 7), the lifetime of S is inversely proportional to y s and M S − M .This feature is evident from the aforementioned figure.We also want the coupling y s to be such that the decay of S starts non-thermally after the freeze out of ψ 0 and for that reason, y s should be 10  The solid and dashed line correspond to M S = 1 TeV and 2 TeV respectively.The shaded region corresponds to where the decay will be active after BBN the decay of S also gets completed before BBN.However, apart from setting the lifetime, y s plays no role in deciding the relic which we will discuss in the following paragraph [88].
We shall now discuss how the model parameters affect the DM number density obtained from the BEQs in Eqs. ( 9) and (10).In Fig. 4, we showcase the evolution of comoving abundances of S and ψ 0 with the dimension less quantity x = M S /T.For illustration, we choose M S = 1000 GeV, M ψ 0 M = 400 GeV and λ S = 0.245.We take y s = 10 −10 and 10 −11 for Fig. 4a and b respectively, to be compatible with the cosmological limitation described above.The comoving abundances of S (Y S ) and ψ 0 (Y ψ 0 ) are denoted by the black dashed and the red solid lines.The orange dotted line corresponds to comoving density leading to observed relic density by PLANCK [4].The green dashed line corresponds to Y ψ 0 in the absence of late-time DM production via non-thermal decay (i.e.y s = 0).We can see from the figures that both S and ψ 0 freeze out from the thermal bath, and after some time, S starts to decay, resulting in additional number densities of ψ 0 .As a result, we see an increase in Y ψ 0 leading to satisfying the observed relic, whereas in the absence of late DM production (green dashed line), it fails to fulfill the observed relic.Comparing Fig. 4a and b, we notice that with decreasing y s the late DM production takes place at a comparably later time (lower T ).However, y s plays no role in deciding the abundance of DM other than determining the lifetime of the S, as S ultimately decays to DM ψ 0 .So, at the later time, we can just assume Y x→∞ To avoid entering in the vicinity of BBN, we take a fixed y s = 10 −10 throughout the analysis.Therefore, with this kind of setup, one can easily accommodate sub-TeV fermion doublet DM with the appropriate choice of M S and λ S .However, despite the fact that the DM ψ 0 can easily fulfill the total observed relic density within the sub-TeV mass range with such minimal setup, the gauge-mediated direct search imposes a strong constraint on DM mass [59].The mass spectrum of and the relevant interaction of direct detection are shown in Fig. 5a and b respectively.The current bound from direct detection experiments by LUX [29], XENON-1T [28] completely rule out the possibility of Dirac DM ψ 0 even up to 2 TeV.Nevertheless, such a scenario with a dark doublet can be revived with pseudo-Dirac DM with the help of an additional scalar triplet as we will discuss in the next section.The elastic direct detection cross-section bound for DM(ψ 0 ) can be easily evaded if the DM turns out to be a pseudo-Dirac state, in which case the Z mediated neutral current vanishes.Such type of scenario can be realized by exploiting the pseudo-Dirac nature of DM where the state ψ 0 splits into two Majorana states.With such motivation, we introduce an additional SU (2) L scalar triplet (with hypercharge, Y = 2) that creates the pseudo-Dirac mass splitting.The relevant parts of the Lagrangian involving are given by, The Yukawa interaction between the vector-like lepton doublet ( ) and the additional scalar triplet( ) is proportional to y ψ and plays a crucial role in generating pseudo-dirac splitting of ψ 0 .However, there are several motivation for introducing a SU (2) L scalar triplet in SM, commonly known as the Type-II seesaw model in literature [13,15,89,90].One of the main motivations for this is to address the explanation of neutrino masses, which have been studied extensively so far [15,89].In this paper, we are only interested in the dark matter analysis.We briefly describe the relevant part of the interaction Lagrangian in this context, L +SM in the Appendix B. The interaction between the dark scalar S and in the modified scenario is described by the Lagrangian where the interaction L S is defined earlier in Eq. ( 4) with v replaced by v d satisfying v 2 d + 2v 2 = 246 GeV.To avoid additional contribution in the thermal abundance of S from , we assume λ S to be zero and M S < M /2 which simplifies our analysis.
The scalar triplet does not acquire any vacuum expectation value (vev).Although electroweak symmetry breaking (EWSB) causes the formation of an induced vev v along the neutral CP even field direction, which can modify the EW parameters.The precision measurements of the EW observ- able constrain the ρ parameters ρ = 1.00038±0.00020[10], putting an upper bound on v as v 2.6 GeV at the 3σ level.There also exists a lower bound on v 10 −9 GeV from lepton flavour violation [91].
Here, the v induces a tiny Majorana mass, y ψ v to the Dirac state ψ 0 thanks to the Yukawa interaction mentioned Eq. ( 12).This Majorana mass term leads to the mass splitting of the ψ 0 state into two physical states, ψ 1 and ψ 2 , after proper diagonalization, as discussed in equation Eq. (C1) 3 .The physical states and their corresponding masses for the lepton doublet are described below: and Here δm is the mass splitting between ψ ± and ψ 0 and is generated due to the quantum correction as previously stated in Eq. (5).The mass splitting between the two pseudo-Dirac states is defined by δm 12 = M ψ 2 − M ψ 1 = 2y ψ v .A typical mass spectrum of the dark fermion doublet is displayed in Fig. 6a for the case y ψ v < δm where The mass spectrum can also be M ψ 2 > M ψ ± > M ψ 1 when y ψ v > δm.In both cases, the lightest neutral pseudo-Dirac state ψ 1 acts as a stable DM candidate.However, for our analysis, we adopt the former hierarchy.With the help of pseudo-Dirac splitting, the Z mediated neutral current interaction turns into and prohibits to write ψ 1 γ μ Z μ ψ 1 interaction (see in Eq. (C3)) as shown in Fig. 6b.The advantage of this setup is that it forbids Z mediated DM-nucleon elastic scattering of ψ 1 and helps to evade the direct detection constraint.We will discuss However, the off-diagonal Z interaction with pseudo-Dirac states,ψ 1,2 in Eq. ( 15) allows for the inelastic DMnucleon scattering, as illustrated in Fig. 7a.In principle, such type of in-elastic scattering kinematically may enable direct search for extremely low mass splitting between two pseudo-Dirac states(ψ 1,2 ), and the condition for such scenario is given as [93], where βc = v DM 650 km/s (escape velocity) and M N is the nucleus mass.For the XENON 1T experiment, Xe is the target nucleus with M N = 130 amu.Following Eq. ( 16), we can derive an upper limit on (δm 12 ) max. as a function of DM mass M ψ 1 .This upper limit sets the threshold below which the Z mediated in-elastic scattering for direct detection is allowed.For DM mass ∼ O(1 TeV), (δm 12 ) max. turns out to be ∼ 250 keV followed by the Eq. ( 16).Therefore the Z-mediated in-elastic scattering of direct search is forbidden with δm 12 250 keV for M ψ 1 ∼ 1 TeV.A detailed discussion on it can be found in the references [93,94].From the definition of δm 12 , we can have a relation between y ψ and v which is depicted by the red dotted line in Fig. 7b.The red shaded region(2y ψ v < (δm 12 ) max. ) below the red dotted line is where inelastic scattering gets allowed.The region above the red dotted line where the Z-mediated inelastic scattering is disallowed.Again the ρ parameter puts an upper bound on v and excluded the region with v > 2.6 GeV, shown by the Gray vertical region.Therefore our point of interest lies within the white region where Z mediated in-elastic scattering is prohibited.
We shall now move to the DM phenomenology of the pseudo-Dirac DM,ψ 1 in the modified scenario with an added scalar triplet .A brief discussion about the interaction Lagrangian of the pseudo-Dirac DM is shown in Appendix C. The thermal abundance of DM ψ 1 is determined by its scalar and gauged mediated annihilation to SM (ψ 1 ψ 1 → SM SM) and co-annihilation to SM (ψ i ψ j ; ψ i ψ ± ; ψ + ψ − → SM SM with i, j = 1, 2).The triplet mediated diagrams depend on the Yukawa coupling y ψ , the masses of the triplet scalars M X , vev of the triplet v and the mixing angle between the CP even scalars sin α.The additional interactions of with the scalar triplet can impact the relic abundance.When the masses of the triplet scalars are smaller than M ψ 1 , the new (co-)annihilation channels, open up and relic density of ψ 1 drops [59].However in our discussion we assume M ψ 1 < M X /2 and small sin α which suppresses the additional interactions due the scalar triplet.Therefore in our analysis we choose a benchmark point(BP) sin α = 10 −4 , M = 3 TeV and v = 0.1 GeV and the masses of the heavy physical states (H 0 , A 0 , H ± , H ++ ) turns out of the order of M (see Eq. (B7)).Hence the additional scalar has no significant role in deciding the relic abundance of DM.
To evaluate the freeze-out abundances of DM (ψ 1 ) with M < M S < M /2 for this setup, we need to solve the coupled Boltzmann equations provided in Eq. (10).And the effective thermal average cross-section, σ v eff in Eq. ( 10) modified as [80,81]: Here . The internal degrees of freedom g ψ 1,2 and g ψ ± are associated with the dark fermion states, ψ 1,2 and ψ ± respectively.Due to the Yukawa coupling of S with ψ 1 , ψ 2 and ψ ± ; the late decay of S will give rise to non-zero abundances of all dark particles ψ 1 , ψ 2 and ψ ± .That ψ ± will eventually decay to the lightest ψ 1 as mentioned in the context of Dirac DM.On the other hand ψ 2 will also decay to ψ 1 promptly due to the strong off diagonal neutral current . In Fig. 8 we plot the comoving abundance of DM, Y ψ 1 with M S /T for the modified scenario.For the plot, we kept fixed M S = 1000 GeV, M = 450 GeV, M = 3 TeV and y S = 10 −10 .We also kept fixed y ψ = 0.1 and v = 0.1 GeV, resulting in pseudo-Dirac mass splitting δm 12 = 20 MeV, consistent with the direct detection requirement.The green dot-dashed line denotes the abundance of DM in the absence of S i.e. no late-time DM production.The cyan dashed line corresponds to the comoving density leading to the correct relic density measured by PLANCK for M = 450 GeV.The blue dashed, red solid and magenta dashed lines correspond to the three different values of λ S = 0.01, 0.23 and 0.5 respectively.From the figure we notice that as the λ S increases, Y ψ ( Y ψ 1 ) decreases.With increase in λ S , the annihilation cross-section of S to SM increases leading to lower freeze out abundance of S, Y x F S .The freeze-out . This property, like the Dirac DM, is easily portrayed in the above-mentioned figure.It is important to note that in the modified scenario, the small pseudo-Dirac mass splitting δm 12 M ψ 1 has no effect on the DM abundance and that DM abundance is almost the same as the Dirac like scenario.As a result, the study of DM relic density under the circumstances of δm 12 M ψ 1 and M ,S < M /2 remains altered for both the scenarios.However, direct detection distinguishes both the Dirac and pseudo-Dirac DM cases as discussed earlier.We now discuss the experimental constraints for the pseudo-Dirac doublet DM scenario.

Direct detection
We discussed in the context of Dirac DM that current bounds from the direct detection experiment rule out Dirac doublet DM for mass even up to 2 TeV.The Majorana nature of the pseudo-Dirac DM easily evades the direct detection constraint.In the modified scenario, the presence of the triplet scalar, which is responsible for pseudo-Dirac splitting, determines the fate of this model in the direct search experiment.Therefore the DM particles can recoil against the nucleus, producing the direct search signature (spin independent) via both tree and loop-level DM-nucleon scattering processes as displayed in Fig. 9.
In Fig. 9 we show the Feynman diagrams for tree-level processes via scalar mixing and loop-level processes via gauge bosons which impact significantly in direct detection crosssection.The loop contributions occur via the gauge mediated diagrams, shown in the Appendix D, which can have a significant contribution to the elastic scattering between DM and detector nucleon.The amplitude for direct detection will be M SI−DD = M tree + M loop .However, for our choice of parameters M = 3 TeV and sin α ∼ O(10 −4 ), the tree level amplitude is mixing and propagator suppressed, while the gauge mediated loop induced amplitude is dominating.A brief discussion on spin-independent direct detection crosssection of the DM ψ 1 followed from [95], can be found in Appendix D.
We found that the loop level contribution(σ SI ∼ 10 −47 cm 2 ) is dominating over the tree level contribution(σ SI ∼ 10 −54 cm 2 ) for M ψ 1 ∼ {100−1000} GeV with the aforementioned benchmark parameters.Therefore the total spinindependent direct detection cross section is well below the existing bounds obtained by the different direct search experiments like XENON-1T [28], Panda-4T [27] and LZ [79].

Indirect detection
Similar to direct detection, the relevant constraint may arise from indirect search experiments like Fermi-LAT [34] and MAGIC [35] by analyzing excess gamma-ray flux.The excess gamma-ray flux can be produced via the production of the SM particles either through DM annihilation or via decay in the local Universe.In this scenario, DM annihilation processes ψ 1 ψ 1 → X X where X = {W − , b, μ − , τ − }; and the subsequent decay X to photons resulted in the production of the gamma ray.Non-observation of DM at the indirect search experiments like Fermi-LAT and MAGIC put an upper bound on the individual thermal averaged annihilation cross-section.For the DM mass M ψ 1 > M W , the most stringent constraint comes from the annihilation process ψ 1 ψ 1 → W − W + .The corresponding Feynman diagram is shown in Fig. 10 which is mediated by ψ ± .
Note the σ v W + W − for pseudo-Dirac DM ψ 1 is smaller than for Dirac DM ψ 0 .The reason for this is the absence of Z mediated s-channel diagram for the process ψ 1 ψ 1 → W + W − in the case of pseudo-Dirac DM.In Fig. 11,

Collider constraint
As stated earlier, due to quantum correction, there exists tiny mass splitting (δm ∼ m π ) between the charged and neutral component of the doublet.For such tiny mass split- ting, the dominate decay mode of the charged fermion is ψ ± → ψ 1 π ± .The corresponding decay width is given by, where m is the mass splitting between ψ ± and ψ 1 , defined as m = δm + δm 12 .Here G F = 1.16638 × 10 −5 , sin θ c = 0.22, f π = 130 MeV and m π = 139.57MeV.If the final decay product is low momentum π ± and stable ψ 1 , it leaves the detector without interaction.Such signature is called the displaced vertex (DV) signatures and LHC has already been constrained such scenario in the context of Higgsino [96,97].
As Higgsino has a similar setup to our doublet fermion DM, we adopt their bound and it turns out for M ψ 1 450 GeV is excluded from the displaced searches at 8 TeV LHC [96].

Results
Finally, we show all the parameter space satisfying the observed relic density via our proposed hybrid setup.In Fig. 12 we do numerical scan with varying M ψ 1 , M S , λ S for a fixed Yukawa coupling y s = 10 −10 .We took a BP with  The color bar resembles the variation in λ S ∈ {0.1 − 1.0}.The grey region is excluded from LHC displaced vertex search.The white region represents the parameter space where the lepton doublet can not be stable observed relic density.In the above mentioned plot we varied M ψ 1 up to 1200 GeV and M S up to 1400 GeV maintaining M S > M ψ 1 GeV so that S can decay on-shell.The color variation in the plot depicts the variation in λ S ∈ {0.1 − 1}.From the figure, we notice that for fixed M S , increase in M ψ 1 calls for an increase in λ S .The reason for such a feature is the fact that a higher value of M ψ 1 leads to a higher Y x F ψ 1 and to meet the observed relic from late-time decay of S, lower value of x F S is required.As Y x F S is inversely proportional to λ S for higher M ψ 1 , the higher value of λ S is needed.This feature is elaborated in the context of Fig. 8.The grey region is excluded from the displaced vertex (DV) signatures [97] which gives the most stringent bound among the experimental constraints.

Conclusion
In this work, we study a minimal setup for lepton doublet dark matter satisfying the observed relic density and direct search constraints in the sub-TeV mass range.The dark matter relic is produced in two different epochs.At first DM abundance is generated thermally from bath particles.At the later epoch, an additional dark sector particle contributes nonthermally to accommodate the total observed DM abundance.The additional dark sector particle was initially in the thermal bath, but it decoupled at the time of non-thermal DM production.This hybrid setup differs from pure WIMP and pure FIMP instances.To realise the scenario, we extend the SM particle content by a SU (2) L lepton doublet and a SM gauge singlet scalar S.An additional Z 2 symmetry has been imposed under which both and S are odd, while all the SM particles are even.In the mass hierarchy M < M S , the neutral component of the lepton doublet, ψ 0 serves as a stable DM candidate.The additional DM production from the late-time non-thermal production can help to compensate for the under-abundance problem in the sub-TeV mass region of .However ψ 0 still faces challenges from direct detection constraints in the sub-TeV mass ranges due to its strong gauge-mediated interactions.To overcome the severe constraint of direct detection (DD), we add a Z 2 even scalar triplet (Y = 2) which helps to generate pseudo-Dirac nature of DM(ψ 1 ), leading to non-diagonal neutral current gauge interaction.Thus using this setup the DM ψ 1 can evade the direct detection constraint.Apart from satisfying the observed DM abundance within the mass range, M ψ 1 1200 GeV, the model has promising detection prospects at the direct, indirect, and collider search experiments.The most stringent limit comes from displaced vertex (DV) signatures [97], which excludes DM masses less than 450 GeV.As a result, the DM mass is consistent with all the constraints within the mass range M ψ 1 ∼ 450−1200 GeV.The DM mass M ψ 1 1200 GeV is ruled out from relic density constraint where DM becomes overabundant.Apart from explaining sub-TeV lepton doublet DM, this minimal framework can also address the non-zero neutrino masses and mixing in presence of the scalar triplet.Thus we present an alternate mechanism where such kind of hybrid set up can provide additional DM abundance.This can help to revive a substantial part of parameter space where the minimal DM models fail to fulfill the observed relic density like in electroweak DM multiplet models.
Acknowledgements The authors thank Dilip Kumar Ghosh for the insightful discussion and helpful suggestions.SJ is supported by CSIR, Government of India, under the NET JRF fellowship scheme with Award file No. 09/080(1172)/2020-EMR-I.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited.[Authors' comment: The experimental data we used in our manuscript for comparison with our estimations are already publicly available.]Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.Funded by SCOAP 3 .SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.and a scalar triplet with hyper-charge Y = 2 reads as [98] where the definition of the co-variant derivative of both the doublet and triplet scalar is defined as, The most general scalar potential including and H can be written as: Note μ 2 > 0 and the does not acquire any vev.However the cubic term in the scalar potential T iσ 2 † leads to the generation of an induce non-vanishing small vev v for the after the EWSB.Then the scalar fields, and can be represented as: In the alignment limit v v d and v = v 2 d + v 2 = 246 GeV.Minimizing the scalar potential at the vacuums (v d and v ) leads to the following conditions . After the EWSB, the two CP even states φ and δ mixed up.The mass matrix can be diagonalized using the orthogonal rotation followed by with the mixing angle The above orthogonal transformation gives rise to two physical states h and H with the physical masses m h ( 125.09 GeV) and m H respectively.Similarly, the CP-odd states are mixed up and lead to one massless Goldstone state eaten by a massive SM Z boson and a massive CP odd eigen state A 0 with mass m A 0 .The orthogonal rotation of the singly charge scalars φ + and + yields one massless Goldstone mode, absorbed by SM W boson and one massive charged eigen state H ± with mass m H ± .The scalar sector also has one massive doubly charged eigen state H ± (≡ ± ) with mass m H ± .The masses of the physical scalars are defined in terms of the couplings, vevs, mixing angle α, and the new mass parameter M as: The Yukawa interaction involving SM lepton and in Eq. (B1), can generate light neutrino masses via Type-II seesaw mechanism [15,89].
for sin α = 10 −4 , the direct detection cross-section is loop dominated.

Fig. 1 Fig. 2
Fig. 1 Cartoon diagrams of relevant interactions active during three different epochs of the Universe to set the relic abundance

Fig. 3
Fig. 3 Variation of τ S with M S − M for different values of y s = 10 −10 , 10 −11 , 10 −12 shown in blue,red and green lines respectively.The solid and dashed line correspond to M S = 1 TeV and 2 TeV respectively.The shaded region corresponds to where the decay will be active after BBN

Fig. 4
Fig. 4 Comoving abundances as a function of inverse temperature for M S = 1000 GeV, M = 400 GeV and λ S = 0.245.We consider y s = 10 −10 in a and y s = 10 −11 in b for demonstration.The black dashed and red solid lines signify Y S and Y ψ 0 respectively.The green dashed dot line correspond to Y ψ 0 in absence of the dark Yukawa coupling i.e. y s = 0.The orange dotted line corresponds to comoving abundances leading to observed relic density for M = 400 GeV

Fig. 5 a
Fig. 5 a Mass spectrum of the Dirac lepton doublet .b Feynman diagram for Z mediated elastic scattering between Dirac DM ψ 0 and nucleon (proton/neutron) N = {p, n}

Fig. 6 a
Fig. 6 a Mass spectrum of in presence of .The mass ordering among charge and neutral states is subject to alert depending on the parameters v and y ψ .b The mass splitting forbids the diagonal neutral current

Fig. 7 a
Fig. 7 a The Feynman diagram for Z mediated in-elastic scattering between ψ 1 and nucleon.b Allowed parameter space in the v vs. y ψ plane.The red shaded region where Z -mediated direct detection is allowed.The gray-shaded vertical region is excluded from ρ parameter.The white region shows the choices of y ψ and v where Z -mediated direct detection limits on lepton doublet DM can be evaded.For our discussion, we chose a benchmark point (BP) which is shown by a blue dot with y ψ = 0.1 and v = 0.1 GeV

Fig. 8
Fig. 8 Comoving abundances of ψ 1 for M S = 1000 GeV, M = 450 GeV with varying λ S ∈ (0.01, 0.23, 0.5) shown by the blue dashed, red solid and magenta dashed lines respectively.The green dashed dot line correspond to Y ψ1 in absence of the dark Yukawa coupling i.e. y s = 0.The thick cyan dotted line corresponds to the comoving abundance leading to observed relic density for M = 450 GeV

Fig. 9
Fig. 9 Elastic scattering of DM ψ 1 with nucleon via a tree level and b loop level process

Fig. 10
Fig. 10 Feynman diagram corresponds to the indirect search mode:ψ 1 ψ 1 → W + W − we plot σ v W + W − as a function of DM mass along with the combined Fermi-LAT and MAGIC exclusion bound.It turns out that the σ v W + W − lies below the bound from indirect detection.Note for the Dirac doublet DM set up the same annihilation channel excluded the DM mass below ∼ 350 GeV.

Fig. 11
Fig. 11 Indirect detection cross-section for DM annihilation to W + W − as a function of DM mass.The blue dashed line corresponds to the pseudo-Dirac DM scenario where DM is ψ 1 and the green solid line represents the Dirac DM where DM is ψ 0 .The combined exclusion bound from Fermi LAT and MAGIC for the DM annihilation channel is shown by cyan region in the same plane

Fig. 12
Fig. 12 Points satisfying observed relic density in M ψ1 vs. M S plane.The color bar resembles the variation in λ S ∈ {0.1 − 1.0}.The grey region is excluded from LHC displaced vertex search.The white region represents the parameter space where the lepton doublet can not be stable 1 GeV and the masses of the heavy scalars (m H 2 ,A 0 ,H + ,H ++ ∼ 3 TeV) such that they have no effect on DM abundance.With such parameter choices, the scalar triplet sector has no such role in deciding Y ψ 1 (x FO ψ 1 ) apart from evading the stringent direct detection bound.Similar to Dirac DM scenario for M ψ 1 < 1200 GeV, Y x→∞ ψ 1 is decided by the gauge interactions and the late time decay of S repopulate Y ψ 1 leading to