Gravitational waves from first-order phase transition in an electroweakly interacting vector dark matter model

Many astrophysical observations show the existence of dark matter (DM). DM constitutes approximately 26% of the energy in the universe Planck:2018vyg . However, the nature of DM remains unclear. Models beyond the standard model (SM) of particle physics often predict DM candidates. In many particle DM models, the freeze-out mechanism Lee:1977ua is used to explain the measured value of the DM energy density. The mechanism requires a pair of DM particles to annihilate into other particles in the thermal bath in the early Universe. The canonical value of the annihilation cross section, which can explain the measured value of the DM energy density, is $\expectationvalue{\sigma v}\simeq 3\times 10^{-26}$ cm³ s${}^{-1}\simeq 1$ pb $c$. This value is of the same order as the cross section of processes by the electroweak interaction and implies that DM particles interact with the SM particles via the electroweak interaction. An example of such DM is the wino DM, which is an SU(2)_L triplet spin-$\frac{1}{2}$ Majorana fermion, and the mass prediction of the thermally produced wino DM is approximately 3 TeV Hisano:2006nn ; Cirelli:2007xd .

The electroweakly interacting vector DM model proposed in Abe:2020mph is one of the attractive DM models. In the model, the electroweak gauge symmetry SU(2)${}_{L}\times$ U(1)_Y is extended into SU(2)${}_{0}\times$SU(2)${}_{1}\times$SU(2)${}_{2}\times$U(1)_Y. The SU(2)${}_{0}\times$SU(2)${}_{1}\times$SU(2)₂ gauge symmetry is spontaneously broken into SU(2)_L by vacuum expectation values of two scalar fields $\Phi_{1}$ and $\Phi_{2}$. An exchange symmetry SU(2)${}_{0}\leftrightarrow$SU(2)₂ is imposed, and it stabilizes linear combinations of gauge fields $V^{a}_{\mu}\equiv(W_{0\mu}^{a}-W_{2\mu}^{a})/\sqrt{2}$, where $W_{j\mu}^{a}$ is the gauge field of SU(2)_j with $j=0,1,2$. After the electroweak symmetry breaking SU(2)${}_{L}\times$U(1)${}_{Y}\to$U(1)${}_{\text{em}}$, quantum corrections make charged component of $V^{a}$, which is denoted as $V^{\pm}$, slightly heavier than the neutral component $V^{0}$. Consequently, $V^{0}$ is a DM candidate. Note that $V^{a}$ arises from the SU(2)₀ and SU(2)₂ gauge fields, and the SU(2)_L gauge field are linear combination of the SU(2)₀, SU(2)₁, and SU(2)₂ gauge fields. Hence, $V^{a}$ directly couples to the SU(2)_L gauge bosons and have the electroweak interaction. Therefore, the model predicts electroweakly interacting vector DM.

Through the electroweak interaction, $V^{a}$ is in equilibrium with the SM particles in the early universe, and the DM abundance is determined by the freeze-out mechanism. In addition to $V^{a}$ and the SM gauge bosons, the model predicts extra heavy gauge bosons, denoted by $W^{\prime\pm}$ and $Z^{\prime}$. They are approximately SU(2)_L triplet and have common mass, $m_{W^{\prime}}\simeq m_{Z^{\prime}}$. Annihilation of a pair of $V^{0}$ depend on $m_{Z^{\prime}}$. For example, if $m_{Z^{\prime}}\lesssim 2m_{V}-m_{W}$, $V^{0}V^{0}\to W^{\prime-}W^{+}$ is open, but otherwise it is kinematically forbidden. Hence, the value of $m_{Z^{\prime}}$ affects the DM relic abundance. The mass of the vector DM that reproduces the measured value of the DM energy density varies from $3$ TeV to $19$ TeV Abe:2020mph .

Other than relic abundance, direct detection of $V^{0}$, $W^{\prime}/Z^{\prime}$ search in collider experiments, and indirect detection of DM were studied in Abe:2020mph ; Abe:2021mry . However, these observables are insufficient to probe all the regions of the model parameter space. Other observables are necessary to test the model comprehensively.

In this work, we focus on gravitational waves (GWs), aiming to extend the region of the parameter space where we can test the electroweakly interacting DM model with future/current experiments. Because the gauge couplings $g_{0}$ and $g_{1}$ are $\sim 1$ to explain the relic abundance Abe:2020mph , the scalar potential is modified by the gauge boson contributions at the loop level. As a result, the scalar potential can experience a first-order phase transition in the early universe. It is well-known that first-order phase transitions generate GWs Grojean:2006bp . To observe the GW spectra, future space-based GW interferometers, such as LISA LISA , DECIGO DECIGO , and BBO BBO , can be used; thus, DM models with a first-order phase transition can be tested through GW observational experiments Schwaller:2015tja ; Chala:2016ykx ; Baldes:2017rcu ; Chao:2017vrq ; Beniwal:2017eik ; Addazi:2017gpt ; Tsumura:2017knk ; Huang:2017rzf ; Huang:2017kzu ; Hektor:2018esx ; Hashino:2018zsi ; Baldes:2018emh ; Madge:2018gfl ; Beniwal:2018hyi ; Bian:2018mkl ; Bai:2018dxf ; Bian:2018bxr ; Shajiee:2018jdq ; Mohamadnejad:2019vzg ; Bertone:2019irm ; Kannike:2019mzk ; Paul:2019pgt ; Croon:2019rqu ; Hall:2019ank ; Chen:2019ebq ; Hall:2019rld ; Barman:2019oda ; Chiang:2019oms ; Borah:2020wut ; Kang:2020jeg ; Pandey:2020hoq ; Hong:2020est ; Alanne:2020jwx ; Bhoonah:2020oov ; Han:2020ekm ; Wang:2020wrk ; Ghosh:2020ipy ; Huang:2020crf ; Deng:2020dnf ; Chao:2020adk ; Azatov:2021ifm ; Zhang:2021alu ; Davoudiasl:2021ijv ; Reichert:2021cvs ; Mohamadnejad:2021tke ; Bian:2021dmp ; Costa:2022oaa ; Liu:2022jdq ; Shibuya:2022xkj ; Costa:2022lpy ; Kierkla:2022odc ; Morgante:2022zvc ; Chakrabarty:2022yzp ; Arcadi:2022lpp ; Frandsen:2022klh . We investigate the GWs generated by the phase transition from SU(2)${}_{0}\times$SU(2)${}_{1}\times$SU(2)₂ to SU(2)_L in the electroweakly interacting vector DM model. As discussed below, it allows us to explore the region of the parameter space where other observables, such as the $W^{\prime}$ search, cannot probe.

The rest of this paper is organized as follows. In Section 2, we briefly introduce the vector DM model with ${\rm SU}(2)_{0}\times{\rm SU}(2)_{1}\times{\rm SU}(2)_{2}\times{\rm U}(1)_{Y}$ gauge symmetry. In Section LABEL:sec:potential, we show the effective potential with finite-temperature effects for this model. We clarify the parameter region with a first-order phase transition, which can produce a detectable GW spectrum. The formula for the GW spectrum from the first-order phase transition is presented in Section LABEL:sec:GW. We discuss in Section LABEL:sec:numerical the testability of the model at the GW observation experiments, such as the LISA, DECIGO, and BBO experiments. Finally, Section 6 presents the conclusions of this study. In appendix A, we give the explict expression of the effective potential at the zero temperature.

In this section, we introduce the electroweakly interacting vector DM model proposed in Abe:2020mph .

The model exhibits gauge symmetry, described as $\text{SU(3)}_{c}\times\text{SU(2)}_{0}\times\text{SU(2)}_{1}\times\text{SU(2)}_{2}\times\text{U(1)}_{Y}$. Here, $\text{SU(3)}_{c}$ and $\text{U(1)}_{Y}$ correspond to the gauge symmetries governing the quantum chromodynamics (QCD) and hypercharge, respectively. Because the QCD sector is the same as that in the SM, we focus on the electroweak sector, denoted as $\text{SU(2)}_{0}\times\text{SU(2)}_{1}\times\text{SU(2)}_{2}\times\text{U(1)}_{Y}$. We use the notation $W^{a}_{j\mu}$ and $B_{\mu}$ for gauge bosons associated with $\text{SU(2)}_{j}$ and $\text{U(1)}_{Y}$, respectively. Here, $j$ can take values of 0, 1, or 2, and $a$ can take values of 1, 2, or 3. $g_{j}$ and $g^{\prime}$ are the gauge couplings for $\text{SU(2)}_{j}$ and $\text{U(1)}_{Y}$, respectively.

We introduce two scalar fields, $\Phi_{1}$ and $\Phi_{2}$, expressed in two-by-two matrices. They transform under gauge transformation as

where $U_{0},U_{1}$, and $U_{2}$ represent two-by-two unitary matrices for $\text{SU(2)}_{0},\text{SU(2)}_{1}$, and $\text{SU(2)}_{2}$, respectively. In addition, we impose the following conditions for $\Phi_{j}$ to reduce their degrees of freedom,

Hence, each $\Phi_{j}$ consists of four real scalar fields.

All other fields remain identical to those in the SM, except that they are charged under $\text{SU(2)}_{1}$ instead of $\text{SU(2)}_{L}$. The charge assignments for the matter fields are summarized in Table 1.

In addition to gauge symmetry, this model exhibits exchange symmetry. The Lagrangian is invariant under the following field transformations:

whereas all other fields remain unchanged. This symmetry is equivalent to the exchange between $\text{SU(2)}_{0}$ and $\text{SU(2)}_{2}$, implying that the gauge couplings of $\text{SU(2)}_{0}$ and $\text{SU(2)}_{2}$ must be identical. It is important to note that under this symmetry, $\frac{W^{a}_{0\mu}-W^{a}_{2\mu}}{\sqrt{2}}$ and $\frac{\Phi_{1}-\Phi_{2}}{\sqrt{2}}$ change sign, whereas $\frac{W^{a}_{0\mu}+W^{a}_{2\mu}}{\sqrt{2}}$, $\frac{\Phi_{1}+\Phi_{2}}{\sqrt{2}}$, and the other fields remain unchanged. Therefore, the symmetry described in Eq. (3) is equivalent to the $Z_{2}$ symmetry commonly used in DM models. The lightest particle among $\frac{W^{a}_{0\mu}-W^{a}_{2\mu}}{\sqrt{2}}$ and $\frac{\Phi_{1}-\Phi_{2}}{\sqrt{2}}$ is stable and is a dark matter candidate for this model.

The Lagrangian of the scalar and the electroweak gauge sectors are described as

where

Some couplings are equal owing to the exchange symmetry described in Eq. (3).

We assume that the scalar fields develop the following vacuum expectation values (VEVs),

These VEVs do not break the exchange symmetry and maintain the $Z_{2}$ symmetry, which stabilizes the DM candidate. We parametrized the component fields of each scalar field as

where $\pi^{\pm}_{3},\pi^{0}_{3},\pi^{\pm}_{j},\pi^{\pm}_{j}$, and $\pi^{0}_{j}$ are would-be Nambu-Goldstone (NG) bosons. Based on the stationary condition, we obtain the followings:

We have studied the GWs originating from the phase transition in the dark sector in the electroweakly interacting vector DM model proposed in Abe:2020mph .

We have calculated the effective potential and investigated the phase transition. At the tree level, the potential has negative curvature at the origin. However, the gauge bosons give positive contributions to the potential at the loop level, as shown in Eq. (LABEL:eq:Veff_T=0_around-origin). For the large gauge couplings, which is typically required to obtain the measured value of the DM energy density, the effective potential has positive curvature at the origin, even at $T=0$. As a result, the phase transition in the dark sector is first order and is strong, $\varphi_{C}/T_{C}\gtrsim 1$, in a wide range of the parameter space. The curvature at the tree level is proportional to $m_{h^{\prime}}^{2}$; thus, the loop contributions are significant for smaller $m_{h^{\prime}}$. Consequently, $\varphi_{C}/T_{C}$ is larger for smaller $m_{h^{\prime}}$ as shown in Fig. LABEL:fig:VCTC. We also have found a lower bound on $m_{h^{\prime}}$ for the phase transition. This is because the too small value of $m_{h^{\prime}}$ makes the tunneling rate from the origin to the true vacuum too small, and then the phase transition does not occur.

We have studied three benchmarks ($m_{V}=7,5$, and 3 TeV) and found that the model predicts a GW spectrum that is detectable in the LISA, DECIGO, and BBO experiments. Each benchmark has a different prediction for the $W^{\prime}$ search in the collider experiments. The heavier $V^{0}$ cases cannot be tested by the $W^{\prime}$ searches at the collider experiments, and thus, the GW detection is important to test the model. For $m_{V}=7$ and 3 TeV, we have found that the GW detection can probe the region of the parameter space where the $W^{\prime}$ searches cannot. For $m_{V}=5$ TeV, the region of the parameter space that the GW detection can probe overlaps with the region accessible by the $W^{\prime}$ searches. However, the former region is narrower, and thus the GW is helpful in specifying the model parameters. Assuming the model explains the measured value of the DM energy density via the freeze-out mechanism, we have found that the model predicts the detectable GW signals if $m_{h^{\prime}}$ is a few TeV. Because it is challenging to search heavy $h^{\prime}$ in collider experiments, utilizing the GW signals in determining $m_{h^{\prime}}$ is crucial.

This work was supported by JSPS KAKENHI Grant Number 19H04615 and 21K03549 [T.A.]. We would like to thank Editage (www.editage.jp) for English language editing.