Hi SPheno users,
I want to generate some spectra of a model linked to the NMSSM. But first, I've tried to generate a simple NMSSM spectrum using a well-known benchmark (https://twiki.cern.ch/twiki/pub/LHCPhys ... p.slha.txt ).
I've downloaded the SARAH Output for SPheno (https://sarah.hepforge.org/trac/wiki/NMSSM) and used it with SPheno 4.0.3. with the input parameters define in the benchmark.
First, I had some errors because the input file generate by SARAH forgot to put the FLAG 3 to 1 in the MODSEL block (corresponding to the NMSSM model while the FLAG 3 = 0 corresponds to the MSSM).
But now I have several tachyonic states that mustn’t be there. (the generate spectrum from another generator is given https://twiki.cern.ch/twiki/pub/LHCPhys ... r.slha.txt ).
Where is the problem ?
I have some problems to put my input file in attachments so here is my input file :
Block MODSEL #
1 1 # 1/0: High/low scale input
2 1 # Boundary Condition
3 1
6 1 # Generation Mixing
Block SMINPUTS # Standard Model inputs
2 1.166370E-05 # G_F,Fermi constant
3 1.187000E-01 # alpha_s(MZ) SM MSbar
4 9.118870E+01 # Z-boson pole mass
5 4.180000E+00 # m_b(mb) SM MSbar
6 1.735000E+02 # m_top(pole)
7 1.776690E+00 # m_tau(pole)
Block MINPAR # Input parameters
1 0.10789E+04 #m0
2 0.98607E+03 #m12
3 0.22663E+01 #tanb
5 -0.11821E+04 #A0
Block EXTPAR # Input parameters
61 0.64379E+00 #lambda
62 0.35121E+00 #kappa
63 -0.31170E+03 #Alambda
64 1E+02 #Akappa
65 0.17864E+03 #Mueff
Block SPhenoInput # SPheno specific input
1 -1 # error level
2 0 # SPA conventions
7 0 # Skip 2-loop Higgs corrections
8 3 # Method used for two-loop calculation
9 1 # Gaugeless limit used at two-loop
10 0 # safe-mode used at two-loop
11 1 # calculate branching ratios
13 1 # 3-Body decays: none (0), fermion (1), scalar (2), both (3)
14 0 # Run couplings to scale of decaying particle
12 1.000E-04 # write only branching ratios larger than this value
15 1.000E-30 # write only decay if width larger than this value
31 -1 # fixed GUT scale (-1: dynamical GUT scale)
32 0 # Strict unification
34 1.000E-04 # Precision of mass calculation
35 40 # Maximal number of iterations
36 5 # Minimal number of iterations before discarding points
37 1 # Set Yukawa scheme
38 2 # 1- or 2-Loop RGEs
50 1 # Majorana phases: use only positive masses (put 0 to use file with CalcHep/Micromegas!)
51 0 # Write Output in CKM basis
52 1 # Write spectrum in case of tachyonic states
55 1 # Calculate loop corrected masses
57 1 # Calculate low energy constraints
65 1 # Solution tadpole equation
75 1 # Write WHIZARD files
76 1 # Write HiggsBounds file
86 0. # Maximal width to be counted as invisible in Higgs decays; -1: only LSP
510 0. # Write tree level values for tadpole solutions
515 0 # Write parameter values at GUT scale
520 1. # Write effective Higgs couplings (HiggsBounds blocks): put 0 to use file with MadGraph!
525 0. # Write loop contributions to diphoton decay of Higgs
530 1. # Write Blocks for Vevacious
EDIT : I forgot to put the content of the Messages.out file so here's it (and sorry for the long post...) :
Problem while reading SPhenoInput, ignoring unknown entry 14 0.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 15 1.0000000000000001E-030
Problem while reading SPhenoInput, ignoring unknown entry 50 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 51 0.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 52 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 55 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 57 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 65 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 75 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 76 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 86 0.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 510 0.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 515 0.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 520 1.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 525 0.0000000000000000
Problem while reading SPhenoInput, ignoring unknown entry 530 1.0000000000000000
Warning from routine SfermionMass1MSSM!
in the calculation of the masses
occurred a negative mass squared!!!
msf2 -109.31373024560503 2928.1202248817544
M_L2, M_R2, Af, mu, vevs, Yuk, T3, Yl, Yr, g, gp
0.0000000000000000 0.0000000000000000 ( 0.0000000000000000 , 0.0000000000000000 ) ( 178.63999999999999 , 0.0000000000000000 ) 97.631103550525737 221.26136997655649 ( 4.1315489095397430E-002, 0.0000000000000000 ) -0.50000000000000000 0.33333333333333331 0.66666666666666663 0.66515045721523214 0.35533125800232790
The error occured in the following chain of routines
SPheno3
Model_NMSSM
TreeMassesNMSSM
SfermionMass3MSSM
SfermionMass1MSSM
Warning from routine SfermionMass1MSSM!
in the calculation of the masses
occurred a negative mass squared!!!
msf2 -109.31373024560503 2928.1202248817544
M_L2, M_R2, Af, mu, vevs, Yuk, T3, Yl, Yr, g, gp
0.0000000000000000 0.0000000000000000 ( 0.0000000000000000 , 0.0000000000000000 ) ( 178.63999999999999 , 0.0000000000000000 ) 97.631103550525737 221.26136997655649 ( 4.1315489095397430E-002, 0.0000000000000000 ) -0.50000000000000000 0.33333333333333331 0.66666666666666663 0.66515045721523214 0.35533125800232790
The error occured in the following chain of routines
SPheno3
Model_NMSSM
TreeMassesNMSSM
SfermionMass1MSSM
Problem in tqli, too many iterations
Routine tqli, too many iterations: 1
Routine ScalarMassNMSSMeff: mP0^2_1 < 0 1
Routine ScalarMassNMSSMeff: mH+^2 < 0 1
Routine SfermionMass1MSSM: negative mass squared 2
Routine SfermionMass3MSSM: negative mass squared 4
Routine TreeMassesNMSSM: negative sneutrino mass squared 3
Best regards,
Robin Ducrocq
Spectrum generation (NMSSM) : tachyonic states with a well-known benchmark
Re: Spectrum generation (NMSSM) : tachyonic states with a well-known benchmark
Dear Robin, I believe that the files that you downloaded from https://sarah.hepforge.org/trac/wiki/NMSSM is pretty much outdated. My recommendation is to regenerate it using the latest SARAH version, that will produce an output for SPheno that works perfectly with SPheno 4.0.3.
Re: Spectrum generation (NMSSM) : tachyonic states with a well-known benchmark
Hi,
These errors look like running SPheno not SPhenoNMSSM. And, of course, as Avelino said, please generate the code using the a newer SARAH version.
Cheers
Florian
These errors look like running SPheno not SPhenoNMSSM. And, of course, as Avelino said, please generate the code using the a newer SARAH version.
Cheers
Florian
Re: Spectrum generation (NMSSM) : tachyonic states with a well-known benchmark
Hi,
Thanks for your answers.
Indeed, I have used SPheno and not SPhenoNMSSM. I have generated a new SPheno Output with the last version of SARAH (4.14.1) and tried to generate a spectrum (with SPhenoNMSSM) with the inputs define during the implementation of the model with SARAH... and it works (no tachyonic states problem).
But when I try to generate a spectrum with the benchmarks define in https://twiki.cern.ch/twiki/bin/view/LH ... markPoints (I'm interested in mSUGRA models so I have only tried benchmarks for mSUGRA), I have a problem of step size :
Problem with RGE running. Errorcode: -12
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP Subroutine TerminateProgram
Normally, those benchmarks must converge.
Here is my Messages.out file :
Problem in rkqs, stepsize underflow
kont -12 1.0000000000000000E-003 -16.811242831518264 0.33622485663036528
t_out 0.0000000000000000 1.0000000000000000E+020
ErrorLevel, Iname: -1 0
The error has occured in the following chain of subroutines:
FirstGuess
runRGE
odeintB
rkqs
rge218
ParametersToG218
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
Hopefully you find the error soon
Best regards,
Robin
Thanks for your answers.
Indeed, I have used SPheno and not SPhenoNMSSM. I have generated a new SPheno Output with the last version of SARAH (4.14.1) and tried to generate a spectrum (with SPhenoNMSSM) with the inputs define during the implementation of the model with SARAH... and it works (no tachyonic states problem).
But when I try to generate a spectrum with the benchmarks define in https://twiki.cern.ch/twiki/bin/view/LH ... markPoints (I'm interested in mSUGRA models so I have only tried benchmarks for mSUGRA), I have a problem of step size :
Problem with RGE running. Errorcode: -12
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP Subroutine TerminateProgram
Normally, those benchmarks must converge.
Here is my Messages.out file :
Problem in rkqs, stepsize underflow
kont -12 1.0000000000000000E-003 -16.811242831518264 0.33622485663036528
t_out 0.0000000000000000 1.0000000000000000E+020
ErrorLevel, Iname: -1 0
The error has occured in the following chain of subroutines:
FirstGuess
runRGE
odeintB
rkqs
rge218
ParametersToG218
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00\00
Hopefully you find the error soon
Best regards,
Robin
Re: Spectrum generation (NMSSM) : tachyonic states with a well-known benchmark
Hi,
there are two possible reasons for that:
1) The BPs are defined for other conventions (for instance, there are two common choices at which scale lambda and kappa are set: m_GUT or m_SUSY)
2) The BP is really on the edge and very fine-tuned, i.e. small differences how things are calculated (e.g. loop order in the matching, SM input parameters,...) can make a difference between valid/not valid. This is very often the case when one uses huge lambda values very close to the perturbativity limit (lambda ~0.65 at m_SUSY). Tweaking the parameters should help in that case.
Cheers,
Florian
there are two possible reasons for that:
1) The BPs are defined for other conventions (for instance, there are two common choices at which scale lambda and kappa are set: m_GUT or m_SUSY)
2) The BP is really on the edge and very fine-tuned, i.e. small differences how things are calculated (e.g. loop order in the matching, SM input parameters,...) can make a difference between valid/not valid. This is very often the case when one uses huge lambda values very close to the perturbativity limit (lambda ~0.65 at m_SUSY). Tweaking the parameters should help in that case.
Cheers,
Florian