Dear experts, I have a model with an extra hidden U(1) gauge group and I'm trying to make Z' a dark matter candidate (odd under ExtraZ2). I have Global[[1]] = {Z[2], Z2}; Global[[2]] = {Z[2], ExtraZ2}; (* Gauge Superfields *) Gauge[[1]]={B, U[1], hypercharge, g1,False, 1, 1}; Gauge[[2]]={WB, SU[2], ...

Thanks Florian. I will check the paper.

Best,
Amin

Hi Florian,

Thanks for clarifying this. So the effect of the term I am looking for is encoded in g_ab which is the matrix {{g_1, \tilde{g}}, {0, g_BL}} you define in Eq. (4.9) for the B-L model in your paper arXiv:1503.04200 ?

Best,
Amin

Hi Florian, When adding an extra vector superfield C to the MSSM and allowing gauge kinetic mixing, I can see SARAH adds a soft term M\lambda_B\lambda_C (where B is the original vector superfield). Now, in the gauge sector of the theory, there are terms as: \delta*(\lambda_C\sigma^{\mu}\partial_{\mu...

Hi Florian, Got it. I will check the rotation matrices for this purpose. I should note that the main SPheno models (not from SARAH) do not have the same nomenclature. I mean, the sleptons, for example, are not named Se_1, Se_2, etc ... but rather ~e_L, ~e_R, etc ... and this eliminates the confusion...

Hi Florian, I noticed that in the output of SPheno, the particles are ranked in increasing order of mass regardless of their PDG codes. For example, left stau (Se_3) PDG code is 1000015 and that for left selectron (Se_1) is 1000011. For my model point, the stau is lighter than the selectron, so inst...

Hi Florian,

The problem turned out to be coming from the DependenceSPheno for ThetaWp, where it was defined as ArcSin instead of ArcCos. I fixed it and everything is working properly now.
Thank you for the help.

Best,
Amin

Hi Florian, Thanks for the tips. I set the matching conditions to tree level and now it works. However, the GUT scale is now 2x10^15 GeV rather than 10^16 GeV. From the theory side of the Stueckelberg mechanism, the effect of this sector is minimal on the MSSM when the gauge kinetic mixing is small....

Hi Florian, Thank you. I have one final question: the error causing the RGEs to stop running is due to SPheno not being able to satisfy g1==g2 at the GUT scale. Looking at the RGEs, g2 is not affected by gauge kinetic mixing; only g1 is. So how come g2 is turning out to be that large? Also, modifyin...

Hi again, In my model, the mixed gauge couplings, g_YB and g_BY, are related to the U(1) gauge coupling g_1 such that g_YB = -g_1*sin(delta), where delta is the coefficient of kinetic mixing. Is there a way to tell SARAH to write RGEs only for g_1, g_2, g_3 and g_p (new gauge coupling of U(1)') sinc...