H-COUP Version 3: A program for one-loop corrected decays of any Higgs bosonsin non-minimal Higgs models

The H-COUP program is provided as a package of Fortran codes, which can compute observables related to Higgs bosons including radiative corrections in various extended Higgs sectors. We give a manual for the latest version of H-COUP (H-COUP_3.0), in which decay rates and branching ratios of all the Higgs bosons can be calculated at one-loop level in EW and Higgs interactions with QCD corrections in the Higgs singlet model, four types of the two Higgs doublet model with a softly-broken Z2 symmetry, and the inert doublet model. The previous version (H-COUP_2.0) can evaluate those only for the standard model like Higgs boson with the mass of 125 GeV (h). In H-COUP_3.0, renormalized quantities are computed based on the gauge-independent on-shell renormalization scheme. The source code of H-COUP_3.0 can be downloaded via the following link: http://www-het.phys.sci.osaka-u.ac.jp/~hcoup. By using H-COUP_3.0, we can compare the precise measurements of the properties of h and direct searches for additional Higgs bosons with their predictions at one-loop level, by which we can reconstruct the structure of the Higgs sector. New version program summary Program Title: H-COUP Version 3 5CPC Library link to program files: https://doi.org/10.17632/f88szmrj5x.3Developer's repository link: http://www-het.phys.sci.osaka-u.ac.jp/~hcoupLicensing provisions: GPLv3Programming language: Fortran90Journal reference of previous version: Comput. Phys. Commum. 257 (2020) 107512Does the new version supersede the previous version?: YesReasons for the new version: The previous version (H-COUP version 2) is a package of Fortran codes to provide one-loop corrected couplings, decay rates and branching ratios of the 125 GeV Higgs boson in the Higgs Singlet Model (HSM), the Two Higgs doublet models (THDMs) with a softly-broken Z2 symmetry, and the Inert Doublet Model (IDM). In the new version (H-COUP version 3), the H-COUP program has been improved to calculate the decay rates and the branching ratios of all the Higgs bosons including one-loop corrections in the above models. These calculations are performed based on the gauge independent on-shell renormalization scheme.Summary of revisions: We have added radiatively corrected decay rates and branching ratios for all the additional Higgs bosons, not only for the 125 GeV Higgs boson, in the extended Higgs models. Furthermore, several choices are provided for renormalization schemes used to evaluate electroweak corrections. The scalar two-point functions are renormalized by on-shell scheme as in the previous version, while the renormalization for the tadpole is performed by either the standard tadpole scheme or the alternative tadpole scheme. The users can choose these renormalization schemes.Nature of problem: Decay rates and branching ratios of all the Higgs bosons can be calculated with one-loop electroweak corrections and QCD corrections in the HSM, four types of the THDMs with the softly-broken Z2 symmetry, and the IDM. In the calculations of additional Higgs boson decays, for the three-body decay processes and loop induce processes of additional Higgs bosons, electroweak corrections are not included. As for the QCD correction, next-to-leading order or next-to-next-to-leading order corrections are included, depending on the decay process.Solution method: Electroweak corrections are evaluated based on the gauge independent on-shell renormalization scheme, while some parameters are renormalized by the MS‾ scheme. QCD corrections are evaluated with the MS‾ scheme.Additional comments including restrictions and unusual features: All functions up to the previous version are included in H-COUP version 3. References [1]S. Kanemura, M. Kikuchi, K. Mawatari, K. Sakurai, K. Yagyu, Comput. Phys. Commun. 257 (2020) 107512 https://doi.org/10.1016/j.cpc.2020.107512 [arXiv:1910.12769 [hep-ph]].[2]S. Kanemura, M. Kikuchi, K. Sakurai, K. Yagyu, Comput. Phys. Commun. 233 (2018) 134–144, https://doi.org/10.1016/j.cpc.2018.06.012 [arXiv:1710.04603 [hep-ph]].[3]M. Aiko, S. Kanemura, K. Sakurai, Nucl. Phys. B 986 (2023) 116047 https://doi.org/10.1016/j.nuclphysb.2022.116047 [arXiv:2207.01032 [hep-ph]].[4]S. Kanemura, M. Kikuchi, K. Yagyu, Nucl. Phys. B 983 (2022) 115906, https://doi.org/10.1016/j.nuclphysb.2022.115906 [arXiv:2203.08337 [hep-ph]].[5]M. Aiko, S. Kanemura, K. Sakurai, Nucl. Phys. B 973 (2021) 115581, https://doi.org/10.1016/j.nuclphysb.2021.115581 [arXiv:2108.11868 [hep-ph]].[6]S. Kanemura, M. Kikuchi, K. Mawatari, K. Sakurai, K. Yagyu, Nucl. Phys. B 949 (2019) 114791, https://doi.org/10.1016/j.nuclphysb.2019.114791 [arXiv:1906.10070 [hep-ph]].[7]S. Kanemura, M. Kikuchi, K. Mawatari, K. Sakurai, K. Yagyu, Phys. Lett. B 783 (2018) 140–149, https://doi.org/10.1016/j.physletb.2018.06.035 [arXiv:1803.01456 [hep-ph]].[8]S. Kanemura, M. Kikuchi, K. Sakurai, K. Yagyu, Phys. Rev. D 96 (2017) (3) 035014, https://doi.org/10.1103/PhysRevD.96.035014 [arXiv:1705.05399 [hep-ph]].[9]S. Kanemura, M. Kikuchi, K. Yagyu, Nucl. Phys. B 917 (2017) 154–177, https://doi.org/10.1016/j.nuclphysb.2017.02.004 [arXiv:1608.01582 [hep-ph]].[10]S. Kanemura, M. Kikuchi, K. Sakurai, Phys. Rev. D 94 (2016) (11), 115011 https://doi.org/10.1103/PhysRevD.94.115011 [arXiv:1605.08520 [hep-ph]].[11]S. Kanemura, M. Kikuchi, K. Yagyu, Nucl. Phys. B 907 (2016) 286–322, https://doi.org/10.1016/j.nuclphysb.2016.04.005 [arXiv:1511.06211 [hep-ph]].[12]S. Kanemura, M. Kikuchi, K. Yagyu, Nucl. Phys. B 896 (2015), 80–137, https://doi.org/10.1016/j.nuclphysb.2015.04.015 [arXiv:1502.07716 [hep-ph]].[13]S. Kanemura, M. Kikuchi, K. Yagyu, Phys. Lett. B 731 (2014) 27–35, https://doi.org/10.1016/j.physletb.2014.02.022 [arXiv:1401.0515 [hep-ph]].[14]S. Kanemura, Y. Okada, E. Senaha, C. P. Yuan, Phys. Rev. D 70 (2004) 115002, https://doi.org/10.1103/PhysRevD.70.115002 [arXiv:hep-ph/0408364 [hep-ph]].