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

Last edited by Martin Gabelmann Jun 28, 2019
Page history

Diphoton_and_digluon_vertices_with_SPheno

Diphoton and digluon vertices with SPheno

For the calculation of the partial width of a neutral scalar Φ decaying into two gluons or two photons SPheno follows closely for the LO and NLO contributions. The partial widths at LO are given by

\Gamma(\Phi \to \gamma \gamma)_{\rm LO} = \frac{G_F \alpha^2(0) m_\Phi^3}{128 \sqrt{2} \pi^3} \Bigg|\sum_f N^f_c Q_f^2 r^\Phi_f A_f(\tau_f) + \sum_s N^s_c r^\Phi_s Q_s^2 A_s(\tau_s) + \sum_v N^v_c r^\Phi_v Q_v^2 A_v(\tau_v) \Bigg|^2, \\ \Gamma(\Phi \to g g)_{\rm LO} = \frac{G_F \alpha_s^2(\mu) m_\Phi^3}{36 \sqrt{2} \pi^3} \Bigg|\sum_f \frac{3}{2} D_2^f r^\Phi_f A_f(\tau_f) + \sum_s \frac{3}{2} D_2^s r^\Phi_s A_s(\tau_s) + \sum_v \frac{3}{2} D_2^v r^\Phi_v A_v(\tau_v) \Bigg|^2.

Here, the sums are over all fermions f, scalars s and vector bosons v which are charged or coloured and which couple to the scalar Φ. Q is the electromagnetic charges of the fields, Nc are the colour factors and D2 is the quadratic Dynkin index of the colour representation which is normalised to \frac12 for the fundamental representation. We note that the electromagnetic fine structure constant α must be taken at the scale μ = 0, since the final state photons are real . In contrast, αs is evaluated at μ = mΦ which, for the case of interest here, is μ = 750 GeV. riΦ are the so-called reduced couplings, the ratios of the couplings of the scalar Φ to the particle i normalised to SM values. These are calculated as

r^\Phi_f = \frac{v}{2 M_f} (C^L_{\bar f f \Phi}+C^R_{\bar f f \Phi}), \\ r^\Phi_s = \frac{v}{2 M^2_s} C_{s s^\* \Phi},\\ r^\Phi_v = -\frac{v}{2 M^2_v} C_{v v^\* \Phi}.

Here, v is the electroweak VEV and C are the couplings between the scalar and the different fields with mass Mi (i = f, s, v). Furthermore,

\tau_x = \frac{m_\Phi^2}{4 m_x^2}

holds and the loop functions are given by

A_f = 2 (\tau + (\tau -1) f(\tau))/\tau^2, \\ A_s = -(\tau-f(\tau))/\tau^2, \\ A_v = -(2 \tau^2 + 3\tau + 3 (2 \tau -1) f(\tau) )\tau^2,

with

f(\tau) = \begin{cases} \text{arcsin}^2 \sqrt{\tau} \hspace{1cm} \text{for} \,\, \tau \le 1,\\ -\frac{1}{4}\left(\log \frac{1+\sqrt{1-\tau^{-1}}}{1-\sqrt{1-\tau^{-1}}} -i\pi \right)^2 \text{for} \,\, \tau gt; 1. \end{cases}

For a pure pseudo-scalar state only fermions contribute, i.e. the LO widths read

\Gamma(A \to \gamma \gamma)_{\rm LO} = \frac{G_F \alpha^2 m_A^3}{32 \sqrt{2} \pi^3} \left|\sum_f N^f_c Q_f^2 r^A_f A^A_f(\tau_f) \right|^2, \\ \Gamma(A \to g g)_{\rm LO} = \frac{G_F \alpha_s^2 m_A^3}{36 \sqrt{2} \pi^3} \left|\sum_f 3 D_2^f r^A_f A^A_f(\tau_f) \right|^2,

where

AfA = f(τ)/τ ,

and rfA takes the same form as rfΦ in , simply replacing Cf̄fΦL, R by Cf̄fAL, R.

These formulae are used by SPheno to calculate the full LO contributions of any CP-even or odd scalar present in a model including all possible loop contributions at the scale μ = mΦ. However, it is well known, that higher order corrections are important. Therefore, NLO, NNLO and even N3LO corrections from the SM are adapted and used for any model under study. In case of heavy colour fermionic triplets, the included corrections for the diphoton decay are

r^\Phi_f \to r_f \left(1 - \frac{\alpha_s}{\pi} \right), \\ r^\Phi_s \to r_s \left(1+\frac{8 \, \alpha_s}{3 \pi} \right).

These expressions are obtained in the limit τf → 0 and thus applied only when mΦ < mf. rfA does not receive any corrections in this limit. For the case mΦ > 100mf, we have included the NLO corrections in the light quark limit given by

r^X_f \to r^X_f \left(1+ \frac{\alpha_s}{\pi} \left[-\frac{2}{3} \log 4\tau + \frac{1}{18} \left(\pi^2 -\log^2 4\tau\right) + 2\log \left( \frac{\mu_{\text{NLO}}^2}{m_f^2} \right) + i \frac{\pi}{3}\left(\frac{1}{3} \log 4\tau +2 \right) \right] \right)

for X = Φ, A. μNLO is the renormalisation scale used for these NLO corrections, chosen to be μNLO = mΦ/2. In the intermediate range of 100mf > mΦ > 2mf, no closed expressions for the NLO correction exist. Our approach in this range was to extract the numerical values of the corrections from HDECAY and to fit them. For the digluon decay rate, the corrections up to N3LO are included and parametrised by

\Gamma(X \to g g) = \Gamma(X \to g g)_{\rm LO} \left(1 + C_X^{\rm NLO} + C_X^{\rm NNLO} + C_X^{\rm \text{N}^3LO} \right)\,,

with

C_\Phi^{\rm NLO} = \left(\frac{95}{4} - \frac76 N_F \right) \frac{\alpha_s}{\pi}\,, \\ C_\Phi^{\rm NNLO} = \Bigg(370.196 + (-47.1864 + 0.90177 N_F) N_F + (2.375 + 0.666667 N_F)\log \frac{m_\Phi^2}{m_t^2}\Bigg) \frac{\alpha^2_s}{\pi^2}\,, \\\ C_\Phi^{\rm \text{N}^3LO} = \left(467.684 + 122.441 \log \frac{m_\Phi^2}{m_t^2} + 10.941 \left(\log \frac{m_\Phi^2}{m_t^2}\right)^2 \right) \frac{\alpha^3_s}{\pi^3}\,,

and

C_A^{\rm NLO} = \left(\frac{97}{4} - \frac76 N_F \right) \frac{\alpha_s}{\pi}\,, \\ C_A^{\rm NNLO} = \left(171.544 + 5 \log \frac{m_\Phi^2}{m_t^2}\right) \frac{\alpha^2_s}{\pi^2}

For pseudoscalar we include only corrections up to NNLO as the \rm \text{N}^3LO are not known for CP-odd scalars.

See also

References

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