Impact of CP-violating interference effects on MSSM Higgs searches

TL;DR

This paper investigates how CP-violating interference and mixing effects between neutral Higgs bosons in the MSSM with complex parameters significantly influence LHC search interpretations, especially near degenerate states.

Contribution

It demonstrates the importance of including interference effects in MSSM Higgs searches, revealing that neglecting them can lead to incorrect exclusion regions.

Findings

Interference effects can be strongly destructive, reducing observable signals.

Mixing increases individual Higgs production cross sections.

Parameter regions previously considered excluded remain viable when interference is considered.

Abstract

Interference and mixing effects between neutral Higgs bosons in the MSSM with complex parameters are shown to have a significant impact on the interpretation of LHC searches for additional Higgs bosons. Complex MSSM parameters introduce mixing between the CP-even and CP-odd Higgs states and generate CP-violating interference terms. Both effects are enhanced in the case of almost degenerate states. Employing as an example an extension of a frequently used benchmark scenario by a non-zero phase $\phi_{A_t}$, the interference contributions are obtained for the production of neutral Higgs bosons in gluon-fusion and in association with b-quarks followed by the decay into a pair of $\tau$-leptons. While the resonant mixing increases the individual cross sections for the two heavy Higgs bosons $h_2$ and $h_3$, strongly destructive interference effects between the contributions involving $h_2$ and $h_3$ leave a considerable parameter region unexcluded that would appear to be ruled out if the interference effects were neglected.

Figures

Equations

$$M_{\tilde{f}}^2$$

$$X=\left(\begin{array}{cc}{M}_2& \sqrt{2}{M}_W{s}_{\beta }\\ \sqrt{2}{M}_W{c}_{\beta }& \mu \end{array}\right).$$

$$Y=\left(\begin{array}{cccc}{M}_1& 0& -M_Z{c}_{\beta }{s}_W& M_Z{s}_{\beta }{s}_W\\ 0& M_2& M_Z{c}_{\beta }{c}_W& -M_Z{s}_{\beta }{c}_W\\ -M_Z{c}_{\beta }{s}_W& M_Z{c}_{\beta }{c}_W& 0& -\mu \\ M_Z{s}_{\beta }{s}_W& -M_Z{s}_{\beta }{c}_W& -\mu & 0\end{array}\right).$$

$${\mathcal{H}}_1$$

$${\mathcal{H}}_2$$

$$\left(\begin{array}{c}h\\ H\\ A\\ G\end{array}\right)=\left(\begin{array}{cccc}-s_{\alpha }& c_{\alpha }& 0& 0\\ c_{\alpha }& s_{\alpha }& 0& 0\\ 0& 0& -s_{{\beta }_n}& c_{{\beta }_n}\\ 0& 0& c_{{\beta }_n}& s_{{\beta }_n}\end{array}\right)\left(\begin{array}{c}{\varphi }_1^0\\ {\varphi }_2^0\\ {\chi }_1^0\\ {\chi }_2^0\end{array}\right),\left(\begin{array}{c}{H}^{\pm }\\ G^{\pm }\end{array}\right)=\left(\begin{array}{cc}-s_{{\beta }_c}& c_{{\beta }_c}\\ c_{{\beta }_c}& s_{{\beta }_c}\end{array}\right)\left(\begin{array}{c}{\varphi }_1^{\pm }\\ {\varphi }_2^{\pm }\end{array}\right),$$

$$g_{i,f_d{\bar{f}}_d}^{\mathrm{tree}}=-\frac{ig{m}_{f_d}}{2M_W}\cdot \left\{-\frac{s_{\alpha }}{c_{\beta }},\frac{c_{\alpha }}{c_{\beta }},i{\gamma }_5\tan \beta \right\},i=h,H,A.$$

$$\text{M}(p^2)=\left(\begin{array}{ccc}{m}_h^2-{\hat{\Sigma }}_{hh}(p^2)& -{\hat{\Sigma }}_{hH}(p^2)& -{\hat{\Sigma }}_{hA}(p^2)\\ -{\hat{\Sigma }}_{Hh}(p^2)& m_H^2-{\hat{\Sigma }}_{HH}(p^2)& -{\hat{\Sigma }}_{HA}(p^2)\\ -{\hat{\Sigma }}_{Ah}(p^2)& -{\hat{\Sigma }}_{AH}(p^2)& m_A^2-{\hat{\Sigma }}_{AA}(p^2)\end{array}\right).$$

$${\hat{\Gamma }}_{ij}(p^2)=i\left[(p^2-m_i^2){\delta }_{ij}+{\hat{\Sigma }}_{ij}(p^2)\right],$$

$${\mathbf{\Delta }}_{hHA}(p^2)=-{\left[{\hat{\mathbf{\Gamma }}}_{hHA}(p^2)\right]}^{-1}=-{\left[i\left(p^2\text{1}-\text{M}(p^2)\right)\right]}^{-1}.$$

$${\Delta }_{ij}(p^2)=\frac{{\hat{\Gamma }}_{ij}{\hat{\Gamma }}_{kk}-{\hat{\Gamma }}_{jk}{\hat{\Gamma }}_{ki}}{{\hat{\Gamma }}_{ii}{\hat{\Gamma }}_{jj}{\hat{\Gamma }}_{kk}+2{\hat{\Gamma }}_{ij}{\hat{\Gamma }}_{jk}{\hat{\Gamma }}_{ki}-{\hat{\Gamma }}_{ii}{\hat{\Gamma }}_{jk}^2-{\hat{\Gamma }}_{jj}{\hat{\Gamma }}_{ki}^2-{\hat{\Gamma }}_{kk}{\hat{\Gamma }}_{ij}^2},$$

$${\Delta }_{ii}(p^2)$$

$$=\frac{i}{p^2-m_i^2+{\hat{\Sigma }}_{ii}^{\mathrm{eff}}(p^2)},$$

$${\hat{\Sigma }}_{ii}^{\mathrm{eff}}(p^2)$$

$$p^2-m_i^2+{\hat{\Sigma }}_{ii}^{\mathrm{eff}}(p^2)=0,$$

$${\hat{Z}}_i^a$$

$$\hat{Z}_{ij}^{a}=\frac{\Delta_{ij}(p^{2})}{\Delta_{ii}(p^{2})}\bigg{\rvert}_{p^{2}=\mathcal{M}^{2}_{a}}$$

$${\hat{\text{Z}}}_{aj}$$

$${\hat{\Gamma }}_{h_a}$$

$${\Delta }_{ij}(p^2)$$

$${\Delta }_a^{\text{BW}}(p^2):=\frac{i}{p^2-{\mathcal{M}}_a^2}=\frac{i}{p^2-M_{h_a}^2+i{M}_{h_a}{\Gamma }_{h_a}}.$$

$$\mathcal{A}$$

$${\mathcal{A}}_{h_a}\equiv {\hat{\Gamma }}_{h_a}^X{\Delta }_a^{\text{BW}}(p^2){\hat{\Gamma }}_{h_a}^Y=\sum _{i,j=h,H,A}{\hat{\Gamma }}_i^X{\hat{\text{Z}}}_{ai}{\Delta }_a^{\text{BW}}(p^2){\hat{\text{Z}}}_{aj}{\hat{\Gamma }}_j^Y,$$

$$b\bar{b}\to$$

$$gg\to$$

$$\displaystyle|\mathcal{A}|^{2}_{\textrm{coh}}=\bigg{\lvert}\sum_{a=1}^{3}\mathcal{A}_{h_{a}}\bigg{\rvert}^{2}\,,~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}|\mathcal{A}|^{2}_{\textrm{incoh}}=\sum_{a=1}^{3}\bigg{\lvert}\mathcal{A}_{h_{a}}\bigg{\rvert}^{2}\,.$$

$$|\mathcal{A}{|}_{\text{int}}^2=|\mathcal{A}{|}_{\text{coh}}^2-|\mathcal{A}{|}_{\text{incoh}}^2=\sum _{a<b}2\text{Re}\left[{\mathcal{A}}_{h_a}{\mathcal{A}}_{h_b}^{\ast }\right].$$

$${\eta }^P:=\frac{{\sigma }_{\mathrm{int}}^P}{{\sigma }_{\mathrm{incoh}}^P}.$$

$${\sigma }^P$$

$$w_{ab}^{P,a}=\frac{{\sigma }_{h_a}^P}{{\sigma }_{h_a}^P+{\sigma }_{h_b}^P}\equiv 1-w_{ab}^{P,b},$$

$${\eta }_a^P$$

Full text

1. 1 Introduction

2. 2 The MSSM with complex parameters at tree level

3. 3 Higher-order $\mathcal{CP}$-violating Higgs mixing

4. 3.1 Propagator corrections and the $\hat{\textbf{Z}}$-matrix

5. 3.2 Resonant mixing-enhancement of cross sections

6. 4 Interference contributions to the processes $\left\{b\bar{b},gg\right\}\rightarrow h_{1},h_{2},h_{3}\rightarrow\tau^{+}\tau^{-}$

7. 5 Impact of interference and mixing on exclusion bounds

8. 5.1 Coherent and incoherent cross sections

9. 5.2 Exclusion bounds with mixing and interference

10. 6 Conclusions

1 Introduction

Despite a rich program of searches for additional Higgs bosons at the LHC, no new scalars beyond the state at about $125\,\textrm{GeV}$ have been found so far. Additional Higgs bosons may be hiding from the searches at the LHC where possible reasons can be reduced couplings, large masses or cancellations such as destructive interference effects.

Extended Higgs sectors are predicted by various models beyond the Standard Model (SM) such as the Two-Higgs-Doublet Model (2HDM), the Minimal Supersymmetric Standard Model (MSSM) and singlet extensions such as the Next-to-Minimal Supersymmetric Standard Model (NMSSM). In this paper we focus on interpretations within the MSSM. The most obvious possibility in this model is to identify the discovered state with a mass of $M_{h}^{\rm{exp}}=125.09\pm 0.24\,\textrm{GeV}$ [1] with the lightest MSSM Higgs boson, which implies the existence of two heavier neutral Higgs boson in the spectrum. The possibility that the discovered state could be the next-to-lightest Higgs boson [2, 3] is highly constrained within the MSSM, but not ruled out (see Ref. [4] for a recent update).

Although the properties of the discovered scalar are so far compatible with the ones predicted for the SM Higgs boson within the present experimental uncertainties, significant deviations from the SM are possible in individual Higgs couplings, cross sections and branching ratios. Concerning the $\mathcal{CP}$ properties of the discovered state, a pure $\mathcal{CP}$-odd nature could be ruled out [5, 6], while only very weak bounds exist so far on an admixture of $\mathcal{CP}$-even and $\mathcal{CP}$-odd components. If a non-zero $\mathcal{CP}$-odd admixture of the state at 125 GeV could be experimentally established, this would be a direct manifestation of the presence of $\mathcal{CP}$ violation in the Higgs sector. In general $\mathcal{CP}$ violation gives rise to a mixing between all neutral Higgs bosons in the Higgs sector. However, in the decoupling region of the MSSM the lightest neutral Higgs boson is SM-like and almost purely $\mathcal{CP}$-even. On the other hand, the two heavy neutral Higgs bosons in this case can very significantly differ from $\mathcal{CP}$ eigenstates and have a large mixing with each other. The search for additional heavy Higgs bosons should therefore take into account the possibility that those heavy Higgs bosons are not necessarily $\mathcal{CP}$ eigenstates.

LHC searches for neutral MSSM Higgs bosons, following on from earlier searches at LEP [7] and the Tevatron [8, 9, 10, 11], are carried out in the inclusive gluon fusion production channel, $gg\rightarrow h_{a}~{}(a=1,2,3)$, and in the production process in association with a pair of bottom quarks. For higher-order calculations of the cross sections of these two production processes within the MSSM see e.g. Refs. [12, 13, 14, 15, 16, 17, 18] (assuming $\mathcal{CP}$ conservation in the Higgs sector) and [19, 20, 21, 22] (including $\mathcal{CP}$ violation). The latter process, which dominates in the MSSM for high $\tan\beta$ values owing to an enhanced bottom Yukawa coupling, is usually written as $b\bar{b}\rightarrow h_{a}$ (according to the five-flavour scheme where the bottom quark is regarded as a parton in the proton) or $gg\rightarrow b\bar{b}h_{a}$ (according to the four-flavour scheme without a bottom parton density distribution in the proton). Searches for neutral MSSM Higgs bosons decaying to down-type fermions have been carried out at the LHC in the $\tau^{+}\tau^{-}$ [23, 24, 25, 26, 27], in the $\mu^{+}\mu^{-}$ [28, 29, 30, 31] and in the $b\bar{b}$ [32, 33] decay channels, where the $\tau^{+}\tau^{-}$ channel provides by far the highest sensitivity. The decay modes of the heavy MSSM Higgs bosons to down-type fermions are enhanced at large $\tan\beta$, whereas the branching ratios of the heavy Higgs bosons into vector bosons vanish in the decoupling limit [34].

The results of the searches are on the one hand reported as nearly model-independent limits on the product of the on-shell production cross section (separately for the inclusive production in gluon fusion and the production in association with bottom quarks) and the considered branching ratio of a single scalar resonance, assuming a narrow width. On the other hand, the search results are also interpreted in model-specific contexts using appropriate benchmark scenarios, for instance the $M_{h}^{\rm{mod+}}$ scenario of the MSSM [35].

In the MSSM the signal is potentially comprised of contributions of all three neutral Higgs bosons. The parameter region in which the masses of the two heavier neutral Higgs bosons of the MSSM are significantly heavier than the mass of the lightest neutral Higgs, where the latter needs to be close to 125 GeV for a phenomenologically viable scenario, corresponds to the decoupling region of the MSSM. The two heavier Higgs states are nearly mass-degenerate over this whole region, with mass splittings that are often below the experimental resolution. In the non-decoupling region even the masses of all three neutral Higgs bosons can be close to each other. In all model-specific interpretations at the LHC it has been assumed so far that the signal contributions from different Higgs bosons can be added incoherently, i.e. no interference contributions have been taken into account. Under the assumption that $\mathcal{CP}$ is conserved in the Higgs sector, which is realised in the MSSM with real parameters, the mass eigenstates are states of definite $\mathcal{CP}$, comprising the light and heavy $\mathcal{CP}$-even states $h$ and $H$, as well as the $\mathcal{CP}$-odd state $A$. In this case interference effects between the heavy Higgs states $H$ and $A$ are absent, while the $\mathcal{CP}$-even states $h$ and $H$ can interfere with each other. If the assumption of $\mathcal{CP}$ conservation in the Higgs sector is dropped, interference effects occur also between the two heavy Higgs states, which can be enhanced by a small mass difference between the two states. In the MSSM $\mathcal{CP}$-violation in the Higgs sector is induced by potentially large loop corrections involving complex parameters. As a consequence, in the general case of complex MSSM parameters the $\mathcal{CP}$-even states $h$ and $H$ mix with the $\mathcal{CP}$-odd state $A$ into the mass eigenstates $h_{1},h_{2},h_{3}$.

The phenomenological consequences of $\mathcal{CP}$-violating effects in the Higgs sector of the MSSM have been investigated for the Higgs searches at LEP [36, 7]. In the analyses using the CPX benchmark scenario [37] the presence of a non-zero phase of the trilinear couplings $A_{t,b}$ and the mixing between the $\mathcal{CP}$-even and $\mathcal{CP}$-odd states were found to have a significant impact on the limits, giving rise to unexcluded parameters regions at much smaller values of the lightest Higgs mass than for the $\mathcal{CP}$-conserving case [36, 7, 38, 39].

In the context of Higgs searches at the LHC, interference effects between the light and the heavy neutral Higgs boson for the case of $\mathcal{CP}$-conservation in an extended Higgs sector have been analysed in Ref. [40] for Higgs production in the MSSM via the decay of a heavy neutralino, in Ref. [41] for $gg\rightarrow h/H\rightarrow VV,\,V=W,Z$ within the 2HDM including background contributions, and in Ref. [42] for the singlet extension of the SM. The interference of two light NMSSM Higgs bosons decaying to two photons has been investigated in Ref. [43] for the $\mathcal{CP}$-violating case. For discussions of interference effects of heavy neutral Higgs bosons with each other and with the background at low $\tan\beta$ in the $t\bar{t}$ final state see Refs. [44, 45, 46] and for interference with the background without $H-A$ interference see Refs. [47, 48]

In the present paper we investigate the impact of mixing and interference effects on the search for heavy Higgs bosons at the LHC in the channels $\left\{b\bar{b},gg\right\}\rightarrow h_{1},h_{2},h_{3}\rightarrow\tau^{+}\tau^{-}$, where the production modes denote the inclusive production in gluon fusion and the production in association with bottom quarks as discussed above. In order to incorporate the interference effects we make use of a generalised narrow-width approximation (NWA) [40, 49], where in contrast to the standard NWA interference contributions are taken into account. Our analysis is carried out in a benchmark scenario that is extended such that it contains a non-zero phase of the trilinear coupling $A_{t}$. The exclusion bounds are evaluated with the help of the program HiggsBounds [50, 51, 52, 53]. We find that the mixing effects between the two heavy Higgs bosons are resonantly enhanced in the parameter regions where the two Higgs bosons are nearly mass-degenerate. The corresponding contributions of the intermediate Higgs states would yield a significant increase of the cross section for the case of complex parameters as compared to the case of real parameters if the interference contributions were neglected. However, the interference effect between the two heavy Higgs bosons gives rise to a very large destructive interference contribution, so that the net effect turns out to be a large suppression of the total cross section times branching ratio in the resonance region (see also Refs. [54, 55, 56]). As a consequence, a parameter region remains unexcluded by LHC searches from Run 1 that would appear to be excluded if interference effects were neglected.

The paper is structured as follows. After fixing the notation for the MSSM with complex parameters at tree level in Sect. 2, we summarise higher-order propagator-mixing in the Higgs sector and discuss the mixing-enhancement of Higgs production cross sections in Sect. 3. In Sect. 4, we quantify and investigate the interference contributions in the processes $\left\{b\bar{b},gg\right\}\rightarrow h_{1},h_{2},h_{3}\rightarrow\tau^{+}\tau^{-}$. Subsequently in Sect. 5, we compare the predicted cross sections with mixing and interference contributions to experimental limits and evaluate their impact on the exclusion bounds in a benchmark scenario with a non-zero phase of the trilinear coupling $A_{t}$. Our conclusions are given in Sect. 6.

2 The MSSM with complex parameters at tree level

Before discussing higher-order mixing of Higgs bosons in Sect. 3, in this section we will specify the notation for the MSSM with complex parameters at tree level, following Ref. [57].

Sfermion sector

Sfermions $\tilde{f}_{L},\tilde{f}_{R}$ mix into the mass eigenstates $\tilde{f}_{1},\tilde{f}_{2}$ within one generation according to the mass matrix

[TABLE]

where $X_{f}:=A_{f}-\mu^{*}\cdot\left\{\cot\beta,\tan\beta\right\}$ for $f$ being an up- or down-type quark, respectively. Besides the trilinear couplings $A_{f}=|A_{f}|e^{i\phi_{A_{f}}}$ also the higgsino mass parameter $\mu=|\mu|e^{i\phi_{\mu}}$ can be complex. Starting from one-loop order, these phases may influence the Higgs sector via sfermion loops.

Gluino sector

The mass of a gluino $\tilde{g}^{a},~{}a=1,2,3,$ is given by $m_{\tilde{g}}=|M_{3}|$, where $M_{3}=|M_{3}|\,e^{i\phi_{M_{3}}}$ is the possibly complex gluino mass parameter. The gluino does not directly couple to Higgs bosons. Hence, the phase $\phi_{M_{3}}$ enters the predictions for the Higgs-boson masses and the wave function normalisation factors for external Higgs bosons only from the two-loop level onwards whereas it has an impact for example on the bottom Yukawa coupling already at one-loop order.

Neutralino and chargino sector

The charginos $\widetilde{\chi}_{i}^{\pm},~{}i=1,2$, are admixtures of the charged winos $\widetilde{W}^{\pm}$ and higgsinos $\widetilde{H}^{\pm}$ via the mass matrix

[TABLE]

Likewise the neutralinos $\widetilde{\chi}_{i}^{0},\,i=1,...,4$, are composed of the neutral electroweak gauginos $\widetilde{B},\,\widetilde{W}^{3}$ and the neutral Higgsinos $\tilde{h}_{d}^{0},\,\tilde{h}_{u}^{0}$:

[TABLE]

Thus, at tree-level, mixing in the chargino sector is governed by the higgsino and wino mass parameters $\mu$ and $M_{2}$, respectively, and in the neutralino sector in addition by the bino mass parameter $M_{1}$. Although all of these three parameters can be complex in principle, only two of the phases are independent, and the choice $\phi_{M_{2}}=0$ is a common convention.

Higgs sector

The two complex scalar Higgs doublets of the MSSM are denoted as

[TABLE]

As a possible relative phase between both doublets vanishes at the minimum of the Higgs potential and the phase of the coefficient of the bilinear term in the Higgs potential can be rotated away, the Higgs sector is $\mathcal{CP}$ conserving at lowest order. The tree level mass eigenstates as eigenstates of $\mathcal{CP}$ result from the diagonalisation of the mass matrices of the neutral and the charged components, respectively,

[TABLE]

with the short-hand notation $s_{x}\equiv\sin x,~{}c_{x}\equiv\cos x$. The mixing angle $\alpha$ is associated with the $\mathcal{CP}$-even Higgs bosons $h,H$; $\beta_{n}$ with the neutral $\mathcal{CP}$-odd Higgs $A$ and Goldstone boson $G$, and $\beta_{c}$ with the charged Higgs $H^{\pm}$ and the charged Goldstone boson $G^{\pm}$. The masses of the $\mathcal{CP}$-odd and the charged Higgs bosons are at tree level related by $m_{H^{\pm}}^{2}=m_{A}^{2}+M_{W}^{2}$. At lowest order, the Higgs sector is fully determined by the two SUSY input parameters (in addition to SM masses and gauge couplings) $\tan\beta\equiv v_{u}/v_{d}$ and $m_{H^{\pm}}$ (or, for conserved $\mathcal{CP}$, equivalently $m_{A}$). We use in this paper a lower-case and upper-case notation in order to indicate tree-level and loop-corrected masses, respectively.

In the decoupling limit of $M_{A}\gg M_{Z}$ the MSSM Higgs sector appears SM-like, and the heavy Higgs bosons are difficult to detect in production and decay channels involving gauge bosons. On the other hand, the couplings to fermions can be either suppressed or enhanced, depending on the angles $\alpha$ and $\beta$. In the processes considered in this work, in particular the couplings of the neutral Higgs bosons $i=h,H,A$ to down-type fermions $f_{d}$ such as $\tau$-leptons and $b$-quarks are involved [58],

[TABLE]

As $c_{\alpha}\rightarrow s_{\beta}$ in the decoupling limit, the couplings of $H$ and $A$ to down-type fermions are enhanced by large values of $\tan\beta$.

3 Higher-order $\mathcal{CP}$-violating Higgs mixing

Higher-order corrections have a sizeable impact on the MSSM Higgs sector. In loop diagrams, particles from all other sectors contribute to Higgs observables such that in particular the trilinear couplings $A_{f}$, the stop and sbottom masses, the gluino mass, and in the sub-leading terms the higgsino mass parameter $\mu$, play – besides $M_{H^{\pm}}$ (or $M_{A}$) and $\tan\beta$ as at lowest order – an important role in the Higgs phenomenology. While $\mathcal{CP}$ is conserved in the Higgs sector at lowest order, those parameters from other sectors that are a priori complex can introduce $\mathcal{CP}$ violation in the Higgs sector at higher orders. In the $\mathcal{CP}$-conserving case, there is only a $2\times 2$ mixing among the two $\mathcal{CP}$-even neutral Higgs bosons, $h$ and $H$. On the contrary, the non-zero phases of complex parameters cause a $\mathcal{CP}$-violating mixing of the scalars $h,H$ and the pseudoscalar $A$ into the mass eigenstates $h_{1},h_{2},h_{3}$ (in addition, there is also mixing with the neutral Goldstone and vector bosons that we will neglect in this work, see Ref. [57] for a discussion). This $3\times 3$ mixing structure is reflected in the mass matrix of the $h,H,A$ system and in the on-shell wave-function normalisation factors, $\hat{\textbf{Z}}$, as we will outline in Sect. 3.1 following Refs. [57, 39, 59]. In Sect. 3.2 we will discuss the consequences of highly admixed mass eigenstates on the production cross sections of each resonance.

For the calculation of radiative corrections, we adopt the hybrid on-shell and $\overline{\rm{DR}}$-renormalisation scheme defined in Ref. [57] such that the masses are renormalised on-shell whereas $\overline{\rm{DR}}$ conditions are employed for the fields and $\tan\beta$. We evaluate the Higgs masses, widths, branching ratios and self-energies with FeynHiggs-2.10.2 [60, 61, 62, 63, 64]111The additional corrections included in more recent versions do not qualitatively alter the effects discussed in this paper. The latest version will be used in an update for the Run 2 results of the LHC searches [65]., including the full momentum-dependent one-loop corrections and leading two-loop contributions.

3.1 Propagator corrections and the $\hat{\textbf{Z}}$-matrix

Propagator matrix

In the general case, all renormalised self-energies $\hat{\Sigma}_{ij}(p^{2})$ of the Higgs bosons $i,j=h,H,A$ are non-zero. Hence, the mass-square matrix M of the neutral Higgs bosons does not only consist of the tree-level masses-squares $m_{i}^{2}$ in the diagonal entries, but also of the momentum-dependent renormalised self-energies $\hat{\Sigma}_{ij}(p^{2})$ on the diagonal and off-diagonal positions. Therefore the matrix of mass-squares becomes

[TABLE]

The self-energies $\hat{\Sigma}_{ij}p^{2}$ also contribute to the renormalised irreducible two-point vertex functions,

[TABLE]

whose elements form the $3\times 3$ matrix $\hat{\boldsymbol{\Gamma}}_{hHA}(p^{2})$ that is related to the mass-square matrix $\textbf{M}(p^{2})$ and to the propagator matrix $\boldsymbol{\Delta}(p^{2})$ via

[TABLE]

The elements of $\boldsymbol{\Delta}(p^{2})$ are the propagators $\Delta_{ij}(p^{2})$ starting as Higgs boson $i$ and ending on $j$ with all possible mixings in between. The off-diagonal propagators ($i\neq j$) have the form

[TABLE]

where the explicit $p^{2}$-dependence of all $\hat{\Gamma}(p^{2})$ has been suppressed in order to simplify the notation. The diagonal propagator can be expressed as

[TABLE]

where in the last step the effective self-energy $\hat{\Sigma}^{\rm{eff}}_{ii}$ is defined as

[TABLE]

starting with the diagonal self-energy $\hat{\Sigma}_{ii}$ at one-loop order and containing the mixing 2-point functions whose products in $\hat{\Sigma}^{\rm{eff}}_{ii}$ arise from two-loop order onwards. The complex poles $\mathcal{M}^{2}_{a}$, where $a=1,2,3$ label the mass eigenstates, are determined as solutions of the equations

[TABLE]

with $i=h,H,A$.

Z-factors

For the correct normalisation of the $S$-matrix elements appropriate on-shell properties of all external particles are required. Pure on-shell renormalisation schemes guarantee these properties, whereas the $\overline{\rm{DR}}$ field renormalisation conditions make the finite on-shell wave-function normalisation factors necessary which ensure that external particles have unit residue and that all mixing contributions vanish on-shell. According to the LSZ formalism [66], they are obtained from the residue of the propagator of the external field, here a neutral Higgs boson $i=h,H,A$, evaluated at one of the complex poles $\mathcal{M}^{2}_{a},~{}a=1,2,3$ [67, 68, 38, 39, 59],

[TABLE]

The normalisation of each external particle $i$ at the mass shell $p^{2}=\mathcal{M}^{2}_{a}$ is accounted for by a factor of $\sqrt{\hat{Z}_{i}^{a}}$. Furthermore, if mixing between $i$ and $j$ occurs on an external line at $\mathcal{M}^{2}_{a}$, the transition factor

[TABLE]

restores the correct normalisation. The products of the overall normalisation factors $\sqrt{\hat{Z}_{a}}$ and transition ratios $\hat{Z}_{aj}$ (note the difference between $\hat{Z}_{aj}$ and $\hat{\textbf{Z}}_{aj}$),

[TABLE]

form the non-unitary matrix $\hat{\textbf{Z}}$. The imaginary parts of the self-energies of unstable particles, evaluated at non-vanishing incoming momentum, cause $\hat{\textbf{Z}}$ to be non-unitary. Accordingly, it does not provide a unitary transformation between the lowest order states $h,H,A$ and the loop-corrected mass eigenstates $h_{1},h_{2},h_{3}$.

The $\hat{\textbf{Z}}$-matrix allows one to express the one-particle irreducible (1PI) vertex functions $\hat{\Gamma}_{h_{a}}$ involving an external loop-corrected mass eigenstate $h_{1},h_{2},h_{3}$ as a linear combination of the 1PI vertex functions of the lowest-order states, $\hat{\Gamma}_{i}$:

[TABLE]

where the ellipsis denotes additional contributions from the mixing of Higgs bosons with Goldstone and vector bosons, which are not contained the $\hat{\textbf{Z}}$-matrix and need to be calculated explicitly if needed. As in the full propagator matrix $\boldsymbol{\Delta}_{hHA}(p^{2})$, if $\mathcal{CP}$ is conserved in the Higgs sector, the $\hat{\textbf{Z}}$-matrix is reduced to the $2\times 2$ mixing among the $\mathcal{CP}$ even states $h$ and $H$, with $\hat{\textbf{Z}}_{3A}=1$, $\hat{\textbf{Z}}_{a3}=\hat{\textbf{Z}}_{3i}=0$ for $a=1,2$ and $i=h,H$. If, on the contrary, complex parameters are present, all elements of $\hat{\textbf{Z}}$ can have non-vanishing values and represent the $3\times 3$ mixing.

Mixing propagators in terms of the mass eigenstates

While the propagators $\Delta_{ij}(p^{2})$ contain the full momentum dependence of the lowest-order propagators and the self-energies, the mixing properties and the leading momentum dependence around the complex pole is adequately approximated by the sum of Breit-Wigner propagators of each mass eigenstate in combination with the on-shell wave-function normalisation factors, $\hat{\textbf{Z}}$. The expansion of the full propagator $\Delta_{ij}(p^{2})$ around all of its complex poles $\mathcal{M}^{2}_{a}=M_{h_{a}}^{2}-iM_{h_{a}}\Gamma_{h_{a}}$ results in [59]

[TABLE]

where the resonant contribution of a mass eigenstate $h_{a},~{}a=1,2,3$, with the loop-corrected mass $M_{h_{a}}$ and the total width $\Gamma_{h_{a}}$ is described by a Breit–Wigner propagator with a constant width,

[TABLE]

Hence, although the $\hat{\textbf{Z}}$ matrix has been introduced for the correct normalisation of external Higgs bosons, it is also applicable to internal Higgs bosons to account for the higher-order mixing properties to a high accuracy. Eq. (20) implies that the amplitude $\mathcal{A}$ of a process with Higgs exchange between the initial $X$ and final state $Y$ can be written as

[TABLE]

with $\hat{\Gamma}_{h_{a}}^{X}=\sum_{i=h,H,A}\hat{\textbf{Z}}_{ai}\hat{\Gamma}_{i}^{X}$ as in Eq. (19). The formulation on the RHS of Eq. (22) allows one to divide the full amplitude into the separate contributions of each resonance $h_{a}$:

[TABLE]

which will be useful in Sect. 4 to compare the coherent and incoherent sums of amplitudes of several mass eigenstates.

3.2 Resonant mixing-enhancement of cross sections

Since the $\hat{\textbf{Z}}$-matrix is non-unitary, cross sections including mixing contributions (with $\hat{\textbf{Z}}$ from Eqs. (16-18)) can be enhanced compared to the case without mixing, which would correspond to a diagonal $\hat{\textbf{Z}}$-matrix (with $\hat{\textbf{Z}}=\mathbb{1}$ for the case where the external states are properly normalised). In particular, in the case of quasi-degenerate states a resonance enhancement is possible as a consequence of significant off-diagonal and imaginary parts of the self-energies and the $\hat{\textbf{Z}}$-factors.

The $3\times 3$ mixing occurring in the general case of complex parameters can change the phenomenology very significantly as compared to the $\mathcal{CP}$-conserving case of real parameters. This is in particular due to the $\mathcal{CP}$-violating mixing of $H,A$ into $h_{2},h_{3}$ since the heavy Higgs bosons $h_{2},h_{3}$ are almost mass-degenerate in the wide range of parameter space characterised by the decoupling limit. This degeneracy and the properties of $h_{2}$ and $h_{3}$ lead to a significant $H,A$-mixing and an enhancement of the squared amplitudes.

For a large mass difference between the two heavy states $H$ and $A$ and the light state $h$, the mixing of $h$ with $H$ and $A$ can be neglected in good approximation. In this limit, the squared amplitudes of the production processes of $h_{2}$ and $h_{3}$, $\left|\mathcal{P}_{h_{a}}\right|^{2}$ for $a=2,3$, scale with $|\hat{\textbf{Z}}_{aH}|^{2}\,|\hat{\Gamma}_{H}|^{2}+|\hat{\textbf{Z}}_{aA}|^{2}\,|\hat{\Gamma}_{A}|^{2}$, which can become significantly larger than $|\hat{\Gamma}_{i}|^{2}$, $i=H,A$. In the case of $|\hat{\Gamma}_{H}|\simeq|\hat{\Gamma}_{A}|$, this scaling factor further simplifies to $|\hat{\textbf{Z}}_{aH}|^{2}+|\hat{\textbf{Z}}_{aA}|^{2}$, which can exceed 1. Hence, the cross sections of $\sigma(pp\rightarrow h_{2})$ and $\sigma(pp\rightarrow h_{3})$ receive an enhancement from the mixing contributions. We are particularly interested here in the production modes in association with a pair of bottom quarks ($b\bar{b}$) and via gluon fusion ($gg$).

As an example case featuring such an enhancement effect, we illustrate the $\mathcal{CP}$-violating mixing in a modified version of the $M_{h}^{\rm{mod+}}$ scenario, originally defined in Ref. [69], where we introduce as in Refs. [55, 59, 56] the complex phase $\phi_{A_{t}}=\frac{\pi}{4}$, set $\phi_{A_{b}}=\phi_{A_{\tau}}=\phi_{A_{t}}$, and increase the value of the higgsino mass parameter to $\mu=1000\,\textrm{GeV}$ (as also proposed in Ref. [69]), which amplifies the $\mathcal{CP}$-violating terms. The input parameters being varied are $M_{H^{\pm}}$ and $\tan\beta$. We refer to this scenario as $\mathbb{C}M_{h}^{\rm{mod+}}$.

Fig. 1 shows the single-particle production cross sections $\sigma(P\rightarrow h_{a})$ for $h_{2}$ and $h_{3}$ in association with a pair of bottom quarks and via gluon fusion. The cross sections for the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario with $\phi_{A_{t}}=\pi/4$, denoted by $3\times 3$, have been normalised in Fig. 1 to the ones for the corresponding $\mathcal{CP}$-conserving scenario with $\phi_{A_{t}}=0$, denoted by $2\times 2$. The predictions for the cross sections have been obtained with FeynHiggs. The displayed results for the ratios $R_{\sigma}(P,h_{a}):=\sigma(P\rightarrow h_{a})_{3\times 3}/\sigma(P\rightarrow h_{a})_{2\times 2}$ show enhancements by factors of more than 5 in certain areas of parameter space where $h_{2}$ and $h_{3}$ are highly admixed states, especially around $M_{H^{\pm}}\simeq 480\,\textrm{GeV}$ and $\tan\beta\simeq 29$. Accordingly, the enhancement can indeed be related to a resonance-type mixing between the $\mathcal{CP}$ eigenstates $H$ and $A$. The observed patterns are similar for $h_{2}$ (Figs. 1, 1) and $h_{3}$ (Figs. 1, 1) and, to a lesser extent, also for production via $b\bar{b}$ (Figs. 1, 1) and $gg$ (Figs. 1, 1). The observed similarities are related to the involved pattern of $\hat{\textbf{Z}}$-factors, as discussed above. The differences between the production via $b\bar{b}$ and $gg$ arise from the structure of the amplitudes. While the production in association with a pair of bottom quarks is dominated by the respective Higgs coupling to bottom quarks (where the phase $\phi_{A_{t}}$ enters via the correction to the relation between the bottom-quark mass and the bottom Yukawa coupling, $\Delta_{b}$, see e.g. Ref. [22]), the loop-induced production via gluon fusion at leading order comprises contributions from the bottom- and top-quark loops as well as from their scalar superpartners (the latter induce an explicit dependence on $\phi_{A_{t}}$).222In the $\mathcal{CP}$-conserving case $A$-boson production in gluon fusion does not receive a squark-loop contribution at leading-order. Those contributions from up-type (s)fermions to gluon fusion have a sizeable impact at lower values of $\tan\beta$. The edge along $\tan\beta\simeq 30$ for $M_{H^{\pm}}\gtrsim 600\,\textrm{GeV}$ in the plots for $gg$-induced production is caused by the swap of the composition of $h_{2}$ and $h_{3}$ being mostly $H$- or $A$-like, respectively.

While so far we have restricted our discussion to the individual cross sections of $h_{2}$ and $h_{3}$, in the next section we will analyse the $\mathcal{CP}$-violating interference contributions to the full process involving production and decay.

4 Interference contributions to the processes

$\left\{b\bar{b},gg\right\}\rightarrow h_{1},h_{2},h_{3}\rightarrow\tau^{+}\tau^{-}$

For the investigation of a full process of Higgs production and decay where the Higgs bosons appear as internal propagators, we calculate the interference of amplitudes with $s$-channel $h_{1},h_{2},h_{3}$ exchange. We consider the production processes via $b\bar{b}$ and $gg$ and the subsequent decay into a pair of $\tau$-leptons, as illustrated in Fig. 2,

[TABLE]

In order to be able to make use of precise predictions for the separate production [22] and decay [61] processes we focus here on the relative interference contribution. For the phenomenological prediction of the complete processes including interference, higher-order corrections and quark- and gluon-luminosities at the LHC, we will combine those relative interference contributions with the computations of hadronic cross sections and branching ratios from the existing literature including available higher-order corrections. Regarding the determination of the relative interference contributions, we restrict ourselves to the computation of the relevant $2\to 2$ partonic processes at leading order (LO). This means that process (24), which involves $b\bar{b}h_{a}$-associated production, is treated at tree-level, while the leading contributions to gluon fusion entering process (25) consist of one-loop diagrams with quarks $q$ and squarks $\tilde{q}$, in particular those of the third generation. Besides, we take propagator-type corrections into account by using Higgs masses, total widths and $\hat{\textbf{Z}}$-factors from FeynHiggs including the full one-loop and dominant two-loop corrections. The procedure to evaluate the relative interference contributions on the basis of leading-order diagrams is motivated by the fact that the vertex corrections to the production and decay factorise (apart from non-factorisable initial-to-final state radiation). It furthermore avoids double-counting of higher-order contributions that are incorporated in the separate production and decay processes. We calculate the LO contributions to the amplitude and the cross section of the full process of production and decay using FeynArts [70, 71, 72, 73, 74] with a model file containing $\hat{\textbf{Z}}$-factors for the Higgs vertices (in this way the Higgs propagators are expressed in terms of $\hat{\textbf{Z}}$-factors and lowest-order Breit–Wigner propagators as given in (20)), FormCalc [75, 76, 77, 78, 79] and, as mentioned above, FeynHiggs for quantities from the Higgs sector.

Quantifying the interference

In order to determine the interference term in each of the two processes (24, 25), we distinguish between the coherent sum of the $2\rightarrow 2$ amplitudes with $h_{1},h_{2},h_{3}$-exchange including the interference and their incoherent sum without the interference,

[TABLE]

The interference term is then obtained from the difference,

[TABLE]

Accordingly, the coherent cross section refers to the cross section based on the coherent sum of amplitudes, $\sigma_{\textrm{coh}}\equiv\sigma\left(|\mathcal{A}|^{2}_{\textrm{coh}}\right)$, and likewise for the incoherent cross section $\sigma_{\textrm{incoh}}$ (omitting the interference term) and the interference part of the cross section $\sigma_{\textrm{int}}$. Furthermore, we define the relative interference term for each production mode $P=b\bar{b},gg$ as

[TABLE]

Beyond this overall interference term, it is also useful to differentiate between the three “squared” terms of each Higgs boson $h_{a}$ in $\sigma_{h_{a}}$ and the three interference terms $\sigma_{\rm{int}_{ab}}$ involving two Higgs bosons $h_{a},h_{b}$ out of the three, respectively, again separately for both production modes $P$,

[TABLE]

In order to express the cross section into separate Higgs contributions (as in the incoherent case), we introduce relative interference factors $\eta_{a}^{P}$ that rescale the individual Higgs rates:

[TABLE]

where we define $\eta_{a}^{P}$ by weighting $\sigma_{\rm{int}_{ab}}^{P}$ according to the $\sigma_{h_{a}},~{}\sigma_{h_{b}}$ contributions with the weight (for $a,b,c=1,2,3$)

[TABLE]

resulting in

[TABLE]

This definition of $\eta_{a}^{P}$ is stable also in the case that one of the individual contributions $\sigma_{h_{a}}^{P}$ is suppressed.

The relative interference contributions $\eta_{a}^{P}$ can then be applied to obtain a result consisting of the separate predictions for the production and decay processes and the respective interference contributions

[TABLE]

In order to confront the theoretical predictions in the considered scenario with the existing experimental bounds we use the program HiggsBounds-4.2.0 [50, 51, 52, 53]. For the incorporation of the interference contributions it is convenient to rescale the ratio of the production cross sections $P\rightarrow h_{a}$ in the MSSM with respect to the SM, which is supplied as an input to HiggsBounds, in the following way,

[TABLE]

The choice of associating the interference contributions with the production processes in HiggsBounds while leaving the branching ratios $\textrm{BR}(h_{a}\rightarrow\tau^{+}\tau^{-})$ unchanged is of course of purely technical nature. One always needs to ensure that interference contributions are properly taken into account for the full process including production and decay, see Ref. [65].

For the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario the results $\eta^{b\bar{b}}$ and $\eta^{gg}$ as defined in (28), given in $\%$, are shown in Fig. 3 within the ($M_{H^{\pm}}$, $\tan\beta$) parameter planes. The interference effect in both processes exhibits a similar pattern. First of all, it is destructive, i.e. $\eta^{P}<0$, throughout the parameter plane (apart from the small dark red regions). Due to the approximate mass degeneracy of $h_{2}$ and $h_{3}$ and the sizeable $H-A$ mixing, the interference contribution becomes in both processes very large in significant parts of the parameter space. At $M_{H^{\pm}}\simeq 480\,\textrm{GeV}$ and $\tan\beta\simeq 29$, the relative interference contribution reaches a minimum of $\eta^{P}\simeq-97\%$ in both processes so that the cross section is almost completely cancelled by the drastic, negative interference term. This strongest interference region is surrounded by a “valley” of also substantial destructive interference.

Conditions for a sizeable interference term

Having observed the occurrence of a large destructive interference contribution at the example of the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario in Fig. 3, we now investigate in more detail under which conditions such an interference can arise. An important criterion is the overlap of two resonances which is determined by the mass splitting and the total widths. In the decoupling limit, $h_{2}$ and $h_{3}$, which are involved in the relevant interference, are quasi degenerate while $h_{1}$ is much lighter. Because of the significant total widths $\Gamma_{h_{2}},\Gamma_{h_{3}}$, the $h_{2}-h_{3}$ overlap

[TABLE]

is also fulfilled for not exactly degenerate states as displayed in Fig. 4.

In order to understand the location of the strongest interference, we examine the couplings that play a role in the interference term compared to those in the incoherent sum of the $b\bar{b}$-initiated process:

[TABLE]

where $g_{h_{a}f\bar{f}}$ with $a=1,2,3,~{}f=\tau,b$ are the tree-level couplings $g_{if\bar{f}}$ from Eq. (7) for $i=h,H,A$, combined with two-loop $\hat{\textbf{Z}}$-factors from FeynHiggs according to Eq. (19):

[TABLE]

Since the masses $m_{\tau},m_{b}$ and other constants cancel out in Eq. (36), the ratio $c_{23}$ is determined by the on-shell $\hat{\textbf{Z}}$-factors and the angles $\cos\alpha$ and $\sin\beta$. Fig. 4 shows that the behaviour of $c_{23}$ largely determines the pattern of the relative interference contribution observed in Fig. 3, whereas effective couplings based on a $p^{2}=0$ approximation of the self-energies or the pure tree-level couplings would yield a completely different pattern, see Ref. [55]. The interference contribution on the squared matrix element level is proportional to [54]

[TABLE]

In the decoupling region of $M_{H}^{\pm}\gg M_{Z}$, the heavy Higgs bosons $h_{2}$ and $h_{3}$ have very similar masses $M_{h_{2}}\simeq M_{h_{3}}$ and widths $\Gamma_{h_{2}}\simeq\Gamma_{h_{3}}$ so that the product $\Delta^{\textrm{BW}}_{2}(s)\Delta^{\textrm{BW*}}_{3}(s)\simeq|\Delta^{\textrm{BW}}_{2}|^{2}$ becomes approximately real. In this limit, the relations $\hat{\textbf{Z}}_{2H}\simeq\hat{\textbf{Z}}_{3A}$, $\hat{\textbf{Z}}_{2A}\simeq-\hat{\textbf{Z}}_{3H}$ and $\cos\alpha\simeq\sin\beta$ simplify Eq. (38) to:

[TABLE]

As expected, the interference term of $h_{2}$ and $h_{3}$ vanishes for the case of $\mathcal{CP}$ conservation in the Higgs sector, in which case $\hat{\textbf{Z}}_{2A}=\hat{\textbf{Z}}_{3H}=0$ holds. In addition to the $\mathcal{CP}$-violating mixing, a non-zero interference term requires non-vanishing imaginary parts of the $\hat{\textbf{Z}}$-factors, which originate from the imaginary parts of the Higgs self-energies. Consequently, replacing the on-shell $\hat{\textbf{Z}}$-factors by real mixing factors U (see e.g. Refs. [57, 59]) in an effective coupling approach with $p^{2}=0$ renders the interference term zero in the decoupling limit even though the U-matrix may contain equally large diagonal and off-diagonal elements. Even if the conditions of $\mathcal{CP}$-violating mixing and the presence of non-vanishing imaginary parts are fulfilled, there might still be a cancellation between the two terms within the bracket in Eq. (39).

Outside the decoupling limit, with unequal masses, widths and mixing properties of $h_{2}$ and $h_{3}$, the full product of angles from the couplings, $\hat{\textbf{Z}}$-factors and complex Breit-Wigner functions has to be taken into account. However, in the relevant part of the considered parameter plane, the decoupling limit is reached. Given the quasi-degeneracy of $M_{h_{2}}$ and $M_{h_{3}}$ shown in Fig. 4, the structure of the $\hat{\textbf{Z}}$-matrix provides in fact a well-suited indication of the relevance of the interference term. In particular, the square of the bracket and the absolute square of the Breit-Wigner function in combination with the overall minus sign in Eq. (39) explain the observed destructive interference effect.

5 Impact of interference and mixing on exclusion bounds

As two consequences of the $\mathcal{CP}$-violating mixing, in the previous sections we have investigated the enhancement of production cross sections for the mass eigenstates and the reduction of the complete processes of production and decay due to destructive interference. In this section we will compare these theory predictions with experimental limits from searches for heavy additional Higgs bosons in the $\tau^{+}\tau^{-}$ final state at Run 1 of the LHC [26, 80]. In the LHC searches for additional Higgs bosons in supersymmetric models so far it has been assumed that $\mathcal{CP}$ is conserved, so that no mixing occurs between $\mathcal{CP}$-even and $\mathcal{CP}$-odd states. Under this assumption, the contributions from $\mathcal{CP}$-even and $\mathcal{CP}$-odd states can be added as an incoherent sum. The question arises in this context in how far the exclusion bounds that have been obtained under the assumption of $\mathcal{CP}$ conservation will be modified in the general case where the possibility of $\mathcal{CP}$-violating mixing and interference effects is taken into account. In the following we will first as an illustration compare single cross section limits with the theory predictions, before performing a more detailed analysis of the exclusion limits in the $(M_{H^{\pm}},\tan\beta)$ parameter plane of the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario using the tool HiggsBounds.

5.1 Coherent and incoherent cross sections

In Fig. 5 we show the theory predictions (turquoise) for $\sigma(pp\rightarrow P\rightarrow h_{2},h_{3}\rightarrow\tau^{+}\tau^{-})$ depending on the mass $M_{\phi}$ of a neutral scalar resonance $\phi$ for each of the production modes $P=b\bar{b},gg$ in comparison with the respective experimental limits on the production of a single resonance (black). The displayed experimental results represent the expected (dotted) and observed (solid) CMS exclusion bound at the $95\%$ CL from Run 1 of the LHC with $\int\mathcal{L}=24.6\,\rm{fb}^{-1}$ reported in Ref. [81] (which is also used in HiggsBounds).

The theory prediction has been calculated using $M_{H^{\pm}}$ and $\tan\beta$ as input parameters and is shown as a function of $M_{\phi}:=M_{h_{3}}$ ($\simeq M_{h_{2}}$ in the relevant region) for representative values of $\tan\beta$. The differences between the incoherent sum restricted to $\mathcal{CP}$-conserving $2\times 2$ mixing for $\phi_{A_{t}}=0$ (dotted), the incoherent sum with $\mathcal{CP}$-violating $3\times 3$ mixing for $\phi_{A_{t}}=\pi/4$ (dashed) and the coherent sum including $3\times 3$ mixing and the interference term for $\phi_{A_{t}}=\pi/4$ (solid) are visible. Figs. 5 and 5 show the $bb$-initiated process, whereas 5 and 5 show the $gg$-initiated one. The two upper plots correspond to $\tan\beta=29$, where the strongest interference is observed, while the lower two plots represent $\tan\beta=25$ and $\tan\beta=19$, respectively, with still a large interference albeit less drastic effect in the $gg$-process. The comparison of the predicted cross sections times branching ratios with and without the interference term to the experimental bound reveals that the presence of a non-vanishing value of the phase of the parameter $A_{t}$ in this scenario has a large impact on the exclusion bounds. In comparison with the case where $\mathcal{CP}$ conservation is assumed ($\phi_{A_{t}}=0$), the cross section in the case of $3\times 3$ mixing would be significantly enhanced in the resonance region if only the incoherent sum of the contributions were taken into account. In contrast, the incorporation of the interference contribution in fact has the overall effect of a reduction of the cross section in the resonance region.

In the comparison of the theory predictions with a single cross section limit considered here, for a given value of $\tan\beta$ one obtains an exclusion at the $95\%$CL for $M_{\phi}$ values where the prediction of $\sigma\times\textrm{BR}$ is larger than the observed limit. Accordingly, the exclusion based on the mixing-enhanced $3\times 3$ incoherent cross section would be stronger than for the case where $\mathcal{CP}$ conservation is assumed ($2\times 2$ with $\phi_{A_{t}}=0$), while the inclusion of the interference contribution in the full result leads to weaker bounds than for the case where $\mathcal{CP}$ conservation is assumed. Therefore in general the lower bound on $M_{\phi}$ for a given value of $\tan\beta$ that can be set by this comparison is reduced as compared to the case where $\mathcal{CP}$ conservation is assumed.

While this shift of the lower bound on $M_{\phi}$ is only moderate if the intersection between the prediction and the observed limit occurs outside of the resonance region, as it is the case in Fig. 5 and in Fig. 5 for the $gg$ production mode (which is less constraining than the $b\bar{b}$ process in the parameter regions considered here), the situation is different if the large destructive interference contribution in the resonance region reduces the cross section such that it falls below the observed limit. This can be seen in Fig. 5 for the $b\bar{b}$ production mode. In addition to the shift in the lower bound on $M_{\phi}$ from about 695 GeV for the case $\phi_{A_{t}}=0$ to about 675 GeV for the case $\phi_{A_{t}}=\pi/4$, in this case a range of unexcluded values of $M_{\phi}$ opens up between 465 GeV and 485 GeV which would appear to be excluded under the assumption of $\mathcal{CP}$ conservation. A similar effect can be seen in Fig. 5 for the $b\bar{b}$ production mode for the smaller value of $\tan\beta$ which is somewhat off the resonance. As a consequence, the resulting interference effect and also the enhancement of the incoherent cross section are less sharp. In this case the lower bound on $M_{\phi}$ is shifted from about 650 GeV for the case $\phi_{A_{t}}=0$ to about 520 GeV for the case $\phi_{A_{t}}=\pi/4$. If instead the interference contribution were neglected for the $3\times 3$ case, the incoherent cross section would indicate an exclusion on $M_{\phi}$ up to about 690 GeV in this case.

In order to display the behaviour in the resonance region in more detail, Fig. 6 shows the results for the coherent cross sections with $\phi_{A_{t}}=\pi/4$ (including the interference contribution) for different values of $\tan\beta$ across the resonance region. Both for the $b\bar{b}$ and $gg$ production modes, the most pronounced interference effect occurs at about $\tan\beta=29$. The interference effect is similarly large for $\Delta\tan\beta=\pm 1$, while it can be seen to be largely washed out and significantly shifted towards smaller (larger) values of $M_{\phi}$ for $\Delta\tan\beta=+4$ ($\Delta\tan\beta=-4$) in this example. As before, the sensitivity of the $gg$ production mode is much smaller in this parameter region, which gives rise to the fact that the main modification of the cross section caused by the interference contribution occurs in the region much below the excluded cross section. For the $b\bar{b}$ production mode, on the other hand, the interference contribution can reduce the cross section to a level below the exclusion limit within the resonance region.

In should be noted in this context that the comparison between the predicted cross section and the observed exclusion limit is of course affected by theoretical (see the discussion in Ref. [15]) and experimental [81] uncertainties. Since the main purpose of our discussion has been the illustration of the qualitative features of the $3\times 3$ mixing and the interference contribution, we do not address the issue of uncertainties in detail here. Assuming comparable experimental uncertainties in the scenario considered in this work and in the standard $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario analysed in Refs. [81, 25], the impact of the mixing and interference effects extends beyond the experimental uncertainty band. Concerning the theoretical uncertainty of the interference contribution, it should be noted that our approach outlined in Sect. 3 involves a systematic treatment of the finite widths of the resonantly produced Higgs bosons, which is crucial for a proper treatment of the interference effects. For a more detailed discussion of the uncertainties related to the incorporation of the interference contribution we refer to Ref. [65].

5.2 Exclusion bounds with mixing and interference

We now turn to an analysis of the exclusion limits in the $(M_{H^{\pm}},\tan\beta)$ parameter plane of the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario using the program HiggsBounds-4.2.0 [50, 51, 52, 53], in which comprehensive information about limits from a variety of searches at LEP, the Tevatron and Run 1 of the LHC in different channels is implemented. Among the channels that are most important for constraining the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario are the ATLAS and CMS searches for neutral heavy Higgs bosons decaying to $\tau\tau$ [25, 26] and to $ZZ(\rightarrow 4l)$ [82, 83], as well as analyses of the neutral Higgs boson at a mass of $125\,\textrm{GeV}$ [84, 85] and searches for charged Higgs bosons [86, 87]. For each parameter point, HiggsBounds determines the search channel at LEP, Tevatron or the LHC with the highest expected sensitivity for an exclusion. For this channel the observed limit is confronted with the predicted cross section in order to determine whether the considered parameter point is excluded at the 95% CL. For the channels $\sigma_{b\bar{b}}(\phi)\,\textrm{BR}(\phi\rightarrow\tau^{+}\tau^{-})$ and $\sigma_{gg}(\phi)\,\textrm{BR}(\phi\rightarrow\tau^{+}\tau^{-})$ the 2-dimensional likelihood information provided by CMS [80, 25] has been implemented333The recent HiggsBounds-5-$\beta$ update [88] which includes also Higgs search results from Run 2 of the LHC will be used in Ref. [65].. For the incorporation of interference effects we have rescaled the ratio of production cross sections used as input for HiggsBounds as described in (33) and (34).

In order to disentangle the interference effects from other effects due to the modified scenario, the result for the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario with $\phi_{A_{t}}=\pi/4$ and $\mu=1000\,\textrm{GeV}$ is shown in Fig. 7 together with the familiar exclusion region in the well-known $M_{h}^{\rm{mod+}}$ scenario [69] with real parameters (in particular $\phi_{A_{t}}=0$) and the default value of $\mu=200\,\textrm{GeV}$. The latter exclusion region is indicated in green in Fig. 7, where the unexcluded parameter region in this type of plot is often denoted as the “LHC wedge” region. The corresponding exclusion region for $\mu=1000\,\textrm{GeV}$ is shown in orange. The larger value of $\mu$ yields stronger exclusion limits since in this case the decay channel of a heavy Higgs boson into higgsino-like neutralinos and charginos is kinematically closed, resulting in an enhancement of the branching ratio into $\tau^{+}\tau^{-}$. The excluded region that one would obtain for the incoherent sum of cross sections of the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario with $\phi_{A_{t}}=\pi/4$ and $\mu=1000\,\textrm{GeV}$ ($3\times 3$ mixing, no interference contribution taken into account) is displayed in blue in the same plot. We find that the incoherent cross section of the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario with $\phi_{A_{t}}=\pi/4$ would result in a stronger exclusion than the corresponding scenario with real parameters and the same value of $\mu$ because of the mixing-enhanced production cross sections as discussed in Sects. 3.2 and 5.1.

In Fig. 7 the impact of the interference contributions on the excluded regions in the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario is investigated. The blue area (the same as in Fig. 7) corresponds to the incoherent sum of cross sections without the interference contribution. In contrast, taking the interference contribution into account in both the $b\bar{b}$- and the $gg$-initiated processes leads to the weaker exclusion limit given by the red area, as a consequence of the destructive interference effect as discussed in Sect. 5.1. The black line indicates the boundary of the excluded region that one would obtain if the interference contribution were taken into account only for the $b\bar{b}$-initiated process, while the grey line indicates the boundary of the excluded region that one would obtain if the interference contribution were taken into account only for the $gg$-initiated process.

The features visible in Fig. 7 can be related to the results of Sect. 5.1, since each horizontal line of constant $\tan\beta$ in Fig. 7 approximately corresponds to one plot as in Fig. 5 for the same $\tan\beta$, separately for each process. It should be noted in this context, however, that in Fig. 7 also contributions from the lightest neutral Higgs boson of the MSSM are taken into account, and the experimental information implemented in HiggsBounds is employed, consisting in particular of a 2-dimensional likelihood distribution [89] for the production via $b\bar{b}$ and $gg$, instead of a single cross-section limit. Furthermore, the quantity on the x-axis of Fig. 7 is $M_{H^{\pm}}$, while the results in Fig. 5 and Fig. 6 are expressed in terms of the mass of a (single) neutral resonance, $M_{\phi}$. While the interference contributions have little impact on the exclusion region for $\tan\beta$ values below about 15, their effect is clearly visible for $\tan\beta\gtrsim 15$ up to the highest values of $\tan\beta$ shown in the plot. The unexcluded “fjord” between about ($M_{H^{\pm}}=480\,\textrm{GeV},\tan\beta=29$) and ($M_{H^{\pm}}=600\,\textrm{GeV},\tan\beta=20$) is a consequence of the destructive interference contributions in the resonance region discussed in Sect. 5.1. It is interesting to note that the full “fjord” only occurs as a consequence of incorporating the interference contributions in both the $b\bar{b}$- and $gg$-initiated processes. If the interference contribution were incorporated only for the $b\bar{b}$-initiated process but neglected for the $gg$-initiated process the resulting exclusion region would have the “bay” shape indicated by the black line in Fig. 7. The “fjord” between $25\lesssim\tan\beta\lesssim 29$ would appear to be excluded in this case. The reason for this relatively large impact of the $gg$-initiated process in a parameter region where it is sub-dominant (see Fig. 5 and Fig. 6) is mainly due to the enhancement of the $gg$-initiated process in the resonance region that would occur for a non-zero phase if the interference contribution were neglected, see Fig. 5. While this parameter region is allowed by the coherent $b\bar{b}$ cross section (see Fig. 5), it would appear to be excluded by the incoherent enhanced $gg$ cross section in this case (see Fig. 5). On the other hand, incorporating the interference contribution only for the $gg$-initiated processes but not for the $b\bar{b}$-initiated one would only have a minor impact, as can be seen by comparing the grey line in Fig. 7 with the boundary of the blue area.

6 Conclusions

In the present paper we have investigated how the limits from the LHC searches for additional heavy SUSY Higgs bosons get modified if instead of the assumption of $\mathcal{CP}$ conservation in the Higgs sector the general case of complex parameters giving rise to $\mathcal{CP}$-violating effects is taken into account. We have considered neutral Higgs boson production in gluon fusion and in association with $b\bar{b}$, and the decay into a pair of $\tau$-leptons in an extension of the well-known $M_{h}^{\rm{mod+}}$ benchmark scenario of the MSSM where $\mathcal{CP}$-violating effects are induced by a non-zero phase of the trilinear coupling $A_{t}$, $\phi_{A_{t}}=\pi/4$.

While in the LHC searches for additional Higgs bosons in supersymmetric models so far it has been assumed that the contributions from $\mathcal{CP}$-even and $\mathcal{CP}$-odd states to the signal cross section can be added as an incoherent sum, such an assumption is not valid in the general case where $\mathcal{CP}$-violating effects in the Higgs sector are taken into account. In the Higgs sector of the MSSM, complex parameters enter via loop corrections, causing $\mathcal{CP}$-violating effects beyond the tree level. As a consequence, the $\mathcal{CP}$ eigenstates $h,H,A$ mix into the mass eigenstates $h_{1},h_{2},h_{3}$. The Higgs sectors of the NMSSM and of a non-supersymmetric 2HDM are in general $\mathcal{CP}$-violating already at lowest order. The most obvious realisation of an extended Higgs sector that is compatible with the present experimental constraints is a scenario where one neutral scalar is SM-like and has a mass of about 125 GeV, while the other neutral Higgs bosons are significantly heavier. In such a decoupling-type scenario the heavy neutral Higgs bosons of the extended Higgs sector are typically nearly mass-degenerate. If $\mathcal{CP}$-violating effects are present, potentially large mixing and interference effects between the heavy neutral Higgs bosons occur as a generic feature of such extended Higgs sectors. These kind of effects should therefore be taken into account in the interpretation of the limits from the searches for additional neutral Higgs bosons.

Using an approximation of the $(3\times 3)$ Higgs propagator matrix in terms of wave function normalisation factors evaluated at the complex poles and Breit–Wigner propagator factors, we have determined the relative interference contributions to the full process of production and decay. We have demonstrated how these interference contributions can be combined with the results for the on-shell production and decay of individual Higgs bosons in order to arrive at a prediction for the full process incorporating interference effects. We have compared this result based on a coherent sum of the relevant contributions with the corresponding incoherent sum that would be obtained if interference contributions were omitted.

Furthermore, we have analytically investigated the relative interference contribution with respect to its dependence on the couplings, wave function normalisation factors, masses and widths. We have pointed out that besides the well-know criterion in terms of the mass difference and sum of the total Higgs widths, the occurrence of a sizeable $\mathcal{CP}$-violating interference effect is also directly related to the presence of imaginary contributions in the propagator matrix. Since the latter can be well approximated by the wave function normalisation factors, we could derive a simple criterion indicating whether interference effects are expected to be relevant.

In the considered scenario we have found that both the mixing effects and the interference contributions are large in the resonance region where the above criteria are fulfilled. In fact, if interference effects were omitted, the mixing effect between the nearly mass-degenerate heavy neutral Higgs bosons $h_{2}$ and $h_{3}$ would give rise to a significant enhancement of the total cross section in this region. However, the large destructive interference contribution, which in the considered scenario reaches a minimum of $-97\%$ at $M_{H^{\pm}}\simeq 480\,\textrm{GeV}$ and $\tan\beta\simeq 29$ in both production processes, overcompensates this effect. The net effect is therefore a large suppression of the total cross section times branching ratio in the resonance region as compared to the case of real parameters.

We have studied the impact of those modifications from mixing and interference effects in the presence of complex parameters on the limits obtained from the LHC searches in comparison with the limits for the $\mathcal{CP}$-conserving case of real parameters. For illustration we have first investigated single cross section limits, followed by a more comprehensive analysis in the $(M_{H^{\pm}},\tan\beta)$ parameter plane of the $\mathbb{C}M_{h}^{\rm{mod+}}$ scenario using the tool HiggsBounds. As expected by our analytical investigation, we have found that both the mixing enhancement and the relative effect of the destructive interference are similar for the two production modes of gluon fusion and associated production with $b\bar{b}$. For moderate and large values of $\tan\beta$, the $b\bar{b}$ cross section is larger than the gluon fusion one, but both processes are relevant for the analysed parameters. We have demonstrated that the $\mathcal{CP}$-violating mixing and interference effects cause a substantial shift of the exclusion bounds in the parameter space along the “LHC wedge” region in comparison to the limits that would be obtained in the corresponding scenario with real parameters. We have shown that in order to obtain those results it is important to incorporate the interference contribution both for the $b\bar{b}$ associated production and the gluon fusion process. Our analysis has indicated that a considerable parameter region remains unexcluded by LHC searches from Run 1 that would appear to be ruled out if interference effects were neglected.

The illustrative study carried out in this paper motivates the analysis of current and future LHC results in scenarios with complex parameters, taking interference effects into account. While in our qualitative investigation we have restricted ourselves to the analysis of the experimental results that were obtained at Run 1 of the LHC and we have used an approximate treatment of the cross sections for the production processes, our results lay the foundation for a more detailed phenomenological analysis incorporating the latest experimental results and the most accurate theoretical predictions. We defer such a more detailed investigation taking into account the latest experimental results from Run 2 of the LHC, employing the recently released tool SusHiMi for the accurate prediction of the Higgs production processes in gluon fusion and $b\bar{b}$ associated production, and studying the variation of the phases of several relevant MSSM parameters, in particular the phase of the gluino mass parameter, to future work [65].

Acknowledgements

We would like to thank Oscar Stål for his help with HiggsBounds and valuable suggestions, and Christian Veelken, Sven Heinemeyer, Shruti Patel and Stefan Liebler for useful discussions. E.F. thanks the DESY theory group where a large part of this work was done. The work of E.F. was partially funded by the German National Academic Foundation and by DESY. The work of G.W. is supported in part by the Collaborative Research Centre SFB 676 of the DFG, “Particles, Strings and the Early Universe”, and by the European Commission through the “HiggsTools” Initial Training Network PITN-GA-2012-316704.