Larbi Rahili
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Research Interests
Larbi Rahili has a total of 43 co-authors. Here are the top co-authors by number of shared publications:
- Abdesslam Arhrib
Université Abdelmalek Essaadi FSTT B. 416 Tangier Morocco
11 shared publications - Mohamed Chabab
Cadi Ayyad University
10 shared publications - Rachid Benbrik
Cadi Ayyad University
8 shared publications - Es-Said Ghourmin
Ibnou Zohr University
6 shared publications - Stefano Moretti
University of Southampton
6 shared publications - Ayoub Hmissou
Faculty of Science Laboratory of Theoretical and High Energy Physics Ibnou Zohr University B.P 8106 Agadir Morocco
5 shared publications
Latest Publications
- Open Access
- 7 July 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 114 (1)
Probing standard-model-like di-Higgs production at photon-photon colliders in the I ( 1 + 2 ) HDM type-I
In this paper, pair production of standard model (SM)-like Higgs bosons, h h , is studied through γ γ scattering at future electron-positron colliders, in the framework of the inert doublet model with two active doublets—i.e., the I( 1 + 2 )HDM for short. The relevance of the process γ γ → h h for such a beyond-the-SM (BSM) scenario stems from the fact that it is a one-loop process at lowest order, wherein inert charged states χ ± contribute alongside W ± , H ± , and heavy fermions (primarily, bottom and top quarks), crucially, at the same perturbative order. Given that χ ± / H ± masses and h S + S − ( S ± = χ ± , H ± ) couplings are very mildly constrained, there exist regions of the parameter space of the I ( 1 + 2 ) HDM where the former can be rather light and the latter rather large. After imposing up-to-date theoretical and experimental constraints on the I ( 1 + 2 ) HDM , it is found that the production rates of such processes at future γ γ machines can be enhanced by up to a factor of ≈ 50 with respect to the SM, significantly exceeding typical yields of conventional two-Higgs-doublet models (2HDMs). Further, thanks to the level of control that one can attain at such facilities on the photon kinematics, leading to excellent invariant mass resolution of the incoming photon pairs, we show how it is possible to extract from this process the value of the χ ± mass (along with those of the active H ± states) with high precision, and whichever of the decays of the h h pair, both with and without beam polarization.
- 7 July 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 114 (1)
Open AccessProbing standard-model-like di-Higgs production at photon-photon colliders in the I ( 1 + 2 ) HDM type-I
In this paper, pair production of standard model (SM)-like Higgs bosons, h h , is studied through γ γ scattering at future electron-positron colliders, in the framework of the inert doublet model with two active doublets—i.e., the I( 1 + 2 )HDM for short. The relevance of the process γ γ → h h for such a beyond-the-SM (BSM) scenario stems from the fact that it is a one-loop process at lowest order, wherein inert charged states χ ± contribute alongside W ± , H ± , and heavy fermions (primarily, bottom and top quarks), crucially, at the same perturbative order. Given that χ ± / H ± masses and h S + S − ( S ± = χ ± , H ± ) couplings are very mildly constrained, there exist regions of the parameter space of the I ( 1 + 2 ) HDM where the former can be rather light and the latter rather large. After imposing up-to-date theoretical and experimental constraints on the I ( 1 + 2 ) HDM , it is found that the production rates of such processes at future γ γ machines can be enhanced by up to a factor of ≈ 50 with respect to the SM, significantly exceeding typical yields of conventional two-Higgs-doublet models (2HDMs). Further, thanks to the level of control that one can attain at such facilities on the photon kinematics, leading to excellent invariant mass resolution of the incoming photon pairs, we show how it is possible to extract from this process the value of the χ ± mass (along with those of the active H ± states) with high precision, and whichever of the decays of the h h pair, both with and without beam polarization.
- Open Access
- 18 March 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (5)
Interpreting the 650 GeV and 95 GeV Higgs anomalies in the next-to-two-Higgs-doublet model
Recent experimental hints from the Large Hadron Collider (LHC) in diphoton and partially in the τ + τ − final states suggest the possible existence of an additional Higgs boson with a mass around 95 GeV. Interestingly, these observations are consistent with earlier results from the Large Electron-Positron collider, which pointed to an excess in b b ¯ final states within a similar mass range. Additionally, CMS has observed an excess in the γ γ b b ¯ final state, indicating a possible resonance near 650 GeV decaying into a pair of Standard Model (SM)-like Higgs bosons or into a SM-like Higgs boson accompanied by a lighter scalar with mass near 95 GeV. In this work, we investigate whether these anomalies can be simultaneously explained within the next-to-two-Higgs-doublet model (N2HDM), an extension of the SM scalar sector featuring two complex Higgs doublets and an additional real singlet. Assuming the existence of a C P -even Higgs state compatible with the 95 GeV excesses (restricted to the γ γ and b b ¯ channels), we analyze the Type-II and Type-Y Yukawa structures, taking the observed 650 GeV resonance to be a C P -even Higgs state. An extensive parameter scan is performed, incorporating the latest constraints from the properties of the observed 125 GeV Higgs boson, direct searches for additional Higgs states, flavor physics data, and electroweak precision observables. Our results show that a heavy C P -even Higgs resonance around 650 GeV, produced predominantly via gluon-gluon fusion and subsequently decaying into a 125 GeV Higgs boson together with another scalar at approximately 95 GeV, can be simultaneously accommodated within both the N2HDM Type-II and Type-Y frameworks in parameter regions that remain consistent with the relevant experimental 2 σ intervals for the reported excesses, once all theoretical and experimental constraints are imposed. This interpretation leads to distinctive and testable predictions for the ongoing LHC Run 3 and the forthcoming high-luminosity LHC phase, in particular through correlated rates in the γ γ b b ¯ , τ + τ − b b ¯ , b b ¯ γ γ , and γ γ τ + τ − final states.
- 18 March 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (5)
Open AccessInterpreting the 650 GeV and 95 GeV Higgs anomalies in the next-to-two-Higgs-doublet model
Recent experimental hints from the Large Hadron Collider (LHC) in diphoton and partially in the τ + τ − final states suggest the possible existence of an additional Higgs boson with a mass around 95 GeV. Interestingly, these observations are consistent with earlier results from the Large Electron-Positron collider, which pointed to an excess in b b ¯ final states within a similar mass range. Additionally, CMS has observed an excess in the γ γ b b ¯ final state, indicating a possible resonance near 650 GeV decaying into a pair of Standard Model (SM)-like Higgs bosons or into a SM-like Higgs boson accompanied by a lighter scalar with mass near 95 GeV. In this work, we investigate whether these anomalies can be simultaneously explained within the next-to-two-Higgs-doublet model (N2HDM), an extension of the SM scalar sector featuring two complex Higgs doublets and an additional real singlet. Assuming the existence of a C P -even Higgs state compatible with the 95 GeV excesses (restricted to the γ γ and b b ¯ channels), we analyze the Type-II and Type-Y Yukawa structures, taking the observed 650 GeV resonance to be a C P -even Higgs state. An extensive parameter scan is performed, incorporating the latest constraints from the properties of the observed 125 GeV Higgs boson, direct searches for additional Higgs states, flavor physics data, and electroweak precision observables. Our results show that a heavy C P -even Higgs resonance around 650 GeV, produced predominantly via gluon-gluon fusion and subsequently decaying into a 125 GeV Higgs boson together with another scalar at approximately 95 GeV, can be simultaneously accommodated within both the N2HDM Type-II and Type-Y frameworks in parameter regions that remain consistent with the relevant experimental 2 σ intervals for the reported excesses, once all theoretical and experimental constraints are imposed. This interpretation leads to distinctive and testable predictions for the ongoing LHC Run 3 and the forthcoming high-luminosity LHC phase, in particular through correlated rates in the γ γ b b ¯ , τ + τ − b b ¯ , b b ¯ γ γ , and γ γ τ + τ − final states.
- Open Access
- 20 January 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (1)
Investigating the 95 GeV Higgs boson excesses within the type-I ( 1 + 2 ) HDM
In this work, we explore how the 2-Higgs doublet model (2HDM) type-I, extended by an inert doublet, can provide an explanation for the recently observed excesses at the Large Hadron Collider (LHC) in the γ γ and τ + τ − final states. Hence, by imposing theoretical constraints and experimental bounds on the model parameter space, our findings show that a light C P -even Higgs boson, h , with a mass around 95 GeV, can account for these anomalies. This result aligns with the excess in b b ¯ signatures reported in earlier data from the Large Electron-Positron (LEP) collider.
- 20 January 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (1)
Open AccessInvestigating the 95 GeV Higgs boson excesses within the type-I ( 1 + 2 ) HDM
In this work, we explore how the 2-Higgs doublet model (2HDM) type-I, extended by an inert doublet, can provide an explanation for the recently observed excesses at the Large Hadron Collider (LHC) in the γ γ and τ + τ − final states. Hence, by imposing theoretical constraints and experimental bounds on the model parameter space, our findings show that a light C P -even Higgs boson, h , with a mass around 95 GeV, can account for these anomalies. This result aligns with the excess in b b ¯ signatures reported in earlier data from the Large Electron-Positron (LEP) collider.
- Open Access
Probing Standard Model-like di-Higgs Production at Photon-Photon Colliders in the I(1+2)HDM Type-I
In this paper, pair production of Standard Model (SM)-like Higgs bosons, $hh$, is studied through $γγ$ scattering at future electron-positron colliders, in the framework of the Inert Doublet Model with two Active Doublets, i.e., the I(1+2)HDM for short. The relevance of the process $γγ\to hh$ for such a Beyond the SM (BSM) scenario stems from the fact that it is a one-loop process at lowest order, wherein inert charged states $χ^\pm$ contribute alongside with $W^\pm$, $H^\pm$ and heavy fermions (primarily, bottom and top quarks), crucially, at the same perturbative order. {Given that $χ^\pm/H^\pm$ masses and $hS^+S^-$ ($S^\pm=χ^\pm, H^\pm$) couplings are very mildly constrained,} there exist regions of the parameter space of the I(1+2)HDM where the former can be rather light and the latter rather large. After imposing up-to-date theoretical and experimental constraints on the I(1+2)HDM, it is found that the production rates of such process at future $γγ$ machines can be enhanced up to a factor of $\approx$ $50$ with respect to the SM, significantly exceeding typical yields of conventional 2-Higgs Doublet Models (2HDMs). Further, thanks to the level of control that one can attain at such facilities on the photon kinematics, leading to excellent invariant mass resolution of the incoming photon pairs, we show how it is possible to extract from this process the value of the $χ^\pm$ mass (along that of the active $H^\pm$ states) with high precision, whichever the decays of the $hh$ pair, both with and without beam polarization.Open AccessProbing Standard Model-like di-Higgs Production at Photon-Photon Colliders in the I(1+2)HDM Type-I
In this paper, pair production of Standard Model (SM)-like Higgs bosons, $hh$, is studied through $γγ$ scattering at future electron-positron colliders, in the framework of the Inert Doublet Model with two Active Doublets, i.e., the I(1+2)HDM for short. The relevance of the process $γγ\to hh$ for such a Beyond the SM (BSM) scenario stems from the fact that it is a one-loop process at lowest order, wherein inert charged states $χ^\pm$ contribute alongside with $W^\pm$, $H^\pm$ and heavy fermions (primarily, bottom and top quarks), crucially, at the same perturbative order. {Given that $χ^\pm/H^\pm$ masses and $hS^+S^-$ ($S^\pm=χ^\pm, H^\pm$) couplings are very mildly constrained,} there exist regions of the parameter space of the I(1+2)HDM where the former can be rather light and the latter rather large. After imposing up-to-date theoretical and experimental constraints on the I(1+2)HDM, it is found that the production rates of such process at future $γγ$ machines can be enhanced up to a factor of $\approx$ $50$ with respect to the SM, significantly exceeding typical yields of conventional 2-Higgs Doublet Models (2HDMs). Further, thanks to the level of control that one can attain at such facilities on the photon kinematics, leading to excellent invariant mass resolution of the incoming photon pairs, we show how it is possible to extract from this process the value of the $χ^\pm$ mass (along that of the active $H^\pm$ states) with high precision, whichever the decays of the $hh$ pair, both with and without beam polarization. - Open Access
- 27 October 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (7)
Could the 650 GeV excess be a pseudoscalar of a three-Higgs-doublet model?
In this study, we propose the interpretation of a 650 GeV excess observed at the Large Hadron Collider (LHC) by the CMS Collaboration in terms of the production of a C P -odd (or pseudoscalar) Higgs boson A , with mass around 650 GeV, decaying into the Standard Model (SM)-like Higgs state h 125 (in turn decaying into γ γ ) and a Z boson (in turn decaying into b b ¯ ), within a three-Higgs doublet model (3HDM) featuring two active and one inert doublet, known as the I ( 1 + 2 ) HDM . This theoretical structure features a spectrum with both the SM-like Higgs boson (with a 125 GeV mass) and a lighter C P -even (or scalar) Higgs state with mass around 95 GeV, h 95 , which is present in this scenario for the purpose of simultaneously explaining anomalies seen in the b b ¯ , γ γ , and τ + τ − final states in searches for additional light Higgs states at the Large Electron-Positron collider and LHC itself. It should be noted that, in the I ( 1 + 2 ) HDM , the inert sector presents loop-induced enhancements to the h 95 → γ γ width via additional scalar charged states, providing a viable mechanism to explain, in particular, the observed (and most significant) diphoton excess at 95 GeV. Taking into account both experimental and theoretical constraints, our results can not only explain the aforementioned anomalies (possibly, aside from the τ + τ − , which is the most marginal one) but also predict, as collateral signals, resonant production of the same C P -odd scalar A followed by the decays: (i) A → h 95 Z , leading to the same γ γ b b ¯ final state displaying the original 650 GeV anomaly and (ii) A → t t ¯ , leading to a well-known and studied signature. Both of these signals are potentially explorable at Run 3 of the LHC and most possibly so at the High-Luminosity LHC, while being consistent with current data at a significance level of 2.5 σ .
- 27 October 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (7)
Open AccessCould the 650 GeV excess be a pseudoscalar of a three-Higgs-doublet model?
In this study, we propose the interpretation of a 650 GeV excess observed at the Large Hadron Collider (LHC) by the CMS Collaboration in terms of the production of a C P -odd (or pseudoscalar) Higgs boson A , with mass around 650 GeV, decaying into the Standard Model (SM)-like Higgs state h 125 (in turn decaying into γ γ ) and a Z boson (in turn decaying into b b ¯ ), within a three-Higgs doublet model (3HDM) featuring two active and one inert doublet, known as the I ( 1 + 2 ) HDM . This theoretical structure features a spectrum with both the SM-like Higgs boson (with a 125 GeV mass) and a lighter C P -even (or scalar) Higgs state with mass around 95 GeV, h 95 , which is present in this scenario for the purpose of simultaneously explaining anomalies seen in the b b ¯ , γ γ , and τ + τ − final states in searches for additional light Higgs states at the Large Electron-Positron collider and LHC itself. It should be noted that, in the I ( 1 + 2 ) HDM , the inert sector presents loop-induced enhancements to the h 95 → γ γ width via additional scalar charged states, providing a viable mechanism to explain, in particular, the observed (and most significant) diphoton excess at 95 GeV. Taking into account both experimental and theoretical constraints, our results can not only explain the aforementioned anomalies (possibly, aside from the τ + τ − , which is the most marginal one) but also predict, as collateral signals, resonant production of the same C P -odd scalar A followed by the decays: (i) A → h 95 Z , leading to the same γ γ b b ¯ final state displaying the original 650 GeV anomaly and (ii) A → t t ¯ , leading to a well-known and studied signature. Both of these signals are potentially explorable at Run 3 of the LHC and most possibly so at the High-Luminosity LHC, while being consistent with current data at a significance level of 2.5 σ .
- Open Access
- 16 September 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (5)
Searching for charged Higgs bosons via at the ILCWe investigate the phenomenology of the charged Higgs boson at the International Linear Collider (ILC) within the framework of the type-X Two-Higgs Doublet Model, where a light charged Higgs boson, with a mass around 200 GeV or even smaller than top quark mass, is still being consistent with flavor physics data as well as with collider experimental data. In the theoretically and experimentally allowed parameter space, the e + e − → H ± W ∓ S (with S = H , A ) production processes can yield signatures with event rates larger than those from e + e − → H + H − and offer sensitivity to the Higgs mixing parameter sin ( β − α ) . We consider the bosonic H ± → W ± S decays, where the neutral scalar S further decays into a pair of tau leptons. We show, through a detector-level Monte Carlo analysis, that the resulting [ τ τ ] [ τ τ ] W W final state could be seen at the ILC with at least 500 GeV center-of-mass energy and 500 fb − 1 of luminosity.
- 16 September 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (5)
Open AccessSearching for charged Higgs bosons via at the ILCWe investigate the phenomenology of the charged Higgs boson at the International Linear Collider (ILC) within the framework of the type-X Two-Higgs Doublet Model, where a light charged Higgs boson, with a mass around 200 GeV or even smaller than top quark mass, is still being consistent with flavor physics data as well as with collider experimental data. In the theoretically and experimentally allowed parameter space, the e + e − → H ± W ∓ S (with S = H , A ) production processes can yield signatures with event rates larger than those from e + e − → H + H − and offer sensitivity to the Higgs mixing parameter sin ( β − α ) . We consider the bosonic H ± → W ± S decays, where the neutral scalar S further decays into a pair of tau leptons. We show, through a detector-level Monte Carlo analysis, that the resulting [ τ τ ] [ τ τ ] W W final state could be seen at the ILC with at least 500 GeV center-of-mass energy and 500 fb − 1 of luminosity.
- Open Access
- 25 June 2025
- Published by Elsevier in Physics Letters B
- Vol. 868
Exploring potential Higgs resonances at 650 GeV and 95 GeV in the 2HDM Type III
- 25 June 2025
- Published by Elsevier in Physics Letters B
- Vol. 868
Open AccessExploring potential Higgs resonances at 650 GeV and 95 GeV in the 2HDM Type III
- Open AccessSearching for charged Higgs bosons via $e^+ e^- \to H^\pm W^\mp S$ at the ILCWe investigate the phenomenology of the charged Higgs boson at the International Linear Collider (ILC) within the framework of the type-X Two-Higgs Doublet Model (2HDM), where a light charged Higgs boson, with a mass around 200 GeV or even smaller than top quark mass, is still being consistent with flavor physics data as well as with the colliders experimental data. In the theoretically and experimentally allowed parameter space, the $e^+ e^- \to H^\pm W^\mp S$ (with $S = H, A$) production processes can yield signatures with event rates larger than those from $e^+ e^- \to H^+ H^-$ and offer sensitivity to the Higgs mixing parameter $\sin(β-α)$. We consider the bosonic $H^\pm \to W^\pm S$ decays, where the neutral scalar $S$ further decays into a pair of tau leptons. We show, through a detector-level Monte Carlo analysis, that the resulting $[ττ][ττ] WW$ final state could be seen at the ILC with at least 500 GeV center-of-mass energy and 500 fb$^{-1}$ of luminosity.Open AccessSearching for charged Higgs bosons via $e^+ e^- \to H^\pm W^\mp S$ at the ILCWe investigate the phenomenology of the charged Higgs boson at the International Linear Collider (ILC) within the framework of the type-X Two-Higgs Doublet Model (2HDM), where a light charged Higgs boson, with a mass around 200 GeV or even smaller than top quark mass, is still being consistent with flavor physics data as well as with the colliders experimental data. In the theoretically and experimentally allowed parameter space, the $e^+ e^- \to H^\pm W^\mp S$ (with $S = H, A$) production processes can yield signatures with event rates larger than those from $e^+ e^- \to H^+ H^-$ and offer sensitivity to the Higgs mixing parameter $\sin(β-α)$. We consider the bosonic $H^\pm \to W^\pm S$ decays, where the neutral scalar $S$ further decays into a pair of tau leptons. We show, through a detector-level Monte Carlo analysis, that the resulting $[ττ][ττ] WW$ final state could be seen at the ILC with at least 500 GeV center-of-mass energy and 500 fb$^{-1}$ of luminosity.
- Open Access
- 22 November 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 110 (9)
Associated charged Higgs boson production within the 2HDM: versus collidersOur goal is to investigate the charged Higgs phenomenology in the framework of 2HDM at the upcoming e + e − and muon colliders. We are primarily concerned with the associated production processes with a fermion pair: ℓ + ℓ − → τ + ν τ H − and ℓ + ℓ − → t b ¯ H − , as well as with the W boson and a neutral Higgs boson: ℓ + ℓ − → W ± H ∓ S ( S = h , H , A ) and ℓ + ℓ − → W ± H ∓ Z . We first update the results for e + e − → { τ + ν τ H − , t b ¯ H − } and then discuss our findings for e + e − → { W ± H ∓ S } for various center-of-mass energies 500 GeV, 1 TeV, 1.5 TeV, and 3 TeV. In the case of muon collider, we show that the new s-channel and t-channel diagrams can increase the cross sections by virtue of their Yukawa couplings. We systematically compare our results for the muon collider with those obtained at the International Linear and Compact Linear colliders. We select benchmark points and conduct signal-background analyses, incorporating detector simulations. For a 3 TeV muon collider, our results show an exclusion region at the 2 σ level and a discovery region at the 5 σ level. Published by the American Physical Society 2024
- 22 November 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 110 (9)
Open AccessAssociated charged Higgs boson production within the 2HDM: versus collidersOur goal is to investigate the charged Higgs phenomenology in the framework of 2HDM at the upcoming e + e − and muon colliders. We are primarily concerned with the associated production processes with a fermion pair: ℓ + ℓ − → τ + ν τ H − and ℓ + ℓ − → t b ¯ H − , as well as with the W boson and a neutral Higgs boson: ℓ + ℓ − → W ± H ∓ S ( S = h , H , A ) and ℓ + ℓ − → W ± H ∓ Z . We first update the results for e + e − → { τ + ν τ H − , t b ¯ H − } and then discuss our findings for e + e − → { W ± H ∓ S } for various center-of-mass energies 500 GeV, 1 TeV, 1.5 TeV, and 3 TeV. In the case of muon collider, we show that the new s-channel and t-channel diagrams can increase the cross sections by virtue of their Yukawa couplings. We systematically compare our results for the muon collider with those obtained at the International Linear and Compact Linear colliders. We select benchmark points and conduct signal-background analyses, incorporating detector simulations. For a 3 TeV muon collider, our results show an exclusion region at the 2 σ level and a discovery region at the 5 σ level. Published by the American Physical Society 2024
- Open AccessAssociated charged Higgs production within the 2HDM: $e^-e^+$ versus $μ^-μ^+$ collidersOur goal is to investigate the charged Higgs phenomenology in the framework of 2HDM at the upcoming $e^+e^-$ and muon colliders. We are primarily concerned with the associated production processes with a fermion pair: $\ell^+ \ell^- \to \tau^+ \nu_{\tau }H^-$ and $\ell^+ \ell^- \to t \bar{b} H^-$, as well as with the W boson and a neutral Higgs boson: $\ell^+ \ell^- \to W^\pm H^\mp S$ ($S=h,\,H,\, A$) and $\ell^+ \ell^- \to W^\pm H^\mp Z$. We first update the results for $e^+e^- \to \{ \tau^+ \nu_{\tau }H^-\ , \ t \bar{b} H^-\}$ and then discuss our findings for $e^+e^- \to \{ W^\pm H^\mp S\}$ for various center of mass energies 500 GeV, 1 TeV, 1.5 TeV and 3 TeV. In the case of muon collider, we show that the new s-channel and t-channel diagrams can increase the cross sections by virtue of their Yukawa couplings. We systematically compare our results for the muon collider with those obtained at the International Linear (ILC) and Compact Linear (CLIC) colliders. We select benchmark points and conduct signal-background analyses, incorporating detector simulations. For a 3 TeV muon collider, our results show an exclusion region at the 2$\sigma$ level and a discovery region at the 5$\sigma$ level.Open AccessAssociated charged Higgs production within the 2HDM: $e^-e^+$ versus $μ^-μ^+$ collidersOur goal is to investigate the charged Higgs phenomenology in the framework of 2HDM at the upcoming $e^+e^-$ and muon colliders. We are primarily concerned with the associated production processes with a fermion pair: $\ell^+ \ell^- \to \tau^+ \nu_{\tau }H^-$ and $\ell^+ \ell^- \to t \bar{b} H^-$, as well as with the W boson and a neutral Higgs boson: $\ell^+ \ell^- \to W^\pm H^\mp S$ ($S=h,\,H,\, A$) and $\ell^+ \ell^- \to W^\pm H^\mp Z$. We first update the results for $e^+e^- \to \{ \tau^+ \nu_{\tau }H^-\ , \ t \bar{b} H^-\}$ and then discuss our findings for $e^+e^- \to \{ W^\pm H^\mp S\}$ for various center of mass energies 500 GeV, 1 TeV, 1.5 TeV and 3 TeV. In the case of muon collider, we show that the new s-channel and t-channel diagrams can increase the cross sections by virtue of their Yukawa couplings. We systematically compare our results for the muon collider with those obtained at the International Linear (ILC) and Compact Linear (CLIC) colliders. We select benchmark points and conduct signal-background analyses, incorporating detector simulations. For a 3 TeV muon collider, our results show an exclusion region at the 2$\sigma$ level and a discovery region at the 5$\sigma$ level.