Entanglement of extremal density matrices of 2-qubit Hamiltonian with Kramers degeneracy arxiv: v1 [quant-ph] 30 Jun 2018

Transcription

1 Entanglement of extremal density matrices of 2-qubit Hamiltonian with Kramers degeneracy arxiv: v [quant-ph] 30 Jun 208 A. Figueroa, O. Castaños, R. López-Peña Instituto de Ciencias Nucleares Universidad Nacional Autónoma de México Apdo. Postal México 0450 D.F. Abstract We establish a novel procedure to analyze the entanglement properties of extremal density matrices depending on the parameters of a finite dimensional Hamiltonian. It was applied to a general 2-qubit Hamiltonian which could exhibit Kramers degeneracy. This is done through the extremal density matrix formalism, which allows to extend the conventional variational principle to mixed states. By applying the positive partial transpose criterion in terms of the Correlation and Schlienz-Mahler matrices on the extremal density matrices, we demonstrate that it is possible to reach both pure and mixed entangled states, changing properly the parameters of the Hamiltonian. For time-reversal invariant Hamiltonians, the extremal pure states can be entangled or not and we prove that they are not time-reversal invariants. For extremal mixed states we have in general 5 possible cases: three of them are entangled and the other two separable. Keywords: Extremal density matrices; 2-qubit system; Entanglement; Kramers degeneracy.

2 . Introduction Recently we have extended the conventional variational method to density matrices of a qudit system. This was done by means of a Lagrange multipliers approach [], using an algebraic procedure [2] and through the geometric formulation of quantum mechanics [3]. We have shown that, for a finite system the mean value of the Hamiltonian operator Ĥ (or any observable achieves its critical values under the condition [ˆρ, Ĥ] = 0, which it is equivalent to the stationary solution of the von Neumann equation. We call extremal density matrices to the states that fulfill this condition. They provide an extremal description of the mean values of the Hamiltonian, and in the case of restricting them to pure states, one recovers the energy spectrum. So, besides of being an alternative tool to find either the eigensystem or detect level crossings in the Hamiltonian without computing its eigenvalues, one obtains information of mixed states which minimize the mean value of the energy. Among other methods to detect entanglement in bipartite systems [4], it has been shown that if the entire system is separable, the partial transposition operation has the property of preserving the positive definiteness of the density matrix. Hence, the necessary condition for separability of a finite dimensional state is the positivity of its partial transposition, also commonly referred as PPT criterion [5]. For the 2-qubit state and the tensorial product of a qubit times a qutrit, the PPT criterion is also a sufficient one [6]. The present article arises from the concern about the entanglement of extremal density matrices, associated to time-reversal invariant Hamiltonian matrices denoting 2-qubit general systems. A Hamiltonian with this type of time-reversal symmetry can be constructed following the procedure indicated by Haake [7]. As we shall shown, the main novelty is that, by applying the extremal density matrix procedure together with the PPT criterion, one can obtain mixed separable or mixed entangled states by changing properly the parameters of the time-reversal invariant Hamiltonian. Extending this result further, at least for 4 and 6 dimensional Hilbert spaces, it will be possible to obtain mixed separable and mixed entangled states tuning the parameters of the Hamiltonian, whether it has a specific symmetry or not. The paper is organized as follows. We start with a summary of the procedure to determine extremal density matrices of a finite dimensional Hamiltonian. In section 3, the Kramers degeneracy or time-reversal invariance is reviewed and the general form of a 4 4 Hamiltonian with this symmetry is given. The PPT criterion is established in section 4, and its connection with the semi-positivity conditions on the extremal density matrices is shown. This connection is given by the Correlation and Schlienz-Mahler matrices [8, 9]. Finally, in section 5, our method is applied to the general four-dimensional Hamiltonian which exhibits Kramers degeneracy. 2

3 2. Extremal density matrices The space of Hermitian matrices can be stratified by means of the Partial Flag Manifolds F (d; m, m 2,..., m k and the number of different strata is equal to the partition function p(d (see 35. If we denote by r the dimension of F (d; m, m 2,..., m k, then any quotient of unitary groups Û F (d; m, m 2,..., m k, is specified by r real parameters (see Table. Consequently, the dimension of Ĥ, i.e., the number of independent real parameters needed to specify Ĥ, is given by the formula [0] ( k dim(ĥ = d2 m 2 j k, ( j= where k is the number of different eigenvalues of Ĥ and {m j} denote their algebraic multiplicities, such that m +m 2 + +m k = d. The quantity in parentheses is the codimension of Ĥ, codim(ĥ, and represents the number of conditions to be fulfilled for a level crossing []. Table : Dimension, Codimensions and Partial Flag Manifolds F (d; m, m 2,..., m k for the unitary orbits of Ĥ. It is supposed that α > β > γ > δ. If they are given explicitly, the dimension dim(ĥ is decreased by k and it is denoted by r, which is also the dimension of the Manifold. Hilbert space Hermitian Eigenvalue Codimension Dimension Partial Flag Manifold dimension diagonal parameters codim(ĥ dim(ĥ Manifold dimension d representation k k j= m2 j k d2 codim(ĥ F (d; m, m 2,..., m k r = dim(ĥ k diag(α, α 3 Point 0 2 diag(α, β U(2/[U( U(] 2 diag(α, α, α 8 Point 0 3 diag(α, β, β U(3/[U( U(2] 4 diag(α, β, γ U(3/[U( U( U(] 6 diag(α, α, α, α 5 Point 0 diag(α, β, β, β U(4/[U( U(3] 6 4 diag(α, α, β, β U(4/[U(2 U(2] 8 diag(α, β, γ, γ U(4/[U( U( U(2] 0 diag(α, β, γ, δ U(4/[U( U( U( U(] 2 The Rayleigh quotient e H (ψ of a Hermitian matrix Ĥ is e H (ψ := ψ Ĥ ψ ψ ψ where ψ is a d-dimensional complex vector. The numerical range W (Ĥ is the set of all possible Rayleigh quotients e H (ψ over the unit vectors. It is a closed interval on the real, 3

4 axis and the eigensystem of Ĥ is associated to the critical points of e H(ψ. Thus W (Ĥ is the convex hull of the eigenvalues [2]. In the density matrix formalism, the numerical range of the Hamiltonian (or any Hermitian operator can be identified with Ĥ = Tr(Ĥ ˆρ and, by restricting to variations along the unitary orbits, it can be shown that under the condition [ˆρ, Ĥ] = 0 the Rayleigh quotient achieves its critical values at Ĥ c = Tr(Ĥ ˆρc, where ˆρ c denotes an extremal density matrix commuting with Ĥ [2, 3]. It is known that two Hermitian operators with vanishing commutator share a common eigenbasis. The difference between degenerate and non-degenerate cases of Ĥ is that, in the former case, not all its eigenvectors are necessarily eigenvectors of ˆρ c. However, it is generally expected that for a given Hamiltonian Ĥ, the dimension of an unknown ˆρc can be determined up to the same Partial Flag Manifold dimension of Ĥ, i.e., from Table, the number of parameters describing ˆρ c is Then the number of free parameters n of ˆρ c is given by dim(ˆρ c = r. (2 n = (d 2 r. (3 An algebraic proof of these results is found by using [Ĥ, ˆρ] = 0, the properties of the generators of the su(d algebra, the tangent vectors of Ĥ at the identity and the Gram matrix G formed with them [2]. Therefore, r = rank G also determines the dimension of the tangent space of the Manifolds (see Table and is an alternative tool to detect level crossings without computing the eigenvalues of Ĥ. 3. Hamiltonians with Time-Reversal Invariance Unitary transformations leaving invariant a Hamiltonian matrix give rise conserved quantities. Anti-unitary transformations sometimes increase the degree of degeneracy as it is the case for the time-reversal invariance, described by an operator ˆT which satisfies ˆT 2 = Î and [ ˆT, Ĥ] = 0. For half-integer total spin S = (2N /2 (N N, evendimensional Hilbert space d = 2N, time-reversal invariance implies Kramers degeneracy, i.e., pairs of energy levels of Ĥ are degenerated [7]. For any state vectors { φ, ψ } and complex numbers {c, c 2 }, ˆT has the properties [5]: i (fermionic condition ˆT 2 = Î, ii (anti-linearity ˆT ( c ψ + c 2 φ = c ˆT ψ + c 2 ˆT φ, iii (anti-unitarity ˆT ψ ˆT φ = ψ φ. 4

5 From these properties, it follows that ˆT has an inverse ˆT, preserves the norm, ˆT ψ is unique and orthogonal to ψ, and consequently, ˆT has no eigenvectors. Three implications can be deduced [6]: All eigenvalues of Ĥ are doubly degenerate (Kramers degeneracy; the corresponding Hilbert space H d can not be decomposed into invariant subspaces with respect to ˆT ; and a symmetry-adapted orthonormal basis for H d exists and it takes the form { ψ, ˆT ψ, ψ 2, ˆT ψ 2,..., ψ N, ˆT ψ N }. In this basis, for N = 2 (d = 4 case, the general form of the traceless Hamiltonian [7], β 0 γ is ɛ iδ Ĥ = 0 β ɛ iδ γ + iσ γ + is ɛ + iδ β 0, (4 ɛ + iδ γ iσ 0 β presents Kramers degeneracy if s = σ and consequently, it is double degenerate. On the other hand, we consider broken the time-reversal invariance by taking s = σ, thus Ĥ is non degenerate. 4. Separability criterion for two mixed qubits A basis for the 2-dimensional matrices is given by the identity and Pauli matrices, ( ( ( ( i 0 ˆσ 0 =, ˆσ 0 =, ˆσ 0 2 =, ˆσ i 0 3 =, 0 with the product rule ˆσ i ˆσ j = i ɛ ijk ˆσ k + δ ij ˆσ 0, and orthogonality relation Tr(ˆσ j ˆσ k = 2δ jk. An arbitrary Hamiltonian Ĥ and single-qubit state ˆρ can be represented as Ĥ = 3 h k ˆσ k, (5 2 ˆρ = 2 k=0 3 r k ˆσ k, (6 with the identifications h k = Tr(Ĥ ˆσ k, r k = Tr(ˆρ ˆσ k and r 0 = Tr(ˆρ =. k=0 Similarly, in the 2-qubit case Ĥ and ˆρ can be parametrized as Ĥ = 3 3 h p q ˆDp,q, (7 4 ˆρ = 4 p=0 q=0 3 p=0 q=0 3 r p q ˆDp,q, (8 5

6 where ˆD p,q = ˆσ p ˆσ q, h p q = Tr(Ĥ ˆD p,q, r p q = Tr(ˆρ ˆD p,q and r 00 = Tr(ˆρ =. In matrix form, ˆρ can be written in terms of 2 2 matrices, ( F G ˆρ =, (9 P Q where P = G and F = ( r03 + r 30 + r 33 + r 0 + r 3 i(r 02 + r 32 4 r 0 + r 3 + i(r 02 + r 32 r 03 + r 30 r 33 + G = ( r0 + r 3 i(r 20 + r 23 r r 22 i(r 2 + r 2 4 r + r 22 + i(r 2 r 2 r 0 r 3 i(r 20 r 23 Q = ( r03 r 30 r 33 + r 0 r 3 i(r 02 r 32 4 r 0 r 3 + i(r 02 r 32 r 03 r 30 + r 33 +,,. The basis matrices { ˆD p,q = ˆσ p ˆσ q } fulfill the orthogonality condition and they comply with the multiplication rules Tr( ˆD j,k ˆDm,n = 4δ jm δ kn, (0 ˆD p,0 ˆD0,q = ˆD 0,q ˆDp,0 = ˆD p,q, ˆD i,0 ˆDj,0 = i ɛ ijk ˆDk,0 + δ ij ˆD0,0, ( ˆD 0,i ˆD0,j = i ɛ ijk ˆD0,k + δ ij ˆD0,0. By means of the previous expressions, the following commutators can be calculated: [ ˆD j,0, ˆDm,n ] = 2i ɛ jnq ˆDm,q, [ ˆD 0,j, ˆDm,n ] = 2i ɛ jmq ˆDq,n, (2 [ ˆD i,j, ˆDp,q ] = 2i (ɛ jql δ ip ˆDl,0 + ɛ ipk δ qj ˆD0,k. The parameterisation (8 is an abbreviation of the Fano form [8, 7], ˆρ = Î 4 + r p 0 ˆDp,0 + r 0 q ˆD0,q + C p q ˆDp,q, (3 4 p= q= p= q= where the two Bloch vectors τ A = (r 0, r 20, r 30, τ B = (r 0, r 02, r 03, (4 6

7 determine the properties of the individual qubits A and B. For s, t =, 2, 3, the matrix C s t = r s t is given by r r 2 r 3 C = r 2 r 22 r 23, (5 r 3 r 32 r 33 from which it is possible to build M s t = C s t r s 0 r 0 t, C r 0 r 0 C 2 r 0 r 02 C 3 r 0 r 03 M = C 2 r 20 r 0 C 22 r 20 r 02 C 23 r 20 r 03. (6 C 3 r 30 r 0 C 32 r 30 r 02 C 33 r 30 r 03 They are known as the Correlation matrix C and the Schlienz-Mahler matrix M, respectively [8], and together they describe the correlations between both subsystems. If C = 0 or M = 0 then the state is separable, while for pure states a good measure of entanglement is β = 4 5 Tr(M T M, (7 where 0 β and, from here on, the upper index T denotes transposition of matrices. The value β = corresponds to the maximal entangled pure state with τ A = τ B = 0, whereas β = 0 to a separable one [9]. Also, β can be interpreted geometrically as a distance between entangled and separable states [8]. Additionally, it was noted in [9, 20] that det C and det M can be written in terms of SU(2 SU(2 polynomial invariants of third and fourth degrees, respectively. In order to extend the applications of the extremal mixed density matrices, we consider the PPT criterion in the four dimensional case given in terms of the positivity conditions of the density matrix together with the Correlation and Schlienz-Mahler matrices. Taking into account (9, the partial transpose of ˆρ with respect to the A subsystem, denoted as ˆρ P T A, and similarly for B, are given by [7] ( ( ˆρ P T A F P =, ˆρ P T B F T G = T G Q P T Q T. (8 The PPT criterion is a necessary condition, for the joint density matrix of the 2-qubits A and B subsystems to be separable. In the 4 and 6 dimensional cases the condition is also sufficient [5], i.e., if the partial transposed density matrix ˆρ P T A (ˆρ P T B is positive definite, the state ˆρ is separable, or else, the state ˆρ is entangled if its partial transposition is not positive definite. Notice that ˆρ P T A and ˆρ P T B share the same characteristic polynomial and 7

8 then the indices A and B can be omitted. This entails that, in the 2-qubit case, the PPT criterion and the C and M matrices are related through the coefficients {a P 2 T, a P 3 T, a P 4 T } of the characteristic polynomial of the 4 4 partial transposed density matrix ˆρ P T from ˆρ as [20] 0 a P T 2 = a 2 3 8, 0 a P T 3 = a det C 6, (9 0 a P T 4 = a det M 256, such that any separable state represented by ˆρ must fulfill those inequalities. Conversely, if an inequality is violated then the state ˆρ is entangled. Hence, given a 4 4 extremal density matrix, one can know immediately whether the states that they represent are separable or not. 5. Entaglement for a 2-qubit extremal density matrices In order to apply the discussed separability criteria together with the extremal density matrix approach, in this section we will consider the general parametrization (3 and two cases in the Hamiltonian (4: The non-degenerate case when s = σ, and the Kramers degeneracy case when s = σ. 5.. Non-degenerate case (s = σ In this case, the matrix Hamiltonian does not commute with the time-reversal antiunitary operator ˆT. Under the condition [ˆρ, Ĥ] = 0, the density matrix commuting with Ĥ takes the following values of r pq r 20 = r ( δ 2 + ɛ 2 δ σ r 0, r 30 = β r 0 γδ γ, r 0 = δ r 02, r 2 = ɛ r ɛ δ, r 2 = δ r 0 σ r, γ r 22 = ɛ (δ r 0 σ r, r 3 = β r γδ γ, r 32 = β ɛ r, r 03 = r 3 = r 23 = r 33 = 0, (20 γδ with 3 free variables r 0, r 02, r. Hence, from Table and Eq. (3, Ĥ is non degenerate. These free variables are determined by establishing the system of polynomial equations (39, where the constants c 2, c 3 and c 4 must lie inside the allowed region exhibited in Fig.. 8

9 For the pure case, associated to c 2 = c 3 = c 4 = 0, the set of solutions for this polynomial systems are denoted as {ˆρ c ±, ˆρc 2± }, whose respective Bloch vectors, Correlation and Schlienz-Mahler matrices are τ± A = (δ, ɛ, 0, τ B δ 2 + ɛ 2 ± = ± ( γ, ( δ E 2 + ɛ 2 + σ, β, + ( γδ γɛ 0 δ ( δ C ± = ± E + δ 2 + ɛ 2 δ 2 + ɛ 2 + σ ɛ 2 + ɛ 2 + σ 0, βδ βɛ 0 τ2± A = (δ, ɛ, 0, τ B δ 2 + ɛ 2 2± = ± ( γ, ( δ E 2 + ɛ 2 σ, β, (2 ( γδ γɛ 0 δ ( δ C 2± = ± E δ 2 + ɛ 2 δ 2 + ɛ 2 σ ɛ 2 + ɛ 2 σ 0, βδ βɛ 0 M ± = M 2± = 0, where E ± = β 2 + γ 2 + δ 2 + σ 2 + ɛ 2 ± 2σ δ 2 + ɛ 2. The respective expectation values of the Hamiltonian are Tr(ˆρ c ± Ĥ = ±E +, Tr(ˆρ c 2± Ĥ = ±E, and it is possible to corroborate that {ˆρ c ±, ˆρc 2± } form a complete set of orthogonal states. They are rank one projectors and ˆρ c + + ˆρc + ˆρc 2+ + ˆρc 2 = Î4. All of them are separable because M ± = M 2± = 0. Due to the convex property of the separable states [2] and the unitarily evolution, it follows that Ĥ does not generate entangled mixed states. In order to study the behavior of the expectation values of the Hamiltonian, we consider the parameters β = γ = σ = ɛ = for two cases: (i the pure case when one has c 2 = c 3 = c 4 = 0, which has four independent solutions for the variables r 0, r 20, r. The energy spectra is a function of the parameter δ and its energy levels are plotted in Fig 2(a. (ii The mixed case is established by taking from the region exhibited in Fig. (c the values c 2 = 59/200, c 3 = 9/400, c 4 = 8/ One has 6 extremal expectation values of the Hamiltonian, two of them correspond to Ĥ c = 0. The results are shown in Fig. 2(b with dotted lines. Notice that the extremal expectation values are contained within the minimum and maximum eigenvalues of the Hamiltonian, as it should be. 9

10 Figure : Solid figure which represents the region of c2, c3, and c4 where the positivity conditions of density matrix are satisfied (see A. The pure case is associated to (c2, c3, c4 = (0, 0, 0 while the maximal mixed state corresponds to (c2, c3, c4 = (3/8, /6, /256. <H>c 0 <H>c d -0-5 (a d -5-0 (b -0 Figure 2: hh ic as a function of δ with β = γ = σ = =. (a Pure case with c2 = c3 = c4 = 0; and (b Mixed case with c2 = 59/200, c3 = 9/400, c4 = 8/ Black continuous lines represent mean values of H in the pure case. The dotted ones correspond to the mixed case. 0

11 5.2. Degenerate case (s = σ Under the condition [ˆρ, Ĥ] = 0, the density matrix commuting with following values for the parameters r pq : Ĥ takes the r 3 = σ r 02 + ɛ r 03 + β r γ, r 32 = σ r 0 δ r 03 + β r 2, r 33 = ɛ r 0 + δ r 02 + β r 3, γ γ r 20 = δ r + ɛ r 2 + σ r 3, r 30 = β r 0 γ γ, r 2 = δ r 0 γ, r 22 = ɛ r 0 γ, r 23 = σ r 0. (22 γ Notice that only 8 parameters were determined and thus one has 7 free variables r 0, r 0, r 02, r 03, r, r 2, r 3. Hence, from Table and Eq. (3, it implies that Ĥ is double degenerate Pure state solution For the pure state case, one has to solve the system of equations ˆρ 2 = ˆρ. Finally one has 2 free parameters r 02 and r 03. One possibility is to solve r 02 = r 03 = 0 and so one gets four solutions for the extremal density matrices denoted by {ˆρ c ±, ˆρc 2± }, whose respective Bloch vectors, Correlation and Schlienz-Mahler matrices are ( γ τ± A = ± E, δ E 2 + ω, β (, τ B 2 ± = E 2 + ω, 0, 0, 2 γ( 2 + ω 2 γδɛ ± βeσ γδσ βeɛ C ± = ± E δ 2 + ω 2 ɛ 2 + ω 2 σ 2 + ω 2, 2 + ω 2 β( 2 + ω 2 βδɛ γeσ βδσ ± γeɛ γ δ ω, (23 2 M ± = C ± ± E 2 + ω 2 ( γ τ2± A = ± E, δ E C 2± = ± E 2 + ω ω, β, τ 2 2± B = E M 2± = C 2± ± E 2 + ω 2 β ( 2 + ω, 0, 0, 2 γ( 2 + ω 2 (γδɛ ± βeσ (γδσ βeɛ δ 2 + ω 2 ɛ 2 + ω 2 σ 2 + ω 2 β( 2 + ω 2 (βδɛ γeσ (βδσ ± γeɛ γ δ 3 2 +ω β 2 0 0,,

12 where we define = β 2 + γ 2, ω = σ 2 + ɛ 2, and E = (det Ĥ 4 = 2 + ω 2 + δ 2. The respective expectation values of the Hamiltonian are Tr(ˆρ c ± Ĥ = ±E, Tr(ˆρc 2± Ĥ = ±E, and it is possible to corroborate that the set {ˆρ c ±, ˆρc 2± } constitutes a complete set of orthogonal rank one projectors because ˆρ c + + ˆρc + ˆρc 2+ + ˆρc 2 = Î4. From (7, the parameter β for each state of {ˆρ c ±, ˆρc 2± } take the same value, β = 6 5 S L( + S L, S L = ω 2 2( 2 + ω 2, (24 indicating entanglement between the qubits for a collection of finite values of {β, γ, σ, ɛ}. The defined S L is precisely the Linear Entropy, S L = ( τ A± 2 = ( τ B± 2. ( One can observe that this quantity is independent of the parameter δ and there is no entanglement if σ = ɛ = Mixed state solution The general extremal mixed state of the Hamiltonian can be written in the form ˆρ mix = P ˆρ c + + P 2 ˆρ c 2+ + P 3 ˆρ c + P 4 ˆρ c 2, (26 where 4 j P j = and 0 P j. By means of the Fano representation (3 and the solutions (23, it is straightforward that ˆρ mix = 3 Î 4 + τp A 4 ˆD p,0 + τq B ˆD 0,q + C p q ˆDp,q, (27 p= q= p= q= 2

13 where τ A = C = ( x γ E, δ z E 2 + ω, x β (, τ B = y 2 E 2 + ω, 0, 0, 2 z γ( 2 + ω 2 z γδɛ + y βeσ z γδσ y βeɛ E x δ 2 + ω 2 x ɛ 2 + ω 2 x σ 2 + ω 2,( ω 2 z β( 2 + ω 2 z βδɛ y γeσ z βδσ + y γeɛ x y γ ω 2 M = C E 2 + ω 2 y z δ 3 2 +ω x y β with the following definitions x = P + P 2 P 3 P 4, y = P P 2 + P 3 P 4 and z = P P 2 P 3 + P 4. The respective determinants of C and M are det C = ω2 2 + ω 2 x y z, (29 ω 2 ( det M = E 2 ( 2 + ω ( 2 + ω 2 x 2 y 2 + δ 2 2 y 2 z 2 E 2 ( 2 + ω 2 xyz. (30 Now we are going to use the PPT criterion by considering the expression (9. Thus in Table 2, for any extremal mixed state (27, the set of values (a P 2 T, a P 3 T, a P 4 T from Eqs. (9 are shown. Therefore we have five possibilities for the extremal mixed density matrices. Three of the cases determine entangled extremal density matrices, which are given for three equal probabilities, i.e., P 2 = P 3 = P 4, two equal probabilities P 3 = P 4, and all probabilities different. The separable cases occur for the maximal mixed state and when one has two equal probabilities. For all the cases we have similar results permutating the values of the probabilities. The PPT criterion can be also applied to the pure case taking P 2 = P 3 = P 4 = 0. For this we have det C = ω 2 /( 2 + ω 2 and det M = ω 2 det C, in agreement with the discussion of the previous subsection., States with Kramers invariance From the 2-qubit degenerate case Hamiltonian studied above, we have found that pure extremal states do not commute with the time-reversal operator ˆT. Due to this observation, we are going to prove that, in general, this occurs for any 2N dimensional Hilbert space. Proposition. If [ ˆT, Ĥ] = 0, extremal pure states do not commute with ˆT. 3

14 Table 2: Values {a P 2 T, a P 3 T, a P 4 T } from Eqs. (9 in terms of the probability coefficients {P, P 2, P 3, P 4 } of the general 2-qubit mixed extremal density matrix (26. Five cases are taken into account according to the different strata shown in the Table for d = 4. (P, P2, P3, P4 a P 2 T a P 3 T a P 4 T (/4, /4, /4, /4 3/8 /6 /256 ( 3 b, b, b, b 3( 2bb (3 8bb 2 ω 2 4( 2 +ω 2 ( 4b3 ( 3bb 3 ω + 2 ( 4b (( 3 4bδ 2 2 ( 2 + ω 2 ((4b 2 + E 2 6 E 2 ( 2 +ω 2 2 ( (/2 b, /2 b, b, b 4 + b 2b2 2 ( 2bb ( b 2 2 b 2 ( ( b 2 c, b, c, c b 2 2bc+ c (c 4bc 2(b b 2c 2 bc 2 ω ( b 2c + 2 (2b + 2c (δ (2b + 2c E 2 ( 2 +ω 2 2 ω b + c(2 3c 2 4( 2 +ω 2 ( 4c(2b + 2c 2 (4c ( 2 + ω 2 ((4c 2 + E 2 ( b c d, b, c, d b 2 d(b + c bc+ b( c d(c + d + cd( c d b 2 (c + d bcd( b c d + ω 2 6 E 2 ( 2 +ω 2 2 ( 2 ( 2 + ω 2 ( 2c 2d 2 ( 2b 2d 2 + b c 2 + c d 2 ω + d 2 4( 2 +ω 2 ( 2b 2c( 2b 2d( 2c 2d δ2 2 ( 2b 2d 2 ( 2b 2c 2 E 2 ( 2 + ω 2 ( 2c 2d( 2b 2d( 2b 2c Proof. Consider that [ ˆT, Ĥ] = 0. Consequently, for k =,..., N, one can construct an orthonormal basis consisting on eigenvectors of Ĥ as { ψ, ˆT ψ, ψ 2, ˆT ψ 2,..., ψ N, ˆT ψ N }, such that the pure states ˆρ k = ˆρ 2 k = ψ k ψ k, ˆη k = ˆη k 2 = ˆT ψ k ψ k ˆT commuting with Ĥ are rank-one orthogonal projectors (Kramers pairs which describe a two-dimensional degenerate space of Ĥ. Suppose now that [ ˆT, ˆρ k ] = 0. Then, it follows that Tr(ˆρ k ˆη k = Tr(ˆρ k ˆρ k =, contrary to the orthogonality of the basis. Because all pure states are unique up to a unitary transformation, the conclusion holds. q.e.d. On the other hand, rank-two projectors constructed by Kramers pairs commute with both Ĥ and ˆT. Proposition 2. If [ ˆT, Ĥ] = 0, rank-two projectors formed by Kramers pairs commute with ˆT. Proof. By means of the Kramers pairs {ˆρ k, ˆη k } defined above, one can construct rank-two projectors as ˆP k = ˆρ k + ˆη k, which satisfies ˆP k 2 = ˆP k. Then, [ ˆT, ˆP k ] = ˆT ˆρ k ˆρ k ˆT + ˆT ( ˆT ˆρk ˆT ( ˆT ˆρ k ˆT ˆT, = ˆρ k ( ˆT + ˆT = 0, where it was used that ˆT 2 = Î plus ˆT ˆT = Î imples ˆT + ˆT = 0. q.e.d. Consequently, in addition of being separable states, the maximal mixed state, mixed 4

15 states with (P, P 2, P 3, P 4 = (/2 b, /2 b, b, b and its permutations, possess Kramers invariance. 6. Summary and Conclusions We provide a self-contained method to determine the extremal density matrices of a finite dimensional time-reversal Hamiltonian. These extremal states commute with the Hamiltonian operator and optimise its mean value, such that the conventional variational principle is extended to mixed states. We also establish a novel procedure to analyze the entanglement of extremal density matrices. It has the advantage of reaching any desirable extremal state, either separable or entangled, by changing the parameters of the Hamiltonian. It was applied for two families of cases of the 2-qubit Hamiltonian, in which, by varying its parameters, it is possible to keep it non degenerate or to acquire Kramers degeneracy. For the non degenerate case, we show that their associated extremal pure and mixed states are separable. When the Hamiltonian exhibits Kramers degeneracy, we have found both possibilities for extremal pure states, depending on the parameters of the Hamiltonian. For the extremal mixed matrices we also have both possibilities classified in five cases according to their eigenvalues degeneracy. The sufficiency of the PPT criterion in the qubit-qubit and qubit-qutrit systems makes that our procedure has no ambiguities and it can be applied to any observable by replacing the Hamiltonian. In higher dimensions, it can be implemented with the respective consideration that the PPT criterion is just a necessary condition. We want to enhance that our procedure encompasses Hamiltonian and states in the same context, thus, it is possible to discuss how the Hamiltonian degeneracy, its symmetry, the purity and entanglement of states in finite dimensional Hilbert space are intertwined. This is neither clear nor direct in the context of the diagonalisation procedure using the secular equation. In comparison with other separability criteria, the advantage of linking the extremal density matrices method with the PPT criterion is the algebraic aspect of the approach, i.e., the posivity conditions on the partial transposed density matrix (inequalities (9 separate in explicit way the regions for which there will be or not entanglement, for both pure and mixed states, by varying the parameters of the Hamiltonian. Consequently, this makes the procedure general and simple, without the need to introduce extra concepts. 5

16 A. Positivity conditions for the Density Operator The characteristic polynomial P d (x for the density matrix acting on a d-dimensional Hilbert space is given by P d (x det(xîd ˆρ = d ( j a j x d j = 0, (3 with the definitions a 0 = a and a d = det ˆρ. For d 2, the real coefficients {a k } are bounded by [22, 23] 0 a k ( d d k, (32 k where ( d k denotes a binomial coefficient. The upper bound defines the most mixed state with maximum von Neumann entropy, while the lower bound specifies pure states which has zero entropy. Additionally, it is known that a k = 0 for all values of k > rank(ˆρ [24]. j=0 The coefficients {a k } can be obtained by means of the Girard-Waring formula [24,28] a k = p(k i= k j= ( (j q ij j q ij q ij! (t j q ij, (33 where t j Tr(ˆρ j, for j =,..., d and {q ij, p(k} denote, respectively, the natural numbers solutions and (the partition function p(k the number of solutions, without regarding to order, of the linear Diophantine equation q i + 2 q i2 + 3 q i3 + + k q ik = k, (34 with i =, 2,..., p(k. For the first four values of k, p(k and {q ij } are given by p(k = = q =, p(k = 2 = 2 { (q, q 2 = (2, 0, (q 2, q 22 = (0, }, { p(k = 3 = 3 p(k = 4 = 5 (q, q 2, q 3 = (3, 0, 0, (q 2, q 22, q 23 = (,, 0, (q 3, q 32, q 33 = (0, 0, }, (35 (q, q 2, q 3, q 4 = (4, 0, 0, 0, (q 2, q 22, q 23, q 24 = (2,, 0, 0, (q 3, q 32, q 33, q 34 = (0, 2, 0, 0, (q 4, q 42, q 43, q 44 = (, 0,, 0, (q 5, q 52, q 53, q 54 = (0, 0, 0,. 6

17 The inverse of the Girard-Waring formula (33 exists and its given by [28] t k = k p(k i= ( M i! k j= ( (j q ij q ij! (a j q ij, (36 where M i k s= q is. The formulas (33 and (36 can also be expressed in Plemelj-Smithies form [29], t a t 2 t a 2 a , t k = a k = k! t k t k 2 t k 3... k t k t k t k 2... t (k a k a k 2 a k 3... ka k a k a k 2... a The density matrix must satisfy the following three properties: (a It has trace one; (b all its eigenvalues are positive or zero; and (c it is Hermitian. Given a monic real polynomial, the inverse problem of deciding when it comes from a density matrix requires these assumptions being translated into polynomial conditions. In other words, the traceone requisite is equivalent to a 0 = a =, the semi-positivity condition implies that expressions (32 must be fulfilled, and the hermiticity condition is taken into account through the Bezoutian matrix B d, i.e., a polynomial with real coefficients has reals roots iff B d is positive definite [25]. In terms of t j Tr(ˆρ j, with j =,..., d, the symmetric Bezoutian matrix is defined by [25 27] d t t 2 t d. t t 2 t.. 3 td B d =. t 2 t.. 3 td+. ( t d t d t d+ t 2(d For a given monic real polynomial in x of degree d, P d (x, having roots {P, P 2,..., P d }, its associated Bezoutian matrix B d has the following properties: a The rank of the Bezoutian equals the number of distinct roots of P d (x [25]. b (Sylvester criteria. The number of real roots of P d (x equals the signature (the difference between positive and negative real roots of its Bezoutian [25]. c (Reality condition. P d (x has all its roots real and distinct iff the Bezoutian matrix is positive definite [25]. 7.

18 d In the case of d = 2, 3, det B d is the only positivity condition of the Bezoutian [30]. e (Degeneracy condition. The discriminant of P (x is equal to the determinant of the Bezoutian, thus the condition for repeated roots of P (x is obtained by the vanishing of det B d [3]. f By means of the Vandermonde matrix, V d = P P 2 P d P 2 P2 2 Pd 2, ( P d P2 d P d d the Bezoutian is obtained as B d = V d V T d. Combining all above results, a monic real polynomial coming from a density matrix must satisfy the following system of d simultaneous polynomial equations: c k = a k, for k = 2, 3,..., d, (39 where constants c k fix the degree of mixing of the system, and the compatible region among them is obtained with the intersection of the semi-positivity conditions of the density matrix from (32 with the respective positivity conditions of the Bezoutian matrix [27]. For instance, taking into account the equality (39, the semi-positivity conditions (32 of the density matrices with dimensions d = 4 are given by and the respective Bezoutian matrix is B 4 = 0 c 2 3 8, 0 c 3 6, (40 0 c 4 256, 4 t 2 t 3 t 2 t 3 t 4 t 2 t 3 t 4 t 5 t 3 t 4 t 5 t 6, 8

19 where, the relations between {a p } with {t k } given in Eq. (33 yields, t 2 = 2 c 2, t 3 = 3c 2 + 3c 3, t 4 = 2(c 2 2c 2 + 4c 3 4c 4 +, (4 t 5 = 5c 2 (c 2 c 3 + 5c 3 + 5c 4 +, t 6 = 9c 2 2 2c 3 2 6(2c 3 + c 4 + c 2 + 3c 3 (c c 4 +. All the positivity conditions on B 4 are ( 2 ( (TrB 4 3 3TrB 4 TrB TrB4 3 6 TrB 4 0, (TrB 4 2 TrB4 2 0, (42 0, det B 4 0, where the last one is the main condition. Nevertheless, the remaining ones are crucial to avoid fake points in the compatible region for {c 2, c 3, c 4 }. Hence, for the set {c 2, c 3, c 4 }, the region which satisfies the inequalities system formed by (40 and (42, is shown in Fig.. Notice that, by making zero c 4, we obtain the d = 3 result, while by making zero two eigenvalues of the density matrix the line associated to the case d = 2 is obtained (c 3 = c 4 = 0. Inside the solid figure (orange color one has the solution for 4 eigenvalues of the density matrix different from zero, whereas the surfaces are associated to 2 degenerated eigenvalues (blue color. The curve for the case with three equal eigenvalues and the other different is also shown (green color. Acknowledgments. This work was partially supported by CONACyT-México (under Project No References [] A. Figueroa, J. López, O. Castaños, R. López-Peña, M. A. Man ko, and V. I. Man ko, Entropy-energy inequalities for qudit states J. Phys. A: Math. Theor (205. [2] O. Castaños, A. Figueroa, J. López, and R. López-Peña, Extremal Density Matrices for the Expectation Value of a Qudit Hamiltonian, J. of Phys: Conf. Ser (207; 9

20 A. Figueroa, J. López, O. Castaños, and R. López-Peña, Extremal density matrices for qudit states, arxiv: , (206. [3] A. Figueroa, O. Castaños, R. López-Peña, and G. Marmo, The variational method for density states a geometrical approach, J. Phys. A: Math. Theor (207. [4] O. Gühne and G. Tóth, Entanglement detection, Phys. Reports 474 (2009. [5] M. Horodecki, P. Horodecki, and R. Horodecki, Separability of mixed states: necessary and sufficient conditions, Phys. Lett. A 223 (996. [6] A. Peres, Separability Criterion for Density Matrices, Phys. Rev. Lett (996. [7] F. Haake, Quantum Signature of Chaos, Springer, (200. [8] J. Schlienz and G. Mahler, Description of entanglement Phys. Rev. A (995. [9] G. Mahler and V. Weberruß, Quantum Networks. Dynamics of Open Nanostructures, Springer, (998. [0] J. Keller, Multiple eigenvalues, Linear Algebra Appl (2008. [] W. J. Caspers, Degeneracy of the eigenvalues of hermitian matrices, J. of Phys.: Conference Series (2008. [2] R. Horn and C. Johnson, Matrix Analysis, Cambridge University Press (203. [3] B. L. van der Waerden, Modern Algebra II, Springer (99. [4] B. Sturmfels, Solving Systems of Polynomial Equations, Number 97, AMS Regional Conference Series (2002. [5] T. F. Jordan, Linear Operators for Quantum Mechanics, Dover (997. [6] N. Rösch, Time-reversal symmetry, Kramers degeneracy and the algebraic eigenvalue problem, Chemical Physics 80 (983. [7] I. Bengtsson and K. Zyczkowski, Geometry of Quantum States. An introduction to Quantum Entanglement, Cambridge University Press (207. [8] P. Aniello, J. Clemente-Gallardo, G. Marmo, and G. F. Volkert, Classical Tensors and Quantum Entanglement I: Pure States, Int. J. Geom. Methods Mod. Phys (200. [9] C. Quesne, SU(2 SU(2 scalars in the enveloping algebra of SU(4, J. Math. Phys. 7, 452 (976. [20] V. P. Gerdt, A. M. Khvedelidze, and Y. G. Palii, Constraints on SU(2 SU(2 invariant polynomials for a pair of entangled qubits, Phys. Atom. Nuclei (20. 20

21 [2] R. A. Bertlmann, H. Narnhofer, and W. Thirring, Geometric picture of entanglement and Bell inequalities, Phys. Rev. A (2002. [22] S. M. Deen, P. K. Kabir, and G. Karl, Positivity Constraints on Density Matrices, Phys. Rev. D (97. [23] M. S. Byrd and N. Khaneja, Characterization of the positivity of the density matrix in terms of the coherence vector representation, Phys. Rev. A (2003. [24] V. Tapia, Invariants and polynomial identities for higher rank matrices, J. Phys. A: Math. Theor (2007; V. Tapia, arxiv:math-ph/ (2002. [25] C. Procesi. Lie Groups: An Approach through Invariants and Representations, Springer (2007. [26] C. Procesi and G. Schwarz, Inequalities defining orbit spaces, Invent. Math (985. [27] V. P. Gerdt, A. M. Khvedelidze, and Y. G. Palii, Describing the Orbit Space of the Global Unitary Actions for Mixed Qudit States, J. Math. Sci (204. [28] H. W. Gould, The Girard-Waring power sum formulas for symmetric functions and Fibonacci sequences, Fibonacci Quarterly (999. [29] I. Gohberg, S. Goldberg, and N. Krupnik, Traces and Determinants of Linear Operators, Birkhäuser Verlag (2000. [30] C. P. Niculescu and L. E. Persson, Convex Functions and their Applications. A Contemporary Approach, Springer (2006. [3] M. Bhattacharya and C. Raman, Detecting level crossings without solving the Hamiltonian. I. Mathematical background, Phys. Rev. A (2007; M. Bhattacharya, How to detect level crossings without looking at the spectrum, Am. J. Phys (