CN2HDM_HDECAY

Released by Margarete Mühlleitner, Jonas Müller, Jonas Wittbrodt, Sophie L. Williamson



A Fortran program for the calculation of the decay widths and branching ratios of the Higgs bosons in the CN2HDM, i.e. the CP-violating 2HDM extended by a complex singlet field.



Description:

The Fortran program CN2HDM_HDECAY calculates the decays of the Higgs bosons in the CN2HDM. The CN2HDM potential is CP-violating and has two Z_2 symmetries one of which is softly broken while the other remains unbroken. The real part of the neutral scalar fields of both the two doublet and the singlet fields are allowed to acquire vacuum expectation values (VEVs), and the then remaining symmetry stabilizes the complex component field of the singlet that becomes the Dark Matter (DM) candidate. The model thus features not only CP violation required for successful baryogenesis but also a DM candidate. The Higgs sector consists of four CP-mixing visible neutral Higgs bosons, a DM candidate and a pair of oppositely charged Higgs bosons. The program includes in the decay widths the most important higher-order QCD effects and consistently neglects all EW corrections. They are not available for the CN2HDM and contrary to the QCD corrections, cannot be taken over from the SM.

A detailed description of the model and the program can be found in M. Mühlleitner, J. Müller, J. Wittbrodt, S.L. Willliamson, arXiv:yymm.wxyz. .


Files:

CN2HDECAY consists of the main program file cn2hdm_hdecay.f. Further linked routines are dmb.f, elw.f, feynhiggs.f, haber.f, hgaga.f, hgg.f, hsqsq.f, susylha.f. The input file is named cn2hdm_hdecay.in. It is the original hdecay.in input file supplemented by additional input values.

When setting in cn2hdm_hdecay.in 2HDM=1 and CN2HDM=1 then the CN2HDM decay widths and branching ratios are calculated. The meaning of the remaining input values is explained in the header of the file cn2hdm_hdecay.f. The output files that are generated are called br.hxy_CN2HDM for the CP-mixing visible neutral Higgs bosons and br.cz_CN2HDM for the charged Higgs bosons. Here x=1,2,3,4 denotes the Higgs boson, i.e. H1, H2, H3, H4, and the files with y=1 include the fermionic decays, with y=2 the decays into gauge bosons and with y=3,4,5 the Higgs-to-Higgs decays, the Higgs-to-gauge-Higgs decays and the total width. For the charged Higgs bosons, the files with z=1,2 include the fermionic decays, and the file with z=3 the decays into a charged W boson plus Higgs boson pair. Note that the Standard Model decay widths and branching ratios are obtained by setting COUPVAR=0, 2HDM=0 and CN2HDM=0 in the input file.