TeX-Output with vector VEVs

Questions concerning the different calculations performed by SARAH and the output of the information into other formats
Post Reply
Martin
Posts: 15
Joined: 10. Feb 2017, 07:31

TeX-Output with vector VEVs

Post by Martin » 6. Jun 2017, 13:15

Hi,

I want to study models where fields which appear in multiple generations receive VEVs. But I do not manage to get a very simple TeX-Output. For example I took the NMSSM model file. There I defined the singlet to appear in three generations

SuperFields[[8]] = {s, 3, sR, 0, 1, 1, RpP, Exp[2*Pi*I/3]};

For the definition of the VEV I did not change anything

{SsR, {vS, 1/Sqrt[2]}, {sigmaS, I/Sqrt[2]},{phiS,1/Sqrt[2]}}

Now when I try to generate some Output using this command

ModelOutput[EWSB, WriteTeX -> True, FeynmanDiagrams -> False,
IncludeVertices -> False, TwoLoopRGEs -> False,
SimplifySums -> False]

it starts to run but after "Writing Mass Matrices to TeX-File" it keeps on running for hours but nothing happens anymore. Looking at the TeX-files they are also pretty much incomplete. The scalar mass matrix stops in the middle.

I would be very glad if someone could help.

Thanks,
Martin

FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: TeX-Output with vector VEVs

Post by FStaub » 6. Jun 2017, 13:43

Hi Martin,

I would need to check what exactly goes on, but have not time before Friday. From your description I guess it comes from the simplification of the mass matrices during the TeX output. My main suspect is this line in Package/Outputs/tex.m

Code: Select all

value = FullSimplify[MMatrizen[[i,j,k]] /. ThetaStep[a___]\[Rule]1] //. subTeXMM //. Delta[ToExpression["gm"<>ToString[k]],ToExpression["gn"<>ToString[j]]]\[Rule]UnitM //. A_[ToExpression["gm"<>ToString[j]],ToExpression["gn"<>ToString[k]]]\[Rule] A //. A_[ToExpression["gn"<>ToString[k]],ToExpression["gm"<>ToString[j]]]\[Rule]Tp[A] //. Conj\[Rule]conj //. a_ conj[x_] + a_ x_ \[Rule]2 a Re[x] //. a_ x_ + a_ conj[x_] \[Rule]2 a Re[x] //. conj[x_] +  x_ \[Rule]2 Re[x]//. x_ + conj[x_]\[Rule] 2 Re[x] //. conj[x_]x_ \[Rule] abs[x] //. x_ conj[x_]\[Rule]abs[x] //. conj[Tp[x_]]\[Rule]Adj[x] //. Tp[conj[x_]]\[Rule]Adj[x] //. conj[Adj[x_]]\[Rule]Tp[x] //. Adj[conj[x_]]\[Rule]Tp[x] //.(a_-b_)(a_+b_)\[Rule](a^2-b^2) //. subRealSum //.RXl[_]\[Rule]0;
You could try to replace all //. just by /.

Cheers,
Florian

Martin
Posts: 15
Joined: 10. Feb 2017, 07:31

Re: TeX-Output with vector VEVs

Post by Martin » 7. Jun 2017, 04:06

Hi Florian,

No worries. It is not super urgent. At the moment I try to guess the right mass matrices from one generation only but it would be much easier of course to check if I did not do any mistake.

I was also thinking that it might be due to some simplification but I hand over the option SimplifySums -> False

Code: Select all

ModelOutput[EWSB, WriteTeX -> True, FeynmanDiagrams -> False, 
 IncludeVertices -> False, TwoLoopRGEs -> False, 
 SimplifySums -> False]
 
In fact, in Package/Outputs/tex.m the code you mention is commented out.

Thanks for your help,
Martin

FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: TeX-Output with vector VEVs

Post by FStaub » 7. Jun 2017, 09:27

Martin wrote: In fact, in Package/Outputs/tex.m the code you mention is commented out.
This part should appear twice, only once commented out as you said.

Martin
Posts: 15
Joined: 10. Feb 2017, 07:31

Re: TeX-Output with vector VEVs

Post by Martin » 7. Jun 2017, 09:36

In my code it appears only once :?

FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: TeX-Output with vector VEVs

Post by FStaub » 7. Jun 2017, 09:45

Ah sorry, the lines are not completely identical. It's easiest to search for subTeXMM" which appears in this context in SARAH-4.11.0 in the lines 1287 and 1537.

Cheers,
Florian

Martin
Posts: 15
Joined: 10. Feb 2017, 07:31

Re: TeX-Output with vector VEVs

Post by Martin » 7. Jun 2017, 10:08

Ok. I have replaced all //. with /. in lines 1287 and 1537 but the problem still persists.

Thanks,
Martin

FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: TeX-Output with vector VEVs

Post by FStaub » 7. Jun 2017, 10:39

Ok. A more radical solution is to change the line to

Code: Select all

value = matIN[[j, k]] /. ThetaStep[a___] -> 1 /. subTeXMM /. 
    Delta[ToExpression["gm" <> ToString[k]], 
      ToExpression["gn" <> ToString[j]]] -> UnitM /. 
   A_[ToExpression["gm" <> ToString[j]], 
     ToExpression["gn" <> ToString[k]]] -> A /. 
  A_[ToExpression["gn" <> ToString[k]], 
    ToExpression["gm" <> ToString[j]]] -> Tp[A]
  
Cheers,
Florian

Martin
Posts: 15
Joined: 10. Feb 2017, 07:31

Re: TeX-Output with vector VEVs

Post by Martin » 7. Jun 2017, 11:46

Ok. Now it runs through but the m_{\phi_S \phi_S} component of the Higgs boson matrix is 5 pages long... :shock:

FStaub
Site Admin
Posts: 822
Joined: 13. Apr 2016, 14:05

Re: TeX-Output with vector VEVs

Post by FStaub » 7. Jun 2017, 17:05

Hi,

that's not unexpected... in the \mu\nuSSM with 3-generations of \nu_R it looks even worse.

Cheers,
Florian

Post Reply