Neutral Higgs-pair production at linear colliders within the general two-Higgs-doublet model: Quantum effects and triple Higgs boson self-interactions

Abstract

The pairwise production of neutral Higgs bosons $({\mathrm{h}}^{0}{\mathrm{A}}^{0},{\mathrm{H}}^{0}{\mathrm{A}}^{0})$ is analyzed in the context of the future linear colliders, such as the ILC and CLIC, within the general two-Higgs-doublet model (2HDM). The corresponding cross sections are computed at the one-loop level, including the full set of contributions at order $\mathcal{O}({\ensuremath{\alpha}}_{\mathrm{ew}}^{3})$ together with the leading $\mathcal{O}({\ensuremath{\alpha}}_{\mathrm{ew}}^{4})$ terms, in full compliance with the current phenomenological bounds and the stringent theoretical constraints inherent to the consistency of the model. We uncover regions across the 2HDM parameter space, mainly for low $\mathrm{tan}\ensuremath{\beta}\ensuremath{\sim}1$ and moderate---and negative---values of the relevant ${\ensuremath{\lambda}}_{5}$ parameter, wherein the radiative corrections to the Higgs-pair production cross section, $\ensuremath{\sigma}({\mathrm{e}}^{+}{\mathrm{e}}^{\ensuremath{-}}\ensuremath{\rightarrow}{\mathrm{A}}^{0}{\mathrm{h}}^{0}/{\mathrm{A}}^{0}{\mathrm{H}}^{0})$, can comfortably reach $|\ensuremath{\delta}\ensuremath{\sigma}|/\ensuremath{\sigma}\ensuremath{\sim}50%$. This behavior can be traced back to the enhancement capabilities of the trilinear Higgs self-interactions---a trademark feature of the 2HDM, with no counterpart in the minimal supersymmetric standard model (MSSM). The corrections are strongly dependent on the actual value of ${\ensuremath{\lambda}}_{5}$ as well as on the Higgs mass spectrum. Interestingly enough, the quantum effects are positive for energies around $\sqrt{s}\ensuremath{\simeq}500\text{ }\text{ }\mathrm{GeV}$, thereby producing a significant enhancement in the expected number of events precisely around the fiducial startup energy of the ILC. The Higgs-pair production rates can be substantial, typically a few tens of femtobarn, therefore amounting to a few thousand events per $500\text{ }\text{ }{\mathrm{fb}}^{\ensuremath{-}1}$ of integrated luminosity. In contrast, the corrections are negative in the highest energy range (viz. $\sqrt{s}\ensuremath{\sim}1\text{ }\text{ }\mathrm{TeV}$ and above). We conclude that a precise measurement of the 2H final states could carry unambiguous footprints of an extended (nonsupersymmetric) Higgs sector. Finally, to better assess the scope of these effects, we compare the exclusive pairwise production of Higgs bosons with the inclusive gauge boson fusion channels leading to $2H+X$ final states, and also with the exclusive triple Higgs boson production. We find that these multiparticle final states can be highly complementary in the overall Higgs bosons search strategy.