s0=
;s1=
;a=1/2;p=500;LaunchKernels[];v[i_]:=Transpose[{UnitVector[p+2,i]}]ρ[r_]:=SumKroneckerProduct[s1,s0,s0,v[n+1]].Transpose[KroneckerProduct[s1,s0,s0,v[n+1]]],{n,0,p}+SumKroneckerProduct[s1,s0,s0,v[n+1]].Transpose[KroneckerProduct[s0,s1,s0,v[n+1]]],{n,0,p}+SumKroneckerProduct[s1,s0,s0,v[n+1]].Transpose[KroneckerProduct[s0,s0,s1,v[n+1]]],{n,0,p}+SumKroneckerProduct[s0,s1,s0,v[n+1]].Transpose[KroneckerProduct[s1,s0,s0,v[n+1]]],{n,0,p}+SumKroneckerProduct[s0,s1,s0,v[n+1]].Transpose[KroneckerProduct[s0,s1,s0,v[n+1]]],{n,0,p}+SumKroneckerProduct[s0,s1,s0,v[n+1]].Transpose[KroneckerProduct[s0,s0,s1,v[n+1]]],{n,0,p}+SumKroneckerProduct[s0,s0,s1,v[n+1]].Transpose[KroneckerProduct[s1,s0,s0,v[n+1]]],{n,0,p}+SumKroneckerProduct[s0,s0,s1,v[n+1]].Transpose[KroneckerProduct[s0,s1,s0,v[n+1]]],{n,0,p}+SumKroneckerProduct[s0,s0,s1,v[n+1]].Transpose[KroneckerProduct[s0,s0,s1,v[n+1]]],{n,0,p}+Sum(n+1)KroneckerProduct[s0,s0,s0,v[n+2]].Transpose[KroneckerProduct[s0,s0,s0,v[n+2]]],{n,0,p}+Sum=ρ[r];=/Tr[];=Chop[Eigenvalues[]];=0;For[i=1,i<Length[ρ[r]]+1,i++,If[Chop[〚i,i〛]==0,=,=-〚i,i〛Log2[〚i,i〛]]];=0;For[j=1,j<Length[ρ[r]]+1,j++,If[〚j〛==0,=,=-〚j〛Log2[〚j〛]]];,{r,0.01,3.5,0.01}];ListLinePlot[Table[(-Re[])/4,{r,0.01,3.5,0.01}],PlotRange{0,1},RotateLabelTrue,FrameTrue,FrameStyleDirective[Black,12,Bold],FrameLabel{Style["r",13,Bold,Black],Style["",13,Bold,Black]},DataRange{0,3.5},PlotRangePadding0,ImageSizeMedium]
1 |
0 |
0 |
1 |
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
2
Cosh[r]
2
a
2n
Tanh[r]
4
Cosh[r]
2
a
2n
Tanh[r]
3
Cosh[r]
(n+1)
KroneckerProduct[s1,s0,s0,v[n+1]].Transpose[KroneckerProduct[s0,s0,s0,v[n+2]]],{n,0,p}+Sum2
a
2n
Tanh[r]
3
Cosh[r]
(n+1)
KroneckerProduct[s0,s1,s0,v[n+1]].Transpose[KroneckerProduct[s0,s0,s0,v[n+2]]],{n,0,p}+Sum2
a
2n
Tanh[r]
3
Cosh[r]
(n+1)
KroneckerProduct[s0,s0,s1,v[n+1]].Transpose[KroneckerProduct[s0,s0,s0,v[n+2]]],{n,0,p}+Sum2
a
2n
Tanh[r]
3
Cosh[r]
(n+1)
KroneckerProduct[s0,s0,s0,v[n+2]].Transpose[KroneckerProduct[s1,s0,s0,v[n+1]]],{n,0,p}+Sum2
a
2n
Tanh[r]
3
Cosh[r]
(n+1)
KroneckerProduct[s0,s0,s0,v[n+2]].Transpose[KroneckerProduct[s0,s1,s0,v[n+1]]],{n,0,p}+Sum2
a
2n
Tanh[r]
3
Cosh[r]
(n+1)
KroneckerProduct[s0,s0,s0,v[n+2]].Transpose[KroneckerProduct[s0,s0,s1,v[n+1]]],{n,0,p}ParallelTable[pf
r
pf
r
pf
r
pf
r
ef
r
pf
r
fdiag
r
pf
r
fdiag
r
fdiag
r
fdiag
r
fdiag
r
pf
r
pf
r
fSp
r
ef
r
fSp
r
fSp
r
fSp
r
fSp
r
ef
r
ef
r
fdiag
r
fSp
r
