SimpDM: SimpDM.fr

File SimpDM.fr, 9.9 KB (added by Wojciech, 3 years ago)
Line 
1(* ************************** *)
2(* ***** Information ***** *)
3(* ************************** *)
4
5M$ModelName = "SimpDM";
6
7M$Information = {
8 Authors -> {"Wojciech Kotlarski"},
9 Version -> "1.0.0",
10 Date -> "29. 05. 2020",
11 Institutions -> {"TU Dresden"},
12 Emails -> {"wojciech.kotlarski@tu-dresden.de"}
13};
14
15FeynmanGauge = False;
16
17(* ***** Change log ***** *)
18
19(* ************************** *)
20(* *** Interaction orders *** *)
21(* *** (as used by mg5) *** *)
22(* ************************** *)
23
24(*
25M$InteractionOrderHierarchy = {
26 {NP, 1}
27};
28*)
29
30(* ************************** *)
31(* **** Particle classes **** *)
32(* ************************** *)
33M$ClassesDescription = {
34
35(* ---------- dark matter ---------- *)
36
37(* Majorana fermion *)
38F[98] == {
39 ClassName -> Xm,
40 SelfConjugate -> True,
41 Mass -> {mXm, 125.0},
42 Width -> {wXm, 0.00407},
43 ParticleName -> "Xm",
44 PDG -> 5000002,
45 PropagatorLabel -> "MFDM",
46 PropagatorArrow -> None,
47 FullName -> "Majorana fermion DM"
48},
49
50(* Dirac fermion *)
51F[99] == {
52 ClassName -> Xd,
53 SelfConjugate -> False,
54 Mass -> {mXd, 125.0},
55 Width -> {wXd,0.00407},
56 PDG -> 5000003,
57 QuantumNumbers -> {LeptonNumber -> 0, Q -> 0},
58 PropagatorLabel -> "DFDM",
59 PropagatorArrow -> None,
60 FullName -> "Dirac fermion DM"
61},
62F[62] == {
63 ClassName -> XdL,
64 QuantumNumbers -> {LeptonNumber -> 0, Y -> 0},
65 Unphysical -> True,
66 SelfConjugate -> False,
67 QuantumNumbers -> {Y -> 0},
68 Definitions -> { XdL[sp1_] :> Module[{sp2}, ProjM[sp1,sp2] Xd[sp2]] }
69},
70
71(* real scalar *)
72S[98] == {
73 ClassName -> Xr,
74 SelfConjugate -> True,
75 Mass -> {mXr, 125.0},
76 Width -> {wXr, 0.00407},
77 PDG -> 5000004,
78 PropagatorLabel -> "RSDM",
79 PropagatorArrow -> None,
80 FullName -> "real scalar DM"
81},
82
83(* complex scalar *)
84S[99] == {
85 ClassName -> Xc,
86 SelfConjugate -> False,
87 Mass -> {mXc, 125.0},
88 Width -> {wXc, 0.00407},
89 PDG -> 5000005,
90 PropagatorLabel -> "CSDM",
91 PropagatorArrow -> None,
92 FullName -> "complex scalar DM"
93},
94
95(* real vector *)
96V[99] == {
97 ClassName -> Xv,
98 SelfConjugate -> True,
99 Mass -> {mXv, 125.0},
100 Width -> {wXv, 0.00407},
101 ParticleName -> "VDM",
102 PDG -> 5000006,
103 PropagatorLabel -> "VDM",
104 PropagatorArrow -> None,
105 FullName -> "VDM"
106},
107
108(* ---------- mediators ---------- *)
109
110(* S-channel *)
111
112(* real scalar *)
113S[90] == {
114 ClassName -> Yr,
115 SelfConjugate -> True,
116 Mass -> {mYr, 1000.},
117 Width -> {wYr, 1},
118 ParticleName -> "YR",
119 PDG -> 5000007,
120 PropagatorLabel -> "CSDM",
121 PropagatorArrow -> None,
122 FullName -> "real scalar mediator"
123},
124
125(* real vector *)
126V[90] == {
127 ClassName -> Yv,
128 SelfConjugate -> True,
129 Mass -> {mYv, 1000.},
130 Width -> {wYv, 1},
131 ParticleName -> "Yv",
132 PDG -> 5000008,
133 PropagatorLabel -> "CSDM",
134 PropagatorArrow -> None,
135 FullName -> "real vector mediator"
136},
137
138(* T-channel *)
139
140(* complex scalar *)
141S[71] == {
142 ClassName -> Tc2,
143 ClassMembers -> {Tc2e, Tc2mu, Tc2ta},
144 Indices -> {Index[Generation]},
145 FlavorIndex -> Generation,
146 SelfConjugate -> False,
147 QuantumNumbers -> {LeptonNumber -> 1, Y -> -1},
148 Unphysical -> True,
149 Definitions -> {Tc2[ff_] :> Tc[ff]}
150},
151S[91] == {
152 ClassName -> Tc,
153 ClassMembers -> {Tce, Tcmu, Tcta},
154 Indices -> {Index[Generation]},
155 FlavorIndex -> Generation,
156 SelfConjugate -> False,
157 Mass -> {mTc, {mTce, 1000}, {mTcmu, 10000}, {mTcta,10000}},
158 QuantumNumbers -> {LeptonNumber -> 1, Q -> -1},
159 Width -> {{wTce, 1}, {wTcmu, 0}, {wTcta, 0}},
160 PDG -> {5000009, 5000109, 5000209},
161 PropagatorLabel -> "CSDM",
162 PropagatorArrow -> None,
163 FullName -> "real scalar mediator"
164},
165
166(* complex vector *)
167
168V[92] == {
169 ClassName -> Tv,
170 SelfConjugate -> False,
171 Mass -> {mTv, 1000},
172 Width -> {wTv, 1},
173 ParticleName -> "Tv",
174 PDG -> 6000010,
175 QuantumNumbers -> {LeptonNumber -> 1, Q -> -1},
176 PropagatorLabel -> "CSDM",
177 PropagatorArrow -> None,
178 FullName -> "real vector mediator"
179}
180};
181
182(* ************************** *)
183(* ***** Parameters ***** *)
184(* ************************** *)
185M$Parameters = {
186
187 gXdXdYr1 == {
188 ParameterType -> External,
189 BlockName -> DMINPUTS,
190 ComplexParameter -> False,
191 OrderBlock -> 1,
192 Value -> 1,
193 InteractionOrder -> {NP,1},
194 Description -> "Fermi constant"
195 },
196 gXdXdYrG5 == {
197 ParameterType -> External,
198 BlockName -> DMINPUTS,
199 ComplexParameter -> False,
200 OrderBlock -> 2,
201 Value -> 1,
202 InteractionOrder -> {NP,1},
203 Description -> "Fermi constant"
204 },
205 gEEYr1 == {
206 ParameterType -> External,
207 BlockName -> DMINPUTS,
208 ComplexParameter -> False,
209 OrderBlock -> 3,
210 Value -> 1,
211 InteractionOrder -> {NP,1},
212 Description -> "Fermi constant"
213 },
214 gEEYrG5 == {
215 ParameterType -> External,
216 BlockName -> DMINPUTS,
217 ComplexParameter -> False,
218 OrderBlock -> 4,
219 Value -> 1,
220 InteractionOrder -> {NP,1},
221 Description -> "Fermi constant"
222 },
223 gXdXdYv1 == {
224 ParameterType -> External,
225 BlockName -> DMINPUTS,
226 ComplexParameter -> False,
227 OrderBlock -> 5,
228 Value -> 1,
229 InteractionOrder -> {NP,1},
230 Description -> "Fermi constant"
231 },
232 gXdXdYvG5 == {
233 ParameterType -> External,
234 BlockName -> DMINPUTS,
235 ComplexParameter -> False,
236 OrderBlock -> 6,
237 Value -> 1,
238 InteractionOrder -> {NP,1},
239 Description -> "Fermi constant"
240 },
241
242 gEEYv1 == {
243 ParameterType -> External,
244 BlockName -> DMINPUTS,
245 ComplexParameter -> False,
246 OrderBlock -> 7,
247 Value -> 1,
248 InteractionOrder -> {NP,1},
249 Description -> "Fermi constant"
250 },
251 gEEYvG5 == {
252 ParameterType -> External,
253 BlockName -> DMINPUTS,
254 ComplexParameter -> False,
255 OrderBlock -> 8,
256 Value -> 1,
257 InteractionOrder -> {NP,1},
258 Description -> "Fermi constant"
259 },
260 gXrXrYr == {
261 ParameterType -> External,
262 BlockName -> DMINPUTS,
263 ComplexParameter -> False,
264 OrderBlock -> 9,
265 Value -> 1,
266 InteractionOrder -> {NP,1},
267 Description -> "Fermi constant"
268 },
269 gXcXcYr == {
270 ParameterType -> External,
271 BlockName -> DMINPUTS,
272 ComplexParameter -> False,
273 OrderBlock -> 10,
274 Value -> 1,
275 InteractionOrder -> {NP,1},
276 Description -> "complex scalar - complex scalar - real scalar coupling"
277 },
278 gXcXcYv == {
279 ParameterType -> External,
280 BlockName -> DMINPUTS,
281 ComplexParameter -> False,
282 OrderBlock -> 11,
283 Value -> 1,
284 InteractionOrder -> {NP,1},
285 Description -> "complex scalar - complex scalar - real vector coupling"
286 },
287 (* all parameters in UFO for Whizard should be real, otherwise one gets a warning like
288 "UFO warning: invalid complex declaration of input parameter 'gEXdTc' ignored!"
289 see https://answers.launchpad.net/whizard/+faq/3056 *)
290 gEXdTcR == {
291 ParameterType -> External,
292 BlockName -> DMINPUTS,
293 OrderBlock -> 12,
294 Value -> 1,
295 InteractionOrder -> {NP,1},
296 Description -> "complex scalar - complex scalar - real vector coupling"
297 },
298 gEXdTcI == {
299 ParameterType -> External,
300 BlockName -> DMINPUTS,
301 OrderBlock -> 13,
302 Value -> 1,
303 InteractionOrder -> {NP,1},
304 Description -> "complex scalar - complex scalar - real vector coupling"
305 },
306 gEXdTv1 == {
307 ParameterType -> External,
308 BlockName -> DMINPUTS,
309 OrderBlock -> 14,
310 Value -> 1,
311 InteractionOrder -> {NP,1},
312 Description -> "complex scalar - complex scalar - real vector coupling"
313 },
314 gEXdTvG5 == {
315 ParameterType -> External,
316 BlockName -> DMINPUTS,
317 OrderBlock -> 15,
318 Value -> 1,
319 InteractionOrder -> {NP,1},
320 Description -> "complex scalar - complex scalar - real vector coupling"
321 },
322 gXmXmYr1 == {
323 ParameterType -> External,
324 BlockName -> DMINPUTS,
325 OrderBlock -> 16,
326 Value -> 1,
327 InteractionOrder -> {NP,1},
328 Description -> "complex scalar - complex scalar - real vector coupling"
329 },
330 gXmXmYvG5 == {
331 ParameterType -> External,
332 BlockName -> DMINPUTS,
333 OrderBlock -> 17,
334 Value -> 1,
335 InteractionOrder -> {NP,1},
336 Description -> "complex scalar - complex scalar - real vector coupling"
337 }
338};
339
340(* ---------- Lagrangian ---------- *)
341
342MyGA5 = ProjP - ProjM;
343
344LBSMIntNoHC = (
345 ebar.(gEEYr1 + I MyGA5 gEEYrG5).e Yr +
346 ebar.Ga[mu].(gEEYv1 + MyGA5 gEEYvG5).e Yv[mu] +
347
348 mubar.(gEEYr1 + I MyGA5 gEEYrG5).mu Yr +
349 mubar.Ga[mu].(gEEYv1 + MyGA5 gEEYvG5).mu Yv[mu] +
350
351 gXrXrYr Xr Xr Yr +
352 gXcXcYr Xcbar Xc Yr +
353 I gXcXcYv (Xcbar del[Xc,mu]-del[Xcbar,mu] Xc) Yv[mu] +
354 Xdbar.(gXdXdYr1 + I MyGA5 gXdXdYrG5).Xd Yr +
355 Xdbar.Ga[mu].(gXdXdYv1 + MyGA5 gXdXdYvG5).Xd Yv[mu] +
356
357 (* Majorana DM *)
358 gXmXmYr1 Xmbar.Xm Yr +
359 gXmXmYvG5 Xmbar.Ga[mu].MyGA5.Xm Yv[mu] +
360
361 (* if written using Tc (and not Tc2) then there are no Tc - gauge boson interactions *)
362 DC[Tc2bar[ff1], mu] DC[Tc2[ff1], mu] +
363 -1/4 FS[Tvbar, mu, nu] FS[Tv, mu, nu]
364 );
365
366LBSMIntHC :=
367 Module[{sp1, sp2, ff1, mu},
368 ExpandIndices[
369 (gEXdTcR + gEXdTcI I) XdLbar[sp1].lR[sp1, ff1] Tc2bar[ff1]
370 ] +
371 Xdbar.Ga[mu].(gEXdTv1 + MyGA5 gEXdTvG5).e Tvbar[mu]
372 ];
373
374LBSM = LBSMIntNoHC + LBSMIntHC + HC[LBSMIntHC];