ewN2HDECAY-A program for the calculation of Electroweak one-loop corrections to Higgs decays in the Next-to-Minimal Two-Higgs-Doublet Model including state-of-the-art QCD corrections

We present in this paper our new program package ewN2HDECAY for the calculation of the partial decay widths and branching ratios of the Higgs bosons of the Next-to-Minimal 2-Higgs Doublet Model (N2HDM). The N2HDM is based on a general CP-conserving 2HDM which is extended by a real scalar singlet field. The program computes the complete electroweak one-loop corrections to all non-loop-induced two-body on-shell Higgs boson decays in the N2HDM and combines them with the state-of-the-art QCD corrections that are already implemented in the existing program N2HDECAY. Most of the independent input parameters of the electroweak sector of the N2HDM are renormalized in an on-shell scheme. The soft-Z2-breaking squared mass scale m2 and the vacuum expectation value vs of the SU(2)L singlet field, however, are renormalized with MS conditions, while for the four scalar mixing angles cei (i = 1, 2, 3) and p of the N2HDM, several different renormalization schemes are applied. By giving out the leading-order and the loop-corrected partial decay widths separately from the branching ratios, the program ewN2HDECAY not only allows for phenomenological analyses of the N2HDM at highest precision, it can also be used for a study of the impact of the electroweak corrections and the remaining theoretical uncertainty due to missing higher-order corrections based on a change of the renormalization scheme. The input parameters are then consistently calculated with a parameter conversion routine when switching from one renormalization scheme to the other. The latest version of the program ewN2HDECAY can be downloaded from the URL https://github.com/ marcel-krauseiewN2HDECAY. Program summary Program Title: ewN2HDECAY Program Files doi: http://dx.doLorg/10.17632/g9wyfrpy68.1 Licensing provisions: GPLv3 Programming language: FORTRAN 77/90; Python 2/3 Nature of problem: Numerical calculation of branching ratios and partial decay widths for the decays of all Higgs bosons in the Next-to-Minimal Two-Higgs-Doublet Model at electroweak one-loop level for non-loop-induced on-shell decays in combination with contributions from the state-of-the-art QCD corrections and off-shell decays already implemented in HDECAY 6.51. Renormalization of the electroweak sector in 10 different renormalization schemes and automatic input parameter conversion for the comparison of numerical results calculated within these different schemes. Solution method: All tree-level and one-loop contributions to the genuine decay amplitudes and for the definition of the renormalization constants are calculated fully analytically and evaluated numerically. The electroweak corrections are consistently combined with the existing corrections implemented in HDECAY 6.51. Additional comments including restrictions and unusual features: For the electroweak corrections, only non-loop-induced on-shell decays are implemented. (C) 2019 Elsevier B.V. All rights reserved.