NMSSMCALC version 5.10, xx October 2022, changes made by Thi Nhung Dao:

complex and real version: implementation of the loop-corrected Higgs self-couplings and of the possibility to compute the loop-corrected Higgs-to-Higgs decays up to O(alpha_s*alpha_t + alpha_t^2).

NMSSMCALC version 5.00, 27 July 2022, changes made by Thi Nhung Dao:

complex and real version: implementation of the electron and muon anomalous magnetic moments (AMM)

NMSSMCALC version 4.00, 15 June 2021: These changes concern both the complex and the real version of NMSSMCALC. They have been made by Thi Nhung Dao.

Added new files for the computation of the O((alpha_t+alpha_lambda+alpha_kappa)^2 + alpha_t alpha_s) corrections: CalcMasses_atlambda.F (computes renormalised Higgs self-energies and two-loop counterterms), CalcMasses_atlambdaH.F (computes unrenormalised Higgs self-energies), CalcMasses_atlambda_W.F (computes renormalised W-boson self-energies), CalcMasses_atlambda_Z.F (computes renormalised Z-bson self-energies); 2LOOPLIB.F (contains one- and two-loop integrals); TLMT4lambda.h (contains the declaration of new variables); changed makefile accordingly.

In the file inp.dat in the block MODSEL the new two-loop corrections have been added.

In the file inp.dat a new block called REGFACTOR has been added to choose the value of the regulator mass MG for the regularisation of the IR divergences appearing in the new two-loop corrections. This is done by setting the value of R, with R=10^3*MG^2/MUR^2, where MUR is the renormalisation scale.

The output of the effective couplings has been adapted to the conventions of HiggsBounds 5.

The output now also constains the Z-factors of the Higgs bosons.

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.

Note that also the labels for the input 7 in the Block MODSEL (concering the renormalization of the top/stop sector) in the input files inp.dat and inp_c.dat changed: 1 now means DRbar scheme w/o gauge running, 2: DRbar running w/ gauge running, 3: OS scheme. For details, we refer to arXiv:1903.11358.

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.