I am trying to redraw the function F as a function of Jt.

Where

$$F = {{(1 + 2{{left| {{D_4}left( t right)} right|}^2})} over 3}$$

I should get the same result in the paper “Annals of Physics 355 (2015) 170–181”, Figure 6 (b) solid-black line. But I get different results!

Is there something wrong?

NB:

D4(t) is defined by equation (12d) and F is defined at the bottom of page 7

```
ClearAll("Global`*")
xx = 0.1;
J = 1.0;
lam = J/xx;
t = tau/J;
k = Sqrt(lam^2 + J^2);
D4 = ((J^22)/(4*(k^2 + lam*k)))*
Exp(-1 I*(lam + k)*t) + ((J^2)/(4*(k^2 - lam*k)))*
Exp(-1 I*(lam - k)*t) - (1/2)*Cos(J*t);
F = (1 + 2*Abs(D4)^2)/3.0;
Plot(F, {tau, 0, 50}, PlotRange -> {0, 1})
```