Skip to content

Commit

Permalink
adding THDM, THDM + Higgsinos + THDM + split models
Browse files Browse the repository at this point in the history 
including 1- and 2-loop boundary conditions to the MSSM from
arxiv:1508.00576 .
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Jan 14, 2016
1 parent b426f01 commit 1faaab7
Show file tree
Hide file tree
Showing 13 changed files with 1,836 additions and 0 deletions.
319 changes: 319 additions & 0 deletions model_files/HGTHDMIIMSSMBC/FlexibleSUSY.m.in
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,319 @@
FSModelName = "@CLASSNAME@";
FSEigenstates = SARAH`EWSB;
AutomaticInputAtMSUSY = False;
FSDefaultSARAHModel = "HGTHDM-II";

(* input parameters *)

MINPAR = {
{3, TanBeta}
};

EXTPAR = {
{0, MSUSY},
{1, MEWSB},
{2, MuInput},
{3, M1Input},
{4, M2Input},
{5, M3Input},
{6, MAInput},
{7, AtInput},
{8, AbInput},
{9, AtauInput},
{100, LambdaLoopOrder}
};

EWSBOutputParameters = { M112, M222 };

(* The high scale where we match to the MSSM *)
HighScale = MSUSY;

HighScaleFirstGuess = MSUSY;

HighScaleInput = {
{g1d , g2 },
{g1dp , GUTNormalization[g1] g1 (* gY *)},
{g2u , g2 },
{g2up , GUTNormalization[g1] g1 (* gY *)},
{Lambda1, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda1th1L + deltaLambda1Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda1th2L)},
{Lambda2, 1/2 (1/4 ( (GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda2th1L + deltaLambda2Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda2th2L)},
{Lambda3, 1/4 (-(GUTNormalization[g1] g1)^2 + g2^2)
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda3th1L + deltaLambda3Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda3th2L},
{Lambda4, -1/2 g2^2
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda4th1L + deltaLambda4Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda4th2L},
{Lambda5, 0
+ UnitStep[THRESHOLD-1] UnitStep[LambdaLoopOrder-1] (deltaLambda5th1L + deltaLambda5Phi1L)
+ UnitStep[THRESHOLD-2] UnitStep[LambdaLoopOrder-2] deltaLambda5th2L}
};

(* The scale where we impose the EWSB conditions
and calculate the spectrum *)
SUSYScale = MEWSB;

SUSYScaleFirstGuess = MEWSB;

SUSYScaleInput = {
{\[Mu] , MuInput },
{MassB , M1Input },
{MassWB , M2Input },
{MassG , M3Input },
{M122 , MAInput^2 Sin[ArcTan[v2/v1]] Cos[ArcTan[v2/v1]]}
};

LowScale = LowEnergyConstant[MT];

LowScaleFirstGuess = LowEnergyConstant[MT];

LowScaleInput = {
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{v1, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Cos[ArcTan[TanBeta]]},
{v2, 2 MZDRbar / Sqrt[GUTNormalization[g1]^2 g1^2 + g2^2] Sin[ArcTan[TanBeta]]}
};

InitialGuessAtLowScale = {
{v1, LowEnergyConstant[vev] Cos[ArcTan[TanBeta]]},
{v2, LowEnergyConstant[vev] Sin[ArcTan[TanBeta]]},
{Yu, Automatic},
{Yd, Automatic},
{Ye, Automatic},
{\[Mu] , MuInput },
{MassB , M1Input },
{MassWB , M2Input },
{MassG , M3Input },
{M122 , MAInput^2 Sin[ArcTan[TanBeta]] Cos[ArcTan[TanBeta]]}
};

DefaultPoleMassPrecision = MediumPrecision;
HighPoleMassPrecision = {hh};
MediumPoleMassPrecision = {};
LowPoleMassPrecision = {};

SMParticles = {
Electron, TopQuark, BottomQuark,
VectorP, VectorZ, VectorG, VectorW, Neutrino
};

(* abbreviations *)
At = AtInput;
Ab = AbInput;
Atau = AtauInput;
Lambda1WagnerLee = 2 Lambda1;
Lambda2WagnerLee = 2 Lambda2;

(* arxiv:1508.00576, Eq. (45) *)
deltaLambda1th1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY
},
(
- kappa/2 ht^4 muMS^4
+ 6 kappa hb^4 AbMS^2 (1 - AbMS^2/12)
+ 2 kappa htau^4 AtauMS^2 (1 - AtauMS^2/12)
+ kappa (g2^2 + gY^2)/4 (3 ht^2 muMS^2 - 3 hb^2 AbMS^2 - htau^2 AtauMS^2)
)
];

(* arxiv:1508.00576, Eq. (46) *)
deltaLambda2th1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY,
AtMS = At / MSUSY
},
(
6 kappa ht^4 AtMS^2 (1 - AtMS^2/12)
- kappa/2 hb^4 muMS^4
- kappa/6 htau^4 muMS^4
- kappa (g2^2 + gY^2)/4 (3 ht^2 AtMS^2 - 3 hb^2 muMS^2 - htau^2 muMS^2)
)
];

(* arxiv:1508.00576, Eq. (47)-(48) *)
deltaLambda3th1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY,
AtMS = At / MSUSY
},
(
kappa/6 muMS^2 (3 ht^4 (3 - AtMS^2)
+ 3 hb^4 (3 - AbMS^2)
+ htau^4 (3 - AtauMS^2))
+ kappa/2 ht^2 hb^2 (3 (AtMS + AbMS)^2
- (muMS^2 - AtMS AbMS)^2
- 6muMS^2)
- kappa/2 (g2^2 - gY^2)/4 (3 ht^2 (AtMS^2 - muMS^2)
+ 3 hb^2 (AbMS^2 - muMS^2)
+ htau^2 (AtauMS^2 - muMS^2))
)
];

(* arxiv:1508.00576, Eq. (49) *)
deltaLambda4th1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY,
AtMS = At / MSUSY
},
(
kappa/6 muMS^2 (3 ht^4 (3 - AtMS^2)
+ 3 hb^4 (3 - AbMS^2)
+ htau^4 (3 - AtauMS^2))
- kappa/2 ht^2 hb^2 (3 (AtMS + AbMS)^2
- (muMS^2 - AtMS AbMS)^2
- 6 muMS^2)
+ kappa/2 g2^2/2 (3 ht^2 (AtMS^2 - muMS^2)
+ 3 hb^2 (AbMS^2 - muMS^2)
+ htau^2 (AtauMS^2 - muMS^2))
)
];

(* arxiv:1508.00576, Eq. (50) *)
deltaLambda5th1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY,
AtMS = At / MSUSY
},
(
- kappa/6 muMS^2 (3 ht^4 AtMS^2 + 3 hb^4 AbMS^2 + htau^4 AtauMS^2)
)
];

(* arxiv:1508.00576, Eq. (53) *)
deltaLambda1Phi1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY
},
(
-kappa/6 (g2^2 + gY^2)/2 (3 ht^2 muMS^2 + 3hb^2 AbMS^2 + htau^2 AtauMS^2)
)
];

(* arxiv:1508.00576, Eq. (54) *)
deltaLambda2Phi1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AtMS = At / MSUSY
},
(
-kappa/6 (g2^2 + gY^2)/2 (3 ht^2 AtMS^2 + 3 hb^2 muMS^2 + htau^2 muMS^2)
)
];

(* arxiv:1508.00576, Eq. (55) *)
deltaLambda3Phi1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
gY = GUTNormalization[g1] g1,
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY,
AtMS = At / MSUSY
},
(
-kappa/6 (g2^2 - gY^2)/4 (3 ht^2 (AtMS^2 + muMS^2) + 3 hb^2 (AbMS^2 + muMS^2) + htau^2 (AtauMS^2 + muMS^2))
)
];

(* arxiv:1508.00576, Eq. (56) *)
deltaLambda4Phi1L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
hb = Yd[3,3]/Cos[ArcTan[v2/v1]],
htau = Ye[3,3]/Cos[ArcTan[v2/v1]],
muMS = \[Mu] / MSUSY,
AbMS = Ab / MSUSY,
AtauMS = Atau / MSUSY,
AtMS = At / MSUSY
},
(
kappa/6 g2^2/2 (3 ht^2 (AtMS^2 + muMS^2) + 3 hb^2 (AbMS^2 + muMS^2) + htau^2 (AtauMS^2 + muMS^2))
)
];

(* arxiv:1508.00576, Eq. (57) *)
deltaLambda5Phi1L = 0;

(* arxiv:1508.00576, Eq. (59) *)
deltaLambda1th2L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
muMS = \[Mu] / MSUSY
},
(
-4/3 kappa^2 ht^4 g3^2 muMS^4
)
];

(* arxiv:1508.00576, Eq. (60) *)
deltaLambda2th2L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
muMS = \[Mu] / MSUSY,
AtMS = At / MSUSY
},
(
16 kappa^2 ht^4 g3^2 (-2 AtMS + 1/3 AtMS^3 - 1/12 AtMS^4)
)
];

(* arxiv:1508.00576, Eq. (61) *)
deltaLambda3th2L = With[{
kappa = 1/(4 Pi)^2,
ht = Yu[3,3]/Sin[ArcTan[v2/v1]],
muMS = \[Mu] / MSUSY,
AtMS = At / MSUSY
},
(
4 kappa^2 ht^4 g3^2 AtMS muMS^2 (1 - 1/2 AtMS)
)
];

deltaLambda4th2L = deltaLambda3th2L;
deltaLambda5th2L = deltaLambda3th2L;
49 changes: 49 additions & 0 deletions model_files/HGTHDMIIMSSMBC/LesHouches.in.HGTHDMIIMSSMBC
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,49 @@
Block MODSEL # Select model
# 12 1000 # DRbar parameter output scale (GeV)
Block FlexibleSUSY
0 1.000000000e-06 # precision goal
1 0 # max. iterations (0 = automatic)
2 0 # algorithm (0 = two_scale, 1 = lattice)
3 1 # calculate SM pole masses
4 2 # pole mass loop order
5 2 # EWSB loop order
6 2 # beta-functions loop order
7 2 # threshold corrections loop order
8 1 # Higgs 2-loop corrections O(alpha_t alpha_s)
9 1 # Higgs 2-loop corrections O(alpha_b alpha_s)
10 1 # Higgs 2-loop corrections O(alpha_t^2 + alpha_t alpha_b + alpha_b^2)
11 1 # Higgs 2-loop corrections O(alpha_tau^2)
12 0 # force output
13 1 # Top quark 2-loop corrections QCD
Block SMINPUTS # Standard Model inputs
1 1.279440000e+02 # alpha^(-1) SM MSbar(MZ)
2 1.166380000e-05 # G_Fermi
3 1.184000000e-01 # alpha_s(MZ) SM MSbar
4 9.118760000e+01 # MZ(pole)
5 4.180000000e+00 # mb(mb) SM MSbar
6 1.733400000e+02 # mtop(pole)
7 1.777000000e+00 # mtau(pole)
8 0.000000000e+00 # mnu3(pole)
9 80.425 # MW pole
11 5.109989020e-04 # melectron(pole)
12 0.000000000e+00 # mnu1(pole)
13 1.056583570e-01 # mmuon(pole)
14 0.000000000e+00 # mnu2(pole)
21 4.750000000e-03 # md(2 GeV) MS-bar
22 2.400000000e-03 # mu(2 GeV) MS-bar
23 1.040000000e-01 # ms(2 GeV) MS-bar
24 1.270000000e+00 # mc(mc) MS-bar
Block MINPAR
3 3 # TanBeta(MZ)
Block EXTPAR # Input parameters
0 1e6 # MSUSY
1 173.34 # MEWSB
2 200 # Mu(MEWSB)
3 200 # M1(MEWSB)
4 200 # M2(MEWSB)
5 200 # M3(MEWSB)
6 200 # MA(MEWSB)
7 0 # At(MSUSY)
8 0 # Ab(MSUSY)
9 0 # Atau(MSUSY)
100 2 # lambda BCs loop order
Loading

0 comments on commit 1faaab7

Please sign in to comment.