# Closing the Window on Light Charged Higgs Bosons in the NMSSM

Authors: Johan Rathsman; Thomas Rössler
Journal: Advances in High Energy Physics
 (2012)

---

## Abstract

In the Next-to-Minimal SuperSymmetric Model (NMSSM) the lightest CP-odd Higgs bosons (a1) can be very light. As a consequence, in addition to the standard charged Higgs boson (h±) decays considered in the MSSM for a light charged Higgs (mh±<mt), the branching fraction for h±→a1W can be dominant. We investigate how this signal can be searched for in tt- production at the Large Hadron Collider (LHC) in the case that (ma1≳2mB) with the a1 giving rise to a single bb--jet and discuss to what extent the LHC experiments are able to discover such a scenario with an integrated luminosity ~20 fb−1. We also discuss the implications of the possible Higgs-signal observed at the LHC.


---

## Body

## 1. Introduction

With the successful start-up of the LHC and the intriguing results from the first year of data-taking at the center of mass energy of 7 TeV and an integrated luminosity close to≈5 fb−1 for each of the ATLAS and CMS experiments, the ongoing run in 2012 is set to be a milestone in particle physics. The possible signal for a Higgs boson around 125 GeV may or may not be confirmed and the search for physics Beyond the Standard Model (BSM) will continue. Contrary to the neutral Higgs boson, which if discovered may need to be analyzed in detail regarding its branching fraction into various channels in order to determine whether it is the Standard Model (SM) Higgs boson or not, the discovery of a charged Higgs boson would be an unmistakable sign of BSM physics.The charged Higgs boson arises in theories with more than one Higgs doublet. The prime example is the MSSM [1, 2] which has two Higgs doublets, leading to two CP-even Higgs bosons (h,H) and one CP-odd (A) in the case of CP-conservation and two charged states (H±) after electroweak symmetry breaking. In this case the two Higgs doublets are required by supersymmetry with one of them giving masses to the uptype fermions and one to the downtype ones. For a complete introduction to the Higgs sector in the MSSM we refer to [3].The main reason for introducing supersymmetry (for a general introduction to supersymmetric theories we refer to [4]) is to solve the so-called hierarchy problem, that is, why the scale of the electroweak interaction is so much smaller than the cut-off of the SM, normally taken to be the Planck mass where gravity becomes the dominant force and the SM breaks down. With the Higgs boson being a scalar particle, the higher-order corrections to its mass are proportional to this cut-off and with the SM only being an effective theory this cut-off dependence cannot be renormalized away. In a supersymmetric theory this problem is essentially solved by the introduction of the fermionic partners of the Higgs fields, which only get logarithmic corrections to their masses and thereby avoids fine-tuning. In turn this means that the Higgs boson masses are also protected from receiving quadratic corrections as long as supersymmetry is not broken or only softly broken.In addition to solving the hierarchy problem supersymmetry also offers a candidate for cold dark matter [5, 6] in the case that R-symmetry is preserved, which in turn is introduced to avoid terms in the Lagrangian that otherwise would mediate proton decay. In this case the lightest supersymmetric particle (LSP) is stable and generically it has the right mass and cross-section to constitute the observed dark matter. Supersymmetry also improves the unification of gauge forces at an hypothesized grand unification scale although the unification is not exact and it does depend on the details of the spectrum of SUSY-particles.As is well known, the MSSM by itself is not without problems. Leaving the question of the precise mechanism for supersymmetry breaking aside, the MSSM faces the so-calledμ-problem. This relates to the magnitude of the dimensionful μ-parameter which couples the two Higgs doublets to each other in the superpotential. In order to avoid large cancellations between this contribution to the Higgs masses and the soft supersymmetry breaking terms as well as having a phenomenologically viable supersymmetric theory (mainly having a large enough chargino mass), the magnitude of μ should be of order the electroweak or supersymmetry breaking scales. The problem is then that there is no a priori reason for this parameter to have any particular value, in principle it could be anything up to the Planck scale, so why is it similar to the electroweak or supersymmetry breaking scales?In the NMSSM theμ-problem is solved by introducing an additional Higgs singlet into the theory. After supersymmetry breaking this field gets a vacuum expectation value (vev) that effectively acts as a μ-term. The original μ-term in the superpotential can then be set to zero without spoiling the viability of the theory. For a more detailed review of the NMSSM we refer to [7, 8].The additional Higgs singlet has important consequences for the phenomenology of the Higgs sector. In the MSSM, the masses of the heavy Higgs bosons (H,A,H±) are closely related to each other as they originate from the same (second) Higgs doublet if viewed in the Higgs basis where only one (the first) of the Higgs doublets has a vev. For example, at tree-level the masses of the CP-odd and charged Higgs bosons are related by mH±2=mA2+mW2. As a consequence the decay H±→AW is typically not open.In the NMSSM the additional Higgs singlet means that there is one more CP-even and one more CP-odd field with a separate mass scale introduced into the Higgs sector. As a consequence, the by now three CP-even and two CP-odd electroweak states will mix into the respective mass eigenstates. Thus the mass-relations from the MSSM will be altered. This is particularly evident in the CP-odd sector where the lightest statea1 may now be much lighter than the charged Higgs boson—even after taking experimental constraints into account as discussed below—opening up the possibility for the h±→a1W decay to be dominant. In turn this means that the search for charged Higgs bosons, in t-quark decays for example, has to be widened also to include this decay channel.The decayh±→a1W has already been considered to different levels of detail [9–11] in the literature and there are constraints from the DELPHI experiment for ma1>12 GeV [12], as well as the CDF experiment for the case a1→τ+τ- [13]. In this paper we want to focus on the region in parameter space where the a1 mass is above the bb- threshold but still so close to it that the two b-quarks will fragment into a single bb--jet. The viability of scenarios with light a1’s has also been considered by [14–18].Our paper is organized as follows. In the next section we give some basic properties of the Higgs sector in the NMSSM that are relevant to our discussion. We then discuss the constraints on the parameter space in Section3, including the latest results from LHC. In Section 4 we illustrate how the signal h±→a1W can be searched for in tt--production taking into account the appropriate backgrounds. Section 5 contains a discussion of the implications of the possible Higgs signal from the ATLAS and CMS experiments and in Section 6 we summarize and conclude.


## 2. Basic Properties of The NMSSM

We consider theZ3-symmetric version of the NMSSM with the superpotential given by
(2.1)WNMSSM=WMSSM+λS^H^u·H^d+κS^3,
where WMSSM is the superpotential of the MSSM with μ set to zero. The soft supersymmetry breaking potential relative to the MSSM is then given by
(2.2)VsoftNMSSM=VsoftMSSM+mS2|S|2+(λAλSHu·Hd+13κAκS3+h.c.),
where the part of VsoftMSSM only depending on the Higgs fields is given by,
(2.3)Vsoft,higgsMSSM=mHu2|Hu|2+mHd2|Hd|2.In additionVsoftMSSM contains all the dependence on the other soft supersymmetry breaking parameters: the gaugino masses M1,M2,M3, the trilinear couplings au,ad,ae, the squark masses MQ,Mu,Md, and finally the slepton masses ML,Me. In the following we will assume minimal flavour violation so that the sfermion mass matrices are diagonal and the trilinear couplings are proportional to the corresponding Yukawa coupling matrices au=Auyu, and so forth.After electroweak symmetry breaking, and assuming that CP is conserved, the Higgs sector will contain three CP-even Higgs bosons (h1,h2,h3), two CP-odd (a1,a2), and one charged (h±), where the states are ordered in terms of increasing mass. In the same way as in the MSSM the minimization conditions for the Higgs potential allow one to trade the mHu,mHd parameters for the doublet vev v≈174 GeV and tanβ=vu/vd. Similarly mS can be expressed in terms of the singlet vev vs which in turn gives rise to the effective μ-parameter, μ=λvs. All in all this leaves us with 6 unknown parameters describing the Higgs sector of NMSSM at tree-level: tanβ,μ,λ,κ,Aλ,Aκ. Below we will trade the latter two parameters for the masses mh± and ma1.As already alluded to the mass-eigenstates are mixtures of the electroweak eigenstates. More specifically, writingSweak=(Re(Hu),Re(Hd),Re(S)) we have hi=SijSjweak. (In the MSSM limit this means that S12=-sinα.) Similarly for the CP-odd states we have ai=AijAjweak with Aweak=(Im(cosβHu+sinβHd),Im(S)). Here the mixing matrix is simply A=(cosθAsinθA-sinθAcosθA). Together with the ratio of the two doublet vevs (or equivalently the rotation angle needed to go to the Higgs basis where only one of the doublets have a vev) tanβ, the mixing matrices S and A specify the reduced couplings to fermions and gauge bosons as given in Table 1.Table 1
Reduced Higgs couplings in the NMSSM compared to the MSSM and the SM (when applicable). Note that the reduced couplings to fermions are identical for all three generations, even if only the third generation is displayed here. The couplings toh3 can be obtained from the h1 ones by the replacements S11→S31 and S12→S32 whereas the couplings to a2 can be obtained from the a1 ones by the replacement cosθA→sinθA.
Vertex
NMSSM
MSSM
SM
h
1
t
t
S
11
/
sin
β
cos
α
/
sin
β
1
h
1
b
b
S
12
/
cos
β
sin
α
/
cos
β
1
h
2
t
t
S
21
/
sin
β
sin
α
/
sin
β
n.a.
h
2
b
b
S
22
/
cos
β
cos
α
/
cos
β
n.a.
a
1
t
t
cot
β
cos
θ
A
cot
β
n.a.
a
1
b
b
tan
β
cos
θ
A
tan
β
n.a.
h
1
V
V
sin
β
S
11
+
cos
β
S
12
sin
(
β
-
α
)
1
h
2
V
V
sin
β
S
21
+
cos
β
S
22
cos
(
β
-
α
)
n.a.
a
1
h
1
Z
(
cos
β
S
11
-
sin
β
S
12
)
cos
θ
A
cos
(
β
-
α
)
n.a.
a
1
h
2
Z
(
cos
β
S
21
-
sin
β
S
22
)
cos
θ
A
sin
(
β
-
α
)
n.a.
h
1
h
+
W
-
cos
β
S
11
-
sin
β
S
12
cos
(
β
-
α
)
n.a.
a
1
h
+
W
-
cos
θ
A
1
n.a.


## 3. Experimental Constraints

In this section we will explore to what extent the process we are interested in is constrained by existing experimental data. Since we are interested in a lighth± (with mh±<mt) and a light a1 there are constraints both from collider experiments as well as low-energy flavour experiments. However, before going into the various constraints we will specify the scenarios that we have considered and then come back to the question of experimental constraints.


### 3.1. Specification of SUSY Scenario Considered

In the following we will consider a variant of the well motivatedmhmax-scenario in the MSSM [19], similarly to what was done in [17]. Thus we will consider a universal scale MSUSY for the sfermion masses at the supersymmetry breaking scale. In other words we assume, as already stated, that the sfermion mass matrices (MQ etc.) are diagonal, and furthermore we assume that all diagonal entries are equal to MSUSY which we keep fixed at 1 TeV. In addition we assume that the gaugino masses are related as in the constrained MSSM, where supersymmetry breaking is assumed to be mediated by gravity, namely, M1=100 GeV, M2=200 GeV, M3=800 GeV. Finally we will assume that At=Ab=Aτ but contrary to what was done in [17] we will let them vary in the range At∈[-5000,5000] GeV so that the amount of mixing between the t~L and t~R is unconstrained.For the Higgs sector we will let all six parameters vary freely. However, as was done in [17], we will trade the Aκ parameter for ma1 and Aλ parameter for mh± using an iterative procedure starting from the tree-level relations:
(3.1)mh±2=2μsin2β(Aλ+κλμ)+mW2-λ2v2ℳP2=(2μsin2β(Aλ+κλμ)λv(Aλ-2κλμ)λv(Aλ-2κλμ)λ2v2sin2β2μ(Aλ+4κλμ)-3κλAκμ),
where the latter gives the masses of the mass-eigenstates a1, a2 after diagonalisation. Thus the parameters we consider with their respective ranges are:
(3.2)tanβ∈[1,60],λ∈[0,0.7],κ∈[-0.7,0.7],μ∈[125,1000]GeV,mh±∈[80,170]GeV,ma1∈[4,150]GeV.The limits for the various parameters have been chosen as follows: Fortanβ, κ and λ we impose perturbativity up to the GUT scale which effectively means that any value out side the above regions is bound to fail. (In addition some points inside these regions also fail because of this requirement.) The lower limits on μ and mh± are dictated by experimental constraints. The upper limit on μ is not a hard one but follows from the implicit assumption that μ should be of order the electroweak scale whereas the upper limit on mh± is given by the condition that the decay t→bh+ should be open. The reason for letting μ vary freely is mainly that this decreases the correlations between the masses of the Higgs bosons as will be discussed more below. Finally the lower limit on ma1 is chosen in order to have a1→τ+τ- open whereas the upper limit follows from having mh±<170 GeV.In order to calculate the resulting models from the inputs we use the package NMSSMTools version 3.2.0 [20, 21] with default settings. Among other things this means that we impose perturbativity of the model up to the GUT scale. Finally, in all scans we generate ~1 M points (with flat priors in the parameters considered) which fulfill the theoretical constraints implemented in NMSSMTools.


### 3.2. Current Experimental Constraints

The most important constraints come from the direct searches for Higgs bosons for which we use the package HiggsBounds version 3.7.0 [22, 23]. In addition there are also constraints from direct searches for supersymmetric particles, various flavour constraints and in principle also the anomalous magnetic moment of the muon as well as the relic density of dark matter. We have not applied the latter two constraints for the following reasons. To investigate the amount of dark matter in the various models one would also need to vary the gaugino masses as was done in [24]. Since we keep these fixed we have not applied the dark matter constraints. On the same vain we have not applied the constraint from the anomalous magnetic moment of the muon, since this depends on the masses of the scalar partners of the muon and the neutrinos, apart from requiring μ to be positive.When it comes to the flavour constraints the situation is more involved in that various constraints have different level of model dependence. On the one hand there are constraints from tree-level mediated processes such asBu→τ+ντ, which only depend on the Higgs sector and on the other hand there are constraints from loop-mediated processes such as Bs→μ+μ- and b→sγ which depend on details of the supersymmetric sector of the model. In the following we will limit ourselves to applying the most severe constraints from Bu→τ+ντ, which limits the available parameter space in [mh±,tanβ] and Bs→μ+μ- which puts limits on [ma1,tanβ]. The latter constraint is especially important since we will consider light a1 which is very constrained by the data from LHCb and CMS [25, 26]. Finally we apply the direct constraints from searches for supersymmetric particles. For all constraints except the Higgs bosons we use NMSSMTools version 3.2.0.The results from the scan are displayed in Figure1 with black points being viable models from a theory point of view but excluded by the direct searches for Higgs bosons and coloured points being allowed by the same constraint. Of primary interest are the allowed regions in [ma1,mh±]. As is clear from the figure there is a distinct region of points with ma1~2mB≈10.6 GeV allowed by all constraints (indicated by green colour) for essentially any value of mh±≳90 GeV. The same is also true in the constrained mhmax scenario where μ=200 GeV and At=-6MSUSY+μcotβ are fixed as was already noted in [17]. However, at difference to the constrained mhmax scenario there are regions with larger ma1 that are allowed also for mh±≲120 GeV.Correlations between the parameters of the scan in the scenario under consideration (see text for details) with the various constraints applied as follows. All points in black are viable model points but are excluded by HiggsBounds. The coloured points are allowed by HiggsBounds and plotted in the following order: red points are excluded byBs→μ+μ-, blue points are excluded by Bu→τ+ντ, yellow points are excluded by direct searches for supersymmetric particles, and finally green points are allowed by all constraints considered.
(a)
(b)
(c)
(d)
(e)Looking at the[mh±,tanβ]-plane, one clearly sees the constraint from Bu→τ+ντ for intermediate tanβ (shown as blue points). For larger tanβ there is a cancellation between the SM and h± contributions, which makes this region allowed by Bu→τ+ντ, but instead the constraints from Bs→μ+μ- come into play (red points). It should be noted that the points excluded by Bu→τ+ντ are plotted on top of the constraints from Bs→μ+μ-. Similarly the constraints from searches for supersymmetric particles are plotted (in yellow) on top of the constraints from B-decays, but with the cut μ>125 GeV there are hardly any points excluded by this constraint. Finally we note that for tanβ≲15 there are allowed points in parameter space for essentially any value of mh±≳90 GeV.Turning to the[ma1,tanβ]-plane we see that the region ma1~2mB is allowed by all constraints up to tanβ≲10, whereas for larger tanβ there are points allowed by direct Higgs boson searches but not allowed by Bs→μ+μ- and Bu→τ+ντ. Given the uncertainties related to the indirect constraints from B-decays we conclude that there is a region in parameter space with ma1~2mB, mh±∈[90,170] GeV, and tanβ∈[1,60] that should be searched for by the ATLAS and CMS experiments. It should also be noted that values both above and below the threshold ma1=2mB are allowed by the constraints.Before turning to the signal of interest, that is,Br(h±→a1W), we also show in Figure 1 the effects of the various constraints when projected onto the [κ,λ] and [At,μ]-planes. From the first of these plots one clearly sees the constraint κ2+λ2≲0.7 which arises from requiring perturbativity up to the GUT scale. From the second we see that the constraints imply |At|≲3500 GeV, which essentially follows from the radiative corrections to the lightest CP-even Higgs becoming small or even negative for large |At| relative to the value MSUSY=1 TeV that we are using. We also see that for μ there are hardly no experimental constraints in the region considered. On the other hand, if we would extend μ to values smaller than 125 GeV then all those points would be excluded by searches for supersymmetric particles.As promised we turn now to the resulting branching ratios for the decayh±→a1W. From Figure 2 we observe the following general feature, the branching ratio can be large as soon as the channel is open (ma1<mh±-mW) except for large tanβ where the decay h±→τντ becomes dominant. Concentrating on those points that pass all the constraints considered and the region ma1~2mB we also see from the lower right plot that Br(h±→a1W)≲0.9 as long as tanβ≲10 and mh±≳100 GeV. Thus we can conclude that not only is the parameter space region ma1~2mB, mh±∈[100,170] GeV, and tanβ∈[1,10] allowed in this region, but the decay h±→a1W is also dominant. In the next section we will exemplify how the so far unexplored region of parameter space with ma1∈[10,12] GeV can be probed by searching for h±→a1W with a1→bb- in tt--production at the LHC.The branching ratio forh±→a1W when scanning over the parameters in the scenario under consideration (see text for details) with the various constraints applied. The colour coding is the same as in Figure 1.
(a)
(b)
(c)
(d)Before ending this section we also show in Figure3 the branching rations for the decay chains of interest, that is, t→bh+, t→bh+→ba1W, and t→bh+→ba1W→bbb-W as a function of tanβ when restricting to parameter space points with mh±<160 GeV. As can be seen from the figure, points with Br(t→bh+) as large as 0.2 are still allowed when the decay h±→a1W is included. This should be compared with the experimental constraints from ATLAS and CMS which so far have assumed that Br(h±→τντ)=1 giving a limit Br(t→bh+)≲ few percent [27, 28].The branching ratios for the decay chainst→bh+, t→bh+→ba1W, and t→bh+→ba1W→bbb-W, respectively, when scanning over the parameters in the scenario under consideration (see text for details) with the various constraints applied and only considering points with mh±<160 GeV. The colour coding is the same as in Figure 1.
(a)
(b)
(c)


## 3.1. Specification of SUSY Scenario Considered

In the following we will consider a variant of the well motivatedmhmax-scenario in the MSSM [19], similarly to what was done in [17]. Thus we will consider a universal scale MSUSY for the sfermion masses at the supersymmetry breaking scale. In other words we assume, as already stated, that the sfermion mass matrices (MQ etc.) are diagonal, and furthermore we assume that all diagonal entries are equal to MSUSY which we keep fixed at 1 TeV. In addition we assume that the gaugino masses are related as in the constrained MSSM, where supersymmetry breaking is assumed to be mediated by gravity, namely, M1=100 GeV, M2=200 GeV, M3=800 GeV. Finally we will assume that At=Ab=Aτ but contrary to what was done in [17] we will let them vary in the range At∈[-5000,5000] GeV so that the amount of mixing between the t~L and t~R is unconstrained.For the Higgs sector we will let all six parameters vary freely. However, as was done in [17], we will trade the Aκ parameter for ma1 and Aλ parameter for mh± using an iterative procedure starting from the tree-level relations:
(3.1)mh±2=2μsin2β(Aλ+κλμ)+mW2-λ2v2ℳP2=(2μsin2β(Aλ+κλμ)λv(Aλ-2κλμ)λv(Aλ-2κλμ)λ2v2sin2β2μ(Aλ+4κλμ)-3κλAκμ),
where the latter gives the masses of the mass-eigenstates a1, a2 after diagonalisation. Thus the parameters we consider with their respective ranges are:
(3.2)tanβ∈[1,60],λ∈[0,0.7],κ∈[-0.7,0.7],μ∈[125,1000]GeV,mh±∈[80,170]GeV,ma1∈[4,150]GeV.The limits for the various parameters have been chosen as follows: Fortanβ, κ and λ we impose perturbativity up to the GUT scale which effectively means that any value out side the above regions is bound to fail. (In addition some points inside these regions also fail because of this requirement.) The lower limits on μ and mh± are dictated by experimental constraints. The upper limit on μ is not a hard one but follows from the implicit assumption that μ should be of order the electroweak scale whereas the upper limit on mh± is given by the condition that the decay t→bh+ should be open. The reason for letting μ vary freely is mainly that this decreases the correlations between the masses of the Higgs bosons as will be discussed more below. Finally the lower limit on ma1 is chosen in order to have a1→τ+τ- open whereas the upper limit follows from having mh±<170 GeV.In order to calculate the resulting models from the inputs we use the package NMSSMTools version 3.2.0 [20, 21] with default settings. Among other things this means that we impose perturbativity of the model up to the GUT scale. Finally, in all scans we generate ~1 M points (with flat priors in the parameters considered) which fulfill the theoretical constraints implemented in NMSSMTools.


## 3.2. Current Experimental Constraints

The most important constraints come from the direct searches for Higgs bosons for which we use the package HiggsBounds version 3.7.0 [22, 23]. In addition there are also constraints from direct searches for supersymmetric particles, various flavour constraints and in principle also the anomalous magnetic moment of the muon as well as the relic density of dark matter. We have not applied the latter two constraints for the following reasons. To investigate the amount of dark matter in the various models one would also need to vary the gaugino masses as was done in [24]. Since we keep these fixed we have not applied the dark matter constraints. On the same vain we have not applied the constraint from the anomalous magnetic moment of the muon, since this depends on the masses of the scalar partners of the muon and the neutrinos, apart from requiring μ to be positive.When it comes to the flavour constraints the situation is more involved in that various constraints have different level of model dependence. On the one hand there are constraints from tree-level mediated processes such asBu→τ+ντ, which only depend on the Higgs sector and on the other hand there are constraints from loop-mediated processes such as Bs→μ+μ- and b→sγ which depend on details of the supersymmetric sector of the model. In the following we will limit ourselves to applying the most severe constraints from Bu→τ+ντ, which limits the available parameter space in [mh±,tanβ] and Bs→μ+μ- which puts limits on [ma1,tanβ]. The latter constraint is especially important since we will consider light a1 which is very constrained by the data from LHCb and CMS [25, 26]. Finally we apply the direct constraints from searches for supersymmetric particles. For all constraints except the Higgs bosons we use NMSSMTools version 3.2.0.The results from the scan are displayed in Figure1 with black points being viable models from a theory point of view but excluded by the direct searches for Higgs bosons and coloured points being allowed by the same constraint. Of primary interest are the allowed regions in [ma1,mh±]. As is clear from the figure there is a distinct region of points with ma1~2mB≈10.6 GeV allowed by all constraints (indicated by green colour) for essentially any value of mh±≳90 GeV. The same is also true in the constrained mhmax scenario where μ=200 GeV and At=-6MSUSY+μcotβ are fixed as was already noted in [17]. However, at difference to the constrained mhmax scenario there are regions with larger ma1 that are allowed also for mh±≲120 GeV.Correlations between the parameters of the scan in the scenario under consideration (see text for details) with the various constraints applied as follows. All points in black are viable model points but are excluded by HiggsBounds. The coloured points are allowed by HiggsBounds and plotted in the following order: red points are excluded byBs→μ+μ-, blue points are excluded by Bu→τ+ντ, yellow points are excluded by direct searches for supersymmetric particles, and finally green points are allowed by all constraints considered.
(a)
(b)
(c)
(d)
(e)Looking at the[mh±,tanβ]-plane, one clearly sees the constraint from Bu→τ+ντ for intermediate tanβ (shown as blue points). For larger tanβ there is a cancellation between the SM and h± contributions, which makes this region allowed by Bu→τ+ντ, but instead the constraints from Bs→μ+μ- come into play (red points). It should be noted that the points excluded by Bu→τ+ντ are plotted on top of the constraints from Bs→μ+μ-. Similarly the constraints from searches for supersymmetric particles are plotted (in yellow) on top of the constraints from B-decays, but with the cut μ>125 GeV there are hardly any points excluded by this constraint. Finally we note that for tanβ≲15 there are allowed points in parameter space for essentially any value of mh±≳90 GeV.Turning to the[ma1,tanβ]-plane we see that the region ma1~2mB is allowed by all constraints up to tanβ≲10, whereas for larger tanβ there are points allowed by direct Higgs boson searches but not allowed by Bs→μ+μ- and Bu→τ+ντ. Given the uncertainties related to the indirect constraints from B-decays we conclude that there is a region in parameter space with ma1~2mB, mh±∈[90,170] GeV, and tanβ∈[1,60] that should be searched for by the ATLAS and CMS experiments. It should also be noted that values both above and below the threshold ma1=2mB are allowed by the constraints.Before turning to the signal of interest, that is,Br(h±→a1W), we also show in Figure 1 the effects of the various constraints when projected onto the [κ,λ] and [At,μ]-planes. From the first of these plots one clearly sees the constraint κ2+λ2≲0.7 which arises from requiring perturbativity up to the GUT scale. From the second we see that the constraints imply |At|≲3500 GeV, which essentially follows from the radiative corrections to the lightest CP-even Higgs becoming small or even negative for large |At| relative to the value MSUSY=1 TeV that we are using. We also see that for μ there are hardly no experimental constraints in the region considered. On the other hand, if we would extend μ to values smaller than 125 GeV then all those points would be excluded by searches for supersymmetric particles.As promised we turn now to the resulting branching ratios for the decayh±→a1W. From Figure 2 we observe the following general feature, the branching ratio can be large as soon as the channel is open (ma1<mh±-mW) except for large tanβ where the decay h±→τντ becomes dominant. Concentrating on those points that pass all the constraints considered and the region ma1~2mB we also see from the lower right plot that Br(h±→a1W)≲0.9 as long as tanβ≲10 and mh±≳100 GeV. Thus we can conclude that not only is the parameter space region ma1~2mB, mh±∈[100,170] GeV, and tanβ∈[1,10] allowed in this region, but the decay h±→a1W is also dominant. In the next section we will exemplify how the so far unexplored region of parameter space with ma1∈[10,12] GeV can be probed by searching for h±→a1W with a1→bb- in tt--production at the LHC.The branching ratio forh±→a1W when scanning over the parameters in the scenario under consideration (see text for details) with the various constraints applied. The colour coding is the same as in Figure 1.
(a)
(b)
(c)
(d)Before ending this section we also show in Figure3 the branching rations for the decay chains of interest, that is, t→bh+, t→bh+→ba1W, and t→bh+→ba1W→bbb-W as a function of tanβ when restricting to parameter space points with mh±<160 GeV. As can be seen from the figure, points with Br(t→bh+) as large as 0.2 are still allowed when the decay h±→a1W is included. This should be compared with the experimental constraints from ATLAS and CMS which so far have assumed that Br(h±→τντ)=1 giving a limit Br(t→bh+)≲ few percent [27, 28].The branching ratios for the decay chainst→bh+, t→bh+→ba1W, and t→bh+→ba1W→bbb-W, respectively, when scanning over the parameters in the scenario under consideration (see text for details) with the various constraints applied and only considering points with mh±<160 GeV. The colour coding is the same as in Figure 1.
(a)
(b)
(c)


## 4.  Search Strategy

In the following section we will perform a signal-to-background analysis for three different charged Higgs masses:mh+=100,130 and 150 GeV, respectively. For definiteness we have used tanβ=50 when simulating the signal but the end results will not depend on this value. The mass of the a1 is set to 11 GeV throughout as an example of a small a1 mass that is just above bb--threshold. The three charged Higgs masses are chosen to illustrate different kinematic properties: at mh±=150 GeV the b-jet from the t→bh+ decay will be rather soft whereas the a1 from the h+ will be harder. For mh±=100 GeV the situation will be the opposite, and in the intermediate case mh±=130 GeV both jets can be relatively hard and in addition the available phase-space will be largest.We aim to reconstruct the signal process where one of thet-quarks decays leptonically via h±→a1W with W→ℓνℓ and the other hadronically via W→jj as illustrated in Figure 4. All cross-sections have been corrected for these enforced W decays as well as a factor of 2 for taking account of the process also with interchanged roles between the t and the t-. Because the a1 is supposed to decay close to threshold to bb-, we aim for a reconstruction where the two b’s from the a1 are clustered together to give a single bb--jet.Figure 4
Illustrations oftt- production in proton-proton collisions with the subsequent decay chain considered in this paper.As backgrounds to the process we consider the irreduciblett-bb- as well as, because of its higher magnitude, tt- with one jet being accidentally b-tagged (weighted with a mistagging probability, assumed to be 0.01 [29, 30]). In order to include also the single top contributions to the background we have simulated the processes t-bW+bb- and t-bW+, respectively, but in the following we will denote them as tt-bb- and tt- for simplicity. For other reducible backgrounds, such as W+b jets, we assume that similar procedures can be applied as in the tt- cross-section determination. For example, requiring two b-tagged jets reduces the W+b jets background to tt- production to about 10% [31]. Several cuts are applied to strengthen the signal and to suppress the background as will be discussed in the following.A center-of-mass energy of 8 TeV at the Large Hadron Collider is assumed throughout the whole analysis. For the generation of the hard matrix elements, we use MadGraph 5 [32] with a fixed renormalization and factorization scale (set as default to the Z mass) and the “CTEQ6L1” parton distribution functions. To supply MadGraph with the proper parameters of the signal we have for simplicity used a simple two Higgs Doublet Model with the masses given above as implemented in the Two Higgs Doublet Model Calculator 2HDMC [33]. The only important difference to the NMSSM then arises from the h±→a1W decay giving and extra factor cos2θA (cf. Table 1). All other steps to generate complete events, such as radiation, underlying events, and hadronization, are carried out using Pythia 8. We start from bare samples of 100000 events for the different signals as well as the tt-bb- background whereas for the tt- background we have 50 times higher statistics to start with. There is no detector simulation included in this exploratory analysis but we have simulated b-tagging in a simplified way as detailed below.For all processes we use the leading order cross-sections obtained from MadGraph. On the one hand this means that there is an overall scale factor which is more or less the same for both signal and background. On the other hand the rates will be lower than what would have been resulted if higher order cross-sections had been used. All in all this means that the signal over background rates we find will be underestimated in this respect. For example we get a LO cross-section forpp→tt- process in the SM of 138 pb to be compared with the NNLL resummed result of 232 pb [34].


### 4.1. Reconstruction of the LeptonicW

To reconstruct the leptonicW, we first need to identify the charged lepton (e or μ) associated to the hard process. The transverse momentum (p⊥) and pseudorapidity (η)-distributions are shown in Figure 5. After applying cuts on the lepton kinematics (p⊥>20 GeV, |η|<2.5), we require the summed p⊥ of the surrounding particles in an (η,ϕ)-cone of size ΔR=Δη2+Δϕ2=0.3 around the lepton to be less than 10 GeV to call it isolated. (The numerical values for the p⊥ cut and the cone size have been optimized by observation of the changes in efficiency and purity when varying the cuts). On the whole event then, we require to have precisely one isolated lepton in the final state.Transverse momentum and pseudorapidity of the charged lepton from theW decay and missing transverse energy (MET) before applying cuts.
(a)
(b)
(c)The next step in reconstructing the leptonicW is to identify the missing energy (MET)/missing p→⊥ of the event with the transverse momentum of the neutrino. Assuming the W to be on mass-shell and using a mass-less neutrino then leaves two possible solutions for the longitudinal momentum of the neutrino pνz=A±D with
(4.1)A=mW2plz+2plz(p→ν⊥p→l⊥)2m⊥l2,D=El2m⊥l4(mW44-(p→ν⊥)2(p→l⊥)2+(p→ν⊥p→l⊥)2+(p→ν⊥p→l⊥)mW2).In order to have a more accurate reconstruction ofpνz we have applied a cut E⊥>30 GeV. In addition, different selection criteria for the choice of the sign have been examined, for example, the invariant mass of the reconstructed h±. Among these, the most viable one turned out to be a simple selection of the smaller |pνz|, which is correct in roughly three quarters of the signal events.


### 4.2. Jet Reconstruction and Tagging

For the reconstruction of the hadronic part of the event, we consider different jet clustering schemes (anti-kT [35], Cambridge/Aachen [36], kT [37, 38]) as well as cone-sizes. All particles, except neutrinos and the isolated lepton, in the rapidity range |η|<4.9 are fed into FastJet [39, 40]. The resulting clustered jets are required to have p⊥>20 GeV. Afterwards, a simplified b-tagging is simulated for all jets in the region |η|<2.5 by comparing the (η,ϕ) of the jets to the b-quarks of the hard process. All jets within ΔR=0.4 are then classified as b-jets.The jet-algorithms require to specify the distance measureR used when calculating the jet measure. Since we want to cluster the bb- pair from the a1 into one jet, it is crucial that the distance between the two subjets is not too large. Figure 6 shows the distance measure between the two b-quarks on parton level for different a1 masses. As can be seen from the figure, for reasonable clustering cone sizes (around 0.4 to 0.6) and with ma1 in the region of interest for this analysis, the two b-quarks will most likely be clustered together as a single jet.Figure 6
The distance between the twob-quarks coming from the a1-decay on parton level for the signal with mh±=130 GeV.As an ideal reconstruction will now give rise to threeb-jets, the correct reconstruction of the h± will be enhanced by the identification of the “wrong” b-jet which comes from the t- together with the hadronically decaying W-. The strategy to achieve this here is to first find the pair of untagged jets that is closest to the W mass and then combine this pair with the b-jet that gives a mass closest to the t mass. This b-jet will then be excluded in the h± reconstruction, reducing the number of b-jets which have to be considered (in an event with a so far correct clustering in the desired way) to two. In addition to this, we put cuts on the quality of the reconstruction by requiring the reconstructed masses shown in Figure 7 to be in the regions mW±20 GeV or mt±30 GeV, respectively. We have checked that the reconstruction of the W and t masses is quite independent of the choice of jet clustering scheme and cone-size as is also shown in Figure 7. Only the 0.3 cone size gives a slightly inferior top reconstruction, but then the a1 will also not give rise to a single b-jet. If not mentioned differently, we thus use the anti-kT algorithm with R=0.5 in the following.Mass reconstruction of the hadronically decayingW (a) and t-quark (b) for the signal with mh±=130 GeV when R=0.5 (a) and using the anti-kT algorithm (b).
(a)
(b)Using the jet-reconstruction outlined above we obtain the jet andb-jet multiplicity, respectively, shown in Figure 8. As is clear from the figure, the number of jets peaks around the expected five and the number of b-jets is typically smaller than ideal which is due to the limited b-tagging region enforced. Also note the large decrease in b-jet multiplicity for the tt- sample from 2 to 3. Here the 3b-jet sample arises from gluon splitting into bb- and we do not take it into account below since that would amount to double counting with the tt-bb- background.Jet (a) andb-jet (b) multiplicity for the signal and backgrounds when using the anti-kT algorithm with R=0.5.
(a)
(b)Given the large background fromtt- for the two b-jet sample we resort to requiring at least three b-jets. Assuming that one of the b-jets has been identified as coming from the t-quark this leaves two b-jets that may originate from the h±. Due to the strong dependence of the hardness of the b-jet from the t→h+b decay on the h± mass, we consider both of these remaining solutions in general as possible candidates for the h± reconstruction. The resulting distributions when combining with the leptonically decaying W are shown in Figure 9. In the h± reconstruction, we thus require 3b-jets for the signals and the tt-bb- background, while we require 2b-jets in the tt- sample. For the latter we then assume that any of the non b-jets inside the b-tagger region (|η|<2.5) can be mistagged with a probability of 0.01 per jet [29, 30].Reconstructed charged Higgs masses, using either the higher-p⊥ (a) or the lower-p⊥ (b) b-jet combined with the leptonically decaying W.
(a)
(b)


### 4.3. Signal Significance and Reach

To estimate the signal reach as well as its significance, we choose a common window of 90 to 160 GeV to integrate the signal as well as the backgrounds. We correct forb-tagging efficiencies which we assume to be 0.6 per b-jet. The resulting cross-sections are shown in Table 2.Table 2
Integrated signal cross-sections (cf. Figure9) for the three different charged Higgs masses as well as the two backgrounds. The integration region is 90 to 160 GeV.
σ
peak [fb]
signal/background
100 GeV
130 GeV
150 GeV
t
t
-
b
b
-
t
t
-
ΣBG
Highp⊥ channel
0.78
7.9
3.7
3.0
0.9
3.9
Lowp⊥ channel
0.97
9.2
4.7
4.8
1.2
6.0S
/
B ratios for the three different mass cases are then calculated for an integrated luminosity of 20 fb−1) and summarized in Table 3. The table also shows the branching ratios at the simulated parameter points, as well as the extrapolations to branching ratios necessary for a 5σ discovery. From the table it is clear that for mh±~130 GeV, the discovery reach is maximal. For smaller or larger mh±, the limitations in phase space will reduce the branching ratio for h±→a1W and t→h±b respectively. In the former case this means that the standard decay channel h±→τντ will also be significant.Table 3
TheS/B ratios obtained for an integrated luminosity of 20 fb−1 for the different signals considered together with the branching ratios for the total decay chain t→bh+→ba1W→bbb-W of the respective parameter points and a linear extrapolation to the critical branching ratio necessary for a 5σ discovery.
m
h
±
100 GeV
130 GeV
150 GeV
b
-jet selection
highp⊥
lowp⊥
highp⊥
lowp⊥
highp⊥
lowp⊥
S
/
B
1.77
1.76
18.0
16.7
8.5
8.5
BR (t→bh+→ba1W→bbb-W)
0.0051
0.022
0.015
BRcrit
0.014
0.014
0.0060
0.0065
0.0085
0.0085Before ending this section we note that similarly to the standard decay modes of the charged Higgs boson it should be possible to use the spin-correlations between the decay products of the two top quarks as a way of enhancing the signal [41–43].


## 4.1. Reconstruction of the LeptonicW

To reconstruct the leptonicW, we first need to identify the charged lepton (e or μ) associated to the hard process. The transverse momentum (p⊥) and pseudorapidity (η)-distributions are shown in Figure 5. After applying cuts on the lepton kinematics (p⊥>20 GeV, |η|<2.5), we require the summed p⊥ of the surrounding particles in an (η,ϕ)-cone of size ΔR=Δη2+Δϕ2=0.3 around the lepton to be less than 10 GeV to call it isolated. (The numerical values for the p⊥ cut and the cone size have been optimized by observation of the changes in efficiency and purity when varying the cuts). On the whole event then, we require to have precisely one isolated lepton in the final state.Transverse momentum and pseudorapidity of the charged lepton from theW decay and missing transverse energy (MET) before applying cuts.
(a)
(b)
(c)The next step in reconstructing the leptonicW is to identify the missing energy (MET)/missing p→⊥ of the event with the transverse momentum of the neutrino. Assuming the W to be on mass-shell and using a mass-less neutrino then leaves two possible solutions for the longitudinal momentum of the neutrino pνz=A±D with
(4.1)A=mW2plz+2plz(p→ν⊥p→l⊥)2m⊥l2,D=El2m⊥l4(mW44-(p→ν⊥)2(p→l⊥)2+(p→ν⊥p→l⊥)2+(p→ν⊥p→l⊥)mW2).In order to have a more accurate reconstruction ofpνz we have applied a cut E⊥>30 GeV. In addition, different selection criteria for the choice of the sign have been examined, for example, the invariant mass of the reconstructed h±. Among these, the most viable one turned out to be a simple selection of the smaller |pνz|, which is correct in roughly three quarters of the signal events.


## 4.2. Jet Reconstruction and Tagging

For the reconstruction of the hadronic part of the event, we consider different jet clustering schemes (anti-kT [35], Cambridge/Aachen [36], kT [37, 38]) as well as cone-sizes. All particles, except neutrinos and the isolated lepton, in the rapidity range |η|<4.9 are fed into FastJet [39, 40]. The resulting clustered jets are required to have p⊥>20 GeV. Afterwards, a simplified b-tagging is simulated for all jets in the region |η|<2.5 by comparing the (η,ϕ) of the jets to the b-quarks of the hard process. All jets within ΔR=0.4 are then classified as b-jets.The jet-algorithms require to specify the distance measureR used when calculating the jet measure. Since we want to cluster the bb- pair from the a1 into one jet, it is crucial that the distance between the two subjets is not too large. Figure 6 shows the distance measure between the two b-quarks on parton level for different a1 masses. As can be seen from the figure, for reasonable clustering cone sizes (around 0.4 to 0.6) and with ma1 in the region of interest for this analysis, the two b-quarks will most likely be clustered together as a single jet.Figure 6
The distance between the twob-quarks coming from the a1-decay on parton level for the signal with mh±=130 GeV.As an ideal reconstruction will now give rise to threeb-jets, the correct reconstruction of the h± will be enhanced by the identification of the “wrong” b-jet which comes from the t- together with the hadronically decaying W-. The strategy to achieve this here is to first find the pair of untagged jets that is closest to the W mass and then combine this pair with the b-jet that gives a mass closest to the t mass. This b-jet will then be excluded in the h± reconstruction, reducing the number of b-jets which have to be considered (in an event with a so far correct clustering in the desired way) to two. In addition to this, we put cuts on the quality of the reconstruction by requiring the reconstructed masses shown in Figure 7 to be in the regions mW±20 GeV or mt±30 GeV, respectively. We have checked that the reconstruction of the W and t masses is quite independent of the choice of jet clustering scheme and cone-size as is also shown in Figure 7. Only the 0.3 cone size gives a slightly inferior top reconstruction, but then the a1 will also not give rise to a single b-jet. If not mentioned differently, we thus use the anti-kT algorithm with R=0.5 in the following.Mass reconstruction of the hadronically decayingW (a) and t-quark (b) for the signal with mh±=130 GeV when R=0.5 (a) and using the anti-kT algorithm (b).
(a)
(b)Using the jet-reconstruction outlined above we obtain the jet andb-jet multiplicity, respectively, shown in Figure 8. As is clear from the figure, the number of jets peaks around the expected five and the number of b-jets is typically smaller than ideal which is due to the limited b-tagging region enforced. Also note the large decrease in b-jet multiplicity for the tt- sample from 2 to 3. Here the 3b-jet sample arises from gluon splitting into bb- and we do not take it into account below since that would amount to double counting with the tt-bb- background.Jet (a) andb-jet (b) multiplicity for the signal and backgrounds when using the anti-kT algorithm with R=0.5.
(a)
(b)Given the large background fromtt- for the two b-jet sample we resort to requiring at least three b-jets. Assuming that one of the b-jets has been identified as coming from the t-quark this leaves two b-jets that may originate from the h±. Due to the strong dependence of the hardness of the b-jet from the t→h+b decay on the h± mass, we consider both of these remaining solutions in general as possible candidates for the h± reconstruction. The resulting distributions when combining with the leptonically decaying W are shown in Figure 9. In the h± reconstruction, we thus require 3b-jets for the signals and the tt-bb- background, while we require 2b-jets in the tt- sample. For the latter we then assume that any of the non b-jets inside the b-tagger region (|η|<2.5) can be mistagged with a probability of 0.01 per jet [29, 30].Reconstructed charged Higgs masses, using either the higher-p⊥ (a) or the lower-p⊥ (b) b-jet combined with the leptonically decaying W.
(a)
(b)


## 4.3. Signal Significance and Reach

To estimate the signal reach as well as its significance, we choose a common window of 90 to 160 GeV to integrate the signal as well as the backgrounds. We correct forb-tagging efficiencies which we assume to be 0.6 per b-jet. The resulting cross-sections are shown in Table 2.Table 2
Integrated signal cross-sections (cf. Figure9) for the three different charged Higgs masses as well as the two backgrounds. The integration region is 90 to 160 GeV.
σ
peak [fb]
signal/background
100 GeV
130 GeV
150 GeV
t
t
-
b
b
-
t
t
-
ΣBG
Highp⊥ channel
0.78
7.9
3.7
3.0
0.9
3.9
Lowp⊥ channel
0.97
9.2
4.7
4.8
1.2
6.0S
/
B ratios for the three different mass cases are then calculated for an integrated luminosity of 20 fb−1) and summarized in Table 3. The table also shows the branching ratios at the simulated parameter points, as well as the extrapolations to branching ratios necessary for a 5σ discovery. From the table it is clear that for mh±~130 GeV, the discovery reach is maximal. For smaller or larger mh±, the limitations in phase space will reduce the branching ratio for h±→a1W and t→h±b respectively. In the former case this means that the standard decay channel h±→τντ will also be significant.Table 3
TheS/B ratios obtained for an integrated luminosity of 20 fb−1 for the different signals considered together with the branching ratios for the total decay chain t→bh+→ba1W→bbb-W of the respective parameter points and a linear extrapolation to the critical branching ratio necessary for a 5σ discovery.
m
h
±
100 GeV
130 GeV
150 GeV
b
-jet selection
highp⊥
lowp⊥
highp⊥
lowp⊥
highp⊥
lowp⊥
S
/
B
1.77
1.76
18.0
16.7
8.5
8.5
BR (t→bh+→ba1W→bbb-W)
0.0051
0.022
0.015
BRcrit
0.014
0.014
0.0060
0.0065
0.0085
0.0085Before ending this section we note that similarly to the standard decay modes of the charged Higgs boson it should be possible to use the spin-correlations between the decay products of the two top quarks as a way of enhancing the signal [41–43].


## 5. Compatibility with Possible Higgs Signal

Recently the ATLAS and CMS experiments announced the combined results of the SM Higgs (ϕ) searches using ≈5 fb−1 of data from 2011 [44, 45]. In short they can be summarized as follows: the CMS experiment has ruled out the region mϕ∈[129,600] GeV at the 95% confidence level as expected whereas in the region [118,129] GeV they have not been able to make any exclusion at all even though they expected to be able to do so at the 95% confidence level, the ATLAS experiment has similarly ruled out the regions mϕ∈[112.9,115.5] GeV, mϕ∈[131,238] GeV, and mϕ∈[251,466] GeV at the 95% confidence level whereas they had expected to rule out the region [125,519] GeV.Instead, both experiments have found an excess of events in the regionsmϕ~126 GeV and mϕ~124 GeV for ATLAS and CMS, respectively, which when corrected for the so-called look-elsewhere-effect and combining the different channels has a statistical significance of about 2σ for each of the two experiments. Although not statistically significant these results have stirred a lot of excitement in the high-energy physics community [24, 46–54]. It will most likely not be possible to draw any final conclusions about whether this is a true signal or not before the end of the run in 2012 which will give another ~20 fb−1 of data.One of the most important properties of the possible signal at the LHC is that theϕ→γγ channel is similar to what is expected from the SM. Therefore we start by considering the would-be signal from the h1 as well as the h2 compared to what is expected in the SM for a Higgs boson (ϕ) with the same mass. Assuming that gluon-gluon fusion dominates the production, which we have verified is always the case in the scenarios we consider, this ratio is given by
(5.1)Rggγγhi=σ(gg→hi)NMSSMσ(gg→ϕ)SMBr(hi→γγ)NMSSMBr(ϕ→γγ)SM≃Γ(hi→gg)NMSSMΓ(ϕ→gg)SMBr(hi→γγ)NMSSMBr(ϕ→γγ)SM,
where in the second equality, following for example [47, 50], we have made the implicit assumption that the difference in radiative corrections for the production and decay processes are canceled in the ratio.The results obtained for this quantity when using the same scan as in Section3 are displayed in Figure 10. From the figure it is clear that if the possible signal seen at the LHC is the h1 then we would have to have ma1≳60 GeV. Even so Rggγγh1 does not become larger than ~0.5 for mh1≳120 GeV in our scan and there appears to be an upper bound mh1≲122 GeV. However, the difference in mass compared to the possible observation is so small that it should be considered to be within the theoretical uncertainty. (For example, including the full one-loop corrections and the two-loop ones from the t and b Yukawa couplings could push the h1 mass bit higher in compliance with the possible signal. Similarly increasing MSUSY would also increase mh1.) In any case it is clear that it is hardly possible to have a light a1 if it is the h1 that has been seen at the LHC.The signal forh1→γγ (left) and h2→γγ (right) relative to the standard model when scanning over the parameters in the scenario under consideration (see text for details) with the various constraints applied. The colour coding is the same as in Figure 1.
(a)
(b)
(c)
(d)
(e)
(f)Next we turn to the possibility that it is theh2 that has been seen by the LHC experiments. The results of the scan are also shown in Figure 10. As can be seen from the figure the results are more promising in this case. There are points with Rggγγh2≳0.1 also for small ma1 and intermediate mh± that have mh2 in the region indicated by the LHC experiments.In order to explore this possibility more closely we show in Figure11 the results when fixing the a1 mass to 11 GeV as was done in Section 4 for the three cases mh±=100,130,160 GeV. As can be seen from the figure, for the intermediate charged Higgs mass it is possible to reach Rggγγh2~0.15. However it should be noted that in this case, as also shown in the figure, the branching fraction for h±→a1W is quite small meaning that the standard decay channel h±→τντ can be used even though with a slightly reduced branching fraction.The signal forh2→γγ relative to the standard model (left) and the branching fraction for h±→a1W (right) in a scan with ma1=11 GeV and mh± fixed at 100 (top), 130 (middle), and 160 (bottom) GeV, respectively. The colour coding for the various constraints is the same as in Figure 1.
(a)
(b)
(c)
(d)
(e)
(f)The last logical possibility would be that it is theh3 that has been observed by the CERN experiments. However, in the scenarios we consider it turns out that Rggγγh3 is always small.The trends seen in Figure11 can be understood on more general grounds from the difficulties of having a light a1 with mass ma1<mhi/2 and at the same time have a large Rggγγhi. The first problem is that unless the triple Higgs coupling ghia1a1 is small the decay hi→a1a1 will become dominant. Looking at the structures of ghia1a1, which for example can be found in [8], this means that both λ and κ typically have to be small. Secondly the decay hi→a1Z will also dominate if ma1<mhi-mZ unless the corresponding reduced coupling given in Table 1 is small. In other words we need (cosβSi1-sinβSi2)cosθA to be small. There are essentially three ways to achieve this. If cosθA is small then a1 will be mainly singlet like and decouple. However, then h±→a1W will also be small. The second possibility is that Si1 and Si2 are small, but then the would-be-signal would be mainly singlet-like and not produced in the first place, so Si3 has to be small. Finally then, the combination (cosβSi1-sinβSi2) could be small but then the complementary combination (sinβSi1+cosβSi2) would have to be large giving an increased coupling hi→VV. All in all this means that is difficult to have a light a1 that it is not decoupled and still have a large Rggγγhi, although it cannot be completely ruled out.


## 6. Summary and Conclusions

We have considered a well motivated class of supersymmetric extensions of the standard model with nonminimal Higgs sector, namely, the CP-conserving NMSSM. In these types of models the additional Higgs singlet can modify the phenomenology of the Higgs sector in many different ways. In this paper we have specifically addressed the possibility of having a light CP-odd Higgs boson close to, but still above thebb- threshold. This in turn means that the light charged Higgs boson, with mass mh±<mt, can decay into a1W in addition to the standard decays h±→τντ, thus invalidating the interpretations made of charged Higgs bosons searches assuming that Br(h±→τντ)=1.When investigating the viability of these types of scenarios we have found that the experimental constraints from direct searches are quite weak even when taking the latest constraints from LHC into account. The constraints from indirect searches inB-decays are more constraining but also more model dependent. Even when including the results from the most important ones, namely Bu→τντ and Bs→μ+μ-, we still find a region of parameter space that is allowed with tanβ∈[1,10]. This is precisely the same region as the one where the decay h±→a1W can be dominant.The phenomenology of these types of scenarios is special in that, due to its low mass, thea1 will decay into a single bb- jet. Even so we have shown that it is possible to reconstruct the decay h±→a1W using standard jet finding algorithms when the W decays leptonically. This requires to use the missing transverse momentum to calculate the four momentum of the W, which can then be combined with the a1-jet to give a mass peak at mh±. The other t quark is assumed to decay hadronically according to t→bjj giving an additional handle to identify the events of interest. The most important background is thus the irreducible one from tt-bb- production but we have also taken into account the tt- background by considering the possibility of mistagging ordinary jets as b-jets.Based on our study we find that with an integrated luminosity of 20 fb−1 it should be possible to discover a charged Higgs bosons in these types of scenarios as long as the combined branching fraction for the decay chain t→bh+→ba1W→bbb-W is larger than ≈0.01.Finally we have also investigated the phenomenological consequences of the possible Higgs signal seen at the LHC on the types of scenarios we consider. We find that it is difficult to have a lighta1 that is not decoupled and at the same time have a combined production and decay into γγ for one of the CP-even Higgs bosons with mass ~125 GeV. As a consequence we have not been able to find regions of parameter space where the h±→a1W decay dominates and at the same time are compatible with the possible Higgs signal. This means that irrespectively of whether the possible Higgs signal is substantiated or not, the LHC experiments should be able to either discover or put very tight constraints on a light charged Higgs bosons also in the NMSSM.



---

*Source: 853706-2012-11-08.xml*