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
  • Matching_to_the_SM_in_SPheno

Last edited by Martin Gabelmann Jun 28, 2019
Page history

Matching_to_the_SM_in_SPheno

Matching to the SM in SPheno

General

The parameters of a model have to be matched to the measured values of masses and couplings in the SM. The procedure to do this differs between a Low or High scale SPheno version:

  1. For a high-scale version one needs to match to the running M**S parameters of the SM
  2. For a low-scale version it is common to match to the measured (pole) values of the SM parameters

Matching to running parameters

Note

The method described here was introduce in SARAH 4.9.0 and replaces previous attempts to define the matching for

  • Models with another gauge group at the SUSY scale
  • Models with several generations of Higgs doublets

General

The SPheno versions by SARAH calculate the running values for the SM gauge (g1SM, g2SM, g3SM) and Yukawa (YeSM, YdSM, YuSM) couplings as well as for the electroweak VEV (vSM) including the full one-loop corrections in the present model. This calculation is done at the scale MZ and the input parameters which are used are

  • The running electroweak and strong coupling constants (α(MZ), αS(MZ))
  • The measured masses of Z- and W- bosons (MZ, MW)
  • The measured masses of all SM fermions (md,mu,ml)
  • The entries of the CKM matrix

The exact calculation of the thresholds is explained here: SPheno_threshold_corrections. The user can then define in the SPheno.m file the matching conditions between the running SM parameters and the parameters in the model via

DEFINITION[MatchingConditions]={
{Model parameters, SM parameter},
...
};

Matching for a One-Higgs doublet model

The simplest matching conditions are those in which the Higgs sector is the same as in the SM:

DEFINITION[MatchingConditions]={
{Ye, YeSM},
{Yd, YdSM},
{Yu, YuSM},
{g1, g1SM},
{g2, g2SM},
{g3, g3SM},
{v, vSM}
};

One can also use as short form

DEFINITION[MatchingConditions]=Default[OHDM];

In this case, SARAH makes use of the knowledge which parameters in the model correspond to the SM parameters.

Matching for a THDM-II-like model

The matching conditions for a model which has a Higgs sector as a THDM-II like the most common SUSY models for instance, would read with this definition

DEFINITION[MatchingConditions]={
{g1, g1SM},
{g2, g2SM},
{g3, g3SM},
{vd, vSM/Sqrt[1+TanBeta^2]},
{vu, TanBeta vd},
{Yu, vSM YuSM/vu},
{Yd, vSM YdSM/vd},
{Ye, vSM YeSM,vd}
};

One can also use as short form

DEFINITION[MatchingConditions]=Default[THDMII];

In this case, SARAH makes use of the knowledge which parameters in the model correspond to the SM parameters.

Matching for Models with extra gauge groups

One can use these conditions also to define the matching conditions in a model with has a S**U(2)R × U(1)B − L gauge group which just gets broken at the SUSY scale to U(1)Y. The matching conditions in this case read

DEFINITION[MatchingConditions]={
 {gR, Sqrt[5/3]*g1RBLFactor*g1SM},
 {gBL, -((-g1SM^2 gRBL + gBLR gR gRBL + Sqrt[g1SM^2 (gBLR - gR)^2 (-g1SM^2 + gR^2 + gRBL^2)])/(g1SM^2 - gR^2))},
 {g2, g2SM},
 {g3, g3SM},
 {vd, vSM/Sqrt[1 + TanBeta^2]},
 {vu, VEVSM1*TanBeta},
 {Yu, YuSM vSM/vu},
 {Yd, YdSM vSM/vd},
 {Ye, YeSM vSM/vd}
};

Here, g1RBLFactor which relates gR and gYS**M is in principle a free parameter which however can be fixed by imposing gauge coupling unification.

Matching for Models with extra Higgs doublets

Even much more complicated matching conditions can be imposed. For instance, the matching in the left-right symmetry model with two generations of bi-doublets are

v_{u,1} = v_L \sin\beta \sin\beta_u\\ v_{d,1} = v_L \cos\beta \sin\beta_d\\ v_{u,2} = v_L \sin\beta \cos\beta_u \\ v_{d,2} = v_L \cos\beta \cos\beta_d\\ Y_{Q_1} =-\frac{Y_d\sqrt{1+\tan\beta_d^2} -Y_u\sqrt{1+ \tan\beta_u^2} }{\tan\beta_d-\tan\beta_u}\\ Y_{Q_2} =\frac{\tan\beta_u Y_d\sqrt{1+\tan\beta_d^2} -Y_u\tan\beta_d\sqrt{1+ \tan\beta_u^2}}{\tan\beta_d-\tan\beta_u}\\ Y_{L_1} =-\frac{Y_e\sqrt{1+\tan\beta_d^2} -Y_\nu\sqrt{1+ \tan\beta_u^2} }{\tan\beta_d-\tan\beta_u}\\ Y_{L_2} =\frac{\tan\beta_u Y_e\sqrt{1+\tan\beta_d^2} -Y_\nu\tan\beta_d\sqrt{1+ \tan\beta_u^2}}{\tan\beta_d-\tan\beta_u}

where three angles are introduced to parametrize all VEVs and the neutrino Yukawa coupling is considered to be given as input. These matching conditions together with the ones in the gauge sector are defined via

(* to get simpler expressions *)
vdAux = vSM/Sqrt[1+TanBeta^2];
vuAux = TanBeta*vSM/Sqrt[1+TanBeta^2];

YuA=Transpose[YuSM]*Sqrt[1+TanBeta^2]/TanBeta
YdA=Transpose[YdSM]*Sqrt[1+TanBeta^2]
YeA=Transpose[YeSM]*Sqrt[1+TanBeta^2]

DEFINITION[MatchingConditions]={
 (* Matching for gauge couplings *)
 {gR, Sqrt[5/3]*g1RBLFactor*g1SM},
 {gBL, (gR g1SM)/Sqrt[gR^2 - g1SM^2]},
 {g2, g2SM},
 {g3, g3SM},

 (* Matching for VEVs *)
 {vPhiu[1] , (TanBeta*TanBetaU*vSM)/(Sqrt[1 + TanBeta^2]*Sqrt[1 + TanBetaU^2])},
 {vPhid[1] , (TanBetaD*vSM)/(Sqrt[1 + TanBeta^2]*Sqrt[1 + TanBetaD^2])},
 {vPhiu[2] , (TanBeta*vSM)/(Sqrt[1 + TanBeta^2]*Sqrt[1 + TanBetaU^2])},
 {vPhid[2] , vSM/(Sqrt[1 + TanBeta^2]*Sqrt[1 + TanBetaD^2])},

 (* Matching for Yukawas *)
 {YQ[1;;3,1;;3,1] , -(vuAux * vPhid[2] * YuA - vdAux * vPhiu[2] * YdA)/(vPhid[2] * vPhiu[1] - vPhid[1] * vPhiu[2])},
 {YQ[1;;3,1;;3,2] ,  (vuAux * vPhid[1] * YuA - vdAux * vPhiu[1] * YdA)/(vPhid[2] * vPhiu[1] - vPhid[1] * vPhiu[2])},
 {YL[1;;3,1;;3,1] ,  ( vdAux * vPhiu[2] *YeA)/(vPhid[2] * vPhiu[1] - vPhid[1] * vPhiu[2])},
 {YL[1;;3,1;;3,2] , (- vdAux * vPhiu[1] * YeA)/(vPhid[2] * vPhiu[1] - vPhid[1] * vPhiu[2])}
};

Matching to on-shell parameters

With version 4.12.0 the matching procedure has been unified for high and low-scale models by introducing a two-scale matching procedure, see 1703.03267. Thus, also for non-SUSY models the full one-loop thresholds are included by default. The definition of the matching conditions is the same as explained above and done via

DEFINITION[MatchingConditions]={
...

(The old approach of writing g1SM, etc. in BoundaryLowScaleInput is no longer recommended!)

However, under some circumstances, the user might be interested in using OS values for the input parameters. That's for instance the case when performing only a tree-level calculation of the scalar masses under the assumption that all higher order corrections can be absorbed into a set of counter-terms. In this case, one can use

Block SPhenoInput
66 0 # Two Scale Matching

in the Les Houches input file. This turns off the new two-scale matching routines. For a low-scale input, this results in a pole mass matching.

Remark

Note, the often used approach to ignore higher order corrections in the scalar sector of non-supersymmetric models should be taken with care as discussed in these papers:

  • N-loop running should be combined with N-loop matching
  • Perturbativity Constraints in BSM Models

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