Remove["Global`*"]
In[]:=
c1=3
10
10
;ω1g=50×2πc1
-15
10
;ω1e=39×2πc1
-15
10
;
ζ
1
=10.0×2πc1
-15
10
;
In[]:=
h=6.62610^-3410^15/(2π);k=1.3806610^-23;T=1.0;β=
1
kT
;Q=Csch[βhω1g/2]/2;
In[]:=
Ag=
ω1g
Sinh[hω1g(β-t/h)]
;M=1;Ae=
ω1e
Sinh[ω1et]
;
In[]:=
Cg=
ω1g
h
Coth[hω1g(β-t/h)/2];Ce=
ω1e
h
Coth[tω1e/2];
In[]:=
dg=
ω1g
h
Tanh[hω1g(β-t/h)/2];
In[]:=
de=
ω1e
h
Tanh[tω1e/2];
In[]:=
K=
h
ω1g
;
In[]:=
Δ=de+dg;ν=dgK;Λ=Ce+Cg;​​Ξ=-KdgK+νν/Δ;
In[]:=
PlotRe
1
Qh
Sqrt
AgAe
ΛΔ
Exp[Ξ]Exp[-
ζ
1
t/2],{t,0,3000},FrameTrue,FrameStyleThickness[0.006],AspectRatio1,BaseStyle{16,FontFamily"TimesRoman"},PlotStyle{Thick,Brown},FrameTicks{True,None},AxesNone,PlotRangeAll
Out[]=
Define your own log range for Log function
In[]:=
myLog[z_]:=If[Arg[z]<0,Log[z]+I2π,Log[z]]
In[]:=
PlotRe
1
Qh
Exp
1
2
myLog
AgAe
ΛΔ
Exp[Ξ]Exp[-
ζ
1
t/2],{t,0.00,3000},FrameTrue,AxesNone,FrameStyleThickness[0.006],AxesNone,AspectRatio1,BaseStyle{16,FontFamily"TimesRoman"},PlotStyle{Thick,Red},FrameTicks{True,None},PlotRangeAll
Out[]=

The above strategy does not work here

In[]:=
Remove["Global`*"]
In[]:=
c1=3
10
10
;ω1g=50×2πc1
-15
10
;ω1e=35×2πc1
-15
10
;
ζ
1
=10.0×2πc1
-15
10
;
In[]:=
h=6.62610^-3410^15/(2π);k=1.3806610^-23;T=1.0;β=
1
kT
;Q=Csch[βhω1g/2]/2;
In[]:=
Ag=
ω1g
Sinh[hω1g(β-t/h)]
;M=1;Ae=
ω1e
Sinh[ω1et]
;
In[]:=
Cg=
ω1g
h
Coth[hω1g(β-t/h)/2];Ce=
ω1e
h
Coth[tω1e/2];
In[]:=
dg=
ω1g
h
Tanh[hω1g(β-t/h)/2];
In[]:=
de=
ω1e
h
Tanh[tω1e/2];
In[]:=
K=
h
ω1g
;
In[]:=
Δ=de+dg;ν=dgK;Λ=Ce+Cg;​​Ξ=-KdgK+νν/Δ;
In[]:=
PlotRe
1
Qh
Sqrt
AgAe
ΛΔ
Exp[Ξ]Exp[-
ζ
1
t/2],{t,0,3000},FrameTrue,FrameStyleThickness[0.006],AspectRatio1,BaseStyle{16,FontFamily"TimesRoman"},PlotStyle{Thick,Brown},FrameTicks{True,None},AxesNone,PlotRangeAll
Out[]=
The Log function does not help here
In[]:=
myLog[z_]:=If[Arg[z]<0,Log[z]+I2π,Log[z]]
In[]:=
PlotRe
1
Qh
Exp
1
2
myLog
AgAe
ΛΔ
Exp[Ξ]Exp[-
ζ
1
t/2],{t,0,3000},FrameTrue,AxesNone,FrameStyleThickness[0.006],AxesNone,AspectRatio1,BaseStyle{16,FontFamily"TimesRoman"},PlotStyle{Thick,Red},FrameTicks{True,None},PlotRangeAll
Out[]=