What can we learn from the total width of the Higgs boson?
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 (). The Higgs boson mass 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 () 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), 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 and colliders Gunion and Haber (1993); Barger et al. (1997).
Assuming the SM-like Higgs couplings, the total width of the Higgs boson can be parametrized as
(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 MeV, and the major branching ratios in the SM are , respectively de Florian et al. (2016). The rescaling factors ’s denote the effects of new physics (NP) in the Higgs boson decay. Since and 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 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 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.

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 . 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 mixes with the neutral CP-even component in the Higgs boson doublet such that the physical Higgs boson we observed can be expressed as:
(2) |
Therefore the couplings of the Higgs boson to other SM particles are scaled by an overall factor . If we assume that, in the minimal setup, there is no hidden sector particles that the Higgs boson 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 . 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 and the mixing angle . 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 .

III Two Higgs Doublet Models
We consider the general CP-conserving 2HDM with soft breaking mass term Branco et al. (2012). The Higgs potential is
(3) | |||||
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 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.
Model | Model | ||||||
---|---|---|---|---|---|---|---|
Type-I | Type-II | ||||||
Type-X | Type-Y |
People usually trade the parameters in the potential with a set of physical parameters: , the electroweak vacuum expectation value (vev); , the ratio of the vevs of two Higgs doublets; , the rotation angle which diagonalizes the mass matrix of the CP-even neutral Higgs sector; , the mass for the tree neutral Higgs, where GeV is identified with the mass of the SM-like Higgs boson, and are masses of the heavy CP even and CP-odd Higgs boson, respectively; , the mass of the charged Higgs boson; , 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 and the charged Higgs boson as .
The Higgs boson width in the 2HDM is dominantly determined by the parameters , , which fix the Higgs couplings to fermions and gauge bosons. The dependence on , and appears when taking into account the partial widths of and . We define the rescaling of the SM Higgs coupling and the to couplings as follows:
(4) | |||||
where the rescaling factors of the Yukawa couplings ’s in each models are summarized in Table 2, and the and are universal for each model,
(5) | |||||
(6) | |||||
Type-I | |||
---|---|---|---|
Type-II | |||
Type-X | |||
Type-Y |
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 and . 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: TeV, . The reason that we chose such a is to satisfy the unitarity constraint in large 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., Chowdhury and Eberhardt (2018). The range of the 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).




In general, if one neglects the partial width contributions from and , then the Higgs-boson total width is uniquely determined by two parameters: and , while the information of other parameters enters via the coupling of the charged Higgs boson to the SM Higgs boson when these two decay channels and 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 combined with a rather small 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 , 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 combined with a relatively small or a negative combined with a relatively large . 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: . In this case, one can work in an effective field theory derived from the Higgs basis, where scales as:
(7) |
where is the coefficient of the which can be either positive or negative and is the coefficient of 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 with different value of assuming a parameter. We plot the contours of the ratio of the Higgs width to the SM value in different models in Fig. 4 with , 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.




IV Georgi-Machacek Model
In the Georgi-Machacek Model, the scalar sector of the SM is extended with a real triplet with hypercharge and a complex triplet () with hypercharge . 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 and the neutral components of the real and complex triplet get an equal vev . The relation between the electroweak vev defined by the Fermi constant and and is given by:
(8) |
The rescaling factor for the gauge boson and fermions are given by:
(9) | |||||
(10) |
where the angle related to the diagonalization of the mass matrix of the neutral scalars, and is defined by the ratio of and . 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 -parameter and . Figure 5 shows the results of the parameter scan in the plane of and (left) and in the plane of and (right). The blue dots represents the ratio within the range and are not be able to discriminant from the SM value. The red and green dots denotes the ratio 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 and . 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 and , and the red dots corresponding to an increased Higgs boson width tend to have smaller ’s.


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 and electroweak gauge bosons 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), 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 are universally shifted Liu et al. (2018) as
(11) |
at the leading order of the chiral expansion, due to Higgs nonlinearity. Here 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 multiplets. In this work, we consider the cases where the top quark and the light fermions are all embedded in the fundamental representation , or the symmetric tensor representation , of the group. Under the paradigm of partial compositeness, there are composite particles mixing with the elementary , which are usually considered as the source of explicit 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
(12) |
while the chirality-flipping Yukawa interactions satisfy the expansion
(13) |
or
(14) |
depending on fermion representations. For the cases considered in this work, and corresponds to the Eq. 13, while corresponds to Eq. 14, as discussed in Ref. Li et al. (2019).
Accordingly, up to the linear order of , the fermion couplings are modified as following:
(15) |
By naive dimensional analysis,
(16) |
where the mass scale and are the mixing parameters, while
(17) |
For the light fermions , due to their small masses , the mixing parameters are expected to be much smaller than the coupling strength if the decay constant is at the order of TeV. As a result, one can neglect the effect of composite resonances in the wave functions of . Furthermore, does not vanish only in the case of .
The effective coupling between the Higgs boson and the gluon is, up to the linear order of ,
(18) |
Note that contribution of the composite resonances in the sectors (but not the light fermions ) are included in the above equation Liu et al. (2017).
Therefore, the Higgs effective couplings in the fermionic sectors are concretely
-
•
:
(19) -
•
and :
(20) -
•
:
(21) where is the ratio of to .
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 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.

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 ), twin weak gauge bosons from the gauge group in the twin sector, and twin gluons Craig et al. (2015a). The Higgs boson is identified as a pseudo Nambu-Goldstone boson from the coset , or equivalently , 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 . For reasonable values of , these twin gauge bosons can be sufficiently heavy. Only the other light degrees of freedom, lighter than , 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 , all the relevant Higgs couplings are as follows: i) for the couplings in the SM sector,
(22) |
ii) for the couplings in the twin sector,
(23) |
assuming the Yukawa couplings are the same, i.e. . 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 does not depend on . Furthermore, the masses of the twin fermions are also connected to the partners in the SM sector, e.g.,
(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 and can be larger or smaller than depending on the scale , which would result in different invisible widths of the Higgs boson.
The particles in the twin sector induce the invisible width as following:
(25) |
where for the -particle and denotes the branching ratio of the decay channel of in the SM.
Note that the running masses of the fermions and in the SM model are roughly GeV and 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 . 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.

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.