HeavyN: heavyN_Dirac.fr

File heavyN_Dirac.fr, 9.0 KB (added by Richard Ruiz, 6 years ago)

Model flie for HeavyN_Dirac

Line 
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(* ************************** *)
26M$ModelName = "SM_HeavyN_Dirac";
27
28M$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};
39FeynmanGauge = 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(* ************************** *)
53M$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
154M$InteractionOrderHierarchy = {
155 {QCD,1},
156 {QED,2}
157};
158
159(* ************************** *)
160(* **** Particle classes **** *)
161(* ************************** *)
162M$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(* ************************** *)
208LNKin := 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 *)
213LNCCbare := 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 );
222LNCC := LNCCbare + HC[LNCCbare];
223
224(* Neutral Current *)
225LNNCBare := 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 );
234LNNC := LNNCBare + HC[LNNCBare];
235
236(* Higgs Interaction *)
237LNHbare := - 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 );
246LNH := LNHbare + HC[LNHbare];
247
248(* Goldstone Interaction *)
249LNGbare := 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);
267LNG := LNGbare + HC[LNGbare];
268
269(* Full N Lagrangian *)
270LN:= LNKin + LNCC + LNNC + LNH + LNG;
271
272(* Combine Lagrangian *)
273LFull := LSM + LN;