two-Higgs-doublet Model Research Articles

Accommodating the LHC charged Higgs boson excess at 130 GeV in the general Two-Higgs Doublet Model

Abstract Charged Higgs bosons are common predictions in most extensions of the Standard Model (SM) Higgs sector. Therefore, their observation would elucidate the nature of the Higgs sector. Motivated by the ATLAS collaboration's latest analysis performed with $139~\text{fb}^{-1}$ of Run 2 data intended to search for charged Higgs boson, produced in top quark decay and subsequently decaying via $H^\pm \rightarrow cb$, where an excess with a local significance of $3\sigma$ is observed at $m_{H^\pm} = 130~\rm{GeV}$, we discuss here the possibility of explaining such excess in the context of the general 2-Higgs Doublet Model (2HDM type-III), after satisfying all theoretical and up-to-date experimental constraints. We also propose phenomenological scenarios to further explore the mass region around 130 GeV in the four Yukawa types of the 2HDM type-III and suggest alternative decay channel $H^\pm \rightarrow cs$ and/or $H^\pm \rightarrow W^{\pm *}h$ to probe the nature of the observed excess (if it is not a statistical fluctuation). Future searches for $H^\pm$ will be critical in confirming or refuting the first hint of a light charged Higgs boson at the LHC.

Two-Higgs-doublet model effective field theory

We construct the general two-Higgs doublet model effective field theory where the effects of additional new physics are parametrized by operators up to mass dimension-six. We further transform this effective theory to the Higgs basis and provide matching of the Wilson coefficients between the two descriptions. We illustrate the advantages of the Higgs basis which include the separation of operators that modify standard model couplings and masses from operators that contribute to scattering processes only, transparent correlations between scattering processes resulting from the same operator, and derivation of correlations between different operators in specific UV completions. For completeness, we also construct specific versions corresponding to four types of two-Higgs-doublet models: type-I, -II, -X, and -Y, distinguished by Z2 symmetries which restrict the couplings of the Higgs doublets to standard model fermions. Furthermore, we derive general vacuum and stability conditions of the scalar potential in the presence of higher-dimensional terms. Published by the American Physical Society 2024

Lessons from LHC on the LFV Higgs decays h→ℓaℓb in the two-Higgs doublet models

Lessons from LHC on the LFV Higgs decays h→ℓaℓb in the two-Higgs doublet models

Neutrino Phenomenology and keV Dark Matter in the Two-Higgs Doublet Model with A4 Symmetry

Abstract We propose a minimal extended seesaw scheme based on the discrete symmetry $A_4\times Z_4\times Z_2\times Z_8$, which can successfully address neutrino phenomenology and keV sterile neutrino dark matter. The lepton mass hierarchy is naturally achieved. Active neutrino mixing angles can reached the best-fit points with the predictive Dirac CP violation phase. The active–sterile mixing matrix elements are small enough to access the observed cosmological dark matter abundance constraint with keV sterile neutrino dark matter. The effective neutrino masses are predicted to be in the ranges of the recent experimental limits.

Maximising CP Violation in naturally aligned Two-Higgs Doublet Models

The Two-Higgs Doublet Model (2HDM) is a well-motivated theoretical framework that provides additional sources of CP Violation (CPV) beyond the Standard Model (SM). After studying the vacuum topology of a general (convex) 2HDM potential, we unambiguously identify three origins of CPV: (i) Spontaneous CPV (SCPV), where the vacuum manifold has at least two degenerate CPV minima disconnected by domain walls, (ii) Explicit CPV (ECPV) with one single CPV ground state, and (iii) Mixed Spontaneous and Explicit CPV (MCPV), where the theory possesses more than one non-degenerate CPV local minimum. Most importantly, we define a novel complex parameter rCP whose norm and phase control the three different realisations of CPV, at least at the tree level. In all these scenarios, only two CPV phases can be made independent, as any third CPV parameter will always be constrained via the CP-odd tadpole condition. Since ECPV vanishes in 2HDMs where SM Higgs alignment is achieved naturally through accidental continuous symmetries, we analyse the possibility of maximising CPV through soft and explicit breaking of these symmetries. We derive upper limits on key CPV parameters that quantify the degree of SM misalignment from constraints due to the non-observation of an electron Electric Dipole Moment (EDM). Finally, we delineate the CP-violating parameter space of the so-constrained naturally aligned 2HDMs that can further be probed at the CERN Large Hadron Collider (LHC).

Impact of new experimental data on the C2HDM: the strong interdependence between LHC Higgs data and the electron EDM

The complex two-Higgs doublet model (C2HDM) is one of the simplest extensions of the Standard Model with a source of CP-violation in the scalar sector. It has a \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathbb{Z}}_{2}$$\\end{document} symmetry, softly broken by a complex coefficient. There are four ways to implement this symmetry in the fermion sector, leading to models known as Type-I, Type-II, Lepton-Specific and Flipped. In the latter three models, there is a priori the surprising possibility that the 125 GeV Higgs boson couples mostly as a scalar to top quarks, while it couples mostly as a pseudoscalar to bottom quarks. This “maximal” scenario was still possible with the data available in 2017. Since then, there have been more data on the 125 GeV Higgs boson, direct searches for CP-violation in angular correlations of τ-leptons produced in Higgs boson decays, new results on the electron electric dipole moment, new constraints from LHC searches for additional Higgs bosons and new results on b → sγ transitions. Highlighting the crucial importance of the physics results of LHC’s Run 2, we combine all these experiments and show that the “maximal” scenario is now excluded in all models. Still, one can have a pseudoscalar component in \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h\ au \\overline{\ au }$$\\end{document} couplings in the Lepton-Specific case as large as 87% of the scalar component for all mass orderings of the neutral scalar bosons.

H-COUP Version 3: A program for one-loop corrected decays of any Higgs bosons in non-minimal Higgs models

The H-COUP program is provided as a package of Fortran codes, which can compute observables related to Higgs bosons including radiative corrections in various extended Higgs sectors. We give a manual for the latest version of H-COUP (H-COUP_3.0), in which decay rates and branching ratios of all the Higgs bosons can be calculated at one-loop level in EW and Higgs interactions with QCD corrections in the Higgs singlet model, four types of the two Higgs doublet model with a softly-broken Z2 symmetry, and the inert doublet model. The previous version (H-COUP_2.0) can evaluate those only for the standard model like Higgs boson with the mass of 125 GeV (h). In H-COUP_3.0, renormalized quantities are computed based on the gauge-independent on-shell renormalization scheme. The source code of H-COUP_3.0 can be downloaded via the following link: http://www-het.phys.sci.osaka-u.ac.jp/~hcoup. By using H-COUP_3.0, we can compare the precise measurements of the properties of h and direct searches for additional Higgs bosons with their predictions at one-loop level, by which we can reconstruct the structure of the Higgs sector. New version program summaryProgram Title: H-COUP Version 3 5CPC Library link to program files:https://doi.org/10.17632/f88szmrj5x.3Developer's repository link:http://www-het.phys.sci.osaka-u.ac.jp/~hcoupLicensing provisions: GPLv3Programming language: Fortran90Journal reference of previous version: Comput. Phys. Commum. 257 (2020) 107512Does the new version supersede the previous version?: YesReasons for the new version: The previous version (H-COUP version 2) is a package of Fortran codes to provide one-loop corrected couplings, decay rates and branching ratios of the 125 GeV Higgs boson in the Higgs Singlet Model (HSM), the Two Higgs doublet models (THDMs) with a softly-broken Z2 symmetry, and the Inert Doublet Model (IDM). In the new version (H-COUP version 3), the H-COUP program has been improved to calculate the decay rates and the branching ratios of all the Higgs bosons including one-loop corrections in the above models. These calculations are performed based on the gauge independent on-shell renormalization scheme.Summary of revisions: We have added radiatively corrected decay rates and branching ratios for all the additional Higgs bosons, not only for the 125 GeV Higgs boson, in the extended Higgs models. Furthermore, several choices are provided for renormalization schemes used to evaluate electroweak corrections. The scalar two-point functions are renormalized by on-shell scheme as in the previous version, while the renormalization for the tadpole is performed by either the standard tadpole scheme or the alternative tadpole scheme. The users can choose these renormalization schemes.Nature of problem: Decay rates and branching ratios of all the Higgs bosons can be calculated with one-loop electroweak corrections and QCD corrections in the HSM, four types of the THDMs with the softly-broken Z2 symmetry, and the IDM. In the calculations of additional Higgs boson decays, for the three-body decay processes and loop induce processes of additional Higgs bosons, electroweak corrections are not included. As for the QCD correction, next-to-leading order or next-to-next-to-leading order corrections are included, depending on the decay process.Solution method: Electroweak corrections are evaluated based on the gauge independent on-shell renormalization scheme, while some parameters are renormalized by the MS‾ scheme. QCD corrections are evaluated with the MS‾ scheme.Additional comments including restrictions and unusual features: All functions up to the previous version are included in H-COUP version 3.

Relevance of one-loop SMEFT matching in the 2HDM

The two-Higgs doublet model (2HDM) is a well understood alternative to the Standard Model of particle physics. If the new particles included in the 2HDM are at an energy scale much greater than the weak scale, the theory can be matched to the Standard Model effective field theory (SMEFT). We compute for the first time the complete one-loop matching at dimension-6. We compare its numerical impact with that of tree-level matching at dimension-8 by performing a global fit to single Higgs and precision electroweak measurements, and we emphasize the importance of comparing one-loop SMEFT results with corresponding one-loop results in the full 2HDM model. In the SMEFT, we consider the relative importance of both one-loop matching and the inclusion of renormalization group evolution. Our results demonstrate the necessity of studying the impact of various expansions to quantify the uncertainties of the SMEFT matching. Published by the American Physical Society 2024

95.4 GeV diphoton excess at ATLAS and CMS

The ATLAS Collaboration has recently reported the results of a low-mass Higgs-boson search in the diphoton final state based on the full Run 2 dataset. The results are based on an improved analysis with respect to the previous analysis, which included a part of the Run 2 data, with a substantially better sensitivity. The “model-dependent” search carried out by ATLAS shows an excess of events at a mass of about 95.4 GeV with a local significance of 1.7σ. The results are compatible with a previously reported excess at the same mass, but with somewhat higher significance of 2.9σ, from the CMS collaboration, also based on the full Run 2 dataset. Combining the two results (neglecting possible correlations), we find a signal strength of μγγATLAS+CMS=0.24−0.08+0.09, corresponding to an excess of 3.1σ. In this work, we investigate the implications of this result, updating a previous analysis based solely on the CMS Run 2 data. We demonstrate that the ATLAS/CMS combined diphoton excess can be interpreted as the lightest Higgs boson in a Two-Higgs doublet model that is extended by a complex singlet (S2HDM) of Yukawa types II and IV, while being in agreement with all other experimental and theoretical constraints. Published by the American Physical Society 2024

The CMS di-photon excess at 95 GeV in view of the LHC Run 2 results

The CMS collaboration has recently reported the results of a low-mass Higgs-boson search in the di-photon final state based on the full Run 2 data set with refined analysis techniques. The new results show an excess of events at a mass of about 95 GeV with a local significance of 2.9σ, confirming a previously reported excess at about the same mass and similar significance based on the first-year Run 2 plus Run 1 data. The observed excess is compatible with the limits obtained in the corresponding ATLAS searches. In this work, we discuss the di-photon excess and show that it can be interpreted as the lightest Higgs boson in the Two-Higgs doublet model that is extended by a complex singlet (S2HDM) of Yukawa types II and IV. We show that the second-lightest Higgs boson is in good agreement with the current LHC Higgs-boson measurements of the state at 125 GeV, and that the full scalar sector is compatible with all theoretical and experimental constraints. Furthermore, we discuss the di-photon excess in conjunction with an excess in the bb¯ final state observed at LEP and an excess observed by CMS in the di-tau final state, which were found at comparable masses with local significances of about 2σ and 3σ, respectively. We find that the bb¯ excess can be well described together with the di-photon excess in both types of the S2HDM. However, the di-tau excess can only be accommodated at the level of 1σ in type IV. We also comment on the compatibility with supersymmetric scenarios and other extended Higgs sectors, and we discuss how the potential signal can be further analyzed at the LHC and at future e+e− colliders.

Muon g-2 and W-mass in a framework of colored scalars: an LHC perspective

A color octet isodoublet can have esoteric origins and it complies with minimal flavour violation. In this study, we take a scenario where the well known Type-X Two-Higgs doublet model is augmented with a color octet isodoublet. We shed light on how such a setup can predict the recently observed value for the W-boson mass. The two-loop Barr-Zee contributions to muon g-2 stemming from the colored scalars are evaluated. It is subsequently found that the parameter space compatible with the observed muon g-2 gets relaxed w.r.t. what it is in the pure Type-X 2HDM by virtue of the contribution from the colored scalars. The extended parameter region therefore successfully accounts for both the W-mass and muon g-2 anomalies simultaneously. Finally, a collider signature leading to a tau ^+ tau ^- b overline{b} final state is explored at the 14 TeV LHC using both cut-based and multivariate techniques. Such a signal can confirm the existence of both colorless as well colored scalars that are introduced by this framework.

Global monopoles in the two-Higgs-doublet-model

We discuss monopoles formed due to the spontaneous breakdown of a global SO(3)HF symmetry within the global two-Higgs doublet model. We explain that the Higgs sector dynamics can be described in terms of two vectors one of which is null, RA=(R0,Ra,R4,R5) for a=1,2,3, with 5 independent components describing the Higgs family symmetry and another, na, with 3 independent components related to the “would-be” Goldstone bosons. When formed from random initial conditions we find that monopoles are formed with a charged vacuum in the centre which couples the two fields together. We find a spherically symmetric solution which is an approximately uniform, unit winding of the sphere in both the Ra and na vectors. These global monopoles are closely related to the Nambu monopole. The additional complexity and structure contained in these monopoles does not appear to prevent the scaling of their density.

Gildener-Weinberg two-Higgs-doublet model at two loops

The Gildener-Weinberg two-Higgs doublet model provides a naturally light and aligned Higgs boson, $H=H(125)$. It has been studied in the one-loop approximation of its effective potential, ${V}_{1}$. An important consequence is that the masses of the model's beyond the Standard Model (BSM) Higgs bosons (${H}^{\ensuremath{'}},A,{H}^{\ifmmode\pm\else\textpm\fi{}}$) are bounded by the sum rule ${({M}_{{H}^{\ensuremath{'}}}^{4}+{M}_{A}^{4}+2{M}_{{H}^{\ifmmode\pm\else\textpm\fi{}}}^{4})}^{1/4}=540\text{ }\text{ }\mathrm{GeV}$. Although they are well within reach of the LHC, searches for them have been stymied by large QCD backgrounds. Another consequence is that $H$ is highly aligned, i.e., $H--{H}^{\ensuremath{'}}$ mixing is small and $H$ has only Standard Model couplings. A corollary of this alignment is that commonly pursued discovery modes such as ${H}^{\ensuremath{'}}$, $A\ensuremath{\leftrightarrow}{W}^{+}{W}^{\ensuremath{-}}$, $ZZ$, $HZ$, and ${H}^{\ifmmode\pm\else\textpm\fi{}}\ensuremath{\leftrightarrow}{W}^{\ifmmode\pm\else\textpm\fi{}}Z,{W}^{\ifmmode\pm\else\textpm\fi{}}H$ are beyond the reach of LHC experiments. To assess the accuracy of the sum rule and Higgs alignment, we study this model in two loops. This calculation is complicated by having many new contributions. We present two formulations of it to calculate the $H--{H}^{\ensuremath{'}}$ mass matrix, its eigenvectors ${H}_{1}$, ${H}_{2}$, and the mass ${M}_{{H}_{2}}$ while fixing ${M}_{{H}_{1}}=125\text{ }\text{ }\mathrm{GeV}$. They give similar results and are in accord with the one-loop results. Requiring ${M}_{A}={M}_{{H}^{\ifmmode\pm\else\textpm\fi{}}}$, we find $180\text{ }\text{ }\mathrm{GeV}\ensuremath{\lesssim}{M}_{A,{H}^{\ifmmode\pm\else\textpm\fi{}}}\ensuremath{\lesssim}380--425\text{ }\text{ }\mathrm{GeV}$ and $550--700\text{ }\text{ }\mathrm{GeV}\ensuremath{\gtrsim}{M}_{{H}_{2}}\ensuremath{\gtrsim}125\text{ }\text{ }\mathrm{GeV}$, with ${M}_{{H}_{2}}$ decreasing as ${M}_{A,{H}^{\ifmmode\pm\else\textpm\fi{}}}$ increase. The corrections to $H$ alignment are below $\mathcal{O}(1%)$. So, the BSM searches above will remain fruitless. Finding the BSM Higgses requires improved sensitivity to their low masses. We discuss three possible searches for this.

Chasing the two-Higgs doublet model in the di-Higgs boson production

We investigate the di-Higgs production at the Large Hadron Collider in the two-Higgs doublet model (2HDM). In particular, we study the production of an extra neutral Higgs boson $\phi$ in association with the Standard Model (SM) Higgs boson $h$ in the Higgs alignment limit. We analyze two scenarios where the additional Higgs $\phi$ is CP-even or -odd state with a large top-Yukawa interaction. The leading contribution of this production comes from the top-quark loop-induced gluon-fusion channel $gg \to h\phi$. The measurement of the $h\phi$ production can probe the quartic couplings in the Higgs potential as well as the top-Yukawa couplings. Imposing both theoretical constraints (from the perturbative unitarity and the vacuum stability bounds) and experimental bounds (from the SM Higgs and flavor physics measurements) on the 2HDM parameter space, we calculate the production cross-section of $gg \to h\phi$. Furthermore, we scrutinize these processes in the parameter spaces where the CMS di-tau and di-photon excesses around 100$\,$GeV, and/or the muon $g-2$ anomaly can be accommodated.

Constraining lepton flavor violating Higgs couplings at the HL-LHC in the vector boson fusion channel

We explore the parameter space of lepton flavor violating (LFV) neutral Higgs Yukawa couplings with the muon and tau leptons that can be probed at the high-luminosity Large Hadron Collider (HL-LHC) via the vector boson fusion (VBF) Higgs production process. Our projected sensitivities for the Standard Model Higgs ($h$) LFV branching ratio $\text{Br}(h\ensuremath{\rightarrow}\ensuremath{\mu}\ensuremath{\tau})$ in the $pp\ensuremath{\rightarrow}hjj\ensuremath{\rightarrow}(h\ensuremath{\rightarrow}\ensuremath{\mu}\ensuremath{\tau})jj$ channel at the HL-LHC are contrasted with the current and future low-energy constraints from the anomalous magnetic moment and electric dipole moment of the muon, as well as with other LFV observables, such as $\ensuremath{\tau}\ensuremath{\rightarrow}3\ensuremath{\mu}$ and $\ensuremath{\tau}\ensuremath{\rightarrow}\ensuremath{\mu}\ensuremath{\gamma}$. We also study the LFV prospects of a generic beyond the Standard Model neutral Higgs boson ($H$) with a mass in the range of ${m}_{H}\ensuremath{\in}[20,800]\text{ }\text{ }\mathrm{GeV}$ and give the projected model-independent upper limits on the VBF production cross section of $Hjj$ times the branching ratio of $H\ensuremath{\rightarrow}\ensuremath{\mu}\ensuremath{\tau}$ at the HL-LHC. We interpret these results in the context of a two-Higgs doublet model as a case study.

Interplay of the charged Higgs boson effects in RD(*) , b→sℓ+ℓ− , and W mass

Current data on semileptonic charged- and neutral-current $B$ decays show deviations from the predictions of the Standard Model. It is well known that a charged Higgs boson, belonging to the two-Higgs doublet model without $Z_2$ symmetry, offers one of the simplest solution to the charged-current $B$ decays. We show that this solution naturally induces a negative shift of $\mathcal{O}(1)$ in the Wilson coefficient ($C_{9\ell}$) of operator $(\bar s_L\gamma_\mu b_L)(\bar \ell\gamma^\mu \ell)$, potentially resolving the tension in neutral-current $B$ decays as well. Interestingly, the lepton universality ratios in $b\to s \ell^+\ell^-$ decays, in tune with the recent LHCb result, remain SM-like. Precision constraints from neutral $B$ and $K$ meson mixing, decays $B_c\to \tau\bar\nu$, $B\to X_s\gamma$, and leptonic decays of $\tau$ and $Z$ can be satisfied. Furthermore, a positive shift in $W$-boson mass, nicely in agreement with the CDF measurement, is also possible, requiring the neutral scalars to be heavier than the charged Higgs but within the sub-TeV region.

Muon g−2 in a type-X 2HDM assisted by inert scalars: A test at the ILC

A two-Higgs doublet model can predict the observed muon $g\ensuremath{-}2$ for an appropriately light pseudoscalar that now faces tight constraints. However, it was shown recently in the companion paper [N. Chakrabarty, Muon $g\ensuremath{-}2$ in a type-X 2HDM assisted by inert scalars: A test at the LHC, preceding paper, arXiv:2112.13126.] that augmenting the two-Higgs doublet model by an additional inert doublet can lead to an explanation to the muon $g\ensuremath{-}2$ anomaly for a much heavier pseudoscalar. In this study, we probe such a framework at the proposed International Linear Collider using beam polarization for $\sqrt{s}=1\text{ }\text{ }\mathrm{TeV}$. Using multivariate techniques, we analyze the signals ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{\ensuremath{\tau}}^{+}{\ensuremath{\tau}}^{\ensuremath{-}}+$ missing transverse energy and ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\rightarrow}{\ensuremath{\mu}}^{+}{\ensuremath{\mu}}^{\ensuremath{-}}+$ missing transverse energy in the lepton- and muon-specific versions of the framework, respectively. Our analysis reveals that the ${e}^{+}{e}^{\ensuremath{-}}$ machine operating at a $3000\text{ }\text{ }{\mathrm{fb}}^{\ensuremath{-}1}$ luminosity predicts a $5\ensuremath{\sigma}$ discovery of a pseudoscalar as heavy as 400 GeV. Comparing with the companion paper, it is concluded that the International Linear Collider is a much more potent machine than the LHC in this regard.

Muon g−2 in a type-X 2HDM assisted by inert scalars: A test at the LHC

A scenario augmenting the well-known type-X two-Higgs doublet model with an additional inert doublet is proposed. The type-X two-Higgs doublet model is known to offer a solution to the muon $g\ensuremath{-}2$ anomaly for a light pseudoscalar. We show that the proposed framework can accommodate a heavier pseudoscalar on account of two-loop Barr-Zee contributions to muon $g\ensuremath{-}2$ stemming from the inert doublet. We subsequently explore an interesting ${\ensuremath{\tau}}^{+}{\ensuremath{\tau}}^{\ensuremath{-}}+$ missing transverse energy signal that can be used to probe the present scenario at the 14 TeV LHC.

Is the new CDF M measurement consistent with the two-Higgs doublet model?

Motivated by the new CDF measurement of the $W$ boson mass reported recently which clearly illustrates a large deviation compared to the Standard Model (SM) prediction. In the present paper, we study the Two-Higgs Doublet Model (2HDM) contributions to $M_W$ and its phenomenological implications in the case where the heavy CP-even $H$ is identified as the observed Higgs boson with a mass of $125$ GeV. Taking into account theoretical and all the available experimental constraints as well as the new CDF measurement, we demonstrate that the 2HDM parameter space can provide a large correction which predicts the W mass close to the new CDF $M_W$ measurement. It is found that $M_{H\pm}=M_A$ is excluded and the splitting of the charged Higgs boson with all other states is positive. We also discuss the consequence on the effective mixing angle $\sin^2\theta_{\text{eff}}$ as well as the phenomenological implications on the charged Higgs and CP-odd decays.

Higgs-like particle decays into γZ and γγ: Fingerprints of some non-supersymmetric models

Recently, ATLAS and CMS experiments at the LHC put on light the relevant results in the measurement precision of the Higgs and BSM. In such a report, where the resonance direct search was made in the γZ channel, a mass adjustment distribution for the reconstructed Z boson and photon was established. Thus, simultaneously with the signal-from-background separation, the number of events has been perfectly described, and an excess signal approximately twice that expected by the Standard Model (SM) has been noticed, which is equivalent to a significance of 2.2 standard deviations. In this study, we examine how any possible new physics models can explain this excess, such as the CP-conserving Two-Higgs doublet model (2HDMs), the Inert doublet model (IDM), and the Higgs triplet model (HTM). While considering the available theoretical and most recent experimental constraints, the phenomenological implications of existing extensions beyond the SM are discussed, and prospects for precision studies of these processes are described. We have found that the excess could be explained in the BSM frameworks studied, depending on the charged and double-charged Higgs boson masses.