This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

What can we learn from the total width of the Higgs boson?

Qing-Hong Cao [email protected] Department of Physics and State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing 100871, China Collaborative Innovation Center of Quantum Matter, Beijing 100871, China Center for High Energy Physics, Peking University, Beijing 100871, China    Hao-Lin Li [email protected] CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, P. R. China    Ling-Xiao Xu [email protected] [email protected] Department of Physics and State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing 100871, China Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy    Jiang-Hao Yu [email protected] CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, P. R. China School of Physical Sciences, University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing 100049, P.R. China Center for High Energy Physics, Peking University, Beijing 100871, China School of Fundamental Physics and Mathematical Sciences, Hangzhou Institute for Advanced Study, UCAS, Hangzhou 310024, China International Centre for Theoretical Physics Asia-Pacific, Beijing/Hangzhou, China
Abstract

As one of the key properties of the Higgs boson, the Higgs total width is sensitive to global profile of the Higgs boson couplings, and thus new physics would modify the Higgs width. We investigate the total width in various new physics models, including various scalar extension, composite Higgs models, and fraternal twin Higgs model. Typically the Higgs width is smaller than the standard model value due to mixture with other scalar if the Higgs is elementary, or curved Higgs field space for the composite Higgs. On the other hand, except the possible invisible decay mode, the enhanced Yukawa coupling in the two Higgs doublet model or the exotic fermion embeddings in the composite Higgs, could enhance the Higgs width greatly. The precision measurement of the Higgs total width at the high-luminosity LHC can be used to discriminate certain new physics models.

I Introduction

After the discovery of the Higgs boson, the next step is to decipher the particle nature of the Higgs boson, e.g., the mass and width of the Higgs boson (hh). The Higgs boson mass mHm_{H} has been measured very precisely at the Large Hadron Collider (LHC), and its value sheds lights on the stability of the vacuum. On the other hand, the width of the Higgs boson (Γh\Gamma_{h}) is not measured yet. We will focus on in this paper the information we can obtain by precision measurements of the Higgs total width.

At the LHC, information of the total width of the Higgs boson can be extracted from off-shell Higgs production Caola and Melnikov (2013); Campbell et al. (2014a, b); Englert et al. (2015); Logan (2015), the global fit results of the on-shell Higgs signals Barger et al. (2012), t¯tH\bar{t}tH and four-top productions Cao et al. (2017). Measurements of the Higgs total width by the ATLAS and CMS collaboration are found in Refs. Aaboud et al. (2018); Sirunyan et al. (2019), respectively, assuming standard model (SM)-like Higgs couplings. More interestingly, it is worthwhile noting that the Higgs width can also be measured through the Higgs line shape at γγ\gamma\gamma and μ+μ\mu^{+}\mu^{-} colliders Gunion and Haber (1993); Barger et al. (1997).

Assuming the SM-like Higgs couplings, the total width of the Higgs boson Γ\Gamma can be parametrized as

Γ\displaystyle\Gamma ΓSM(cb2Brbb+cτ2Brττ+cc2Brcc\displaystyle\simeq\Gamma_{\text{SM}}\left(c_{b}^{2}\text{Br}_{bb}+c_{\tau}^{2}\text{Br}_{\tau\tau}+c_{c}^{2}\text{Br}_{cc}\right.
+cg2Brgg+cW2BrWW+cZ2BrZZ\displaystyle\ \ \ \ \ \ +\ c_{g}^{2}\text{Br}_{gg}+c_{W}^{2}\text{Br}_{WW}+c_{Z}^{2}\text{Br}_{ZZ}
+cγ2Brγγ+cγZ2BrγZ)+Γinvisible,\displaystyle\left.\ \ \ \ \ \ +\ c_{\gamma}^{2}\text{Br}_{\gamma\gamma}+c_{\gamma Z}^{2}\text{Br}_{\gamma Z}\right)+\Gamma_{\text{invisible}}, (1)

where the major decay modes of the SM Higgs boson, and the invisible width, are included. The numerical values of the SM total width is ΓSM4.1\Gamma_{\text{SM}}\simeq 4.1 MeV, and the major branching ratios in the SM are Brbb0.584,Brττ0.0627,Brcc0.029,Brgg0.0856,BrWW0.214,BrZZ0.0262\text{Br}_{bb}\simeq 0.584,\text{Br}_{\tau\tau}\simeq 0.0627,\text{Br}_{cc}\simeq 0.029,\text{Br}_{gg}\simeq 0.0856,\text{Br}_{WW}\simeq 0.214,\text{Br}_{ZZ}\simeq 0.0262, respectively de Florian et al. (2016). The rescaling factors cic_{i}’s denote the effects of new physics (NP) in the Higgs boson decay. Since Brγγ\text{Br}_{\gamma\gamma} and BrγZ\text{Br}_{\gamma Z} are highly suppressed in the SM, we neglect these two decay channels when considering the total width of the Higgs boson. The invisible decay of the Higgs boson takes place in many NP models, therefore, we also include the invisible width Γinvisible\Gamma_{\text{invisible}} in our study. We refer to Ref. Brivio et al. (2019) for a recent analysis of the Higgs width in the framework of SM effective field theory (EFT).

It remains unknown whether the Higgs boson is fundamental or composite, while the couplings of the Higgs boson to gauge bosons often tend to be smaller than the SM values regardless of this nature. In the case of a fundamental Higgs boson, the decrease of the Higgs-gauge-boson couplings is usually caused by mixing between the Higgs and the new scalar, whereas the exception exist for the Georgi-Machacek model Georgi and Machacek (1985); Chanowitz and Golden (1985) where cW,Zc_{W,Z} can be enhanced. On the other hand, when the Higgs boson is composite and is the pNGB comes from the global symmetry breaking, the Higgs boson to gauge boson couplings decrease due to the misalignment between the gauged direction and the true vacuum direction where the physical Higgs boson fluctuates around. The above arguments for the decrease of the Higgs boson to gauge boson couplings sometimes also hold for the the Higgs boson to fermion couplings, which leads to a decreased Higgs-boson total width if no invisible Hidden sector decay considered. However, in the two Higgs doublet models (2HDMs) Branco et al. (2012), the couplings of the Higgs boson to fermions can be enhanced when deviating from the alignment limit, while in the minimal composite Higgs boson model, an exotic fermion embedding can also leads to the enhanced Yukawa couplings of the top quark to the Higgs boson. Therefore we investigate several popular models that can modify the Higgs-boson total width significantly.

Given the measurement of the total Higgs width, it is possible to discriminate certain classes of new physics models. For the models we studied in this paper, an increased Higgs decay width indicating either an enhanced Yukawa coupling or the existence of invisible decay in the extended scalar models when the Higgs is an elementary particles. When the Higgs is an pNGB, an increased Higgs width usually needs an exotic fermion setup where the SM fermions are embedded in higher dimensional representations. On the other hand, If one observes a smaller Higgs width compared to the SM, a heavy scalar particle that mixes with Higgs is expected to be found for the Higgs as an elementary particle, such a particle may or may not exist in the pNGB Higgs scenario depending on whether the UV is strong coupled or not. In Fig. 1, we present that how different kinds of new physics models are classified based on the Higgs width modification and nature of the Higgs boson.

Refer to caption
Figure 1: Various kinds of new physics models are classified in the flowchart, based on the nature of the Higgs boson, and future measurement of the total Higgs width.

The paper is organized as follows. In Sec. II we discuss the simplest case of the scalar extension, real singlet scalar model Profumo et al. (2007). In Sec. III we focus on four types of two Higgs doublet models (2HDMs), which correspond to four different assignments of Yukawa couplings between fermions and Higgs doublets Branco et al. (2012). In Sec. V we focus on minimal composite Higgs boson model Agashe et al. (2005), in which the Higgs boson emerges as a pseudo Nambu-Goldstone boson (PNGB) from the coset SO(5)/SO(4)SO(5)/SO(4). In Sec. VI we turn our attention to the twin Higgs boson paradigm Chacko et al. (2006). In contrast to the original mirror twin Higgs boson model Chacko et al. (2006), the fraternal twin Higgs-boson model Craig et al. (2015a) is considered, in which only the third generation of the SM fermions have twin partners. Finally, we conclude.

II Real Singlet Model

The real singlet scalar model is the most simplest extension to the SM scalar sector that is possible to generate strong first order electroweak phase transition. It has been extensively studied in the literature O’Connell et al. (2007); Profumo et al. (2007). In this model, after electroweak symmetry breaking, the real singlet scalar SS mixes with the neutral CP-even component h0h_{0} in the Higgs boson doublet such that the physical Higgs boson hh we observed can be expressed as:

h=h0cosθ+Ssinθ.\displaystyle h=h_{0}\cos\theta+S\sin\theta. (2)

Therefore the couplings of the Higgs boson hh to other SM particles are scaled by an overall factor cosθ\cos\theta. If we assume that, in the minimal setup, there is no hidden sector particles that the Higgs boson hh can decays to and the mass of the other singlet-like scalar is larger than half of the SM-like Higgs mass, then the Higgs-boson total width is the SM value scaled by cosθ\cos\theta. In this case, any confirmation of an enhanced Higgs total width can help to ruled out this minimal setup. This conclusion could be extended to extended scalar models, such as the general model setup considered in Ref. Ramsey-Musolf et al. (2021).

Figure 2 displays the ratio of the Higgs boson width to its SM value in the plane of the invisible decay branching ratio BinvB_{inv} and the mixing angle sin2θ\sin^{2}\theta. The chosen ranges of the two axes are within the current LHC constraints  Aad et al. (2019). The pink region denotes the parameter space where the deviation of the total width is within 5% of the SM value such that it cannot be discriminated with the SM with the future high luminosity LHC (HL-LHC). However, if the Higgs boson width is more than the SM value by 25%, then it is likely that there is an invisible decays of the Higgs boson in the singlet model, regardless of the mixing angle θ\theta.

Refer to caption
Figure 2: Contour plots on the total width of the Higgs boson in the singlet scalar extended model in the plane of the mixing angle sin2θ\sin^{2}\theta and the invisible decay branching ratio BinvB_{inv}. The parameter space in the pink region the SM value cannot be differentiated with the SM with the future HL-LHC.

III Two Higgs Doublet Models

We consider the general CP-conserving 2HDM with soft Z2Z_{2} breaking mass term Branco et al. (2012). The Higgs potential is

V(ϕ1,ϕ2)\displaystyle V(\phi_{1},\phi_{2}) =\displaystyle= m112(ϕ1ϕ1)+m222(ϕ2ϕ2)(m122(ϕ1ϕ2)+h.c.)\displaystyle m_{11}^{2}(\phi_{1}^{\dagger}\phi_{1})+m_{22}^{2}(\phi_{2}^{\dagger}\phi_{2})-\left(m_{12}^{2}(\phi_{1}^{\dagger}\phi_{2})+{\rm h.c.}\right) (3)
+λ12(ϕ1ϕ1)2+λ22(ϕ2ϕ2)2+λ3(ϕ1ϕ1)(ϕ2ϕ2)\displaystyle+\frac{\lambda_{1}}{2}(\phi_{1}^{\dagger}\phi_{1})^{2}+\frac{\lambda_{2}}{2}(\phi_{2}^{\dagger}\phi_{2})^{2}+\lambda_{3}(\phi_{1}^{\dagger}\phi_{1})(\phi_{2}^{\dagger}\phi_{2})
+λ4(ϕ1ϕ2)(ϕ2ϕ1)+12[λ5(ϕ1ϕ2)2+h.c.],\displaystyle+\lambda_{4}(\phi_{1}^{\dagger}\phi_{2})(\phi_{2}^{\dagger}\phi_{1})+\frac{1}{2}\left[\lambda_{5}(\phi_{1}^{\dagger}\phi_{2})^{2}+{\rm h.c.}\right]\ ,

where all the parameters in the above potential are assumed to be real. Four types of Yukawa interactions can be introduced with different assignments of the Z2Z_{2} charge of the fermion fields to prevent the tree-level flavor changing neutral current (FCNC). In each model, different types of the right handed fermion fields couple to different Higgs doublets; see Table 1.

Table 1: Models without tree-level FCNC where the uRiu_{R}^{i}, dRid_{R}^{i} and eRie_{R}^{i} represents the right-handed up-type quark, down-type quark and charged lepton fields, respectively. In different types of models they couple to different Higgs doublets.
Model uRiu_{R}^{i} dRid_{R}^{i} eRie_{R}^{i} Model uRiu_{R}^{i} dRid_{R}^{i} eRie_{R}^{i}
Type-I ϕ2\phi_{2} ϕ2\phi_{2} ϕ2\phi_{2} Type-II ϕ2\phi_{2} ϕ1\phi_{1} ϕ1\phi_{1}
Type-X ϕ2\phi_{2} ϕ2\phi_{2} ϕ1\phi_{1} Type-Y ϕ2\phi_{2} ϕ1\phi_{1} ϕ2\phi_{2}

People usually trade the parameters in the potential with a set of physical parameters: vv, the electroweak vacuum expectation value (vev); tanβ\tan\beta, the ratio of the vevs of two Higgs doublets; α\alpha, the rotation angle which diagonalizes the mass matrix of the CP-even neutral Higgs sector; m1(2,3)m_{1(2,3)}, the mass for the tree neutral Higgs, where m1=125.1m_{1}=125.1 GeV is identified with the mass of the SM-like Higgs boson, m2m_{2} and m3m_{3} are masses of the heavy CP even and CP-odd Higgs boson, respectively; mH±m_{H^{\pm}}, the mass of the charged Higgs boson; m122m_{12}^{2}, the parameters that is sensitive to several the theoretical bonds of the theory (perturbative unitarity, stability) Kling et al. (2016). In addition, we denotes the physical SM-like Higgs boson as hh and the charged Higgs boson as H±H^{\pm}.

The Higgs boson width in the 2HDM is dominantly determined by the parameters α\alpha, β\beta, which fix the Higgs couplings to fermions and gauge bosons. The dependence on m122m_{12}^{2}, m2m_{2} and mH±m_{H^{\pm}} appears when taking into account the partial widths of hγγh\to\gamma\gamma and hZγh\to Z\gamma. We define the rescaling of the SM Higgs coupling and the hh to H±H^{\pm} couplings as follows:

h\displaystyle{\mathcal{L}}_{h} =\displaystyle= f=u,d,mfvcff¯fhV=Z,W(1+δV)mV2vcVVμVμh\displaystyle-\sum_{f=u,d,\ell}\frac{m_{f}}{v}c_{f}{\overline{f}}fh-\sum_{V=Z,W}\frac{(1+\delta_{V})m_{V}^{2}}{v}c_{V}V^{\mu}V_{\mu}h (4)
gcvhH+H,\displaystyle-\frac{g_{c}}{v}hH^{+}H^{-},

where the rescaling factors of the Yukawa couplings cfc_{f}’s in each models are summarized in Table 2, and the cVc_{V} and gcg_{c} are universal for each model,

cV\displaystyle c_{V} =\displaystyle= sin(βα),\displaystyle\sin(\beta-\alpha), (5)
gc\displaystyle g_{c} =\displaystyle= 2(m122m122sin2β)cot2βcos(βα)\displaystyle 2\left(m_{1}^{2}-\frac{2m_{12}^{2}}{\sin 2\beta}\right)\cot 2\beta\cos(\beta-\alpha) (6)
+\displaystyle+ (m12+2mH±24m122sin2β)sin(βα).\displaystyle\left(m_{1}^{2}+2m_{H^{\pm}}^{2}-\frac{4m_{12}^{2}}{\sin 2\beta}\right)\sin(\beta-\alpha).
Table 2: The rescaling factor of the Yukawa couplings to the SM value.
cuc_{u} cdc_{d} clc_{l}
Type-I cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha} cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha} cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha}
Type-II cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha} sβαcβαtβs_{\beta-\alpha}-c_{\beta-\alpha}t_{\beta} sβαcβαtβs_{\beta-\alpha}-c_{\beta-\alpha}t_{\beta}
Type-X cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha} cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha} sβαcβαtβs_{\beta-\alpha}-c_{\beta-\alpha}t_{\beta}
Type-Y cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha} sβαcβαtβs_{\beta-\alpha}-c_{\beta-\alpha}t_{\beta} cβα/tβ+sβαc_{\beta-\alpha}/t_{\beta}+s_{\beta-\alpha}

We use the public code 2HDMC-1.7.0 Eriksson et al. (2010) to calculate the total width of the Higgs boson numerically. Figure 3 plots the contours of the ratio of the total width to the SM value in four types of models in the plane of cos(βα)\cos(\beta-\alpha) and tanβ\tan\beta. Different color lines represent the contours of the ratio, the gray region is excluded by the global fit of the single Higgs boson production, which is obtained directly from Fig. 1 in Ref. Chowdhury and Eberhardt (2018). In each model, we have set mass parameters as: mH±=m2=m3=1m_{H^{\pm}}=m_{2}=m_{3}=1 TeV, m122=m22sinβcosβm_{12}^{2}=m_{2}^{2}\sin\beta\cos\beta. The reason that we chose such a m122m_{12}^{2} is to satisfy the unitarity constraint in large tanβ\tan\beta regions Kling et al. (2016). The mass difference between the charged Higgs boson and the neutral heavy Higgs boson should be less than 300 GeV, i.e., |mH±m2,3|300GeV|m_{H^{\pm}}-m_{2,3}|\leq 300~{}{\rm GeV} Chowdhury and Eberhardt (2018). The range of the tanβ\tan\beta is taken to be (0.3, 32), this region is well consistent with the region allowed by 95.4% confidential level when taking into account all the available constraints Chowdhury and Eberhardt (2018).

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 3: Contour of the ratio of the total width of the Higgs boson to the SM value in various 2HDMs in the cos(βα)\cos(\beta-\alpha) vs tanβ\tan\beta plane. We have set the mass parameters in each types of model as: mH±=m2=m3=1m_{H^{\pm}}=m_{2}=m_{3}=1 TeV, m122=1m_{12}^{2}=1 TeV2.

In general, if one neglects the partial width contributions from hγγh\to\gamma\gamma and hZγh\to Z\gamma, then the Higgs-boson total width is uniquely determined by two parameters: cos(βα)\cos(\beta-\alpha) and tanβ\tan\beta, while the information of other parameters enters via the coupling of the charged Higgs boson to the SM Higgs boson gcg_{c} when these two decay channels hγγh\to\gamma\gamma and hZγh\to Z\gamma are taken into account.

Figure 3 shows that ratio of the Higgs-boson total width to the SM value in the 2HDMs roughly ranges from 0.8 to 1.25, and only the Type-I model can have the ratio less than 0.85; on the other hand, only the Type-X model can have the ratio larger than 1.15 . Therefore, all the four types of models are likely to be excluded if the LHC find a ratio larger than 1.25 or less than 0.8 . If one find the ratio is in the range of (1.15,1.25) ((0.8,0.85)), then it must be Type-X (Type-I) model provided the CP conserving 2HDM were realized in the nature.

Figure 3 also demonstrates that in the Type-I model only a positive cos(βα)\cos(\beta-\alpha) combined with a rather small tanβ\tan\beta can give an enhanced Higgs total width, the reason is simply that all the fermion couplings are enhanced in this region. However, in the Type-II and Type-X model, an enhanced Higgs total width can only obtained by a negative cos(βα)\cos(\beta-\alpha), the main reason is that the couplings to the bottom quarks are enhanced. In Type-X model, one can obtain an enhanced Higgs total width with either a positive cos(βα)\cos(\beta-\alpha) combined with a relatively small tanβ\tan\beta or a negative cos(βα)\cos(\beta-\alpha) combined with a relatively large tanβ\tan\beta. In the former case, the increase of the Higgs width is mainly due to enhanced couplings to the top quarks which further increases the decay width to gluons, while for the later case, the increase of the Higgs width is mainly caused by the enhanced decay width to bottom quarks. In addition, in Type-II, Type-X and Type-Y models, one can find a nearly vertical line corresponding to the ratio of the Higgs width equal to one, the reason is that the increase of decay widths of some channels due to increased fermions couplings is exactly compensated by the reduction in decay widths of some other channels caused by decreased Higgs couplings to other fermions and gauge bosons.

In the decoupling limit, where we assume the masses of the three Higgs are equal and much larger than the electroweak scale: m2=m3=mH±Λvm_{2}=m_{3}=m_{H^{\pm}}\sim\Lambda\gg v. In this case, one can work in an effective field theory derived from the Higgs basis, where cos(βα)\cos(\beta-\alpha) scales as:

cos(βα)=Z6v2Λ2,\displaystyle\cos(\beta-\alpha)=-Z_{6}\frac{v^{2}}{\Lambda^{2}}, (7)

where Z6Z_{6} is the coefficient of the |H1|2H1H2|H_{1}|^{2}H_{1}^{\dagger}H_{2} which can be either positive or negative and Λ2\Lambda^{2} is the coefficient of |H2|2|H_{2}|^{2} in the Higgs basis Bélusca-Maïto et al. (2017). Then one can translate the bond on the Higgs width onto the new physics scale Λ\Lambda with different value of tanβ\tan\beta assuming Z6Z_{6} a 𝒪(1){\cal O}(1) parameter. We plot the contours of the ratio of the Higgs width to the SM value in different models in Fig. 4 with Z6=1Z_{6}=1, the gray regions are excluded by the Higgs global fit of the LHC Run1 and Run2 results as above Chowdhury and Eberhardt (2018), the pink region represents the parameter space in which the deviation of the Higgs total width is within 5% of the SM value so that the new physics effects cannot be probed by the measurement of the Higgs total width at the HL-LHC.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 4: Contour of the ratio of the total width of the Higgs boson to the SM value in various 2HDMs with the parameter Z6=1Z_{6}=1. The gray regions are excluded by the Higgs global fit, the pink regions are the parameter spaces such that the deviation of the Higgs total width from the SM value are too small to probe at future HL-LHC.

IV Georgi-Machacek Model

In the Georgi-Machacek Model, the scalar sector of the SM is extended with a real triplet (ξ+,ξ0,ξ)(\xi^{+},\xi^{0},\xi^{-}) with hypercharge Y=0Y=0 and a complex triplet (χ++,χ+,χ0\chi^{+}+,\chi^{+},\chi^{0}) with hypercharge Y=2Y=2. A good review of the model can be found in Ref Hartling et al. (2014). After electroweak symmetry breaking, the Higgs doublet gets a vev vHv_{H} and the neutral components of the real and complex triplet get an equal vev vχv_{\chi}. The relation between the electroweak vev vv defined by the Fermi constant GFG_{F} and vHv_{H} and vχv_{\chi} is given by:

v2=vH2+8vχ2.\displaystyle v^{2}=v_{H}^{2}+8v_{\chi}^{2}. (8)

The rescaling factor for the gauge boson and fermions are given by:

cW,Z\displaystyle c_{W,Z} =\displaystyle= e22sW2cW2(cosαsinβ263sinαcosβ)\displaystyle\frac{e^{2}}{2s_{W}^{2}c_{W}^{2}}\left(\cos\alpha\sin\beta-\frac{2\sqrt{6}}{3}\sin\alpha\cos\beta\right) (9)
cf\displaystyle c_{f} =\displaystyle= cosαsinβ,\displaystyle\frac{\cos\alpha}{\sin\beta}, (10)

where the angle α\alpha related to the diagonalization of the mass matrix of the neutral scalars, and cosβ\cos\beta is defined by the ratio of vHv_{H} and vv. From Eq. 10, one can find that the Higgs boson couplings of fermions and gauge bosons can be larger than 1, which may lead to an enhanced Higgs boson width.

To analyze the Higgs boson width in the model, we use GMCALC-1.4.1 Hartling et al. (2014) to scan the allowed parameter space. The scanning method and the parameter ranges are given in Ref. Hartling et al. (2014). Each parameter point in the scan is required to satisfy various constraints: the correct electroweak vev, the stability of the scalar potential, the tree-level unitarity, the experimental bounds on the SS-parameter and bsγb\to s\gamma. Figure 5 shows the results of the parameter scan in the plane of sinβ\sin\beta and sinα\sin\alpha (left) and in the plane of cfc_{f} and cVc_{V} (right). The blue dots represents the ratio Γ/ΓSM\Gamma/\Gamma_{SM} within the range (0.95,1.05)(0.95,1.05) and are not be able to discriminant from the SM value. The red and green dots denotes the ratio Γ/ΓSM\Gamma/\Gamma_{SM} larger than 1.05 and smaller 0.95, respectively, which are possible to be probed by the future experiments. We also find that enhanced fermion couplings strongly correlate with an enlarged Higgs width, while an enlarged Higgs width does not necessarily corresponds to enhanced gauge boson couplings. Especially for a Higgs boson width enlarged by a factor of 1.5 or more, a decreased gauge boson coupling is observed. We also plot the scanned points in the plane of sinβ\sin\beta and sinα\sin\alpha. One can find from the plot that these two angle parameters almost determine the value of the Higgs boson width, as the parameter points with different colors are well separated. As one expected, the SM value corresponds to the point with sinβ=1\sin\beta=1 and cosα=0\cos\alpha=0, and the red dots corresponding to an increased Higgs boson width tend to have smaller sinβ\sin\beta’s.

Refer to caption
Refer to caption
Figure 5: The scattered plot of the scanned points on the sinβ\sin\beta vs sinα\sin\alpha plane (left) and cfc_{f} vs cVc_{V} plane (right), each points satisfies the theoretical bounds and indirect experimental bounds. The blue dots are those with Γ/ΓSM\Gamma/\Gamma_{SM} within [0.95, 1.05], which cannot be distinguished by the future HL-LHC experiment. The red and green dots corresponds to Γ/ΓSM>1.05\Gamma/\Gamma_{SM}>1.05 and Γ/ΓSM<0.95\Gamma/\Gamma_{SM}<0.95 respectively and are possible to be probed by the future HL-LHC experiments.

V Minimal Composite Higgs Models

If the Higgs boson is a pseudo Nambu-Goldstone boson (PNGB) emerging from strong dynamics at TeV scale, its coupling to light fermions b,c,τb,c,\tau and electroweak gauge bosons W±,ZW^{\pm},Z would be modified with respect to their SM values. In this case, the sign and magnitude of Higgs coupling modifications are mainly dictated by two effects: the Higgs nonlinearity and the composite resonances Li et al. (2019).

The Higgs nonlinearity denotes the nontrivial curvature of the coset space from which the PNGB Higgs boson emerges Alonso et al. (2016). To be concrete, the compact cosets have positive curvatures, while the non-compact ones have negative curvatures. Among various cosets Bellazzini et al. (2014), SO(5)/SO(4)SO(5)/SO(4) receives the most attention, and it is also known as the minimal coset that accommodates the custodial symmetry. Therefore this is called as the minimal composite Higgs model Agashe et al. (2005). As a result, the Higgs coupling to the electroweak gauge bosons W±,ZW^{\pm},Z are universally shifted Liu et al. (2018) as

cW=cZ=κhWWκhWWSM=κhZZκhZZSM=1ξ\displaystyle c_{W}=c_{Z}=\frac{\kappa_{hWW}}{\kappa^{\text{SM}}_{hWW}}=\frac{\kappa_{hZZ}}{\kappa^{\text{SM}}_{hZZ}}=\sqrt{1-\xi} (11)

at the leading order of the chiral expansion, due to Higgs nonlinearity. Here ξv2/f2\xi\equiv v^{2}/f^{2} is defined as the ratio of the electroweak scale and the decay constant of the PNGB Higgs boson.

In the fermionic sector, modifications of Higgs boson couplings also depend on fermion embeddings into SO(5)SO(5) multiplets. In this work, we consider the cases where the top quark and the light fermions b,c,τb,c,\tau are all embedded in the fundamental representation 55, or the symmetric tensor representation 1414, of the SO(5)SO(5) group. Under the paradigm of partial compositeness, there are composite particles mixing with the elementary t,b,c,τt,b,c,\tau, which are usually considered as the source of explicit SO(5)SO(5) breaking in the fermionic sector; otherwise the Higgs boson would be a massless exact Goldstone boson when these mixings (as well as gauge couplings) are turned off. After integrating out these composite fermions, they would contribute to the low energy effective couplings of the Higgs boson, through both the chirality-conserving wave functions and the chirality-flipping Yukawa vertices of the elementary fermions. Following Ref. Li et al. (2019), the chirality-conserving wave functions are universally expanded as

ΠfL=Π0fL+Π1fLsh2+Π2fLsh4+,ΠfR=Π0fR+Π1fRsh2+Π2fRsh4+;\displaystyle\begin{aligned} \Pi_{f_{L}}&=\Pi_{0f_{L}}+\Pi_{1f_{L}}\ s^{2}_{h}+\Pi_{2f_{L}}\ s^{4}_{h}+\cdots,\ \\ \Pi_{f_{R}}&=\Pi_{0f_{R}}+\Pi_{1f_{R}}\ s^{2}_{h}+\Pi_{2f_{R}}\ s^{4}_{h}+\cdots;\end{aligned} (12)

while the chirality-flipping Yukawa interactions satisfy the expansion

ΠtLtR\displaystyle\Pi_{t_{L}t_{R}} =Π1tLtRchsh+Π2tLtRchsh3+,\displaystyle=\Pi_{1t_{L}t_{R}}\ c_{h}s_{h}+\Pi_{2t_{L}t_{R}}\ c_{h}s^{3}_{h}+\cdots\ , (13)

or

ΠtLtR\displaystyle\Pi_{t_{L}t_{R}} =Π1tLtRsh+Π2tLtRsh3+,\displaystyle=\Pi_{1t_{L}t_{R}}\ s_{h}+\Pi_{2t_{L}t_{R}}\ s^{3}_{h}+\cdots\ , (14)

depending on fermion representations. For the cases considered in this work, (t,b,c,τ)L,R5(t,b,c,\tau)_{L,R}\subset 5 and (t,b,c,τ)L,R14(t,b,c,\tau)_{L,R}\subset 14 corresponds to the Eq. 13, while (t,b,c,τ)L5,(t,b,c,τ)R1(t,b,c,\tau)_{L}\subset 5,(t,b,c,\tau)_{R}\subset 1 corresponds to Eq. 14, as discussed in Ref. Li et al. (2019).

Accordingly, up to the linear order of ξ\xi, the fermion couplings are modified as following:

cf\displaystyle c_{f} κfκfSM\displaystyle\equiv\frac{\kappa_{f}}{\kappa_{f}^{\text{SM}}}
=132ξξ(Π1fLΠ0fL+Π1fRΠ0fR)+2ξΠ2fLfRΠ1fLfR.\displaystyle=1-\frac{3}{2}\xi-\xi\left(\frac{\Pi_{1f_{L}}}{\Pi_{0f_{L}}}+\frac{\Pi_{1f_{R}}}{\Pi_{0f_{R}}}\right)+2\xi\frac{\Pi_{2f_{L}f_{R}}}{\Pi_{1f_{L}f_{R}}}. (15)

By naive dimensional analysis,

Π1fLyL2f2m2,Π1fRyR2f2m2,Π0fLΠ0fR1,\displaystyle\Pi_{1f_{L}}\sim\frac{y_{L}^{2}f^{2}}{m_{*}^{2}},\ \Pi_{1f_{R}}\sim\frac{y_{R}^{2}f^{2}}{m_{*}^{2}},\ \Pi_{0f_{L}}\sim\Pi_{0f_{R}}\sim 1, (16)

where the mass scale mgfm_{*}\sim g_{*}f and yL,Ry_{L,R} are the mixing parameters, while

Π1fLfRyLyRf2m,Π2fLfRyLyRf2m.\displaystyle\Pi_{1f_{L}f_{R}}\sim\frac{y_{L}y_{R}f^{2}}{m_{*}},\ \Pi_{2f_{L}f_{R}}\sim\frac{y_{L}y_{R}f^{2}}{m_{*}}. (17)

For the light fermions b,c,τb,c,\tau, due to their small masses mb,c,τyLyRgfm_{b,c,\tau}\sim\frac{y_{L}y_{R}}{g_{*}}f, the mixing parameters yL,Ry_{L,R} are expected to be much smaller than the coupling strength gg_{*} if the decay constant ff is at the order of 11 TeV. As a result, one can neglect the effect of composite resonances in the wave functions of b,c,τb,c,\tau. Furthermore, Π2fLfR\Pi_{2f_{L}f_{R}} does not vanish only in the case of (t,b,c,τ)L,R14(t,b,c,\tau)_{L,R}\subset 14.

The effective coupling between the Higgs boson and the gluon is, up to the linear order of ξ\xi,

cg\displaystyle c_{g} κhggκhggSM\displaystyle\equiv\frac{\kappa_{hgg}}{\kappa^{\text{SM}}_{hgg}}
=132ξ+2ξΠ2tLtRΠ1tLtR+f=b,c,τξ(Π1fLΠ0fL+Π1fRΠ0fR).\displaystyle=1-\frac{3}{2}\xi+2\xi\frac{\Pi_{2t_{L}t_{R}}}{\Pi_{1t_{L}t_{R}}}+\sum_{f=b,c,\tau}\xi\left(\frac{\Pi_{1f_{L}}}{\Pi_{0f_{L}}}+\frac{\Pi_{1f_{R}}}{\Pi_{0f_{R}}}\right)\ . (18)

Note that contribution of the composite resonances in the b,c,τb,c,\tau sectors (but not the light fermions b,c,τb,c,\tau) are included in the above equation Liu et al. (2017).

Therefore, the Higgs effective couplings in the fermionic sectors are concretely

  • (t,b,c,τ)L,R5(t,b,c,\tau)_{L,R}\subset 5:

    cb,c,τ=cg132ξ.\displaystyle c_{b,c,\tau}=c_{g}\simeq 1-\frac{3}{2}\xi\ . (19)
  • (t,b,c,τ)L5(t,b,c,\tau)_{L}\subset 5 and (t,b,c,τ)R1(t,b,c,\tau)_{R}\subset 1:

    cb,c,τ=cg112ξ.\displaystyle c_{b,c,\tau}=c_{g}\simeq 1-\frac{1}{2}\xi\ . (20)
  • (t,b,c,τ)L,R14(t,b,c,\tau)_{L,R}\subset 14:

    cb,c,τ=cg132ξ+2r12ξ,\displaystyle c_{b,c,\tau}=c_{g}\simeq 1-\frac{3}{2}\xi+2\ r_{12}\ \xi\ , (21)

    where r12r_{12} is the ratio of Π2fLfR\Pi_{2f_{L}f_{R}} to Π1fLfR\Pi_{1f_{L}f_{R}}.

With the above results, we obtain the total width of the Higgs boson in minimal composite Higgs models, as shown in Fig. 6. We see explicitly that ξ0\xi\to 0 corresponds to the decoupling limit where one cannot distinguish a composite Higgs boson from an elementary one. For small fermion representation, the total width of the Higgs boson tends to be smaller than the value in the SM; while for symmetric tensor representation, the total width could be larger.

Refer to caption
Figure 6: The total width of the Higgs boson in minimal composite Higgs models, in which the left-handed and right-handed fermions are embedded in different representations of the SO(5)SO(5) group.

VI Fraternal Twin Higgs Model

It is fascinating to consider models in which the Higgs boson could decay into invisible particles. Typically invisible decay exist in neutral naturalness models Chacko et al. (2006); Craig et al. (2015b); Xu et al. (2020). One of these typical model is the twin Higgs scenario Chacko et al. (2006), where the naturalness problem is addressed with some particles which are not charged under the SM gauge groups.

In this section we consider the fraternal twin Higgs model Craig et al. (2015a). In contrast to the original mirror twin Higgs setup Chacko et al. (2006), the fraternal model is more minimal and cosmologically safe. The model ingredients of the twin sector consists of one additional Higgs doublet, twin fermions (including t~R,b~R,τ~R,Q~L,L~L\widetilde{t}_{R},\widetilde{b}_{R},\widetilde{\tau}_{R},\widetilde{Q}_{L},\widetilde{L}_{L}), twin weak gauge bosons from the SU(2)SU(2) gauge group in the twin sector, and twin gluons g~\widetilde{g} Craig et al. (2015a). The Higgs boson is identified as a pseudo Nambu-Goldstone boson from the coset SU(4)/SU(3)SU(4)/SU(3), or equivalently SO(8)/SO(7)SO(8)/SO(7), depending whether the custodial symmetry is included in the unbroken group. After the global symmetry breaking, three of the Goldstone bosons would be eaten by the twin weak gauge bosons, which makes their masses roughly at the order of symmetry-breaking scale ff. For reasonable values of ff, these twin gauge bosons can be sufficiently heavy. Only the other light degrees of freedom, lighter than mh/2m_{h}/2, could contribute to the invisible width of the Higgs boson. On top of this, all the SM couplings are modified due to the Higgs nonlinearity.

Up to the linear order of ξ=v2/f2\xi=v^{2}/f^{2}, all the relevant Higgs couplings are as follows: i) for the couplings in the SM sector,

cb,c,τ\displaystyle c_{b,c,\tau} =cg112ξ,\displaystyle=c_{g}\simeq 1-\frac{1}{2}\xi\ ,
cW\displaystyle c_{W} =cZ=1ξ;\displaystyle=c_{Z}=\sqrt{1-\xi}; (22)

ii) for the couplings in the twin sector,

cb~,τ~=cg~κhb~b~,hτ~τ~,hg~g~κhbb,hττ,hggSMξ,\displaystyle c_{\widetilde{b},\widetilde{\tau}}=c_{\widetilde{g}}\equiv\frac{\kappa_{h\widetilde{b}\widetilde{b},h\widetilde{\tau}\widetilde{\tau},h\widetilde{g}\widetilde{g}}}{\kappa^{\text{SM}}_{hbb,h\tau\tau,hgg}}\simeq\sqrt{\xi}, (23)

assuming the Yukawa couplings are the same, i.e. yt,b,τ=yt~,b~,τ~y_{t,b,\tau}=y_{\widetilde{t},\widetilde{b},\widetilde{\tau}}. All the above Higgs boson couplings in the twin sector are correlated to their counterparts in the SM sector, in order to realize the naturalness condition Li et al. (2019); as one can see that cb~,τ~,g~2+cb,τ,g2c^{2}_{\widetilde{b},\widetilde{\tau},\widetilde{g}}+c^{2}_{b,\tau,g} does not depend on ξ\xi. Furthermore, the masses of the twin fermions are also connected to the partners in the SM sector, e.g.,

mt~,b~,τ~mt,b,τ1ξξfv,\displaystyle\frac{m_{\widetilde{t},\widetilde{b},\widetilde{\tau}}}{m_{t,b,\tau}}\simeq\frac{\sqrt{1-\xi}}{\sqrt{\xi}}\sim\frac{f}{v}, (24)

assuming the Yukawa couplings are the same in the SM and the twin sectors. The twin gluons are massless. They could form glueballs by themself or other bound states with twin fermions below the scale of twin confinement Craig et al. (2015a). Note that mb~m_{\widetilde{b}} and mτ~m_{\widetilde{\tau}} can be larger or smaller than mh/2m_{h}/2 depending on the scale ff, which would result in different invisible widths of the Higgs boson.

The particles in the twin sector induce the invisible width as following:

ΓinvisibleΓSMξ\displaystyle\Gamma_{\text{invisible}}\simeq\Gamma_{\text{SM}}\xi [Brgg\displaystyle\Bigg{[}~{}~{}\text{Br}_{gg}
+Brττ(1xτ~1xτ)32Θ(mh2mτ~)\displaystyle+\text{Br}_{\tau\tau}\left(\frac{1-x_{\widetilde{\tau}}}{1-x_{\tau}}\right)^{\frac{3}{2}}\Theta(m_{h}-2m_{\widetilde{\tau}})
+Brbb(1xb~1xb)32Θ(mh2mb~)],\displaystyle+\text{Br}_{bb}\left(\frac{1-x_{\widetilde{b}}}{1-x_{b}}\right)^{\frac{3}{2}}\Theta(m_{h}-2m_{\widetilde{b}})\Bigg{]}, (25)

where xi4mi2/mh2x_{i}\equiv 4m_{i}^{2}/m_{h}^{2} for the ii-particle and Brjj\text{Br}_{jj} denotes the branching ratio of the decay channel of hjjh\to jj in the SM.

Note that the running masses of the fermions bb and τ\tau in the SM model are roughly mb(MZ)2.9m_{b}(M_{Z})\simeq 2.9 GeV and mτ(MZ)1.746m_{\tau}(M_{Z})\simeq 1.746 GeV, respectively. We calculate the total width of the Higgs boson in the fraternal twin Higgs model from Eq. 24 and Eq. 25. Figure 7 plots the ratio of the total width of the Higgs boson to the SM value as a function of ξ\xi. We notice that the total width is moderately smaller than the SM value, therefore, it is very hard to distinguish the fraternal twin Higgs model from the SM with only the information of Higgs total width.

Refer to caption
Figure 7: The total width of the Higgs boson in fraternal twin Higgs model.

VII Conclusion

The nature of the Higgs boson remains as one unsolved puzzle nowadays. Motivated from the future precise measurement of the Higgs width in the future HL-LHC, we investigated that if the Higgs is elementary or composite, how the model setup could affect the Higgs total width. Typically the Higgs total width is expected to be smaller than the standard model values because of the mixture with other scalars or curved Higgs field space for the composite Higgs. If there is invisible decay channels, the Higgs width could be enhanced, as expected. We further showed that there are also cases that the Higgs width is enhanced due to model setup.

Depending on the model setup, the following results are in order:

  • For the singlet extended model, and typical scalar extension models other than the doublet and GM models, the Higgs width is smaller than the SM if no other invisible decay presents.

  • In the general two Higgs doublet models, the Higgs width can be larger or smaller than the SM value, depending on values of the model parameters. It depends on the concrete value of the Higgs width measured to determine whether a given type of the two Higgs doublet models are excluded or not.

  • In the Georgi-Machacek model, one will find that an increased Higgs boson coupling always accompanies an enhance Higgs to fermion couplings.

  • For the case of a pseudo-Nambu-Goldstone Higgs boson, including both the minimal composite Higgs model and the twin Higgs model, the Higgs total width is strongly preferred to be smaller than the SM prediction except for the case where the left handed and right handed top quark are both embedded in the 14 representation of SO(5). Therefore, only if the Higgs width is measured to be smaller than the SM value, the minimal composite Higgs model and the twin Higgs model are favored.

Overall, if the Higgs total width is measured to be significantly larger or smaller than the SM value, it is possible for us to falsify several new physics models considered above, although the caveat still exists, in which new physics might contain some new hidden sector we never consider and thus the SM Higgs can decay invisibly.

We note that the Higgs total width is sensitive to almost all kinds of Higgs couplings. On the other hand, the Higgs production and partial decay width provide different kinds of information on the Higgs coupling. We expect that combining the Higgs coupling information and the Higgs total width information at the HL-LHC, the Higgs sector could be better understood in near future.

Acknowledgements.
The work is supported in part by the National Science Foundation of China under Grant Nos. 11725520, 11675002, 11635001, 12022514, 11875003 and 12047503. JHY is also supported by the National Key Research and Development Program of China under Grant No. 2020YFC2201501.

References

  • Caola and Melnikov (2013) F. Caola and K. Melnikov, Phys. Rev. D88, 054024 (2013), eprint 1307.4935.
  • Campbell et al. (2014a) J. M. Campbell, R. K. Ellis, and C. Williams, JHEP 04, 060 (2014a), eprint 1311.3589.
  • Campbell et al. (2014b) J. M. Campbell, R. K. Ellis, and C. Williams, Phys. Rev. D89, 053011 (2014b), eprint 1312.1628.
  • Englert et al. (2015) C. Englert, Y. Soreq, and M. Spannowsky, JHEP 05, 145 (2015), eprint 1410.5440.
  • Logan (2015) H. E. Logan, Phys. Rev. D92, 075038 (2015), eprint 1412.7577.
  • Barger et al. (2012) V. Barger, M. Ishida, and W.-Y. Keung, Phys. Rev. Lett. 108, 261801 (2012), eprint 1203.3456.
  • Cao et al. (2017) Q.-H. Cao, S.-L. Chen, and Y. Liu, Phys. Rev. D95, 053004 (2017), eprint 1602.01934.
  • Aaboud et al. (2018) M. Aaboud et al. (ATLAS), Phys. Lett. B786, 223 (2018), eprint 1808.01191.
  • Sirunyan et al. (2019) A. M. Sirunyan et al. (CMS), Phys. Rev. D99, 112003 (2019), eprint 1901.00174.
  • Gunion and Haber (1993) J. F. Gunion and H. E. Haber, Phys. Rev. D48, 5109 (1993).
  • Barger et al. (1997) V. D. Barger, M. S. Berger, J. F. Gunion, and T. Han, Phys. Rept. 286, 1 (1997), eprint hep-ph/9602415.
  • de Florian et al. (2016) D. de Florian et al. (LHC Higgs Cross Section Working Group) (2016), eprint 1610.07922.
  • Brivio et al. (2019) I. Brivio, T. Corbett, and M. Trott (2019), eprint 1906.06949.
  • Georgi and Machacek (1985) H. Georgi and M. Machacek, Nucl. Phys. B262, 463 (1985).
  • Chanowitz and Golden (1985) M. S. Chanowitz and M. Golden, Phys. Lett. 165B, 105 (1985).
  • Branco et al. (2012) G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher, and J. P. Silva, Phys. Rept. 516, 1 (2012), eprint 1106.0034.
  • Profumo et al. (2007) S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, JHEP 08, 010 (2007), eprint 0705.2425.
  • Agashe et al. (2005) K. Agashe, R. Contino, and A. Pomarol, Nucl. Phys. B719, 165 (2005), eprint hep-ph/0412089.
  • Chacko et al. (2006) Z. Chacko, H.-S. Goh, and R. Harnik, Phys. Rev. Lett. 96, 231802 (2006), eprint hep-ph/0506256.
  • Craig et al. (2015a) N. Craig, A. Katz, M. Strassler, and R. Sundrum, JHEP 07, 105 (2015a), eprint 1501.05310.
  • O’Connell et al. (2007) D. O’Connell, M. J. Ramsey-Musolf, and M. B. Wise, Phys. Rev. D75, 037701 (2007), eprint hep-ph/0611014.
  • Ramsey-Musolf et al. (2021) M. J. Ramsey-Musolf, J.-H. Yu, and J. Zhou (2021), eprint 2104.10709.
  • Aad et al. (2019) G. Aad et al. (ATLAS) (2019), eprint 1909.02845.
  • Kling et al. (2016) F. Kling, J. M. No, and S. Su, JHEP 09, 093 (2016), eprint 1604.01406.
  • Eriksson et al. (2010) D. Eriksson, J. Rathsman, and O. Stal, Comput. Phys. Commun. 181, 189 (2010), eprint 0902.0851.
  • Chowdhury and Eberhardt (2018) D. Chowdhury and O. Eberhardt, JHEP 05, 161 (2018), eprint 1711.02095.
  • Bélusca-Maïto et al. (2017) H. Bélusca-Maïto, A. Falkowski, D. Fontes, J. C. Romão, and J. P. Silva, Eur. Phys. J. C77, 176 (2017), eprint 1611.01112.
  • Hartling et al. (2014) K. Hartling, K. Kumar, and H. E. Logan (2014), eprint 1412.7387.
  • Li et al. (2019) H.-L. Li, L.-X. Xu, J.-H. Yu, and S.-H. Zhu, JHEP 09, 010 (2019), eprint 1904.05359.
  • Alonso et al. (2016) R. Alonso, E. E. Jenkins, and A. V. Manohar, Phys. Lett. B756, 358 (2016), eprint 1602.00706.
  • Bellazzini et al. (2014) B. Bellazzini, C. Csaki, and J. Serra, Eur. Phys. J. C74, 2766 (2014), eprint 1401.2457.
  • Liu et al. (2018) D. Liu, I. Low, and Z. Yin, Phys. Rev. Lett. 121, 261802 (2018), eprint 1805.00489.
  • Liu et al. (2017) D. Liu, I. Low, and C. E. M. Wagner, Phys. Rev. D96, 035013 (2017), eprint 1703.07791.
  • Craig et al. (2015b) N. Craig, S. Knapen, and P. Longhi, Phys. Rev. Lett. 114, 061803 (2015b), eprint 1410.6808.
  • Xu et al. (2020) L.-X. Xu, J.-H. Yu, and S.-H. Zhu, Phys. Rev. D 101, 095014 (2020), eprint 1810.01882.