





title: TwoLoop SelfEnergies and Tadpoles



permalink: /TwoLoop_SelfEnergies_and_Tadpoles/










[Category:Calculations](/Category:Calculations "wikilink")






General










SARAH can calculate twoloop contributions to the selfenergies of real scalars. There are two equivalent approaches implemented in the SPheno interface of SARAH to perform these calculations






1. an effective potential approach



2. a diagrammatic approach with vanishing external momenta






Because of the very complicated form of the results there is no output of the corresponding expressions in the Mathematic or LaTeX format but the results are just included in the [SPheno](/SPheno "wikilink") output. Since both calculations are based on a completely independent implementation and use a different approach they are very useful to perform cross checks. For the MSSM and NMSSM both calculations reproduce exactly the results obtained by widely used routines implemented in other tools. However, for nonminimal SUSY models are no references available to compare with. Thus, the only possibility to cross check the results is within SPheno and comparing the two different methods.






### Using the different methods






It is explained in [Using SPheno for twoloop masses](/Using_SPheno_for_twoloop_masses "wikilink") how the user can make use of the different possibilities to calculate twoloop masses.






Effective potential calculation










### Selfenergies and tadpoles from the effective potential in general






We consider a set of real scalars {*ϕ*<sub>*k*</sub>} which are diagonalised to physical states {*h*<sub>*k*</sub>}. The scalar potential at treelevel is *V*<sup>*T*</sup>({*ϕ*<sub>*k*</sub>}). The pole masses of *h*<sub>*i*</sub> are in general the eigenvalues of the loop corrected mass matrix given by






$\\frac{\\partial^2 V^T}{\\partial \\phi_i \\partial \\phi_j} + \\sum_n \\Pi^{(n)}_{ij}(p^2)$






Here, *Π*<sup>(*n*)</sup>(*p*<sup>2</sup>) is the selfenergy at the *n*loop level which usually depends on the external momenta. The parameters appearing in *V*<sup>*T*</sup> should be chosen to minimise the effective potential. For a pole mass *m*<sub>*h*<sub>*i*</sub></sub> the momentum is fixed to be *p*<sup>2</sup> = *m*<sub>*h*<sub>*i*</sub></sub><sup>2</sup>. To include this momentum dependence of *Π*<sup>(*n*)</sup> a diagrammatic calculation is necessary. However, things become significantly easier if one considers the limit of *p*<sup>2</sup> = 0. In this limit *Π*<sup>(*n*)</sup>(0) is equivalent to the second derivative of the effective potential at the *n*loop level,






$\\label{eq:effpotApprox}



\\Pi_{ij}^{(n)}(0) = \\frac{\\partial^2 V^{(n)}}{\\partial \\phi_i \\partial \\phi_j} \\, .$






SARAH/SPheno uses a hybrid method for the calculation of the scalar masses at the twoloop level in the *D**R*′ renormalization scheme: while the oneloop corrections are calculated including the full momentum dependence, the twoloop corrections are derived in the effective potential approximation.






### Calculation of the twoloop self energies in the effective potential approximation






CPviolating is neglected in the following. Therefore, a set of neutral, complex scalars *H*<sub>*i*</sub> are decomposed after symmetry breaking as






$H_i = \\frac{1}{\\sqrt{2}}\\left(\\phi_i + i \\sigma_i + v_i \\right) \\, .$






*ϕ*<sub>*i*</sub> are the CPeven components, *σ*<sub>*i*</sub> the CPodd ones and *v*<sub>*i*</sub> are the vacuum expectation values (VEVs). Under this assumption, the twoloop corrections to the mass matrix of real scalars in the effective potential approach are given by






$\\label{eq:Pi2L}



\\Pi^{(2)}_{ij}(0) = \\frac{\\partial^2 V^{(2)}}{\\partial v_i \\partial v_j}\\,.$






In addition, one requires the tadpole contributions which are the first derivatives of the effective potential,






$\\label{eq:Ti2L}



\\delta t_i^{(2)} = \\frac{\\partial V^{(2)}}{\\partial v_i} \\,.$






The tadpoles then also contribute to the scalar masses by shifting the parameters; for example if we treat the vevs *v*<sub>*i*</sub> as fixed, “allloop” correct values, then one can exchange them for scalar masssquared at a given loop order via the tadpole equations, and in turn find compact equations for the total shift in the masssquareds of






$\\begin{aligned}



\\big( \\Delta \\mathcal{M}^2_S \\big)_{ij}^{\\rm eff} =& \\frac{\\partial^2 V^{(2)}}{\\partial v_i \\partial v_j} \\frac{\\delta_{ij}}{v_i} \\frac{\\partial V^{(2)}}{\\partial v_i },\\end{aligned}$ However, this is not done by SARAH which instead solves the tadpole equations and then uses the parameters derived from these in the treelevel mass calculation. The reason is that SARAHallows for a more general choice of variables to solve for via the tadpole equations.






*V*<sup>(2)</sup> receives contributions from the possible topologies:






[rahmenloszentriert](/Datei:Diagrams.png "wikilink")






The generic results for all of these diagrams in Landau gauge have been presented in Ref. and SARAH makes heavily use of these results. However, not all topologies are implemented so far but *V**V*, *V**V**S*, *V**V**V*, *S**V*, and *G**G**V* are neglected. The reason is that these only lead to nonzero contributions if massive vector bosons are involved. Working in the gaugeless limit, all contributions stemming from broken gauge groups are neglected so far because of the many complications that they entail.






The remaining topologies are those which have been considered so far in the MSSM to obtain the dominant twoloop corrections and which are implemented in public computer tools. The impact on the Higgs mass is considered to be moderate when neglecting contributions propotional to the electroweak gauge couplings in models with the SM gauge sector, since they are significantly smaller than the strong and top Yukawa coupling – and the diagrams involving these couplings are of typically lower multiplicity. In fact, by inspecting the form of the loop functions it is evident that the topologies *S**V*, *V**V*, *V**V**V*, *V**V**S* < *m**a**t**h* > *a**n**d* < *m**a**t**h* > *G**G**V* lead to contributions proportional to *g*<sub>*E**W*</sub><sup>4</sup>*v*<sup>2</sup> < *m**a**t**h* > *w**h**e**r**e* < *m**a**t**h* > *g*<sub>*E**W*</sub> stands for the electroweak gauge coupling; these are hence subdominant to the contributions from $SSV, FFV, \\overline{FF}V$ diagrams involving the electroweak gauge bosons (which we are also neglecting) which would be proportional to *g*<sub>*E**W*</sub><sup>2</sup>*v*<sup>2</sup> < *m**a**t**h* > .*N**o**t**e*, *f**o**r**c**o**n**s**i**s**t**e**n**c**y*, *S**A**R**A**H**a**l**s**o**s**e**t**s**a**l**l**g**a**u**g**e**c**o**u**p**l**i**n**g**s**a**r**i**s**i**n**g**i**n**t**h**e* < *m**a**t**h* > *D*terms of broken groups in the vertices to zero. The masses used in the loops are treelevel masses calculated from running *D**R*′ parameters. In this context, there are two possibilities to treat *D*term contributions to the treelevel masses: either one can work in the gaugeless limit where these contributions are also put to zero, or one can work with the full treelevel masses as this is also done in diagrammatic calculations. SARAH offers both possibilities. In the gaugeless limit the *S**U*(2)<sub>*L*</sub> × *U*(1)<sub>*Y*</sub> gauge symmetries become global symmetries. Therefore, it is obvious that no gauge dependence has been introduced by including only Goldstone diagrams but no diagrams of the corresponding vector bosons. Of course, when including *D*terms in the mass matrices the derivatives of the *D*terms are still forced to vanish.






For the calculation of the twoloop effective potential the expressions of Ref. given in twocomponent notation have been translated into fourcomponent language. All necessary generic expressions have been implemented in SARAH. SARAHuses these expressions to generate `Fortran`code for all twoloop diagrams which are possible in the considered model assuming the topologies *S**S*, *S**S**V*, *S**S**S*, *F**F**S*, $\\overline{FF}S$, *F**F**V* and $\\overline{FF}V$.






### Numerical derivation






As soon as the twoloop effective potential is calculated, one can obtain the twoloop corrections to the Higgs mass by performing the derivative to calculate the tadpoles and selfenergies numerically. The numerical derivation in SPhenois done by using Ridders’ method of polynomial extrapolation with dynamical stepsize. SARAH provides two different methods to perform this derivation:






1. *'Purely numerical derivation:* in this approach the entire potential is derived with respect to the VEVs.



2. **Semianalytical derivation:** in this approach the derivatives of all masses and couplings with respect to the VEVs are calculated separately also in a numerical way. However, all one and twoloop derivatives of the loopfunctions with respect to their arguments have been calculated analytically and implemented in the output SPhenocode. The derivatives of the potential are then calculated combining both results using the chain rule. This can be easily done because every contribution to the potential is a product of couplings (*c*<sub>1</sub>, *c*<sub>2</sub>), masses (*m*<sub>*i*</sub>) and a loop function *f*<sub>*X*</sub> with a coefficient *k*:






$\\begin{aligned}



V^{(2)}_X &=& k\\cdot (c_1 c_2) \\cdot f_{X}(m_1^2,m_2^2,m_3^2), \\quad \\text{for } X=FFS,FFV,SSV\\\\



V^{(2)}_X &=& k\\cdot (c_1 c_2) \\cdot m_{F1} m_{F2} \\cdot f_{X}(m_{F1}^2,m_{F2}^2,m_3^2), \\quad \\text{for } X=\\overline{FFS},\\overline{FFV}\\\\



V^{(2)}_X &=& k\\cdot (c_1) \\cdot f_{X}(m_1^2,m_2^2), \\quad \\text{for } X=SS\\end{aligned}$






The second method is numerically slightly more expensive but it is also more stable. In particular, in the presence of large hierarchies in the VEVs the purely numerical method could become inaccurate.






Diagrammatic calculation










A fully diagrammatic calculation for twoloop contributions to scalar selfenergies with SARAH –SPheno became available with Ref. . In this setup a new set of generic expressions is used. The following twoloop diagrams are included in the limit *p*<sup>2</sup> = 0:






[zentriert](/Datei:2loop_diagrammatic.png "wikilink")






These are again the diagrams which don’t vanish in general in the gaugeless limit. The results of Ref. have the advantage that the expressions which are derived from the effective potential are much simpler than taking the limit *p*<sup>2</sup> → 0 in other twoloop functions available in literature. The diagrammatic method gives completely equivalent results to the effective potential calculation but is usually numerically more robust and it is now the default method. [One can switch in SPheno easily between the diagrammatic and effective potential calculation.](/Using_SPheno_for_twoloop_masses "wikilink")






### CP Violation






The diagrammatic calculation can also be used to get the twoloop corrections to real scalars which don't receive a VEV and for models with CP violation in the Higgs sector . No modifications are required: for the computation of masses for CPodd scalars in CPconserving theories and for models with CP violation exactly the same routines apply because the results of Ref. are given in terms of real scalars, and the CPodd scalars are just CPeven scalars with different labels.






Handling of the Goldstones










It is well known that the derivatives of the effective potential in the Landau gauge may suffer from divergences due to massless Goldstone bosons : the derivatives of some loop functions have infrared singularities. In the gaugeless limit of the MSSM, this problem is circumvented because the masses of the pseudoscalars become independent of the Higgs vevs, and so the derivatives of the effective potential do not contain any singular functions. However, once one goes beyond the MSSM (even to the NMSSM) this problem reappears. The singularities in the first derivative of the potential may be tamed by resummation methods, but the second derivatives may remain singular and it has been suggested that the problem could be resolved by passing to the pole mass calculation instead . However, working in the gaugeless limit this problem is usually not present: as stated above, the treelevel masses entering the calculation are calculated at the minimum of (full) treelevel potential. Thus, using the Lagrangian parameters at this minimum in the mass matrices expressed in the gaugeless limit, the Goldstone masses are usually nonzero and no divergences in the numerical evaluation show up. This is one strong motivation for working in the gaugeless limit: it reduces the dependence of the treelevel Goldstone \`\`masses'' that enter into the calculation on the renormalisation scale.






See also










 [Loop functions](/Loop_functions "wikilink")



 [Calculation of the mass spectrum with SPheno](/SPheno_mass_calculation "wikilink")



 [Using SPheno for twoloop masses](/Using_SPheno_for_twoloop_masses "wikilink")






References



 


\ No newline at end of file 