Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
SARAH SARAH
  • Project overview
    • Project overview
    • Details
    • Activity
  • Packages & Registries
    • Packages & Registries
    • Container Registry
  • Analytics
    • Analytics
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Members
    • Members
  • Activity
Collapse sidebar
  • GOODSELL Mark
  • SARAHSARAH
  • Wiki
  • Using_SPheno_for_two loop_masses

Last edited by Martin Gabelmann Jun 28, 2019
Page history

Using_SPheno_for_two loop_masses

Using SPheno for two-loop masses

SARAH automatically writes the necessary Fortran routines in the SPhenooutput to calculate the two-loop Higgs corrections. There are five flags which can be used in the Les Houches input file to adjust the properties of the two-loop calculation

Block SPhenoInput #
...
  7  ... # Skip two loop masses
  8  ... # Choose two-loop method
  9  ... # Gaugeless limit
 10  ... # Safe mode
...      #
400  ... # Step-size for purely-numerical method
401  ... # Step-size for semi-analytical method

The following values are possible:

  • SPhenoInput[7]:
    • 0: Don’t skip two-loop masses
    • 1: Skip two-loop masses
  • SPhenoInput[8]:
    • 1: Effective potential two-loop calculation with purely numerical derivation
    • 2: Effective potential two-loop calculation with analytical derivation of loop functions
    • 3: Diagrammatic two-loop calculation functions (default)
    • 9: Use hardcoded routines available for some models
  • SPhenoInput[9]:
    • 0: Turn off gauge-less limit
    • 1: Use gauge-less limit (default)
  • SPhenoInput[10]:
    • 0: Turn off the safe-mode (default)
    • 1: Use safe-mode
  • SPhenoInput[400]: a real number (default: 0.5)
  • SPhenoInput[401]: a real number (default: 0.001)

Note that the hard-coded two-loop routines are not included by default in the SPhenooutput of SARAH. To include them, in the SPheno.m file the flag

Use2LoopFromLiterature = True;

has to be set.

The flags SPhenoInput[400] and SPhenoInput[401] can be used to check numerical stability of the derivation. If the step size is choosen to be too small or large the numerical derivation might suffer from some instabilites. We found that the initial step size for derivation with the fully numerical method usually needs a larger initial step-size for the considered VEV especially for heavy SUSY spectra. The reason is that the potential is of O(MSUS**Y4) and the overall value only changes slightly when the VEVs are varied. The second method usually operates acceptably with a smaller initial step size because objects of at most order O(MSUS**Y2, MZ2) are derived numerically. In addition, we make the approximation that in the purely numerical approach mass squareds in the loop which are smaller than 10−5 times the largest mass squared in the loop is taken to be zero. For the semi-analytical approach we take a limit of 10−8.

Safe-Mode and beyond

For the effective potential calculation, there is a “safe mode” for SPheno via flag 10: in this case SPhenostarts with a large initial step size which is decreased automatically. It checks for what range of the initial step size the results are numerically stable by comparing the results obtained with different inital step sizes. If no stable range is found an error is returned. In the unlikely case that both methods suffer from numerical instabilities the chances are good that the diagrammatic works. Another, less recommended, possibility is to increase the numerical precision by passing from double to quadruple precision. For this purpose, the Makefile located in SPheno-3.3.0/src must be changed. The line

PreDef = -DGENERATIONMIXING  -DONLYDOUBLE

should be replaced by

PreDef = -DGENERATIONMIXING  -DQUADRUPOLE

Afterwards the entire SPhenocode must be recompiled via

 make cleanall
 make Model=$MODEL

This will slow down the numerical evaluation significantly. However, have not found an example for resonable SUSY masses where this was necessary.

See also

Clone repository

Home

Index

  • Additional terms in Lagrangian
  • Advanced usage of FlavorKit
  • Advanced usage of FlavorKit to calculate new Wilson coefficients
  • Advanced usage of FlavorKit to define new observables
  • Already defined Operators in FlavorKit
  • Already defined observables in FlavorKit
  • Auto-generated templates for particles.m and parameters.m
  • Automatic index contraction
  • Basic definitions for a non-supersymmetric model
  • Basic definitions for a supersymmetric model
  • Basic usage of FlavorKit
  • Boundary conditions in SPheno
  • CalcHep CompHep
  • Calculation of flavour and precision observables with SPheno
  • Checking the particles and parameters within Mathematica
  • Checks of implemented models
  • Conventions
  • Decay calculation with SPheno
  • Defined FlavorKit parameters
  • Definition of the properties of different eigenstates
  • Delete Particles
  • Different sets of eigenstates
  • Diphoton and digluon vertices with SPheno
  • Dirac Spinors
  • FeynArts
  • Fine-Tuning calculations with SPheno
  • Flags for SPheno Output
  • Flags in SPheno LesHouches file
  • FlavorKit
  • FlavorKit Download and Installation
  • Flavour Decomposition
  • GUT scale condition in SPheno
  • Gauge Symmetries SUSY
  • Gauge Symmetries non-SUSY
  • Gauge fixing
  • Gauge group constants
  • General information about Field Properties
  • General information about model implementations
  • Generating files with particle properties
  • Generic RGE calculation
  • Global Symmetries SUSY
  • Global Symmetries non-SUSY
  • Handling of Tadpoles with SPheno
  • Handling of non-fundamental representations
  • HiggsBounds
  • Higher dimensionsal terms in superpotential
  • Input parameters of SPheno
  • Installation
  • Installing Vevacious
  • LHCP
  • LHPC
  • LaTeX
  • Lagrangian
  • Loop Masses
  • Loop calculations
  • Loop functions
  • Low or High scale SPheno version
  • Main Commands
  • Main Model File
  • Matching to the SM in SPheno
  • MicrOmegas
  • ModelOutput
  • Model files for Monte-Carlo tools
  • Model files for other tools
  • Models with Thresholds in SPheno
  • Models with another gauge group at the SUSY scale
  • Models with several generations of Higgs doublets
  • More precise mass spectrum calculation
  • No SPheno output possible
  • Nomenclature for fields in non-supersymmetric models
  • Nomenclature for fields in supersymmetric models
  • One-Loop Self-Energies and Tadpoles
  • One-Loop Threshold Corrections in Scalar Sectors
  • Options SUSY Models
  • Options non-SUSY Models
  • Parameters.m
  • Particle Content SUSY
  • Particle Content non-SUSY
  • Particles.m
  • Phases
  • Potential
  • Presence of super-heavy particles
  • RGE Running with Mathematica
  • RGEs
  • Renormalisation procedure of SPheno
  • Rotations angles in SPheno
  • Rotations in gauge sector
  • Rotations in matter sector
  • SARAH in a Nutshell
  • SARAH wiki
  • SLHA input for Vevacious
  • SPheno
  • SPheno Higgs production
  • SPheno Output
  • SPheno and Monte-Carlo tools
  • SPheno files
  • SPheno mass calculation
  • SPheno threshold corrections
  • Setting up SPheno.m
  • Setting up Vevacious
  • Setting up the SPheno properties
  • Special fields and parameters in SARAH
  • Superpotential
  • Support of Dirac Gauginos
  • Supported Models
  • Supported gauge sectors
  • Supported global symmetries
  • Supported matter sector
  • Supported options for symmetry breaking
  • Supported particle mixing
  • Tadpole Equations
  • The renormalisation scale in SPheno
  • Tree-level calculations
  • Tree Masses
  • Two-Loop Self-Energies and Tadpoles
  • UFO
  • Usage of tadpoles equations
  • Using SPheno for two-loop masses
  • Using auxiliary parameters in SPheno
  • VEVs
  • Vertices
  • Vevacious
  • WHIZARD