Wolfram Cloud Document

Remove["Global`*"]

In[]:=

cp=310^10;

ω

1

=25×2Picp10^-15;

ω

2

=60×2Picp10^-15;γ1=20×2Picp10^-15;z=3×2Picp10^-15;ℏ=6.62610^-3410^15/(2π);θ=1.3806610^-23;d=2.5;β=

1

θT

;T=2;

In[]:=

A=(

2

ω

1

-

2

ω

2

)(2

2

ω

2

);r=

ω

1

/

ω

2

;η=-I

ω

2

(A+1)t;ϵ=

ω

2

(-βℏ+It);φ=-I

ω

2

At;Q=

-βℏ

ω

2

/2



1-

-βℏ

ω

2



;σ=

2

η

-

2

φ

;x=η-φ;y=η+φ;ϰ1=

kSinh[

xy

]

xy

-

k-1+Cosh[

x

y

]

y

;ϰ2=

kSinh[

xy

]

xy

+

k-1+Cosh[

x

y

]

y

;ϰ3=-

1

4

5/2

(xy)

2

k

x3x

xy

+y

xy

+4

3/2

(xy)

-4x

xy

Cosh[

xy

]+(x-y)

xy

Cosh[2

xy

]-4xySinh[

xy

];v=



2

ω

1

-

2

ω

2

(1-

-2t

ω

1



)

-2t

ω

1



2

(

ω

1

-

ω

2

)

-

2

(

ω

1

+

ω

2

)

;m=v;k=-I

d

2

ω

1

2

ω

2

t;g=Log

2

ω

1

ω

2

2Cos[t

ω

1

]

ω

1

ω

2

+Sin[t

ω

1

]

2

ω

1

+

2

ω

2



;

In[]:=

lin=

ϵ





2g+ϵ



ϰ1(vϰ1-2ϰ2)+2

g



ϰ1ϰ2+

ϵ



m

2

ϰ2



2-4

g+ϵ



+2

2ϵ





2g



-mv



π

1-(2

g+ϵ



-

2(g+ϵ)



+

2ϵ



mv)

;corrfunc=

1

πQ

Expϰ3+

g

2

+

ϵ

2

-

I

ω

2

(rd)

2

t-3zt2lin//N;

In[]:=

Plot[corrfunc//Re,{t,0,3000},PlotRangeAll]

Out[]=

In[]:=

corrfuncII=

1

πQ

Expϰ3+

g

1

+

ϵ

2

-

I

ω

2

(rd)

2

t-3zt2lin//N;

Session II

In[]:=

Plot[g//Re,{t,0,3000},PlotRangeAll]

Out[]=

In[]:=

Plot[

g

//Re,{t,0,3000},PlotRangeAll]

Out[]=