Spectrum generation (NMSSM) : tachyonic states with a well-known benchmark
Posted: 16. Apr 2019, 14:13
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
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