Effective Operator

[Geilson]

Dear Staub,

I was trying to write a neutrino mass in a susy version of 331-v2 through effective operator l.rho.l.rho (as l.Hu.l.Hu). However, the operator displays an error:

Calculate superpotential: 1/1(All Done)
Part::pkspec1: The expression a cannot be used as a part specification. >>
Part::pkspec1: The expression a cannot be used as a part specification. >>
checking charge conservation: local symmetries okay. global symmetries okay.

I tried some more operators with 4 triplets and always the same error show up. Doesn't SARAH support this kind of operators?

Cheers,
Geilson

[FStaub]

Hi,

so far the option to define effective superpotential terms A.B.C.D was only tested for case where A.B and C.D are already gauge invariances. Since this is not the case for your terms I guess the problem might originate from that. I'll try to check the next days if this can be easily generalised.

Best,
Florian

[FStaub]

Sorry, forget my last reply. I was thinking in terms of the SM gauge group, ie. combining triplets/doublets under SU(2). So, I tested the following toy model:

Code: Select all

(* Gauge Superfields *)

Gauge[[1]]={B,  U[1],  xcharge,     g1,False, 0, 1};
Gauge[[2]]={WB, SU[3], left,        g2,True,  0, 1};
Gauge[[3]]={G,  SU[3], color,       g3,False, 0, 1};


(* Chiral Superfields *)

SuperFields[[1]] = {l, 3, {RL, -vL, eL},  -1, -3, 1};
SuperFields[[2]] = {e, 3, conj[eR],     1, 1,  1};
SuperFields[[3]] = {R, 3, conj[RR],     2, 1,  1};

SuperFields[[4]] =  {rho,  1, {rhopp, rho0, rhop},  1, 3, 1};
SuperFields[[5]] =  {eta,  1, {eta2p, eta1m, eta0},  0, 3, 1};
SuperFields[[6]] =  {chi,  1, {chi0, chimm, chim}, -1, 3, 1};


(*----------------------------------------------*)
(*   DEFINITION                                 *)
(*----------------------------------------------*)

NameOfStates={GaugeES};

SuperPotential = Wop l.rho.l.rho + Yl l.e.eta + Ye l.R.chi

SARAH runs for me without any error messages. Thus, if you have the problems with the newest version of SARAH, I would need your model files to see why it doesn't work for you.

Cheers,
Florian

[Geilson]

Hi,

The problem was the Mathematica version. I work in version 10 of Mathematica and using the version 8 of my workmate everything run well. Thanks.


Best,
Geilson