Clear["`*"];​​Do​​​​
ρφAB
r,α
=
2n
Tanh[r]
2
Cosh[r]
2
α
0
0
2n
Tanh[r]
2
Cosh[r]
α
1-
2
α
n+1
Cosh[r]
0
2n
Tanh[r]
2
Cosh[r]
n(1-
2
α
)
2
Cosh[r]
2
Tanh[r]
0
0
0
0
2n
Tanh[r]
2
Cosh[r]
2
α
2
Tanh[r]
0
2n
Tanh[r]
2
Cosh[r]
α
1-
2
α
n+1
Cosh[r]
0
0
2n
Tanh[r]
2
Cosh[r]
(n+1)(1-
2
α
)
2
Cosh[r]
;​​
TρφAB
r,α
=
ρφAB
r,α
[[1,1]]
ρφAB
r,α
[[2,1]]
ρφAB
r,α
[[1,3]]
ρφAB
r,α
[[2,3]]
ρφAB
r,α
[[1,2]]
ρφAB
r,α
[[2,2]]
ρφAB
r,α
[[1,4]]
ρφAB
r,α
[[2,4]]
ρφAB
r,α
[[3,1]]
ρφAB
r,α
[[4,1]]
ρφAB
r,α
[[3,3]]
ρφAB
r,α
[[4,3]]
ρφAB
r,α
[[3,2]]
ρφAB
r,α
[[4,2]]
ρφAB
r,α
[[3,4]]
ρφAB
r,α
[[4,4]]
;​​
TρφAB
r,α
=
TρφAB
r,α
Tr[
TρφAB
r,α
];​​​​
NeABC
rb,rc
=0.5*(Total[Abs[Eigenvalues[
TρφAB
r,α
]]]-1);;​​​​,{r,0,3,0.01},{α,0,1,0.01}​​​​ListPlot3D[Table[
Ne
r,α
,{r,0,3,0.01},{α,0,1,0.01}],AxesLabel{Text["α"],Text["r"],Text["Ne"]},DataRange{{0,1},{0,3}},PlotRange{0,0.6}]
Power
:Infinite expression
1
0.
encountered.
Infinity
:Indeterminate expression 0ComplexInfinity encountered.
Eigenvalues
:Input matrix contains an indeterminate entry.
Sum
:Sum does not converge.
Power
:Infinite expression
1
0.
encountered.
Infinity
:Indeterminate expression 0ComplexInfinity encountered.
Eigenvalues
:Input matrix contains an indeterminate entry.
Sum
:Sum does not converge.
Power
:Infinite expression
1
0.
encountered.
General
:Further output of Power::infy will be suppressed during this calculation.
Infinity
:Indeterminate expression 0ComplexInfinity encountered.
General
:Further output of Infinity::indet will be suppressed during this calculation.
Eigenvalues
:Input matrix contains an indeterminate entry.
General
:Further output of Eigenvalues::mindet will be suppressed during this calculation.
Sum
:Sum does not converge.
General
:Further output of Sum::div will be suppressed during this calculation.
​
​