Hello SARAH staff,

I've been having trouble when running MakeSPheno[] in a customized version of the MRSSM. The tadpole calculation takes forever (I have not waited until it ends), which must be because in SPheno.m I solve for vevs of the triplet T and singlet S

ParametersToSolveTadpoles = {mHd2,mHu2,vT, vS};

I'm aware that this is usually harder than solving for soft masses {mS2, mT2}. But in a .nb file I have been able to get analytical solutions for vS and vT by eliminating mS2 and mT2 via TadpoleEquation[]. Is there a way of telling SARAH to use my tadpole solutions and proceed to finish execution of MakeSPheno[]? My custom SPheno.m is pasted below.

Thanks,

===============================================================================
(*A customized SPheno.m that incorporates B-terms for the adjoints and solves for vS and vT*)

OnlyLowEnergySPheno = True;

(*Free MSSM parameters*)
MINPAR={
{3,TanBeta}
};

(*New free parameters absent/different than in the MSSM*)
EXTPAR = {
(*Dirac gaugino masses*)
{1, MDBinput},
{2, MDWinput},
{3, MDOinput},

(*R soft masses*)
(*Free parameters because <R>=0, i.e. no min. conds. *)
{10, mRd2input},
{11, mRu2input},

(*New mu terms*)
{40, MuDinput},
{41, MuUinput},

(*B-terms*)
{45, BmuInput},
{46, BsInput}, (* ! added by me *)
{47, BtInput}, (* ! added by me *)

(*SUSY R-H-adj Yukawas*)
{80, LSDinput},
{81, LSUinput},
{82, LTDinput},
{83, LTUinput},

(*Adjoint soft masses *)
{91, mS2input},
{92, mT2input},
{93, moc2input},

(*Sfermion soft masses*)
{101, mq2input},
{102, mq2entry33input},
{103, ml2input},
{104, ml2entry33input},
{105, mu2input},
{106, mu2entry33input},
{107, md2input},
{108, md2entry33input},
{109, me2input},
{1010, me2entry33input}
};

RealParameters = {TanBeta};

AssumptionsTadpoleEquations = {conj -> Identity};
ParametersToSolveTadpoles = {
mHd2,
mHu2,
vT, (* ! added by me *)
vS (* ! added by me *)
};

RenormalizationScaleFirstGuess = 1000.0^2;
RenormalizationScale = 1000.0^2;

(*Forbidden by R-symmetry*)
BoundaryHighScale = {
{B[MuD], 0},
{B[MuU], 0},
{\[Mu], 0}
};

UseBoundarySUSYatEWSB = True;

BoundaryEWSBScale={

{vd,Sqrt[4 mz2/(g1^2+g2^2)]*Cos[ArcTan[TanBeta]]},
{vu,Sqrt[4 mz2/(g1^2+g2^2)]*Sin[ArcTan[TanBeta]]},

(*Copied from Master.nb*)
{ vS, -((Sqrt[2]*(LamSD*LamTD*vd^2 + LamSU*LamTU*vu^2)*(LamTD*MuD*vd^2 + LamTU*MuU*vu^2 + g2*MDWBT*(-vd^2 + vu^2)) +
(g1*MDBS*(vd^2 - vu^2) - Sqrt[2]*(LamSD*MuD*vd^2 + LamSU*MuU*vu^2))*(16*MDWBT^2 + 4*mT2 + LamTD^2*vd^2 + LamTU^2*vu^2 + 8*B[superT]))/
((LamSD*LamTD*vd^2 + LamSU*LamTU*vu^2)^2 - (8*MDBS^2 + 2*mS2 + LamSD^2*vd^2 + LamSU^2*vu^2 + 4*B[superS])*
(16*MDWBT^2 + 4*mT2 + LamTD^2*vd^2 + LamTU^2*vu^2 + 8*B[superT]))) },

{ vT,(Sqrt[2]*g1*MDBS*(-vd^2 + vu^2) + 2*(LamSD*MuD*vd^2 + LamSU*MuU*vu^2) -
(Sqrt[2]*(8*MDBS^2 + 2*mS2 + LamSD^2*vd^2 + LamSU^2*vu^2 + 4*B[superS])*
(Sqrt[2]*(LamSD*LamTD*vd^2 + LamSU*LamTU*vu^2)*(LamTD*MuD*vd^2 + LamTU*MuU*vu^2 + g2*MDWBT*(-vd^2 + vu^2)) +
(g1*MDBS*(vd^2 - vu^2) - Sqrt[2]*(LamSD*MuD*vd^2 + LamSU*MuU*vu^2))*(16*MDWBT^2 + 4*mT2 + LamTD^2*vd^2 + LamTU^2*vu^2 +
8*B[superT])))/((LamSD*LamTD*vd^2 + LamSU*LamTU*vu^2)^2 - (8*MDBS^2 + 2*mS2 + LamSD^2*vd^2 + LamSU^2*vu^2 + 4*B[superS])*
(16*MDWBT^2 + 4*mT2 + LamTD^2*vd^2 + LamTU^2*vu^2 + 8*B[superT])))/(LamSD*LamTD*vd^2 + LamSU*LamTU*vu^2) }

};

(*Assigns SPheno variables to the model parameters*)
BoundarySUSYScale={

(*Dirac gaugino masses*)
{MDBS, MDBinput},
{MDWBT,MDWinput},
{MDGoc,MDOinput},

(*R soft masses*)
{mRd2, mRd2input},
{mRu2, mRu2input},

(*New mu terms*)
{MuD, MuDinput},
{MuU, MuUinput},

(*B-terms*)
{B[\[Mu]], BmuInput},
{B[superS], BsInput}, (* ! added by me *)
{B[superT], BtInput}, (* ! added by me *)

(*SUSY R-H-adj Yukawas*)
{LamSD, LSDinput},
{LamSU, LSUinput},
{LamTD, LTDinput},
{LamTU, LTUinput},

(*Adjoint soft masses*)
{mS2, mS2input}, (* ! added by me*)
{mT2, mT2input}, (* ! added by me*)
{moc2, moc2input},

(*Sfermion soft masses*)
(*Set by hand*)
{mq2, DIAGONAL mq2input}, (* ! added by me*)
{mq2[3,3], mq2entry33input}, (* ! added by me*)

{ml2, DIAGONAL ml2input}, (* ! added by me*)
{ml2[3,3], ml2entry33input}, (* ! added by me*)

{mu2, DIAGONAL mu2input}, (* ! added by me*)
{mu2[3,3], DIAGONAL mu2entry33input}, (* ! added by me*)

{md2, DIAGONAL md2input}, (* ! added by me*)
{md2[3,3], DIAGONAL md2entry33input}, (* ! added by me*)

{me2, DIAGONAL me2input}, (* ! added by me*)
{me2[3,3], DIAGONAL me2entry33input}, (* ! added by me*)

(*Forbidden by R-symmetry*)
{B[MuD],0},
{B[MuU],0},
{\[Mu], 0}

};

ListDecayParticles = Automatic;
ListDecayParticles3B = Automatic;

====================================================================

Hi,

The problem might come from the loop corrections, see for instance http://www.stauby.de/sarah_userforum/vi ... ?f=4&t=246
for an explanation how these are treated by SARAH.

Maybe, it helps to work with assumptions, see http://stauby.de/sarah_wiki/index.php?t ... ith_SPheno

Cheers
Florian

Thanks for the response, Florian

Indeed the by-hand implementation of tadpole solutions is explained somewhere in the manual.