1 | (* ************************************************************************* *)
|
---|
2 | (* ***** ***** *)
|
---|
3 | (* ***** FeynRules model file: VLferm ***** *)
|
---|
4 | (* ***** Authors: J. Fuentes-Martin and D. Faroughy ***** *)
|
---|
5 | (* ***** ***** *)
|
---|
6 | (* ************************************************************************* *)
|
---|
7 |
|
---|
8 | (* ************************** *)
|
---|
9 | (* ***** Information ***** *)
|
---|
10 | (* ************************** *)
|
---|
11 |
|
---|
12 | M$ModelName = "VLferm";
|
---|
13 |
|
---|
14 | M$Information = {Authors -> {"J. Fuentes-Martin", "D. Faroughy"},
|
---|
15 | Version -> "1.0",
|
---|
16 | Date -> "21.06.2022",
|
---|
17 | Institutions -> {"University of Granada", "University of Zurich"},
|
---|
18 | Emails -> {"javier.fuentes@ugr.es", "faroughy@physik.uzh.ch"}
|
---|
19 | };
|
---|
20 |
|
---|
21 | (* ************************** *)
|
---|
22 | (* *** Interaction orders *** *)
|
---|
23 | (* ************************** *)
|
---|
24 |
|
---|
25 | M$InteractionOrderHierarchy = { {QCD,1}, {NP,1}, {QED,2} };
|
---|
26 |
|
---|
27 |
|
---|
28 | (* ************************** *)
|
---|
29 | (* **** Particle classes **** *)
|
---|
30 | (* ************************** *)
|
---|
31 |
|
---|
32 | M$ClassesDescription = {
|
---|
33 | (* Vector-like fermions: physical fields *)
|
---|
34 | F[5] == {
|
---|
35 | ClassName -> Ln,
|
---|
36 | SelfConjugate -> False,
|
---|
37 | Mass -> {MLn, 1000},
|
---|
38 | Width -> {WLn, 10},
|
---|
39 | QuantumNumbers -> {LeptonNumber -> 1},
|
---|
40 | PropagatorLabel -> "Ln",
|
---|
41 | PropagatorType -> Straight,
|
---|
42 | PropagatorArrow -> Forward,
|
---|
43 | PDG -> 72,
|
---|
44 | ParticleName -> "Ln",
|
---|
45 | AntiParticleName -> "Ln~",
|
---|
46 | FullName -> "Neutral vector-like lepton"
|
---|
47 | },
|
---|
48 | F[6] == {
|
---|
49 | ClassName -> Le,
|
---|
50 | SelfConjugate -> False,
|
---|
51 | Mass -> {MLe, 1000},
|
---|
52 | Width -> {WLe, 10},
|
---|
53 | QuantumNumbers -> {Q -> -1, LeptonNumber -> 1},
|
---|
54 | PropagatorLabel -> "Le",
|
---|
55 | PropagatorType -> Straight,
|
---|
56 | PropagatorArrow -> Forward,
|
---|
57 | PDG -> 71,
|
---|
58 | ParticleName -> "Le-",
|
---|
59 | AntiParticleName -> "Le+",
|
---|
60 | FullName -> "Charged vector-like lepton"
|
---|
61 | },
|
---|
62 | F[7] == {
|
---|
63 | ClassName -> Qu,
|
---|
64 | SelfConjugate -> False,
|
---|
65 | Indices -> {Index[Colour]},
|
---|
66 | Mass -> {MQu, 2000},
|
---|
67 | Width -> {WQu, 10},
|
---|
68 | QuantumNumbers -> {Q -> 2/3},
|
---|
69 | PropagatorLabel -> "Qu",
|
---|
70 | PropagatorType -> Straight,
|
---|
71 | PropagatorArrow -> Forward,
|
---|
72 | PDG -> 62,
|
---|
73 | ParticleName -> "Qu",
|
---|
74 | AntiParticleName -> "Qu~",
|
---|
75 | FullName -> "Up-type vector-like quark"
|
---|
76 | },
|
---|
77 | F[8] == {
|
---|
78 | ClassName -> Qd,
|
---|
79 | SelfConjugate -> False,
|
---|
80 | Indices -> {Index[Colour]},
|
---|
81 | Mass -> {MQd, 2000},
|
---|
82 | Width -> {WQd, 10},
|
---|
83 | QuantumNumbers -> {Q -> -1/3},
|
---|
84 | PropagatorLabel -> "Qd",
|
---|
85 | PropagatorType -> Straight,
|
---|
86 | PropagatorArrow -> Forward,
|
---|
87 | PDG -> 61,
|
---|
88 | ParticleName -> "Qd",
|
---|
89 | AntiParticleName -> "Qd~",
|
---|
90 | FullName -> "Down-type vector-like quark"
|
---|
91 | },
|
---|
92 |
|
---|
93 | (* Vector-like fermions: unphysical fields *)
|
---|
94 |
|
---|
95 | F[16] == {
|
---|
96 | ClassName -> chiL,
|
---|
97 | Unphysical -> True,
|
---|
98 | Indices -> {Index[SU2D]},
|
---|
99 | FlavorIndex -> SU2D,
|
---|
100 | SelfConjugate -> False,
|
---|
101 | QuantumNumbers -> {Y -> -1/2},
|
---|
102 | Definitions -> {
|
---|
103 | chiL[sp1_,1] :> Ln[sp1],
|
---|
104 | chiL[sp1_,2] :> Le[sp1]
|
---|
105 | }
|
---|
106 | },
|
---|
107 | F[17] == {
|
---|
108 | ClassName -> chiQ,
|
---|
109 | Unphysical -> True,
|
---|
110 | Indices -> {Index[SU2D], Index[Colour]},
|
---|
111 | FlavorIndex -> SU2D,
|
---|
112 | SelfConjugate -> False,
|
---|
113 | QuantumNumbers -> {Y -> 1/6},
|
---|
114 | Definitions -> {
|
---|
115 | chiQ[sp1_,1,cc_] :> Qu[sp1,cc],
|
---|
116 | chiQ[sp1_,2,cc_] :> Qd[sp1,cc]
|
---|
117 | }
|
---|
118 | }
|
---|
119 | };
|
---|
120 |
|
---|
121 | (* ************************** *)
|
---|
122 | (* ***** Parameters ***** *)
|
---|
123 | (* ************************** *)
|
---|
124 |
|
---|
125 | M$Parameters = {
|
---|
126 | (* External Parameters *)
|
---|
127 | betaLQL == {
|
---|
128 | ParameterType -> External,
|
---|
129 | Value -> 1.0,
|
---|
130 | InteractionOrder -> {NP,0},
|
---|
131 | TeX -> Subsuperscript[\[Beta],L,QL],
|
---|
132 | BlockName -> NPLQCOUP,
|
---|
133 | OrderBlock -> 8,
|
---|
134 | Description -> "LQ Q L left-handed coupling"
|
---|
135 | },
|
---|
136 | betaRQL == {
|
---|
137 | ParameterType -> External,
|
---|
138 | Value -> 1.0,
|
---|
139 | InteractionOrder -> {NP,0},
|
---|
140 | TeX -> Subsuperscript[\[Beta],R,QL],
|
---|
141 | BlockName -> NPLQCOUP,
|
---|
142 | OrderBlock -> 9,
|
---|
143 | Description -> "LQ Q L right-handed coupling"
|
---|
144 | },
|
---|
145 | betaL3L == {
|
---|
146 | ParameterType -> External,
|
---|
147 | Value -> 1.0,
|
---|
148 | InteractionOrder -> {NP,0},
|
---|
149 | TeX -> Subsuperscript[\[Beta],L,3L],
|
---|
150 | BlockName -> NPLQCOUP,
|
---|
151 | OrderBlock -> 10,
|
---|
152 | Description -> "LQ q3 L coupling"
|
---|
153 | },
|
---|
154 | betaLQ3 == {
|
---|
155 | ParameterType -> External,
|
---|
156 | Value -> 1.0,
|
---|
157 | InteractionOrder -> {NP,0},
|
---|
158 | TeX -> Subsuperscript[\[Beta],L,Q3],
|
---|
159 | BlockName -> NPLQCOUP,
|
---|
160 | OrderBlock -> 11,
|
---|
161 | Description -> "LQ Q l3 coupling"
|
---|
162 | },
|
---|
163 | betaL2L == {
|
---|
164 | ParameterType -> External,
|
---|
165 | Value -> 0.2,
|
---|
166 | InteractionOrder -> {NP,0},
|
---|
167 | TeX -> Subsuperscript[\[Beta],L,2L],
|
---|
168 | BlockName -> NPLQCOUP,
|
---|
169 | OrderBlock -> 12,
|
---|
170 | Description -> "LQ q2 L coupling"
|
---|
171 | },
|
---|
172 | betaLQ2 == {
|
---|
173 | ParameterType -> External,
|
---|
174 | Value -> 0.2,
|
---|
175 | InteractionOrder -> {NP,0},
|
---|
176 | TeX -> Subsuperscript[\[Beta],L,Q2],
|
---|
177 | BlockName -> NPLQCOUP,
|
---|
178 | OrderBlock -> 13,
|
---|
179 | Description -> "LQ Q l2 coupling"
|
---|
180 | },
|
---|
181 |
|
---|
182 | zetaLQQ == {
|
---|
183 | ParameterType -> External,
|
---|
184 | Value -> 1.0,
|
---|
185 | InteractionOrder -> {NP,0},
|
---|
186 | TeX -> Subsuperscript[\[Zeta],L,QQ],
|
---|
187 | BlockName -> NPZPCOUP,
|
---|
188 | OrderBlock -> 13,
|
---|
189 | Description -> "Zp Q Q left-handed coupling"
|
---|
190 | },
|
---|
191 | zetaRQQ == {
|
---|
192 | ParameterType -> External,
|
---|
193 | Value -> 1.0,
|
---|
194 | InteractionOrder -> {NP,0},
|
---|
195 | TeX -> Subsuperscript[\[Zeta],R,QQ],
|
---|
196 | BlockName -> NPZPCOUP,
|
---|
197 | OrderBlock -> 14,
|
---|
198 | Description -> "Zp Q Q right-handed coupling"
|
---|
199 | },
|
---|
200 | zetaLQ2 == {
|
---|
201 | ParameterType -> External,
|
---|
202 | Value -> 0.2,
|
---|
203 | InteractionOrder -> {NP,0},
|
---|
204 | TeX -> Subsuperscript[\[Zeta],L,Q2],
|
---|
205 | BlockName -> NPZPCOUP,
|
---|
206 | OrderBlock -> 15,
|
---|
207 | Description -> "Zp Q q2 left-handed coupling"
|
---|
208 | },
|
---|
209 | zetaLLL == {
|
---|
210 | ParameterType -> External,
|
---|
211 | Value -> 1.0,
|
---|
212 | InteractionOrder -> {NP,0},
|
---|
213 | TeX -> Subsuperscript[\[Zeta],L,LL],
|
---|
214 | BlockName -> NPZPCOUP,
|
---|
215 | OrderBlock -> 16,
|
---|
216 | Description -> "Zp L L left-handed coupling"
|
---|
217 | },
|
---|
218 | zetaRLL == {
|
---|
219 | ParameterType -> External,
|
---|
220 | Value -> 1.0,
|
---|
221 | InteractionOrder -> {NP,0},
|
---|
222 | TeX -> Subsuperscript[\[Zeta],R,LL],
|
---|
223 | BlockName -> NPZPCOUP,
|
---|
224 | OrderBlock -> 17,
|
---|
225 | Description -> "Zp L L right-handed coupling"
|
---|
226 | },
|
---|
227 | zetaLL2 == {
|
---|
228 | ParameterType -> External,
|
---|
229 | Value -> 0.2,
|
---|
230 | InteractionOrder -> {NP,0},
|
---|
231 | TeX -> Subsuperscript[\[Zeta],L,L2],
|
---|
232 | BlockName -> NPZPCOUP,
|
---|
233 | OrderBlock -> 18,
|
---|
234 | Description -> "Zp L l2 left-handed coupling"
|
---|
235 | },
|
---|
236 |
|
---|
237 | kappaLQQ == {
|
---|
238 | ParameterType -> External,
|
---|
239 | Value -> 1.0,
|
---|
240 | InteractionOrder -> {NP,0},
|
---|
241 | TeX -> Subsuperscript[\[Kappa],L,QQ],
|
---|
242 | BlockName -> NPGPCOUP,
|
---|
243 | OrderBlock -> 10,
|
---|
244 | Description -> "Gp Q Q left-handed coupling"
|
---|
245 | },
|
---|
246 | kappaRQQ == {
|
---|
247 | ParameterType -> External,
|
---|
248 | Value -> 1.0,
|
---|
249 | InteractionOrder -> {NP,0},
|
---|
250 | TeX -> Subsuperscript[\[Kappa],R,QQ],
|
---|
251 | BlockName -> NPGPCOUP,
|
---|
252 | OrderBlock -> 11,
|
---|
253 | Description -> "Gp Q Q left-handed coupling"
|
---|
254 | },
|
---|
255 | kappaLQ2 == {
|
---|
256 | ParameterType -> External,
|
---|
257 | Value -> 0.2,
|
---|
258 | InteractionOrder -> {NP,0},
|
---|
259 | TeX -> Subsuperscript[\[Kappa],L,Q2],
|
---|
260 | BlockName -> NPGPCOUP,
|
---|
261 | OrderBlock -> 12,
|
---|
262 | Description -> "Gp Q q2 left-handed coupling"
|
---|
263 | },
|
---|
264 |
|
---|
265 | (* Internal Parameters *)
|
---|
266 | betaLchiQ == {
|
---|
267 | ParameterType -> Internal,
|
---|
268 | Indices -> {Index[Generation]},
|
---|
269 | Definitions -> {betaLchiQ[1] -> 0},
|
---|
270 | Value -> {betaLchiQ[2] -> betaLQ2, betaLchiQ[3] -> betaLQ3},
|
---|
271 | InteractionOrder -> {NP,0},
|
---|
272 | TeX -> Subsuperscript[\[Beta],L,Ql],
|
---|
273 | Description -> "Q-SM-LQ coupling vector"
|
---|
274 | },
|
---|
275 | betaLchiL == {
|
---|
276 | ParameterType -> Internal,
|
---|
277 | Indices -> {Index[Generation]},
|
---|
278 | Definitions -> {betaLchiL[1] -> 0},
|
---|
279 | Value -> {betaLchiL[2] -> betaL2L, betaLchiL[3] -> betaL3L},
|
---|
280 | InteractionOrder -> {NP,0},
|
---|
281 | TeX -> Subsuperscript[\[Beta],L,qL],
|
---|
282 | Description -> "SM-L-Zp coupling vector"
|
---|
283 | },
|
---|
284 | zetaLchiQ == {
|
---|
285 | ParameterType -> Internal,
|
---|
286 | Indices -> {Index[Generation]},
|
---|
287 | Definitions -> {zetaLchiQ[1] -> 0, zetaLchiQ[3] -> 0},
|
---|
288 | Value -> {zetaLchiQ[2] -> zetaLQ2},
|
---|
289 | InteractionOrder -> {NP,0},
|
---|
290 | TeX -> Subsuperscript[\[Zeta],L,Qq],
|
---|
291 | Description -> "SM-Q-Zp coupling vector"
|
---|
292 | },
|
---|
293 | zetaLchiL == {
|
---|
294 | ParameterType -> Internal,
|
---|
295 | Indices -> {Index[Generation]},
|
---|
296 | Definitions -> {zetaLchiL[1] -> 0, zetaLchiL[3] -> 0},
|
---|
297 | Value -> {zetaLchiL[2] -> zetaLL2},
|
---|
298 | InteractionOrder -> {NP,0},
|
---|
299 | TeX -> Subsuperscript[\[Zeta],L,Ll],
|
---|
300 | Description -> "SM-L-Zp coupling vector"
|
---|
301 | },
|
---|
302 | kappaLchiQ == {
|
---|
303 | ParameterType -> Internal,
|
---|
304 | Indices -> {Index[Generation]},
|
---|
305 | Definitions -> {kappaLchiQ[1] -> 0, kappaLchiQ[3] -> 0},
|
---|
306 | Value -> {kappaLchiQ[2] -> kappaLQ2},
|
---|
307 | InteractionOrder -> {NP,0},
|
---|
308 | TeX -> Subsuperscript[\[Kappa],L,Qq],
|
---|
309 | Description -> "SM-Q-Gp coupling vector"
|
---|
310 | }
|
---|
311 | };
|
---|
312 |
|
---|
313 |
|
---|
314 | (* ************************** *)
|
---|
315 | (* ***** Lagrangian ***** *)
|
---|
316 | (* ************************** *)
|
---|
317 |
|
---|
318 | LVLfermkin := Block[{mu},
|
---|
319 | ExpandIndices[
|
---|
320 | I*(chiQbar.Ga[mu].DC[chiQ, mu] + chiLbar.Ga[mu].DC[chiL, mu]),
|
---|
321 | FlavorExpand->{SU2W,SU2D}]
|
---|
322 | ];
|
---|
323 |
|
---|
324 | LU1VL := Block[{ff1,ff2,s1,s2,s3,cc1,ii1,mu},
|
---|
325 | ExpandIndices[
|
---|
326 | gU/Sqrt[2] * VLQ[mu,cc1] *
|
---|
327 | (
|
---|
328 | betaLQL * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiL[s3,ii1]
|
---|
329 | + betaRQL * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiL[s3,ii1]
|
---|
330 | + betaLchiL[ff1] * newCKM[ff2,ff1] * uqbar[s1,ff2,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * Ln[s3]
|
---|
331 | + betaLchiL[ff1] * dqbar[s1,ff1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * Le[s3]
|
---|
332 | + betaLchiQ[ff1] * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * LL[s3,ii1,ff1]
|
---|
333 | )
|
---|
334 | ,FlavorExpand-> SU2D]
|
---|
335 | ];
|
---|
336 |
|
---|
337 | LZpVL := Block[{s1,s2,s3,mu,cc1,ii1,ff1},
|
---|
338 | ExpandIndices[
|
---|
339 | gZp/(2 * Sqrt[6]) * Zp[mu] *
|
---|
340 | (
|
---|
341 | zetaLQQ * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiQ[s3,ii1,cc1]
|
---|
342 | + zetaRQQ * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiQ[s3,ii1,cc1]
|
---|
343 | -3 * zetaLLL * chiLbar[s1,ii1] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiL[s3,ii1]
|
---|
344 | -3 * zetaRLL * chiLbar[s1,ii1] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiL[s3,ii1]
|
---|
345 | )
|
---|
346 | ,FlavorExpand-> SU2D]
|
---|
347 | ];
|
---|
348 |
|
---|
349 | LZpVLMix := Block[{s1,s2,s3,mu,cc1,ii1,ff1,ff2},
|
---|
350 | ExpandIndices[
|
---|
351 | gZp/(2 * Sqrt[6]) * Zp[mu] *
|
---|
352 | (
|
---|
353 | zetaLchiQ[ff1] * HC[newCKM[ff2,ff1]] * Qubar[s1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * uq[s3,ff2,cc1]
|
---|
354 | + zetaLchiQ[ff1] * Qdbar[s1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * dq[s3,ff1,cc1]
|
---|
355 | -3 * zetaLchiL[ff1] * chiLbar[s1,ii1] * Ga[mu,s1,s2] * ProjM[s2,s3] * LL[s3,ii1,ff1]
|
---|
356 | )
|
---|
357 | ,FlavorExpand-> SU2D]
|
---|
358 | ];
|
---|
359 |
|
---|
360 | LGpVL := Block[{s1,s2,s3,aa1,cc1,cc2,mu,ii1,ff1},
|
---|
361 | ExpandIndices[
|
---|
362 | gGp * Gp[mu,aa1] *
|
---|
363 | (
|
---|
364 | kappaLQQ * chiQbar[s1,ii1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiQ[s3,ii1,cc2]
|
---|
365 | + kappaRQQ * chiQbar[s1,ii1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiQ[s3,ii1,cc2]
|
---|
366 | )
|
---|
367 | ,FlavorExpand-> SU2D]
|
---|
368 | ];
|
---|
369 |
|
---|
370 |
|
---|
371 | LGpVLMix := Block[{s1,s2,s3,aa1,cc1,cc2,mu,ff1,ff2},
|
---|
372 | ExpandIndices[
|
---|
373 | gGp * Gp[mu,aa1] *
|
---|
374 | (
|
---|
375 | kappaLchiQ[ff1] * HC[newCKM[ff2,ff1]] * Qubar[s1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * uq[s3,ff2,cc2]
|
---|
376 | + kappaLchiQ[ff1] * Qdbar[s1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * dq[s3,ff1,cc2]
|
---|
377 | )
|
---|
378 | ,FlavorExpand-> SU2D]
|
---|
379 | ];
|
---|
380 |
|
---|
381 | LVLferm := LVLfermkin + LU1VL + HC[LU1VL] + LZpVL + LZpVLMix + HC[LZpVLMix] + LGpVL + LGpVLMix + HC[LGpVLMix];
|
---|