One-Loop Corrections to the Two-Body Decays of the Neutral Higgs Bosons in the Complex NMSSM

TL;DR

This paper calculates precise one-loop electroweak corrections to the two-body decays of neutral Higgs bosons in the complex NMSSM, enhancing the accuracy of theoretical predictions for collider searches.

Contribution

It provides the complete one-loop electroweak corrections to neutral Higgs decays in the CP-violating NMSSM, implemented in the new NMSSMCALCEW code, improving prediction precision.

Findings

Enhanced precision in Higgs decay predictions.

Implementation of electroweak corrections in NMSSMCALCEW.

Supports collider searches for new physics.

Abstract

Since no direct signs of new physics have been observed so far indirect searches in the Higgs sector have become increasingly important. With the discovered Higgs boson behaving very Standard Model (SM)-like, however, indirect new physics manifestations are in general expected to be small. On the theory side, this makes precision predictions for the Higgs parameters and observables indispensable. In this paper, we provide in the framework of the CP-violating Next-to-Minimal Supersymmetric extension of the SM (NMSSM) the complete next-to-leading order (SUSY-)electroweak corrections to the neutralHiggs boson decays that are on-shell and non-loop induced. Together with the also provided SUSY-QCD corrections to colored final states, they are implemented in the Fortran code NMSSMCALC which already includes the state-of-the art QCD corrections. The new code is called NMSSMCALCEW. This way we provide the NMSSM Higgs boson decays and branching ratios at presently highest possible precision and thereby contribute to the endeavor of searching for New Physics at present and future colliders.

Tables

Table 1. Table 1: Input parameters for the NMSSM scan. All parameters have been varied independently between the given minimum and maximum values.

	$t_{\beta }$	$\lambda$	$\kappa$	$M_1$	$M_2$	$M_3$	$A_t$	$m_{{\tilde{Q}}_3}$	$m_{{\tilde{t}}_R}$	$M_{H^{\pm }}$	$A_{\kappa }$	$\left|{\mu }_{\text{eff}}\right|$
				in TeV
min	1	0	-0.7	0.5	0.5	1.8	-3	0.6	1	0.5	-2	0.2
max	20	0.7	0.7	1	1	2.5	3	3	3	1.5	2	1

Table 2. Table 2: Scenario1: Masses and main components of the neutral Higgs bosons at two-loop order 𝒪 ​ ( α t ​ α s + α t 2 ) 𝒪 subscript 𝛼 𝑡 subscript 𝛼 𝑠 superscript subscript 𝛼 𝑡 2 {\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2}) , using DR ¯ ¯ DR \overline{\text{DR}} renormalization in the top/stop sector.

$\overline{\text{DR}}$	$H_1$	$H_2$	$H_3$	$H_4$	$H_5$
two-loop $𝒪\left({\alpha }_t{\alpha }_s+{\alpha }_t^2\right)$	125.14	698.05	813.53	1391.6	1392.56
main component	$h_u$	$a_s$	$h_s$	$a$	$h_d$

Table 3. Table 3: Scenario1: Masses and main components of the neutralinos and charginos at tree and one-loop level in the three renormalization schemes OS1, OS2 and DR ¯ ¯ DR \overline{\text{DR}} .

		$M_{{\tilde{\chi }}_1^0}$	$M_{{\tilde{\chi }}_2^0}$	$M_{{\tilde{\chi }}_3^0}$	$M_{{\tilde{\chi }}_4^0}$	$M_{{\tilde{\chi }}_5^0}$	$M_{{\tilde{\chi }}_1^{+}}$	$M_{{\tilde{\chi }}_2^{+}}$
OS1	tree-level	265.97	276.05	565.2	730.92	920.76	270.72	730.83
	one-loop	273.22	282.48	565.2	730.78	914.39	278.02	730.83
OS2	tree-level	265.97	276.05	565.2	730.78	920.76	270.72	730.69
	one-loop	273.22	282.48	565.2	730.78	914.39	278.02	730.83
$\overline{\text{DR}}$	tree-level	265.5	276.2	563.47	694.32	920.76	270.38	694.18
	one-loop	273.21	282.44	565.21	730.01	914.39	278.01	729.86
	main component	${\tilde{H}}_d^0$	${\tilde{H}}_u^0$	$\tilde{B}$	${\tilde{W}}_3$	$\tilde{S}$	${\tilde{H}}^{+}$	${\tilde{W}}^{+}$

Table 4. Table 4: Scenario 1: Loop-corrected decay widths Γ l superscript Γ 𝑙 \Gamma^{l} , relative loop corrections δ 𝛿 \delta and loop-corrected branching ratios BR l of all kinematically allowed decays into chargino and neutralino pairs of heavy Higgs bosons in the OS1 and DR ¯ ¯ DR \overline{\text{DR}} renormalization schemes. The results in the OS2 scheme are nearly the same as in OS1 scheme.

		${\Gamma }^l\left[\text{MeV}\right]$	$\delta [\%]$	BRl [$\%$]		${\Gamma }^l\left[\text{MeV}\right]$	$\delta [\%]$	BRl [$\%$]
$H_2\to {\tilde{\chi }}_1^0{\tilde{\chi }}_1^0$	OS1	405.67	0.08	25.14	$H_2\to {\tilde{\chi }}_1^0{\tilde{\chi }}_2^0$	0.43	-0.08	0.03
	$\overline{\text{DR}}$	407.74	8.4	25.15		0.4	-16.34	0.02
$H_2\to {\tilde{\chi }}_2^0{\tilde{\chi }}_2^0$	OS1	404.44	0.07	25.06	$H_2\to {\tilde{\chi }}_1^{+}{\tilde{\chi }}_1^{-}$	802.3	0.08	49.72
	$\overline{\text{DR}}$	404.84	7.23	24.98		807.18	9.03	49.8
$H_3\to {\tilde{\chi }}_1^0{\tilde{\chi }}_1^0$	OS1	321.53	0.08	22.01	$H_3\to {\tilde{\chi }}_1^0{\tilde{\chi }}_2^0$	0.78	-0.11	0.05
	$\overline{\text{DR}}$	322.45	8.13	22.04		0.73	-16.24	0.05
$H_3\to {\tilde{\chi }}_2^0{\tilde{\chi }}_2^0$	OS1	308.94	0.07	21.15	$H_3\to {\tilde{\chi }}_1^{+}{\tilde{\chi }}_1^{-}$	633.52	0.08	43.38
	$\overline{\text{DR}}$	309.0	7.08	21.12		635.37	8.79	43.42
$H_4\to {\tilde{\chi }}_1^0{\tilde{\chi }}_1^0$	OS1	102.86	0.02	0.8	$H_5\to {\tilde{\chi }}_1^0{\tilde{\chi }}_1^0$	151.31	-0.01	1.18
	$\overline{\text{DR}}$	89.4	-11.31	0.69		132.77	-13.56	1.04
$H_4\to {\tilde{\chi }}_1^0{\tilde{\chi }}_2^0$	OS1	17.87	0.08	0.14	$H_5\to {\tilde{\chi }}_1^0{\tilde{\chi }}_2^0$	9.29	0.24	0.07
	$\overline{\text{DR}}$	14.59	-11.68	0.11		7.65	1.94	0.06
$H_4\to {\tilde{\chi }}_1^0{\tilde{\chi }}_3^0$	OS1	395.9	-5.5	4.59	$H_5\to {\tilde{\chi }}_1^0{\tilde{\chi }}_3^0$	423.04	0.02	3.3
	$\overline{\text{DR}}$	451.86	6.59	3.5		448.9	7.83	3.5
$H_4\to {\tilde{\chi }}_1^0{\tilde{\chi }}_4^0$	OS1	1448.21	-0.03	11.22	$H_5\to {\tilde{\chi }}_1^0{\tilde{\chi }}_4^0$	1161.18	-0.02	9.05
	$\overline{\text{DR}}$	1438.6	-3.68	11.15		1157.98	-2.81	9.03
$H_4\to {\tilde{\chi }}_1^0{\tilde{\chi }}_5^0$	OS1	194.66	-0.02	1.51	$H_5\to {\tilde{\chi }}_1^0{\tilde{\chi }}_5^0$	338.35	-0.02	2.64
	$\overline{\text{DR}}$	195.02	-1.63	1.51		341.12	-1.58	2.66
$H_4\to {\tilde{\chi }}_2^0{\tilde{\chi }}_2^0$	OS1	27.86	-0.13	0.22	$H_5\to {\tilde{\chi }}_2^0{\tilde{\chi }}_2^0$	34.6	-0.1	0.27
	$\overline{\text{DR}}$	25.54	-20.08	0.2		32.72	-15.35	0.26
$H_4\to {\tilde{\chi }}_2^0{\tilde{\chi }}_3^0$	OS1	425.3	0.02	3.29	$H_5\to {\tilde{\chi }}_2^0{\tilde{\chi }}_3^0$	413.74	0.01	3.22
	$\overline{\text{DR}}$	450.13	7.77	3.49		438.11	6.48	3.42
$H_4\to {\tilde{\chi }}_2^0{\tilde{\chi }}_4^0$	OS1	1125.46	-0.02	8.72	$H_5\to {\tilde{\chi }}_2^0{\tilde{\chi }}_4^0$	1379.31	-0.03	10.75
	$\overline{\text{DR}}$	1116.5	-3.34	8.66		1366.13	-3.91	10.66
$H_4\to {\tilde{\chi }}_2^0{\tilde{\chi }}_5^0$	OS1	382.4	-0.02	2.96	$H_5\to {\tilde{\chi }}_2^0{\tilde{\chi }}_5^0$	180.21	-0.02	1.4
	$\overline{\text{DR}}$	381.89	-2.35	2.96		180.34	-1.86	1.41
$H_4\to {\tilde{\chi }}_3^0{\tilde{\chi }}_3^0$	OS1	1.2	0.19	0.01	$H_5\to {\tilde{\chi }}_3^0{\tilde{\chi }}_3^0$	1.19	0.18	0.01
	$\overline{\text{DR}}$	1.29	28.87	0.01		1.29	27.81	0.01
$H_4\to {\tilde{\chi }}_3^0{\tilde{\chi }}_4^0$	OS1	3.18	0.05	0.02	$H_5\to {\tilde{\chi }}_3^0{\tilde{\chi }}_4^0$	1.48	0.06	0.01
	$\overline{\text{DR}}$	2.95	-2.79	0.02		1.43	0.89	0.01
$H_4\to {\tilde{\chi }}_1^{+}{\tilde{\chi }}_1^{-}$	OS1	393.59	-0.01	3.05	$H_5\to {\tilde{\chi }}_1^{+}{\tilde{\chi }}_1^{-}$	257.68	-0.03	2.01
	$\overline{\text{DR}}$	343.95	-13.61	2.67		215.16	-19.09	1.68
$H_4\to {\tilde{\chi }}_1^{+}{\tilde{\chi }}_2^{-}$	OS1	2440.51	-0.01	18.91	$H_5\to {\tilde{\chi }}_1^{+}{\tilde{\chi }}_2^{-}$	2482.01	-0.01	19.34
	$\overline{\text{DR}}$	2461.41	-0.61	19.0		2498.8	-0.9	19.49
$H_4\to {\tilde{\chi }}_2^{+}{\tilde{\chi }}_1^{-}$	OS1	2440.51	-0.01	18.91	$H_5\to {\tilde{\chi }}_2^{+}{\tilde{\chi }}_1^{-}$	2482.01	-0.01	19.34
	$\overline{\text{DR}}$	2461.41	-0.61	19.0		2498.8	-0.9	19.49

Table 5. Table 5: Scenario2 : The tree-level and one-loop corrected stop and sbottom masses in the DR ¯ ¯ DR \overline{\text{DR}} and OS schemes.

		$m_{{\tilde{t}}_1}\left[\text{GeV}\right]$	$m_{{\tilde{t}}_2}\left[\text{GeV}\right]$	$m_{{\tilde{b}}_1}\left[\text{GeV}\right]$	$m_{{\tilde{b}}_2}\left[\text{GeV}\right]$
OS	tree	334	1166	1122	2297
	1loop	334	1166	1125	2297
$\overline{\text{DR}}$	tree	585	1216	1181	3000
	1loop	341	1152	1152	2294

Table 6. Table 6: Scenario2 : Mass values in GeV and main components of the neutral Higgs bosons at two-loop order 𝒪 ​ ( α t ​ α s + α t 2 ) 𝒪 subscript 𝛼 𝑡 subscript 𝛼 𝑠 superscript subscript 𝛼 𝑡 2 {\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2}) obtained by using the DR ¯ ¯ DR \overline{\text{DR}} renormalization in the top/stop sector.

	$H_1$	$H_2$	$H_3$	$H_4$	$H_5$
two-loop $𝒪\left({\alpha }_t{\alpha }_s+{\alpha }_t^2\right)$	123.63	1621.65	1865.39	1895.83	2538.29
main component	$h_u$	$a_s$	$h_d$	$a$	$h_s$

Equations

$$W_{\text{NMSSM}}=W_{\text{MSSM}}-{ϵ}_{ab}\lambda \hat{S}{\hat{H}}_d^a{\hat{H}}_u^b+\frac{1}{3}\kappa {\hat{S}}^3,$$

$$\displaystyle W_{\text{MSSM}}=-\epsilon_{ab}\bigl{(}y_{u}\hat{H}_{u}^{a}\hat{Q}^{b}\hat{U}^{c}-y_{d}\hat{H}_{d}^{a}\hat{Q}^{b}\hat{D}^{c}-y_{e}\hat{H}_{d}^{a}\hat{L}^{b}\hat{E}^{c}\bigr{)}\;,$$

$${\mathcal{L}}_{\text{NMSSM}}^{\text{soft}}={\mathcal{L}}_{\text{MSSM}}^{\text{soft}}-m_S^2|S{|}^2+({ϵ}_{ab}{A}_{\lambda }\lambda S{H}_d^a{H}_u^b-\frac{1}{3}{A}_{\kappa }\kappa S^3+h.c.).$$

$${\mathcal{L}}_{\text{MSSM}}^{\text{soft}}$$

$$+{ϵ}_{ab}(y_u{A}_u{H}_u^a{\tilde{Q}}^b{\tilde{u}}_R^{\ast }-y_d{A}_d{H}_d^a{\tilde{Q}}^b{\tilde{d}}_R^{\ast }-y_e{A}_e{H}_d^a{\tilde{Q}}^b{\tilde{e}}_R^{\ast }+h.c.)$$

$$-\frac{1}{2}(M_1\tilde{B}\tilde{B}+M_2{\tilde{W}}_i{\tilde{W}}_i+M_3\tilde{G}\tilde{G}+h.c),$$

$$H_d=\left(\begin{array}{c}\frac{1}{\sqrt{2}}(v_d+h_d+i{a}_d)\\ h_d^{-}\end{array}\right),H_u=e^{i{\phi }_u}\left(\begin{array}{c}{h}_u^{+}\\ \frac{1}{\sqrt{2}}(v_u+h_u+i{a}_u)\end{array}\right),S=\frac{e^{i{\phi }_s}}{\sqrt{2}}(v_s+h_s+i{a}_s).$$

$${\mathcal{L}}_{\text{neutral}}^m=\frac{1}{2}{\varphi }^T{M}_{\varphi \varphi }\varphi .$$

$$\text{diag}(m_{h_1}^2,m_{h_2}^2,m_{h_3}^2,m_{h_4}^2,m_{h_5}^2,0)$$

$$(h_1,h_2,h_3,h_4,h_5,G{)}^T$$

$${\mathcal{L}}_{\text{charged}}^m=\left(\begin{array}{c}{h}_d^{+},h_u^{+}\end{array}\right)M_{h^{+}{h}^{+}}\left(\begin{array}{c}{h}_d^{-}\\ h_u^{-}\end{array}\right),$$

$$\displaystyle M_{h^{+}h^{+}}=\frac{1}{2}\begin{pmatrix}{t_{\beta}}&1\\ 1&1/{t_{\beta}}\end{pmatrix}\bigg{[}$$

$$\displaystyle-\frac{2|\lambda|^{2}M_{W}^{2}s_{\theta_{W}}^{2}}{e^{2}}s_{2\beta}\bigg{]}\;,$$

$$\tan \beta =v_u/v_d.$$

$$M_{H^{\pm }}^2=M_W^2+\frac{|\lambda |v_s}{s_{2\beta }\cos ({\phi }_{\lambda }+{\phi }_u+{\phi }_s)}\left(\sqrt{2}\text{Re}{A}_{\lambda }+|\kappa |v_s\cos ({\phi }_{\kappa }+3{\phi }_s)\right)-\frac{2|\lambda {|}^2{M}_W^2{s}_{{\theta }_W}^2}{e^2}.$$

$$M_N=$$

$$\left(\begin{array}{ccccc}{M}_1& 0& -c_{\beta }{M}_Z{s}_{{\theta }_W}& M_Z{s}_{\beta }{s}_{{\theta }_W}{e}^{-i{\phi }_u}& 0\\ 0& M_2& c_{\beta }{M}_W& -M_W{s}_{\beta }{e}^{-i{\phi }_u}& 0\\ -c_{\beta }{M}_Z{s}_{{\theta }_W}& c_{\beta }{M}_W& 0& -\lambda \frac{v_s}{\sqrt{2}}{e}^{i{\phi }_s}& -\frac{\sqrt{2}{M}_W{s}_{\beta }{s}_{{\theta }_W}\lambda e^{i{\phi }_u}}{e}\\ M_Z{s}_{\beta }{s}_{{\theta }_W}{e}^{-i{\phi }_u}& -M_W{s}_{\beta }{e}^{-i{\phi }_u}& -\lambda \frac{v_s}{\sqrt{2}}{e}^{i{\phi }_s}& 0& -\frac{\sqrt{2}{M}_W{c}_{\beta }{s}_{{\theta }_W}\lambda }{e}\\ 0& 0& -\frac{\sqrt{2}{M}_W{s}_{\beta }{s}_{{\theta }_W}\lambda e^{i{\phi }_u}}{e}& -\frac{\sqrt{2}{M}_W{c}_{\beta }{s}_{{\theta }_W}\lambda }{e}& \sqrt{2}\kappa v_s{e}^{i{\phi }_s}\end{array}\right)$$

$$\displaystyle\tilde{\chi}_{i}^{0}=\left(\begin{array}[]{c}\chi_{i}^{0}\vskip 3.0pt plus 1.0pt minus 1.0pt\\ \overline{\chi_{i}^{0}}\end{array}\right)\quad\ \text{with}\quad\ \chi^{0}_{i}=N_{ij}\psi^{0}_{j},\quad\ i,\,j=1,\dots,5~{},$$

$$\overline{{\chi }_i^0}=i{\sigma }_2{\chi }_i^{0\ast }.$$

$${\psi }_R^{-}=\left(\begin{array}{c}{\tilde{W}}^{-}\\ {\tilde{H}}_d^{-}\end{array}\right)\text{mbox}and{\psi }_L^{+}=\left(\begin{array}{c}{\tilde{W}}^{+}\\ {\tilde{H}}_u^{+}\end{array}\right)$$

$$\mathcal{L}=({\psi }_R^{-}{)}^T{M}_C{\psi }_L^{+}+h.c.,$$

$$M_C=\left(\begin{array}{cc}{M}_2& \sqrt{2}{s}_{\beta }{M}_W{e}^{-i{\phi }_u}\\ \sqrt{2}{c}_{\beta }{M}_W& \lambda \frac{v_s}{\sqrt{2}}{e}^{i{\phi }_s}\end{array}\right).$$

$$\text{diag}(m_{{\tilde{\chi }}_1^{\pm }},m_{{\tilde{\chi }}_2^{\pm }})=U^{\ast }{M}_C{V}^{†},$$

$${\tilde{\chi }}_L^{+}=V{\psi }_L^{+}\text{mbox}and{\tilde{\chi }}_R^{-}=U{\psi }_R^{-},$$

$$\displaystyle\tilde{\chi}^{+}_{i}=\left(\begin{array}[]{c}\tilde{\chi}^{+}_{L_{i}}\vskip 3.0pt plus 1.0pt minus 1.0pt\\ \overline{\tilde{\chi}^{-}_{R_{i}}}\end{array}\right)$$

$$\mathcal{L}=$$

$$+\overline{{\tilde{\chi }}_k^0}/p{P}_L{\tilde{\chi }}_k^0+\overline{{\tilde{\chi }}_k^0}/p{P}_R{\tilde{\chi }}_k^0-\overline{{\tilde{\chi }}_k^0}{\left[N^{\ast }{M}_N{N}^{†}\right]}_{kl}{P}_L{\tilde{\chi }}_l^0-\overline{{\tilde{\chi }}_k^0}{\left[N{M}_N^{†}{N}^T\right]}_{kl}{P}_R{\tilde{\chi }}_l^0,$$

$$M_{\tilde{t}}=\left(\begin{array}{cc}{m}_{{\tilde{Q}}_3}^2+m_t^2+M_Z^2{c}_{2\beta }(\frac{1}{2}-\frac{2}{3}{s}_{{\theta }_W}^2)& m_t\left(A_t^{\ast }{e}^{-i{\phi }_u}-{\mu }_{\text{eff}}/t_{\beta }\right)\\ m_t\left(A_t{e}^{i{\phi }_u}-{\mu }_{\text{eff}}^{\ast }/t_{\beta }\right)& m_t^2+m_{{\tilde{t}}_R}^2+\frac{2}{3}{M}_Z^2{c}_{2\beta }{s}_{{\theta }_W}^2\end{array}\right),$$

$$M_{\tilde{b}}=\left(\begin{array}{cc}{m}_{{\tilde{Q}}_3}^2+m_b^2+M_Z^2{c}_{2\beta }(-\frac{1}{2}+\frac{1}{3}{s}_{{\theta }_W}^2)& m_b\left(A_b^{\ast }-e^{i{\phi }_u}{\mu }_{\text{eff}}{t}_{\beta }\right)\\ m_b\left(A_b-e^{-i{\phi }_u}{\mu }_{\text{eff}}^{\ast }{t}_{\beta }\right)& m_b^2+m_{{\tilde{b}}_R}^2-\frac{1}{3}{M}_Z^2{c}_{2\beta }{s}_{{\theta }_W}^2\end{array}\right),$$

$${\mu }_{\text{eff}}=\frac{\lambda v_s{e}^{i{\phi }_s}}{\sqrt{2}}.$$

$$\text{diag}(m_{{\tilde{q}}_1}^2,m_{{\tilde{q}}_2}^2)=U^{\tilde{q}}{M}_{\tilde{q}}{U}^{\tilde{q}†},\left(\begin{array}{c}{\tilde{q}}_1\\ {\tilde{q}}_2\end{array}\right)=U^{\tilde{q}}\left(\begin{array}{c}{\tilde{q}}_L\\ {\tilde{q}}_R\end{array}\right),q=t,b,$$

$$m_{H_d}^2,m_{H_u}^2,m_S^2,M_W^2,M_Z^2,e,\tan \beta ,v_s,{\phi }_s,{\phi }_u,|\lambda |,{\phi }_{\lambda },|\kappa |,{\phi }_{\kappa },\text{Re}{A}_{\lambda },\text{Im}{A}_{\lambda },\text{Re}{A}_{\kappa },\text{Im}{A}_{\kappa }.$$

$$t_{\varphi }\equiv {\frac{\partial V_{\text{Higgs}}}{\partial \varphi }|}_{\text{Min.}}\stackrel{!}{=}0.$$

Full text

1 Introduction

The discovery of a scalar particle by the LHC experiments ATLAS [1] and CMS [2] and the subsequent investigation of its properties revealed a Higgs boson that behaves very Standard Model (SM)-like. Also years after its discovery there are no evidences for new physics from direct searches. In this situation the precise investigation of the Higgs sector plays an important role. Indirect effects of physics beyond the SM (BSM) might show up in the properties of the discovered Higgs boson. With a mass of 125.09 GeV [3] it does not exclude the possibility for the Higgs boson of a supersymmetric (SUSY) extension of the SM, like the minimal (MSSM) or the next-to-minimal (NMSSM) ones. Supersymmetry certainly belongs to the best motivated and most intensively studied BSM extensions, and the NMSSM, with a Higgs sector consisting of seven Higgs bosons arising after electroweak symmetry breaking (EWSB) from the two doublet and singlet fields of the Higgs sector, provides a rich phenomenology [4, 5]. The experimental limits strongly restrict possible new physics effects in the Higgs sector and call for precision in the theory predictions for the Higgs boson observables. This is also necessary in order to be able to distinguish between new physics extensions in case of discovery.

In this paper, we concentrate on the NMSSM Higgs boson decays. While the (SUSY-)QCD corrections can be taken over from the MSSM case with the appropriate modifications and a minimum of effort, this is not the case for the electroweak (EW) corrections. In the recent years there has been some progress on this subject. In the CP-conserving NMSSM, members of our group computed the next-to-leading order (NLO) SUSY-EW and SUSY-QCD corrections to the decays of CP-odd NMSSM Higgs bosons into stop pairs and found that the both the EW and the SUSY-QCD corrections are significant and can be of opposite sign [6]. The authors of [7, 8] provided in the framework of the CP-conserving NMSSM its full one-loop renormalization and the two-body Higgs decays at one-loop order in the on-shell (OS) renormalization scheme. A generic calculation of the two-body partial decays widths at full one-loop level was provided in [9] in the $\overline{\text{DR}}$ scheme. In [10] the full one-loop corrections for the neutral CP-violating NMSSM Higgs bosons were calculated to their decays into fermions and gauge bosons and combined with the leading QCD corrections. For the Higgs-to-Higgs decays, we provided in previous papers the complete one-loop [11] and the order ${\cal O}(\alpha_{t}\alpha_{s})$ two-loop [12] corrections in the CP-conserving and CP-violating NMSSM, respectively.

In this work, we compute, in the framework of the CP-violating NMSSM, the complete next-to-leading order (SUSY-)electroweak corrections to the neutral NMSSM Higgs boson decays into all tree-level induced SM final states, i.e. into fermion and massive gauge boson pairs, but also into non-SM pairs, namely gauge and Higgs boson final states, chargino and neutralino pairs, and into squarks. Where applicable we combine our corrections with the already available (SUSY-)QCD corrections. We furthermore include the complete one-loop corrections to the decays into Higgs boson pairs, cf. Refs. [11, 12]. For the loop-induced decays into gluon and photon pairs as well as a photon and a $Z$ boson no corrections are provided as they would be of two-loop order. For the first time, we present the one-loop corrections to the electroweakino, stop and sbottom masses in the context of the CP-violating NMSSM, by applying both OS and $\overline{\text{DR}}$ schemes. We have implemented our corrections in our original code NMSSMCALC [13], which calculates, based on a mixed OS-$\overline{\text{DR}}$ scheme, the NMSSM Higgs mass corrections and decays in both the CP-conserving and CP-violating case. This way we provide the NMSSM Higgs boson decays and branching ratios at presently highest possible precision including the state-of-the-art (SUSY-)QCD and the computed (SUSY-)EW corrections. In the EW higher-order corrections we not only include the NLO vertex corrections but also take into account the proper on-shell conditions of the external decaying Higgs bosons up to two-loop order ${\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2})$. This is the order up to which the mass corrections for the NMSSM Higgs bosons both in the CP-conserving [14] and CP-violating case [15, 16, 17] have been implemented in NMSSMCALC. The new program is called NMSSMCALCEW can be obtained at the url: https://www.itp.kit.edu/$\sim$maggie/NMSSMCALCEW/. Here also a detailed description of the program and its structure are given, instructions on how to compile and run it as well as information on modifications, which is constantly updated. A brief description of the code is given in Appendix B.

The paper is organized as follows. In Section 2 we introduce the NMSSM sectors at tree level, that are relevant for our computation, and set our notation before we move on to the NMSSM at one-loop level in Section 3. We here describe the renormalization of the Higgs, chargino/neutralino, and squark sectors as well as the loop corrections to the Higgs boson masses and mixings, to the neutralino and chargino masses, and finally to the squark masses and their mixings. Section 4 is devoted to the detailed presentation of our calculation of the one-loop corrections to the neutral non-loop induced Higgs boson decays into on-shell final states, namely into fermion pairs, massive gauge boson pairs, final states with one gauge and one Higgs boson, neutralino and chargino pairs, and squark final states. In Section 5 we present the numerical analysis of the one-loop corrections to the Higgs boson branching ratios into SM and SUSY final states, where we discuss in particular the size of the newly implemented corrections to both the branching ratios and to the electroweakino and third generation squark masses. Our conclusions are given in Section 6. Explicit expressions of the counterterm couplings for the decays of neutral Higgs bosons into a squark pair are displayed in Appendix A.

2 The NMSSM at Tree Level

We are working in the complex NMSSM with a preserved $\mathbb{Z}_{3}$ symmetry. The Lagrangian of the NMSSM can be divided into the Lagrangian of the MSSM and the additional part coming from the NMSSM. For convenience of the reader and to set our notation, we give here the parts of the Lagrangian that are relevant for our calculations. For the Higgs sector we need the NMSSM Higgs potential. It is derived from the NMSSM superpotential $W_{\text{NMSSM}}$, the corresponding soft SUSY-breaking terms, and the $D$-term contributions. With the Higgs doublet superfields $\hat{H}_{u}$ and $\hat{H}_{d}$ coupling to the up- and down-type quark superfields, respectively, and the singlet superfield $\hat{S}$, we have for the NMSSM superpotential

[TABLE]

where $a,b=1,2$ are the indices of the fundamental $SU(2)_{L}$ representation and $\epsilon_{ab}$ is the totally antisymmetric tensor with $\epsilon_{12}=\epsilon^{12}=1$ . The MSSM part reads

[TABLE]

in terms of the left-handed quark and lepton superfield doublets $\hat{Q}$ and $\hat{L}$ and the right-handed up-type, down-type, and electron-type superfield singlets $\hat{U}$, $\hat{D}$, and $\hat{E}$, respectively. Charge conjugation is denoted by the superscript $c$, and color and generation indices have been omitted. The NMSSM superpotential contains the coupling $\kappa$ of the self-interaction of the new singlet superfield and the coupling $\lambda$ for the $\hat{S}$ interaction with the two Higgs doublet superfields. Both couplings are complex. The quark and lepton Yukawa couplings $y_{d},$ $y_{u}$, and $y_{e}$ are in general complex. However, in case of no generation mixing, as assumed in this paper, the phases of the Yukawa couplings can be absorbed through a redefinition of the quark fields, so that the phases can be chosen arbitrarily without changing the physical meaning [18]. The soft SUSY-breaking NMSSM Lagrangian in terms of the component fields $H_{u},H_{d}$ and $S$ reads

[TABLE]

It contains two more complex parameters specific to the NMSSM, the soft SUSY-breaking trilinear couplings $A_{\lambda}$ and $A_{\kappa}$. The soft SUSY-breaking MSSM contribution can be cast into the form

[TABLE]

where the SM-type and SUSY fields corresponding to a superfield (denoted with a hat) are represented by a letter without and with a tilde, respectively. The indices $\tilde{Q}$ ($\tilde{L}$) of the soft SUSY-breaking masses denote, exemplary for the first generation, the left-handed quark (lepton) doublet component fields of the corresponding quark and lepton superfields, and $\tilde{u}_{R},\tilde{d}_{R},\tilde{e}_{R}$ the right-handed component fields for the up-type and down-type quarks, and charged leptons, respectively,. The trilinear couplings $A_{u}$, $A_{d}$ and $A_{e}$ of the up-type and down-type quarks and charged leptons are in general complex, whereas the soft SUSY-breaking mass terms $m_{x}^{2}$ ($x=S,H_{u},H_{d},\tilde{Q},\tilde{u}_{R},\tilde{d}_{R},\tilde{L},\tilde{e}_{R}$) are real. The soft SUSY-breaking mass parameters of the gauginos, $M_{1}$, $M_{2}$, $M_{3}$, for the bino, the winos, and the gluinos, $\tilde{B}$, $\tilde{W}_{i}$ ($i=1,2,3$), and $\tilde{G}$, corresponding to the weak hypercharge $U(1)$, the weak isospin $SU(2)$, and the color $SU(3)$ symmetry, are in general complex. The $R$-symmetry can be exploited to choose either $M_{1}$ or $M_{2}$ to be real. In this paper we keep both $M_{1}$ and $M_{2}$ complex.

Expanding the scalar Higgs fields about their vacuum expectation values (VEVs) $v_{u}$, $v_{d}$, and $v_{s}$, two further phases, $\varphi_{u}$ and $\varphi_{s}$, are introduced which describe the phase differences between the VEVs,

[TABLE]

For vanishing phases, the fields $h_{i}$ and $a_{i}$ with $i=d,u,s$ correspond to the CP-even and CP-odd part of the neutral entries of $H_{d}$, $H_{u}$ and $S$. The charged components are denoted by $h_{d,u}^{\pm}$. In this paper, we set the phases of the Yukawa couplings to zero. We furthermore re-phase the left- and right-handed up-quark fields as $u_{L}\to e^{-i\varphi_{u}}u_{L}$ and $u_{R}\to e^{i\varphi_{u}}u_{R}$, so that the quark and lepton mass terms yield real masses.

After electroweak symmetry breaking (EWSB) the six Higgs interaction states mix and in the basis $\phi=(h_{d},h_{u},h_{s},a_{d},a_{u},a_{s})$ the mass term is given by

[TABLE]

The mass matrix $M_{\phi\phi}$ is obtained from the second derivative of the Higgs potential with respect to the Higgs fields in the vacuum. The explicit expression of the mass matrix $M_{\phi\phi}$ can be found in Ref. [15]. The transformation into mass eigenstates at tree level can be performed with orthogonal matrices $\mathcal{R},\mathcal{R}^{G}$,

[TABLE]

where the matrix $\mathcal{R}^{G}$ is used first to single out the Goldstone boson $G$. The tree-level Higgs mass eigenstates are denoted by the small letter $h$, and their masses are ordered as $m_{h_{1}}\leq...\leq m_{h_{5}}$.

The mass term for the charged components of the Higgs doublets in

[TABLE]

is given by

[TABLE]

where $M_{W}$ is the mass of the $W$ boson, $\theta_{W}$ the electroweak mixing angle, $e$ the electric charge and $\varphi_{\lambda},$ $\varphi_{\kappa}$ the complex phases of $\lambda$ and $\kappa$, respectively. The angle $\beta$ is defined as

[TABLE]

Here and in the following we use the short hand notation $c_{x}=\cos x$, $s_{x}=\sin x$ and $t_{x}=\tan x$. Diagonalizing this mass matrix by a rotation matrix with the angle $\beta_{c}$, for which at Born level $\beta_{c}=\beta$, one obtains the charged Higgs boson mass as

[TABLE]

The charged Goldstone boson $G^{\pm}$, on the other hand, is massless.

The fermionic superpartners of the neutral Higgs bosons, $\tilde{H}_{d}^{0}$, $\tilde{H}_{u}^{0}$, $\tilde{S}$, and of the neutral gauge bosons, $\tilde{B}$, $\tilde{W}_{3}$, mix, and in the Weyl spinor basis $\psi^{0}=(\tilde{B},\tilde{W}_{3},\tilde{H}^{0}_{d},\tilde{H}^{0}_{u},\tilde{S})^{T}$ the neutralino mass matrix $M_{N}$ is given by

[TABLE]

after EWSB, where $M_{Z}$ is the $Z$ boson mass. The symmetric neutralino mass matrix can be diagonalized by a $5\times 5$ matrix $N$, yielding $\text{diag}(m_{\tilde{\chi}^{0}_{1}},m_{\tilde{\chi}^{0}_{2}},m_{\tilde{\chi}^{0}_{3}},m_{\tilde{\chi}^{0}_{4}},m_{\tilde{\chi}^{0}_{5}})=N^{*}M_{N}N^{\dagger}$, where the absolute mass values are ordered as $|m_{\tilde{\chi}^{0}_{1}}|\leq...\leq|m_{\tilde{\chi}^{0}_{5}}|$. The neutralino mass eigenstates $\tilde{\chi}^{0}_{i}$, expressed as Majorana spinors, can then be obtained by

[TABLE]

where, in terms of the Pauli matrix $\sigma_{2}$,

[TABLE]

The fermionic superpartners of the charged Higgs and gauge bosons are given in terms of the Weyl spinors $\tilde{H}_{d}^{\pm}$, $\tilde{H}_{u}^{\pm}$, $\tilde{W}^{-}$ and $\tilde{W}^{+}$. With

[TABLE]

the mass term for these spinors is of the form

[TABLE]

where

[TABLE]

The chargino mass matrix $M_{C}$ can be diagonalized with the help of two unitary $2\times 2$ matrices, $U$ and $V$, yielding

[TABLE]

with $m_{\tilde{\chi}^{\pm}_{1}}\leq m_{\tilde{\chi}^{\pm}_{2}}$. The left-handed and the right-handed part of the mass eigenstates are

[TABLE]

respectively, with the mass eigenstates ($i=1,2$)

[TABLE]

written as Dirac spinors. In summary, the bilinear terms in the chargino and neutralino mass eigenstates are given by

[TABLE]

where the left- and right-handed projectors are defined as $P_{L/R}=\left(1\mp\gamma_{5}\right)/2$ and $i,j=1,2$ and $k,l=1,\ldots,5$.

The scalar partners of the left- and right-handed quarks are denoted by $\tilde{q}_{L}$ and $\tilde{q}_{R}$, respectively. The mixing matrix for the top squark is given by

[TABLE]

while the bottom squark matrix reads

[TABLE]

where

[TABLE]

The mass eigenstates are obtained by diagonalizing these squark matrices with the unitary transformations

[TABLE]

with the usual convention $m_{\tilde{q}_{1}}\leq m_{\tilde{q}_{2}}$.

3 The NMSSM at One-Loop Level

3.1 Renormalization

3.1.1 The Higgs Sector

For the Higgs sector we follow the mixed on-shell OS-$\overline{\text{DR}}$ renormalization scheme described and applied in Refs. [14, 15, 16, 17]. We do not repeat all details here but quote the most important formulae. There are eighteen parameters entering the Higgs sector at tree level,

[TABLE]

Note that for the sake of convenience we decompose the complex trilinear couplings $A_{\lambda}$ and $A_{\kappa}$ into a real and imaginary part in contrast to Ref. [15] where the absolute values and complex phases were used. It was found in Ref. [15] that the four complex phases $\varphi_{s},$ $\varphi_{u},$ $\varphi_{\lambda}$ and $\varphi_{\kappa}$ do not need to be renormalized. We verified this statement and will discard them in our renormalization procedure.

In our introduction of the NMSSM Higgs sector in Sec. 2 we have already replaced the $U(1)$ and $SU(2)$ gauge couplings $g^{\prime}$ and $g$ and the VEV $v$ by the physical observables $M_{W}$, $M_{Z}$ and $e$. It is convenient to further convert, where possible, the input parameters of the set Eq. (3.31) to parameters that can be interpreted more easily in terms of physical quantities. Thus we trade the three soft SUSY-breaking mass parameters $m_{H_{d}}^{2},m_{H_{u}}^{2},m_{S}^{2}$ as well as $\text{Im}\,A_{\lambda}$ and $\text{Im}\,A_{\kappa}$ for the tadpole parameters $t_{\phi}$ ($\phi=h_{d},h_{u},h_{s},a_{d},a_{s}$). These coefficients of the terms of the Higgs potential $V_{\text{Higgs}}$ are linear in the Higgs boson fields and have to vanish, in order to ensure the minimum at non-vanishing VEVs $v_{u},$ $v_{d}$, $v_{s}$,

[TABLE]

It is debatable whether the tadpole parameters can be called physical quantities, but certainly their introduction is motivated by physical interpretation. In the same way, in a slight abuse of the language, we will call the renormalization conditions for the tadpole parameters on-shell. With the new set of input parameters, we allow for two possible renormalization schemes in our Higgs mass calculation. The difference between the two schemes relates to the treatment of the charged Higgs mass. In the first scheme the charged Higgs mass is an OS input parameter,

[TABLE]

while in the second scheme $\text{Re}\,A_{\lambda}$ is an input parameter renormalized as a $\overline{\text{DR}}$ parameter and the charged Higgs mass is a derived quantity,

[TABLE]

For the definition of the one-loop OS counterterms we refer the reader to Ref. [15]. The counterterms of the $\overline{\text{DR}}$ parameters do not contribute to the final physical results but they are kept to ensure UV finiteness. We define these counterterms solely in the Higgs sector by requiring that all renormalized self-energies of the Higgs bosons be finite. This is different with respect to the definition in Ref. [15] where the chargino and neutralino sectors were used. The numerical results between the two definitions are identical, however. We renormalize the Higgs fields in the $\overline{\text{DR}}$ scheme as described in Ref. [15] at one-loop level, and in Refs. [16, 17] at two-loop order ${\cal O}(\alpha_{t}\alpha_{s})$ and ${\cal O}(\alpha_{t}^{2})$, respectively.

3.1.2 The Chargino and Neutralino Sector

The chargino and neutralino sectors are described by fourteen real parameters: $M_{W},M_{Z},$ $\tan\beta,v_{s},\varphi_{s},\varphi_{u},\lvert\lambda\rvert,$ $\varphi_{\lambda},\lvert\kappa\rvert,\varphi_{\kappa},|M_{1}|,\varphi_{M_{1}},|M_{2}|,\varphi_{M_{2}}$. Since the first ten of these already appear in the Higgs sector, there remains to define the renormalization conditions for the four parameters $|M_{1}|,\varphi_{M_{1}},|M_{2}|,\varphi_{M_{2}}$.444In the case of the MSSM one needs to renormalize three complex parameters $M_{1},M_{2},\mu$. There are no physical renormalization conditions to fix the counterterms of the phases $\varphi_{M_{1}},\varphi_{M_{2}}$. It has been found in Ref. [15] that the complex phases of $M_{1}$ and $M_{2}$ do not need to be renormalized. We verified this statement in our computation.555The same holds true in the complex MSSM [19]. The three complex phases of $M_{1}$, $M_{2}$ and $\mu$ do not need to be renormalized in order to render all Green functions finite. In addition, we have to renormalize the chargino and neutralino fields in order to obtain finite self-energies. In the literature there exist two descriptions for the introduction of wave function renormalization (WFR) constants. In the Espriu-Manzano-Talavera (EMT) description two independent renormalization constants were introduced for incoming and outgoing fermions [20, 21, 19]. Thanks to more degrees of freedom one can keep contributions arising from absorptive parts of the loop integrals and eliminate completely the mixing self-energies thereby fulfilling the standard OS conditions. However, the hermicity of the renormalized Lagrangian is not satisfied anymore. In the Denner description [22], one WFR constant was used instead. It preserves the hermicity constraint, but the absorptive part of the loop integral must be eliminated. We want to investigate the effect of the absorptive part and therefore apply both descriptions. In the following we will derive the formulae in the EMT method. From these formulae, one can easily obtain the ones in the Denner description. The bare parameters and fields are replaced by the renormalized ones and the corresponding counterterms as

[TABLE]

where $i,j=1,2$ and $k,l=1,\ldots,5$. Since the neutralinos are Majorana fermions we have

[TABLE]

Note that we do not need to renormalize the rotation matrices $U,$ $V$ and $N$ because their counterterms are redundant. They always appear in combination with WFR constants. One can therefore always redefine the WFR constants to absorb the counterterms of the rotation matrices, as shown in Ref. [23]. In general, the renormalized self-energies $\hat{\Sigma}$ of the fermions can be cast into the following form, cf. Ref. [22],

[TABLE]

with

[TABLE]

where $\Sigma_{ij}(p)$ denotes the unrenormalized self-energy of the transition $\tilde{\chi}^{+}_{i}\to\tilde{\chi}^{+}_{j},~{}i,j=1,2$, for the charginos and $\tilde{\chi}^{0}_{i}\to\tilde{\chi}^{0}_{j},~{}i,j=1,\cdots,5$, for the neutralinos. For the charginos, the tree-level mass matrix $M^{\text{tree}}_{\tilde{\chi}}$ and its counterterm $\delta M^{\text{tree}}_{\tilde{\chi}}$ are given by

[TABLE]

and for the neutralinos by

[TABLE]

In the following, we will discuss the OS conditions for the general fermion fields $\tilde{\chi}_{i}$ having the tree-level masses $m_{\tilde{\chi}_{i}}$. The renormalized fermion propagator matrix is given by

[TABLE]

where the renormalized one-particle irreducible (1PI) two-point functions $\hat{\Gamma}$ are related to the renormalized self-energies as

[TABLE]

The propagator matrix has complex poles at ${\cal M}^{2}_{\tilde{\chi}_{i}}=M_{\tilde{\chi}_{i}}^{2}-iM_{\tilde{\chi}_{i}}\Gamma_{\tilde{\chi}_{i}}$, where $\Gamma_{\tilde{\chi}_{i}}$ denotes the decay widths. In the OS scheme we require that the tree-level masses are equal to the physical masses, the mixing terms are vanishing at the poles and that the residues of the propagators are unity666In the Denner description, $\hat{\Sigma}_{ii}$ in Eq. (3.51), $\hat{\Sigma}_{ij}$ in Eq. (3.52) and $\hat{\Gamma}_{ii}$ in Eq. (3.53) are replaced by $\widetilde{\text{Re}}\hat{\Sigma}_{ii}$, $\widetilde{\text{Re}}\hat{\Sigma}_{ij}$ and $\widetilde{\text{Re}}\hat{\Gamma}_{ii}$, respectively. $\widetilde{\text{Re}}$ means that one takes only the real part of the loop integrals but leaves the couplings unaffected.,

[TABLE]

In addition, we require the chiral structure to vanish in the OS limit777If we use the Denner OS conditions in Ref. [22] then these relations are automatically satisfied for the real parts only.,

[TABLE]

Applying the decomposition in Eq. (3.42) and the tree-level relations

[TABLE]

one obtains the mass counterterms $\delta m_{\tilde{\chi}_{i}}=\mbox{Re}\,(\delta M_{\tilde{\chi}}^{\text{tree}})_{ii}$, with

[TABLE]

the off-diagonal wave function renormalization constants, $\delta Z_{L/R,ij},\delta\bar{Z}_{L/R,ij}$,888In the Denner description $\Sigma_{ij}^{L/R/Ls/Rs}$ are replaced by $\widetilde{\text{Re}}\Sigma_{ij}^{L/R/Ls/Rs}$.,999Note that the $\delta\bar{Z}$ are obtained from $\delta Z$ by interchanging the $m_{\tilde{\chi}}^{2}$ of the electroweakinos, but not the $m_{\tilde{\chi}}$.

[TABLE]

and the diagonal wave function renormalization constants101010In the Denner description $a=b=0$ and the $\Sigma^{L/R/Ls/Rs}_{ii}$ are replaced by $\widetilde{\text{Re}}\Sigma^{L/R/Ls/Rs}_{ii}$.

[TABLE]

where

[TABLE]

Our results coincide with those given in Ref. [19]. The above wave function renormalization constants have been chosen such that for all fermions the mixing terms are cancelled and the correct propagators are produced at the tree-level mass values. Note that in case we do not have enough parameters to renormalize all fermions on-shell, only some of these fermions satisfy OS conditions, i.e. their tree-level masses are equal to the physical masses. The remaining fermions have loop-corrected masses. This is the case for the electroweakino sector. Given the fact that we have already fixed the renormalization scheme for the Higgs sector, only the two gaugino masses $M_{1}$ and $M_{2}$ remain to be renormalized, while we have 7 masses (5 neutralino and 2 chargino masses), so that only two of them can be set OS. The remaining 5 particles receive loop-corrections to their masses. At loop level the mixing between fermions is in general not vanishing any more, and the residues of the propagators are not unity. These effects should be taken into account if the loop corrections to the masses are large. This is not the case for the renormalization schemes chosen here, therefore we neglect these effects.

In the chargino and neutralino sector hence only the two counterterms $\delta\lvert M_{1}\rvert$ and $\delta\lvert M_{2}\rvert$ remain to be determined. There are 21 different ways to choose two out of the 7 masses for the OS conditions. We will consider here two different schemes. In the first scheme (OS1), we require the masses of the wino-like chargino $\tilde{\chi}_{i}^{+}$ and the bino-like neutralino $\tilde{\chi}_{k}^{0}$ to be OS. The bino-like neutralino is sensitive to $M_{1}$ while the wino-like chargino is sensitive to $M_{2}$. Note that we do not choose the chargino and neutralino by referring to a fixed index order since they may not be sensitive to $M_{1}$ or $M_{2}$. This can then lead to numerical instability, as was found in the MSSM [23, 24, 25, 26] and in [7, 8] for the NMSSM. We denote the tree-level masses for the neutralinos (charginos) by a small letter $m_{\tilde{\chi}^{0(+)}_{i}}$, and the loop corrected masses by a capital letter $M_{\tilde{\chi}^{0(+)}_{i}}$. In the OS scheme the tree-level masses are equal to the loop-corrected ones. We define the counterterm mass matrices of the chargino and neutralino sector in the interaction basis, $\delta M_{N}$ and $\delta M_{C}$, through

[TABLE]

with the neutralino mass matrix given in Eq. (2.13) and the chargino mass matrix in Eq. (2.20). The counterterms for $M_{1}$ and $M_{2}$ are then given by

[TABLE]

In the second scheme (OS2), we use the masses of the bino-like neutralino, denoted by $\tilde{\chi}_{k}^{0}$, and of the wino-like neutralino, denoted by $\tilde{\chi}_{l}^{0}$, as inputs. The renormalization conditions for their OS masses are given by

[TABLE]

This results in the two solutions for the counterterms $\delta|M_{1}|$ and $\delta|M_{2}|$,

[TABLE]

with

[TABLE]

For the field renormalization constants of the charginos and neutralinos, we impose the OS conditions for the tree-level masses. Besides the two OS schemes, we will also adopt the $\overline{\text{DR}}$ renormalization scheme for $M_{1}$ and $M_{2}$, while for the field renormalization constants we use the OS conditions.

3.1.3 The Squark Sector

We consider here only the third-generation squarks. The results for the first- and second-generation squarks are obtained analogously. There are seven parameters to be renormalized in this sector,

[TABLE]

where $A_{t},A_{b}$ are complex and the mass terms are real. We denote their corresponding counterterms as

[TABLE]

and define the squark-mass counterterm matrices as

[TABLE]

with $M_{\tilde{q}}$ given in Eq. (2.27) for the stops and in Eq. (2.28) for the sbottoms. The renormalization of the remaining parameters appearing in the squark mass matrices has been specified in the renormalization of the Higgs sector, more specifically see Refs. [14, 15, 16, 17].

We have to renormalize the squark fields in order to make the squark self-energies finite. Here we use both the EMT and the Denner description. For the EMT description we have to introduce two separate WFR constants, one for the particle and one for the anti-particle. We introduce the squark WFR constants for the particles and anti-particles in the mass eigenstate basis as111111In the Denner description, we have $\delta\bar{Z}_{\tilde{q}}=\delta Z_{\tilde{q}}^{\dagger}$.

[TABLE]

The renormalized self-energies in the mass eigenstate basis are given by ($i,j=1,2$)

[TABLE]

where we denote by $\Sigma^{\tilde{q}}_{ij}$ the unrenormalized self-energies for the $\tilde{q}_{i}^{*}\to\tilde{q}_{j}^{*}$ transition.121212Note that in the real NMSSM the unrenormalized self-energies for the $\tilde{q}_{i}^{*}\to\tilde{q}_{j}^{*}$ transition and for the $\tilde{q}_{i}\to\tilde{q}_{j}$ transition are identical. They are different, however, in the complex case. In the following, we give the OS counterterms. The $\overline{\text{DR}}$ counterterms are then easily obtained by taking only the divergent parts of the corresponding OS counterterms.

Applying the decomposition of the fermionic self-energies as given in Eq. (3.42), the mass counterterm in the OS scheme for the top and bottom quark, respectively, is given by $(q=t,b)$131313In the Denner description, Re  is replaced by $\widetilde{\text{Re}}$.

[TABLE]

The OS conditions for the scalar renormalized self-energies are ($i,j=1,2$)141414In the Denner description, $\text{Re}\,\hat{\Sigma}^{\tilde{q}}_{ii}$ is replaced by $\widetilde{\text{Re}}\hat{\Sigma}^{\tilde{q}}_{ii}$ in Eqs. (3.82) and (3.84) while in Eq. (3.83) $\text{Re}\,\hat{\Sigma}^{\tilde{q}}_{ij}$ is replaced by $\widetilde{\text{Re}}\hat{\Sigma}^{\tilde{q}}_{ij}$.

[TABLE]

We apply the conditions in Eqs. (3.82) and (3.83) for the top squark to determine $\delta A_{t}$, $\delta m_{\tilde{Q}_{3}}^{2}$ and $\delta m_{\tilde{t}_{R}}^{2}$,151515In the Denner description, $\text{Re}\,\Sigma^{\tilde{t}}_{ii}$ in Eqs. (3.88) and (3.89) is replaced by $\widetilde{\text{Re}}\Sigma^{\tilde{t}}_{ii}$. Note, however, that Eq. (3.90) is the same in both the EMT and the Denner description. We use $\widetilde{\text{Re}}$ in the definition of $\delta Y_{\tilde{q}}$ ($q=t,b$) so that $\delta A_{q}=(\delta A_{q}^{*})^{*}$. The contribution from the imaginary part of the loop integrals is then moved into $\delta Z_{ij}^{\tilde{q}}$ and $\delta\bar{Z}_{ij}^{\tilde{q}}$.

[TABLE]

where

[TABLE]

There remain two parameters from the bottom squark sector to be determined, $A_{b}$, $m_{\tilde{b}_{R}}^{2}$. We choose the OS scheme where the bottom squark with the dominant contribution from the right-handed sbottom, which we denote by $\tilde{b}_{i_{R}}$, is OS and the mixing between the two bottom squark states vanishes. The three counterterms $\delta m_{\tilde{b}_{R}}^{2},\text{Re}\,\delta A_{b},\text{Im}\,\delta A_{b}$ are then obtained by solving three linear equations

[TABLE]

where $(x,y,z)=(\delta m_{\tilde{b}_{R}}^{2},\text{Re}\,\delta A_{b},\text{Im}\,\delta A_{b})$ and161616In the Denner description, $\text{Re}\,\Sigma^{\tilde{b}}_{i_{R}i_{R}}$ in Eq. (3.94) is replaced by $\widetilde{\text{Re}}\Sigma^{\tilde{b}}_{i_{R}i_{R}}$.

[TABLE]

where we have introduced the notation $\delta Y_{\tilde{b}}$ for later use. The other bottom squark mass gets loop corrections.171717Note that in principle all four masses of the stops and sbottoms can be renormalized on-shell simultaneously by adapting the input parameters appropriately. Its loop-corrected mass $M_{\tilde{b}_{j}}$ is obtained by solving iteratively the following equation181818In the Denner description $\text{Re}\,\hat{\Sigma}^{\tilde{b}}_{jj}$ is replaced by $\widetilde{\text{Re}}\hat{\Sigma}^{\tilde{b}}_{jj}$.

[TABLE]

We stop the iteration when the difference between two consecutive solutions is less than $10^{-5}$.

The OS wave function renormalization constants for the squarks are given by191919In the Denner description, $\text{Re}\,\Sigma^{\tilde{q}}_{ii}$ in Eq. (3.98) is replaced by $\widetilde{\text{Re}}\Sigma^{\tilde{q}}_{ii}$ and $\text{Re}\,\Sigma^{\tilde{q}}_{12},\text{Re}\,\Sigma^{\tilde{q}}_{21}$ in Eqs. (3.99),…,(3.102) are replaced by $\widetilde{\text{Re}}\Sigma^{\tilde{q}}_{12},\widetilde{\text{Re}}\Sigma^{\tilde{q}}_{21}$, respectively.

[TABLE]

where the $\delta Y_{\tilde{q}}$ are given in Eq. (3.90) and Eq. (3.96) for $q=t$ and $q=b$, respectively. Besides the OS scheme we also provide the option to use the $\overline{\text{DR}}$ scheme for all parameters and the wave function renormalization constants. In the $\overline{\text{DR}}$ scheme, all squarks receive loop-corrections to their masses. We will discuss this issue in Subsection 3.4.

3.2 Loop-corrected Higgs Boson Masses and Mixings

Since we use the mixed OS-$\overline{\text{DR}}$ renormalization scheme for the Higgs sector parameters together with the $\overline{\text{DR}}$ scheme for the Higgs fields, all Higgs bosons are mixed and receive loop corrections to their masses. For the evaluation of the loop-corrected Higgs boson masses and the Higgs mixing matrix, we use the numerical results obtained from NMSSMCALC [13]. In this code the two-loop corrected Higgs boson masses are obtained by determining the zeros of the determinant of the two-point function $\hat{\Gamma}^{h}(p^{2})$ with

[TABLE]

where $m_{h_{i}}$ are the tree-level masses and $\hat{\Sigma}_{ij}^{h}(p^{2})$ is the renormalized self-energy of the $h_{i}\to h_{j}$ transition at $p^{2}$. In NMSSMCALC, we have included in the renormalized Higgs self-energies the complete one-loop contributions with full momentum dependence [14, 15] and the two-loop contributions of ${\cal O}(\alpha_{t}\alpha_{s})$ [12] and of ${\cal O}(\alpha_{t}^{2})$ [17] in the gaugeless limit at zero momentum. The loop-corrected masses of the Higgs bosons are then sorted by ascending masses202020We denote loop-corrected Higgs mass eigenstates by capital letters $H_{i}$ ($i=1,...,5$).

[TABLE]

We have improved the stability of the determination of the Higgs boson masses in NMSSMCALC by implementing two-point loop integrals with complex momentum. In the old version of NMSSMCALC, in order to take into account the contribution of the imaginary part of the complex momentum we expanded the renormalized Higgs self-energies around the real part of the complex momentum as

[TABLE]

Note that this was done only for the one-loop correction with full momentum dependence. This expansion is not good when $\text{Re}\,p^{2}$ is close to threshold regions in which $\partial\hat{\Sigma}_{ij}^{h}(\text{Re}\,p^{2})/\partial\text{Re}\,p^{2}$ contains threshold singularities. To overcome this problem one can use complex masses for the loop particles or complex momenta. Using complex masses requires the decay widths of the particles. These have to be obtained in an iterative procedure which is very time consuming. We have decided to use the complex momenta and to keep the masses real. We have implemented the two-point loop integrals with complex momenta and therefore do not use any more the mentioned approximation. We have confirmed that the evaluation of the Higgs masses is stable in the singularities region and the differences between Higgs masses using the complex momentum and the expansion in Eq. (3.105), defined as $(M_{H_{i}}^{\text{expansion}}-M_{H_{i}}^{\text{complex~{}}p^{2}})/M_{H_{i}}^{\text{expansion}}$ ($i=1,\ldots,5$) are of per mille level.

In processes with external Higgs bosons finite wave-function renormalization factors ${\bf Z}^{H}$ have to be taken into account in order to ensure the on-shell properties of these Higgs bosons. The wave-function renormalization factor matrix performing the rotation to the OS states is given by [27]

[TABLE]

where

[TABLE]

with

[TABLE]

Here prime denotes the derivative with respect to $p^{2}$.

3.3 Loop-corrected Neutralino and Chargino Masses

Within the OS and $\overline{\text{DR}}$ schemes defined for the neutralino and chargino sector, the electroweakinos cannot all be renormalized OS, and there remain neutralinos and charginos that receive loop corrections to their masses. In the following we define our procedure to determine the loop-corrected masses for fermions in the general case where mixing contributions are present212121In the literature there exist many papers that deal with loop-corrected masses for neutralinos and charginos in the MSSM such as Refs. [28, 29, 23] to name a few of them. For the definition of the loop-corrected masses applied there, OS conditions for the field renormalization constants were applied to eliminate mixing effects between fermions at tree-level mass values.. This procedure will be used for both the OS and the $\overline{\text{DR}}$ scheme. To give an intuitive definition of the loop-corrected fermion masses, we express the propagator matrix in terms of left- and right-handed Weyl spinors, $\psi_{D}=(\psi_{L},\psi_{R})^{T}$. In this basis, the tree-level propagator matrix is given by

[TABLE]

with

[TABLE]

where $\vec{\sigma}=(\sigma^{1},\sigma^{2},\sigma^{3})$ denotes the three Pauli matrices. The mass matrix is given by

[TABLE]

The loop-corrected propagator matrix in the Weyl basis $(\psi_{L},\psi_{R})^{T}$ is defined as

[TABLE]

The loop-corrected mass $M$ is determined from the real pole $p^{2}=M^{2}$ of the propagator matrix satisfying the equation

[TABLE]

The solution of Eq. (3.113) is given by

[TABLE]

which can be solved iteratively. When the fermion is OS, $p^{2}=m^{2}$, the above relation is nothing else but the OS condition obtained from Eq. (3.51). For the case of $n$ Dirac spinors, the 1PI two-point function in the basis $(\psi_{L}^{1},\psi_{L}^{2},\cdots,\psi_{L}^{n},\psi_{R}^{1},\psi_{R}^{2},\cdots,\psi_{R}^{n})$ is a $2n\times 2n$ matrix

[TABLE]

The matrices $a,b,c,d$ are $2\times 2$ matrices in case of charginos and $5\times 5$ matrices for neutralinos. With $m_{\tilde{\chi}_{i}}$ generically denoting the mass of an electroweakino with index $i$, the matrices are given by

[TABLE]

with $i,j=1,2$ for charginos and $i,j=1,...,5$ for neutralinos.

The poles of the propagator matrix are the solutions of the equation

[TABLE]

This is equivalent to [30],

[TABLE]

In practice, we solve the equation numerically through iteration together with the diagonalization of the mass matrix $M_{\tilde{\chi}}=c(p^{2})b^{-1}(p^{2})d(p^{2})a^{-1}(p^{2})$ to obtain the complex poles. The loop-corrected masses are then obtained from the real parts of these complex poles. This procedure is applied for the calculation of the loop-corrected masses for neutralinos using the OS definition of the neutralino WFR constants. However, for the chargino sector the mass matrix $M_{\tilde{\chi}^{\pm}}$ contains infrared (IR) divergences at arbitrary momentum. We have implemented the approximation used in Ref. [23] and calculate the loop-corrected chargino masses by using the formula

[TABLE]

3.4 Loop-corrected Squark Masses and Mixings

In our $\overline{\text{DR}}$ scheme, the counterterms

[TABLE]

contain only the UV divergent parts. For the renormalization of the squark fields we use a modified OS scheme. In the following, we will describe the details of this scheme. We first compute the $\overline{\text{DR}}$ squark WFR constants by taking the UV-divergent parts of the OS counterterms, defined in Eqs. (3.98) to (3.102), with tree-level mass values. Using these $\overline{\text{DR}}$ squark WFR constants, we then compute the loop-corrected squark masses $M_{\tilde{q}_{i}}$ that are obtained by solving iteratively the equations

[TABLE]

We have assumed here that the off-diagonal renormalized self-energies vanish and that the residues of the propagators are unity at the loop-corrected masses. This is equivalent to redefining the diagonal squark WFR constants at loop-corrected masses as

[TABLE]

and the off-diagonal WFR constants as

[TABLE]

In the above equations the renormalized self-energies are computed with the $\overline{\text{DR}}$ squark WFR constants. We keep also the imaginary part of the two-point loop integrals in the renormalized self-energies. Note that the WFR constants $\delta Z^{\tilde{q}}_{ij}(M_{\tilde{q}_{k}}^{2})$ will enter the evaluation of the decay width. The diagonal WFR constants $\delta Z^{\tilde{q}}_{ii}(M_{\tilde{q}_{k}}^{2})$ contain IR divergences evaluated at the loop-corrected masses. These IR divergences will cancel exactly with those arising from the virtual part and the real radiation contributions which are also evaluated at the loop-corrected masses. We have verified that this statement is true for both the EW and the QCD corrections.

4 Higher-Order Corrections to the Two-Body Decays of the

Neutral Higgs Bosons

In this section, we present those two-body decay channels that we have improved by including the missing NLO EW and QCD corrections. These channels are the decays into OS SM fermion pairs, OS massive gauge bosons, into a pair of Higgs and gauge bosons, into chargino or neutralino pairs and into top or bottom squark pairs. We will not discuss decays into gluon pairs, photon pairs or $Z\gamma$ which can be found in Ref. [13]. For these decays, NLO EW corrections are of two-loop order as the leading order (LO) decay widths are already loop-induced. The inclusion of the NLO EW corrections to Higgs-to-Higgs decays on the other hand have been presented in Ref. [11] and the dominant two-loop corrections of the ${\cal O}(\alpha_{t}\alpha_{s})$ have been provided in Ref. [12].

For our computation we have used several programs. The generation of the amplitudes was done by FeynArts [31, 32] using a model file created by SARAH [33, 34, 35, 36]. The output amplitudes were further processed using FeynCalc [37, 38] for the simplification of the Dirac matrices and for the tensor reduction. The one-loop integrals were evaluated with the help of LoopTools [39].

4.1 Higgs Boson Decays into Fermion Pairs

In order to make use of the published results of higher-order corrections in the literature for CP-even and CP-odd Higgs bosons, it is convenient to write the interaction vertex of the complex NMSSM Higgs boson $h_{i}$ $(i=1,\ldots,5)$ and quarks as

[TABLE]

where the scalar and pseudoscalar coupling coefficients for the up- and down-type quarks at tree-level are given by

[TABLE]

where ${\cal R}_{ij}$ ($i,j=1,5$) denotes the matrix elements of the mixing matrix rotating the tree-level Higgs gauge eigenstates to the mass eigenstates, see Eq. (2.8).

Following the prescription outlined in our publication [13], we improve the widths of the Higgs boson decays into quark pairs by including the missing SUSY–QCD and SUSY–EW corrections. We decompose the EW corrections into the known QED corrections arising from a virtual photon exchange and a real photon emission and the remaining unknown EW corrections from the genuine EW one-loop diagrams.

The one-loop SUSY–QCD corrections originate from loop diagrams with the exchange of a gluino, while the SUSY–EW corrections stem from loop diagrams with weak gauge bosons $W,Z$, fermions, Higgs bosons and their superpartners in the internal lines. They are both IR finite quantities. The computation of the Higgs boson decays into a bottom quark pair shows that the bottom quark mass counterterm contains terms proportional to ${t_{\beta}}$. This contribution is large in the large-${t_{\beta}}$ regime and universal. In many cases, this contribution is the leading part of the SUSY–QCD and SUSY–EW corrections and can be absorbed into an effective bottom quark Yukawa coupling. This can be done by using an effective Lagrangian formalism [40, 41, 42]. In Ref. [13], we have presented the effective bottom Yukawa couplings in the real and complex NMSSM. We do not repeat every detail here but only quote the relevant formulae. In Eq. (4.128) we have given the tree-level scalar and pseudoscalar coupling coefficients appearing in the Feynman rule for the CP-violating Higgs bosons $h_{i}$ to a bottom-quark pair. The Feynman rule for the effective coupling including the leading SUSY–QCD and SUSY–EW corrections [40, 41, 42, 43, 44, 45, 46, 47, 48, 49] (denoted by a tilde to mark the inclusion of the corrections) is also decomposed into a scalar and a pseudoscalar part and reads [13]

[TABLE]

with

[TABLE]

where

[TABLE]

The correction $\Delta_{b}$ including the leading SUSY–QCD and SUSY–EW corrections can be cast into the form

[TABLE]

with the one-loop corrections given by

[TABLE]

where $\alpha_{t}=y_{t}^{2}/(4\pi)$ with $y_{t}=\sqrt{2}m_{t}/(v\sin\beta)$ is the top-Yukawa coupling and $C_{F}=4/3$. The generic function $I$ is defined as

[TABLE]

Note that the scale of $\alpha_{s}$ in the SUSY–QCD corrections has been set equal to $\mu_{R}=(m_{\tilde{b}_{1}}+m_{\tilde{b}_{2}}+|M_{\tilde{g}}|)/3$, while in the SUSY–EW corrections it is $\mu_{R}=(m_{\tilde{t}_{1}}+m_{\tilde{t}_{2}}+|\mu|)/3$. The strong coupling constant $\alpha_{s}$ is evaluated with five active flavors.

The decay width of the CP-violating NMSSM Higgs bosons $H_{i}$ into $q\bar{q}$, including the QCD, SUSY–QCD, QED, EW and SUSY–EW corrections, can then be cast into the form

[TABLE]

where $x_{q}=m_{q}^{2}/M_{H_{i}}^{2}$, and

[TABLE]

and

[TABLE]

In the numerical analysis presented in Sec. 5 we will use the quantity $\Gamma^{\text{SEW(+QCD)}}$ for the decays into fermion pairs to denote the partial decay widths including the SUSY-EW and (for the quarks) SUSY-QCD corrections, i.e. exactly the partial decay width as defined in Eq. (4.138) with the loop-corrected $\Gamma^{S}$ and $\Gamma^{P}$ given in Eqs. (4.139) and (4.140), respectively. In contrast, we will denote by $\Gamma^{\text{tree}}$ the partial decay widths that only include the $\Delta_{b}$ corrections, i.e. Eq. (4.138) but with $\Gamma^{S}$ and $\Gamma^{P}$ given by the first lines in Eqs. (4.139) and (4.140), respectively.

Note that in $\delta{\cal M}^{\text{rem},S/P}$, $\delta^{S/P}_{\text{sub}}$ and $\delta{\cal M}_{GZ,\text{mix}}$ (which will be explained below) we use the tree-level Higgs couplings $g_{h_{k}q\bar{q}}$ to the quarks. We take the occasion to remind the reader that tree-level mass eigenstates are always denoted by $h_{i}$ and loop-corrected ones by $H_{i}$. Unless stated otherwise, this means that we use tree-level couplings for external $h_{i}$ but with loop-corrected masses, and for particles inside loop diagrams we always use tree-level masses and tree-level couplings. We comment on the various terms appearing in Eqs. (4.138), (4.139) and (4.140) one by one. The one-loop QED corrections, denoted by $\Delta_{\text{QED}}$ have been known in the SM for a long time, cf. Refs. [50, 51, 52, 53, 54, 55]. In the limit $m_{q}\ll M_{H_{i}}$ they are given by Ref. [56]222222Note that actually in the code we have programmed the QCD corrections for the completely massive case at next-to-leading order, translated to the $\overline{\mbox{MS}}$ scheme, and interpolated with the massless expression for large Higgs masses, according to the implementation in HDECAY [57, 58]. ,

[TABLE]

where $Q_{q}$ denotes the electric charge of the quark $q$. The one-loop SM QCD corrections are similar to the QED corrections, with the replacement of $Q_{q}^{2}\,\alpha$ by $(4/3)\alpha_{s}(M_{H_{i}}^{2})$. It is well-known that the SM QCD corrections are rather large. The large logarithmically enhanced part has been absorbed in Eq. (4.138) into the running $\overline{\text{MS}}$ quark mass $\overline{m}_{q}(M_{H_{i}}^{2})$ at the corresponding Higgs mass scale $M_{H_{i}}$ to improve the convergence of the perturbative expansion. The QCD corrections can be taken over from the MSSM case by adapting the Higgs couplings [50, 51, 52, 59, 53, 60, 61, 62, 63, 64, 65, 66, 67, 68]. After subtracting the enhanced part, the remaining QCD correction $\Delta_{\text{QCD}}$ reads

[TABLE]

where $N_{F}=5$ active flavors are taken into account. In the CP-conserving case we also include top quark induced corrections $\Delta_{t}^{S/P}$ by adding them to $\Delta_{\text{QCD}}$. They can be taken over from the MSSM case and read [50, 51, 52, 59, 53, 60, 61, 62, 63, 64, 65, 66, 67],

[TABLE]

In the decay into a bottom quark pair, the large universal corrections proportional to ${\cal O}(\alpha_{s}{t_{\beta}},\alpha_{b}{t_{\beta}})$ are resummed into the effective bottom Yukawa couplings $\tilde{g}_{h_{i}b\bar{b}}^{S,P}$, given in Eqs. (4.131) and (4.132), while in the decay into top quarks we use the tree-level values of the effective couplings, i.e.

[TABLE]

with $g_{h_{i}t\bar{t}}$ given in Eq. (4.129). The remaining SUSY–QCD and SUSY–EW corrections are collected in $\delta{\cal M}^{\text{rem},S}$, $\delta{\cal M}^{\text{rem},P}$ and $\delta{\cal M}_{GZ,\text{mix}}$, where

[TABLE]

with $\delta{\cal M}^{\text{counter},S/P}$ denoting the counterterm contributions. Since in the decay into bottom quarks we have resummed the dominant corrections into the effective couplings, in the remaining SUSY–QCD and SUSY–EW correction we have to subtract these corrections to avoid double counting by adding appropriate counterterms. This is taken care of by the last terms in Eqs. (4.139) and (4.140), respectively, to which we will come back below.

The term $\delta{\cal M}_{GZ,\text{mix}}$ is the sum of the contributions from the mixing of the CP-odd component of the Higgs bosons with the neutral Goldstone boson $G$ and with the $Z$ boson, respectively. We use the tree-level masses for the Higgs bosons in the loops in order to ensure the proper cancellation of the UV-divergent pieces, but we use the loop-corrected Higgs boson masses for the external particles in the evaluation of the wave-function renormalization factors, of the amplitudes, and of the decay widths. It is well-known that the use of loop-corrected Higgs boson masses for the external particles violates the Slavnov-Taylor identity of the amplitude $\delta{\cal M}_{GZ,\text{mix}}$ [11, 69, 27]. To restore the gauge symmetry one should use the tree-level masses for both the external and the internal Higgs bosons. This causes a mismatch with the phase-space factor (where we use the loop-corrected Higgs boson masses) and with the evaluation of the other amplitudes. We therefore use the loop-corrected masses for the external particles also in these contributions, which are then computed in the unitary gauge. Note that we apply the same method also for the other decays that contain the contribution $\delta{\cal M}_{GZ,\text{mix}}$. In Ref. [11], $\delta{\cal M}_{GZ,\text{mix}}$ was computed only for the $i$th tree-level mass eigenstate of the decay $H_{i}\to q\bar{q}$. This may cause instabilities in the case of large mixing between Higgs boson mass eigenstates at loop-level. We avoid this by multiplying it also with the WFR factor ${\bf Z}^{H}$.

The remaining SUSY–QCD and SUSY–EW corrections are computed in the Feynman diagrammatic approach. The corrections consist of the contributions from genuine one-loop diagrams and the counterterms. The counterterms are given by

[TABLE]

with the expressions for the scalar and pseudoscalar parts, $\delta\lambda_{h_{i}q\bar{q}}^{S/P}$, reading232323Note, that the angle $\beta$ in the sense of a mixing angle does not get renormalized.

[TABLE]

with

[TABLE]

The counterterm $\delta v$ for the VEV is given by

[TABLE]

The electric coupling $e$, the $W$ and $Z$ boson masses and $\tan\beta$ are renormalized according to the Higgs sector. The top and bottom quarks are renormalized OS using both the EMT and the Denner descriptions. The terms being related to left-handed and right-handed OS wave function renormalization constants for the quarks ($\delta Z_{qL/R}$) and anti-quarks ($\delta\bar{Z}_{qL/R}$) ($q=t,b$) are242424In the Denner description, $\Sigma_{q}$ is replaced by $\widetilde{\text{Re}}\Sigma_{q}$.

[TABLE]

The $\overline{\text{DR}}$ wave function renormalization constants for the Higgs bosons are denoted by $\delta Z_{h_{i}h_{k}}$. We have checked the UV-finiteness of the SUSY–QCD and SUSY–EW corrections to the decay amplitude.

As mentioned above, in the decay into bottom quarks we have to take care to avoid double counting after resumming the dominant part of the SUSY–QCD and SUSY–EW corrections into the effective bottom coupling. To subtract these contributions we add in the decays into a $b$-quark pair the following counterterms to Eq. (4.139) and Eq. (4.140),

[TABLE]

where $\Delta_{b}$ equals $\Delta_{b}^{\text{QCD}}-\Delta_{1}$ and $\Delta^{\text{elw}}_{b}$ for the SUSY–QCD and SUSY–EW corrections, respectively. For the decays into a top-quark pair, these contributions are

[TABLE]

In the decays into strange quarks we also include the one-loop SUSY–QCD corrections. They are obtained after substituting $\Delta_{b}$ as given in Eq. (4.133) with

[TABLE]

The decays into charm quarks are treated as the decays into top quarks, with the appropriate replacements.

The decays into lepton final states $l=e,\mu,\tau$ do not receive QCD corrections. Their SUSY–EW corrected decay width is given by

[TABLE]

where $x_{l}=m_{l}^{2}/M_{H_{i}}^{2}$. The $\Delta_{\text{QED}}^{S,P}$ are given by Eq. (4.141) after replacing $(Q_{q},m_{q})$ with $(Q_{l},m_{l})$. Furthermore, we resum the dominant SUSY–EW corrections into the effective couplings $\tilde{g}_{h_{i}ll}^{S,P}$. They are obtained from the effective couplings $\tilde{g}_{h_{i}b\bar{b}}$ in Eqs. (4.131) and (4.132) after replacing $\Delta_{b}$ with $\Delta_{l}$, where $\Delta_{l}$ in the complex case is given by

[TABLE]

The contributions $\Gamma^{S,P}_{H_{i}\to l\bar{l}}$ are obtained from the ones given in Eqs. (4.139) and (4.140) after replacing $q\to l$.

4.2 Higgs Boson Decays into $W^{+}W^{-}$ and $ZZ$

We now address the higher-order corrections to the Higgs boson decays into gauge boson pairs. We consider corrections only for on-shell decays. Off-shell decays are still treated at tree-level as done in NMSSMCALC [13]. The one-loop corrected decay amplitude for the decay of a CP-violating NMSSM Higgs boson $H_{i}$ ($i=1,...,5$) into a pair of massive gauge bosons $V=Z,W^{\pm}$, $H_{i}(p)\to V(k_{1})V(k_{2})$, is given by

[TABLE]

where $\epsilon_{\mu}(k_{1})$ and $\epsilon_{\nu}(k_{2})$ are the polarization vectors of the two gauge bosons with four-momenta $k_{1}$ and $k_{2}$, respectively. Note that the $GZ,\text{mix}$ contribution vanishes. The tree-level amplitudes for the two final-state pairs are

[TABLE]

with

[TABLE]

And the tensor structure of the NLO corrections is given by

[TABLE]

The genuine one-loop triangle diagram contributions are denoted by ${\cal M}_{\text{1L}}^{(i)}$ ($i=1,2,3$) and the counterterm contribution by ${\cal M}_{\text{1L}}^{\text{CT}}$. The term ${\cal M}_{\text{1L}}^{(3)}$ vanishes in the CP-conserving case.

The decay width for the decay $H_{i}\to VV$, including the NLO corrections, is given by

[TABLE]

where $R=1/2$ for $V=Z$ and $R=1$ for $V=W$. The improved tree-level decay width reads

[TABLE]

with

[TABLE]

and

[TABLE]

The NLO partial width for the $Z$-boson pair final state contains only virtual contributions,

[TABLE]

For the $W$-boson pair final state it contains both virtual and real radiation contributions,

[TABLE]

The virtual part is given by

[TABLE]

with

[TABLE]

The formulae for ${\cal M}_{\text{1L}}^{(1)}$ and ${\cal M}_{\text{1L}}^{(2)}$ are quite lengthy and we do not display them explicitly here. The counterterm contributions for $V=W$ and $Z$, respectively, read

[TABLE]

For the decay $H_{i}\to W^{+}W^{-}$ we have to include to the contribution from the radiation of a real photon in order to get an infrared-finite result. This contribution is given by

[TABLE]

where

[TABLE]

The formula is in agreement with the result given in Ref. [70]. Here, we have neglected the arguments of the Bremsstrahlung integrals $I_{i_{1},\ldots,i_{2}}^{j_{1},\ldots,j_{2}}(M_{H_{i}},M_{W},M_{W})$ for the sake of readability. In terms of the photon momentum $q$, the $W^{+}$ momentum $k_{1}$ and the $W^{-}$ momentum $k_{2}$, these integrals are defined as

[TABLE]

where $k_{0}$ is the four-momentum of $H_{i}$ and $i_{l},j_{k}=0,1,2$. The plus signs correspond to $k_{1}$, $k_{2}$, the minus sign to $k_{0}$. Their analytic expressions are given in Ref. [22].

We have checked the UV finiteness of the NLO decay widths of both the $H_{i}\to ZZ$ and the $H_{i}\to W^{+}W^{-}$ decays. The IR divergence in the decay $H_{i}\to W^{+}W^{-}$, however, is more demanding. At strict one-loop level, *i.e. *one must use the tree-level Higgs boson mass for the external Higgs boson and the unity WFR factor matrix, the IR finiteness if fulfilled. However, the use of the loop-corrected Higgs masses and the WFR factors ${\bf Z}^{H}_{ij}$ breaks the IR finiteness, because different orders of perturbation theory are mixed in this case. At tree level there exists a relation between the coupling of the neutral Higgs boson with the charged Goldstone bosons and the coupling of the neutral Higgs boson with the $W$ bosons. Defining the Lagrangian for the interaction between Higgs and Goldstone bosons by

[TABLE]

it is given by

[TABLE]

with the tree-level Higgs boson mass $m_{h_{j}}$. In order to obtain an IR-finite result while using the loop-corrected mass $M_{H_{i}}$, we chose to modify the coupling $g_{h_{j}G^{+}G^{-}}$ as

[TABLE]

where the tree-level mass $m_{h_{j}}^{2}$ has been replaced by the loop-corrected mass $M_{H_{i}}^{2}$ of the decaying Higgs boson $H_{i}$. We verified that the modification of this coupling ensures IR finiteness while not affecting UV finiteness. The same method has been used in Ref. [71]. While taking the loop-corrected mass for the external decaying Higgs boson ensures compatibility with the observation of a 125 GeV SM-like Higgs boson, this approach breaks gauge invariance, however. For more details on this issue, we refer to an upcoming publication [72].

For the non-SM-like neutral Higgs bosons, the tree-level coupling $g_{h_{i}VV}$ is in general suppressed, in particular in case the Higgs boson with mass 125.09 GeV behaves very SM-like. In this case, the one-loop corrected decay width $\Gamma^{\text{1L}}(H_{i}\to VV)$ can be even larger than the tree-level improved one $\Gamma^{\text{tree}}(H_{i}\to VV)$. This becomes a problem when the one-loop correction is negative, as then the one-loop corrected partial decay width becomes negative. In this case, we have to include the one-loop squared term, which is formally of higher order. For the decay $H_{i}\to ZZ$, we will include the one-loop squared contribution.252525Note, however, that the thus obtained result has to be taken with caution. The complete two-loop calculation contributes further terms that might cause the complete two-loop result to differ considerably from the result obtained in the here applied approximation. Moreover, the inclusion of (part of) the two-loop corrections explicitly includes a dependence on the renormalization scheme chosen at one-loop order that would need to be cancelled by transforming the input parameters appropriately so as not to become inconsistent. We still use this approach in order to obtain physical, i.e. positive, partial decay widths and hence physical branching ratios. Since the partial decay width is suppressed here anyway, the effect of the difference between the approximation and the full two-loop result on the branching ratio is expected to be subleading. Still, the code NMSSMCALC will print out a warning to make the user aware of this issue. In particular, the decay width now is given by

[TABLE]

where

[TABLE]

where ${\cal M}_{\text{1L}}^{(1)+\text{CT}}$ is the sum of ${\cal M}_{\text{1L}}^{(1)}+{\cal M}_{\text{1L}}^{\text{CT}}$. For the decay $H_{i}\to W^{+}W^{-}$, the form factor ${\cal M}_{\text{1L}}^{(1)+\text{CT}}$ contains IR divergences so that we cannot treat it in the same way as in the decay $H_{i}\to ZZ$. Note, however, that the 1L decay width $\Gamma^{\text{1L}}(H_{i}\to WW)$ can always be divided into three parts that are separately UV and IR finite: the (s)fermion contribution arising from loops containing SM model fermions and their superpartners, the chargino/neutralino contribution from loops with internal charginos and neutralinos and the gauge/Higgs contribution from loops with gauge and Higgs particles. In many cases the dominant contribution is the (s)fermion part. That was also observed in the MSSM case [73, 71]. We therefore add to the one-loop corrected decay with $H_{i}\to W^{+}W^{-}$ the one-loop squared contribution from the (s)fermion part only. This part is IR finite as it solely involves fermions and sfermions but no photons. Both for the decays into $ZZ$ and into $WW$ we include the respective one-loop squared terms in case the one-loop contribution is larger than 80% of the tree-level decay width.

4.3 Higgs Boson Decays into a $Z$ Boson and a Higgs Boson

The one-loop corrected amplitude for the decay of a heavy Higgs boson $H_{i}$ with four-momentum $p$ into a light Higgs boson $H_{j}$ and a $Z$ boson, with four-momenta $k_{1}$ and $k_{2}$, respectively, $H_{i}(p)\to H_{j}(k_{1})Z(k_{2})$, can be written as

[TABLE]

where

[TABLE]

The tree-level expression ${\cal M}^{(0)}_{h_{i}h_{j}Z}$ reads

[TABLE]

and the one-loop term ${\cal M}^{(1)}_{h_{i}h_{j}Z}$ consists of the genuine one-loop diagram contribution and the counterterm part given by

[TABLE]

Also here, the contribution from the one-loop diagrams with the transition $h_{i}\to Z(G)$, ${\cal M}^{GZ,\text{mix}}_{h_{i}h_{j}Z}$, is calculated in the unitary gauge. The improved tree-level decay width is given by

[TABLE]

and the NLO decay width by

[TABLE]

with the 2-particle phase-space factor

[TABLE]

Since the formulae for the 1-loop amplitudes are quite lengthy we do not display them explicitly here. Note that, as in the decay into massive gauge bosons, in Eq. (4.191) we also included, keeping in mind the caveat mentioned there, one-loop contributions squared as the one-loop corrections can be large and negative.

4.4 Higgs Boson Decays into Neutralinos and Charginos

The couplings of a neutral Higgs boson $h_{i}$ with the electroweakinos can be defined as

[TABLE]

where $\tilde{\chi}$ stands generically for the neutralinos and charginos and $g_{h_{i}\tilde{\chi}_{k}\tilde{\chi}_{j}}^{R}=\left(g_{h_{i}\tilde{\chi}_{k}\tilde{\chi}_{j}}^{L}\right)^{*}$. At tree level, the left- and right-handed coefficients for the Higgs-chargino couplings are given in terms of [13]

[TABLE]

where $i=1,\ldots,5$, $j,k=1,2$, and for the Higgs-neutralino couplings we have

[TABLE]

with $l,m=1,\ldots,5$. The decay width for the decay of a Higgs boson $H_{i}$ into a neutralino pair or a chargino pair including higher order corrections is given by

[TABLE]

where $R=1/2$ for identical final states and $R=1$ otherwise. The improved tree-level decay width reads

[TABLE]

where the 2-body phase space factor is

[TABLE]

in terms of the improved tree-level amplitude

[TABLE]

The one-loop decay width for the decay into a neutralino pair is given by

[TABLE]

and for the decay into a chargino pair it is

[TABLE]

where the virtual contribution can be cast into the form

[TABLE]

with the left- and right-handed one-loop amplitudes containing genuine triangle, counterterm and ’$GZ,\text{mix}$’ contributions,

[TABLE]

We do not display explicitly here the lengthy expressions for the triangle and ’$GZ,\text{mix}$’ contributions. The explicit expressions of the counterterm amplitudes for the decays into neutralinos are

[TABLE]

and for the decays into charginos they read

[TABLE]

The right-handed counterterm amplitudes are equal to the complex conjugate of the corresponding left-handed parts after interchanging the indices of the charginos and neutralinos in the final state. The real photon contribution for the decays into a chargino pair is expressed in terms of the Bremsstrahlung integrals as

[TABLE]

where the arguments of the Bremsstrahlung integrals $I_{i_{1},\ldots,i_{n}}^{j_{1},\ldots,j_{m}}(M_{H_{i}},M_{\tilde{\chi}^{+}_{j}},M_{\tilde{\chi}^{-}_{k}})$ have been neglected. Note that we use the loop-corrected masses for the external Higgs boson and the external charginos and neutralinos in the tree-level, virtual and real contributions. However, for particles inside loops we use the tree-level masses and tree-level couplings. This does not affect the UV-finiteness but can break the IR-finiteness in the decay into a pair of charginos. We overcome this problem by replacing the tree-level mass of the chargino in the loop diagrams with a photon by the corresponding loop-corrected chargino mass. Our treatment is different from Ref. [9] where the authors define an IR divergent counterterm to cancel the mismatch between the real and virtual contributions. Note finally, that in case the NLO decay width into neutralino final states becomes negative, the improved tree-level decay width is calculated instead in NMSSMCALCEW, including the ${\bf Z}^{H}$ factor.

4.5 Higgs Boson Decays into Squark Pairs

The NLO corrections to the decay of a neutral Higgs boson into a squark-antisquark pair consist of the QCD and EW corrections. In the CP-conserving NMSSM, the NLO corrections to the decay of a CP-odd Higgs boson into a stop pair have been calculated and discussed in Ref. [6]. We extend this computation to the CP-violating case and include also the decay into a sbottom-antisbottom pair in this paper. The NLO QCD corrections are positive and large. They can be larger than 100% as observed in Ref. [6] while the EW correction are negative and can be of up to $-40\%$. In our calculation, we have implemented both the OS and the $\overline{\text{DR}}$ scheme. We have three options here. First, the seven parameters are renormalized in the OS scheme. Second, the parameters of the stop sector, $m_{t},m_{\tilde{Q}_{3}},m_{\tilde{t}_{R}},A_{t}$, are renormalized in the OS scheme while the remaining parameters, $m_{b},m_{\tilde{b}_{R}},A_{b}$, are renormalized in the $\overline{\mbox{DR}}$ scheme. Third, all parameters are renormalized in the $\overline{\mbox{DR}}$ scheme. The loop-corrected decay width is decomposed into the improved tree-level, one-loop QCD and one-loop EW decay widths,

[TABLE]

Denoting the color factor by $N_{F}$, with $N_{F}=3$, the improved tree-level decay width is given by

[TABLE]

in terms of the improved tree-level amplitude

[TABLE]

The tree-level Higgs-squark-squark couplings are given by [13]

[TABLE]

with

[TABLE]

The one-loop QCD and EW contributions to the decay width are given by the sum of the virtual and real contributions, respectively,

[TABLE]

For the virtual QCD contribution we have

[TABLE]

with the 2-body phase space factor $R_{2}$ defined in Eq. (4.198). The expression for the virtual EW contribution is different from the QCD one due to an extra contribution containing the transition $h_{i}\to G,Z$. Explicitly, we have

[TABLE]

The explicit expressions for the counterterm contributions are quite lengthy and given in Appendix A. We do not display, however, the more cumbersome amplitudes of the virtual QCD and EW contributions, ${\cal M}^{\Delta,\text{QCD}}_{h_{i^{\prime}}\tilde{q}_{j}\tilde{q}_{k}^{*}}$ and ${\cal M}^{\Delta,\text{EW}}_{h_{i^{\prime}}\tilde{q}_{j}\tilde{q}_{k}^{*}}$, respectively.

The real photon radiation contribution in the EW corrections is expressed in terms of the Bremsstrahlung integrals as

[TABLE]

As usual, we have neglected the arguments of the Bremsstrahlung integrals $I_{l}(M_{H_{i}}^{2},M_{\tilde{q}_{j}}^{2},M_{\tilde{q}_{k}}^{2})$ and $I_{lm}(M_{H_{i}}^{2},M_{\tilde{q}_{j}}^{2},M_{\tilde{q}_{k}}^{2})$ ($l,m=1,2$ and $j,k=1,2$). The real gluon radiation contribution in the QCD corrections can be obtained from the EW real photon radiation contribution by replacing $Q_{q}^{2}\alpha$ with $C_{F}\alpha_{s}^{2}$, where $C_{F}=4/3$ for $SU(3)_{C}$. We have checked the UV and IR finiteness of the EW and QCD corrections. We have compared numerically with the NLO EW and QCD corrections in the OS scheme for the decay $A_{2}\to\tilde{t}_{1}\tilde{t}_{2}$ [6] using their description in the real NMSSM and found full agreement.

5 Numerical Results

To illustrate the importance of the higher-order corrections to the decays of the light and heavy neutral Higgs bosons and to test the stability of the NLO results in various regions of the parameter space we have performed a scan in the NMSSM parameter space. The parameter points are checked against compatibility with the experimental constraints from the Higgs data by using the programs HiggsBounds5.3.2 [74, 75, 76] and HiggsSignals2.2.3 [77]. These programs require as input the effective couplings of the Higgs bosons, normalized to the corresponding SM values, as well as the masses, the widths and the branching ratios of the Higgs bosons. These have been obtained for the SM and NMSSM Higgs bosons from the Fortran code NMSSMCALC [78, 13]. One of the neutral CP-even Higgs bosons is identified with the SM-like Higgs boson – it will be called $h$ from now on – and its mass is required to lie in the range

[TABLE]

For the SM input parameters we use the following values  [79, 80]

[TABLE]

Concerning the NMSSM sector, we follow the SUSY Les Houches Accord (SLHA) format [81] in which the soft SUSY breaking masses and trilinear couplings are understood as $\overline{\text{DR}}$ parameters at the scale

[TABLE]

This is also the renormalization scale that we use in the computation of the higher-order corrections. Note that we chose the charged Higgs boson mass as an OS input parameter. The computation of the ${\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2})$ corrections to the Higgs boson masses is done in the $\overline{\text{DR}}$ renormalization scheme of the top/stop sector. We have included the contribution of the gauge parameters $g_{1},g_{2}$ into the conversion from pole to $\overline{\text{DR}}$ top masses. In Table 1 we summarize the ranges applied in our parameter scan. In order to ensure perturbativity we apply the rough constraint

[TABLE]

The remaining mass parameters of the third generation sfermions that are not listed in the table are chosen as

[TABLE]

The mass parameters of the first and second generation sfermions are set to

[TABLE]

We have performed two scans. In the first (smaller) scan we took care to select only such scenarios where the lightest CP-even Higgs boson $H_{1}$ is singlet-like and the second lightest CP-even Higgs boson is the SM-like Higgs boson. We refer to this scan as scan1 in the following. In the second (larger) scan, called scan2 in the following, we only retained scenarios where the SM-like Higgs boson is the lightest CP-even Higgs boson. Both scans allow for points that have a $\chi^{2}$ computed by HiggsSignals-2.2.3 that is consistent with an SM $\chi^{2}$ within $2\sigma$. All the branching ratios shown in the following have been calculated by implementing the here presented higher-order corrections to the various decay widths in NMSSMCALC. In this way the new EW corrections are combined with the state-of-the-art higher-order QCD corrections already implemented in NMSSMCALC. Note, however, that the EW 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, are taken into account in the computation of the total decay width and branching ratios.

5.1 Decays into SM Fermion Pairs

In the old implementation in NMSSMCALC the tree-level couplings entering the various decay widths were improved by including loop effects in the Higgs mixing matrix elements. Thus, the tree-level rotation matrix ${\cal R}$ was replaced by the loop-corrected rotation matrix

[TABLE]

evaluated at zero external momentum both at one-loop and at two-loop order to ensure unitarity.262626We remind the reader, that in contrast the one-loop corrected masses are obtained at non-vanishing external momenta and the two-loop corrections at zero external momenta. The implementation here differs by the fact that in the computation of ${\bf Z}^{H}$ we include the momentum dependence at one-loop order and we do not apply the approximation of Ref. [14] to deduce ${\bf Z}^{H}$ but proceed as described in Eqs. (3.106)-(3.108). In the following, we call the couplings where we apply ${\cal R}^{l}$ as obtained from Eq. (5.224) with zero external momentum and by applying the approximation of Ref. [14] ’effective tree-level couplings’ while those with ${\bf Z}^{H}$ calculated according to Eqs. (3.106)-(3.108) including the momentum dependence at one-loop order are denominated ’improved couplings’.

The decays into SM fermion pairs in the old implementation in NMSSMCALC were calculated using the loop-corrected rotation matrix, ${\cal R}^{\text{l}}$, evaluated at zero external momenta and by including the $\Delta_{b}$ corrections272727For simplicity, we collectively call them $\Delta_{b}$ corrections although we also include the corresponding corrections in the decays into strange quarks and into leptons. into the effective tree-level couplings, as specified in Ref. [13]. The thus obtained ’effective couplings’ are given by282828We call them ’effective couplings’ and not ’effective tree-level couplings’ as they also contain the $\Delta_{b}$ corrections.

[TABLE]

Beyond the $\Delta_{b}$ approximation no further SUSY-EW nor SUSY-QCD corrections were included. To quantify the difference between the branching ratio computed in this paper and the old implementation in NMSSMCALC we introduce the relative change in the branching ratio for the decay $H_{i}\to X_{j}X_{k}$ as

[TABLE]

with $X_{j}X_{k}\equiv f\bar{f}$ for the decays into fermions. Here the branching ratio $\text{BR}_{{\bf Z}^{H}}^{\text{SEW(+SQCD)}}(H_{i}\to f\bar{f})$ means that we include the SUSY-EW corrections (and SUSY-QCD corrections for the decays into quarks) together with the wave-function renormalization factor into the decay width of the decay $H_{i}\to f\bar{f}$. The formulae are given by Eq. (4.138) for the decays into quarks and by Eq. (4.159) for the decays into leptons together with the definitions Eqs. (4.139) and (4.140). The branching ratio in the old implementation in NMSSMCALC is denoted by $\text{BR}_{{\cal R}^{l}}^{\text{tree}}(H_{i}\to f\bar{f})$ (although it also includes the $\Delta_{b}$ corrections where applicable). The quantity $\Delta_{\text{BR}}$ hence gives information on the importance of the improvement of the branching ratios by the ${\bf Z}^{H}$ factor and the SEW(+SQCD) corrections. This quantity will also be used in the investigations of the decays into gauge boson pairs and into a pair of $Z$ and Higgs bosons.

The SM-like Higgs boson is given by the $h_{u}$-like Higgs state292929As the SM-like Higgs boson has to comply with the experimentally measured Higgs rates and for small values of $\tan\beta$, as preferred by the NMSSM, is dominantly produced through gluon fusion it needs a substantial coupling to top quarks so that it is the $h_{u}$-dominated Higgs state that turns out to be SM-like., and in our scans we found valid scenarios where this can be the lightest or the second lightest of the CP-even Higgs bosons. We first consider only the parameter points where the lightest CP-even Higgs boson $H_{1}$ is the singlet-like state, i.e. has a large $h_{s}$ component. These are points obtained in the above described scan1. Here and in the following we denote a Higgs boson $H_{i}$ to be dominantly $x$-like ($x=h_{s},h_{u},h_{d},a_{s},a_{d}$) if the corresponding matrix element squared $|{\cal R}^{l}_{ix}|^{2}$ exceeds 80%. When $H_{1}$ is $h_{s}$-like, the question which final state has the largest decay width strongly depends on the amount of admixture of $h_{d}$ and $h_{u}$ components to the singlet-state. In Fig. 1 we show, for all parameter points that pass our constraints, the scatter plot of the $H_{1}$ branching ratios into SM particles against its $h_{d}$ component represented by the element ${\cal R}^{l}_{1h_{d}}$ of the loop-corrected Higgs rotation matrix303030Note, that ${\cal R}^{l}_{1h_{d}}$ is the $(1,h_{d})$-component of the mixing matrix ${\cal R}^{l}$ given by Eq. (5.224), evaluated at zero external momentum both at one- and at two-loop order and where for the computation of ${\bf Z}^{H}$ the approximation of Ref. [14] is used. In the computation of the loop-corrected branching ratios, however, we of course use the new implementation described at the beginning of this subsection.. The mass of $H_{1}$ lies between 70 and 118 GeV for these points. As can be inferred from the plot, the dominant decays are those into $b\bar{b}$, $c\bar{c}$, $\tau\bar{\tau}$ and $gg$. In most cases the branching ratio into a bottom-quark pair is dominant followed by the decay into $\tau\tau$. However, when the $h_{d}$ component of $H_{1}$ is very small the branching ratios into $gg$ and $c\bar{c}$ become competitive and can even be larger than those for the decay into $b\bar{b}$ with values beyond 60% for the $gg$ final state and of up to 35-39% for $c\bar{c}$ in some of the scenarios. In this case, i.e. for |{\cal R}^{l}_{1h_{d}}|\raisebox{-3.69899pt}{~{}\shortstack{<$\\ [-0.07cm]$\sim}}~{}0.02, also the branching ratios into $\gamma\gamma$ and into the off-shell final state $W^{+*}W^{*-}$ increase and can reach up to about 30% in the latter and about 2% in the former case. The branching ratio into the off-shell $Z^{*}Z^{*}$ final state, which also increases then, is about one order of magnitude smaller than the one into $W^{+*}W^{-*}$. But already for |{\cal R}^{l}_{1h_{d}}|\raisebox{-3.69899pt}{~{}\shortstack{>$\\ [-0.07cm]$\sim}}~{}0.02 the decay into $b\bar{b}$ takes over again and reaches branching ratio values of up to 90% followed by the branching ratio into $\tau\tau$ with values of up to 10%.

In order to investigate the importance of the higher-order corrections we define for our new implementation the relative correction of the partial width for the decay $H_{i}\to X_{j}X_{k}$ as

[TABLE]

with the higher-order decay widths for the decays $H_{i}\to q\bar{q}$ into quarks given in Eq. (4.138) and the higher-order decay widths for the decays $H_{i}\to l\bar{l}$ into leptons given in Eq. (4.159) and with the tree-level decay width $\Gamma^{\text{tree}}_{{\bf Z}^{H}}$ including only the $\Delta_{b}$ corrections. The tree-level and higher-order decay widths are both evaluated with the new implementation of ${\bf Z}^{H}$. Note that the quantity $\delta$ gives information on the importance of the SEW(+SQCD) corrections in the decay width alone as the factor ${\bf Z}^{H}$ cancels in the ratio. In Fig. 2 we show the scatter plot of the relative change of the branching ratios, $\Delta_{\text{BR}}(H_{1}f\bar{f})$, $f=b,c,\tau$, for all the parameter points passing the constraints, against $\text{BR}_{{\bf Z}^{H}}^{\text{SEW(+SQCD)}}(H_{1}\to f\bar{f})$.

The color code in Fig. 2 as well as in Figs. 3–6 denotes the sizes of the relative corrections of the partial decay widths. The points where the absolute value of $\delta^{\text{SEW(+SQCD)}}$ exceeds 10% are marked in blue, those with $\left|\delta^{\text{SEW(+SQCD)}}\right|$ in the [5,10]% range in black and those with relative corrections less than 5% in red. For Figs. 5b and 6b we distinguish two regimes for the larger corrections, in blue where $\left|\delta^{\text{SEW+SQCD}}\right|$ is in the [10,20]% range and in green where $\left|\delta^{\text{SEW+SQCD}}\right|$ is in the [20,40]% range; in Fig. 4b we also add two other categories of points, in cyan where $\left|\delta^{\text{SEW+SQCD}}\right|$ is in the [40,60]% range and in pink where $\left|\delta^{\text{SEW+SQCD}}\right|$ is in the [60,80]% range. Note that for the $\tau$ decays there are no SUSY-QCD corrections. The ballpark of the relative change $\Delta_{\text{BR}}({H_{1}b\bar{b}})$ in the branching ratios between the old and the new implementation ranges below about 30% with vertex corrections $|\delta|$ smaller than 5%. There are some very rare scenarios where $\left|\Delta({H_{1}b\bar{b}})\right|$ exceeds 50% and where at the same time the relative vertex corrections are between 5 and 10%. We investigated these cases and observed that there is an accidental cancellation either in the effective tree-level or in the improved couplings. These parameter points lead to similar results for the $\tau$ final states, i.e. $|\Delta_{\text{BR}}({H_{1}\tau\bar{\tau}})|>50$% and at the same time $|\delta|$ between 5 and 10%. The cancellation results in a suppression of the branching ratio, to less than 4% for the $\tau\bar{\tau}$ final state and 10% at maximum for the $b\bar{b}$ final state. In most of the cases, the large $\Delta_{\text{BR}}({H_{1}f\bar{f}})$ is due to the use of the wave-function renormalization factor ${\bf Z}^{H}$, however. There are also cases with a cancellation between the SUSY-EW/SUSY-QCD corrections and the wave-function renormalization factor ${\bf Z}^{H}$ correction. This results in $\Delta_{\text{BR}}({H_{1}c\bar{c}})$ being less than 1%.

We have performed the same analysis for the heavier Higgs bosons, using the full set of points from our scan2. Figure 3 is the scatter plot of the $\Delta_{\text{BR}}({hf\bar{f}})$, $f=b,\tau,c,t$, against $\text{BR}_{{\bf Z}^{H}}^{\text{SEW(+SQCD)}}(h\to f\bar{f})$ for an SM-like Higgs boson $h$, while Figures 4, 5, and 6 are the scatter plots of $\Delta_{\text{BR}}({H_{i}f\bar{f}})$, $f=b,t$, against $\text{BR}_{{\bf Z}^{H}}^{\text{SEW(+SQCD)}}(H_{i}\to f\bar{f})$ for a heavy $a_{s}$-, $a$- and $h_{d}$-like Higgs boson i.e. $H_{i}=H_{a_{s}}$, $H_{a}$, $H_{h_{d}}$, respectively. As before, the SM-like Higgs boson $h$ is always $h_{u}$-like and decays dominantly into a bottom-quark pair with a branching ratio of about 60%, as expected, followed by the decay into a $\tau$ pair and the decay into a $c$-quark pair. As can be inferred from Fig. 3, the relative changes between the old and the new implementation are much smaller than for the singlet-like lightest Higgs boson and amount only to a few percent. The relative vertex corrections $|\delta|$ are below 10% for the $b$-quark pair final state and below 5% for the decays both into $\tau\bar{\tau}$ and $c\bar{c}$.

For the heavy Higgs bosons, the decay into a top quark pair can become kinematically possible. We start by discussing the decay pattern of the heavy singlet state $H_{a_{s}}$, with a mass between 120 GeV and 1.7 TeV, into the $b$-quark and $t$-quark final states, presented in Fig. 4.313131Since we will not gain much new information, for $H_{a_{s}}$, $H_{a}$ and $H_{h_{d}}$ we do not show the corresponding plots into $\tau\bar{\tau}$ and $c\bar{c}$. For the $b\bar{b}$ final states the relative change in the branching ratios due to the new implementation is mostly between -20% and 20%. We also find points where the relative change is close to 100%, in particular for branching ratios close to 100%. Most points exhibit small relative vertex corrections (see, red points in Fig. 4a), so that the large changes of $\Delta_{\text{BR}}(H_{a_{s}}b\bar{b})$ are due to the implementation of ${\bf Z}^{H}$. This is especially the case for large branching ratios close to 100%. There are a few points where the relative vertex corrections lie between 5 and 10% (black) and even above 10% (blue). This happens for the cases where the effective tree-level couplings $H_{a_{s}}f\bar{f}$ are suppressed. The relative change $\Delta_{\text{BR}}$ can still be very small when the effects from the ${\bf Z}^{H}$ factor and the vertex corrections cancel. The decay pattern for the $t\bar{t}$ channel, finally, is displayed in Fig. 4b. The branching ratio takes all values between almost 0 and 100%. The relative changes $\Delta_{\text{BR}}(H_{a_{s}}t\bar{t})$ are mostly between -20% and 20% and close to 0% for large branching ratios above about 60%. We also observe large $\Delta_{\text{BR}}$, in particular for branching ratios close to zero. This is mainly due to very suppressed effective tree-level $H_{a_{s}}t\bar{t}$ couplings corresponding to the regions where the branching ratio $\mathrm{BR}(H_{a_{s}}\to b\bar{b})$ is enhanced. These regions correspond to large values of $\tan\beta$ close to the upper bound in our scan, or to smaller mass values of $H_{a_{s}}$ with not sufficient phase space to decay into an on-shell top-quark pair. In these regions the relative corrections $|\delta|$ are most of the time below 40%, and for cases where $|\delta|<5\%$ the large changes in $\Delta_{\text{BR}}$ are mostly due to the use of the wave-function renormalization factor ${\bf Z}^{H}$. For larger branching ratios the relative corrections $|\delta|$ are mostly below 10% (red and black points). Some rare scenarios display corrections above 40% and up to 80% (in cyan and in pink), again mostly in regions with lower branching ratios.

Similar observations can be made for the other heavier Higgs states $H_{a}$ (with a mass between 539 GeV and 2 TeV) and $H_{h_{d}}$ (with a mass between 548 GeV and 2 TeV), with the notable exception that the relative changes $\Delta_{\text{BR}}({H_{a/h_{d}}}X_{i}X_{j})$ are more reduced and never reach 100%. The relative changes $\Delta_{\mathrm{BR}}(H_{a/h_{d}}b\bar{b})$ are most of the time positive and below 40% as seen in Figs. 5a and 6a. The decays into top-quark pairs can be dominant where the decays into $b\bar{b}$ are suppressed, and the relative changes $\Delta_{\text{BR}}$ between the old and new implementation are close to zero when $\mathrm{BR}(H_{a/h_{d}}\to t\bar{t})\to 100\%$ as seen in Figs. 5b and 6b. For some rare scenarios the relative vertex corrections $|\delta|$ can reach 40%, depicted in green in the figures. Note that $\mathrm{BR}(H_{a/h_{d}}\to b\bar{b})$ can reach 90%, corresponding to regions where the effective tree-level coupling $H_{a/h_{d}}b\bar{b}$ is strongly enhanced due to large values of $\tan\beta$ while at the same time the effective tree-level coupling $H_{a/h_{d}}t\bar{t}$ is strongly suppressed.

5.2 Decays into a Massive Gauge Boson Pair

In the CP-conserving case, the heavy Higgs boson that can decay into two on-shell massive gauge bosons is $h_{d}$-like. The tree-level coupling of a Higgs boson $H_{i}$ to $VV$ ($V=W,Z$) is proportional to

[TABLE]

Due to the SM-like (i.e. $h_{u}$-like) Higgs boson coupling with almost SM-strength to the massive gauge bosons the tree-level coupling of the $h_{d}$-like heavy Higgs boson to $VV$ is almost zero because of sum rules. This leads to very suppressed tree-level partial decay widths $\Gamma(H_{h_{d}}\to VV)$.

In order to compare the results obtained in this paper with the old implementation in NMSSMCALC using the tree-level coupling together with the loop-corrected rotation matrix ${\cal R}^{l}$, we show in Fig. 7a the relative change $\Delta_{\text{BR}}(H_{h_{d}}WW)$ of the branching ratio into $WW$ between the old and the new implementation including the NLO-EW vertex corrections as described in Sec. 4.2 and the improvement with the ${\bf Z}^{H}$ factor, as a function of the loop-corrected branching ratio ${\rm BR}_{{\bf Z}^{H}}^{\rm SEW}(H_{h_{d}}\to W^{+}W^{-})$. The plotted points are those of our scan that pass the constraints we have applied. We display in Fig. 7b the same but for the decay into $ZZ$. The color and symbol code denotes the magnitude of the relative NLO electroweak vertex corrections alone, with $\delta$ defined in Eq. (5.227). As can be inferred from the plots both the ${\bf Z}^{H}$ factor and the NLO electroweak corrections can be responsible for the large relative changes in the branching ratios. This is in particular reflected by the black points for which the vertex corrections are below 20% and at the same time the relative changes $\Delta_{\text{BR}}$ can reach up to 100%. In some cases there is a cancellation between the two contributions (the vertex corrections and the ${\bf Z}^{H}$ factor) leading to relatively small relative changes in the branching ratios. These cases are the pentagon-marked full (red) points in Fig. 7a for the decay $H_{h_{d}}\to W^{+}W^{-}$, which are mostly located in regions where |\Delta_{\text{BR}}(H_{h_{d}}WW)|\raisebox{-3.69899pt}{~{}\shortstack{<$\\ [-0.07cm]$\sim}}~{}25\% while the relative vertex correction $|\delta|$ is at least 100%. In the case of the decay into a $Z$ boson pair, however, the bulk of these pentagon-marked full (red) points, indicating again a relative vertex correction $|\delta|$ of at least 100%, induces large relative changes of $\Delta_{\text{BR}}(H_{h_{d}}ZZ)$ close to 100%. These points correspond to a region which is discussed in more detail in the next paragraph.

We also note that there are two regions concentrating many points for the decay into a $Z$ boson pair, the region for which $\Delta_{\text{BR}}(H_{h_{d}}ZZ)\simeq 0$% and the one for which $\Delta_{\text{BR}}(H_{h_{d}}ZZ)\simeq 100$%. This is in contrast to the decay into a $W$ boson pair which is mostly centered around $\Delta_{\text{BR}}(H_{h_{d}}WW)\simeq 0$% for vertex corrections $|\delta|<80$% and much more scattered for the points where $80<|\delta|<100$%, displayed with cyan-pentagon-marked points (for $|\delta|\geq 100$% the above described cancellation takes place in the decay $H_{h_{d}}\to W^{+}W^{-}$). This presence of the second region in the $Z$ boson final state, for which $\Delta_{\text{BR}}$ is close to 100% can be explained by the occurrence of many parameter points having a very suppressed tree-level coupling $H_{h_{d}}ZZ$. They also correspond to regions where the loop-corrected partial decay width $\Gamma(H_{h_{d}}\to ZZ)$ is higher, up to 1 GeV, while the decay width is at most 5 MeV for the region centered around $\Delta_{\text{BR}}(H_{h_{d}}ZZ)=0$. Note, that while the tree-level couplings $H_{h_{d}}ZZ$ and $H_{h_{d}}WW$ are the same, the loop-corrected decay widths differ by the fact that the decay into $WW$ bosons receives real corrections and that in the one-loop squared contributions to the decay width $H_{h_{d}}\to W^{+}W^{-}$ we only include the (s)fermion contributions in contrast to the decay $H_{h_{d}}\to ZZ$.323232We remind the reader that we take into account this part of the two-loop corrections in case the one-loop corrected partial decay width becomes negative, see also Eq. (4.2).

5.3 Decays into a $Z$ Boson and a Higgs Boson

In the searches for heavy pseudoscalars, this decay can be an important search channel [82, 83]. We are interested here in how large the branching ratio can be and how important are the newly included higher-order corrections, in the case for which on-shell decays are possible. With the obtained valid set of parameter points, we present in Fig. 8a scatter plots of the relative changes of the branching ratios between the old and new implementation for the decay of a heavy pseudoscalar-like Higgs boson $H_{a}$ into $ZH_{1}$ and in Fig. 8b for the decay into $ZH_{2}$, against the respective loop-corrected branching ratios. The color and marker codes denote the relative sizes of the one-loop vertex corrections within specific ranges, identical to the ranges used in the previous sub-section for the decays into gauge boson pairs. The mass values of the individual involved CP-even Higgs bosons in the final state range in $123~{}\text{GeV}\leq m_{H_{1}}\leq 127~{}\text{GeV}$ and $463~{}\text{GeV}\leq m_{H_{2}}\leq 1.73~{}\text{TeV}$, while the mass of the decaying Higgs boson ranges in $539~{}\text{GeV}\leq m_{H_{a}}\leq 2.0~{}\text{TeV}$ for the on-shell decay into $ZH_{1}$ pairs and in $713~{}\text{GeV}\leq m_{H_{a}}\leq 2.0~{}\text{TeV}$ for the on-shell decay into $ZH_{2}$ pairs. As can already be inferred from the mass ranges, the $H_{1}$ state is the SM-like Higgs boson $h$, while the $H_{2}$ state is the singlet-like scalar Higgs boson $H_{h_{s}}$.

We observe that the branching ratios into the $ZH_{1}$ final state remain very small, below 0.4%, while those of the decay into $ZH_{2}$ can reach 11%. This is due to the nature of the $H_{1}$ Higgs boson that is SM-like, with very suppressed tree-level $H_{1}H_{a}Z$ couplings. The relative changes $\Delta_{\text{BR}}$ of the branching ratio for the decay $H_{a}\to ZH_{1}$ are mostly between 0 and -75% corresponding to relative vertex corrections $|\delta^{\text{SEW}}|$ being at most 60% (black, blue, and pink points), while a few points corresponding to higher vertex corrections up to more than 100% (green, cyan, and red points) can reach $\Delta_{\text{BR}}=\pm 100$%. These extreme points correspond to very small values for the branching ratios themselves which explains in turn the very large relative corrections $|\delta|$ that we observe. Note that in these decays we take into account the one-loop squared term as described in Eq. (4.191) which makes up for the main contribution to the very large relative corrections. As for the decay $H_{a}\to ZH_{2}$ the relative correction $\Delta_{\text{BR}}$ is most of the time between 0 and $\pm 25$%, corresponding to points where the relative vertex corrections $|\delta|$ are below 20%. A few points display larger $\Delta_{\text{BR}}$ values, and also larger relative vertex corrections $|\delta|$ that can reach 100% and even beyond, again for points that display very small branching ratios, below about $10^{-4}\,\%$. Note that there are points for which the correction $|\Delta_{\text{BR}}|$ is rather limited, below 25%, while the relative vertex correction can reach 40% (for one scenario even more than 100%). This can be explained by a sign compensation between the ${\bf Z}^{H}$ factor and the vertex correction.

The corresponding results for the heavy singlet-like Higgs boson $H_{h_{s}}$ decaying into $ZA_{1}$ and $ZA_{2}$ are shown in Figs. 9a and Figs. 9b, respectively. We see that the maximum achieved branching ratios for the decay $H_{h_{s}}\to ZA_{1}$ are below 20% and for a large number of parameter points are tiny. In the case of the decay $H_{h_{s}}\to ZA_{2}$ the branching ratio can reach around 15%. In most of the cases $A_{1}$ is singlet-like, corresponding to points where the branching ratio is small (below 10%), while $A_{2}$ is doublet-like. The mass values of the individual involved Higgs bosons in the final state range in $120~{}\text{GeV}\leq m_{A_{1}}\leq 1.50~{}\text{TeV}$ and $562~{}\text{GeV}\leq m_{A_{2}}\leq 1.63~{}\text{TeV}$, while the mass of the decaying Higgs boson ranges in $464~{}\text{GeV}\leq m_{H_{h_{s}}}\leq 1.75~{}\text{TeV}$ for the on-shell decay into $ZA_{1}$ pairs and in $696~{}\text{GeV}\leq m_{H_{h_{s}}}\leq 1.75~{}\text{TeV}$ for the on-shell decay into $ZH_{2}$ pairs. The cases with larger branching ratios for the decay $H_{h_{s}}\to ZA_{1}$ (larger than 10%) correspond mostly to the few $A_{1}$ pseudoscalar Higgs bosons with doublet-like admixture and mass values above $400$ GeV.

The relative changes $\Delta_{\mathrm{BR}}$ in the branching ratios are mostly positive, and can reach values of 100%. For some very rare scenarios we get $\Delta_{\mathrm{BR}}(H_{h_{s}}ZA_{1})$ close to -100%, while $\Delta_{\mathrm{BR}}(H_{h_{s}}ZA_{2})$ is not below -15%. The relative vertex corrections are moderate for the decay $H_{h_{s}}\to ZA_{1}$, mostly $|\delta|\leq 20$% (black points). This means that the large changes in $\Delta_{\mathrm{BR}}(H_{h_{s}}ZA_{1})$ are mostly due to the ${\bf Z}^{H}$ factor. For very small branching fractions below $10^{-4}\%$ larger vertex corrections are possible, mainly because the denominator in the definition of $\delta$ is very small in these regions and can lead to sharp changes in $\delta$. Note that the bulk of the changes between the old and the new implementation in these cases stems from the vertex corrections. In the case of the decay $H_{h_{s}}\to ZA_{2}$ the relative vertex corrections are mostly small, with values $|\delta|<5$% (red triangle-marked points). For large relative changes $\Delta_{\text{BR}}(H_{h_{s}}ZA_{2})$ the ${\bf Z}^{H}$ factor is responsible then.

5.4 Decays into Charginos and Neutralinos

We start by investigating the loop corrections to the masses of the charginos and neutralinos using the three renormalization schemes OS1, OS2 and $\overline{\text{DR}}$, imposed on the two gaugino masses $M_{1}$ and $M_{2}$, as defined in Section 3.1.2. According to the SLHA format that we apply in our code, $M_{1},M_{2}$ are $\overline{\text{DR}}$ parameters given at the scale $M_{\text{SUSY}}=\sqrt{m_{\tilde{Q}_{3}}m_{\tilde{t}_{R}}}$. When we use the OS schemes we have to translate the $\overline{\text{DR}}$ input parameters to the OS values by applying the approximate transformation formulae

[TABLE]

where $\delta M_{1/2}^{\text{finOSi}}$ are the finite parts of the $M_{1/2}$ counterterms computed in the OSi ($i=1,2$) renormalization scheme. Since the finite parts $\delta M_{1/2}^{\text{finOSi}}$ should be computed with OS input parameters we have used an iterative method to obtain these. For all parameter points in our scan, the size of the loop corrections to the neutralino and charginos masses, quantified by $\Delta M^{\tilde{\chi}}_{i}=M_{\tilde{\chi}_{i}}^{\text{loop}}-m_{\tilde{\chi}_{i}}^{\text{tree}}$, with $i=1,...,5$ for the neutralinos and $i=1,2$ for the charginos, never exceeds 46 GeV.

We exemplary present here a particular point, called scenario1, with the soft SUSY breaking masses and trilinear couplings given by

[TABLE]

and the remaining input parameters set to333333The imaginary part of $A_{\kappa}$ is obtained from the tadpole condition.

[TABLE]

The Higgs boson masses and their main composition in terms of singlet/doublet and scalar/pseudoscalar components at two-loop order ${\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2})$ for $\overline{\text{DR}}$ renormalization in the top/stop sector computed by NMSSMCALC, are summarized in Table 2.

For scenario1, we present in Table 3 the tree-level and loop-corrected masses of the neutralinos and charginos in the three different renormalization schemes and for the Denner description. As expected the wino-like neutralino and the wino-like chargino which couple to the electroweak gauge bosons, get significant loop corrections in the $\overline{\text{DR}}$ scheme. The one-loop corrected masses themselves, however, barely differ in the three renormalization schemes so that the remaining theoretical error due to missing higher-order corrections is very small.

We vary the phases of the gaugino mass parameters $M_{1}$ and $M_{2}$, $\varphi_{M_{1}}$ and $\varphi_{M_{2}}$, in order to study their effect on the loop-corrected neutralino and chargino masses. Note that these complex phases have negligible impact on the Higgs sector [15, 17]. We use the $\overline{\text{DR}}$ scheme for this analysis and show in the upper plots of Fig. 10 the loop corrections to the electroweakino masses,

[TABLE]

as function of $\varphi_{M_{1}}$ (left) and $\varphi_{M_{2}}$ (right). We apply a subtraction $\delta M_{x}$ to the mass corrections $\Delta M_{x}$ of the different electroweakinos that allows us to show all corrections, which can be very different in size, in one plot. In the lower plots of Fig. 10 we show the differences

[TABLE]

with $x$ denoting any electroweakino, between the loop corrected electroweakino masses computed using the Denner and the EMT descriptions presented in Section 3.1.2. As can be inferred from the upper plots, the wino-like neutralino $\tilde{\chi}_{4}^{0}$ (cyan line) and wino-like chargino $\tilde{\chi}_{2}^{+}$ (black line) receive the largest loop corrections of about 35 GeV in absolute values. The corrections to the Higgsino-like neutralinos $\tilde{\chi}_{1}^{0}$ (blue line) and $\tilde{\chi}_{2}^{0}$ (green line), the singlino-like neutralino $\tilde{\chi}_{0}^{5}$ (red line) and the Higgsino-like chargino $\tilde{\chi}_{1}^{+}$ (pink line) range around 6-8 GeV. The correction to the bino-like neutralino $\tilde{\chi}_{3}^{0}$ (orange line) is somewhat smaller with values around 1.5-1.8 GeV. A small difference between the Denner and the EMT descriptions of about 1-2 MeV is observed for the Higgsino-like neutralinos, cf. lower plots. We do not see any difference, however, for the chargino masses. This is because we have used the approximation in Eq. (3.122) for the loop-corrected chargino masses. There is a compensation between $\text{Re}\,\Sigma^{\tilde{\chi}^{+},Ls}_{ii}(m_{\tilde{\chi}_{i}^{\pm}}^{2})$ and $\text{Re}\,\Sigma^{\tilde{\chi}^{+},Rs}_{ii}(m_{\tilde{\chi}_{i}^{\pm}}^{2})$ that kills the effect of the imaginary part of the loop integral function.

In order to study the loop corrections on the decay widths, we computed the tree-level and loop-corrected decay widths, defined in Section 4.4, for the three different renormalization schemes OS1, OS2 and $\overline{\text{DR}}$ for the scenario1. Note that we use the loop-corrected masses for external Higgs bosons, charginos and neutralinos not only in the loop-corrected decay widths but also in the tree-level ones. For illustration, we present in Table 4 for the decays of all Higgs bosons in all possible electroweakino final states the loop-corrected decay widths $\Gamma^{l}\equiv\Gamma_{{\bf Z}^{H}}^{\text{SEW}}(H_{i}\to\tilde{\chi}_{j}\tilde{\chi}_{k})$, the relative corrections $\delta(H_{i}\tilde{\chi}_{j}\tilde{\chi}_{k})$ as defined in Eq. (5.227), and the loop-corrected branching ratios $\mbox{BR}^{l}\equiv\mbox{BR}_{{\bf Z}^{H}}^{\text{SEW}}(H_{i}\to\tilde{\chi}_{j}\tilde{\chi}_{k})$ using scenario1, for the OS1 and $\overline{\text{DR}}$ renormalization schemes. We found that $\Gamma^{l}$ is almost identical in the two OS schemes. The relative size of the loop corrections is below 10% in the OS scheme. The relative corrections in the $\overline{\text{DR}}$ scheme are always larger than in the OS schemes. For some channels with small decay widths, we see significant corrections in the $\overline{\text{DR}}$ scheme. For example in the decay channel $H_{4}\to\tilde{\chi}_{2}^{0}\tilde{\chi}_{2}^{0}$, the relative loop correction is $-0.13\%$ in the OS scheme while it is $-20\%$ in the $\overline{\text{DR}}$ scheme. Based on our investigation, we conclude that it is better to use the OS scheme in the decays of the neutral Higgs bosons into electroweakinos. The largest uncertainty due to missing higher-order corrections that we estimate from the variation of the renormalization schemes is found to be 17% in the decay $H_{5}\to\tilde{\chi}_{1}^{0}\tilde{\chi}_{2}^{0}$. We also studied the difference between the Denner and EMT descriptions and did not observe any significant difference. Defining the difference as $(\Gamma^{\text{D}}_{x}-\Gamma^{\text{EMT}}_{x})/\Gamma^{\text{D}}_{x}$ with $\Gamma_{x}$ being the loop-corrected decay width for some Higgs decay into an electroweakino final state, we see that the differences are of per mille level for all investigated decays.

5.5 Decays into a Squark Pair

We start by discussing the top and bottom squark masses in the OS and $\overline{\text{DR}}$ renormalization schemes defined in Section 3.1.3. We follow the SLHA convention where the input parameters $m_{\tilde{Q}_{3}},m_{\tilde{t}_{R}},m_{\tilde{b}_{R}},A_{t},A_{b}$ are $\overline{\text{DR}}$ parameters at the scale $M_{\text{SUSY}}$. When we apply the OS scheme, these parameters must be translated into OS parameters by applying the approximate transformation formula ($i=1,2$)

[TABLE]

with $X=m_{\tilde{Q}_{3}},m_{\tilde{t}_{R}},m_{\tilde{b}_{R}},A_{t},A_{b}$ and the finite part of their OS counterterms denoted as $\delta X^{\text{finOSi}}$. We have used an iterative method to obtain a stable value of $\delta X^{\text{finOSi}}$. Note that we include both the NLO QCD and the full NLO EW contribution in the conversion Eq. (5.234).343434This is a bit different from the Higgs mass calculation in NMSSMCALC where we include the NLO QCD correction and the NLO Yukawa correction of order ${\cal O}(\alpha_{t})$ to the conversion in the OS renormalization scheme of the top/stop sector. Using these OS parameters together with the OS top mass in the tree-level mass matrices, we obtain the top and bottom squark masses. When we apply the $\overline{\text{DR}}$ renormalization scheme, the top pole mass has to be translated to the $\overline{\text{DR}}$ top mass for which we follow the description in appendix C of Ref. [17]. The $\overline{\text{DR}}$ top and bottom masses together with the $\overline{\text{DR}}$ parameters $m_{\tilde{Q}_{3}},m_{\tilde{t}_{R}},m_{\tilde{b}_{R}},A_{t},A_{b}$ are then used in the tree-level mass matrices to get the tree-level rotation matrices. They are subsequently used in the computation of the renormalized self-energies of the top and bottom squarks to obtain the loop-corrected squark masses as described in Section 3.4. In principle, we would expect that the loop-corrected masses computed in the $\overline{\text{DR}}$ scheme are closer to the OS masses in the OS description if one includes more higher order corrections. We consider here a parameter point (scenario2) given by the following soft SUSY breaking masses and trilinear couplings353535This parameter point is allowed by HiggsBounds5.3.2 and its $\chi^{2}$ computed by HiggsSignals2.2.3 is consistent with an SM $\chi^{2}$ less than 1$\sigma$.

[TABLE]

with the CP phases given by

[TABLE]

The remaining input parameters have been set to

[TABLE]

With the given $\overline{\text{DR}}$ parameters of the squark sector, we obtain their corresponding OS parameters

[TABLE]

The tree-level and loop-corrected masses of the stops and sbottoms in the OS and $\overline{\text{DR}}$ scheme are shown in Table 5. We see that for the $\overline{\text{DR}}$ scheme there are large changes between the tree-level and one-loop masses, in particular for the lightest stop $\tilde{t}_{1}$. The loop-corrected masses, however, are then closer to each other in both schemes, as expected. The maximum difference if found for the light sbottom $\tilde{b}_{1}$ mass, where the OS and $\overline{\text{DR}}$ results differ by 27 GeV at one-loop order (compared to 59 GeV at tree level). The two-loop corrected neutral Higgs boson masses at ${\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2})$ together with their respective main component are displayed in Table 6. The SM-like Higgs boson mass is around 124 GeV while the remaining Higgs spectrum is quite heavy with masses above 1.6 TeV.

We now turn on the complex phase of $A_{t}$. In the left plot of Fig. 11, we show the tree-level (black), NLO EW (green), NLO QCD (blue), and full, i.e. NLO QCD+EW, (red) corrections to the partial width of the decay $H_{4}\to\tilde{t}_{1}\tilde{t}_{1}^{*}$ as function of the phase $\varphi_{A_{t}}$ in both the OS (full lines) and the $\overline{\text{DR}}$ (dashed lines) schemes while their corresponding branching ratios are depicted in the right plot. The decay $H_{4}\to\tilde{t}_{1}\tilde{t}_{1}^{*}$ vanishes in the CP-conserving limit where $H_{4}$, which is $a$-like in scenario2, is a CP-odd Higgs boson. (Note that CP-odd Higgs bosons at tree-level only couple to two different stops.) In the OS scheme, the relative EW corrections $\delta$ (see Eq. (5.227) for the definition) vary in the range (-6%,-10 %) and the QCD corrections in the range $(-4\%,-8\%)$ depending on the phase $\varphi_{A_{t}}$ that is varied from zero to $\pm\pi/2$. In the $\overline{\text{DR}}$ scheme, the relative EW corrections are of order $-10\%$ and the relative QCD corrections of $22\%$ and depend slightly on the phase $\varphi_{A_{t}}$. We define the relative differences between the OS and $\overline{\text{DR}}$ decay widths and branching ratios, as

[TABLE]

and

[TABLE]

respectively. For $\varphi_{A_{t}}=-\pi/2$, the relative difference $\Delta_{\Gamma}$ of the partial decay widths between the OS and $\overline{\text{DR}}$ schemes is then about 40% at tree-level while it reduces to $4\%$ when both QCD and EW corrections are included, so that at one-loop level we clearly see a reduction of the theoretical error due to missing higher-order corrections For the relative error in the branching ratios we find values between 32% and 27% at tree level and between 0.3% and 3% at one-loop order including both the EW and QCD corrections while the phase $\varphi_{A_{t}}$ is varied from $\pm\pi/2$ to zero.

For the decay $H_{4}(\equiv H_{a})\to\tilde{t}_{1}\tilde{t}_{2}^{*}$, we show in the upper panels of Fig. 12 the partial decay width (left) and branching ratio (right) at tree-level (black), NLO EW (green), NLO QCD (blue), and NLO EW+QCD (red) as a function of $\varphi_{A_{t}}$, both for the OS (full) and $\overline{\text{DR}}$ (dashed) scheme . In the middle panels we show the relative NLO EW, NLO QCD and NLO EW+QCD corrections which are defined as

[TABLE]

and

[TABLE]

respectively. The lower panels display the relative differences between the OS and $\overline{\text{DR}}$ decay widths and branching ratios, $\Delta_{\Gamma}$ and $\Delta_{\text{BR}}$, as defined in Eq. (5.239) and Eq. (5.240), respectively. The corrections vary slightly with the phase $\varphi_{A_{t}}$. The EW corrections are negative in both schemes while the QCD corrections are positive and of the same order of magnitude. This shows the importance to include both types of corrections to make reliable predictions. Overall, the relative corrections $\delta$ in the $\overline{\text{DR}}$ scheme are larger than in the OS scheme. As can be inferred from the bottom left panel of Fig. 12, for $\varphi_{A_{t}}=0$, the relative difference in the partial width, $\Delta_{\Gamma}$, between the OS and $\overline{\text{DR}}$ scheme at tree level is about $37\%$ and decreases dramatically to less than $12\%$ when both the EW and the QCD corrections are included.363636When we only include the EW corrections the scheme dependence increases when going from tree- to one-loop level. Overall, the behavior is as expected, however, when the full set of corrections is included. This shows that care has to be taken, when estimating the theoretical uncertainty due to missing higher-order corrections based on a change of the renormalization schemes, if not all corrections of a given loop order are included. See also Ref. [17] for a similar discussion. Similar results are found for the branching ratios, presented on the right plots of Fig. 12, with smaller values of $21\%$ at tree level and $7\%$ at full one-loop order. Note that in the right hand side plots we treated the decays $H_{a}\to\tilde{t}_{1}\tilde{t}_{2}^{*}$ and $H_{a}\to\tilde{t}_{2}\tilde{t}_{1}^{*}$ at the same level of precision while all other decays are computed at the highest possible precision.

In the CP-invariant scenario the decay width of decay $H_{a}\to\tilde{t}_{1}\tilde{t}_{2}^{*}$ is equal to the one of its charge conjugate decay $H_{a}\to\tilde{t}_{2}\tilde{t}_{1}^{*}$. For non-vanishing $\varphi_{A_{t}}$, however, the CP asymmetry, defined as

[TABLE]

is non-zero. In Fig. 13 we show the CP asymmetry as a function of $\varphi_{A_{t}}$. We see that the CP asymmetry appears already at tree-level, which results from the imaginary part of the WFR factor ${\bf Z}^{H}$ and the imaginary part of the tree-level couplings $g_{h_{i}\tilde{q}_{j}\tilde{q}_{k}^{*}}$. The relative change of the asymmetry due to loop corrections ranges between 18% and -9% in the OS scheme while in the $\overline{\text{DR}}$ scheme it is about 8% when the phase $\varphi_{A_{t}}$ is varied from $-\pi/4$ to $\pi/4$.

We present in Fig. 14 the same plots for the decay widths and branching ratios as in Fig. 12 but for the decay of $H_{3}\to\tilde{t}_{1}\tilde{t}_{2}^{*}$ which is the dominant decay channel of $H_{3}$. In scenario2 $H_{3}$ is $h_{d}$-like. For both the OS and the $\overline{\text{DR}}$ scheme the NLO EW corrections to the decay width are negative and the relative corrections $\delta_{\Gamma}$ are around -3% in the OS and -8% in the $\overline{\text{DR}}$ scheme. The NLO QCD corrections on the other hand are positive and their relative size can reach $8\%$ in the OS scheme and around $37\%$ in the $\overline{\text{DR}}$ scheme. For $\varphi_{A_{t}}=0$, the difference between the decay widths in the OS and the $\overline{\text{DR}}$ scheme, $\Delta_{\Gamma}$, is about $36\%$ at tree-level and reduces to 12% at full one-loop level. The corresponding numbers for the branching ratios are similar.

We have observed that for this parameter point the EW and QCD corrections in the OS scheme are smaller in the $\overline{\text{DR}}$ scheme. In the OS scheme we have seen that this is due to a cancellation between the genuine triangle diagram contribution and the counterterm contribution while this does not happen in the $\overline{\text{DR}}$ scheme. Overall the inclusion of the QCD corrections in addition to the EW corrections reduces the difference between the OS and $\overline{\text{DR}}$ results.

6 Conclusions

In this paper, we have presented in the framework of the CP-violating NMSSM our calculation of the complete (SUSY-)EW to the on-shell Higgs boson decays into fermion pairs, massive gauge boson final states, gauge and Higgs boson final state pairs, electroweakino and stop and sbottom pairs. Where applicable we have included the SUSY-QCD corrections. We have implemented these new corrections into NMSSMCALC, a Fortran code for the computation of the Higgs mass spectrum up to presently two-loop order ${\cal O}(\alpha_{t}\alpha_{s}+\alpha_{t}^{2})$ and the calculation of their branching ratios. The code already included in the branching ratios the state-of-the-art QCD corrections and the $\Delta_{b}$ corrections as well as decays into off-shell massive gauge boson pairs and decays with off-shell heavy quarks in the final state. The new code is called NMSSMCALCEW.

The consistent implementation of our newly computed corrections provides the presently highest level of precision in the calculation of the NMSSM Higgs boson decays. In contrast to the previous NMSSMCALC version, we have included the ${\bf Z}^{H}$ factor with full momentum dependence in order to render the loop-corrected masses on-shell. For the decays into the electroweakinos and into squark pairs different renormalization schemes were implemented. The numerical analysis has demonstrated that the relative change in the branching ratios due to this new treatment and the newly implemented corrections is significant. The analysis of the decays into chargino/neutralino and squark pairs for different renormalization schemes has shown that the one-loop corrections reduce the theoretical uncertainty due to missing higher-order corrections. The new code NMSSMCALCEW can be obtained at the url:

https://www.itp.kit.edu/$\sim$maggie/NMSSMCALCEW/.

Appendix A Counterterm contribution to the decay of a neutral Higgs boson into a squark pair

In this appendix, we give explicit expressions of the counterterm couplings for the decays of a neutral Higgs boson into a squark pair. Note that we have used the redefined WFR constants for the squarks given in Eqs. (3.125) and (3.126). For the EW corrections, the counterterm entering Eq. (4.216) for the Higgs decay into a stop pair is given by

[TABLE]

[TABLE]

and for the decay into a sbottom pair

[TABLE]

where

[TABLE]

For the QCD corrections, the counterterm couplings are obtained from ${\cal M}^{\text{CT,EW}}_{h_{i}\tilde{q}_{j}\tilde{q}_{k}^{*}}$ by setting the set of EW counterterms, $\delta Z_{e},\delta M_{Z}^{2},\delta v,\delta v_{s},\delta\lambda,\delta s_{\theta_{W}},\delta c_{\theta_{W}},\delta{s_{\beta}},\delta{c_{\beta}}$ and $\delta Z_{h_{i}h_{i^{\prime}}}$ to zero. Note that the counterterm $\delta\mu_{\text{eff}}$ is given in terms of $\delta\lambda$ and $\delta v_{s}$ as

[TABLE]

Appendix B The Code NMSSMCALCEW

We here give a brief introduction into our new code NMSSMCALCEW that includes the newly calculated and here presented corrections to the decay widths of the neutral NMSSM Higgs bosons in the CP-violating NMSSM, as well as the newly calculated and here presented one-loop corrections to the chargino, neutralino, stop, and sbottom masses. It is derived from the code NMSSMCALC, which is described in detail in Ref. [13]. We here concentrate on the new features in NMSSMCALCEW with respect to NMSSMCALC.

NMSSMCALCEW requires ’LoopTools’ version 2.14 (or higher) [84, 85] to work with the EW corrections in the decay part. If ’LoopTools’ is not installed yet, it can be obtained from the url: http://www.feynarts.de/looptools/

In order to generate the executable, download and unpack the tar file ’nmssmcalcew.tar.gz’. It contains two subdirectories called ’nmssmcalc_rew_alphat2-master’ and ’nmssmcalc_cew_alphat2-master’ for the real and complex NMSSM, respectively. Go to the subdirectory of the version in which you want to work in. Open in a text editor the file ’makefile’ and in line 31 provide the absolute path to the ’LoopTools’ binary directory located in the main directory of LoopTools. Modify also the line 66 (to make sure it refers to the correct ’lib’ sub-directory within the ’LoopTools’ binary directory). In case the package is compiled without the EW corrections in the decay widths, the flag ’yesEW’ can be switched to ’FALSE’ on line 19 and ’LoopTools’ is not needed anymore. Subsequently, all files are compiled by typing ’make’. An executable ’run’ is created. By typing’./run’ the executable is run.

For the code to be run, the user has to provide the input files for ’CalcMasses.F’ (default name ’inp.dat’) and for bhdecay(_c).f (to be named ’bhdecay.in’). The user also has the choice to provide in the command line the names of the input and output files for ’CalcMasses.F’ (first and second argument) and the name of the output file provided by the decay routine (third argument). Hence the command will be ’run name_file1 name_file2 name_file3’ in this case. Sample input files ’inp.dat’ and ’bhdecay.in’ are included in the .tar files. By typing ’make clean’ the executable as well as the object files generated in the ’obj’ directory are removed.

In ’bhdecay.in’ that is used by NMSSMCALCEW, new options have been included. They are

• •

’ischhXX’ to choose the renormalization scheme for the loop corrected electroweakino masses. The options are 1 (OS1), 2 (OS2) and 3 (DRbar). The two OS schemes are specified in Section 3.1.2.

• •

’ischhst’ to choose the renormalization scheme for the stop sector. The options are 1 (OS NLO-EW, OS NLO-QCD), 2 (OS NLO-EW, DRbar NLO-QCD), 3 (DRbar NLO-EW, DRbar NLO-QCD).

• •

’ischhsb’ to choose the renormalization scheme for the sbottom sector. The options are 1 (OS NLO-EW, OS NLO-QCD), 2 (OS NLO-EW, DRbar NLO-QCD), 3 (DRbar NLO-EW, DRbar NLO-QCD).

• •

’iewh’ to choose the levels of NLO SUSY-EW (SUSY-QCD) corrections that are included. The options are 0 (as in NMSSMCALC 3.00), 1 (decays as implemented in NMSSMCALC, but with the ${\bf Z}^{H}$ factor), 2 (full NLO corrections as described in this paper). Both for option 1 and 2, the loop-corrected electroweakino and stop/sbottom masses are used in the phase space factor.

Further information on the organization of the files for the code and their functionalities as well as modifications on the code (which are constantly updated) can be found at the webpage of NMSSMCALCEW. The code has been tested on a Linux machine.

Acknowledgments

JB acknowledges the support from the Carl-Zeiss foundation. TND and MM acknowledge financial support from the DFG project Precision Calculations in the Higgs Sector - Paving the Way to the New Physics Landscape (ID: MU 3138/1-1). TND’s work is funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.01-2017.78.

Note added While this paper was being completed, Ref. [86] appeared, where the impact of the EW corrections to heavy Higgs boson decays into fermion pairs was discussed. It was shown that these contributions are dominated by Sudakov double logarithms. The authors furthermore stressed the relevance of three-body decays for a consistent evaluation of the total widths and branching ratios at complete one-loop order.