Friday, 8 December 2017

quantum mechanics - Linear purity and von Neumann entropy


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.




Answer



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ϕρ12eiϕ(1ρ11))

here unknown real variables 0ρ111, ρ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)(12ρ11)2+4(ρ12)2tanh1((12ρ11)2+4(ρ12)2).
After a bit of algebra it can be shown that : Trρlogρ=12log(1Tr[ρ2]2)2Tr[ρ2]1tanh1(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) : ˜ρ=(p001p)
where 0p1. Now calculating the quantities under consideration : Tr[ρ2]=Tr[˜ρ2]=p2+(1p)2,
observe that Tr[ρ2] is a function of single variable, which can be inverted (assuming we knowTr[ρ2]) to get p=12±2Tr[ρ2]12
out of which we choose only physical branch (satisfying the condition 0p1 02Tr[ρ2]11). 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]11)log(12Tr[ρ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 : ˜ρ=(p1pn),
here 0pk1 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]=nk=1p2k

and Tr[ρ.logρ]=Tr[˜ρ.log˜ρ]=nk=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

Understanding Stagnation point in pitot fluid

What is stagnation point in fluid mechanics. At the open end of the pitot tube the velocity of the fluid becomes zero.But that should result...