NMSSMTools

# NMHDECAY, NMSPEC and NMGMSB

NMHDECAY
The Fortran code NMHDECAY computes the masses, couplings and decay widths of all Higgs bosons of the NMSSM, and the masses of all sparticles, in terms of its parameters at the electroweak (or susy breaking) scale: the Yukawa couplings lambda and kappa, the soft trilinear terms A_lambda and A_kappa, and tan(beta) and mu_eff = lambda*S. (Instead of A_lambda, the MSSM-like parameter M_A can also be used as input.) The computation of the Higgs spectrum includes leading electroweak corrections, two loop terms and propagator corrections. The computation of the decay widths is carried out as in HDECAY, but momentarily without three body decays. 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. B physics constraints from b -> s gamma, Delta M_q, B -> mu+mu- and B+ -> tau+ nu_tau are included as in ref. [4] below and updated occasionally. The dark matter relic density, direct and indirect detection cross sections can be computed via a link to a NMSSM version of the MicrOMEGAs code [3]. SLHA conventions for input and output are used. NMHDECAY is part of the NMSSMTools package that can be downloaded below.

## Authors: Ulrich Ellwanger*, John F. Gunion**, Cyril Hugonie***

* Laboratoire de Physique Theorique, Universite de Paris XI, Batiment 210, F-91405 Orsay, France.
** Davis Institute for High Energy Physics, University of California, Davis, California 95616, USA.
*** Laboratoire Physique Theorique et Astroparticules, Univ. Montpellier II, F-34095 Montpellier, France

## References:

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

NMSPEC
NMSPEC allows to compute the sparticle and Higgs masses, Higgs decay widths and couplings in the NMSSM, with soft Susy breaking terms specified at the GUT scale. Exceptions are the soft singlet mass and kappa, that are both determined in terms of the other parameters through the minimization equations of the Higgs potential. The soft Higgs and gaugino masses at the GUT scale can be chosen as non-universal, if desired. NMSPEC is included in the NMSSMTools package below.

## Authors: Ulrich Ellwanger*, Cyril Hugonie**

* Laboratoire de Physique Theorique, Universite de Paris XI, Batiment 210, F-91405 Orsay, France.
** Laboratoire Physique Theorique et Astroparticules, Univ. Montpellier II, F-34095 Montpellier, France

## References:

U. Ellwanger and C. Hugonie, "NMSPEC: A Fortran code for the sparticle and Higgs masses in the NMSSM with GUT scale boundary conditions", Comput.Phys.Commun.177:399-407,2007, arXiv:hep-ph/0612134

NMGMSB
NMGMSB allows to compute the sparticle and Higgs masses, Higgs decay widths and couplings in the NMSSM, where soft Susy breaking terms are specified at by gauge mediated supersymmetry breaking, plus terms originating from couplings of the Singlet to messengers. The boundary conditions at the messenger scale are taken from the paper below, and described in the file README. Either the soft Singlet mass or a Susy tadpole term XIS are determined in terms of the other parameters through the minimization equations of the Higgs potential. Sample input- and output files are included. It is part of the NMSSMTools package described below.

## Authors: Ulrich Ellwanger*, C.-C. Jean-Louis*

* Laboratoire de Physique Theorique, Universite de Paris XI, Batiment 210, F-91405 Orsay, France

## References:

U. Ellwanger, C.-C. Jean-Louis, A.M. Teixeira, "Phenomenology of the General NMSSM with Gauge Mediated Supersymmetry Breaking",
JHEP 0805:044 (2008), arXiv:0803.2962.
Boundary conditions following A. Delgado, G.F. Giudice and P. Slavich, Dynamical mu term in gauge mediation,'' Phys. Lett. B653 (2007) 424, arXiv:0706.3873:

## Author: Cyril Hugonie*

* Laboratoire Physique Theorique et Astroparticules, Univ. Montpellier II, F-34095 Montpellier, France

## References:

B. Allanach, M. Badziak, C. Hugonie and R. Ziegler, Gauge Mediation in the NMSSM with a Light Singlet: Sparticles within the Reach of LHC Run II,'' arXiv:1510.03143 .

NMSDECAY
NMSDECAY allows to compute sparticle widths and branching ratios. It is based on a generalization of SDECAY, 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.

## Authors: Debottam Das*, Ulrich Ellwanger*, Ana M. Teixeira**

* Laboratoire de Physique Theorique, Universite de Paris XI, Batiment 210, F-91405 Orsay, France
** Laboratoire de Physique Corpusculaire, CNRS/IN2P3 -- UMR 6533, Campus des Cezeaux, 24 Av. des Landais, F-63171 Aubiere Cedex, France

## References:

Debottam Das, Ulrich Ellwanger, Ana M. Teixeira,
"NMSDECAY: A Fortran Code for Supersymmetric Particle Decays in the Next-to-Minimal Supersymmetric Standard Model" arXiv:1106.5633 , Comput.Phys.Commun. 183 (2012) 774-779

SHOULD BE CITED TOGETHER WITH

M. Muhlleitner (Karlsruhe, Inst. Technol.), A. Djouadi (Orsay, LPT & CERN, Theory Division), Y. Mambrini (Orsay, LPT),
"SDECAY: A Fortran code for the decays of the supersymmetric particles in the MSSM", Comput.Phys.Commun.168:46-70 (2005), arXiv:hep-ph/0311167.

NMHDECAY-CPV
NMHDECAY-CPV allows to treat the CP-violating NMSSM allowing for complex input parameters at the SUSY scale, following essentially JHEP 1506 (2015) 052, arXiv:1503.07087 , by Florian Domingo.

## References:

NMSSMTools
NMSSMTools contains NMHDECAY, NMSPEC, NMGMSB, NMSDECAY and NMHDECAY-CPV.

## Authors: Ulrich Ellwanger*, Cyril Hugonie**

* Laboratoire de Physique Theorique, Universite de Paris XI, Batiment 210, F-91405 Orsay, France.
** Laboratoire Physique Theorique et Astroparticules, Univ. Montpellier II, F-34095 Montpellier, France

HOW TO USE NMSSMTOOLS: README (plain text)
Directory NMSSMTools_5.0.0 (gzipped and tared), contains NMHDECAY, NMSDECAY, NMSPEC, NMGMSB, NMHDECAY-CPV, micromegas_4.3, EXPCON for experimental constraints, and sample input- and output files.

## How to use NMSSMTools:

• After the download of NMSSMTools_x.tgz (where x is the version number), type "tar zxvf NMSSMTools_x.tgz". This will create a directory named NMSSMTools_x. Go into this directory. There you find the directories
• "sources", "main" and "micromegas", which contain source files;
• "EXPCON", which contains data files corresponding to experimental constraints;
• "SAMPLES", which contains sample input and output files;

### COMPILATION:

To compile, type first "make init", then "make". A first compilation may take a while, since most subroutines of micromegas are compiled. The following eight executable routines are created in the directory "main": nmhdecay, nmhdecay_rand, nmhdecay_grid, nmspec, nmspec_rand, nmspec_grid, nmgmsb and nmgmsb_rand. If a subroutine in the directory "sources" was modified, one has to type "make init" and "make" again. If a routine in the directory "main" was modified, it suffices to type "make" again. To delete all the already compiled codes type "make clean".

### INPUT FILES:

From the version 2.0.0 onwards, the task to be performed by an input file is independent from its name. Any name is allowed, provided it contains the three letters "inp"; it can be of the general form PREFIXinpSUFFIX where PREFIX and SUFFIX can contain dots etc..

The input file can be located in any directory specified by a PATH.

To run any input file PREFIXinpSUFFIX, type "run PATH/PREFIXinpSUFFIX". (PATH is optional; if absent, the input file has to be located in the same directory as the script file "run".)

The output files are located in the directory specified by PATH. They have the following format:
If one single point in the parameter space is evaluated:
PREFIXspectrSUFFIX (includes the previous PREFIXdecaySUFFIX), PREFIXomegaSUFFIX (if the relic density is computed)
If scans are performed:
PREFIXerrSUFFIX as well as PREFIXoutSUFFIX
However, the task to be performed by an input file must be specified in the BLOCK MODSEL at the beginning (see the SLHA2 conventions in B. Allanach et al., SUSY Les Houches Accord 2, Comp.Phys.Commun.180:8-25,2009, arXiv:0801.0045.

For more details see the description in README .

# Latest update: 19.10.2016

## History:

• 19.10.2016: version 5.0.0
- Treatment of the CP-violating NMSSM made possible allowing for complex input parameters at the SUSY scale, following essentially JHEP 1506 (2015) 052, arXiv:1503.07087 , by Florian Domingo.
- Updated constraints on the SM-like Higgs couplings from the ATLAS and CMS combination in JHEP 1608 (2016) 045, arXiv:1606.02266 , with thanks to Jeremy Bernon
- New constraints from SM-like Higgs decaying to 2 CP-odd Higgs states with 4 fermions in the final state from ATLAS ( arXiv:1505.01609 , Phys.Rev. D92 (2015) no.5, 052002) and CMS ( arXiv:1510.06534 , JHEP 1601 (2016) 079), arXiv:1506.00424 , Phys.Lett. B752 (2016) 146-168), CMS-PAS-HIG-14-019, CMS-PAS-HIG-14-022, CMS-PAS-HIG-14-041, CMS-PAS-HIG-15-011).
- New and updated constraints on dark matter direct detection (both spin dependent and spin independent) from LUX ( arXiv:1602.03489 , Phys.Rev.Lett. 116 (2016) no.16, 161302, and arXiv:1608.07648 ), PICO ( arXiv:1510.07754 , Phys.Rev. D93 (2016) no.5, 052014, and arXiv:1601.03729 , Phys.Rev. D93 (2016) no.6, 061101) and PandaX-II ( arXiv:1607.07400 ).
- Interface with the version 4.3 of micromegas (including improved computation of the freeze-out temperature), with thanks to Alexander Pukhov.
- Bug corrected in subroutine decay.f (affecting the Higgs to bb partial width).
- Corrected Z3-breaking NMSSM-specific contributions to the Higgs masses from charginos/neutralinos (following a report from Pawel Szczerbiak).
- Corrections of the Z3-breaking terms in the constraints from false minima (subroutine checkmin.f).
- Corrections of bugs in the subroutine DMSQUARK_YUK in msferm.f (affecting sfermion masses), with thanks to Florian Domingo.
- Corrections of bugs in the subroutine bsg.f (B physics constraints) thanks to Florian Domingo and Jonathan Da Silva.
- Numerous corrections of small bugs like COMMON BLOCK size not matching, array indices out of bounds, implicit types ill-defined, possible divisions by zero, overflows/underflows. Possible compiler-dependent effects.
• 30.06.2016: version 4.9.3
- Minor update of bsg.f, with thanks to Florian Domingo.
• 17.05.2016: version 4.9.2
- Constraints from CMS-PAS-HIG-13-010 on H_125->AA->4mu corrected (the bounds were too low by a factor 1/100), with thanks to Robin Cameron Aggleton.
• 18.03.2016: version 4.9.1
- Gravitino relic density for nmspec_rand and nmspec_mcmc and the output of XIFGUT, XISGUT in nmspec_grid.f, nmspec_rand.f and nmspec_mcmc.f corrected, with thanks to Yong-Gon Yi.
• 04.03.2016: version 4.9.0
- New tests of K-physics: K+-> Pi+ nu nubar, KL-> Pi0 nu nubar, DMK, epsK, thanks to Florian Domingo.
- B-physics: Updates of previous and addition of new processes: B->Xs mu+mu-, B->Xs mu+mu-, b -> d gamma, Bd->mu+mu-, B-> Xs nu nubar, B+-> K+ nu nubar, B-> Ks nu nubar, BR[B+ -> D tau+ nu_tau]/BR[B+ -> D l+ nu_l], BR[B+ -> D* tau+ nu_tau]/BR[B+ -> D* l+ nu_l], thanks to Florian Domingo.
- Tests of SM Higgs couplings: Updates of parameters (based on Lilith-1.1.3 with database 15.09, thanks to J. Bernon).
- Decays of light scalars: Muon/strange quark loops added to H/A -> gammagamma (but BRs not reliable for M_H/A < 1 GeV).
- msferm.f: (minor) bugs corrected, thanks to Florian Domingo. <
- The gravitino mass and, for flag I7=2, branching ratios into the gravitino are computed.
- Gauge mediation: for flag I3=/=0, the non-thermal relic density due to the gravitino is computed by MicrOmegas from the gravitino/NLSP mass ratio.
• 01.10.2015: version 4.8.2
- In bsg.f: Delta M_d,s also upgraded to NLO, numerically unstable functions corrected (thanks to Daniele Barducci for spotting), with thanks to Florian Domingo, see arXiv:1512.02091.
• 23.09.2015: version 4.8.1
- Added missing files to the update of bsg.f in 4.8.0.
- Missing (rare) Higgs-to-Higgs decays are printed in the output, with thanks to Marcin Badziak.
• 09.09.2015: version 4.8.0
- Major update of bsg.f, with thanks to Florian Domingo, see arXiv:1512.02091: Most B-physics observables are upgraded to full NLO; recent SM results of 1503.01789 (B -> Xs gamma), 1311.0903 (Bs -> mu+mu-) and 1503.04849 (B -> Xsl+l-) have been taken into account. Experimental and lattice input from current HFAG and FLAG data.
• 13.08.2015: version 4.7.1
- Bug in lhchig.f concerning chi2gam/chi2zz in case of two scalars near 125 GeV corrected, with thanks to Nils-Erik Bomark
• 02.07.2015: version 4.7.0
- For Higgs mass precision >=1, running stop/sbottom masses and trilin. couplings at the scale QSTSB are computed without threshold corrections such that they are defined in DR_bar (in order to comply with SLHA).
- The SUSY scale can now be defined in BLOCK EXTPAR, item 0 (in order to comply with SLHA).
- Yukawa corrections ~ht, hb to the stop/sbottom pole masses added
- Lmunu in runpar takes into account the Z_3 violating contribution MUP to the singlino mass (with thanks to F. Domingo).
- Chi2 of the signal rates of the SM-like Higgs boson are computed for the CP-even Higgs with mass closest to 125.1 GeV, indep. of its mass
- Gauge and SM Yukawa couplings at the SUSY scale are computed in an extra routine "getsusycoup" attached to runpar.f
• 22.04.2015: version 4.6.0
- NMSDECAY can be called for GMSB models, decays into a gravitino are included
- More general GMSB models of the type discussed by A. Delgado, G. F. Giudice and P. Slavich in PLB 653, 424 (2007) arXiv:0706.3873 can be defined.
- Correction of bugs concerning selectrons/smuons in NMSDECAY (thanks to Daniele Barducci and Nishita Desai for spotting)
- Bugs in decays for very light (< 1 GeV) Higgs bosons corrected
- NMSPEC improved for very large (> 10 TeV) masses at the GUT scale
- Optional output without privately defined SLHA-BLOCK's (and some additional BLOCK's); useful as param_card.dat for MadGraph as indicated in README
• 25.02.2015: version 4.5.1
- Correction of a bug in the compilation of micrOMEGAs on Mac (with thanks to A. Pukhov)
- Improvement of the convergence algorithm in NMSPEC.
- Contributions of two distinct Higgs bosons to the gammagamma/ZZ signal rates are summed only if the mass difference is below 3 GeV.
- Constraint on the invisible Z-width reduced to 0.5 MeV due to update of the SM contribution in A. Freitas, 1401.2447.
• 18.12.2014: version 4.5.0
- Calculation of relic density in MicrOmegas made possible for LSP masses below 1 GeV (with thanks to A. Pukhov)
- Scalar and pseudoscalar Higgs masses below 1 GeV made possible (branching ratios not quite reliable, however)
- The switch 1 in BLOCK MODSEL allows to insist on the Z_3-invariant NMSSM for the values 3,4,5
• 21.11.2014: version 4.4.1
- Hadronic parameters for B-physics are updated, with thanks to F. Domingo. For details see the file bsg.f.
- Calculation of the relic density in MicrOmegas corrected for non-Z_3-invariant terms, with thanks to J. Gainer et al. for raising the issue, and A. Semenov and A. Pukhov for the corrections.
- Calculation of the relic density in MicrOmegas improved in cases of s-channel resonances, with thanks to A. Pukhov.
- Minor bug (sign) in mgluino.f corrected, impact on the gluino mass <~2 GeV.
• 15.09.2014: version 4.4.0
- Constraints from B-physics are updated using values from the HFAG web page. For details see the file bsg.f.
- The default value of the SM-like Higgs mass is 125.1 GeV. Constraints on signal strengths are updated, with thanks to J. Bernon, B. Dumont and S. Kraml.
- Constraints on H->tautau from ATLAS-CONF-2014-049 added to the ones from CMS.
- Constraints from CMS-PAS-HIG-13-010 on H_125->AA->4mu added.
- Regularisation of IR logarithm in mhiggs.f for very a light singlino and Higgs precision 0, with thanks to M. Cahill-Rowley, J. Gainer, J. Hewett and T. Rizzo.
- Regularisation of some IR logarithms for a very light pseudoscalar, and of some potential divergencies for very heavy squarks for Higgs precision 1-2.
- The use of Higgs precision 1-2 is blocked in the presence of Z_3 violating terms, since this part of the code does not include Z_3 violating terms.
• 26.06.2014: version 4.3.0
- Subroutines for Markov Chain Monte Carlo scans in parameter space are provided, separately for the general NMSSM (NMHDECAY_MCMC), the NMSSM with parameters defined at the GUT scale (NMSPEC_MCMC) and the NMSSM with gauge mediated supersymmetry breaking (NMGMSB_MCMC). For applications see README and the files provided in the directory SAMPLES.
- corrections in the Higgs mass calculation (Higgs mass precision 1 and 2):
- correction of the two loop term in the MSbar-DRbar transition of h_top (with thanks to M. Muehlleitner and P. Slavich)
- the light Higgs contribution to the transition from the top pole mass to the running top mass is expressed in terms of h_top at m_top and not M_SUSY; increases the SM-like Higgs mass by 0.2-0.3 GeV (with thanks to P. Slavich)
• 08.01.2014: Commented checks for PROBs in nmhdecay.rand uncommented (no change of version number)

• 10.12.2013: version 4.2.1
- BLOCKs in the files *spectr* for input for HiggsBounds corrected
- (minor) bug in ftpar.f corrected, with thanks to S. Trojanowski
• 6.12.2013: version 4.2.0
- LUX constraints on the direct detection cross section implemented
- Constraints from CMS-PAS-HIG-12-050 on H->tautau implemented (quite conservative for MA<~250 GeV, but realistic for MA>~300 GeV)
- Higgs decays into two off-shell Z/W-bosons can be included via a flag (from Hdecay, http://people.web.psi.ch/spira/hdecay/)
- Bug in the reduced couplings of CP-odd Higgses to gluons/photons corrected (the branching fractions were not affected)
- The reduced couplings of Higgses to b-quarks include now the DELMB corrections (also used in loop induced couplings)
- The default value of the b-quark pole mass in init.f is set to 4.54 GeV
- Higgs signal rates into invisible final states (pairs of neutralinos) indicated in the output
- In the output *spectr*: BLOCKs with parameters at the GUT scale renamed
- The content of the previous output files *decay* is appended to the output files *spectr*
- The output files *spectr* can be used as input for HiggsBounds from http://higgsbounds.hepforge.org/ (BLOCKs added)
- For Mac users: new compiler (clang) used for compatibility with OSX 10.9.
• 18.10.2013: version 4.1.2
- bug in sources/fullhiggs.f corrected. Leads to an increase of the SM-like Higgs mass for PFLAG=1 or 2 and MA >> mtop.
- minor corrections in MicrOmegas and the files sources/full* in order to avoid NaNs for unphysical parameter sets.
• 27.08.2013: version 4.1.1
- bug in the total width of stau2 in NS_sleptons.f corrected, with thanks to L. Duflot.
- bug in ftpar.f ~MUP corrected, with thanks to Doyoun Kim.
• 21.07.2013: version 4.1.0
- constraints on a SM-like Higgs mass (125.7 \pm 3 GeV) and its signal rates (from arXiv 1306.2941 ) added; modifiable in the subroutine Higgs_CHI2 at the end of sources/lhchig.f.
- guess for the initial values of the parameters in nmspec_rand.f corrected; with thanks to Jiang Yun.
• 23.05.2013: version 4.0.0
- The general NMSSM with non-Z_3 invariant parameters either at the (low) SUSY scale or at the (high) GUT scale is made possible, as well as the nMSSM with kappa=0. (See README for the allowed options; now most parameters not given in input files are assumed to be = 0).
- NMSSM specific radiative corrections to CP-odd and charged Higgs masses from chargino/neutralino loops have been added for "Higgs precision = 0".
- Chargino/neutralino 1-loop corrections to the CP-even Higgs masses use parameters at the scale QSTSB (resumming some large logarithms), and some NMSSM specific large logarithms have been added for "Higgs precision = 0".
- The relic density and the direct detection cross sections are computed with the new version micromegas 3.0 (as of May 2013), which includes virtual contributions to neutralino annihilation cross sections (e.g. neutralino-neutralino -> WW) and loop contributions to neutralino+neutralino -> gamma+gamma. (Hence a first compilation of the code and a calculation of the relic density takes longer than before. In relden.F, virtual decays can be switched off setting setVVdecay(0) instead of the default setVVdecay(1).)
- The direct detection cross section is computed in terms of sigma_PiN and sigma_S instead of sigma_PiN and sigma_0; now sigma_PiN and sigma_S can be specified in init.f.
- The exp. constraint on dark matter relic density is now 0.1187 +/- 0.0017 (arXiv:1303.5076), but one has to allow for theoretical errors which are difficult to estimate. The default window in init.f is 0.1187 +/- 0.0085, but can be modified by the user.
- New interface between NMSSMTools and\A0 micromegas 3.0 using the full Higgs effective lagrangian with parameters computed at the same precision level as for Higgs masses and couplings. These parameters can be found in the\A0BLOCK EFFECTIVE_COUPLINGS at the end of the output file #spectr#.
- With "Higgs precision = 1 or 2", the relic density and the direct detection cross sections are computed by micromegas with the correct Higgs masses (previously micromegas used always the Higgs masses as computed with "Higgs precision = 0").
- The diagonal chargino loop contributions to H/A->Z+gamma are included (with thanks to Peiwen Wu for discussions).
• 30.01.2013: version 3.2.4
- Factor 2 in the charged Higgs loop contribution to H->gammagamma and H->Zgamma corrected (with thanks to Peiwen Wu, L. Sluter and P. Richardson; impact on the partial widths <2% for a charged Higgs mass >~ 300 GeV)
• 11.01.2013: version 3.2.3
- XENON100 bounds corrected, with thanks to pwwu@itp.ac.cn
- lhchig.f updated: more reduced significances for Higgs bosons at the LHC added, but LHC bounds omitted (left to the user); notation of the reduced significances in the ouput files corrected (with thanks to C. Beskidt)
• 21.12.2012: version 3.2.2
- Constraints from B-physics updated
- Constraints from XENON100 updated
- Bugs in tdecay concerning t -> H+ + b (decays near the threshold, and exp. constraints) corrected (with thanks to Neil Christensen and Zhen Liu)
- Convergence of the charged Higgs mass calculation in fullhig insured (with thanks to Jinmian Li)
• 12.7.2012: version 3.2.1
- Provisional suppression of LHC constraints on Higgs signal cross sections (to be updated)
• 23.3.2012: version 3.2.0
- In NMSPEC: improved convergence of RGEs with boundary conditions at the Susy and GUT scales for large Yukawa couplings
- In NMSPEC: mu_eff and kappa at the Susy scale can be given as input, then MHD and MHU at the GUT scale become outputs
- bounds from XENON100 on the direct detection cross section can be applied
- LHCb constraints on BR(Bs->mu+mu-) applied
- LHC constraints on Higgs signal cross sections updated
- code of M3H^2 in the output corrected (with thanks to B. Fuks)
• 16.1.2012: version 3.1.0
- Constraints from single (not combined) Higgs search channels H->tautau, H->bb, H->ZZ, H->WW, top->bH+ ->tau+nu (CMS) and H->2gamma (ATLAS) added
- Compilation improved to be compatible with Fedora
- Improvement of the numerical integration of the Boltzmann equations inside MicrOmegas
- CP-odd Higgs mixing corrected (for PFLAG > 0 only)
- The fine-tuning with respect to the input parameters is computed following arXiv:1107.2472 (The command "CALL FTPAR(PAR,...)" in MAIN is commented for scans)
• 7.10.2011: version 3.0.2
- NMSDECAY: a warning appears on the screen, of QCD corrections to colored sparticle decay widths are negative and larger than the tree level widths. Then the relative corrections "corr", leading to a negative width when added according to (1+corr), are replaced by 1/(1-corr).
- the measured (not computed) W- and Z-widths are added to the output of NMSDECAY
- the BLOCKs for parameters at the GUT-scale in the SLHA-output of NMHDECAY and NMSPEC are renamed
- Higgs-significances at the LHC are computed only if a corresponding flag "7" in the input file is switched on
• 12.9.2011: version 3.0.1
- NMSDECAY corrected, with special thanks to Oscar Stal: negative chargino masses now considered correctly; output corrected
- RGE for the tadpole term xi_S corrected, with special thanks to Renato Fonseca
- constraints from B_s -> mu mu updated (comb. LHCb + CMS)
- constraints from H/A -> tau tau updated according to CMS-PAS-HIG-11-020
• 29.6.2011: Release of version 3.0.0
- NMSDECAY is added, which allows to compute sparticle widths and branching ratios (see above). It is based on a generalization of SDECAY, including the corresponding QCD corrections and 3-body decay modes. With special thanks to M. Muhlleitner, A. Djouadi and Y. Mambrini.
- gfortran is set as default compiler for LINUX. f77 is no longer recommended due to potential problems with complex functions in NMSDECAY.
• 17.5.2011: version 2.3.7
- charged Higgs mass with complete 1-loop and Yukawa induced 2-loop corrections (if improved corrections to Higgs masses are switched on), with special thanks to K.H. Phan and P. Slavich
• 11.5.2011: bug in H->bb corrected
• 22.4.2011: version 2.3.6
- LHC constraints on H/A -> 2 tau included
- constraints from ALEPH on H->AA->4tau removed for MA > 9.4 GeV (due to A-eta_b mixing, see arXiv:1105.1722)
- decays H/A -> bb allowed only for MH/MA > 2*MB ~ 10.6 GeV
• 4.2.2011: version 2.3.5
- top quark decays into charged Higgs and top squarks + chi_0 included
- Tevatron constraints on t -> b + H+, H+ -> c+s or tau+nu_tau or A1+A1 included
• 4.1.2011: version 2.3.4
- updated constraints on B_s -> mu+ mu- and B+ -> tau+ nu_tau
• 13.12.2010: version 2.3.3
- corrected constraints from ALEPH (JHEP 1005:049,2010, arXiv:1003.0705)
• 14.10.2010: version 2.3.2
- the improved corrections to the Higgs masses as in G. Degrassi, P. Slavich, Nucl.Phys.B825:119-150,2010, arXiv 0907.4682 , can be switched on (without or with pole mass corrections), with special thanks to P. Slavich
- improved calculation of h_top(Q_Susy)/Z-factors for rad. corrections
- bug in the 2-loop RGEs for A_lambda corrected (with thanks to F. Staub)
- bug in the BRs(A -> sfermions) corrected (with thanks to L. Suter)
- normalisation of the partial width H -> W+H- corrected (with thanks to L. Suter)
- slight correction of the function AQCD in decay.f
- constraints from ALEPH (JHEP 1005:049,2010, arXiv:1003.0705) on H->AA->4tau included
- corrections in the output SLHA: Higgs vev v(Q) properly normalised, smuons indexed by L,R
• 18.11.2009: version 2.3.1
- benchmark points corrected
- small correction in the radiative contributions in decay.f (for M_Higgs > M_stop)
- output format of the GUT-/messenger scale corrected
• 10.11.2009: Release of version 2.3.0
- constraints from B_s -> X_s mu+ mu- added (with thanks to F. Domingo)
- algorithm for the iterative determination of mu_eff in nmspec improved
- conventions for the Susy and soft singlet masses mu' and m_S^2' in nmgmsb modified according to the new version of SLHA2 (additional factors 1/2)
- small corrections in the radiative contributions in decay.f, lowmuk.f, checkmin.f and the two loop RGE for h_t (-> slightly modified benchmark points)
• 31.3.2009: Release of version 2.2.0
- Two loop RGEs for all couplings and soft terms including the NMSSM-specific contributions completed, pole mass corrections for the CP-odd and charged Higgs scalars improved.
- Corresponding small changes of the numerical values of the benchmark points.
• 6.2.2009: version 2.1.1
- Version "grid" added to NMGMSB
- Benchmark point P1 in BMPGMSB corrected
• 10.12.2008: Release of version 2.1.0
- Significances for the detection of all neutral Higgs bosons at the LHC (at low luminosity, 30 fb-1, and at high luminosity, 300 fb-1) can be computed, following U. Ellwanger, J.F. Gunion and C. Hugonie, JHEP {0507} (2005) 041, arXiv:hep-ph/0503203, and G. Belanger, C. Hugonie and A. Pukhov, arXiv:0811.3224 [hep-ph].
- New version of micrOMEGAS (2.2) included, which allows to compute the LSP annihilation cross section, the resulting photon spectrum for indirect detection and spin dependent/spin independent LSP-nucleon cross sections for direct detection, following G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, arXiv:0803.2360 [hep-ph], and G. Belanger, C. Hugonie and A. Pukhov, arXiv:0811.3224 [hep-ph].
- Constraints from radiative Upsilon decays and the eta_b mass measurement are implemented as in F. Domingo, U. Ellwanger, E. Fullana, C. Hugonie and M.A. Sanchis-Lozano, arXiv:0810.4736 [hep-ph].
- Initial guess for the parameters at the scale Q2 in the routines NMGMSB corrected (with thanks to Hisano, Junji)
• 30.5.2008: version 2.0.1
- Bug in the rad. corr. to M2 in mneu.f corrected (with thanks to Junjie Cao)
- Output of the BLOCK HMIX corrected in order to coincide with SLHA2 conventions (with thanks to A. Pukhov)
- Inappropriate floats in bsg.f corrected (with thanks to A. Pukhov)
- Warnings in case of inappropriate input files modified
• 13.05.2008: Release of version 2.0.0
- The architecture of the code has been modified: From the version 2.0.0 onwards, the task to be performed by an input file is independent from its name. However, the task to be performed by an input file must be specified in the BLOCK MODSEL at the beginning, see the new description README.
- The routines NMGMSB have been added, which allow for boundary conditions specified by gauge mediated supersymmetry breaking (see above).
- SM error bars for the muon anomalous magnetic moment have been updated. Still: the SM hadronic contribution to (g-2)_muon is estimated from e+e- data only, data from tau decays is ignored; consequently the exp. BNL result differs from the SM by ~3 sigma. If the line "CALL MAGNMU(PAR,PROB)" is present in the main routines, it is checked whether a point in parameter space differs by at most by 2 sigma from the exp. result. Hence, points in parameter space for which the Susy contribuion is too small (or vanishing or of the wrong sign) generate warnings or are considered as phenomenologically excluded. If this constraint is considered as too strong by the user it is enough to comment the line "CALL MAGNMU(PAR,PROB)" in the main program.
- Exp. constraints from the Tevatron on sbottoms updated.
• 17.03.2008: Release of version 1.3.0
- The muon anomalous magnetic moment is computed (to two loops) according to F. Domingo and U. Ellwanger, JHEP 0712:090 (2007), arXiv:0710.3714 [hep-ph]; a warning is issued if the susy contribution incl. theoret. errors is less than ~10^{-9} or larger than ~5*10^{-9}
- The input of the soft parameters of the first two generations is simplified for the versions nmhdecay
- The versions nmspec_grid and _rand allow to scan over non-universal soft parameters at the GUT scale
- A script "run" is added which allows to run user-defined input files
• 25.02.2008: version 1.2.5, 2 loop RGE for A_muon corrected
• 20.02.2008: version 1.2.4
- Bug in electroweak corrs. in bsg.f corrected (effect essentially on B->mumu) -> the Benchmark point 3 in BMP.tgz above has been updated
- The Smuon masses (and the trilinear coupling A_muon) are treated separately from the Selectron masses
- Small bug in the error management (if PROB(I) < 0) corrected
• 08.02.2008: version 1.2.3
- Bug in msferm corrected (slepton masses of the two first generations)
- Small bug in mhiggs corrected (bosonic loop corrs.; effect <~1 MeV on Mh_1)
- Floats rewritten as "D0"; Fortran compiler g77 used if available
• 05.02.2008: version 1.2.2
- Treatment of floats in bsg.f improved, consistency of the output format with SLHA2 improved
- Benchmark points from the paper "Benchmark scenarios for the NMSSM", A. Djouadi et al., arXiv:0801.4321 [hep-ph], made available (see Downloads above)
• 03.12.2007: version 1.2.1, small bug in the RGE for lambda corrected
• 30.10.2007: Release of version 1.2.0:
- inclusion of constraints from B physics according to F. Domingo and U. Ellwanger, JHEP 0712:090 (2007), arXiv:0710.3714 [hep-ph]
- non-universal boundary conditions at the GUT scale in NMSPEC are made possible
- improved convergence of NMSPEC for large lambda, improved numerical agreement between NMHDECAY and NMSPEC
- the scale of mu_eff in the subroutine msferm.f corrected (with thanks to Junjie Cao)
• 12.07.2007: version 1.1.1:
- improved convergence of nmspec near focus point regions
- bug in MA-dependent threshold corrections for g1 and g2 in rges.f corrected (affects only the calculation of MGUT)
• 11.06.2007: Bug in the large tan(beta) corrections to the branching ratios (in decay.f) corrected (with thanks to I. Rottlaender)
• 30.03.2007: Bug in the error management in nmhdecay (in case of negative Higgs masses squared) corrected
• 27.03.2007: Release of version 1.1.0:
- the MSSM-like parameter MA at M_SUSY can be used as input instead of A_lambda (within nmhdecay, nmhdecay_scan or nmhdecay_grid)
- improved diagonalisation of Higgs mass matrices (avoids problems in cases with huge mass differences)
- bug in the calculation of the Higgs wave function renormalisation constants Z in the case MA >> M_susy corrected
- if soft Higgs masses - as computed internally in nmhdecay - are much larger than the SUSY scale, a warning is indicated
- the micrOMEGAs version included in the NMSSMTools package has been updated:
- bugs in the calculation of the Higgs widths in micrOMEGAs corrected
- if the Higgs widths computed by micrOMEGAs are still unreliable (which could affect the relic density if the LSP annihilates through a Higgs resonance) a warning is indicated in the output file
• 31.01.2007: version 1.0.2: bug in the Higgs loop corrections to M_higgs corrected,
• 08.01.2007: version 1.0.1: bug in the Susy threshold corrections in the running of h_top and h_bottom, and the electroweak two loop corrections corrected (with thanks to P. Slavich)
• 12.12.2006: Release of NMSSMTools, version 1.0.0