Dear Florian, All,

I am a new user of SARAH and I am currently trying to generate an MSSM model with high scale type 1 seesaw. I started from the Seesaw1 model in the SARAH package, with minor changes (I modified the left sneutrinos sector, so that now they are written in the same manner as in the MSSM).

However, I get some error when performing this step (see ERROR_SARAH2.png in attachment) : Writing SPheno two loop diagrammatic Higgs mass
The detailed error is then given in SARAH_ERROR.png (in attachment)

It seems that this error gets translated into Fortran code for SPheno in line 6105 of 2LPole_Seesaw1_MSSM_SU5.txt (in attachment) the error in the fortran code reads "coup3 = cplSvSvcSvcSv({{}, {Sv, i4}, {conj[Sv], i2}, {}}[[1,2]],i4,i2,{{}, {Sv, i4}, {conj[Sv], i2}, {}}[[4,2]])"

When sending this coupling to 0, SPheno compile otherwise it does not understand what is this line. It looks like the Mathematica code is pasted directely intro Fortran, likely due to the empty {} arrays (see error above).

I provided in attachment all the SARAH model related files.

Many thanks for your help,
Best,
Jordan

Dear Jordan,

The problem is that you have removed an important part of the Seesaw1 model. By removing the line

{SvL, {0,0}, {sigmaL, \[ImaginaryI]/Sqrt[2]},{phiL,1/Sqrt[2]}}

from the Seesaw1 model definition and the references to sigmaL, phiL, SARAH no-longer treats the sneutrinos correctly. This is because after EWSB the sneutrinos are split into pairs of real fields with different masses, because there is no longer any symmetry keeping them complex. So in the two loop routines SARAH doesn't know what to do with the quartic sneutrino coupling: they are particularly pernickety, I am afraid.

Since you have added a right-handed sneutrino to the model, if you want to be completely correct you ought to do the same for the right-handed sneutrinos and then get the components to mix with the left-handed ones.

Best,
Mark