zero SI and SD cross section

[waqasahmed]

Dear Florian,

Hope you are fine and doing well. I have one question related to SI and SD cross sections. In SSP scan, I open a Flag 50->0 as you mention in the manual when we include dark matter calculation. But I found in my data many points have zero crossection printed. I want to know what is the reason for this. I paste blow the output of one point.

Thanking you,
best regards,
Waqas




Block DARKMATTER # MicrOmegas results
1 0.00515 # relic density
100 0.185054 # ~su1 ~su2 -> u1 u2
101 0.122968 # ~su1 ~Su1 -> g g
102 0.12231 # ~su2 ~Su2 -> g g
103 0.082034 # ~su1 ~su1 -> u1 u1
104 0.081588 # ~su2 ~su2 -> u2 u2
105 0.041001 # ~sd1 ~su1 -> d1 u1
106 0.040948 # ~sd2 ~su2 -> d2 u2
107 0.028747 # ~sd2 ~su1 -> d2 u1
108 0.028743 # ~sd1 ~su2 -> d1 u2
109 0.024552 # ~su1 ~Su1 -> g Z
110 0.024416 # ~su2 ~Su2 -> g Z
111 0.022962 # ~sd1 ~Su1 -> g Wm
112 0.022932 # ~sd2 ~Su2 -> g Wm
113 0.019208 # ~sd3 ~su1 -> d3 u1
114 0.019206 # ~sd3 ~su2 -> d3 u2
115 0.014403 # ~su1 ~Su1 -> g A
116 0.014323 # ~su2 ~Su2 -> g A
201 0. #
202 0. #
203 0. #
204 0. #
301 0. #
302 0. #
303 0. #
304 0. #

[FStaub]

At least this point has obviously a squark LSP which can't be DM.

[waqasahmed]

Dear Florian,
Thanks for your guidance. I have another question related to NHSSM model given in Sarah. I work at low scale. When I consider trilinear terms > 1 Gev than i have error like that

Calculating mass spectrum
1 .-iteration
Loop corrections to Yukawa couplings at renormalisation scale too large!
Note: The following floating-point exceptions are signalling: IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP Subroutine TerminateProgram


could you please tell me how to fixed that one.

here I attached the model files.

best regards,
Waqas

[FStaub]

Hi,

there is nothing really to be fixed. It just tells you that for your chosen parameters perturbation theory is no longer trustworthy. You can try what happens when you turn off these checks, but I can't guarantee that the output is trustworthy.
If you want to try anyway, remove the following from Boundary_MODEL.f90

Code: Select all

!-------------------------------------------------- 
!Either we have run into a numerical problem or 
!perturbation theory breaks down 
!-------------------------------------------------- 
If ((Minval(Abs(mf_l_Q/mf_l)).Lt.0.1_dp)&
&.Or.(Maxval(Abs(mf_l_Q/mf_l)).Gt.10._dp)) Then
Iname=Iname-1
kont=-405
Call AddError(405)
    Write(*,*) " Loop corrections to Yukawa couplings at EW scale too large!" 
    Call TerminateProgram
Else If ((Minval(Abs(mf_d_Q/mf_d)).Lt.0.1_dp)&
&.Or.(Minval(Abs(mf_d_Q/mf_d)).Gt.10._dp)) Then
Iname=Iname-1
kont=-406
Call AddError(406)
    Write(*,*) " Loop corrections to Yukawa couplings at EW scale too large!" 
    Call TerminateProgram
Else If ((Minval(Abs(mf_u_Q/mf_u)).Lt.0.1_dp)&
&.Or.(Minval(Abs(mf_u_Q/mf_u)).Gt.10._dp)) Then
Iname=Iname-1
kont=-407
Call AddError(407)
    Write(*,*) " Loop corrections to Yukawa couplings at EW scale too large!" 
    Call TerminateProgram
End If