ZeeBabu: zeebabu_NLO_public.fr

File zeebabu_NLO_public.fr, 10.9 KB (added by Richard Ruiz, 2 years ago)

FeynRules model file for Zee-Babu Lagrangian

Line 
1(* ************************************************************************* *)
2(* FeynRules Model file for a Zee-Babu extension of the Standard Model. *)
3(* Contains two new scalar fields: k--(Y=-2) and h-(Y=-1) that couple *)
4(* to all SM leptons through Yukawa couplings, Z/A bosons through gauge *)
5(* couplings, and the SM Higgs through dimensionless couplings. *)
6(* Neutrinos are massless and aligned with flavor states in this model file. *)
7(* The Lagrangian is written with Goldstone bosons. *)
8(* *)
9(* Contact author: R. Ruiz [rruiz [at] ifj.edu.pl] *)
10(* *)
11(* The original Zee-Babu papers are: *)
12(* 1. Zee, PLB ('85) [Zee:1985rj] *)
13(* 2. Zee, NPB ('86) [Zee:1985id] *)
14(* 3. Babu, PLB ('88) [Babu:1988ki] *)
15(* At LO and NLO in QCD, this model was implemented in: *)
16(* 1. R. Ruiz [arXiv:2206.14833], *)
17(* and extends the default FeynRules SM model file. *)
18(* Please cite accordingly. *)
19(* *)
20(* ************************************************************************* *)
21
22(* ************************** *)
23(* ***** Information ***** *)
24(* ************************** *)
25M$ModelName = "SM_ZeeBabu";
26
27M$Information = {
28 Authors -> {"R. Ruiz"},
29 Version -> "1.0",
30 Date -> "2022 June 30",
31 Institutions -> {"IFJ PAN"},
32 Emails -> {"rruiz@ifj.edu.pl"},
33 References -> {"R. Ruiz [arXiv:2206.14833];",
34 "Zee, PLB ('85) [Zee:1985rj];",
35 "Zee, NPB ('86) [Zee:1985id];",
36 "Babu, PLB ('88) [Babu:1988ki]"
37 },
38 URLs -> {"https://feynrules.irmp.ucl.ac.be/wiki/ZeeBabu"}
39};
40FeynmanGauge = True;
41
42
43(* ************************** *)
44(* ***** Change log ***** *)
45(* ************************** *)
46(* v1.1: Public release of NLO model file *)
47
48(* ************************** *)
49(* ***** Parameters ***** *)
50(* ************************** *)
51M$Parameters = {
52 (* External Parameters *)
53
54 (* ***** Scalar Sector ***** *)
55
56 muZB == {
57 ParameterType -> External,
58 BlockName -> ZeeBabuPotential,
59 OrderBlock -> 1,
60 Value -> 1000.0,
61 ComplexParameter -> False,
62 InteractionOrder -> {NP,1},
63 TeX -> Subscript[mu,ZB],
64 Description -> "LNV mass scale; h-h-k ZB coupling"
65 },
66
67 lamhZB == {
68 ParameterType -> External,
69 BlockName -> ZeeBabuPotential,
70 OrderBlock -> 2,
71 Value -> 1.0,
72 ComplexParameter -> False,
73 InteractionOrder -> {NP,1},
74 TeX -> Subscript[lambda,hZB],
75 Description -> "h-h-h-h ZB coupling"
76 },
77
78 lamkZB == {
79 ParameterType -> External,
80 BlockName -> ZeeBabuPotential,
81 OrderBlock -> 3,
82 Value -> 1.000000001,
83 ComplexParameter -> False,
84 InteractionOrder -> {NP,1},
85 TeX -> Subscript[lambda,kZB],
86 Description -> "k-k-k-k ZB coupling"
87 },
88
89 lamhZBkZB == {
90 ParameterType -> External,
91 BlockName -> ZeeBabuPotential,
92 OrderBlock -> 4,
93 Value -> 1.000000002,
94 ComplexParameter -> False,
95 InteractionOrder -> {NP,1},
96 TeX -> Subscript[lambda,hZBkZB],
97 Description -> "h-h-k-k ZB coupling"
98 },
99
100 lamhZBH == {
101 ParameterType -> External,
102 BlockName -> ZeeBabuPotential,
103 OrderBlock -> 5,
104 Value -> 1.000000003,
105 ComplexParameter -> False,
106 InteractionOrder -> {{NP,1},{QED,1}},
107 TeX -> Subscript[lambda,hZBH],
108 Description -> "h-h-H(-H) ZB-SM coupling"
109 },
110
111 lamkZBH == {
112 ParameterType -> External,
113 BlockName -> ZeeBabuPotential,
114 OrderBlock -> 6,
115 Value -> 1.000000004,
116 ComplexParameter -> False,
117 InteractionOrder -> {{NP,1},{QED,1}},
118 TeX -> Subscript[lambda,kZBH],
119 Description -> "k-k-H(-H) ZB-SM coupling"
120 },
121
122 (* ***** Symmetric Couplings ***** *)
123
124 gee == {
125 ParameterType -> External,
126 BlockName -> ZeeBabuSymmetric,
127 OrderBlock -> 1,
128 Value -> 1.000000010,
129 ComplexParameter -> False,
130 InteractionOrder -> {NP,1},
131 TeX -> Subscript[g,ee],
132 Description -> "ee coupling to k--"
133 },
134
135 gemu == {
136 ParameterType -> External,
137 BlockName -> ZeeBabuSymmetric,
138 OrderBlock -> 2,
139 Value -> 1.000000011,
140 ComplexParameter -> False,
141 InteractionOrder -> {NP,1},
142 TeX -> Subscript[g,emu],
143 Description -> "emu coupling to k--"
144 },
145
146 getau == {
147 ParameterType -> External,
148 BlockName -> ZeeBabuSymmetric,
149 OrderBlock -> 3,
150 Value -> 1.000000012,
151 ComplexParameter -> False,
152 InteractionOrder -> {NP,1},
153 TeX -> Subscript[g,etau],
154 Description -> "etau coupling to k--"
155 },
156
157 gmumu == {
158 ParameterType -> External,
159 BlockName -> ZeeBabuSymmetric,
160 OrderBlock -> 4,
161 Value -> 1.000000013,
162 ComplexParameter -> False,
163 InteractionOrder -> {NP,1},
164 TeX -> Subscript[g,mumu],
165 Description -> "mumu coupling to k--"
166 },
167
168 gmutau == {
169 ParameterType -> External,
170 BlockName -> ZeeBabuSymmetric,
171 OrderBlock -> 5,
172 Value -> 1.000000014,
173 ComplexParameter -> False,
174 InteractionOrder -> {NP,1},
175 TeX -> Subscript[g,mutau],
176 Description -> "mutau coupling to k--"
177 },
178
179 gtautau == {
180 ParameterType -> External,
181 BlockName -> ZeeBabuSymmetric,
182 OrderBlock -> 6,
183 Value -> 1.000000015,
184 ComplexParameter -> False,
185 InteractionOrder -> {NP,1},
186 TeX -> Subscript[g,tautau],
187 Description -> "tautau coupling to k--"
188 },
189
190 (* ***** Anti-symmetric Couplings ***** *)
191
192 femu == {
193 ParameterType -> External,
194 BlockName -> ZeeBabuAntisymmetric,
195 OrderBlock -> 1,
196 Value -> 1.000000020,
197 ComplexParameter -> False,
198 InteractionOrder -> {NP,1},
199 TeX -> Subscript[f,emu],
200 Description -> "emu coupling to h-"
201 },
202
203 fetau == {
204 ParameterType -> External,
205 BlockName -> ZeeBabuAntisymmetric,
206 OrderBlock -> 2,
207 Value -> 1.000000021,
208 ComplexParameter -> False,
209 InteractionOrder -> {NP,1},
210 TeX -> Subscript[f,etau],
211 Description -> "etau coupling to h-"
212 },
213
214 fmutau == {
215 ParameterType -> External,
216 BlockName -> ZeeBabuAntisymmetric,
217 OrderBlock -> 3,
218 Value -> 1.000000022,
219 ComplexParameter -> False,
220 InteractionOrder -> {NP,1},
221 TeX -> Subscript[f,mutau],
222 Description -> "mutau coupling to h-"
223 },
224
225
226 (* Internal Parameters *)
227 mhZBTree == {
228 ParameterType -> Internal,
229 Value -> Sqrt[mhZB^2 - lamhZBH*vev^2 /2],
230 ComplexParameter -> False,
231 Description -> "tree-level mass of h- before EWSB"
232 },
233
234 mkZBTree == {
235 ParameterType -> Internal,
236 Value -> Sqrt[mkZB^2 - lamhZBH*vev^2 /2],
237 ComplexParameter -> False,
238 Description -> "tree-level mass of k-- before EWSB"
239 },
240
241
242(* ************************** *)
243(* ***** Containers ***** *)
244(* ************************** *)
245
246 gll == {
247 ParameterType ->Internal,
248 Indices ->{Index[Generation],Index[Generation]},
249 ComplexParameter->False,
250 Unitary ->False,
251 InteractionOrder->{NP,1},
252 Value ->{
253 gll[1,1]->gee, gll[1,2]->gemu, gll[1,3]->getau,
254 gll[2,1]->gemu, gll[2,2]->gmumu, gll[2,3]->gmutau,
255 gll[3,1]->getau, gll[3,2]->gmutau, gll[3,3]->gtautau},
256 TeX ->Subscript[g,ll],
257 Description ->"symmetric ll coupling matrix to k--"
258 },
259
260 fll == {
261 ParameterType ->Internal,
262 Indices ->{Index[Generation],Index[Generation]},
263 ComplexParameter->False,
264 Unitary ->False,
265 InteractionOrder->{NP,1},
266 Value ->{
267 fll[1,1]->0, fll[1,2]->femu, fll[1,3]->fetau,
268 fll[2,1]->-femu, fll[2,2]->0, fll[2,3]->fmutau,
269 fll[3,1]->-fetau, fll[3,2]->-fmutau, fll[3,3]->0},
270 TeX ->Subscript[f,ll],
271 Description ->"anti-symmetric ll coupling matrix to h-"
272 }
273
274};
275
276
277M$InteractionOrderHierarchy = {
278 {QCD,1},
279 {QED,2},
280 {NP,1}
281};
282
283
284(* ************************** *)
285(* **** Particle classes **** *)
286(* ************************** *)
287M$ClassesDescription = {
288 (*Zee-Babu Scalars*)
289 S[138] == {
290 ClassName -> hZB,
291 SelfConjugate -> False,
292 Mass -> {mhZB,300.},
293 Width -> {whZB,1},
294 QuantumNumbers -> {Q->-1,LeptonNumber -> 2},
295 PropagatorLabel -> "h-",
296 PropagatorType -> ScalarDash,
297 PropagatorArrow -> True,
298 ParticleName -> "h-",
299 AntiParticleName-> "h+",
300 PDG -> 38,
301 FullName -> "h-"
302 },
303
304 S[161] == {
305 ClassName -> kZB,
306 SelfConjugate -> False,
307 Mass -> {mkZB,500.},
308 Width -> {wkZB,1},
309 QuantumNumbers -> {Q->-2,LeptonNumber -> 2},
310 PropagatorLabel -> "k--",
311 PropagatorType -> ScalarDash,
312 PropagatorArrow -> True,
313 ParticleName -> "k--",
314 AntiParticleName-> "k++",
315 PDG -> 61,
316 FullName -> "k--"
317 },
318
319 S[238] == {
320 ClassName -> hGauge,
321 Unphysical -> True,
322 SelfConjugate -> False,
323 QuantumNumbers-> {Y->-1,LeptonNumber -> 2},
324 Definitions -> { hGauge -> hZB }
325 },
326
327 S[261] == {
328 ClassName -> kGauge,
329 Unphysical -> True,
330 SelfConjugate -> False,
331 QuantumNumbers-> {Y->-2,LeptonNumber -> 2},
332 Definitions -> { kGauge -> kZB }
333 }
334
335};
336
337(* ************************** *)
338(* ***** Lagrangian ***** *)
339(* ************************** *)
340
341(* hZB and kZB covariant derivatives and masses before EWSB *)
342LZBKin := DC[hGaugebar,mu] DC[hGauge,mu] + DC[kGaugebar,mu] DC[kGauge,mu] - mhZBTree^2 hGaugebar hGauge - mkZBTree^2 kGaugebar kGauge;
343
344(* hZB and kZB self-couplings *)
345LZBScalarSelf := lamhZB * (hGaugebar hGauge)^2 + lamkZB * (kGaugebar kGauge)^2 + lamhZBkZB * (hGaugebar hGauge) (kGaugebar kGauge);
346
347(* hZB and kZB couplings to SM Higgs *)
348LZBHiggs := lamhZBH (hGaugebar hGauge) * Phibar[ii] * Phi[ii] + lamkZBH (kGaugebar kGauge) * Phibar[ii] * Phi[ii];
349
350(* hZB and kZB LNV coupling *)
351LlnvBare := muZB hGauge hGauge kGaugebar;
352Llnv := LlnvBare + HC[LlnvBare];
353
354(* kZB couplings to leptons *)
355LkllBare := kGaugebar * CC[lRbar][sp1,ff1].gll[ff1,ff2].lR[sp1,ff2];
356Lkll := LkllBare + HC[LkllBare];
357
358(* hZB couplings to leptons *)
359(* ii,jj = isospin; ff1,ff2 = flavor; sp = spin *)
360LhllBare := hGaugebar * CC[LLbar][sp, ii, ff1].fll[ff1,ff2].Eps[ii, jj].LL[sp, jj, ff2];
361Lhll := LhllBare + HC[LhllBare];
362
363(* Full ZB Lagrangian *)
364LScalar := LZBScalarSelf + LZBHiggs + Llnv;
365LZB:= LZBKin + Lkll + Lhll - LScalar;
366
367(* Full SM+ZB Lagrangian *)
368LFull := LSM + LZB;