NMSSMCALC version 3.00, 27 March 2019: These changes concern both the complex and the real version of NMSSMCALC. They have been made by Thi Nhung Dao.

  • Added new files CalcMasses_atat.F, TLMT4.h to compute the O(alpha_t^2) corrections to the Higgs boson masses of the complex NMSSM and changed makefile accordingly.
  • CalcMasses.F has been changed to add the O(alpha_t^2) corrections.
  • The file conversion_polemass_DRbarmass.f has been added to perform the conversion from the top-quark pole mass to the DRbar mass taking into account the O(alpha_t^2) contributons and to compute the strong coupling g3 at the SUSY scale. The makefile has been changed accordingly.
  • Changes in CalcMasses.F to include also the EW corrections to the conversion of the DRbar parameters At, M_Qtilde_3 and M_tR to the OS parameters.
  • Change of g3 in CalcMasses2loop.F to be in accordance with our O(alpha_t^2) corrections.
  • This change only concerns the complex version: Changed the electron EDM bound to the new experimental value of November 2018.


  • NMSSMCALC version 2.00, 14 January 2019

  • Correction of definition of absAu(3) and phiAu(3) in CalcMasses.F (real and complex case).
  • Correction of definition of ReAlambdaDR and ImAlambdaDR in CalcMasses2loop.F (real and complex case).
  • Correction of decay width H+ -> W H3 in bhdecay_c.f (so only complex case): amh2 changed to amh3.


  • NMSSMCALC version 2.00, 28 December 2018

  • Added new file Zfactors_New.F and changed makefile accordingly. With this the two-point functions and hence the Z-factor and one-loop Higgs self-energies can be computed with complex momentum. This gives stable results near thresholds. (The changes are denoted by 'Nhung changed/added ...'.) In accordance with this change, there are changes in CalcMasses.F. (They are marked by '(end) Maggie changed 21 Aug 2018'.)
  • Changes in CalcMasses.F in order to account for the fact that the SLHA2 requires Re(Alambda) and Im(Alambda) as input. (The changes are marked by '(end) MMM changed 22 Aug 2018'.)
  • Changes in CalcMasses2loop.F in the computation of Sqrt(EPhiu). It is replace by Cos(phiu/2d0)+(0d0,1d0)*Sin(phiu/2d0) (The changes are marked by '(end) MMM changed 21 Aug 2018'.)
  • These changes concern both the real and the complex NMSSMCALC programs.


  • NMSSMCALC version 2.00, 8 October 2017

  • CalcMasses2loop.F: Added 'dmtep = 0' to set the epsilon part of the counterterm to zero in accordance with our OS renormalization scheme.


  • NMSSMCALC version 2.00, 22 April 2016

  • CalcMasses2loop.F: Bug fixed in the 2-loop counterterm for the charged Higgs mass, when A_lambda is given as input. Thanks to Florian Staub


  • NMSSMCALC version 2.00, 13 August 2015

  • CalcMasses.F and nmssmcalc.f have been changed as now in the complex NMSSM also the electric dipole moments (EDMs) are calculated. The computation is performed with the files EDMs.F and EDMloopfunc.F, which are included in the tar file for the complex NMSSM.
  • bhdecay.f, bhdecay_c.f and susylha_nmssm.f changed to write out now *all* parameters in the block EXTPAR. This is also the block that is used in inp.dat, respectively, inp_c.dat. Furthermore, now all switches, that have been set by the user, are given out.


  • NMSSMCALC version 1.03, 7 June 2015

  • bhdecay.f and bhdecay_c.f changed to write out now also the Higgs coupling blocks read in by HiggsBounds.
  • bhdecay_c.f and bhdecay_f: If there are problems in the Higgs-to-Higgs off-shell decays, their decay width is set to zero.
  • bhdecay_c.f and bhdecay.f: The total off-shell decay width has been initialized to 0: hitohhoff(i)=0, to avoid problems when NMSSMCALC is run in a loop over parameters.
  • bhdecay_c.f: Bug fix in Hi-> bb, i.e. ratcoup = ghit(i)/cdabs(xghigb(i)) corrected to ratcoup = ghit(i)/cdabs(xghib(i)). Thanks to Shoaib Munir.


  • NMSSMCALC version 1.03, 11 February 2015

  • Switched from the MSbar running value of alphas to the DRbar running value.
  • Bug that caused wrong output values in 'Block GAUGE Q=...' fixed.
  • Bug that caused program to crash for At=0 fixed.
  • Bug that caused program to crash for M3=M_Q3 fixed.


  • NMSSMCALC version 1.03, 12 January 2015

  • Bug fixed in CalcMasses.F and CalcMasses2loop.F which occurs for degenerate sbottom masses. Thanks to Florian Staub and Shoaib Munir for pointing this out to us.


  • NMSSMCALC version 1.03, 5 January 2015

  • CalcMasses now version 2.00: Implementation of the two-loop O(alpha_t alpha_s) corrections to the NMSSM Higgs boson masses in the real and the complex NMSSM.
  • Due to this implementation there is a new additional file CalcMasses2loop.F.
  • In inp.dat the Block MODSEL has been extended, so that the user can choose the loop level for the calculation of the loop corrected masses and the renormalisation scheme in the top/stop sector, i.e. DR_bar or OS.
  • The implementation of the two-loop corrections furthermore lead to small changes in the files susylha_nmssm.f and in the Makefile.
  • Bug fix in hsqsq_nmssmcalc.f: xlamvs/sqrt(2) changed to xlamvs.
  • Changes in bhdecay.f and bhdecay_c.f so that warnings are given out in case there are problems in the off-shell decays.
  • nmssmcalc.f: Changes in print statements on the screen.


  • NMSSMCALC version 1.02, 24 August 2014

  • nmssmcalc.f: Changes in print statements on the screen.


  • Changes for NMSSMCALC version 1.02, 30 June 2014

  • CalcMasses.F: in subroutine neutralinomasses: fixed a bug in an if-statement.
  • CalcMasses.F: changed interpretation of a negative mu-value for the real limit.
  • bhdecay.f: fixed a definition in the common block of elemhtt_NMSSM.


  • NMSSMCALC version 1.02, 11 April 2014

  • Off-shell Higgs-to-Higgs decays stabilized in real and complex NMSSM. (Thanks to S. Munir)
  • Warning is issued if QCD corrections to Higgs decays into stops/sbottoms become too large in some extreme scenarios. In this case the tree-level decay widths are calculated. This is only relevant for the real NMSSM. (Thanks to S. Munir)
  • tan(beta) is now consistently read in from EXTPAR. (Thanks to the CPC referee.)
  • Consistent use of VEV.


    NMSSMCALC version 1.01, based now on HDECAY 6.10, 22 January 2014

  • Typo corrected in a do loop in CalcMasses.F (thanks to S. Munir)
  • If the mass calculation does not converge sufficiently fast, a warning is issued and the last masses are given out by CalcMasses.F.
  • Changes in bhdecay.f: signs in Delta_tau/mu corrected, elw. corrections extended, elw. corrections to Delta_b/s added; pdg code for H1 in output file corrected (thanks to S. Munir); warning included in case of problems in off-shell decays.
  • Changes in bhdecay_c.f: signs in Delta_tau/mu corrected; pdg code for H1 in output file corrected (thanks to S. Munir); warning included in case of problems in off-shell decays.