1 | (* ************************** *)
|
---|
2 | (* ***** Information ***** *)
|
---|
3 | (* ************************** *)
|
---|
4 |
|
---|
5 | M$ModelName = "SimpDM";
|
---|
6 |
|
---|
7 | M$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 |
|
---|
15 | FeynmanGauge = False;
|
---|
16 |
|
---|
17 | (* ***** Change log ***** *)
|
---|
18 |
|
---|
19 | (* ************************** *)
|
---|
20 | (* *** Interaction orders *** *)
|
---|
21 | (* *** (as used by mg5) *** *)
|
---|
22 | (* ************************** *)
|
---|
23 |
|
---|
24 | (*
|
---|
25 | M$InteractionOrderHierarchy = {
|
---|
26 | {NP, 1}
|
---|
27 | };
|
---|
28 | *)
|
---|
29 |
|
---|
30 | (* ************************** *)
|
---|
31 | (* **** Particle classes **** *)
|
---|
32 | (* ************************** *)
|
---|
33 | M$ClassesDescription = {
|
---|
34 |
|
---|
35 | (* ---------- dark matter ---------- *)
|
---|
36 |
|
---|
37 | (* Majorana fermion *)
|
---|
38 | F[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 *)
|
---|
51 | F[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 | },
|
---|
62 | F[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 *)
|
---|
72 | S[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 *)
|
---|
84 | S[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 *)
|
---|
96 | V[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 *)
|
---|
113 | S[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 *)
|
---|
126 | V[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 *)
|
---|
141 | S[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 | },
|
---|
151 | S[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 |
|
---|
168 | V[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 | (* ************************** *)
|
---|
185 | M$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 |
|
---|
342 | MyGA5 = ProjP - ProjM;
|
---|
343 |
|
---|
344 | LBSMIntNoHC = (
|
---|
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 |
|
---|
366 | LBSMIntHC :=
|
---|
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 |
|
---|
374 | LBSM = LBSMIntNoHC + LBSMIntHC + HC[LBSMIntHC];
|
---|