Next-to-leading order NMSSM decays with CP-odd Higgs bosons and stops

We compute the full next-to-leading order supersymmetric (SUSY) electroweak (EW) and SUSY-QCD corrections to the decays of CP-odd NMSSM Higgs bosons into stop pairs. In our numerical analysis we also present the decay of the heavier stop into the lighter stop and an NMSSM CP-odd Higgs boson. Both the EW and the SUSY-QCD corrections are found to be significant and have to be taken into account for a proper prediction of the decay widths.

After electroweak symmetry breaking (EWSB) the Higgs sector of the NMSSM consists of seven physical Higgs bosons. In the CP-conserving case, which we assume to be valid here, these are three neutral CP-even, two neutral CP-odd and two charged Higgs bosons. The discovery of all Higgs particles is challenging though not impossible at the high-energy option of the LHC. In [48] we investigated the discovery prospects for the NMSSM Higgs The NMSSM Higgs potential is obtained from the NMSSM superpotential, which we assume to be scale invariant, the soft SUSY breaking terms and the D-term contributions, which are the same as in the MSSM. In terms of the superfieldsĤ u andĤ d , coupling to the up-and down-type quarks, respectively, and the singlet superfieldŜ the superpotential reads where i, j = 1, 2 are the SU(2) L indices and we have introduced the totally antisymmetric tensor ǫ ij with ǫ 12 = 1. Working in the CP-invariant NMSSM, the dimensionless parameters λ and κ are chosen to be real. The MSSM superpotential W MSSM is given by with the soft SUSY breaking MSSM Lagrangian The gaugino fields are denoted byB,W k (k = 1, 2, 3) andG, and the left-handed squarks and sleptons are arranged in doublets denoted byQ = (ũ L ,d L ) T ,L = (ν L ,ẽ L ) T while the right-handed fields are denoted byũ R ,d R andẽ R . The soft SUSY breaking mass parameters m 2 X (R) of the scalar fields X = S, H d , H u ,Q,Ũ ,D,L,Ẽ are real, while the gaugino mass parameters M 1 , M 2 and M 3 and the soft SUSY breaking trilinear couplings A Y (Y = λ, κ, U, D, E) are in general complex. In the CP-conserving case assumed here, they are, however, real. Again, the respective quark and lepton superfields are understood to refer to all three fermion generations. Note that we have set soft SUSY breaking terms linear and quadratic in the singlet field S to zero.
After expanding the Higgs fields around their vacuum expectation values (VEVs) v u , v d and v s , chosen to be real and positive,

JHEP10(2015)024
the Higgs mass matrices for the three scalar, two pseudoscalar and the charged Higgs bosons can be derived from the tree-level scalar potential. The mass matrix decomposes into two mass matrices for the CP-even and the CP-odd Higgs fields. The squared 3×3 mass matrix M 2 S for the CP-even Higgs fields can be diagonalized through a rotation matrix R S which yields the CP-even mass eigenstates H i (i = 1, 2, 3) as (2.6) The H i are ordered by ascending mass, In order to obtain the CPodd mass eigenstates A 1 , A 2 and the massless Goldstone boson G first a rotation R G to separate G is applied, and then a rotation R P to obtain the mass eigenstates which are ordered such that The minimization of the Higgs potential V requires the terms linear in the Higgs fields to vanish in the vacuum. The corresponding coefficients, which are called tadpoles, therefore have to be zero. The tadpole conditions for the CP-even fields can be exploited to replace m 2 Hu , m 2 H d and m 2 S by the tadpole parameters t h d , t hu and t hs . Replacing the SU(2) L and U(1) Y gauge couplings g and g ′ and the VEVs v u and v d by the electric charge e, the gauge boson masses M W , M Z and by tan β, the tree-level NMSSM Higgs sector can then be parameterized by the twelve parameters The VEVs of the neutral components of the Higgs fields are denoted by the brackets around the corresponding fields. The sign conventions for λ and tan β are chosen such that they are positive. The κ, A λ , A κ and µ eff on the other hand can have both signs. Note also, that the parameter A λ can be traded for the charged Higgs boson mass M H ± , which we will do in the following. From now on we will drop the subscript 'eff'. Note that the inclusion of higher order corrections requires also the soft SUSY breaking mass terms for the scalars and the gauginos as well as the trilinear soft SUSY breaking couplings.

The tree-level decay width
We start by discussing the tree-level decay width of a pseudoscalar Higgs boson A i (i = 1, 2) into a pair of stops. The stop mass matrix in the interaction basis (t L ,t R ) reads

JHEP10(2015)024
in terms of the soft SUSY breaking mass parameters mQ and mt R , the soft SUSY breaking trilinear coupling A t , the higgsino mixing parameter µ, the top and the Z boson masses m t and M Z , the mixing angle β and the Weinberg angle θ W . The µ parameter is generated dynamically in the NMSSM and given by The stop mass matrix is diagonalized by yielding the stop mass eigenstatest i (i = 1, 2) as where s = L, R and for the squark masses we have mt 1 < mt 2 by convention. The mixing angle θt and the squark masses are given by and In case of the pseudoscalar only the coupling to two different stop mass eigenstates is non-vanishing. For the tree-level decay width Γ LO we have where λ(x, y, z) = (x − y − z) 2 − 4yz is the two-body phase space function and the coupling G 12 A j (j = 1, 2) of the pseudoscalar A j to the stops reads where v is the VEV given by v 2 = v 2 u + v 2 d and R P jk are the elements of the rotation matrix defined in eq. (2.7). In particular, they are the tree-level mixing matrix elements. In the kinematics of the decay, however, i.e. for the external Higgs field, we use the two-loop corrected Higgs boson masses at order O(α t α s ), which include the full EW corrections at one-loop order. The renormalization of the Higgs fields and the computation of the mass corrections have been described in refs. [103,121,122]. We follow the conventions of these papers, to which we refer the reader for more details. In order to ensure the on-shell properties of the external Higgs field, which in the calculation of the Higgs mass corrections JHEP10(2015)024 has been renormalized in the mixed on-shell-DR scheme, the finite wave function renormalization factors Z [124] have to be taken into account. The application of the factor Z to the tree-level matrix R P (in G 12 A j ) leads to the rotation matrix R P,l which, modulo the Goldstone boson G, rotates the interaction eigenstates a and a s to the loop corrected mass eigenstates A 1 and A 2 , cf. [103], These and the loop-corrected masses are taken from the Fortran code NMSSMCALC [120], in which we choose the on-shell (OS) renormalization for the top/stop sector in order to be in accordance with the renormalization scheme chosen later on both in the SUSY-QCD corrections and in the electroweak corrections. In eq. (3.10) we have summed over both possible final statest 1t * 2 andt * 1t 2 . In the MSSM, i.e. leaving out the singlet contribution ∝ λ in the coupling G 12 For small values of tan β and not too heavy pseudoscalars, the decay into stops can compete with and even dominate over the decays into top quarks and into charginos and neutralinos. In the NMSSM, this statement has to be taken with caution, however, as the singlet component in the coupling G 12 A i , depending on the scenario, can come with both signs and hence increase or decrease the decay width.
The importance of the pseudoscalar Higgs decays into lighter stops depends on the size of the pseudoscalar Higgs coupling to stop pairs and the available phase space as well as the magnitude of other possible pseudoscalar decay widths. Due to the large amount of parameters entering already the tree-level Higgs sector, the identification of regions in the multidimensional parameter space that lead to large branching ratios is a non-trivial task. In general, the pseudoscalar decays discussed in this paper are important for NMSSM setups with a heavy pseudoscalar, not too heavy stops and where the decays of pseudoscalars into SM particles and electroweakinos are suppressed. Such mass configurations appear in scenarios with large charged Higgs masses and stop sectors with large mixing and/or small soft SUSY breaking stop masses. Furthermore, small values of tan β increase the Higgs couplings to stop pairs. The investigated stop decays on the other hand are interesting for scenarios with large stop mass splitting, not too heavy pseudoscalars and suppressed decays with electroweakinos in the final states. Light pseudoscalars can be realized through small values of κA κ .

SUSY-QCD corrections
The SUSY-QCD corrections for the NMSSM pseudoscalar A i decay width differ from the ones of the MSSM [83][84][85] solely in the tree-level coupling to the stops G 12 A i . We shortly repeat them here for completeness and in order to introduce our renormalization scheme.
The virtual corrections at order O(α s ) to the pseudoscalar Higgs decays into stops consist of loop diagrams with a gluon, respectively, gluino exchanged in the A it1t2 vertex and of a contribution involving the four-squark vertex, cf. figure 1 (upper). Note that the mixing contributions due to off-diagonal self-energies are absent in the case of pseudoscalar Higgs decays, as A i only couples to different stop mass eigenstates. The computation of the virtual diagrams leads to ultraviolet (UV) and infrared (IR) divergences. We work in dimensional reduction [125,126], which preserves SUSY at the one-loop level. The fields and couplings are then treated in 4 dimensions while the loop integral is performed in D = 4 − 2ǫ dimensions. The UV divergences appear as poles in ǫ and are canceled by the wave function counterterms and the counterterm renormalizing the A it1t2 interaction. The IR divergences are regularized by the introduction of a fictitious gluon mass ζ. 2 The IR divergences left over after renormalization are canceled after adding the real corrections. These consist of the radiation of an additional gluon off the final state stops and are shown in figure 1 (lower).
The one-loop corrected decay amplitude can be written as where again the Z factors appear to ensure the on-shell properties of the external loopcorrected Higgs field. The ∆ QCD A k are given by the sum of the virtual, real and counterterm From now on a factor C F α s /(4π) with C F = 4/3 is factorized out and already included in eq. (4.2). The virtual corrections receive contributions ∆ g A k from the gluon exchange JHEP10(2015)024 diagram, ∆g A k from the gluino exchange diagram and ∆ 4t from the diagram involving the 4-squark vertex. They are given by for the gluon exchange, for the gluino exchange, and for the 4-squark vertex diagram, where B 0 and C 0 are the Passarino-Veltman scalar twoand three-point functions [127,128], cf. appendix A for their definitions.
The counterterm corrections consist of the renormalization of the external squark wave functions Zt jtj (j = 1, 2) and the renormalization of the A kt1t2 interaction vertex. Note that the wave function renormalization of the pseudoscalar does not contribute at order O(α s ) and hence to the SUSY QCD corrections. The parameters λ, v, v s , tan β, M W , M Z and e are not renormalized by the strong interaction, so that the counterterm reads leading to The counterterm δA t of the trilinear coupling is given by the quark and squark mass counterterms, δm t and δmt 1,2 , and the mixing angle counterterm δθt, (4.9) In the (s)quark sector we adopt OS renormalization with the quark and squark masses defined as the poles of their respective propagators and the squark wave function renormalization constants defined such that the residues of the poles are equal to one. Defining the following structure for the quark self-energy, where P L,R denote, respectively, the leftand right-chiral projector, we have for the top quark mass counterterm The squark mass and wave function counterterms are given by (j = 1, 2) In eq. (4.12) the Σt jtj denote the diagonal parts of the squark self-energies. The diagrams contributing at order O(α s ) to the squark and quark self-energies are depicted in figure 2 (upper) and (lower), respectively.
The mass counterterms read Here B ′ 0 (k 2 ; m 1 , m 2 ) denotes the derivative with respect to k 2 . The mixing angle counterterm is renormalized as where Σt itj denotes the respective squark self-energies, so that The real corrections finally in terms of dilogarithms read where in the two-body phase space function λ( we have neglected the arguments for better readability. We have furthermore introduced We explicitly checked that the dependence on the fictitious gluon mass drops out in the final result, as it should be.
We have cross-checked our results against the ones available in the literature for the MSSM [83,84], both on the analytical and on the numerical level. For the former we performed the MSSM limit of our formulae, and for the latter we took the MSSM limit of the input parameters to evaluate the decay widths.

JHEP10(2015)024 5 The one-loop electroweak corrections
The NLO electroweak corrections consist of the virtual corrections to the vertex and the counterterm contributions to cancel the UV divergences. The top and stop fields as well as the parameters specified below need to be renormalized. For the external Higgs boson we use the two-loop corrected Higgs boson mass at order O(α t α s ), including the full EW corrections at one-loop order, so that we also need to renormalize the Higgs field. In the loops, however, the tree-level masses for the Higgs bosons have to be used so that the UV divergences are canceled properly.
For the NLO EW corrections, in addition the mixings of the decaying CP-odd Higgs boson A i with the Z boson and the Goldstone boson G have to be included. The matrix element for the EW corrected pseudoscalar Higgs decay into stops hence reads where ∆ EW A j represents the sum of the 1-particle irreducible (1PI) diagrams ∆ V,EW A j contributing to the EW virtual corrections of the vertex and of the counterterms ∆ CT,EW Due to massless photons in the loops we also encounter IR divergences. These are canceled by adding the real corrections ∆ R,EW , where a photon is radiated off the final state stop lines. We hence have for the EW corrected decay amplitude Again we have dropped the arguments in the two-body phase space function λ. In the following we will discuss the individual contributions. The virtual corrections consist of the 1PI diagrams given by the triangle diagrams with scalars, fermions and gauge bosons in the loops, as shown in the first two rows of figure 3, and of the diagrams involving fourparticle vertices, cf. figure 3 (last row). For better readability, for the scalars S appearing in the loops we only list the particle types but not the combination of scalars that are allowed by the theory for the various vertices. Let us remark, however, that in the four- compared to the MSSM. The former is due to the singlet admixture in H l , the latter is proportional to the NMSSM specific coupling λ. The diagrams have been generated with the Mathematica package FeynArts 3.6 [129,130] and evaluated with FormCalc 7.3 [131,132] in two independent calculations. The integrals have been computed with LoopTools 2.7 [131,132]. The numerical results of both calculations agree and have been cross-checked numerically against a third calculation, that did not use any of the tools to evaluate and simplify the amplitudes, and which takes the loop functions from HDECAY [133][134][135] and SDECAY [136,137]. The UV divergences encountered in the computation of the virtual corrections are canceled by the counterterms that are the sum of the stop wave function corrections and of the counterterm renormalizing the A jt1t2 interaction, Because of the antisymmetric G 12 A j coupling the stop wave function corrections are given by The stops are renormalized on-shell, with the renormalization conditions given in eqs. (4.12) and (4.13). The diagrams, that contribute to the here required electroweak self-energies are displayed in figure 4 (upper two rows). For the one-loop EW corrections the vertex counterterm reads The individual counterterms are derived from the renormalization of the input parameters. We follow ref. [ The coupling g and the VEV v appearing in eq. (5.7) are given in terms of these parameters by Additionally, for the renormalization of the Higgs sector in the computation of the higher order corrections to the Higgs boson masses, we have the tadpole parameters, the mass of the charged Higgs boson M H ± , the NMSSM parameter κ and the trilinear coupling Aκ, which need to be renormalized at loop level, cf. [103]. Figure 3. Generic diagrams contributing to the electroweak corrections of the decay Γ(A i →t 1t2 ) with i, j, k, p = 1, 2, l = 1, 2, 3 and m, n = 1, . . . , 5. Figure 4. Diagrams contributing to the electroweak squark self-energies (first two rows) and quark self-energies (last row) with i, j, k, p = 1, 2, l = 1, 2, 3 and m = 1, . . . , 5. Figure 5. Generic one-loop diagrams contributing to the mixing δM G,Z mix,i of a pseudoscalar A i with the Z and Goldstone boson. from which their counterterms can be derived. The details of the renormalization of the counterterms for the first six input parameters can be found in [103], so that they are not repeated here. The formulae for the OS renormalization of the top and stop masses are given in eqs. (4.11)-(4.13). The squark and quark self-energies for the EW one-loop corrections are depicted in figure 4. At EW one-loop order the counterterm for A t now reads

JHEP10(2015)024
The counterterm for the mixing angle θt is renormalized as in eq. (4.19), however with the self-energies given by the diagrams shown in figure 4 (first two rows). The diagrams for the contributions to the electroweak corrections stemming from the mixings of the pseudoscalar A i with the Z boson and the Goldstone boson, δM G,Z mix,i , are shown in figure 5. Using the Slavnov-Taylor identity [138,139], on can obtain the mixing contributions through whereΣ A i Z denotes the renormalized A i -Z mixing self-energy and G 12 G is the Goldstone boson coupling to the stops, Note, in particular that the external momenta have to be set to the tree-level mass M (0) A i in order to ensure gauge invariance. The last piece which is missing in the decay, in order to get also an IR finite result, is the real corrections term ∆ R,EW A k . This is the same as for the QCD corrections, but with the gluon replaced by the photon. In the formula for the QCD corrections, eq. (4.21), this means that the coupling and color factors have to be replaced accordingly, i.e.
with ∆ R A k given in eq. (4.21).

JHEP10(2015)024
The full NLO decay width including the SUSY-QCD and -EW corrections is then given by where Γ NLO EW , defined in eq. (5.4), includes the leading order decay width and Γ QCD has been defined in eq. (4.2).

Stop decays at NLO SUSY-QCD and SUSY-EW
The calculation of the higher order corrections to the pseudoscalar Higgs decays into a stop pair can easily be transferred to the decay of a heavy stop into a pseudoscalar and the light stop,t The LO decay width is obtained from eq. (3.10) by adapting the kinematic factor and dividing by the color factor 3 and the factor 2 due to the summation of the two charge conjugated stop pair final states, The NLO SUSY-QCD decay width can be written as with the correction factor ∆ QCD The virtual corrections 4 and the counterterms are given by the same expressions as for the pseudoscalar decay presented in section 4, i.e.
In the real corrections, though, the roles of A k andt 2 have to be interchanged,

JHEP10(2015)024
with ∆ R A k given in eq. (4.21). The NLO SUSY-EW corrections are composed of the same electroweak virtual corrections to the vertex, 5 the same counterterms and the same mixing contributions of the pseudoscalar with the Z and Goldstone boson, ∆ EW A k and δM G,Z mix,i , as in the pseudoscalar decay. In the real corrections, however, again M 2 A i and m 2 t 2 have to be interchanged. We hence have and ∆ EW A k and δM G,Z mix,i given in section 5. The full NLO decay width including SUSY-QCD and -EW corrections is given by QCD given by eqs. (6.4)-(6.7) and Γ NLO EW by eq. (6.8).

Numerical analysis
For our numerical analysis we first perform a scan in the NMSSM parameter space in order to find scenarios that are in accordance with the LHC Higgs and SUSY data. The compatibility with the LHC Higgs data has been checked by using the programs HiggsBounds [140][141][142] and HiggsSignals [143]. The effective couplings of the NMSSM Higgs bosons, normalized to the corresponding SM values, as well as the masses, the widths and the branching ratios of the NMSSM Higgs bosons, which are required as inputs for these programs, have been obtained from the Fortran code NMSSMCALC [120]. The loop induced Higgs coupling to gluons normalized to the corresponding coupling of a SM Higgs boson with same mass is obtained by taking the ratio of the partial widths for the Higgs decays into gluons in the NMSSM and the SM, respectively. These include the QCD corrections up to next-to-nextto-next-to leading order in the limit of heavy quarks [144][145][146][147][148][149][150][151][152][153] and squarks [154,155], taken over from the SM, respectively, MSSM case, while the EW corrections are unknown for the SUSY case and hence consistently neglected also in the SM decay width. The stop mass

JHEP10(2015)024
values have been chosen such that they are not excluded by present ATLAS [156][157][158][159][160][161][162][163] and CMS [164][165][166][167][168][169][170] searches. The squark masses of the first two generations are heavy enough not to be in conflict with LHC data. The SM input parameters that we use are [171,172] α We choose the same renormalization scale µ R as the one used in the calculation of the loop-corrected Higgs boson masses, that enter the decay widths. It is given by the SUSY scale M s , which we set

Pseudoscalar Higgs boson decays into stop pairs
In this subsection we present the impact of the SUSY-QCD and -EW corrections on the decay of a heavy pseudoscalar Higgs boson into a pair of stop quarks. We will investigate the corrections as a function of mt 1 and mg. They have the largest impact on the genuine corrections to the decay widths. The impact of the other parameters is either tiny or they enter the Higgs sector already at tree level with a large impact on the pseudoscalar mass value, so that the genuine effect of the higher order corrections on the decay widths cannot be disentangled any more. The parameter point, which we have chosen from the set of parameter points that survive the LHC constraints, is given by the soft SUSY breaking masses and trilinear couplings mũ R ,c R = md R ,s R = mQ 1,2 = mL 1,2 = mẽ R ,μ R = 3 TeV , mt R = 536.43 GeV , mQ 3  We follow the SLHA format [173,174], in which the parameters λ, κ, A κ , µ eff , tan β as well as the soft SUSY breaking masses and trilinear couplings are understood as DR parameters at the scale µ R = M s , 6 whereas the charged Higgs mass is an OS parameter. As input for our computation we take the soft SUSY breaking trilinear stop coupling A t , however, consistently as an OS parameter. The conversion from the DR to the OS scheme is done within NMSSMCALC and yields 7 A OS t = 1435 GeV . corrections taken into account, as a function of mt 1 , which is varied around the parameter point defined in eqs. (7.1)-(7.5) with mt 1 ≈ 281 GeV. 8 Note, that stops can still be rather light [162,163,170,175,176], down to about 240 GeV for arbitrary neutralino masses [163] and even lower when taking into account the actualt 1 branching ratio [176]. The figure illustrates the effect of the higher order corrections, although the thus obtained parameter configurations are not all in accordance with the applied constraints anymore. The lower plot displays the relative corrections This plot demonstrates that both the inclusion of the EW and the QCD corrections is required in order to properly predict the decay width. In figure 7 we show the branching ratios corresponding to the widths of figure 6. They have been obtained by replacing in NMSSMCALC the corresponding width with our loop corrected width. 9 The higher order corrections to the remaining decay modes needed in the computation of the total width for the branching ratio, are the ones as implemented in NMSSMCALC. In particular they include the dominant higher order QCD corrections. The decays into a bottom quark pair furthermore take into account the higher order SUSY-QCD and the approximate SUSY-EW corrections up to one-loop accuracy. The decays into a strange quark pair include the dominant resummed SUSY-QCD corrections and the ones into a τ pair the dominant resummed SUSY-EW corrections. For details, we refer the reader to ref. [120] and the references therein. The branching ratio at LO of our investigated parameter point amounts to BR LO (A 2 →t 1t2 ) = 40.8% .
(7.11) 8 The variation of mt 1 between 170 and 300 GeV corresponds to a variation of A OS t between 1371 and 1721 GeV. 9 In NMSSMCALC the SUSY QCD corrections to the decays into squarks, as derived from [86], are taken into account. These include improvements in the decays into sbottoms, which are required in parts of the parameter space, that are not relevant for us. Furthermore, we include the EW corrections. We therefore consistently turned off the corrections implemented in NMSSMCALC in the decays into squarks and included instead our corrections. The formulae for the NLO SUSY-QCD and SUSY-EW corrections to the decays of CP-odd NMSSM Higgs bosons into stops shall be included in NMSSMCALC and will be made publicly available in a future release of the program. The net effect of the NLO EW and QCD corrections is an increase of the branching ratio by ∆BR QCD+EW (A 2 →t 1t2 ) = +20.8% . (7.12) In the plot of figure 7 we again vary mt 1 around the chosen parameter point, illustrated by the pink line in the plot. As can be read off the lower plot the total NLO corrections increase the branching ratio by up to a bit more than 40% in the shown mt 1 range. Thus this decay remains the most important one also after the inclusion of the NLO corrections. Figure 8 finally shows the dependence of the higher order corrections on the gluino mass. The EW corrections of course do not depend on mg, while the QCD corrections show a very mild dependence on the gluino mass, apart from the region around mg ≈ 535 GeV. The kink that appears here, arises in thet 2 self-energy at the threshold where mt 2 = mg + m t .
The size of the higher order corrections sensitively depends on the scenario. Thus we find for the scenario defined by This can also be inferred from figure 9 which shows the NLO corrections to the decay widths as a function of mt 1 . As demonstrated in the lower plot, the NLO QCD corrections are of the order of 20-25%, while the EW corrections range between about -17% and -8%, leading to an overall increase of the cross section between 3% and 17% due to the combined NLO corrections.

Stop decays into a pseudoscalar
Decays of the heavy stop into a pseudoscalar and the light stop can occur and become important when there is a large mass splitting between the two stop mass eigenstates.
Stop decays into a SM-like Higgs boson final state have recently been discussed in [177] JHEP10 (2015) Figure 9. Same as figure 6, but now for the initial scenario (marked by the pink line in the plot) given by eqs. (7.13)-(7.14). and the production of NMSSM Higgs bosons in squark cascade decays in [76,178]. In order to show the impact of the SUSY-EW and -QCD corrections on the stop decay width we chose the following parameter set, which leads to an NMSSM Higgs and SUSY spectrum in accordance with the LHC data:  The search for New Physics is one of the main tasks at the LHC. In the absence of any direct sign of new resonances so far, the precise investigation of the Higgs sector becomes more and more important. Physics beyond the SM might reveal itself in modified Higgs decay rates compared to the SM expectations or in the discovery of additional Higgs bosons, unambiguous sign of a non-SM Higgs sector. In view of the complexity of the experimental analyses and the plethora of still possible New Physics extensions, it is evident that the success of this research program depends on the precise predictions of parameters and observables from the theory side. In this paper we have calculated the NLO SUSY-EW and -QCD corrections to the decays of a pseudoscalar NMSSM Higgs boson into stop pairs and of the heavier stop into the lighter stop and a pseudoscalar Higgs boson. Both processes can become important in certain regions of the parameter space and hence contribute to the discovery channels of either the pseudoscalar Higgs boson and/or the stop quarks. The NLO corrections turn out to be important and, depending on the scenario, the EW corrections can be of same size as the QCD corrections and also come with opposite sign. Therefore not only the inclusion of higher order corrections is important but also the consideration of both the QCD and the EW corrections is indispensable for making reliable predictions.

JHEP10(2015)024
with p, p 1 , p 2 denoting the external momenta, which are taken as incoming, and m, m 1 , m 2 , m 3 the masses of the loop particles.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.