The Fortran code NMHDECAY [1,2] computes the masses, couplings and decay widths of all Higgs bosons as well as the masses of all sparticles of the NMSSM in terms of its parameters at the electroweak or susy breaking scale. The computation of the Higgs spectrum includes leading electroweak corrections, two loop terms and propagator corrections as in [3]. (We are grateful to P. Slavich for providing the corresponding routines.) The computation of the decay widths is carried out as in HDECAY [4]. The code is modified in the case of light Higgs bosons in order to take into account interactions with the strong-coupling sector beyond the partonic approximation. In the chiral limit, the decays of a CP-even Higgs are implemented following the description in [5], and those of a CP-odd Higgs follow [6]. At quark thresholds, Higgs mixings with the CP-even or CP-odd quarkonia are considered as described in [7,8,9]. Each point in parameter space is checked against negative Higgs boson searches at LEP, and negative sparticle searches at LEP and the Tevatron, including unconventional channels relevant for the NMSSM. LHC constraints on the couplings of the observed Higgs boson with mass ~125 GeV are checked and updated regularly (cf. History). B and K physics constraints are included as in ref. [10,11] and also updated occasionally. Limits from CLEO and BABAR on radiative decays into leptons, quarks or invisible final states are also considered [12,13]. A calculation of the muon anomalous magnetic moment in the NMSSM is also provided, following [14]. The dark matter relic density, direct and indirect detection cross sections can be computed via a link to a NMSSM version of the micrOMEGAs package [15]. Constraints on dark matter direct detection (both spin dependent and spin independent) are also checked with up to date experimental data. SLHA conventions [16] for input and output are used. NMHDECAY is included in the NMSSMTools package that can be downloaded. Sample input and output files are included.

* Laboratoire de Physique Theorique, Universite de Paris-Sud and Univ. Paris-Saclay, F-91405 Orsay, France.

** Davis Institute for High Energy Physics, University of California, Davis, California 95616, USA.

*** Laboratoire Univers et Particules de Montpellier, Universite de Montpellier, F-34095 Montpellier, France.

The Fortran code NMSPEC [17] proceeds like NMHDECAY, but with soft Susy breaking terms specified at the GUT scale. The soft Susy breaking terms at the GUT scale can be chosen as non-universal, if desired. SLHA conventions [16] for input and output are used. NMSPEC is included in the NMSSMTools package that can be downloaded. Sample input and output files are included.

* Laboratoire de Physique Theorique, Universite de Paris-Sud and Univ. Paris-Saclay, F-91405 Orsay, France.

** Laboratoire Univers et Particules de Montpellier, Universite de Montpellier, F-34095 Montpellier, France.

The Fortran code NMGMSB allows for gauge mediated supersymmetry breaking terms, including couplings of the singlet to messengers. The user can either assume general boundary conditions at the messenger scale, as described in [18], or boundary conditions at the GUT scale following [19,20]. The gravitino mass is estimated and, if it is the LSP, micrOMEGAs is used to compute the NLSP relic density, which is then rescaled by mass(gravitino)/mass(NLSP). The dark matter direct and indirect detection cross sections are not computed. SLHA conventions [16] for input and output are used. NMGMSB is included in the NMSSMTools package that can be downloaded. Sample input and output files are included.

* Laboratoire de Physique Theorique, Universite de Paris-Sud and Univ. Paris-Saclay, F-91405 Orsay, France.

** Laboratoire Univers et Particules de Montpellier, Universite de Montpellier, F-34095 Montpellier, France.

The Fortran code NMHDECAY_CPV allows to treat the CP-violating NMSSM allowing for complex input parameters at the SUSY scale following [21]. The same collider limits as in the CP-conserving case are included, with the addition of limits from electric dipole moments. The dark matter relic density, direct and indirect detection cross sections are not computed. SLHA conventions [16] for input and output are used. NMHDECAY_CPV is included in the NMSSMTools package that can be downloaded. Sample input and output files are included.

* Instituto de Fisica Teorica (UAM-CSIC), E-28049 Madrid, Spain.

** Laboratoire Univers et Particules de Montpellier, Universite de Montpellier, F-34095 Montpellier, France.

The Fortran code NMSDECAY [22] allows to compute sparticle widths and branching ratios. It is based on a generalization of SDECAY [23], including the corresponding QCD corrections and 3-body decay modes. Slepton 3-body decays, possibly relevant in case of a singlino-like LSP, have been added. Decays to gravitinos are included for GMSB models. The output follows SLHA conventions [16]. NMSDECAY is included in the NMSSMTools package that can be downloaded.

* Laboratoire de Physique Theorique, Universite de Paris-Sud and Univ. Paris-Saclay, F-91405 Orsay, France.

** Laboratoire de Physique Corpusculaire, CNRS/IN2P3, F-63171 Aubiere, France.

*** Laboratoire Univers et Particules de Montpellier, Universite de Montpellier, F-34095 Montpellier, France.

The code micrOMEGAS [15,24,25,26], which explores the properties of dark matter in generic models, is included in the NMSSMTools package. It computes the relic density of the LSP in the NMSSM. In GMSB models, where the gravitino can be the true LSP, micrOMEGAs computes the NLSP relic density, which is then rescaled as explained above. All annihilation, coannihilation and semi-annihilation channels are included. The code also computes several observables such as the cross-sections for both spin dependent and spin independent interactions of the (neutralino) LSP on protons as well as the rates for its scattering on nuclei. Version 5.0 allows to compute the relic density of feebly interacting dark matter candidates via the freeze-in mechanism [27].

* Laboratoire d'Annecy-le-Vieux de Physique Théorique, Université de Savoie, F-74941 Annecy-le-Vieux, France.

** Skobeltsyn Institute of Nuclear Physics, Moscow State University, Moscow 119992, Russia.

*** Joint Institute of Nuclear research, JINR, 141980 Dubna, Russia.

- [1] U. Ellwanger, J. F. Gunion and C. Hugonie, JHEP 0502 (2005) 006, arXiv:hep-ph/0406215

- [2] U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 175 (2006) 290-303, arXiv:hep-ph/0508022

- [3] G. Degrassi, P. Slavich, Nucl. Phys. B 825 (2010) 119-150, arXiv:hep-ph/0907.4682

- [4] A. Djouadi, J. Malinowski, M. Spira, Comput. Phys. Commun. 108 (1998) 56-74 arXiv:hep-ph/9704448

- [5] S. Narison, Phys. Lett. B 236 (1990) 474

- [6] F. Domingo, JHEP 1703 (2017) 052, arXiv:hep-ph/1612.06538

- [7] M. Drees and K. i. Hikasa, Phys. Rev. D 41 (1990) 1547

- [8] F. Domingo, U. Ellwanger and M. A. Sanchis-Lozano, Phys. Rev. Lett. 103 (2009) 111802, arXiv:hep-ph/0907.0348

- [9] F. Domingo and U. Ellwanger, JHEP 1106 (2011) 067, arXiv:hep-ph/1105.1722

- [10] F. Domingo, U. Ellwanger, JHEP 0712 (2007) 090, arXiv:hep-ph/0710.3714

- [11] F. Domingo, Eur. Phys. J. C 76 (2016) no.8, 452, arXiv:hep-ph/1512.02091

- [12] F. Domingo, U. Ellwanger, E. Fullana, C. Hugonie and M. A. Sanchis-Lozano, JHEP 0901 (2009) 061, arXiv:hep-ph/0810.4736

- [13] F. Domingo, JHEP 1104 (2011) 016, arXiv:hep-ph/1010.4701

- [14] F. Domingo and U. Ellwanger, JHEP 0807 (2008) 079, arXiv:hep-ph/0806.0733

- [15] G. Belanger, F. Boudjema, C. Hugonie, A. Pukhov, A. Semenov, JCAP 0509 (2005) 001, arXiv:hep-ph/0505142

- [16] B. Allanach et al., Comput. Phys. Commun. 180 (2009) 8-25, arXiv:0801.0045

- [17] U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 177 (2007) 399-407, arXiv:hep-ph/0612134

- [18] U. Ellwanger, C.-C. Jean-Louis, A. M. Teixeira, JHEP 0805 (2008) 044, arXiv:0803.2962

- [19] B. Allanach, M. Badziak, C. Hugonie and R. Ziegler, Phys. Rev. D 92, 015006 (2015), arXiv:1510.03143

- [20] M. Badziak, N. Desai, C. Hugonie and R. Ziegler, Eur. Phys. J. C 79 (2019) no.1, 67, arXiv:1810.05618

- [21] F. Domingo, JHEP 1506 (2015) 052, arXiv:1503.07087

- [22] D. Das, U. Ellwanger, A. M. Teixeira, Comput. Phys. Commun. 183 (2012) 774-779, arXiv:1106.5633

- [23] M. Muhlleitner, A. Djouadi, Y. Mambrini, Comput. Phys. Commun. 168 (2005) 46-70, arXiv:hep-ph/0311167

- [24] G. Bélanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 174 (2006) 577–604, arXiv:hep-ph/0405253

- [25] G. Bélanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 176 (2007) 367–382, arXiv:hep-ph/0607059

- [26] G. Bélanger, F. Boudjema, A. Pukhov, and A. Semenov, Comput. Phys. Commun. 185 (2014) 960–985, arXiv:1305.0237

- [27] G. Bélanger, F. Boudjema, A. Goudelis, A. Pukhov, B. Zaldivár, Comput. Phys. Commun. 231 (2018) 173-186, arXiv:1801.03509