Di-higgs enhancement by neutral scalar as probe of new colored sector

TL;DR

This paper explores how an extra neutral scalar, produced via gluon fusion through new colored particles, can enhance Higgs pair production and potentially explain an observed excess at the LHC.

Contribution

It proposes a model where a neutral scalar enhances di-Higgs production via colored particles and analyzes its experimental constraints and future collider prospects.

Findings

Potential explanation for the ATLAS $b\bar{b}\gamma\gamma$ excess as a 300 GeV scalar resonance.

Enhanced di-Higgs production signals at future colliders.

Constraints from current LHC data on the model parameters.

Abstract

We study a class of models in which the Higgs pair production is enhanced at hadron colliders by an extra neutral scalar. The scalar particle is produced by the gluon fusion via a loop of new colored particles, and decays into di-Higgs through its mixing with the Standard Model Higgs. Such a colored particle can be the top/bottom partner, such as in the dilaton model, or a colored scalar which can be triplet, sextet, octet, etc., called leptoquark, diquark, coloron, etc., respectively. We examine the experimental constraints from the latest Large Hadron Collider (LHC) data, and discuss the future prospects of the LHC and the Future Circular Collider up to 100TeV. We also point out that the 2.4$\sigma$ excess in the $b\bar b\gamma\gamma$ final state reported by the ATLAS experiment can be interpreted as the resonance of the neutral scalar at 300GeV.

Tables

Table 1. Table 1: Colored particles that may run in the loop represented by the blob in Fig.  1 , and their possible parameters. We assume that they are S ​ U ​ ( 2 ) L 𝑆 𝑈 subscript 2 𝐿 SU(2)_{L} singlets. The electromagnetic charge Q 𝑄 Q is fixed to allow a mixing with either top or bottom quark for the Dirac spinor and a Yukawa coupling with a pair of SM fermions for the complex scalar; see Appendix  C . In the last row, F 𝐹 F stands for T 𝑇 T or B 𝐵 B .

	Dirac spinor			complex scalar
field	$T$	$B$	…	${\varphi }_{\mathrm{𝟑}}$	${\varphi }_{\mathrm{𝟔}}$	${\varphi }_{\mathrm{𝟖}}$	…
$SU{\left(3\right)}_C$	$\mathrm{𝟑}$	$\mathrm{𝟑}$	…	$\mathrm{𝟑}$	$\mathrm{𝟔}$	$\mathrm{𝟖}$	…
$Q$	$\frac{2}{3}$	$-\frac{1}{3}$	…	$-\frac{1}{3}$, $-\frac{4}{3}$	$\frac{1}{3}$, $-\frac{2}{3}$, $\frac{4}{3}$	0, $-1$	…
$\Delta b_g$	$\frac{2}{3}$	$\frac{2}{3}$	…	$\frac{1}{6}$	$\frac{5}{6}$	$1$	…
$\Delta b_{\gamma }$	$\frac{16}{9}$	$\frac{4}{9}$	…	$\frac{1}{9}$, $\frac{16}{9}$	$\frac{2}{9}$, $\frac{8}{9}$, $\frac{32}{9}$	0, $\frac{8}{3}$	…
$\eta$		$y_F{N}_F\frac{v}{M_F}$		${\kappa }_{\varphi }{N}_{\varphi }\frac{fv}{M_{\varphi }^2}$

Equations

$$H^0$$

$$S$$

$$⟨H^0⟩$$

$$\Delta \mathcal{L}$$

$$\Gamma \left(s\to hh\right)$$

$$\Gamma \left(s\to VV\right)$$

$$\Gamma \left(s\to t\bar{t}\right)$$

$$M_T$$

$${\mathcal{L}}_{\text{eff}}$$

$$\frac{1}{g_s^2}$$

$${\mathcal{L}}_{\text{eff}}^{hgg}$$

$${\mathcal{L}}_{\text{eff}}^{sgg}$$

$${\mathcal{L}}_{\text{eff}}^{h\gamma \gamma }$$

$${\mathcal{L}}_{\text{eff}}^{s\gamma \gamma }$$

$$\eta$$

$$M_{\varphi }^2$$

$${\mathcal{L}}_{\text{eff}}$$

$$\Gamma \left(s\to gg\right)$$

$$b_g$$

$$b_{\gamma }$$

$$\hat{\sigma }\left(gg\to s\right)$$

$${\sigma }_s$$

$$\sigma \left(pp\to s\right)$$

$$\frac{\text{d}{\mathcal{L}}^{gg}}{\text{d}\tau }$$

$${\hat{\sigma }}_{\text{LO}}^{\text{SM}}\left(gg\to hh\right)$$

$${\mu }_{hhh}$$

$$\frac{N_T{y}_T^2{\sin }^2\theta }{y_t^2/2}\frac{y_T^2{f}^2}{M_T^2}.$$

$$\frac{{\mu }_{\text{eff}}{M}_T}{32\pi v^2}{\sin }^3\theta$$

$${\mu }_{\text{ggF}}$$

$${\mu }_{h\to \gamma \gamma }$$

$${\mu }_{h\to gg}$$

$${\mu }_{h\to f\bar{f},WW,ZZ}$$

$$\left(\frac{{\Gamma }_h}{{\Gamma }_h^{\text{SM}}}\right)$$

$$M_T,M_B\gtrsim 800\text{GeV}.$$

Full text

1 Introduction

The di-Higgs production will continue to be one of the most important physics targets in the Large Hadron Collider (LHC) and beyond, since its observation leads to a measurement of the tri-Higgs coupling, and will provide a test if it matches with the Standard Model (SM) prediction [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. Since its production in the SM is destructively interfered by the top-quark box-diagram contribution, sizable production of di-Higgs directly implies a new physics signature [12].

It is important to examine in what kind of a model the di-Higgs signal is enhanced. Indeed the enhancement has been pointed out in the models with two Higgs doublets [13, 14, 15, 16, 17, 18, 19, 20, 21], type-II seesaw [22], light colored scalars [23], heavy quarks [24], effective operators [25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39], dilaton [40], strongly interacting light Higgs and minimal composite Higgs [41, 42, 43, 44], little Higgs [45, 46, 47], twin Higgs [48], Higgs portal interactions [40, 49, 50, 51, 52, 53, 54, 55, 56, 57], supersymmetric partners [58, 59, 5, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71], and Kaluza-Klein graviton [72]. Other related issues are discussed in Refs. [73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90]. The triple Higgs productions at the LHC and the future circular collider (FCC) are also discussed in Refs. [91, 92, 93].

In this paper, we study a class of models in which the di-Higgs process is enhanced by a resonant production of an extra neutral scalar particle. Its production is radiatively induced by the gluon fusion via a loop of new colored particles. Its tree-level decay is due to the mixing with the SM Higgs boson. As concrete examples of the new colored particle that can decay into SM ones in order not to spoil cosmology, we examine the top/bottom partner, such as in the dilaton model, and the colored scalar which are triplet (leptoquark), sextet (diquark), and octet (coloron).

We are also motivated by the anomalous result reported by the ATLAS Collaboration: the 2.4 $\sigma$ excess in the search of di-Higgs signal using $b\bar{b}$ and $\gamma\gamma$ final states with the $m_{(b\bar{b})(\gamma\gamma)}(=m_{hh})$ invariant mass at around 300 GeV [15]. The excess in $m_{(\gamma\gamma)}$ distribution is right at the SM Higgs mass on top of both the lower and higher mass-side-band background events. The requested signal cross section roughly corresponds to 90 times larger than what is expected in the SM. Thus the enhancement, if from new physics, should be dramatically generated via e.g. a new resonance at $300\,\text{GeV}$.

This paper is organized as follows. In Sec. 2, we present the model. In Sec. 3, we show how the di-Higgs event is enhanced. In Sec. 4, we examine the constraints on the model from the latest results from the ongoing LHC experiment. In Sec. 5, we present a possible explanation for the 2.4 $\sigma$ excess. In Sec. 6, we summarize our result and provide discussion. In Appendix A, we show how the effective interaction between the new scalar and Higgs is obtained from the original Lagrangian. In Appendix B, we give a parallel discussion for the $Z_{2}$ model. In Appendix C, we spell out the possible Yukawa interactions between the colored scalar and the SM fields.

2 Model

We consider a class of models in which the di-Higgs ($hh$) production is enhanced by the schematic diagram depicted in Fig. 1, where $s$ denotes the new neutral scalar and the blob generically represents an effective coupling of $s$ to the pair of gluons via the loop of the extra heavy colored particles. We assume that $h$ and $s$ are lighter and heavier mass eigenstates obtained from the mixing of the neutral component of the $SU(2)_{L}$-doublet $H$ and a real singlet $S$ that couples to the extra colored particles:

[TABLE]

where $\theta$ is the mixing angle and $v$ and $f$ denote the vacuum expectation values (VEVs):

[TABLE]

with $v\simeq 246\,\text{GeV}$ and $m_{h}=125\,\text{GeV}$. We phenomenologically parametrize the effective $shh$ interaction as

[TABLE]

where $\mu_{\text{eff}}$ is a real parameter of mass dimension unity, whose explicit form in terms of original Lagrangian parameters is given in Appendix A. We note that the parameter $\mu_{\text{eff}}$ is a purely phenomenological interface between the experiment and the underlying theory in order to allow a simpler phenomenological expression for the tree-level branching ratios; see Sec. 2.1. We note also that the $\theta$-dependent $\mu_{\text{eff}}\!\left(\theta\right)$ goes to a $\theta$-independent constant in the small mixing limit $\theta^{2}\ll 1$; see Appendix A for detailed discussion. In Sec. 4, it will indeed turn out that only the small, but non-zero, mixing region is allowed in order to be consistent with the signal-strength data of the 125 GeV Higgs at the LHC.

The extra colored particle that runs in the loop, which has been generically represented by the blob in Fig. 1, can be anything that couples to $S$. It should be sufficiently heavy to evade the LHC direct search and decay into SM particles in order not to affect the cosmological evolution. In this paper, we consider the following two possibilities: a Dirac fermion that mixes with either top or bottom quark and a scalar that decays via a new Yukawa interaction with the SM fermions. For simplicity, we assume that the new colored particles are singlet under the SU(2)L in both cases.

In Table 1, we list the colored particles of our consideration. The higher rank representations of $SU(3)_{C}$ for the colored scalars are terminated at $\boldsymbol{8}$ in order not to have too higher dimensional Yukawa operators.111 The ultraviolet completion of the higher dimensional operator requires other new colored particles. We assume that their contributions are subdominant. E.g. they do not contribute to the effective $ggs$ vertex if they do not have a direct coupling to $S$.

The triplet $\phi_{\boldsymbol{3}}$ is nothing but the leptoquark. It is worth noting that the leptoquark with $Y=-1/3$ may account for $R_{D^{(*)}}$, $R_{K}$, and $\left(g-2\right)_{\mu}$ anomalies simultaneously [94].

2.1 Tree-level decay

The scalar $s$ may dominantly decay into di-Higgs at the tree level due to the coupling (4):

[TABLE]

For $m_{s}>2m_{Z}$, the partial decay rate into the pair of vector bosons $s\to VV$ with $V=W,Z$ are

[TABLE]

where $\delta_{Z}=1$, $\delta_{W}=2$, and $x_{V}=m_{V}^{2}/m_{s}^{2}$; see e.g. Ref. [95]. Similarly for $m_{s}>2m_{t}$, the partial decay width into a top-quark pair is

[TABLE]

Note that the tree-level branching ratios become independent of $\theta$ thanks to the parametrization (4).

The total decay width $\Gamma_{\text{total}}$ is the sum of the above rates at the tree level. In the small mixing limit $\theta^{2}\ll 1$, the tree-level decay width becomes small and the loop level decay, which is described in Sec. 2.3, can be comparable to it. The diphoton constraint is severe in this parameter region, as will be discussed in Sec. 4.

In Fig. 2, we plot the tree-level branching ratios in the $\mu_{\text{eff}}$ vs $m_{s}$ plane. Note that the $\theta$-dependence drops out of the tree-level branching ratios when we use $\mu_{\text{eff}}$ as a phenomenological input parameter as in Eq. (4) because then all the decay channels have the same $\theta$ dependence $\propto\sin^{2}\theta$.

2.2 Effective coupling to photons and gluons

We first consider the vector-like top-partner $T$ as the colored particle running in the loop that is represented as the blob in Fig. 1. The bottom-partner $B$ can be treated in the same manner, as well as the colored scalars.

The mass of the top partner is given as

[TABLE]

where $m_{T}$ and $y_{T}$ are the vector-like mass of $T$ and the Yukawa coupling between $T$ and $S$, respectively. The top-partner $T$ mixes with the SM top quark. We note that limit $m_{T}\to 0$ corresponds to an effective dilaton model.222 The particular dilaton model in Ref. [96] corresponds to the identification of the lighter 125 GeV scalar to be an $S$-like one, contrary to this paper.

Given the kinetic term of gluon that is non-canonically normalized,

[TABLE]

the effective coupling after integrating out the top and $T$ can be obtained by the replacement $\left\langle S\right\rangle\to S$ and $\left\langle H^{0}\right\rangle\to H^{0}$ in the running coupling; see e.g. Refs. [97, 96]:

[TABLE]

where $b_{g}^{\text{top}}$ and $\Delta b_{g}$ are the contributions of top and $T$ to the beta function, respectively. To use this formula, we need to assume the new colored particles are slightly heavier than the neutral scalar. For a Dirac spinor in the fundamental representation, $b_{g}^{\text{top}}=\Delta b_{g}={1\over 2}\times{4\over 3}={2\over 3}$. The resultant effective interactions for the canonically normalized gauge fields are

[TABLE]

where $F_{\mu\nu}$ being the (canonically normalized) field strength tensor of the photon, $\alpha_{s}$ and $\alpha$ denoting the chromodynamic and electromagnetic fine structure constants, respectively, $N_{c}=3$, $b_{\gamma}^{\text{SM}}\simeq-6.5$ and

[TABLE]

with $N_{T}$ being the number of $T$ introduced. The values $\Delta b_{g}={1\over 2}\times{4\over 3}={2\over 3}$ and $\Delta b_{\gamma}=N_{c}Q_{T}^{2}\times{4\over 3}={16\over 9}$ are listed in Table 1.

The bottom partner $B$ can be treated exactly the same way. According to Table 1, $\Delta b_{\gamma}$ becomes one fourth compared to the above.

For the colored scalar $\phi$, its diagonal mass is given as

[TABLE]

where we have assumed the $Z_{2}$ symmetry $S\to-S$ for simplicity; $m_{\phi}$ is the original diagonal mass in the Lagrangian; and $\kappa_{\phi}$ is the quartic coupling between $S$ and $\phi$.333 The three point interaction between the neutral and the colored scalar can be introduced. If the sign of the three and the four point couplings are opposite, $\eta$ can be enhanced in some parameter region.

The possible values of the electromagnetic of $\phi$ are $Q=-1/3$ and $-4/3$ for the leptoquark $\phi_{\boldsymbol{3}}$; $Q=1/3$, $-2/3$, and $4/3$ for the color-sextet $\phi_{\boldsymbol{6}}$; and $Q=0$ and $-1$ for the color-octet $\phi_{\boldsymbol{8}}$; see Appendix C. Correspondingly the values of $\Delta b_{g}$ are ${1\over 2}\times{1\over 3}={1\over 6}$, $\left({N_{c}\over 2}+1\right)\times{1\over 3}={5\over 6}$, and $N_{c}\times{1\over 3}=1$, and $\Delta b_{\gamma}$ are $Q^{2}$, $2Q^{2}$, and ${8\over 3}Q^{2}$. Again the effective interactions are obtained as in Eqs. (11)–(14) from the replacement (10) with the substitution $y_{T}/M_{T}\to\kappa_{\phi}f/M_{\phi}^{2}$, where $f$ has been the VEV of $S$; see Eq. (3). Note that the expression for $\eta$ is now $\eta=\kappa_{\phi}N_{\phi}{fv/M_{\phi}^{2}}$, where $N_{\phi}$ is the number of $\phi$ introduced. We list all these parameters in Table 1.

2.3 Loop-level decay

No direct contact to the gauge bosons are allowed for the singlet scalar $S$, and the tree-level decay of $s$ into a pair of gauge bosons is only via the mixing with the SM Higgs boson. Therefore the decay of $s$ to $gg$ and $\gamma\gamma$ are only radiatively generated. Given the effective operators from the loop of heavy colored particle

[TABLE]

the partial decay widths are

[TABLE]

where the factor 8 difference comes from the number of degrees of freedom of gluons in the final state. Concretely,

[TABLE]

If we go beyond the scope of this paper and allow the particles in the loop to be charged under $SU(2)_{L}$, then the loop contribution to the decay channels to $Z\gamma$, $ZZ$ and $W^{+}W^{-}$ might also become significant; see e.g. Ref. [98].

3 Production of singlet scalar at hadron colliders

We calculate the production cross section of $s$ via the gluon fusion with the narrow width approximation:444 The colored particles running in the blob in Fig. 1 might also have a direct coupling with the quarks in the proton, and possibly change the production cross section of $s$ if it is extremely large. In this paper we assume that this is not the case.

[TABLE]

where

[TABLE]

Therefore, we reach the expression with the gluon parton distribution function (PDF) for the proton $g(x,\mu_{F})$:

[TABLE]

where $\tau:=m_{s}^{2}/s$ and

[TABLE]

is the luminosity function, in which the factorization scale $\mu_{F}$ is taken to be $\mu_{F}=\sqrt{\tau s}$.555 Notational abuse of $s$ for the singlet scalar field and for the Mandelstam variable of $pp$ scattering should be understood.

Using the leading order CTEQ6L [99] PDF, we plot in Fig. 3 the production cross section $\sigma\!\left(pp\to s\right)$ as a function of $m_{s}$ for a phenomenological benchmark setting $\left|b_{g}\right|={\Delta b_{g}\over 2}{v\over m_{s}}$ with $\Delta b_{g}={2\over 3}$ (top/bottom partner). Other particles just scale as $\sigma\!\left(pp\to s\right)\propto\left(\Delta b_{g}\right)^{2}$. The value $\sqrt{s}=14\,\text{TeV}$ is motivated by the High-Luminosity LHC; 28 TeV and 33 TeV by the High-Energy LHC (HE-LHC); and 75 TeV, and 100 TeV by the Future Circular Collider (FCC) [100, 101, 102].

We see that typically the top/bottom partner models give a cross section $\sigma\!\left(pp\to s\right)\gtrsim 1\,\text{fb}$, which could be accessed by a luminosity of $\mathcal{O}\!\left(\text{ab}^{-1}\right)$, for the scalar mass $m_{s}\lesssim 1.3\,\text{TeV}$, 2 TeV, and 4 TeV at the LHC, HE-LHC, and FCC, respectively.

Several comments are in order:

• •

Our setting corresponds to putting $M_{T}=y_{T}N_{T}m_{s}$ in Eq. (15) in order to reflect the naive scaling of $\eta\sim v/f$ with $f\sim m_{s}$; recall that we need $M_{T}\gtrsim m_{s}$ to justify integrating out the top partner to write down the effective interactions (11)–(14).

• •

Here we have used the leading order parton distribution function. The higher order corrections may be approximated by multiplying an overall factor $K$, the so-called $K$-factor, which takes value $K\simeq 1.6$ for the SM Higgs production at LHC; see e.g. Ref. [95].

• •

The SM cross section for $pp\to hh$ is of the order of 10 fb and $10^{3}\,\text{fb}$ for $\sqrt{s}=8\,\text{TeV}$ and 100 TeV, respectively [12]. We are interested in the on-shell production of $s$, and the non-resonant SM background can be discriminated by kinematical cuts. The detailed study is beyond the scope of this paper and will be presented elsewhere.

• •

When we consider the new resonance with a narrow width (21), we can neglect the box contribution from the extra colored particles as the box contribution gets a suppression factor666 In the SM, the $gg\to hh$ cross section takes the following form at the leading order [12]:

$\displaystyle\hat{\sigma}^{\text{SM}}_{\text{LO}}\!\left(gg\to hh\right)$ $\displaystyle=\int_{\hat{t}_{-}}^{\hat{t}_{+}}\text{d}\hat{t}{G_{\text{F}}^{2}\alpha_{s}^{2}\over 256\left(2\pi\right)^{3}}\left[\left|{\mu_{hhh}v\over\left(\hat{s}-m_{h}^{2}\right)+im_{h}\Gamma_{h}}F^{\text{SM}}_{\bigtriangleup}+F^{\text{SM}}_{\Box}\right|^{2}+\left|G^{\text{SM}}_{\Box}\right|^{2}\right],$

where $G_{\text{F}}$ is the Fermi constant; $\mu_{hhh}=3m_{h}^{2}/v$ is the $hhh$ coupling in the SM; and $F^{\text{SM}}_{\bigtriangleup}$, $F^{\text{SM}}_{\Box}$, and $G^{\text{SM}}_{\Box}$ are the triangular and box form factors, approaching $F^{\text{SM}}_{\bigtriangleup}\to 2/3$, $F^{\text{SM}}_{\Box}\to-2/3$, and $G^{\text{SM}}_{\Box}\to 0$ in the large top-quark-mass limit. A large cancellation takes place between $F^{\text{SM}}_{\bigtriangleup}$ and $F^{\text{SM}}_{\Box}$ as is well known. For the on-shell resonance production of $s$, on the other hand, the triangle contribution from the fermion loops dominates over the box loop contribution: The new triangle contribution for $s$ can be well approximated by replacing the expression for the SM as

$\displaystyle\mu_{hhh}$ $\displaystyle\to\mu_{\rm eff}\sin\theta,$ $\displaystyle\,m_{h}$ $\displaystyle\to m_{s},$ $\displaystyle\Gamma_{h}$ $\displaystyle\to\Gamma_{s},$ $\displaystyle F^{\text{SM}}_{\bigtriangleup}$ $\displaystyle\to\Delta b_{g}\eta\cos\theta+b_{g}^{\rm top}\sin\theta,$

and the new box contribution of the top partner can be obtained from that of the SM-top quark with the multiplicative factor

$\displaystyle{N_{T}y_{T}^{2}\sin^{2}\theta\over y_{t}^{2}/2}\,{y_{T}^{2}f^{2}\over M_{T}^{2}}.$

Finally, taking the ratio of the size of the box contribution and the triangle contribution with $\Delta b_{g}=2/3$ and $\eta=y_{T}N_{T}v/M_{T}\sim N_{T}v/M_{T}$, $y_{T}\sim y_{t}$, and $m_{s}\Gamma_{s}\sim\mu_{\text{eff}}^{2}\sin^{2}\theta/32\pi$, we get the result in Eq. (25).

[TABLE]

4 LHC constraints

We examine LHC constraints on the model for various $m_{s}$. That is, we verify constraints from 125 GeV Higgs signal strength, from $s\to ZZ\to 4l$ search, from $s\to\gamma\gamma$ search, and from the direct search of the colored particles running in the blob in Fig. 1.

4.1 Bound from Higgs signal strength

We first examine the bound on $\theta$ and $\eta$ from the Higgs signal strengths in various channels. The “partial signal strength” for the Higgs production becomes

[TABLE]

where ggF, VBF, VH, and ttH are the gluon fusion, vector-boson fusion, associated production with vector, and that with a pair of top quarks, respectively; see e.g. Ref. [103] for details. Similarly, the partial signal strength for the Higgs decay is

[TABLE]

where the ratio of the total widths is given by

[TABLE]

with $\text{Br}_{h\to\text{SM others}}^{\text{SM}}=0.913$, $\text{Br}_{h\to\gamma\gamma}^{\text{SM}}=0.002$ and $\text{Br}_{h\to gg}^{\text{SM}}=0.085$. We compare these values with the corresponding constraints given in Ref. [103]. Results are shown in Fig. 4 for the matter contents summarized in Table 1. We note that the region near $\theta\simeq 0$ is always allowed by the signal strength constraints, though it is excluded by the di-photon search as we will see.

4.2 Bound from $s\to ZZ\to 4l$

One of the strongest constraints on the model comes from the heavy Higgs search in the four lepton final state at $\sqrt{s}=13\,\text{TeV}$ at ATLAS [104]. Experimentally, an upper bound is put on the cross section $\sigma\!\left(pp\to s\to ZZ\to 4l\right)$, with $l=e,\mu$, for each $m_{s}$. Its theoretical cross section is obtained by multiplying the production cross section (23) by the branching ratio $\operatorname{BR}\!\left(s\to ZZ\right)=\Gamma\!\left(s\to ZZ\right)/\Gamma\!\left(s\to\text{all}\right)$ and $\left(\operatorname{BR}_{\text{SM}}\!\left(Z\to ee,\mu\mu\right)\right)^{2}\simeq\left(6.73\%\right)^{2}$; see Sec. 2.1.

In Fig. 6, we plot $2\sigma$ excluded regions on the $\mu_{\text{eff}}$ vs $m_{s}$ plane with varying $b_{g}$ from 0 to 1 with incrementation 0.2. The weakest bound starts to exist on the plane from $b_{g}=0.2$. $K$-factor is set to be $K=1.6$. The experimental bound becomes milder for large $\mu_{\text{eff}}$ because the di-Higgs channel dominate the decay of the neutral scalar. The large fluctuation of the bound is due to the statistical fluctuation of the original experimental constraint.

We note that though we have focused on the strongest constraint at the low $m_{s}$ region, the other decay channels of $WW\to l\nu qq$ and of $ZZ\to\nu\nu qq$ and $ll\nu\nu$ may also become significant at the high mass region $m_{s}\gtrsim 700\,\text{GeV}$.

4.3 Bound from $s\to\gamma\gamma$

A strong constraint comes from the heavy Higgs search in the di-photon final state at $\sqrt{s}=13\,\text{TeV}$ at ATLAS [105]. Experimentally, an upper bound is put on the cross section $\sigma\!\left(pp\to s\to\gamma\gamma\right)$ for each $m_{s}$. Its theoretical cross section is obtained by multiplying the production cross section (23) by the branching ratio $\operatorname{BR}\!\left(s\to\gamma\gamma\right)=\Gamma\!\left(s\to\gamma\gamma\right)/\Gamma\!\left(s\to\text{all}\right)$; see Sec. 2.1. Since this constraint is strong in the small mixing region, where the loop-level decay is comparable to the tree-level decay, we include the loop-level decay channels into $\Gamma\!\left(s\to\text{all}\right)$ for this analysis; see Sec. 2.3.

In Fig 7, we plot the $2\sigma$-excluded regions on $\mu_{\text{eff}}$ vs $m_{s}$ plane for $\sin\theta=$0.01, 0.03, 0.05, and 0.1, with varying $b_{g}b_{\gamma}$ from 0 to 2 with incrementation 0.2. $K$-factor is set to be $K=1.6$. If $\sin\theta=0.01$, broad region is excluded for $b_{g}b_{\gamma}=0.4$. On the other hand, the experimental bound is negligibly weak in the case of $\sin\theta=0.1$. The large fluctuation of the bound is due to the statistical fluctuation of the original experimental constraint.

In Fig. 8, we plot the same $2\sigma$-excluded regions on the $\sin\theta$ vs $\eta$ plane for $m_{s}=300\,\text{GeV}$, $600\,\text{GeV}$, $900\,\text{GeV}$, $1200\,\text{GeV}$, and $1500\,\text{GeV}$. In the left and right panels, we set $\mu_{\text{eff}}=1\,\text{TeV}$ and $\mu_{\text{eff}}=\sqrt{3}m_{s}^{2}/v$. The latter corresponds to $\Gamma\!\left(s\to hh\right)=\sum_{V=W,Z}\Gamma\!\left(s\to VV\right)$ which is chosen such that there are sizable di-Higgs event and that $\mu_{\text{eff}}$ is not too large. $K$-factor is set to be $K=1.6$. We emphasize that the small mixing limit $\sin\theta\to 0$ is always excluded by the di-photon channel in contrast to the other bounds, though it cannot be seen in Fig. 8 in the small $\eta$ region due to the resolution.

The bound from $s\to Z\gamma$ is weaker and we do not present the result here.

4.4 Bound from direct search for colored particles

We first review the mass bound on the extra colored particles. For the $SU(2)_{L}$ singlet $T$ and $B$ [106, 107],

[TABLE]

The mass bound for the leptoquark $\phi_{\boldsymbol{3}}$, diquark $\phi_{\boldsymbol{6}}$, and coloron $\phi_{\boldsymbol{8}}$ are given in Refs. [108, 109], [110], and [111] as

[TABLE]

respectively, depending on the possible decay channels.

For the top-partner $M_{T}\gtrsim 800\,\text{GeV}$ with $\theta\simeq 0$, we get $\eta\lesssim 0.3y_{T}N_{T}$. Therefore, we need rather large Yukawa coupling $y_{T}\simeq 2.2$ for $N_{T}=1$ in order to account for Eq. (33) by Eq. (35).777 Strictly speaking, the bound on $M_{T}$ slightly changes when $N_{T}\geq 2$, and hence the bound for $y_{T}N_{T}$ could be modified accordingly.

The same argument applies for the bottom partner since it has the same $\Delta b_{g}=2/3$.

Similarly for a colored scalar with $M_{\phi}\gtrsim 0.7$, 1.1, 5.5, and 7 TeV, we get $\eta\lesssim\kappa_{\phi}N_{\phi}{f\over 2\,\text{TeV}}$, $\kappa_{\phi}N_{\phi}{f\over 4.9\,\text{TeV}}$, $\kappa_{\phi}N_{\phi}{f\over 123\,\text{TeV}}$, and $\kappa_{\phi}N_{\phi}{f\over 200\,\text{TeV}}$, respectively. For $\theta\simeq 0$, the value of $b_{g}$ is suppressed or enhanced by extra factors ${1\over 6}/{2\over 3}=1/4$, ${5\over 6}/{2\over 3}={5\over 4}$, and $1/{2\over 3}=3/2$, respectively, compared to the top partner. Therefore, from Eq. (36), we need $\kappa_{\phi}N_{\phi}f\gtrsim 5$–13 TeV, 106 TeV, and 54 TeV for $\phi_{\boldsymbol{3}}$, $\phi_{\boldsymbol{6}}$, and $\phi_{\boldsymbol{8}}$, respectively, in order to account for the 2.4$\sigma$ excess at $\theta^{2}\ll 1$.

5 Accounting for 2.4$\sigma$ excess of $b\bar{b}\gamma\gamma$ by $m_{s}=300\,\text{GeV}$

It has been reported by the ATLAS Collaboration that there exist 2.4$\sigma$ excess of $hh$-like events in the $b\bar{b}\gamma\gamma$ final state [15]. This corresponds to the extra contribution to the SM cross section888 At $\sqrt{s}=8\,\text{TeV}$, $\sigma_{\text{SM}}\!\left(pp\to hh\right)=9.2\,\text{fb}$. The expected number of events are $1.3\pm 0.5$, $0.17\pm 0.04$, and $0.04$ for the non-$h$ background, single $h$, and the SM $hh$ events, respectively. Since the observed number of events is 5, excess is $5-1.3-0.17=3.5$, which is $3.5/0.04=87.5$ times larger than the SM $hh$ events. Therefore, the excess corresponds to $9.2\,\text{fb}\times 87.5=0.8\,\text{pb}$.

[TABLE]

In Fig. 9, we plot the branching ratio at $m_{h}=300\,\text{GeV}$ as a function of $\mu_{\text{eff}}$.

5.1 Signal

With $m_{s}=300\,\text{GeV}$, we get the luminosity functions

[TABLE]

That is,

[TABLE]

In Fig. 10, we plot the preferred contour to explain the 2.4$\sigma$ excess at $m_{s}=300\,\text{GeV}$, where the shaded region is excluded at the 95% C.L. by the $\sigma\!\left(pp\to s\to ZZ\to 4l\right)_{13\,\text{TeV}}$ constraint that has been discussed in Sec. 4.2. We have assumed the $K$-factor $K=1.6$.

We see that at the benchmark point $\theta\simeq 0$, the lowest and highest possible values of $\mu_{\text{eff}}$ and $\eta$ are, respectively,

[TABLE]

in order to account for the cross section (33). The ratio of the upper bound on $\eta$ is given by the scaling $\propto\left(\Delta b_{g}\right)^{2}$.

5.2 Constraints

When $m_{s}=300\,\text{GeV}$, the $95\%$ C.L. upper bound at $\sqrt{s}=13\,\text{TeV}$ is $\sigma\!\left(s(\text{ggF})\to ZZ\to 4l\right)_{13\,\text{TeV}}\lesssim 0.8\,\text{fb}$ [104]; see also Fig. 6. The corresponding excluded region is plotted in Fig. 10.

Currently, the strongest direct constraint on the di-Higgs resonance at $m_{s}=300\,\text{GeV}$ comes from the $\sqrt{s}=8\,\text{TeV}$ data in the $b\bar{b}\gamma\gamma$ final state at CMS [112] and in $b\bar{b}\tau\tau$ at ATLAS [113]:

[TABLE]

at the 95% C.L. The preferred value (33) is still within this limit.

We note that the current limit for the $m_{s}=300\,\text{GeV}$ resonance search at $\sqrt{s}=13\,\text{TeV}$ is from the $b\bar{b}\gamma\gamma$ final state at ATLAS [113] and from $b\bar{b}bb$ at CMS [112]:

[TABLE]

at the 95% C.L. This translates to the $\sqrt{s}=8\,\text{TeV}$ cross section:

[TABLE]

This is weaker than the direct 8 TeV bound (37).

The branching ratio for $s\to\gamma\gamma$ is999 The power of $m_{s}$ dependence is valid in the limit $m_{s}\gg 2m_{h}$.

[TABLE]

We see that the loop suppressed decay into diphoton is negligible compared to the tree-level decay via the interaction (4). For $m_{s}=300$ GeV, the cross section at $\sqrt{s}=13\,\text{TeV}$ is

[TABLE]

We see that the loop-suppressed $\Gamma\!\left(s\to\gamma\gamma\right)$ becomes the same order as $\Gamma\!\left(s\to hh\right)$ when $\theta\lesssim 10^{-3}$ and that the region $\theta\lesssim 10^{-2}$ is excluded by the diphoton search, $\sigma\!\left(pp\to s\to\gamma\gamma\right)_{\text{13\,\text{TeV}}}\lesssim 10\,\text{fb}$ [105], for a typical set of parameters that explains the 300 GeV excess; see also Sec. 4.3.

We comment on the case where the neutral scalar is charged under the $Z_{2}$ symmetry, $S\to-S$, or is extended to a complex scalar charged under an extra U(1), $S\to e^{i\varphi}S$. In such a model, the effective coupling in the small mixing limit becomes

[TABLE]

see Appendix B. That is, for a given $m_{s}$, there is an upper bound on the product $\mu_{\text{eff}}\,\eta$: $\mu_{\text{eff}}\,\eta\lesssim m_{s}$. On the other hand, the production cross section and the di-Higgs decay rate of $s$ are proportional to $\eta^{2}$ and $\mu_{\text{eff}}^{2}$, and hence there is a preferred value of $\mu_{\text{eff}}\,\eta$ in order to account for the 2.4 $\sigma$ excess by $m_{s}=300\,\text{GeV}$; see Fig. 10. In the $Z_{2}$ model and the $U(1)$ model, this preferred value exceeds the above upper bound. That is, they cannot account for the excess. More rigorous proof can be found in Appendix B.

On the other hand, a singlet scalar that does not respect additional symmetry does not obey this relation (42). Because of this reason, a singlet scalar without $Z_{2}$ symmetry is advantageous to enhance the di-Higgs signal in general, and can explain the excess by $m_{s}=300\,\text{GeV}$.

6 Summary and discussion

We have studied a class of models in which the di-Higgs production is enhanced by the $s$-channel resonance of the neutral scalar that couples to a pair of gluons by the loop of heavy colored fermion or scalar. As such a colored particle, we have considered two types of possibilities:

• •

the vector-like fermionic partner of top or bottom quark, with which the neutral scalar may be identified as the dilaton in the quasi-conformal sector,

• •

the colored scalar which is either triplet (leptoquark), sextet (diquark), or octet (coloron).

We have presented the future prospect for the enhanced di-Higgs production in the LHC and beyond. Typically, the top/bottom partner models give a cross section $\sigma\!\left(pp\to s\right)\gtrsim 1\,\text{fb}$, which could be accessed by a luminosity of $\mathcal{O}\!\left(\text{ab}^{-1}\right)$, for the scalar mass $m_{s}\lesssim 1.3\,\text{TeV}$, 2 TeV, and 4 TeV at the LHC, HE-LHC, and FCC, respectively.

We have examined the constraints from the direct searches for the di-Higgs signal and for a heavy colored particle, as well as the Higgs signal strengths in various production and decay channels. Typically small and large mixing regions are excluded by the diphoton resonance search and by the Higgs signal strength bounds, respectively. Region of small $\mu_{\text{eff}}$ is excluded by the diphoton search as well as by the $s\to ZZ\to 4l$ channel.

We also show a possible explanation of the 2.4$\sigma$ excess of the di-Higgs signal in the $b\bar{b}\gamma\gamma$ final state, reported by the ATLAS experiment. We have shown that the $Z_{2}$ model explained in Appendix B cannot account for the excess, while the general model in Appendix A can. A typical benchmark point which evades all the bounds and can explain the excess is

[TABLE]

For the top/bottom partner $T,B$, the required value to explain the 2.4$\sigma$ excess for the Yukawa coupling is rather large $y_{F}N_{F}\gtrsim 2.2$, where $N_{F}$ is the number of $F=T,B$ introduced. For the colored scalar $\phi$, required value of the neutral scalar VEV, $f=\left\langle S\right\rangle$, are

[TABLE]

where $\kappa_{\phi}$ and $N_{\phi}$ are the quartic coupling between the colored and neutral scalars and the number of colored scalar introduced, respectively.

In this paper, we have restricted ourselves to the case where the colored particle running in the blog in Fig. 1 are $SU(2)_{L}$ singlet. Cases for doublet, triplet, etc., which could be richer in phenomenology, will be presented elsewhere. We have assumed $M_{F},M_{\phi}\gtrsim m_{s}$ to justify integrating out the colored particle. It would be worth including loop functions to extend the region of study toward $M_{F},M_{\phi}\lesssim m_{s}$. A full collider simulation of this model for HL-LHC and FCC would be worth studying. A theoretical background of this type of the neutral scalar assisted by the colored fermion/scalar is worth pushing, such as the dilaton model and the leptoquark model with spontaneous $B-L$ symmetry breaking.

Acknowledgement

The work of K. Nakamura and K.O. are partially supported by the JSPS KAKENHI Grant No. 26800156 and Nos. 15K05053, 23104009, respectively. S.C.P. and Y.Y. are supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP) (No. 2016R1A2B2016112). K. Nishiwaki would like to thank David London and the Group of Particle Physics of Université de Montréal for the kind hospitality at the final stage of this work.

Appendix

Appendix A General scalar potential

We write down the most general renormalizable potential including the SM Higgs $H$ and the singlet $S$:

[TABLE]

with101010 Our $\lambda_{H}$ differs from the conventionally used $\lambda$ by $\lambda_{H}=2\lambda$, with $\lambda=m_{h}^{2}/2v^{2}\simeq 0.13$ in the SM; see e.g. Refs. [114, 115].

[TABLE]

where $m_{S}^{2}$ and $m_{H}^{2}$ are (potentially negative) mass-squared parameters; $\lambda_{S}$, $\kappa$, $\lambda_{H}$ are dimensionless constants; $\mu_{S}$ and $\mu$ are real parameters of the mass dimension unity; and the tadpole term of $S$ is removed by the field redefinition $S\to S+\text{const.}$ The $Z_{2}$ model corresponds to setting $\mu_{S}=\mu=0$ which are prohibited by the $Z_{2}$ symmetry: $S\to-S$.

The vacuum condition reads

[TABLE]

Using this vacuum condition, and putting Eqs. (1) and (2), we can always rewrite $m_{H}^{2}$ and $m_{S}^{2}$ in terms of $v$, $f$, and other parameters. The mixing angle can be written as

[TABLE]

Now the effective coupling in Eq. (4) is written as

[TABLE]

In the small mixing limit $\theta^{2}\ll 1$, we obtain

[TABLE]

We note that the first term also goes to a constant for fixed $v,f$ because of Eq. (51): $\left(\kappa f+\mu\right)\propto\theta$. More explicitly,

[TABLE]

as $\theta\to 0$. That is, the $shh$ coupling vanishes in the small mixing limit: $\mu_{\text{eff}}\sin\theta\to 0$. Let us emphasize that this is a general feature since the $shh$ coupling necessarily requires the non-zero mixing term $v\,sh$ that is obtained by the replacement $h\to v$. In order to take this feature into account, we have parametrized the effective coupling as in Eq. (4).

The mass eigenvalues satisfy the relations,

[TABLE]

where we suppose $m_{s}>m_{h}\simeq 125\,\text{GeV}$. The tachyon free condition is that the right hand sides of Eqs. (55) and (56) are positive. Also, from the condition that the quartic terms are positive in the large field limit for any linear combination of two fields, we obtain $\lambda_{S}>0$, $\lambda_{H}>0$, and ${\left({\lambda_{H}\lambda_{S}\over 3}-\kappa^{2}\right)\left({\lambda_{H}\over 2}+{\lambda_{S}\over 3!}-\kappa\right)}>0$.111111 When we allow higher dimensional operators such as $S^{6}$, this vacuum stability condition can be violated. In this analysis, we restrict ourselves within the potential up to quartic order terms, and assumes that this condition is met.

In the model without the $Z_{2}$ symmetry, we can remove the parameters $\mu$ and $\mu_{S}$ using Eqs. (55) and (56). Then the mixing angle (51) may be rewritten as

[TABLE]

Such a solution for $\lambda_{H}>0$ exists when only when

[TABLE]

We see that the small mixing limit corresponds to $\lambda_{H}\searrow m_{h}^{2}/v^{2}$. Also one may remove $\mu,\mu_{S}$ from the small mixing limit (54):

[TABLE]

where we used Eqs. (55) and (56) in the first step, and substitued the $\lambda_{H}\searrow m_{h}^{2}/v^{2}$ limit in the next step. We see that the Higgs-singlet mixing $\kappa$ remains to be a free parameters even in the small mixing limit.

If we want to explain the $b\bar{b}\gamma\gamma$ excess [15], we set $m_{s}\simeq 300\,\text{GeV}$, and get

[TABLE]

Even in the small mixing limit, $\mu_{\text{eff}}$ in Eq. (59) with $m_{s}=300\,\text{GeV}$ can be as large as $\mu_{\text{eff}}\simeq 1\,\text{TeV}$ ($2\,\text{TeV}$) for $\kappa=-1$ ($-3$), which is well within the current experimental bound; see Fig. 10; if we are happy with an extremely large value, say, $\kappa=-4\pi$, we may push it up to $\mu_{\text{eff}}\simeq 6.7\,\text{TeV}$.

Appendix B $Z_{2}$ model

We consider the $Z_{2}$ model with $\mu=\mu_{S}=0$. The discussion is parallel to Appendix A. The mixing angle reads

[TABLE]

Especially in the limit $v\ll f$, we get $\tan 2\theta\to{6\kappa\over\lambda_{S}}{v\over f}$. Eqs. (55) and (56) may be solved e.g. as

[TABLE]

For $\lambda_{H}>0$, the solution with $m_{s}>m_{h}>0$ again exists when and only when the condition (58) is met. This condition also ensures $\lambda_{S}$ to be positive. Putting Eq. (62) into Eq. (61), we again obtain Eq. (57).

Finally, the small mixing limit of the effective coupling becomes

[TABLE]

If we want to set $m_{s}=300\,\text{GeV}$, we get $\mu_{\text{eff}}\simeq 490\,\text{GeV}$ in the small mixing limit $\theta^{2}\ll 1$, which is already excluded by the $s\to ZZ\to 4l$ search; see Fig. 10. The $Z_{2}$ model cannot explain the 2.4 $\sigma$ excess reported by ATLAS. For larger values of $m_{s}$, the $Z_{2}$ model is still viable.

Appendix C Yukawa interaction between colored scalar and SM particles

For the scalar in the fundamental representation $\phi_{\boldsymbol{3}}$, the possible Yukawa interactions are

[TABLE]

depending on the hypercharge of $\phi_{\boldsymbol{3}}$: $-1/3$, $-1/3$, and $-4/3$, respectively. The superscript c denotes the charge conjugation.

We note that we can in principle write down the following diquark interactions:

[TABLE]

depending on the hypercharge of $\phi_{\boldsymbol{3}}$: $-1/3$, $-4/3$, $2/3$ and $-1/3$, respectively, where $a,b,c$ and $i,j$ represent the indices of the $SU(3)_{C}$ and $SU(2)_{L}$ fundamental representations, respectively, and $\epsilon$ is the totally antisymmetric tensor. The coexistence of the leptoquark and the diquark interactions leads to rapid proton decay. Since the diquark interactions are strongly restricted compared with the leptoquark in direct searches in hadron colliders, we focus on the situation that only the leptoquark interactions are switched on. The diquark interactions can be forbidden e.g. by the $B-L$ symmetry.

For the symmetric scalar $\phi_{\boldsymbol{6}}$, a possible Yukawa is either one of

[TABLE]

depending on the hypercharge of $\phi_{\boldsymbol{6}}$: $4/3$, $1/3$, $-2/3$, and $1/3$, respectively.

For adjoint scalar, a possible lowest-dimensional Yukawa is either one of

[TABLE]

depending on the hypercharge of $\phi_{\boldsymbol{8}}$: [math], $-1$, $-1$, and [math], respectively, where we have assigned $Y_{H}=+1/2$ and $\Lambda$ denotes an ultraviolet cutoff scale.