Given a density matrix ρ, its linear purity is Trρ2, its von Neumann entropy is −Trρlogρ. Knowing one how do I calculate the other?
Edit:
My thought is that the linear purity and the von Neumann entropy are a parallel set of benchmarks for mixedness. Their inter-conversion should be like unit conversion.
See it this way, as a state goes from a pure state to a maximally mixed state, Trρ2 goes from 1 to 0, and −Trρlogρ goes from 0 to logd, like two rulers. A point on one should correspond to a point on the other.
Updated answer :
Summary : It is not always that knowing Trρ2 is enough to calculate −Trρlogρ. i.e., −Trρlogρ cannot be determined solely from the knowledge of Trρ2 alone. For a qubit case (2-level system with 2×2 density matrix) it is always possible. For general n-level systems (n>2), more information is required and hence −Trρlogρ is not always an unique function of Trρ2. This will be shown below. Hence the claim made (on intuitive grounds) by OP is only true for a qubit system (nevertheless an useful case), in general it is not feasible. The general reason will be provided towards the end after treating the problem in a hard way for qubit case.
Qubit case (demonstrating possibility) :
Hard way : Consider an arbitary 2×2 density matrix (Positive semi-definite Hermitian matrix with unit trace) : (ρ11ρ12eiϕρ12e−iϕ(1−ρ11))
here unknown real variables
0≤ρ11≤1,
ρ12 and
ϕ.
Consider the explicit expression for Trρ2 : Trρ2=(ρ11)2+(1−ρ11)2+2(ρ12)2,
and an explicit expression for
−Trρlogρ :
−Trρlogρ=−12log(−ρ211+ρ11−(ρ12)2)−√(1−2ρ11)2+4(ρ12)2tanh−1(√(1−2ρ11)2+4(ρ12)2).
After a bit of algebra it can be shown that :
−Trρlogρ=−12log(1−Tr[ρ2]2)−√2Tr[ρ2]−1tanh−1(√2Tr[ρ2]−1).
Hence for the qubit case : −Trρlogρ is an unique function of Trρ2.
Simple way : Since
ρ is a (Positive definite) Hermitian matrix (with unit trace) it can always be diagonalised with an unitary matrix to get (the eigevalues of
ρ are enough to calculate the observables under consideration as they involve
Tr of functions of
ρ alone) :
˜ρ=(p001−p)
where
0≤p≤1. Now calculating the quantities under consideration :
Tr[ρ2]=Tr[˜ρ2]=p2+(1−p)2,
observe that
Tr[ρ2] is a function of single variable, which can be inverted (assuming we know
Tr[ρ2]) to get
p=12±√2Tr[ρ2]−12
out of which we choose only physical branch (satisfying the condition
0≤p≤1 ⇒ 0≤2Tr[ρ2]−1≤1).
Note : The solution switches between the above two around p=12. This when used in (irrespective of the branch) we get
−Tr[ρ.logρ]=−Tr[˜ρ.log˜ρ]=12((√2Tr[ρ2]−1−1)log(1−√2Tr[ρ2]−1)−(√2Tr[ρ2]−1+1)log(√2Tr[ρ2]−1+1)+log(4)).
Again demonstrating the fact that for qubit case OP's proposition is feasible. (Both the results obtained here and above for
Tr[ρ.logρ] are equivalent.)
General n-level system case (impossibility argument) :
For general
n×n case
ρ (again a
n×n positive semi-definite Hermitian matrix with unit trace, can be diagonalized and quantities under consideration can be calculated in eigenbasis), general density matrix in its eigenbasis is given as :
˜ρ=(p1⋱pn),
here
0≤pk≤1 for all
k=1,...,n along with the constraint
∑nk=1pk=1.
Now calculating the quantities of interest (in basis diagonalizing ρ) : Tr[ρ2]=Tr[˜ρ2]=n∑k=1p2k
and
−Tr[ρ.logρ]=−Tr[˜ρ.log˜ρ]=−n∑k=1pklogpk.
Since
Tr[ρ2] (given) imposes only one constraint apart form another constraint
∑nk=1pk=1,
{pk} is an under determined set and cannot be uniquely expressed in terms of
Tr[ρ2] (except for
n=2 qubit case which is concretely demonstrated above explicitly), hence
−Tr[ρ.logρ] cannot be uniquely determined in terms of
Tr[ρ2] alone (more information is needed through imposition of some constraints).
Graphical demonstartion of non uniqueness for 3-level system :
To conclude, OP's proposed claim, although interesting, holds feasible only for the simple
(of couse relevant) case of a qubit. In general n-level (n>2) system case von Neumann
entropy cannot be uniquely determined by linear purity.
No comments:
Post a Comment