What is the entropy of a density matrix which is the sum of two unitarily equivalent projectors? - MathOverflow most recent 30 from http://mathoverflow.net 2013-06-20T09:23:56Z http://mathoverflow.net/feeds/question/28525 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/28525/what-is-the-entropy-of-a-density-matrix-which-is-the-sum-of-two-unitarily-equival What is the entropy of a density matrix which is the sum of two unitarily equivalent projectors? Jess Riedel 2010-06-17T16:06:34Z 2010-06-21T19:38:45Z <h2>Construction</h2> <p>Suppose I have a density matrix $\rho$ which is proportional to a projector $P$ formed by tensoring together $N$ small projectors $P^{(i)}$ of rank 2:</p> <p>$P^{(i)} = |a\rangle_i\langle a| + |b\rangle_i\langle b|$</p> <p>$P=P^{(1)}\otimes \cdots \otimes P^{(N)} = \bigotimes_{i=1}^N P^{(i)}$</p> <p>$\rho = P/2^n \quad , \quad \mathrm{Tr} \rho = 1 $.</p> <p>(I'm a physicists, and I don't really know how to express $P^{(i)}$ explicitly without using bra-ket notation.) The small projectors $P^{(i)}$ act on their own (potentially very large) Hilbert space $H^{(i)}$ spanned by, say, $\{|a\rangle_i, |b\rangle_i, |c\rangle_i,...\}$, and $P$ acts on $H = \bigotimes_{i=1}^N H^{(i)}$.</p> <p>Then since the $2^N$ non-zero eigenvalues of $\rho$ are all the same, the entropy is obviously</p> <p>$H[\rho] = -\sum_{s} \lambda_s \ln \lambda_s = N \ln 2 $.</p> <p>Now, suppose I split the density matrix in half, and operate on each half with either a <em>local</em> set of unitaries (all identical except that the act on different $H^{(i)}$) or their inverses. That is, suppose</p> <p>$U = \bigotimes_{i=1}^N U^{(i)}$</p> <p>$U^{(i)} |a\rangle_i = |a+\rangle_i \quad , \quad {U^{(i)}}^{-1} |a\rangle_i = |a-\rangle_i \quad , \quad U^{(i)} |b\rangle_i = |b+\rangle_i \quad , \quad {U^{(i)}}^{-1} |b\rangle_i = |b-\rangle_i $</p> <p>with a new density matrix $\eta$:</p> <p>$\eta = (U P U^{-1} + U^{-1} P U)/2^{N+1}$.</p> <p>Assume we know the inner products</p> <p>$\langle a-|b+\rangle \quad, \quad \langle a+|b-\rangle \quad, \quad \langle a+|a-\rangle \quad, \quad \langle b+|b-\rangle \quad ,$</p> <p>(which are the same for all $i$ because the $U^{(i)}$ are all identical). </p> <h2>Question</h2> <p>Is it possible to calculate the entropy $H[\eta]$? </p> <p>If useful, it's fine to assume that $U^{(i)}$ is infitesimal so that </p> <p>$\langle a-|b+\rangle \sim 0 \quad, \quad \langle a+|b-\rangle \sim 0 \quad, \quad \langle a+|a-\rangle \sim 1 \quad, \quad \langle b+|b-\rangle \sim 1 \quad .$</p> <h2>Simpler case</h2> <p>It might be helpful to see the solution to the (much) simpler case where the $P^{(i)}$ are rank 1. In this case, we can write</p> <p>$P^{(i)} = |a\rangle_i\langle a|$</p> <p>$P = |a\rangle_1\langle a| \otimes \cdots \otimes |a\rangle_N\langle a| = |A\rangle\langle A|$</p> <p>$\rho = P \quad , \quad \mathrm{Tr} \rho = 1$</p> <p>$\eta = (U P U^{-1} + U^{-1} P U)/2 $</p> <p>$\quad = (|A+\rangle\langle A+| + |A-\rangle\langle A-|)/2$</p> <p>and get that the two non-zero eigenvalues are</p> <p>$\lambda_\pm = [1 \pm |\langle A+|A-\rangle_i|]/2$</p> <p>where</p> <p>$\langle A+|A-\rangle = \prod_{i=1}^N \langle a+|a-\rangle_i = \langle a+|a-\rangle^N$.</p> <h2>Elaboration of Answer</h2> <p>Steve Flammia graciously provided the answer below, which hinges on the relationship between the eigenvalues of the two operators $P+Q$ and $PQP$, where $P$ and $Q$ are orthogonal projectors. (The eigenvalues of $PQP$ are easy to compute because the matrix is diagonal in the same basis as $P$ and the matrix elements are expressed in terms of the known inner products.) As an exercise for myself, and to help out anyone else reading this question, the relationship is proved here.</p> <p>Consider two orthogonal projectors $P$ and $Q$ with respective images $M$ and $N$. We can assume the intersections</p> <p>$ \quad M \cap N \quad , \quad M \cap N^\perp \quad , \quad M^\perp \cap N \quad , \quad M^\perp \cap N^\perp$</p> <p>are all {0}, i.e. $M$ and $N$ are in "generic position". (If not, then the Hilbert space can be decomposed into $H=H' \oplus H_0$, where $M$ and $N$ leave both $H'$ and $H_0$ invariant, $M$ and $N$ are simultaneously diagonal on $H_0$, and the above intersections are trivial on $H'$. We then restrict our attention $H'$.) As shown by Halmos [1], for any pair of orthogonal projectors $P$ and $Q$ of equal rank we can decompose the Hilbert space into two subspaces of equal dimension, $H=K \oplus K$, in which</p> <p>$P = \pmatrix{1&amp; 0 \\ 0&amp; 0}$</p> <p>$Q = \pmatrix{C^2 &amp; CS \\ CS &amp; S^2}$</p> <p>where $C$ and $S$ are positive matrices on $K$ which commute and for which $C^2 + S^2 =1$. If we let $c_i$ be the eigenvalues of $C$, then the non-zero eigenvalues of </p> <p>$PQP = \pmatrix{C^2&amp; 0 \\ 0&amp; 0}$</p> <p>are $c_i^2$ and the eigenvalues of</p> <p>$P+Q=\pmatrix{1 + C^2&amp; CS \\ CS &amp; S^2}$</p> <p>can be shown to be $1 \pm c_i$.</p> <h2>References</h2> <p>[1] Halmos, P.R. <a href="http://www.jstor.org/stable/1995288" rel="nofollow">"Two Subspaces"</a>. Transactions of the American Mathematical Society, Vol. 144 (Oct., 1969), pp. 381-389. </p> http://mathoverflow.net/questions/28525/what-is-the-entropy-of-a-density-matrix-which-is-the-sum-of-two-unitarily-equival/28537#28537 Answer by Steve Flammia for What is the entropy of a density matrix which is the sum of two unitarily equivalent projectors? Steve Flammia 2010-06-17T17:13:46Z 2010-06-17T17:13:46Z <p>First observation, you can unitarily rotate the first term without changing the spectrum so you can just consider a single unitary with out loss of generality. Call the rotated term $Q$.</p> <p>Next, it helps to know about the canonical form for two projectors. Given two projectors P and Q in general position, you can always compute the eigenvalues of P+Q as follows. They are simply $1 \pm x_j$ where the $x_j$ are the positive square roots of the eigenvalues of the operator $PQP$ (or $QPQ$... it doesn't matter). </p> <p>For your problem, because of the tensor product structure, this is actually quite easy to do. The eigenvalues $x_j$ will be all possible products of the local eigenvalues for $P^{(i)} Q^{(i)} P^{(i)}$, with a square root. The entropy is then a sum over these configurations, suitably normalized. </p>