One-loop Corrections to the Two-Body Decays of the Charged Higgs Bosons in the Real and Complex NMSSM

We evaluate the full next-to-leading order supersymmetric (SUSY) electroweak and SUSY-QCD corrections to the on-shell two-body decays of the charged Higgs bosons in the framework of the CP-conserving and CP-violating Next-to-Minimal Supersymmetric extension of the Standard Model (NMSSM). Our corrections are implemented in the code NMSSMCALCEW in order to compute the branching ratios of the charged Higgs boson where we also take into account the state-of-the-art QCD corrections already included in the code. We investigate the impact of the NLO corrections for each decay mode in a wide range of the parameter space that is allowed by the theoretical and experimental constraints. The new version of NMSSMCALCEW is made publicly available.


Introduction
While the Higgs boson with a mass of 125 GeV [1] discovered by the LHC experiments ATLAS [2] and CMS [3] behaves very Standard Model (SM)-like [4,5] the pending open questions within the SM call for extensions of the Higgs sector. The lacking direct discovery of any new physics sign so far forces us to focus more and more on the indirect discovery of physics beyond the SM which in turn requires precise predictions of the considered observables.
The very SM-like nature of the discovered Higgs boson calls for sophisticated experimental techniques together with precise theoretical predictions in order to reveal new physics signs. In the context of the NMSSM, various radiative corrections to Higgs observables have been calculated. More specifically in the context of Higgs boson decay widths, the next-to-leading order (NLO) SUSY-electroweak (EW) and SUSY-QCD corrections to the decays of CP-odd Higgs bosons of the CP-conserving NMSSM into stop pairs have been computed in [36]. The full one-loop renormalization of the CP-conserving NMSSM has been worked out in [37,38] together with the computation of the one-loop two-body Higgs decays in the on-shell (OS) renormalization scheme. In [39] a generic calculation of the two-body decays widths at full one-loop level was provided in the DR scheme. The full one-loop corrections for the neutral Higgs decays into fermions and gauge bosons was performed in [40] in the framework of the CPviolating NMSSM and combined with the leading QCD corrections. More recently, the effects of Sudakov logarithms on fermionic decays of heavy Higgs bosons, which appear though radiative corrections of electroweak gauge bosons, have been studied in Ref. [41]. In [42], members of our group completed the evaluation of the NLO SUSY-EW and SUSY-QCD corrections to the full set of two-body on-shell decays of the neutral Higgs bosons in the CP-violating NMSSM, including the state of the art QCD corrections. For the Higgs-to-Higgs decays, the complete one-loop corrections [43] and two loop corrections of order O(α s α t ) [44] have been calculated in both the CP-conserving and the CP-violating NMSSM. As for the charged Higgs bosons, authors of our group have computed the SUSY-EW corrections to the charged Higgs two-body decays into a neutral Higgs boson and a W + boson and studied the gauge dependence arising from the mixing of different loop orders due to the inclusion of mass corrections to the external Higgs boson [45] in order to comply with phenomenological constraints. The authors of [46] discussed this problem further and proposed two strategies to preserve or restore gauge invariance.
In this work, we complete the one-loop higher-order corrections to the on-shell charged Higgs boson decays and compute the NLO SUSY-EW and NLO SUSY-QCD corrections to all charged Higgs boson decays into two-particle final states in the framework of the CP-conserving and CP-violating NMSSM, with the exception of Higgs pair final states that are left for future work. More specifically, we evaluate the corrections to the decays into SM fermions, electroweakinos, sleptons and squarks. The one-loop corrections are based on the renormalization schemes that we have introduced and applied in Refs. [42,44,[47][48][49]. Our newly computed higher-order corrections to the charged Higgs decays as well as those obtained previously in [45] are included in the Fortran code NMSSMCALCEW [42]. The program is based on NMSSMCALC [50] that has been derived by extending the Fortran code HDECAY [51,52] to the NMSSM and as such contains the state-of-the-art QCD corrections and relevant off-shell decays as well as higher-order SUSY corrections through effective couplings. In the first version of NMSSMCALCEW we already included the SUSY-EW corrections to the neutral NMSSM Higgs boson decays of the CP-violating NMSSM published in Ref. [42] as well as their SUSY-QCD corrections to the coloured final states. The code NMSSMCALCEW with now also the higher-order corrections to the charged Higgs boson decays is available at the url: http.//www.itp.kit.edu/∼maggie/NMSSMCALCEW/ With our implemented corrections in NMSMMCALCEW we compute the decay widths and branching ratios of the charged Higgs bosons in the CP-conserving and CP-violating NMSSM and investigate the impact of the SUSY-EW and SUSY-QCD corrections, taking into account the recent theoretical and experimental constraints.
The paper is organized as follows. In section 2, we introduce the NMSSM at tree level to fix our notation for the higher-order computations. In section 3, we briefly discuss the renormalization of the Higgs sector, the electroweakino and the squark sector. In section 4 we move on to the calculations of the radiative corrections to the decay widths of the on-shell two-body decays of the charged Higgs boson. Our numerical analysis is presented in section 5. Conclusions are given in section 6.

The NMSSM at Tree Level
In this section, we briefly introduce the model and set up the notation used throughout the paper. Our notation and conventions follow those already used and explained in detail in our previous works [42-45, 47-49, 53]. We work in the framework of the scale-invariant NMSSM, in which the superpotential is constrained by a Z 3 symmetry. The NMSSM is obtained from the MSSM by adding a gauge-singlet chiral superfieldŜ to the MSSM field content. Since the new additional field is not charged under the SU (2) L × U (1) Y symmetry, the modification of the supersymmetric Lagrangian with respect to the MSSM only arises in the superpotential, whereĤ u ,Ĥ d are the Higgs doublet superfields and ij (i, j = 1, 2) is the totally antisymmetric tensor, with 12 = 12 = 1, and i, j denote the indices of the fundamental SU (2) L representation.
Here and in the following we sum over repeated indices. The dimensionless NMSSM-specific parameters λ and κ are complex in the CP-violating NMSSM. The MSSM superpotential W MSSM is written in terms of the quark and lepton chiral superfieldsQ,L,Û ,D andÊ as where, for simplicity, we neglect colour and generation indices. Assuming flavour conservation, the Yukawa couplings y u , y d and y e are diagonal 3×3 matrices in flavour space. In accordance with the applied Z 3 symmetry, the MSSM parameter µ has been set to zero.
In the NMSSM, the soft supersymmetry (SUSY) breaking Lagrangian reads and H u are the complex scalar components of the chiral superfieldsQ,L,Ĥ d andĤ u , respectively. Similarly,ũ R ,d R andẽ R denote the complex scalar components of the right-handed quark and lepton chiral superfields. Moreover,B,W i (i = 1, 2, 3) andG represent the bino, wino and gluino fields, with masses M j (j = 1, 2, 3), respectively. Finally, the parameters m 2 X are the squared soft SUSY-breaking masses of the fields X = S, H d , H u ,Q,ũ R ,d R ,L,ẽ R . The terms A α (α = u, d, e, λ, κ) are the soft SUSY-breaking trilinear couplings. In the CP-violating NMSSM, the trilinear couplings and the gaugino mass parameters M j (j = 1, 2, 3) are complex.
The Higgs Sector The tree-level Higgs potential derived from the F -and D-terms in the supersymmetric Lagrangian, and the soft SUSY-breaking Lagrangian, Eq.(3), reads where g 1 and g 2 are the gauge couplings of the U (1) Y and SU (2) L symmetry, respectively. The complex scalar Higgs doublet fields H u and H d and the complex singlet field S are expressed in terms of the component fields and vacuum expectation values (VEVs) as where v u , v d and v s are the VEVs of H u , H d , and S, respectively. The two CP-violating phases ϕ u and ϕ s describe the phase differences between the VEVs. Following Refs. [44,[47][48][49], we set the phases of the Yukawa couplings to zero and rephase the left-and right-handed up-quark fields as u L → e −iϕu u L and u R → e iϕu u R , so that the quark and lepton mass terms yield real masses. Inserting Eq. (5) into the gauge sector of the Lagrangian, we obtain the masses the charged gauge bosons W ± and the neutral Z boson, respectively, with the SM VEV v 246 GeV being related to v d and v u as where the three Lagrangian parameters g 1 , g 2 and v have been traded for the three physical observables M W , M Z and the electric coupling e. In the following, the three soft SUSY-breaking mass parameters m 2 H d , m 2 Hu , m 2 S as well as Im A λ and Im A κ will be replaced by the tadpole parameters t φ i (φ i = h d , h u , h s , a d , a s ) 1 . At tree level these tadpoles vanish at the minimum of the Higgs potential. They are, however, non zero at one-loop level, and hence kept for the renormalization procedure. There is also the option to replace the parameter Re A λ by the charged Higgs boson mass.
The Electroweakino Sector The fermionic partners of the neutral Higgs bosons, the neutral higgsinosH u ,H d and the singlinoS, mix with the neutral gauginosB andW 3 , resulting in five neutralinos. The mass term of the neutralinos is given in the basis of the Weyl spinor field with the 5 × 5 symmetric mass matrix M N , By introducing the 5 × 5 orthogonal matrix N , the mass matrix is diagonalized by performing the rotation 2 diag(mχ0 where the matrix N transforms the fields ψ 0 i ≡ (ψ 0 ) i into the mass eigenstates χ 0 i ≡ (χ 0 ) i (i = 1, ..., 5), i.e., The neutralino mass eigenstates are given by the Majorana fieldsχ 0 i = (χ 0 i , χ 0 i ) T (i = 1, ..., 5). By convention, the mass ordering of theχ 0 i is chosen as mχ0 1 ≤ ... ≤ mχ0 5 . We choose M 1 and M 2 as input parameters, so that the masses mχ0 i are derived quantities.
Similarly, the charged higgsinosH ± d andH ± u mix with the charged gauginosW ± resulting in the charginos. In the basis of the spinors u ,H ± d ,W ± , the chargino mass terms are expressed as with 1 The tadpole parameter ta u is not an independent quantity since it is related to ta d as ta u = ta d t β . 2 We choose the mass eigenvalues to be positive, possible complex phases are absorbed into the rotation matrix.
The spinors ψ + L , ψ − R are rotated into the mass eigenstatesχ where U and V are 2 × 2 unitary matrices. Thereby, the mass matrix for the charginos can be diagonalized as with the charginos described by the Dirac spinors (i = 1, 2) The convention used for the mass ordering is the same as that for the neutralinos, namely, mχ± In total, the electroweakino sector requires four additional input parameters which are the absolute values of the gaugino masses and their corresponding two complex phases, hence The Sfermion Sector In the sfermion sector, we only consider the third generation which is the relevant one for our calculation of the charged Higgs boson decays. The mass terms for the left-handed and right-handed sfermions are derived from the soft-breaking Lagrangian Eq. (3) as well as from the F -and D-terms of the supersymmetric Lagrangian. This yields the mass matrices for the stops and sbottoms, with the isospins I t , I b and the electric charges Q t , Q b given by I t = 1 2 , The mass matrix which mixes left-and right-handed staus reads with Q τ = −1, I τ = −1/2. The sneutrino masses are given by with the generation indices i = 1, 2, 3. In our setting, only left-handed neutrinos exist and hence the neutrinos are massless.
Using a 2 × 2 rotation matrix for the sfermions, Uf , which relatesf L ,f R to the mass eigenstatesf 1 andf 2 , we get The sfermion mass matrices are diagonalized as (f =t,b,τ ) In the calculation of the mass matrices Mt, Mb and Mτ , the soft SUSY-breaking masses m 2Q and the trilinear soft SUSY-breaking couplings A t , A b , A τ are chosen as input parameters, so that the mass eigenvalues are outputs, for which the convention mf 1 ≤ mf 2 is used as before.
In addition to the input parameters of the Higgs and electroweakino sectors, we have seven more input parameters for the third generation of the squarks, and four parameters for the third-generation of the sleptons

Renormalization of the NMSSM
In order to obtain UV-finite results at one-loop level, the renormalization of the parameters and external fields is mandatory. In particular, the bare parameters p 0 of the Lagrangian are replaced by the corresponding renormalized parameters, p, and the counterterms, δp, as and the bare fields φ 0 are expressed via the renormalized fields φ and the wave-function renormalization constants (WFRCs) Z φ as In our previous studies [42,44,[47][48][49], we have established the renormalization schemes in the complex NMSSM. We will adopt these procedures for the computation of the NLO corrections to the two-body decays of the charged Higgs bosons and summarize here the main points.
• The tadpole counterterms are chosen in such way that the minimum of the Higgs potential does not change at one-loop level, leading to where t (1) φ i are the tadpole contributions at one-loop level.
• The SM electroweak parameters e, M W , M Z , inspired by their experimental measurements, are renormalized in the on-shell (OS) scheme. The W and Z mass counterterms are given by and the electric charge counterterm reads where Σ T V V , V, V = γ, W, Z, denote the transverse part of the respective self-energies. The notation 'Re' means that we take only the real part of the loop integrals. The above definition of the electric charge counterterm has been chosen to avoid the dependence of the results on large logarithms ln m f from the light fermion f = t contributions [54]. It corresponds to the input of the fine structure constant at the Z boson mass. Therefore, in Eq. (39) the photon self-energy Σ light,T γγ includes only the light fermion contributions.
• The charged Higgs mass in the calculation of the NLO decay widths is always renormalized in the OS scheme and its counterterm is defined via the charged Higgs self-energy, Note that in the computation of the loop-corrected Higgs masses in NMSSMCALCEW if the charged Higgs mass is an input parameter it is renormalized in the OS scheme. In case the DR parameter ReA λ is an input parameter, however, the charged Higgs mass is a derived quantity. The derived value will then be used as input in the decay part as an OS parameter.
• The parameters |λ|, |κ|, v s , tan β, and ReA κ are renormalized in the DR scheme. The counterterms for |λ| ,|κ| and v s can be defined from the electroweakino sector [47,53] or from the neutral Higgs boson sector. Both derivations lead to the same results which has been confirmed explicitly by our previous calculations.
• The gaugino mass parameters |M 1 |, |M 2 | can be renormalized in the OS or the DR scheme. In Ref. [42] we worked with three possible renormalization schemes which we called OS1, OS2 and DR. In OS1, the OS conditions are imposed on the wino-like chargino and the bino-like neutralino, while in OS2 they are applied to the wino-and the bino-like neutralinos. All three schemes lead to small loop corrections to the masses of the remaining neutralinos and charginos, cf. Ref. [42].
• The renormalization of the squark sector, i.e. of the parameters m t , m b , m 2Q and A b can be performed in the OS or the DR scheme. In Ref [42] we provided the counterterms for both schemes and we also presented our evaluation of the one-loop SUSY-QCD and SUSY-EW corrections to the masses of the squarks. We observed large loop corrections to the squark masses and the decay widths of the heavy neutral Higgs bosons into squarks [36,42]. The counterterms of the slepton sector can be derived similarly to the squark sector and are implemented both for the OS and the DR scheme.
• The neutral Higgs fields are renormalized in the DR scheme. The OS property of an external neutral Higgs boson is ensured by using the wave-function renormalization factors (WFRFs) Z H . Thereby potentially large corrections arising from external lines with neutral Higgs bosons can be resummed into the decay widths via WFRFs. The expressions for Z H can be found in e.g. Ref. [45] for the real NMSSM and Ref. [42] for the complex case.
• The charged Higgs fields are also renormalized in the DR scheme. Their WFRCs are related to the WFRCs of the neutral Higgs fields as they are in the same doublet [48,49]. Unlike the neutral Higgs field, we cannot use the resummed WFRFs for external charged Higgs bosons since they contain an infrared (IR) divergence due to the contribution from the massless photon. We therefore expand the charged Higgs WFRFs to take only into account the pure one-loop term which is equivalent to the use of an OS charged Higgs wave function renormalization constant δZ H ± , • The SM fermionic fields and the neutralino and chargino fields are renormalized in the OS scheme defined at the respective tree-level masses of the fields. The counterterm expressions for the neutralinos and charginos can be found in Ref [42]. The OS property of the external fermionic line is justified if the corresponding mass of the fermion is also renormalized in the OS scheme. This is the case for the SM fermionic fields. However, in case of the neutralinos and charginos, some or all of them 3 cannot be renormalized OS. As mentioned in the sixth point in this list, the loop corrections to the masses of the neutralinos and charginos in the studied renormalization schemes, OS1, OS2 and DR, are small. This allows us to use the OS WFRCs defined at the tree-level masses in the decay processes. In the cases where loop-corrections to particles on the external line are large, i.e. squarks, we have to use the modified squark WFRCs which are defined at the loop-corrected masses, [42], for the diagonal WFRCs and for the off-diagonal WFRCs. Here Σq ij denotes the unrenormalized self-energy for theq * i → q * j transition and theΣq stands for the renormalized one. The capital letter Mq is used for the loop-corrected mass while the small letter mq for the tree-level mass. The superscript 'div' means that we take only the UV divergent part. Mq k stands for the loop-corrected mass of the external squark whereq k can be eitherq i orq j . The renormalization of the slepton fields is performed in the OS scheme in accordance with their mass renormalization. Given our flexibility in choosing different renormalization schemes, many particles may get loop-corrections to their masses, e.g. Higgs bosons, neutralinos, charginos, squarks, and sleptons. Loop corrections to the Higgs boson masses have been computed at the two-loop order O(α t α s + α 2 t ) [48,49] while for the neutralinos, charginos, squarks and sleptons loop corrections are calculated at one-loop level [42]. These corrections, depending on the chosen renormalization scheme, have been implemented in the code NMSSMCALCEW [42]. We use the loop-corrected masses in the initial and in the final states of the decay processes. However, we use the tree-level masses and tree-level couplings, if not specified otherwise, in the internal lines of the loop diagrams to ensure the cancellation of the UV divergences.

Higher-Order Corrections to the Two-Body Decays of the Charged Higgs Bosons
In this section we describe the computation of the NLO EW, QCD as well as SUSY-EW, and SUSY-QCD corrections to various two-body decays of the positively charged Higgs boson H + . The NLO corrections for the decay widths of the negatively charged Higgs boson are derived in a similar way. In the real NMSSM the decay widths of the H + and H − are the same, however, in the presence CP-violating phases they may not be equal. In the following, we will discuss the decay modes H + → qq , H + → ν, H + →χ +χ0 , H + →qq , and H + →˜ ν . The NLO SUSY-EW corrections to Higgs boson decays into neutral Higgs bosons and a W + boson as well as their gauge dependences have been discussed in Ref. [45]. Higher-order corrections to the two-body decays of the neutral Higgs bosons in the CP-conserving and CP-violating NMSSM have been presented in Refs. [42][43][44].
We used the public codes FeynArts [55,56], FormCalc [57] and FeynCalc [58,59] in order to generate the Feynman diagrams and calculate the squared amplitudes. The NMSSM model file for FeynArts has been generated with the help of SARAH [60][61][62][63]. We have performed two independent calculations for all processes studied here and the results are in full agreement.
The UV divergent diagrams are regularized by the dimensional reduction scheme [64] which is equivalent to the constrained differential renormalization scheme [65,66] at one loop level as was shown in Ref. [57] and preserves SUSY at this level.

Charged Higgs Boson Decays into Fermions
In our previous implementation in NMSSMCALC, the charged Higgs boson decay width includes higher-order QCD corrections for the quark final states 4 while the decays into lepton final states were evaluated at tree level. We have taken into account, however, the universal large corrections proportional to tan β that arise from the SUSY-QCD (for the top-bottom final state) and the SUSY-EW (for the top-bottom and the τ -neutrino final states) corrections by absorbing them into effective Yukawa couplings as described in our publication [50]. In this section we improve the calculation of the decay widths by including the missing one-loop SUSY-QCD and SUSY-EW corrections, i.e. their finite remainders. We present here the decay into the third generation fermions. Other decays into the first and second generation are kept as implemented in NMSSMCALC, i.e. apart from the absorption into effective couplings where applicable no further one-loop SUSY-QCD and SUSY-EW corrections have been included, since they are suppressed by the smallness of couplings due to the light quark masses and/or small CKM matrix elements. We note that in the setting of NMSSMCALC the CKM matrix is set to unity in the calculation of the loop corrections. This setting is valid for the third generation since the mixing with the first and second generations is small. We, however, keep the corresponding CKM matrix element in front of the charged Higgs decay widths as a factor and do not renormalize it.
The Yukawa interaction terms of the charged Higgs boson with third generation quarks and leptons read where the projection operators P L/R are given by P L/R = (1 ∓ γ 5 )/2 and V tb is the CKM matrix element. The explicit expressions of the tree-level couplings are given by In NMSSMCALC we included higher-order QCD corrections to the decay width H + → tb based on HDECAY [51,52,70]. In HDECAY, there is an interpolation between the decay width that is valid We include the newly computed one-loop SUSY-QCD and SUSY-EW corrections in the latter case. Our implemented loop-corrected decay width is hence given by , with m t/b being pole masses and denotes the two-body phase space factor. Explicit expressions for δ + tb , δ + bt and δ − tb can be found in Ref. [70]. The universal SUSY-QCD and SUSY-EW corrections that become large in the large tan β regime are resummed into effective bottom Yukawa couplings. These ∆ b corrections are incorporated into the decay width through the R factor, where the explicit expression for ∆ b is given by [70][71][72][73][74][75][76][77][78] with the one-loop corrections for the complex NMSSM [50] is the top-Yukawa coupling and C F = 4/3. The generic function I is defined as The remaining one-loop SUSY-QCD and SUSY-EW corrections are denoted by Γ SUSYQCD (47), respectively. The SUSY-EW correction can be expressed as The ∆ L/R,vert SUSYEW contributions are computed from the corresponding one-loop triangle diagrams. 5 The counterterm contributions read where δv is given in terms of the counterterms for e, M W and M Z as, The contributions from the transition of the charged Higgs boson on the external line to the charged W + boson and the charged Goldstone boson, respectively, do not vanish. They are denoted by ∆ L/R,GWmix SUSYEW and given by The last term in Eq.(56) represents the subtraction terms which are needed to avoid double counting arising from the ∆ b corrections. They read with ∆ elw b given in Eq. (52). The contribution from real photon emission, Γ H + →tbγ in Eq. (55), is needed to cancel the infrared divergences (IR) from the triangle contributions. The decay width for the real emission of a photon reads where we have dropped the arguments (M 2 H ± , m 2 t , m 2 b ) of the I functions for better readability. In the above formulae, the I functions are the Bremsstrahlung integrals defined as [54] where k 0 , k 1 and k 2 stand for the incoming momentum of the charged Higgs boson and the outgoing momentum of top and bottom quark, respectively, while k γ is the outgoing momentum of the photon. For the sign of the momentum products k γ · k in and k γ · k jm , the minus sign is chosen only when the indices i n , j m are 0. In the I functions defined above, the integrals I 00 , I 01 , I 02 , I 11 , I 22 , and I 12 are IR divergent. We checked that the analytical expression Eq. (64) is consistent with the generic result for the real corrections to S → FF (S denotes a scalar, F a fermion) calculated in Ref. [39].
The SUSY-QCD contribution, Γ SUSYQCD , can be computed in a similar way as the SUSY-EW corrections, where the left-and right-handed form factors of the renormalized vertex g L/R,SUSYQCD H + bt are given by The contribution ∆ L/R,vert SUSYQCD comes from one-loop diagrams containing a virtual gluino. The counterterm contributions ∆ L/R,CT SUSYQCD are given by Eqs. (57) and (58) by setting the counterterms δv, δt β , δs β , δc β , δZ H ± to zero and including in the computation of δm b,t and δZ t,b L,R instead of the SUSY-EW the SUSY-QCD corrections. The subtraction terms of the SUSY-QCD contributions read with ∆ QCD b and ∆ 1 given in Eqs. (51) and (53), respectively.
The higher-order corrections for the charged Higgs decays into leptons and neutrinos can be expressed in the same way as the decays into quarks, H + → tb. Since the left-handed form factor of the vertex H + lν vanishes the formula for the NLO width is given by where the effective tree-level coupling resumming the large ∆ τ corrections reads and The vertex correction ∆ vert R is given by the corresponding genuine one-loop triangle diagrams. The counterterm contribution, which is needed to cancel the UV divergences arising from the genuine one-loop triangle contribution, is denoted by ∆ CT R and reads For the one-loop H ± -W/G mixing contributions ∆ GWmix R to the effective coupling we get and the subtraction term to avoid double counting reads The real photon emission Γ H + →ντ γ can be cast into the form where the arguments (M 2 H ± , 0, m 2 τ ) of the bremsstrahlung functions have been dropped.

Charged Higgs Boson Decays into Electroweakinos
In this section, we present the results for charged Higgs bosons decaying into a pair of electroweakinos. The Lagrangian describing the interaction between a charged particle H + k (k = 1 for the charged Goldstone boson G + and k = 2 for the charged Higgs boson H + ) and a neutralinõ χ 0 i (i = 1, ..., 5) and a charginoχ − j (j = 1, 2) reads The tree-level couplings are given by with the matrix elements In the previous version of NMSSMCALC, we included only the tree-level decay width which reads Γ tree 6 We now include the NLO SUSY-EW corrections to the decay width as Γ NLO with Γ SUSYEW 6 Note that by capital letters we denote the loop-corrected final state particle masses while the masses of the particles running in the loops are tree-level masses.
The left-and right-handed form factors ∆ L/R consist of the genuine triangle one-loop diagram contributions ∆ vert L/R , the counterterms ∆ CT L/R and the contributions ∆ GWmix L/R from the transitions of H + → G + /W + on the external line, The counterterm contributions read can be expressed in terms of the renormalized self-energyΣ where The IR divergences which appear in the ∆ vert L/R vertex corrections are removed by the contribution from the real photon emission. It is given by

Charged Higgs Decays into Squarks and Sleptons
In this section, we give the expressions for the loop-corrected decays of the charged Higgs boson into squarks and sleptons. The Lagrangian describing the interactions between the charged particle H + k (with H 1 = G + and H 2 = H + ) and the squarks and sleptons reads with i, j = 1, 2 and the tree-level couplings with The fermion Yukawa couplings are defined as The NLO decay width for the process H + →t ib * j consists of the tree-level decay width and the NLO SUSY-QCD 7 and SUSY-EW corrections, where the tree-level decay width was implemented in the NMSSMCALC as Γ tree with the two-body phase space factor The newly computed NLO SUSY-QCD and SUSY-EW corrections are given by The ∆ The contributions to ∆ vert SUSYQCD (∆ vert SUSYEW ) arise from one-loop triangle diagrams containing at least a gluon or a gluino (all possible particles except for gluons and gluinos). The SUSY-EW counterterm contribution can be expressed as where the (δC H + 2bt * ) ab are obtained from Eqs. (96)-(99) by differentiating these expressions with respect to their parameters except for the rotation matrix Z H ± , resulting in with The SUSY-QCD counterterm ∆ CT SUSYQCD is similar to the counterterm ∆ CT SUSYEW , with the modification that the counterterms δg 2 , δv, δt β , δµ eff do not receive QCD contributions and are hence set to zero then, and the remaining counterterms are obtained from SUSY-QCD loop corrections instead of SUSY-EW corrections. The mixing contributions with a W + boson and a charged Goldstone boson are given by Finally the real photon and gluon emission contributions, necessary to cancel the IR divergences arising from the vertex corrections, are given by The NLO decay width for the decay H + →τ * iν (i = 1, 2) is composed of the tree-level contribution and the SUSY-EW corrections and given by with the tree-level decay width Γ tree where the ∆ SUSYEW H +ν ττ * i are given by While ∆ vert SUSYEW denotes the contributions from the one-loop triangle diagrams of the loopcorrected decay H + →τ * iν τ , the counterterm contribution ∆ CT SUSYEW reads where The contribution of the G/W mixing reads and the real photon emission is given by where again the arguments (M 2 H ± , M 2 τ i , M 2 ντ ) of the I functions have been dropped.

Numerical Results
In the following we will discuss the impact of the computed higher-order corrections on the charged Higgs boson decays and branching ratios. In order to get an overall picture we performed a scan in the NMSSM parameter space and only retained those data sets whose phenomenology is in accordance with the most recent experimental results. For this purpose, the parameter points were checked against compatibility with the experimental constraints from the Higgs data by using the programs HiggsBounds [79][80][81] and HiggsSignals [82]. The effective couplings of the Higgs bosons normalized to the corresponding SM values, that are required as input for these programs, have been obtained with the Fortran code NMSSMCALCEW [42]. One of the neutral CP-even Higgs bosons is identified with the SM-like Higgs boson and will be called h from now on. Its mass is required to lie in the range The SM input parameters have been chosen as [83,84] For the NMSSM sector we follow the SUSY Les Houches Accord (SLHA) format [85] in which the soft SUSY breaking masses and trilinear couplings are understood as DR parameters at the scale This is also the renormalization scale that we use in the computation of the higher-order corrections. The code NMSSMCALCEW provides the option to choose either A λ or M H ± as input parameter. We adopted the latter choice and used the charged Higgs boson mass as an OS input parameter. The computation of the O(α t α s + α 2 t ) corrections to the Higgs boson masses is done in the DR renormalization scheme of the top/stop sector [48,49]. In Table 1 we summarize the ranges applied in the parameter scans. In order to ensure perturbativity we apply the rough constraint The bottom trilinear coupling has been fixed to The mass parameters of the first and second generation sfermions are chosen as From the scan we retain those points that have a χ 2 computed by HiggsSignals-2.2.3 that is consistent with an SM χ 2 within 2σ. 8 For the scan and the results presented in the first part of the numerical analysis we keep the CP-violating phases equal to zero. In the second part we turn on various CP-violating phases in order to study their impact individually. Note that in all our CP-conserving scenarios it is the lightest CP-even Higgs boson H 1 that is SM-like and has a mass around 125 GeV.
All the branching ratios shown in the following have been calculated by implementing the higher-order corrections to the various charged Higgs boson decay widths in NMSSMCALCEW. In this way the newly computed corrections are combined with the state-of-the-art higher-order QCD corrections already included in NMSSMCALCEW. Note, however, that the (SUSY-)EW and SUSY-QCD corrections are only taken into account if the respective decay is kinematically allowed. Otherwise, the corresponding decay width without the higher-order corrections discussed in this paper, which only apply for on-shell decays, is used in the computation of the total decay width and branching ratios.
For the computation of radiative corrections we used the following renormalization schemes unless stated differently (see also Section 3): In the electroweakino sector we used the OS1 renormalization scheme. For the SUSY-EW corrections to the decays into stop-sbottom and stau-sneutrino pairs the DR scheme was applied for the stop, sbottom sector and the OS scheme for the stau sector. For details on the definition of the schemes, we refer to Ref. [42].

The Branching Ratios
In Fig. 1 we show as function of the charged Higgs boson mass the NLO branching ratios of the charged Higgs boson decays into the various possible final states, namely into the SM-like final states tb, τ + ν τ , µ + ν µ and cs, and the new physics final states W + H,χ + 1χ 0 ,χ + 2χ 0 ,tb * andτ * ν τ . In the decays into a charged W + boson plus Higgs final state we have summed up all W + H i final states, so that the branching ratio BR(H + → W + H) is given by 9 Analogously, in the branching ratios into the electroweakino final states we have summed over the neutralinos, hence Also the decays into the sfermion final states are summed over so that The NLO branching ratios include the higher-order corrections to the Higgs decays widths as presented in the draft, namely the SUSY-EW corrections as well as the QCD and SUSY-QCD corrections for the coloured final states. More specifically, the formulae for the loop-corrected decay widths are given in Eq. (47) for the decay in the top-bottom final state, in Eq. (70) for the decay into τ + ν τ , in Eqs. (83) for the decays into electroweakino pairs, in Eqs. (101) for the decays into stop-sbottom pairs, and in Eq. (116) for those into stau-sneutrino pairs. The implemented higher-order corrections to the decays into charged W + boson plus Higgs final states have been described in Ref. [45] for the CP-even Higgs bosons in the CP-conserving NMSSM. We have extended this to the CP-violating case. The decays into the SM-like final states include the higher-order QCD and resummed SUSY corrections as specified in the manual for HDECAY [51,52] and extended to the NMSSM case in [50]. As mentioned above the SUSY-EW and SUSY-QCD corrections are only included in on-shell decays. Otherwise, (where applicable) only QCD corrections and resummed corrections through effective couplings are included. We furthermore include in the decays with on-shell neutral Higgs bosons in the external states the resummed Z H factors, cf. Sec. 3.
As can be inferred from the plots, the largest branching ratios are given by the decays into top-bottom final states (gray circles) with values of up to almost 100% for charged Higgs mass values below 2 TeV. For larger values of the charged Higgs mass the decays into electroweakinos become dominant. For M H ± < 1.2 TeV, the decays into a charged W plus a neutral Higgs boson provide the second largest branching ratio for some parameter points. As stated above we show here the sum over all possible neutral Higgs bosons. The resulting branching ratio, indicated by the orange circles, can reach up to 98%. The decays into the electroweakinos (cyan diamonds for χ + 1χ 0 , violet diamonds forχ + 2χ 0 ) can reach up to 54% (55%) forχ + 1χ 0 (χ + 2χ 0 ) when summed up. The summed-up branching ratios into stau-sneutrino pairs (blue diamonds) can go up to 17% for charged Higgs masses below 1.4 TeV and specific parameter configurations, whereas the decays into stop-sbottom pairs (black diamonds) become important for large charged Higgs masses and can have branching ratios of up to 29% in their sum. The branching ratios into τ + ν τ (pink circles) reach 20%. The branching ratios into µ + ν µ (red circles) and cs (green circles) attain at most 7 · 10 −2 % and 4 · 10 −2 %, respectively. The comparison of these scatter plots with the corresponding ones for the leading order (LO) branching ratios shows that the overall pattern of the distribution of the branching ratios does not change when NLO corrections are included. For individual parameter points the changes can be substantial, however. In the following, we will discuss the impact of the higher-order corrections for the various final states separately.

Impact of Higher-Order Corrections
For the discussion of the impact of the NLO corrections on the decay width of the decay H + → XY we introduce the relative correction of the partial width as We furthermore define the relative change in the branching ratio for the decay H + → XY as .
This quantity allows us to directly identify large corrections in the branching ratios that are not 'artificially' enhanced because of tiny LO branching ratios.
We have to specify what we mean by LO widths and branching ratios. They are the LO quantities calculated with 'Higgs effective tree-level couplings', which means that the Higgs tree-level rotation matrix elements have been replaced by the loop-corrected ones. For decays with neutral Higgs bosons in the final state, additionally the improved resummed Z H factor as described in Ref. [42] is included in the LO decay widths and branching ratios. Note, that these 'LO' quantities also include the QCD corrections and resummed SUSY-EW and SUSY-QCD corrections in effective quark couplings as already implemented in the first release of NMSSMCALC [50] and described there. In fact the use of the word 'LO' in essence means that we thereby refer to the old implementation in NMSSMCALC without the genuine SUSY-EW and SUSY-QCD vertex corrections. This means that the definitions Eqs. (134) and (135) give us information on the effects of the newly computed corrections, namely the SUSY-EW and SUSY-QCD vertex corrections, respectively their finite remainders, compared to the previous implementation in NMSSMCALC which only uses the improved LO decay widths as defined here. Figure 2 displays the relative, ∆ BR , due to the impact of the SUSY-QCD and SUSY-EW corrections on the branching ratio BR(H + → tb) (left) and the impact of the SUSY-EW corrections on BR(H + → τ + ν τ ) (right) as a function of their respective NLO branching ratios. The color code indicates the respective relative corrections to the partial width in per cent for the newly computed SUSY-QCD and SUSY-EW corrections. It shows that the impact of the corrections on the partial width for the decay into tb is of moderate size, ranging between -20% to +2% with a few outliers going down to -29%. The relative change in the branching ratio is of similar size with values between -20% and +8% and a few outliers going down to about -30%. Splitting up the contributions, we find that apart from a few outliers the relative SUSY-EW corrections to the partial width (branching ratio) lie between -16% and -2% (-12% and +6%), whereas the relative SUSY-QCD corrections range between -11% and +7% (-10% and +4%). The SUSY-EW corrections on the width are negative and of comparable size as the SUSY-QCD ones which underlies the importance of including both types of corrections. Large relative corrections to the decay widths of up to about -30% arise from the sum of SUSY-QCD and SUSY-EW corrections with same sign. The dominant contributions to both the SUSY-EW corrections and the SUSY-QCD corrections stem from one-particle irreducible triangle diagrams. The SUSY-EW corrections to the decay width into τ + ν τ mostly lie between -17% and +7% and between -10% and +15% for ∆ BR for the bulk of the points.  135) for the H + decays into tb (left) and and τ + ν τ (right) as a function of the respective NLO-corrected (SUSY-QCD and SUSY-EW for the former and SUSY-EW for the latter decay) branching ratio for all scan parameter points passing the constraints. The color code indicates the relative correction in the partial width as defined in Eq. (134).

Decays into Gauge plus Higgs Boson Pairs
In Fig. 3 we show the relative change in the branching ratios due to our newly computed genuine SUSY-EW corrections and as colour code the relative correction for the partial decay widths of the H + decays into charged W + boson plus Higgs final states as a function of the corresponding NLO branching ratio. We do not classify the Higgs final states by the mass eigenstates but by the gauge eigenstates, i.e. we show final states with mostly h u , h s and a s final states (decays into W + h d and W + a d are kinematically closed). Mostly j-like (j = h u , h s , a s ) means that the mixing matrix element squared |R ij | 2 of the Higgs eigenstate H i exceeds 0.5. Note, that the h u -like state corresponds to the SM-like Higgs boson as compatibility with the Higgs data requires a maximum coupling to the top-quark. As can be inferred from the upper left plot the bulk of the corrections to the branching ratio (decay width) for the decay into an h u -like, i.e. SM-like, Higgs boson together with the charged W + boson lies between -42% and +10% (-35% and +3%). There are a few outliers with somewhat larger corrections. We note that the branching ratio into W + h u always remains below 0.75% and is hence rather unimportant for charged Higgs decays.
The branching ratios into singlet-like CP-even Higgs final states, W + h s , with up to 25% reach larger values than those into W + h u . The relative corrections to the decay widths are moderate and range between -18% and +4% with the change in the branching ratio being between -14% and +10%. There are two outliers with ∆ BR reaching up to -20%. Similarly, the corrections to the decay width in the CP-odd singlet-like Higgs state are moderate with corrections to the decay width between -10% and +15% and to the branching ratio between mostly -2% and +20%. A few outliers can reach corrections to the branching ratio of up to -23%.
Note that for all scattering plots here we do not consider branching ratios that are smaller than 10 −4 , since they are not phenomenologically interesting. We can have large relative corrections in these cases because of the suppression of the tree-level couplings.

Decays into Electroweakinos
We now turn to the impact of the SUSY-EW corrections on the decays into electroweakino pairs. The relative changes in the branching ratios and the relative corrections of the decay widths are shown for the final states in the gauge basis, in Fig. 4 for the charged wino and in Fig. 5 for the charged higgsino final state, respectively, together with a neutral electroweakino as specified in the figure labels. In the plots we show results for the parameter points passing our constraints after applying the following cuts: We cut the ratio to lie in the range This ensures that there are no large hierarchies among the left-handed couplings 10 of the charged Higgs to a neutralino-chargino pair, cf. Eq. (80), which would otherwise blow up the NLO corrections compared to the LO width. If r ∈ [0.5... 1.5], there exist cancellations in the treelevel couplings of the decay in question, that lead to a suppression of the tree-level decay width. Contributions coming from a neutralino close in mass with an enhanced tree-level coupling will dominate and can be very large. In NMSSMCALCEW we print out a warning if this case occurs. Furthermore, for the charged Higgs decay H + →χ + iχ 0 j we impose the following cuts on the mass differences to avoid large mixing effects between the close-in-mass electroweakino masses that also induce huge NLO corrections. Since we use the OS scheme for the WFR constants of the electroweakinos, if two neutralinos (charginos) are degenerate then the WFR constant contributions are dominant and huge. 11 In the code we print out a warning if degenerate cases are involved. In our study, we fix the renormalization scheme for all points to be OS thus encountering about one hundred points 12 the corrections have the typical size of EW corrections that we comment on in the following. In case of unnaturally large loop corrections we recommend to change the renormalization scheme.
The maximum relative corrections do not differ much for the charged wino and charged Higgsino final states when comparing the final states with the same neutral electroweakino. The smallest corrections are found for neutral down-type Higgsino production together with a charged wino or Higgsino, i.e.W +H 0 d andH + uH 0 d . The relative corrections of the partial width lie between about -16% and +24% for the former and -37% and +2% for the latter, for the bulk of the points. A few outliers involve also corrections of up to -55% to theH + uH 0 d decay. The relative corrections of the branching ratios lie between about -4% and +20% forW +H 0 d production and -30% and 0% forH + The relative corrections to the decay widths lie between about -35% and +40% (depending on the specific final state) for the bulk of the parameter points, and those to the branching ratios between about -30% and +30% apart from a few outliers. The largest corrections are found forW +B ,H + uB production where the relative corrections to the partial widths range between -34% (-34%) and +77% (+57%) for W +B (H + uB ) and to the branching ratios between -30% and +40% barring a few outliers that can go up to -100%. These are found, however, for small LO widths and branching ratios so that the relative correction easily gets enhanced.
Let us also briefly comment on the size of the branching ratios. The largest branching ratios are obtained for theH + uW 3 andW +W 3 final states with 37% and 30%, respectively, followed byH + uB (35%) andW +B (25%) production. The branching ratios intoW +H 0 d ,W +H 0 u ,W +S , H + uH 0 u ,H + uS reach maximum values between 10 and 20%. The smallest branching ratio is found for theH + uH 0 d final state with at most 0.5%.

Decays into Sfermions
Like for the electroweakino decays also in the decays into sfermions we apply some cuts to get rid of artificially enhanced corrections due to specific corners in the parameter space. Thus we require for the decays into stop-sbottom pairs, H +t ibj , that the involved masses fulfill so that large corrections originating from threshold effects in triangle loop diagrams do not play a role. Furthermore, we apply the following cuts on the mass differences 13 This way large corrections arising from degenerate squark mixing contributions are removed. In case of large loop corrections, we recommend the user to change the renormalization scheme of 12 These points are in the decays of H + intoW +B ,W +S ,H + uB andH + uS . with relative corrections beyond 100% among 10,000 allowed points. Without the applied cuts the scattering plots will involve large scales so that it becomes difficult to read off the corrections of most of the points that have mild corrections. After applying these cuts 13 We remind the reader that the final state particle masses receive loop corrections as they are renormalized DR whereas the masses of the particles in the loop are taken at tree level. Capital letters denote loop-corrected, small letters denote tree-level masses.  are shown in Fig. 6 for the stop-sbottom final state on the left and for the stau-sneutrino final state on the right. Note that we sum over all possible final statest ib * j (i, j = 1, 2) in the left and τ * iν τ (i = 1, 2) in the right plot. The reason why we have less points in the stop-sbottom than in the stau-sneutrino case is simply because in a lot of the scenarios of our scan compatible with the applied constraints the former decays are kinematically closed.
The decays into squarks receive SUSY-EW and SUSY-QCD corrections. For most of the points the relative SUSY-QCD corrections to the decay widths are distributed between 10% and 90%, while the typical size for the SUSY-EW corrections is between -30% and +22%. The combination of both corrections alters the decay width between -37% and +92%. Some parameter points can reach corrections of up to more than 100% which is due to the renormalization of the squark wave functions in the SUSY-QCD corrections and calls for future improvements through the inclusion of higher-order corrections or resummation. Barring these outliers the relative corrections in the branching ratios range between -10% and +50%.
The relative corrections to the decays into stau-sneutrino, stemming only from SUSY-EW corrections, lie between -64% and +72%, ranging for the bulk of the corrections only between -31% and +39%, however. For ∆ BR we find values between -25% and 30%.

The Impact of CP Violation
For the investigation of the impact of CP violation on the higher-order corrections to the charged Higgs decays we chose the following parameter point from the set of valid scenarios obtained in our scan, The remaining parameters and phases are fixed as where for simplicity we choose the notation ϕ µ eff ≡ ϕ µ .
In the following we vary the three phases ϕ µ , ϕ M 2 , ϕ At individually away from their benchmark value zero while keeping all the other phases fixed to zero in order to quantify the effect induced by CP violation through the respective non-zero phase. For the phase of λ, we set ϕ λ = 2ϕ µ /3 so that we do not encounter CP violation at tree level in the Higgs sector. Note, that we vary the phases almost up to their maximum values leading to scenarios that are already excluded by the EDM constraints. For illustrative purposes, we still allow for these variations, however. We stop all plots at ±0.47π and not at ±π/2 because for phases |ϕ µ | > 0.47π the singlet-like Higgs boson has a negative mass. In Fig. 7 (left) we show the tree-level and loop-corrected decay width of the charged Higgs decay into the top-bottom final state as function of a variation of either ϕ µ (red), ϕ At (blue) or ϕ M 2 (green) while keeping all other phases to zero and ensuring a vanishing tree-level CPviolating phase in the Higgs sector. The lower panels show the relative corrections of the partial decay width as defined in Eq. (134) as a function of the respective non-zero phase. Despite the vanishing CP-violating phase at tree level we still see a small dependence of the tree-level decay width on the CP-violating phase ϕ µ and ϕ At , respectively. This is due to the ∆ b corrections included in our definition of the tree-level decay width into quarks, which depends on the phases of µ eff and A t . As for the loop corrections, the relative correction which is δ Γ = −19% for the chosen parameter point, barely changes when ϕ M 2 is turned on, whereas for non-zero ϕ µ it varies from -19% to −11% for ϕ µ = −0.47π and −12.5% for ϕ µ = +0.47π. The dependence on ϕ At is given by values ranging from -19% at vanishing phase to −14% at ϕ At = ±0.47π.  is the SM-like Higgs boson, cf. Fig. 7 (right). For ϕ µ and ϕ At each, the relative corrections vary from about -10% to at −0.47π to about +40% at +0.47π. The dependence on ϕ M 2 is very weak, however. The LO dependence on ϕ µ and ϕ At is due to the inclusion of the two-loop corrections into the final state Higgs boson masses, where the dependence on ϕ µ is largest. For the W + H 2 final state ϕ µ and ϕ M 2 change δ Γ from 3% at zero phase to 4.5% (1% for ϕ µ and 1.5% for ϕ M 2 ) at 0.47π (−0.47π). The dependence on ϕ At is very small. The loop corrections to H + → W + A 1 show a smaller dependence on ϕ µ , ϕ M 2 and again the dependence on ϕ At is almost negligible. The other decays into charged W + boson plus Higgs final states are kinematically closed.
The charged Higgs couplings to the electroweakinos depend on the phases ϕ µ and ϕ M 2 so that the LO decay widths already show a dependence on these CP-violating phases. We exemplary show the effect on the decay H + →χ + 1χ 0 2 in Fig. 8 (left). In this case theχ + 1 is a higgsino-like chargino while theχ 0 2 is a higgsino-like neutralino. (The CP-violating impact is found to be less important for theχ + 1χ 0 1 /χ 0 3 final states of this benchmark point.) The relative correction δ Γ shows a substantial dependence on all three phases.
The tree-level charged Higgs couplings to squarks depend on the phases ϕ µ and ϕ At so that the LO width shows a substantial dependence on these phases that is also translated to NLO, as can be inferred from Fig. 8 (right), which shows the decay H + →t 1b * 1 as a function of the CP-violating phases at LO and NLO. (The decays into other stop-sbottom final states are kinematically closed.) The dependence on ϕ M 2 is only radiatively induced and very weak.

The Impact of the Renormalization Scheme
We now want to discuss the impact of the renormalization scheme. The renormalization scheme dependence of the decay widths also gives us a possibility to roughly estimate the remaining theoretical error due to missing higher-order corrections. The code NMSSMCALCEW that we use for the computation of the decay widths and branching ratios follows the SLHA conventions where the soft SUSY breaking parameters are understood to be DR input parameters at the scale M SUSY which per default is given as in Eq. (126). Consequently, depending on the chosen     Table 4: The tree-level and one-loop corrected stau and tau sneutrino masses in the DR and the OS scheme. renormalization scheme for the computation of the higher-order corrections the input parameters have to be converted to the applied scheme where necessary. For example, if we use OS renormalization in the stop sector then the soft SUSY-breaking parameters (m 2 Q 3 ,m 2 t R , A t ) affecting the stop sector must be converted from the SLHA DR parameters to OS parameters. With these converted parameters the NLO width is calculated then. In Ref. [42] we outlined in detail the procedure applied in NMSSMCALCEW to convert the input parameters. In the following we show the change in the LO and NLO widths for the electroweakino and sfermion decays when applying different renormalization schemes after consistently converting the input parameters. The benchmark point used in the plots is the same as for the investigation of the impact of CP-violation in the previous subsection 5.7, we only vary tan β in the following. In order to quantify the renormalization scheme dependence we introduce where Γ OS (Γ DR ) denotes partial width evaluated in the OS (DR) scheme. It should be noted that in both renormalization schemes we use the loop-corrected masses for the external lines only while we use tree-level masses and tree-level couplings for particles inside loops. For our chosen parameter point, we present the tree-level and loop-corrected masses for the electroweakinos in Table 2, for the stops/sbottoms in Table 3, and for the staus/tau-sneutrino in Table 4.  , respectively, for OS1 renormalization 14 (red) and for DR renormalization (blue) at LO (dashed) and NLO (full) as a function of tan β. The lower insert displays ∆Γ where in the definition Eq. (143) OS has to be replaced by OS1. As we can inferred from the plots for all three decays the dependence on the renormalization scheme decreases when going from LO to NLO, as expected. For theχ + 1χ 0 1 andχ + 1χ 0 3 final states, already at LO the dependence is very small with ∆ Γ ≈ 1.5% for the former and 0.32% for the latter, almost independently of tan β. It gets reduced to close to 0% at NLO. For theχ + 1χ 0 2 final state the LO difference between OS1 and DR renormalization is larger than in the other decays, as is the relative NLO correction to the decay width in the OS1 scheme. For this final state the LO dependence on the renormalization scheme gets reduced from about 18% to 3% at NLO. These results are in accordance with the observed small dependence of the loop-corrected final state masses on the renormalization scheme, cf. Table 2. Figure 10: The partial decay width for the decay H + →t 1b * 1 as a function of tan β for OS (red) and DR (blue) renormalization in the stop sector at LO (dashed) and NLO (full) for the SUSY-EW (left) and the SUSY-QCD corrections (right). Lower insert: The relative difference in the widths due to different renormalization schemes at LO (dashed) and NLO (full).
In Fig. 10 we display the LO and the NLO widths for the charged Higgs decay intot 1b * 1 for OS renormalization (red) and for DR renormalization in the stop sector (blue) at LO (dashed) and NLO (full) for the SUSY-EW (left) and the SUSY-QCD corrections (right), as a function of tan β. The lower insert displays ∆Γ as defined in Eq. (143). As expected, the plots show that the LO dependence on the renormalization scheme gets reduced when going to NLO, both for the SUSY-EW and the SUSY-QCD corrections, namely from about 38% (39%) to 35% (37%) for the SUSY-EW corrections for tan β = 4.5 (6.5). For the SUSY-QCD corrections the reduction is more important, going down to about 30%. The overall larger dependence on the renormalization schemes compared to the electroweakino final states is also reflected in the larger dependence of the final state squark masses on the renormalization scheme, cf. Table 3. Figure 11 shows the H + →τ * ν τ decay (summing over the two possibleτ i (i = 1, 2) states) for OS (red) and the DR (blue) renormalization at LO (dashed) and SUSY-EW NLO (full) as a function of tan β. For the chosen parameter point the SUSY-EW corrections are rather small, and the rather small dependence on the renormalization scheme at LO with ∆ Γ = 5% gets reduced to 1% at NLO.

Conclusions
In this paper we complete the computation of the NLO SUSY-EW and SUSY-QCD corrections to the charged Higgs boson decays in the CP-conserving and CP-violating NMSSM that was started in a previous paper with the corrections to the charged W + plus Higgs boson final states. We provide the missing corrections to the on-shell two-body decays into the SM fermion, the Figure 11: The partial decay width for the decays into stau-sneutrino pairs as a function of tan β for OS (red) and DR (blue) renormalization at LO (dashed) and NLO (full). Lower insert: The relative difference in the widths due to different renormalization schemes at LO (dashed) and NLO (full). electroweakino and the sfermion final states. For the decays into electroweakinos and sfermions we provide the corrections for different renormalization schemes, chosen to be OS and DR. This allows us to roughly estimate the remaining theoretical uncertainty due to missing higher-order corrections. All corrections have been implemented in the code NMSSMCALCEW thereby combining them with the already incorporated state-of-the-art higher-order QCD corrections.
In our numerical analysis we find that the newly computed SUSY-EW and SUSY-QCD corrections are significant and need to be included for meaningful predictions of the decay branching ratios. In specific corners of the parameter space with small LO widths and/or large mixing effects the corrections can become very large. As for the dependence on the CP-violating phases of the various parameters the effects are of typical size for radiatively induced CP violation. The investigation of the renormalization scheme dependence shows a good perturbative convergence of the higher-order corrections. Our results contribute to the improvement of the precision on the predictions for Higgs boson observables in beyond-the-SM extensions that is required for the correct interpretation of new physics effects being looked for at the LHC.