1 | (* ************************************************************************************ *)
|
---|
2 | (* FeynRules Model file for a heavy neutrino extension of the Standard Model. *)
|
---|
3 | (* Contains three massive Dirac neutrinos that couples to all SM leptons through *)
|
---|
4 | (* mixing between the active (flavor) eigenstates and mass eigenstates. *)
|
---|
5 | (* The Lagrangian is written in the Feynman Gauge. *)
|
---|
6 | (* *)
|
---|
7 | (* Contact author: R. Ruiz [richard.ruiz [at] uclouvain.be] *)
|
---|
8 | (* *)
|
---|
9 | (* The Lagrangian is based on Atre, et. al. [arXiv:0901.3589]. *)
|
---|
10 | (* The model was first implemented in Pascoli, et. al. [arXiv:1805.09335]. *)
|
---|
11 | (* Model files are extensions of default FeynRules SM model file *)
|
---|
12 | (* and the HeavyN libraries by Degrande, et al. [arXiv:1602.06957] *)
|
---|
13 | (* Please cite accordingly. *)
|
---|
14 | (* *)
|
---|
15 | (* Neutrino mass (mN) and mixing parameters (VlN) between heavy mass eigenstate and *)
|
---|
16 | (* flavor eigenstates are taken to be independent, phenomenological parameters. *)
|
---|
17 | (* Mixing parameters (VeN, VmuN, VtaN) are taken to be real. *)
|
---|
18 | (* This allows for maximum flexibility and model independence when calculating rates. *)
|
---|
19 | (* *)
|
---|
20 | (* Feynman gauge is set to true. *)
|
---|
21 | (* ************************************************************************************ *)
|
---|
22 |
|
---|
23 | (* ************************** *)
|
---|
24 | (* ***** Information ***** *)
|
---|
25 | (* ************************** *)
|
---|
26 | M$ModelName = "SM_HeavyN_Dirac";
|
---|
27 |
|
---|
28 | M$Information = {
|
---|
29 | Authors -> {"R. Ruiz"},
|
---|
30 | Version -> "3.1",
|
---|
31 | Date -> "14 November 2018",
|
---|
32 | Institutions -> {"CP3 / Universite Catholique de Louvain"},
|
---|
33 | Emails -> {"richard.ruiz@uclouvain.be"},
|
---|
34 | References -> {"Atre, et al, JHEP0905, 030 (2009) [arXiv:0901.3589]",
|
---|
35 | "Degrande, et al, PRD94, 053002 (2016) [arXiv:1602.06957]",
|
---|
36 | "Pascoli, et al, PLB786, 106 (2018) [arXiv:1805.09335]"},
|
---|
37 | URLs -> {"https://feynrules.irmp.ucl.ac.be/wiki/HeavyN"}
|
---|
38 | };
|
---|
39 | FeynmanGauge = True;
|
---|
40 |
|
---|
41 |
|
---|
42 | (* ************************** *)
|
---|
43 | (* ***** Change log ***** *)
|
---|
44 | (* ************************** *)
|
---|
45 | (* v1.1: Public release of LO model file *)
|
---|
46 | (* v2.1: Added Goldstone couplings for Feynman Gauge and NLO implementation *)
|
---|
47 | (* v2.2: Corrected relative sign between Yukawa and gauge couplings. *)
|
---|
48 | (* v3.1: Public release of Dirac NLO model file. *)
|
---|
49 |
|
---|
50 | (* ************************** *)
|
---|
51 | (* ***** Parameters ***** *)
|
---|
52 | (* ************************** *)
|
---|
53 | M$Parameters = {
|
---|
54 | (* External Parameters *)
|
---|
55 |
|
---|
56 | VeN1 == {
|
---|
57 | ParameterType -> External,
|
---|
58 | BlockName -> NUMIXING,
|
---|
59 | OrderBlock -> 1,
|
---|
60 | Value -> 1.0,
|
---|
61 | ComplexParameter -> False,
|
---|
62 | TeX -> Subscript[V,eN1],
|
---|
63 | Description -> "Mixing between ve flavor/gauge state and N1 mass state"
|
---|
64 | },
|
---|
65 |
|
---|
66 | VeN2 == {
|
---|
67 | ParameterType -> External,
|
---|
68 | BlockName -> NUMIXING,
|
---|
69 | OrderBlock -> 2,
|
---|
70 | Value -> 0.0,
|
---|
71 | ComplexParameter -> False,
|
---|
72 | TeX -> Subscript[V,eN2],
|
---|
73 | Description -> "Mixing between ve flavor/gauge state and N2 mass state"
|
---|
74 | },
|
---|
75 |
|
---|
76 | VeN3 == {
|
---|
77 | ParameterType -> External,
|
---|
78 | BlockName -> NUMIXING,
|
---|
79 | OrderBlock -> 3,
|
---|
80 | Value -> 0.0,
|
---|
81 | ComplexParameter -> False,
|
---|
82 | TeX -> Subscript[V,eN3],
|
---|
83 | Description -> "Mixing between ve flavor/gauge state and N3 mass state"
|
---|
84 | },
|
---|
85 |
|
---|
86 | VmuN1 == {
|
---|
87 | ParameterType -> External,
|
---|
88 | BlockName -> NUMIXING,
|
---|
89 | OrderBlock -> 4,
|
---|
90 | Value -> 0.0,
|
---|
91 | ComplexParameter -> False,
|
---|
92 | TeX -> Subscript[V,muN1],
|
---|
93 | Description -> "Mixing between vm flavor/gauge state and N1 mass state"
|
---|
94 | },
|
---|
95 |
|
---|
96 | VmuN2 == {
|
---|
97 | ParameterType -> External,
|
---|
98 | BlockName -> NUMIXING,
|
---|
99 | OrderBlock -> 5,
|
---|
100 | Value -> 1.0,
|
---|
101 | ComplexParameter -> False,
|
---|
102 | TeX -> Subscript[V,muN2],
|
---|
103 | Description -> "Mixing between vm flavor/gauge state and N2 mass state"
|
---|
104 | },
|
---|
105 |
|
---|
106 | VmuN3 == {
|
---|
107 | ParameterType -> External,
|
---|
108 | BlockName -> NUMIXING,
|
---|
109 | OrderBlock -> 6,
|
---|
110 | Value -> 0.0,
|
---|
111 | ComplexParameter -> False,
|
---|
112 | TeX -> Subscript[V,muN3],
|
---|
113 | Description -> "Mixing between vm flavor/gauge state and N3 mass state"
|
---|
114 | },
|
---|
115 |
|
---|
116 | VtaN1 == {
|
---|
117 | ParameterType -> External,
|
---|
118 | BlockName -> NUMIXING,
|
---|
119 | OrderBlock -> 7,
|
---|
120 | Value -> 0.0,
|
---|
121 | ComplexParameter -> False,
|
---|
122 | TeX -> Subscript[V,taN1],
|
---|
123 | Description -> "Mixing between vt flavor/gauge state and N1 mass state"
|
---|
124 | },
|
---|
125 |
|
---|
126 | VtaN2 == {
|
---|
127 | ParameterType -> External,
|
---|
128 | BlockName -> NUMIXING,
|
---|
129 | OrderBlock -> 8,
|
---|
130 | Value -> 0.0,
|
---|
131 | ComplexParameter -> False,
|
---|
132 | TeX -> Subscript[V,taN2],
|
---|
133 | Description -> "Mixing between vt flavor/gauge state and N2 mass state"
|
---|
134 | },
|
---|
135 |
|
---|
136 | VtaN3 == {
|
---|
137 | ParameterType -> External,
|
---|
138 | BlockName -> NUMIXING,
|
---|
139 | OrderBlock -> 9,
|
---|
140 | Value -> 1.0,
|
---|
141 | ComplexParameter -> False,
|
---|
142 | TeX -> Subscript[V,taN3],
|
---|
143 | Description -> "Mixing between vt flavor/gauge state and N3 mass state"
|
---|
144 | },
|
---|
145 |
|
---|
146 | (* Internal Parameters *)
|
---|
147 | gN == {
|
---|
148 | ParameterType -> Internal,
|
---|
149 | Definitions -> {gN->ee/sw},
|
---|
150 | InteractionOrder -> {NP,1},
|
---|
151 | TeX -> Subscript[g,N]}
|
---|
152 | };
|
---|
153 |
|
---|
154 | M$InteractionOrderHierarchy = {
|
---|
155 | {QCD,1},
|
---|
156 | {QED,2}
|
---|
157 | };
|
---|
158 |
|
---|
159 | (* ************************** *)
|
---|
160 | (* **** Particle classes **** *)
|
---|
161 | (* ************************** *)
|
---|
162 | M$ClassesDescription = {
|
---|
163 | (*Dirac Neutrino*)
|
---|
164 | F[131] == {
|
---|
165 | ClassName -> N1,
|
---|
166 | SelfConjugate -> False,
|
---|
167 | QuantumNumbers -> {LeptonNumber -> 1},
|
---|
168 | Mass -> {mN1,300.},
|
---|
169 | Width -> {WN1,0.303},
|
---|
170 | PropagatorLabel -> "N1",
|
---|
171 | PropagatorType -> Straight,
|
---|
172 | PropagatorArrow -> True,
|
---|
173 | ParticleName -> "N1",
|
---|
174 | PDG -> {9990012},
|
---|
175 | FullName -> "N1"},
|
---|
176 |
|
---|
177 | F[132] == {
|
---|
178 | ClassName -> N2,
|
---|
179 | SelfConjugate -> False,
|
---|
180 | QuantumNumbers -> {LeptonNumber -> 1},
|
---|
181 | Mass -> {mN2,500.},
|
---|
182 | Width -> {WN2,1.50},
|
---|
183 | PropagatorLabel -> "N2",
|
---|
184 | PropagatorType -> Straight,
|
---|
185 | PropagatorArrow -> True,
|
---|
186 | ParticleName -> "N2",
|
---|
187 | PDG -> {9990014},
|
---|
188 | FullName -> "N2"},
|
---|
189 |
|
---|
190 | F[133] == {
|
---|
191 | ClassName -> N3,
|
---|
192 | SelfConjugate -> False,
|
---|
193 | QuantumNumbers -> {LeptonNumber -> 1},
|
---|
194 | Mass -> {mN3,1000.},
|
---|
195 | Width -> {WN3,12.3},
|
---|
196 | PropagatorLabel -> "N3",
|
---|
197 | PropagatorType -> Straight,
|
---|
198 | PropagatorArrow -> True,
|
---|
199 | ParticleName -> "N3",
|
---|
200 | PDG -> {9990016},
|
---|
201 | FullName -> "N3"}
|
---|
202 | };
|
---|
203 |
|
---|
204 |
|
---|
205 | (* ************************** *)
|
---|
206 | (* ***** Lagrangian ***** *)
|
---|
207 | (* ************************** *)
|
---|
208 | LNKin := I N1bar[s1].Ga[v,s1,s2].del[N1[s2],v] - mN1 N1bar[s1]N1[s1] \
|
---|
209 | + I N2bar[s1].Ga[v,s1,s2].del[N2[s2],v] - mN2 N2bar[s1]N2[s1] \
|
---|
210 | + I N3bar[s1].Ga[v,s1,s2].del[N3[s2],v] - mN3 N3bar[s1]N3[s1];
|
---|
211 |
|
---|
212 | (* Charge Current *)
|
---|
213 | LNCCbare := gN/Sqrt[2] *( VeN1 * N1bar.W[m].ProjM[m].e \
|
---|
214 | + VmuN1 * N1bar.W[m].ProjM[m].mu \
|
---|
215 | + VtaN1 * N1bar.W[m].ProjM[m].ta ) \
|
---|
216 | + gN/Sqrt[2] *( VeN2 * N2bar.W[m].ProjM[m].e \
|
---|
217 | + VmuN2 * N2bar.W[m].ProjM[m].mu \
|
---|
218 | + VtaN2 * N2bar.W[m].ProjM[m].ta ) \
|
---|
219 | + gN/Sqrt[2] *( VeN3 * N3bar.W[m].ProjM[m].e \
|
---|
220 | + VmuN3 * N3bar.W[m].ProjM[m].mu \
|
---|
221 | + VtaN3 * N3bar.W[m].ProjM[m].ta );
|
---|
222 | LNCC := LNCCbare + HC[LNCCbare];
|
---|
223 |
|
---|
224 | (* Neutral Current *)
|
---|
225 | LNNCBare := 1/2 * gN/cw *( VeN1 * N1bar.Z[m].ProjM[m].ve \
|
---|
226 | + VmuN1 * N1bar.Z[m].ProjM[m].vm \
|
---|
227 | + VtaN1 * N1bar.Z[m].ProjM[m].vt ) \
|
---|
228 | + 1/2 * gN/cw *( VeN2 * N2bar.Z[m].ProjM[m].ve \
|
---|
229 | + VmuN2 * N2bar.Z[m].ProjM[m].vm \
|
---|
230 | + VtaN2 * N2bar.Z[m].ProjM[m].vt ) \
|
---|
231 | + 1/2 * gN/cw *( VeN3 * N3bar.Z[m].ProjM[m].ve \
|
---|
232 | + VmuN3 * N3bar.Z[m].ProjM[m].vm \
|
---|
233 | + VtaN3 * N3bar.Z[m].ProjM[m].vt );
|
---|
234 | LNNC := LNNCBare + HC[LNNCBare];
|
---|
235 |
|
---|
236 | (* Higgs Interaction *)
|
---|
237 | LNHbare := - gN*mN1/(2*MW) *( VeN1 * N1bar.ProjM.ve H \
|
---|
238 | + VmuN1 * N1bar.ProjM.vm H \
|
---|
239 | + VtaN1 * N1bar.ProjM.vt H ) \
|
---|
240 | - gN*mN2/(2*MW) *( VeN2 * N2bar.ProjM.ve H \
|
---|
241 | + VmuN2 * N2bar.ProjM.vm H \
|
---|
242 | + VtaN2 * N2bar.ProjM.vt H ) \
|
---|
243 | - gN*mN3/(2*MW) *( VeN3 * N3bar.ProjM.ve H \
|
---|
244 | + VmuN3 * N3bar.ProjM.vm H \
|
---|
245 | + VtaN3 * N3bar.ProjM.vt H );
|
---|
246 | LNH := LNHbare + HC[LNHbare];
|
---|
247 |
|
---|
248 | (* Goldstone Interaction *)
|
---|
249 | LNGbare := I *gN*mN1/(2*MW) *( VeN1 * vebar.ProjP.N1 G0 \
|
---|
250 | + VmuN1 * vmbar.ProjP.N1 G0 \
|
---|
251 | + VtaN1 * vtbar.ProjP.N1 G0) \
|
---|
252 | + I *gN*mN2/(2*MW) *( VeN2 * vebar.ProjP.N2 G0 \
|
---|
253 | + VmuN2 * vmbar.ProjP.N2 G0 \
|
---|
254 | + VtaN2 * vtbar.ProjP.N2 G0) \
|
---|
255 | + I *gN*mN3/(2*MW) *( VeN3 * vebar.ProjP.N3 G0 \
|
---|
256 | + VmuN3 * vmbar.ProjP.N3 G0 \
|
---|
257 | + VtaN3 * vtbar.ProjP.N3 G0) \
|
---|
258 | + I *gN*mN1/(Sqrt[2]*MW) *( VeN1 * ebar.ProjP.N1 GPbar \
|
---|
259 | + VmuN1 * mubar.ProjP.N1 GPbar \
|
---|
260 | + VtaN1 * tabar.ProjP.N1 GPbar) \
|
---|
261 | + I *gN*mN2/(Sqrt[2]*MW) *( VeN2 * ebar.ProjP.N2 GPbar \
|
---|
262 | + VmuN2 * mubar.ProjP.N2 GPbar \
|
---|
263 | + VtaN2 * tabar.ProjP.N2 GPbar) \
|
---|
264 | + I *gN*mN3/(Sqrt[2]*MW) *( VeN3 * ebar.ProjP.N3 GPbar \
|
---|
265 | + VmuN3 * mubar.ProjP.N3 GPbar \
|
---|
266 | + VtaN3 * tabar.ProjP.N3 GPbar);
|
---|
267 | LNG := LNGbare + HC[LNGbare];
|
---|
268 |
|
---|
269 | (* Full N Lagrangian *)
|
---|
270 | LN:= LNKin + LNCC + LNNC + LNH + LNG;
|
---|
271 |
|
---|
272 | (* Combine Lagrangian *)
|
---|
273 | LFull := LSM + LN;
|
---|