NMSSMTools

TOOLS FOR THE CALCULATION OF THE


HIGGS AND SPARTICLE SPECTRUM IN THE NMSSM:


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. The dark matter relic density 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:



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:




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 (amongst others for the points P1 - P5 in the paper) 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 [hep-ph].



NMSDECAY
NMSDECAY allows to compute sparticle widths and branching ratios (not yet for NMGMSB). 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 [hep-ph].

    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.



NMSSMTools
NMSSMTools contains NMHDECAY, NMSPEC, NMGMSB and NMSDECAY.

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

Downloads:
    HOW TO USE NMSSMTOOLS: README (plain text)
    Directory NMSSMTools_3.1.0 (gzipped and tared), contains NMHDECAY, NMSDECAY, NMSPEC, NMGMSB, micromegas_2.2, EXPCON for experimental constraints, and sample input- and output files (amongst others the benchmark points from in the paper "Benchmark scenarios for the NMSSM", A. Djouadi et al., arXiv:0801.4321 [hep-ph].)

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_2.2", which contain source files;
  • "EXPCON", which contains data files corresponding to experimental constraints;
  • "SAMPLES", which contains sample input and output files;
  • "BMPSUGRA" contains mSUGRA-like benchmark points as in A. Djouadi et al., arXiv:0801.4321 [hep-ph];
  • "BMPGMSB" contains GMSB-like benchmark points as in U. Ellwanger et al., arXiv:0803.2962 [hep-ph].

    COMPILATION:

    To compile, type first "make init", then "make". A first compilation may take a while, since most subroutines of micromegas_2.2 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, PREFIXdecaySUFFIX, PREFIXlhcsigSUFFIX (if Higgs significances are computed), 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, arXiv:0801.0045 [hep-ph]).

    For more details see the description in README .


    Latest update: 16.1.2012

    History:
  • 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 [hep-ph] (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 [hep-ph])
      - 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