Complex neutrino Yukawa

[avirupghosh93]

Thanks Florian,
For your reply. But there is a different error with real Yukawa. I have tried it with linear-Seesaw model files given in SARAH.....the same problem persists there even in real Yukawa case.
In the real case the error is:-
"Problem in RGE Running. Errorcode: -1006"
While in case of imaginary coupling :-
"NaN appearing in CalculateMFv"
I have tested the used set of Yukawa coupling to give finite and positive values for light neutrino masses in Mathematica.
Thanks in advance,
Avirup

[FStaub]

The error '1006' points towards a problem with QP precision. Did you really change the makefile and re-compiled everything afterwards?

[avirupghosh93]

Thanks Florian,
For your reply. I had compiled the linear seesaw after editing src/makefile .
One interesting observation is that for MR set to zero the code is running for YVIN=diag{1.0E-01,1.0E-01,1.0E-01}.
But when it is added to IMYVIN =diag{1.0E-01,1.0E-01,1.0E-01} Spheno is terminating with the error "Problem in RGE running: -1006". It is not also running when MR is set non-zero. Three ellustrative values of MR are 1E+03,1E+14,1E+16.
Other parameters like YSL and BMR are kept zero throughout.
Thanks,
Avirup

[FStaub]

Can you try more realistic numbers:
1) Imaginary part smaller than real
2) all entries of MR not bigger than 10^4

[avirupghosh93]

Hi Florian,
Is there any such constraint that real part of Yukawa has to be larger than imaginary part?
I have tried with following Yukawa's:-
(1) Yv=diag(1E-01, i 1E-05, 1E-01)
The error is "Nan appearing in CalculateMFv"
(2) Yv=diag(1E-00, i 1E-05, 1E-01)
The error is "Problem with RGE running.Errorcode : -12"
In both cases MR=diag(100,100,100)
Does NaN implies negative mass squared?
Thanks
Avirup

[FStaub]

This input works fine for me:

Code: Select all

Block MRIN    #  
1 1   1.0000000E+03         # MR(1,1)
1 2   0.000000E+00         # MR(1,2)
1 3   0.000000E+00         # MR(1,3)
2 1   0.000000E+00         # MR(2,1)
2 2   1.0000000E+03         # MR(2,2)
2 3   0.000000E+00         # MR(2,3)
3 1   0.000000E+00         # MR(3,1)
3 2   0.000000E+00         # MR(3,2)
3 3   1.0000000E+03         # MR(3,3)
Block YSLIN    #  
1 1   0.1000000E+00         # YSL(1,1)
1 2   0.000000E+00         # YSL(1,2)
1 3   0.000000E+00         # YSL(1,3)
2 1   0.000000E+00         # YSL(2,1)
2 2   0.1000000E+00         # YSL(2,2)
2 3   0.000000E+00         # YSL(2,3)
3 1   0.000000E+00         # YSL(3,1)
3 2   0.000000E+00         # YSL(3,2)
3 3   0.1000000E+00         # YSL(3,3)
Block YVIN    #  
1 1   1.0000000E-01         # Yv(1,1)
1 2   0.0000000E+00         # Yv(1,2)
1 3   0.000000E+00         # Yv(1,3)
2 1   0.000000E+00         # Yv(2,1)
2 2   1.0000000E-01         # Yv(2,2)
2 3   0.000000E+00         # Yv(2,3)
3 1   0.000000E+00         # Yv(3,1)
3 2   0.000000E+00         # Yv(3,2)
3 3   1.0000000E-01         # Yv(3,3)
Block IMYVIN    #  
1 1   1.0000000E-02         # Yv(1,1)
1 2   0.000000E+00         # Yv(1,2)
1 3   0.000000E+00         # Yv(1,3)
2 1   0.000000E+00         # Yv(2,1)
2 2   1.0000000E-02         # Yv(2,2)
2 3   0.000000E+00         # Yv(2,3)
3 1   0.000000E+00         # Yv(3,1)
3 2   0.000000E+00         # Yv(3,2)
3 3   1.0000000E-02         # Yv(3,3)

So, either it's a problem with your point or with your installation. in the second case, a

Code: Select all

make cleanall

might help.

Florian

[avirupghosh93]

Thanks Florian,
It worked after doing "make cleanall"......
Thanks,
Avirup

[avirupghosh93]

Hi Florian,
SPheno works perfectly well for Yukawa couplings for which the eigenvalues of squared mass matrix of neutrino has positive eigenvalues. What modifications do I need to do for the set of Yukawas for which the mass squared eigenvalues are negative.....
Also for certain set of Yukawas the mass squared matrix has positive eigenvalues but SPheno stops running because of Problem in RGE running.....does such problem also occurs for only low scale input i.e if I use
LesHouches.in.LinSeesaw_low?
Any suggestions will be really helpful.
Thanks,
Avirup

[FStaub]

What modifications do I need to do for the set of Yukawas for which the mass squared eigenvalues are negative.....

Nothing. That's an ill defined point with tachyonic states.

does such problem also occurs for only low scale input i.e if I use
LesHouches.in.LinSeesaw_low?

You'll know when you try...

[SDrau]

Dear Florian,

I'll join in on questioning about QP precision. In the wiki it is said, that one needs to enable QP precision in the model file and I can't seem to find that statement anywhere else. Is it enough to delete the -DONLYDOUBLE flag and recompile, or should one also use, for example, QuadruplePrecision={Fv}; in the model file?

Regards,

S.Drauksas