Gravitational wave effects and phenomenology of a two-component dark matter model

We study an extension of the Standard Model (SM) which could have two candidates for dark matter (DM) including a Dirac fermion and a vector dark matter (VDM) under a new $U(1)$ gauge group in the hidden sector. The model is classically scale-invariant and the electroweak symmetry breaks because of loop effects. We investigate the parameter space allowed by current experimental constraints and phenomenological bounds. We probe the parameter space of the model in the mass range $1<M_V<5000$ GeV and $1<M_{\psi}<5000$ GeV. It has been shown that there are many points in this mass range that are in agreement with all phenomenological constraints. The electroweak phase transition has been discussed and it has been shown that there is region in the parameter space of the model consistent with DM relic density and direct detection constraints that, at the same time, can lead to first order electroweak phase transition. The gravitational waves produced during the phase transition could be probed by future space-based interferometers such as LISA and BBO.


Introduction
The Standard Model (SM) has been the most effective way to describe the functioning of the world around us, but is incomplete because of challenges like matter-antimatter asymmetry, hierarchy problem, and dark matter (DM).DM is estimated to make up approximately 27 % of our universe, as indicated by a lot of astrophysical and cosmological evidence [1].One of the main goals of particle physicists is to predict and find a particle that can satisfy the properties of DM, which can be a window to physics beyond the standard model (BSM).
Weakly interacting massive particles (WIMPs) are the most popular candidate for DM, with the freeze-out scenario being the most popular choice [2].WIMP paradigm is essential background for almost any discussion of particle DM and the triple coincidence of motivations from particle theory, particle experiment, and cosmology is known as the WIMP miracle.However, no trace of DM has so far been found in direct detection experiments.Due to strong constraints on direct detection experiments in one-component DM models, multi-component DM models seem more appropriate in some ways .
In the SM, electroweak phase transition is of the second order [54,55] and does not generate the gravitational wave (GW) signal (for a recent review see [56]).A first-order phase transition can be caused by certain extensions of the SM and the DM candidate, leading to the creation of GWs .In the early universe when two local minima of free energy (potential) co-exist for some range of temperatures (critical temperature), strongly first-order electroweak phase transition can take place.After that, the relevant scalar fields can quantum-mechanically tunnel into the new phase and through the nucleation of bubbles and collide with each other to cause a significant background of GWs [86][87][88][89][90].The discovery of GWs resulting from the first-order phase transition can be the consequence of physics BSM, which can be a supplement to ground experiments like those conducted using the Large Hadron Collider (LHC).Unlike GWs from strong astrophysical sources [91], these waves have a range between millihertz and decihertz [92].
The Laser Interferometer Space Antenna (LISA) [93] and Big Bang Observer (BBO) [94] are two space-based GW detectors which are expected to observe GWs resulting from cosmological phase transitions in future years.On the other hand, one of the Sakharov conditions [95] which explains the matter-antimatter asymmetry in universe is the thermal imbalance that occurs in first-order phase transitions.
As mentioned, one of the fundamental challenges of particle physics is the hierarchy problem.
A potential solution to this problem is to drop the Higgs mass term in the potential.SM without the Higgs mass term is scale-invariant.In this paper, we present a classically scale-invariant extension of the SM where all the particle masses are generated using the Coleman-Weinberg mechanism [96].The model includes three new fields: a fermion, a complex singlet scalar and a vector field with U D (1) gauge symmetry.We examine two scenarios in this paper.In scenario A, we consider only the fermionic field as DM.In scenario B, both fermionic and vector fields are considered as DM.We probe the parameter space of the model according to constraints from relic density and direct detection.DM relic density is reported by the Planck Collaboration [97] and DM-Nucleon cross section is constrained by XENONnT experiment results [98].We investigate the possibility of the electroweak phase transition with respect to the bounded parameter space, where we use the effects of the effective potential of the finite temperature.We probe the parameter space of the model which is consistent with the said phenomenological constraints and also leads to a strong first-order electroweak phase transition.The GW signal resulting from this phase transition also has been studied in the LISA and BBO detectors.
Here is the organization of the paper.In the next section, we introduce the model.In Sect.3, we study the phenomenology of the Scenario A including relic density, direct detection, invisible Higgs decay and the resulting GWs.Section 4 is dedicated to the phenomenology of the Scenario B and its GW signals.Finally, our conclusion comes in Sect. 5.

The model
In this section, we consider an extension of the SM to explain DM phenomenology.In this regard, the model contains three new fields in which a vector dark matter V µ and a Dirac fermion field ψ can play the role of DM.A complex scalar,S, mediates between SM and the dark Field Table 1: The charges of the dark sector particles under the new U (1) D symmetry.
sector.In the model V µ , ψ and S are charged under a new dark U (1) D gauge group.All of these fields are singlet under SM gauge groups.We suppose the mass of the fermion was produced by breaking of dark U D (1) gauge symmetry, and so was constrained by other parameters of the model.However, in the model, the dark sector is invariant under the transformations of the U D (1) gauge group: where the U (1) D charge of the new particles, Q l,r,s , are given in Table 1.
The Lagrangian for the model is given by the following renormalizable interactions, where L SM is the SM Lagrangian without the Higgs potential term, The covariant derivative is , and On the imposition of dark charge conjugation symmetry, we do not assume a kinetic mixing term between V µ and the U (1) Y gauge boson of the SM.
The most general scale-invariant potential V (H, S) that is renormalizable and invariant under gauge symmetry is Note that the quartic portal interaction, λ SH (S * S)(H † H), is the only connection between the dark sector and the SM.
SM Higgs field H as well as dark scalar S can receive VEVs, respectively breaking the electroweak and U D (1) symmetries.In unitary gauge, the imaginary component of S can be absorbed as the longitudinal part of V µ .In this gauge, we can write where h 1 and h 2 are real scalar fields which can receive VEVs.Now, the tree-level potential becomes There is a Z 2 symmetry for ψ, making it a stable particle.In addition, if the mass of V µ is less than two times of the mass of ψ, then both V µ and ψ are viable DM candidates.
For the Hessian matrix, we define: Necessary and sufficient conditions for local minimum of V tree in which vacuum expectation values h 1 = ν 1 and h 2 = ν 2 , have be written as: where det(H(ν 1 , ν 2 )) is determinant of the Hessian matrix.Condition (2.8) for non-vanishing VEVs leads to the following constraints Conditions (2.8) and (2.9) require λ H > 0, λ S > 0, and λ SH < 0. However, condition (2.10) will not be satisfied, because det(H(ν 1 , ν 2 )) = 0.When the determinant of the Hessian matrix is zero, the second derivative test is inconclusive, and the point (ν 1 , ν 2 ) could be any of a minimum, maximum or saddle point.However, in the model, constraint (2.11) is defined as a flat direction, in which V tree = 0. Thus, it is the stationary line or a local minimum.
The important point is that for other directions V > 0, and the tree level potential only vanishes along the flat direction.Thus, the full potential of the theory will be dominated by higher-loop contributions along flat direction and specifically by the one-loop effective potential.Considering one-loop effective potential, , can lead to a small curvature in the flat direction which picks out a specific value along the ray as the minimum with V 1−loop ef f < 0 and vacuum expectation value ν 2 = ν 2 1 + ν 2 2 characterized by a RG scale Λ.Since at the minimum of the one-loop effective potential V tree 0 and along the flat direction (where V tree = 0) is a global minimum of the full potential, and so spontaneous symmetry breaking takes place.As a result, we suppose h 1 → ν 1 + h 1 and h 2 → ν 2 + h 2 and the electroweak symmetry breaks with value ν 1 = 246 GeV.In tree level potential since h 1 and h 2 mix with each other, they can be rewritten by the mass eigenstates H 1 and H 2 as where H 2 is along the flat direction, thus M H 2 = 0, and H 1 is perpendicular to the flat direction which we identify it as the SM-like Higgs observed at the LHC with M H 1 = 125 GeV.After the where M ψ and M V are the masses of vector and fermion fields after symmetry breaking.Conditions (2.13) constrain free parameters of the model up to three independent parameters.We choose M V , M ψ and g v as the independent parameters of the model.
Since in tree level, M H 2 = 0, and the elastic scattering cross section of DM off nuclei becomes severely large, the model actually is excluded by direct detection experiments.However, the radiative corrections give a mass to the massless eigenstate H 2 .One-loop corrections to the potential, via the Gildener-Weinberg formalism [99], shift the scalon mass to the values that can be even higher than SM Higgs mass.Along the flat direction, the one-loop effective potential where a and b are the dimensionless constants that given by k . (2.15) In (2.15), M k and g k are the tree-level mass and the internal degrees of freedom of the particle k, respectively (In our convention g k takes positive values for bosons and negative values for fermions).
Minimizing (2.14) shows that the potential has a non-trivial stationary point at a value of the RG scale Λ, given by Eq. (2.16) can now be used to find the form of the one-loop effective potential along the flat direction in terms of the one-loop VEV ν It is remarkable that the scalon does not remain massless beyond the tree approximation.Re- , M H 2 will be Considering (2.15), the scalon mass can be expressed in terms of other particle masses where M W,Z,t are the masses of W, Z gauge bosons, and top quark, respectively.As mentioned before, M H 1 = 125 GeV and be a minimum, it must be less than the value of the potential at the origin, hence it must be negative.From (2.19), we have a constraint on the parameter space of our model where M H 2 > 0.
Note that according to (2.19) and (2.13), M H 2 is completely determined by the independent parameters of the model, i.e., M V , M ψ and the coupling g v .These constraints are due to the scale invariance conditions which were imposed on the model.Depending on the new particle masses of our model, we examine two different scenarios.In scenario A, the ψ field is considered as DM.In scenario B, M V < 2M ψ and both V µ and ψ fields are considered as DM.In the following, we examine the phenomenology of each scenario separately.
3 Scenario A 3.1 DM phenomenology

Relic density
In this case, since we do not assume M V < 2M ψ , the only candidate for DM is ψ.In the WIMP scenario, at first the early universe is hot and very dense and all particles are in thermal equilibrium.Then the universe cools until its temperature falls below the mass of DM particles, and the amount of DM becomes Boltzmann suppressed, dropping exponentially as e −m X /T .As the universe expands, the DM particles are diluted and can no longer find each other until they are annihilated and are out of thermal equilibrium with the SM particles.Then the DM particles freeze out and their number asymptotically reaches a constant value as their thermal relic density.The evolution of the number density of DM particle (n ψ ) with time is governed by the Boltzmann equation: where H is the Hubble parameter and n eq ψ ∼ (m ψ T ) 3/2 e −m ψ /T is the particle density before particles get out of equilibrium.The relevant Feynman diagrams for DM production are shown in the Fig. 1.We calculate the relic density numerically for the ψ particle by implementing the model into micrOMEGAs [100].Figure 2 shows the parameter space of the model in agreement with the observed density relic [97].As can be seen, there is agreement for 400 < M V < 5000 GeV , 20 < M ψ < 2500 GeV and 0.1 < g v < 6.

Direct detection
WIMPs may be detected by the scattering off normal matter through processes XSM →XSM .
Given a typical WIMP mass of m X ∼ 100 GeV and WIMP velocity υ ∼ 10 −3 , the deposited recoil energy is limited to ∼ 100 keV, so detection requires highly-sensitive, low-background and deep-site detectors.Such detectors are insensitive to very strongly-interacting DM, which would   be stopped in the atmosphere or earth and would be undetectable underground.The spinindependent direct detection(DD) cross sections of ψ were obtained using the micrOMEGAs package [100].
In Fig 3, the parameter space of the model is drawn in agreement with the limits of relic density, XENONnT and neutrino floor.For 20 < M ψ < 1000 GeV, there will be points of the parameter space that fall below the XENONnT limit.As can be seen from Fig. 3, for M ψ < M H 1 /2 there will be points that allow the investigation of the invisibility Higgs decay, which will be investigated in the next section.

Invisible Higgs decay
As mentioned, for a parameter space consistent with ψ relic density and DD, SM Higgs(H 1 ) can only kinematically decay into a pair of ψ.Therefore, H 1 can contribute to the invisible decay mode with a branching ratio: where Γ(h) SM = 4.15 [MeV] is total width of Higgs boson [101].The partial width for process H 1 → ψψ is given by: The SM prediction for the branching ratio of the Higgs boson decaying to invisible particles which coming from process h → ZZ * → 4ν [102], [103], [104], [105] is, 1.2 × 10 −3 .CMS Collaboration has reported the observed (expected) upper limit on the invisible branching fraction of the Higgs boson to be 0.18(0.10)at the 95% confidence level, by assuming the SM production cross section [106].A similar analysis was performed by ATLAS collaboration in which an observed : The cross points depict allowed region that is consistent with invisible Higgs decay at [107].
upper limit of 0.145 is placed on the branching fraction of its decay into invisible particles at a 95% confidence level [107].
Figure 4, shows the allowed range of parameters by considering ATLAS [107] upper limit for invisible Higgs mode.

Finite temperature potential
In addition to the 1-loop zero-temperature potential (2.17), we can also consider the 1-loop corrections at finite temperature in the effective potential, which is with thermal functions The above functions can be expanded in terms of modified Bessel functions of the second kind, K 2 (x) [71], The contribution of resummed daisy graphs is also as follows [109] V daisy (H 2 , T ) = where the sum runs only over scalar bosons and longitudinal degrees of freedom of the gauge bosons.Thermal masses, Π k (T ), are given by Finally, the one-loop effective potential including both one-loop zero temperature 2.17 and finite temperature 3.4 and 3.7 corrections is given by In order to get V ef f (0, T ) = 0 at all temperatures, We make the following substitution: Now we are ready to study the phase transition and the resulting gravitational waves.

Gravitational waves
The characteristic of first-order phase transitions is the existence of a barrier between the symmetric and broken phases.The electroweak phase transition takes place after the temperature of the universe drops below the critical temperature(T C ).At this temperature, effective potential(3.9)has two degenerate minimums, one in H 2 = 0 and the other in H 2 = ν C = 0 : By solving these two equations, one can obtain ν C and T C .If this phase transition is strongly first order, it can satisfy the condition of departure from thermal equilibrium, which is one of Sakharov conditions for creating baryonic asymmetry in the universe.There is a criteria for strongly electroweak phase transition [95,110] which is as follows (3.12) The transition from the false to the true vacuum proceeds via thermal tunneling at finite temperature.This concept can be grasped in the context of formation of bubbles of the broken phase in the sea of the symmetric phase.Once this has happened, the bubble spreads throughout the universe, converting false vacuums into true ones.Bubbles formation starts at the nucleation temperature T N , where one can estimate T N by the condition S 3 (T N )/T N ∼ 140 [111].The function S 3 (T ) is the three-dimensional Euclidean action for a spherical symmetric bubble given by where H 2 satisfies the differential equation which minimizes S 3 : with the boundary conditions: In order to solve Eq. 3.14 and find the Euclidean action 3.13, we used the AnyBubble package [112].In the following, we will show that the nucleation temperature(T N ) will be much lower than the critical temperature (T C ), indicating a very strong phase transition.
GWs resulting from the strong first-order electroweak phase transitions are have three causes, which are as follows: • collisions of bubble walls and shocks in the plasma, • sound waves to the stochastic background after collision of bubbles but before expansion has dissipated the kinetic energy in the plasmal, • turbulence forming after bubble collisions.
These three processes may coexist, and each one contributes to the stochastic GW background: There are four thermal parameters that control the above contributions: • T N : the nucleation temperature, • α: the ratio of the free energy density difference between the true and false vacuum and the total energy density, where ρ * is • β: the inverse time duration of the phase transition, • υ ω : the velocity of the bubble wall which is anticipated to be close to 1 for the strong transitions [113].
Isolated spherical bubbles cannot be used as a source of GWs, and these waves arise during the collision of the bubbles.The collision contribution to the spectrum is given by [114] Ω coll (f )h 2 = 1.67 where S coll parametrises the spectral shape and is given by where Hz. (3.22) The collision of bubbles produces a massive movement in the fluid in the form of sound waves that generate GWs.This is the dominant contribution to the GW signal, and is given by [115] Ω The spectral shape of S sw is where Hz. (3.25) Plasma turbulence can also be caused by bubble collisions, which is a contributing factor to the GW spectrum and is given by [116] Ω where In computing GW spectrum we have used [117,118] where the parameters κ, κ υ , and κ turb denote the fraction of latent heat that is transformed into gradient energy of the Higgs-like field, bulk motion of the fluid, and MHD turbulence, respectively.
To investigate the GWs resulting from the first-order electroweak phase transition, we choose three benchmark points.These points are presented in table 2. Figure 5 shows the potential behavior for both critical and nucleation temperatures.Also, S 3 /T changes in terms of temperature are shown.In table 2 all relevant quantities, including independent parameters of the model, DM properties, and phase transition parameters, are given.The benchmarks 1 and 2 are consistent with direct detection constraint while the benchmark 3 is outside the range of XENONnT and is placed under the neutrino floor.The GW spectrum for these benchmark points is depicted in figure 6.The GW spectrum for these benchmarks 1,(2,3) falls within the observational window of LISA(BBO).Therefore, for benchmark 3, GW can be a special way to probe it.
), ( ), ( where j runs over SM massive particles, H 1 and H 2 .In σ ab→cd υ all annihilations are taken into account except σ ψV →ψV υ which does not affect the number density.In the above relations, for simplicity in the writing, only the annihilation and conversion contributions are shown in the equations.But in practice, all contributions, even semi-annihilations are included in the micrOMEGAs package to solve these equations [119]. )], (4.3) )], (   2.
where g is the degrees of freedom parameter and M pl is the Planck mass.It is clear from the above equations that there are new terms in the Boltzmann equations that describe the conversion of two DM particles into each other.Because these two cross sections are also described by the same matrix element, we expect σ V V →ψψ υ and σ ψψ→V V υ are not independent and their relation is: The interactions between the two DM components take place by exchanging two scalar mass eigenstates H 1 and H 2 where the coupling of V to H 1 is suppressed by sin α.For this reason, it usually is the H 2 -mediated diagram that gives the dominant contribution.We also know that the conversion of the heavier particle into the lighter one is relevant.The relic density for any DM candidate associated with the Y at the present temperature is given by the following relation: where h is the Hubble expansion rate at present time in units of 100(km/s)/M pc.We used the micrOMEGAs package [119] to numerically solve coupled Boltzmann differential equations.
According to the data from the Planck Collaboration [97], the DM constraint in this model reads We also define the fraction of the DM density of each component by, In Fig. 8, the parameter space consistent with DM relic density is obtained.As can be seen, there is an agreement with the relic density observed for 300 < M V < 5000 GeV, 200 < M ψ < 5000 GeV and 0.1 < g v < 6.Of course, it is necessary to mention that the contribution of semi-annihilations in the model is important.For example, for benchmark 2 in Table 3, using the micrOMEGAs package, we found that the share of the cross section of the process of V ψ → ψ H 2 includes 24 % of the cross section of all processes.

Direct detection
We investigate constraints on parameters space of the model which are imposed by searching for scattering of DM-nuclei.The spin-independent direct detection(DD) cross sections of V and ψ are determined by H 1 and H 2 exchanged diagrams [30,70]: where M N is the nucleon mass and f N ≃ 0.3 parameterizes the Higgs-nucleon coupling.
Various DD experiments have placed constraints on DM-Nucleon spin independent cross section, such as LUX [120], PandaX-II [121], XENON1T [122],LZ [123] and XENONnT [98].Of course, these experiments are gradually approaching what is called the neutrino floor [124], which is a the irreducible background coming from scattering of SM neutrinos on nucleons.We use the XENONnT [98] experiment results to constrain the parameter space of our model.In this experiment there is a minimum upper limit on the spin-independent WIMP-nucleon cross section of 2.58 × 10 −47 cm 2 for a WIMP mass of 28 GeV .In order to study the effect of the direct detection experiment on the model, we use rescaled DM-Nucleon cross section ξ V σ V and ξ ψ σ ψ .
In Fig.
Consequently, the effective 5(6)-dimensional interaction terms for DM-quark interactions at low energies will be AB sin α cos α tive coupling between DM and quarks approaches zero, resulting in a dip in the DM-nucleon cross-section.In [125], a study has been done on degenerate Higgs scenario.In [126], it is shown by using the high resolution of the diphoton channel of the Higgs boson decays, the mass difference between the two degenerate states ∆m 3 GeV is disfavored at the 2σ level from the LHC Run-I data.In order to test a degenerate Higgs scenario, a possible proposal is consider- shown [125] that for ∆m 1 GeV, it is possible to distinguish between two Higgs.

Invisible Higgs decay
In this scenario, according to Fig. 9, because there is no point where M ψ,V,H 2 < M H1 /2, it is not necessary to check the invisible Higgs decay.

Electroweak phase transition and gravitational waves
To investigate the phase transition and the resulting GWs, we follow the procedure of Sec.3.2.
We select three benchmark points as shown in Table 3. Benchmarks 1 and 2 are consistent with direct detection constraint, while benchmark 3 is placed under the neutrino floor.Figure 10 shows the changes of potential and S 3 /T in terms of temperature.The GW spectrum for these benchmark points is depicted in Fig. 11.The GW spectrum for these benchmarks falls within the observational window of BBO.For benchmark 1, gravitational wave peak will also be in LISA's range.For the benchmark 3, which is below the neutrino floor limit, future gravitational wave detectors can provide a special way to discover and investigate.

Conclusion
We have considered an extension of the SM with a new U (1) symmetry in the dark part.3.
collider constraints (invisible Higgs decay in scenario A).A three-dimensional parameter space acquisition was completed.
We focused our attention on the phase transition dynamics after presenting the model and exploring DM phenomenology.The full finite-temperature effective potential of the model at the one-loop level was obtained to investigate the nature and strength of the electroweak phase transition, with the aim of exploring its nature and strength.A first-order electroweak phase transition can exist when there is a barrier between the broken and symmetric phases.It was demonstrated that the finite-temperature effects induce such a barrier and thereby give rise to a phase transition which can generate GWs.
After studying the phase transition, we investigated the resulting GWs.The parameters required to investigate GWs can all be calculated from our presented model and are a function of the independent parameters of the model.We have demonstrated that the model can survive

Figure 1 :
Figure 1: The relevant Feynman diagrams for DM relic density production cross section.

Figure 2 :
Figure 2: The allowed range of parameter space consistent with DM relic density.

Figure 3 :
Figure 3: The allowed range of parameter space consistent with DM relic density and DD.

Figure 4
Figure4: The cross points depict allowed region that is consistent with invisible Higgs decay at[107].
3.27, h * is the value of the inverse Hubble time at GW production, redshifted to today, h * = 1.65 × 10 −5 T N 100 In this scenario, M V < 2M ψ and both V µ and ψ fields are considered as DM.The evolution of the number density of DM particles with time are governed by the Boltzmann equation.The coupled Boltzmann equations for fermion ψ and vector DM are given by: The relevant Feynman diagrams for DM production are shown in Fig 7.By choosing x = M/T and Y = n/s, , where T and s are the photon temperature and the entropy density, respectively, one can rewrite the Boltzmann equations in terms of Y = n/s:

Figure 5 :Figure 6 :
Figure 5: In (a), (b) and (c) Potential behavior are given for critical temperature and nucleation temperature.In (d), (e) and (f) S 3 /T changes in terms of temperature are also given for all three benchmarks.

Figure 7 :
Figure 7: The relevant Feynman diagrams for DM relic density cross section including: a(annihilation), b(conversion) and c(semi-annihilation).

Figure 8 :
Figure 8: The allowed range of parameter space consistent with DM relic density.

Figure 9 :
Figure 9: The allowed range of parameter space consistent with DM relic density and DD.In (a) ξ V σ V VS M V and in (b) ξ ψ σ ψ VS M ψ has shown.

Figure 10 :
Figure 10: In (a), (b) and (c) Potential behavior are given for critical temperature and nucleation temperature.In (d), (e) and (f) S 3 /T changes in terms of temperature are also given for all three benchmarks.

Figure 11 :
Figure 11: GW spectrum for benchmark points of the Table3.

Table 2 :
Three benchmark points with DM and phase transition parameters.
9, rescaled DM-Nucleon cross sections(ξ V σ V and ξ ψ σ ψ ) are depicted for the parameters that are in agreement with the relic density.It is clear from the figure that there are some points between the XENONnT direct detection bound and the neutrino floor which can be probed in future direct detection experiments.As is expected from 4.10, by reducing the mass difference between H 1 and H 2 , the cross section decreases and therefore allowable points increase.In our model, DM interacts with nucleons through H 1 and H 2 mediators.The relevant terms in the Lagrangian are Ah 1 [qq] and Bh 2 [V µ V µ ( ψψ)] for vector (spinor) DM, where A and B are some constants.Therefore, theH 1 mediator involves A cos αH 1 [qq] − B sin αH 1 [V µ V µ ( ψψ)].Similarly, for the H 2 mediator, the terms are A sin αH 2 V h 2 Ω ψ h 2 Ω DM h 2 ξ V σ V (cm 2 ) ξ ψ σ ψ (cm 2 ) 1 1.17

Table 3 :
Three benchmark points with DM and phase transition parameters.