arxiv: v1 [math-ph] 10 Apr 2012

Transcription

1 Simple Derivation of the Linbla Equation Philip Pearle 1, 1 Department of Physics, Hamilton College, Clinton, NY 1333 (Date: April 11, 01) The Linbla equation is an evolution equation for the ensity matrix in quantum theory. It is the general linear, Markovian, form which ensures that the ensity matrix is hermitian, trace 1, positive an completely positive. Some elementary examples of the Linbla equation are given. The erivation of the Linbla equation presente here is simple in that all it uses is the expression of a hermitian matrix in terms of its orthonormal eigenvectors an real eigenvalues. Thus, it is appropriate for stuents who have learne the algebra of quantum theory. Where helpful, arguments are first given in a two-imensional hilbert space. arxiv: v1 math-ph 10 Apr 01 PACS numbers: I. INTRODUCTION The ensity matrix ρ is a useful operator for quantum mechanical calculations. For a given system, one may be unsure about what is the state vector. If the possible state vectors an their associate probabilities are { ψ i, p i }, one creates the proper1 ensity matrix ρ i p i ψ i ψ i. (1) It is hermitian: ρ = ρ. It is trace 1: 1 = T rρ j φ j ρ φ j, where { φ j } are an arbitrary complete orthonormal set of vectors. It is positive: v ρ v 0 for an arbitrary vector v. All these properties can easily be verifie from Eq.(1). One can use the ensity matrix to conveniently calculate probabilities or mean values. If a measurement is set up to result in one of the eigenstates φ i of an operator O, that outcome s probability is φ i ρ φ i an the mean eigenvalue of O is T roρ. Because each iniviual state vector evolves unitarily uner the system Hamiltonian H (assume for simplicity here to be time inepenent), ψ, t = exp ith ψ, 0, the ensity matrix in Eq.(1) satisfies the evolution equation t ρ(t) = ih, ρ(t) = ρ(t) = e iht ρ(0)e iht. () The operator acting on ρ(0) is often calle a superoperator since it escribes a linear transformation on an operator: it operates on both sies of ρ, so to speak. The case sometime arises where the system S uner consieration is a subsystem of a larger system S + S, an S is not measure. The pure (so-calle because it is forme from a single state vector) ensity matrix for the joint system is R(t) = im C im (t) φ i χ m jn Electronic aress: ppearle@hamilton.eu C jn (t) φ j χ n, where { φ i } { χ m } are orthonormal bases for S, S respectively, an im C im(t) = 1. By taking the trace of R with respect to S, one arrives at an improper1 ensity matrix for S from which preictions can be extracte: ρ(t) = ijm C im (t)c jm (t) φ i φ j. One can easily see that this is hermitian, trace 1 an positive. However, while the ensity matrix of S + S evolves unitarily, the ensity matrix of the subsystem S evolving uner the influence of S generally oes not evolve unitarily. Nonetheless, sometimes ρ(t)/t can be written in terms of ρ for a range of times earlier than t. Sometimes that range is short compare to the time scale of evolution of ρ so that one may make an approximation whereby ρ(t)/t epens linearly just on ρ(t). This is quite useful, an is what shall be consiere in this paper. In this case, the evolution equation is highly constraine by the requirements on ρ(t), to be satisfie at all times: hermiticity, trace 1 an positivity. (The latter proves too general to simply implement, so a stronger requirement is impose, calle complete positivity see Section IV).) The result, for an N-imensional Hilbert space, is the Linbla3 (or Linbla-Gorini- Kossakowsky-Suarshan4) evolution equation for the ensity matrix: ρ(t) = ih, ρ(t) t 1 1 L L ρ(t) + ρ(t)l L L ρ(t)l. (3) In Eq.(3), the hamiltonian H is an arbitrary hermitian operator, but the Linbla operators {L } are completely arbitrary operators. Actually, as shall be shown, there nee be no limitation on the number of terms in the sum in Eq.(3), but this can always be reuce to a sum of 1 terms. It is not a necessary conition, but if the equation is to be time-translation-invariant, the operators are timeinepenent.

2 II. LINDBLAD EXAMPLES Before eriving Eq.(3), we give a few examples of the non-unitary evolutions it escribes. Since unitary evolution is well known, we shall let H = 0. It shall be seen that relaxation to some equilibrium (constant) ensity matrix is reaily escribe. For simplicity, four of the five examples shall be in a N = hilbert space (a restriction that is reaily lifte). A. Ranom Phases Consier a state vector written in a basis whose phase factors unergo ranom walk. Given an initial state vector ψ = a φ 1 + b φ ( a + b = 1), suppose at time t, it has evolve to with probability ψ, t = ae iθ1 φ 1 + be iθ φ P (θ 1, θ )θ 1 θ = θ 1 θ 1 1 πλ1 t πλ t e The ensity matrix is ρ(t) = θ 1 θ P (θ 1, θ ) ψ ψ = a φ 1 φ 1 + b φ φ λ 1 t θ 1 e 1 λ t θ. +e 1 tλ1+λ ab φ 1 φ + a b φ φ 1. We see that the off-iagonal elements ecay at a fixe rate while the iagonal elements remain constant. It satisfies t ρ(t) = 1 λ 1 + λ ρ(t) Q i ρ(t)q i, i=1 where the projection operator Q i φ i φ i. To see that this is a Linbla equation, note that Q i = Q i an i=1 Q i Q i = 1. There are two Linbla operators: ientify L i (λ 1 + λ )Q i in Eq.(3). B. Unitary Jump Suppose in time t, a state vector ψ, t has probability λt of changing to exp ig ψ, t (probability 1 λt of being unchange), where G is a hermitian operator an exp ig 1. The ensity matrix at time t + t is therefore ρ(t + t) = (1 λt)ρ(t) + λte ig ρ(t)e ig, so its evolution equation is t ρ(t) = λρ(t) e ig ρ(t)e ig. This is of the Linbla form, with one Linbla operator L λ exp ig. In the basis where G is iagonal with elements (g 1, g ), we get ρ ii /t = 0 an t ρ 1(t) = λρ 1 (t) 1 e i(g g1). So, again, its iagonal elements remain constant. Its offiagonal elements ecay at the fixe rate λ1 cos(g g 1 ) an their phases change. C. Ranom Unitary Transformation Suppose in time t a state vector ψ, t unergoes a unitary transformation to e igθ ψ, t = 1 igθ 1 (Gθ) +... ψ, t with probability P (θ)θ = θ 4πλt e 1 4λt θ. The ensity matrix at t + t (neglecting terms of orer higher than t) is given by ρ(t + t) = P (θ)θe igθ ρ(t)e igθ = P (θ)θ ρ(t) θ 1 G ρ(t) + ρ(t) 1 G Gρ(t)G = ρ(t) λt giving the Linbla equation G ρ(t) + ρ(t)g Gρ(t)G t ρ(t) = λ G ρ(t)+ρ(t)g Gρ(t)G = λ GG, ρ(t) with one Linbla operator L λg. In the basis where G is iagonal with elements (g 1, g ), t ρ ij(t) = λ (g i g j ) ρ ij (t). So, again, its iagonal elements remain constant but its off-iagonal elements ecay at a rate etermine by the ifference in eigenvalues. D. State Exchange Suppose, in time t, with probability λt, a state vector ψ, t exchanges its basis states φ 1, φ, becoming σ ψ, t, where φ i σ φ j is the Pauli matrix with iagonal elements 0 an off-iagonal elements 1. It is easy to see that the ensity matrix evolution equation is

3 3 t ρ(t) = λ ρ(t) σρ(t)σ an is of the Linbla form, with one Linbla operator L λσ. The ensity matrix elements therefore satisfy t ρ 11(t) = t ρ (t) = λρ 11 (t) ρ (t), t ρ 1(t) = t ρ 1(t) = λρ 1 (t) ρ 1 (t). The iagonal ensity matrix elements change in this example, ecaying to 1/. The off-iagonal matrix elements keep their real parts while the imaginary parts ecay to 0. E. State Transitions Here we consier arbitrary N. Suppose in time t, a state vector ψ, t makes a transition to state with probability n ψ, t m n ψ, t P mn p m λt ψ, t n (p m 0, N m=1 p m = 1). The probability of all such transitions is mn P mn = λt, so the state vector is unchange with probability 1 λt. Define Q mn m n. Note that N m,n=1 p mq mnq mn = 1. The ensity matrix at time t + t is ρ(t + t) = (1 λt)ρ(t) + λt p m Q mn ρ(t)q nm, so its evolution equation is t ρ(t) = λ ρ(t) m,n=1 m,n=1 p m Q mn ρ(t)q nm. This is of the Linbla form, with (one more than the necessary maximum!) Linbla operators L mn λpm Q mn : The matrix elements of the ensity matrix obey t ρ rs(t) = λ ρ rs (t) p r δ rs ρ nn (t) = λ ρ rs (t) p r δ rs. The off-iagonal elements ecay at a uniform rate. The iagonal elements o not remain constant. They ecay to preetermine values p r : ρ rr (t) = ρ rr (0)e λt + p r 1 e λt. This might be useful in moeling the approach to thermal equilibrium, where the states m are energy eigenstates an p r is the Boltzmann probability Z 1 exp E r /kt. n III. APPLICATION OF CONSTRAINTS We now turn to eriving the Linbla equation as the most general equation satisfying the constraints. While the hilbert space iscusse here shall be assume of imension N, N may be allowe to go to to infinity an, also, the argument may reaily be extene to a continuum basis. To make the argument easier to follow, examples of how its steps apply to a two-imensional hilbert space shall occasionally be inserte. The Markov constraint is that the ensity matrix ρ(t ) ρ at a later time t, epens only upon the ensity matrix ρ(t) ρ at an earlier time t, not upon the ensity matrix over a range of earlier times. The linearity constraint, combine with the Markov constraint, is that the matrix elements of ρ can be written as the sum of constants multiplying the matrix elements of ρ rather than, say, powers of the matrix elements of ρ or any other kin of function of these matrix elements: ρ ij = A ir,js ρ rs. (4) Here ρ ij φ i ρ φ j with φ j some convenient orthonormal basis an, similarly, ρ is expresse in the same basis. The constants A ir,js can be functions of t, t. There are N 4 constants, an each can be complex, so there are N 4 real constants involve in Eq.(4). The hermiticity constraint ρ ij ρ ji = ρ ij, applie to Eq.(4), results in A js,ir A ir,js ρ rs. (5) 1. Two-imensional space: hermiticity Suppose we have the equation B 11 ρ 11 + B 1 ρ 1 + B 1 ρ 1 + B ρ = T rbρ = 0, (the B ij are constants) which hols for all possible ensity matrices. Then, one can see B = 0 as follows. First choose the ensity matrix ρ 11 = 1, with all other elements vanishing: thus, B 11 = 0. Similarly, one shows B = 0. Next, employ the ensity matrix ρ ij = 1/, which results in B 1 + B 1 = 0 Finally, use the ensity matrix ρ 11 = ρ = 1/, ρ 1 = ρ 1 = i/, which results in B 1 B 1 = 0 an so B 1 = B 1 = 0. Therefore, B = 0 The four ensity matrices use here, 1 (1 + σ3 ), 1 (1 σ3 ), 1 (1 + σ1 ), 1 (1 + σ ), (written in terms of the Pauli matrices) we shall call the ensity matrix basis. Any matrix can be written

4 4 as a linear sum with constant (complex) coefficients of these four matrices. More than that, they form a matrix basis for hermitian matrices, in that any hermitian matrix can be written as a linear sum with constant (real) coefficients of these four matrices. More than that, an this is the reason for their eployment here, they form a matrix basis for ensity matrices, in that any ensity matrix can be written as a linear sum with constant positive real coefficients of these four matrices such that the sum of the coefficients a up to 1. This basis is to be istinguishe from another basis, the Pauli matrices plus the ientity matrix, which we shall call the Pauli+1 basis. This is also a matrix basis for hermitian matrices but it is not a ensity matrix basis.. N-imensional space: hermiticity Generalizing, if we have an equation B sr ρ rs T rbρ = 0 (6) for a matrix B, which hols for all vali ρ, then B = 0. This can be seen by using an -size ensity matrix basis, (generalizing the -size ensity matrix basis of the previous section). First choose ρ kk = 1 with all other elements vanishing, which implies B kk = 0. Then for particular values of k, l, choose ρ kk = ρ ll = ρ kl = ρ lk = 1/ with all other elements vanishing, from which one fins B kl + B lk = 0. Finally, choose ρ kk = ρ ll = iρ kl = iρ lk = 1/, from which one fins B kl B lk = 0, so B kl = B lk = 0. Letting k, l range over all possible pairs of inices results in B = Two-imensional space: evolution equation an trace constraint It therefore follows from Eq.(5) that A js,ir = A ir,js where each inex can take on the values 1 or. A matrix B for which Bm,n = B n,m is a hermitian matrix. Therefore, A is a hermitian matrix, where we regar the number pairs 11, 1, 1, as four ifferent inices. That is, A is a 4 4 imensional matrix. The most general 4 4 hermitian matrix is characterize by 16 real numbers (the four real iagonal matrix elements an the six complex matrix elements above the iagonal). Since there are 3 real numbers which characterize the most general superoperator in a two imensional space, the conition of hermiticity of the ensity matrix has cut that number in half. A hermitian matrix can be written in terms of its orthonormal eigenvectors an eigenvalues, an that ecomposition shall prove very useful here. There are four real eigenvalues, λ, where = 1,, 3, 4. Corresponing to each eigenvalue is an eigenvector E in the four imensional complex vector space. The four complex components of Eir make 3 real numbers, but they are constraine. Each eigenvector is normalize to 1: i,r=1 E ir E ir = 1 + i0 provies 8 constraints, lowering the number of free components to 4. The orthogonality of E 1 to the other three vectors provies 6 constraints, the orthogonality of E to the remaining two vectors provies 4 constraints, an the orthogonality of E 3 to E 4 provies two constraints. Thus, there are 1 constraints on the 4 free components, so the eigenvectors contain 1 free components. These, together with the four eigenvalues, comprise the 16 real numbers characterizing A. An example of such an orthonormal basis is given by 1/ multiplying the Pauli+1 basis. If we write the four components of E as a four imensional vector with components E11, E1, E1, E, then E 1 1/ σ 1 has components 1/ 0, 1, 1, 0, E 1/ σ has components 1/ 0, i, i, 0, E 3 1/ σ 3 has components 1/ 1, 0, 0, 1, E 4 1/ 1 has components 1/ 1, 0, 0, 1, It is easy to verify that this is an orthonormal set of vectors. Although each E is a vector in a four imensional space, with four components Eij, E can also be regare as an operator in the two-imensional hilbert space with four matrix elements Eij. This leas to a neat way of writing the orthogonality relations for these eigenvectors. Instea of N N i=1 r=1 E ir Eβ ir = δβ, we can write T re E β = δ β, where E β is the hermitian conjugate (complex conjugate transpose) of E β. It is easy to see how this works for the example where E is 1/ the Pauli+1 basis. The expression for the components of A written in terms of its eigenvectors an eigenvalues is A ir,js = 4 λ E ire js. Putting this into Eq.(4) results in the evolution equation ρ ij = 4 λ E irρ rs E js = 4 λ E ρe. Now, lets impose the trace constraint, i.e., i=1 ρ ii = 1. In terms of components this says 1 = 4 λ i=1 E irρ rs E is. Writing 1 = i=1 ρ ii = δ srρ rs, the trace constraint can be written as 4 λ E si E ir δ sr ρ rs = 0 i=1

5 5 or in matrix notation as 4 T r λ E E 1 ρ = 0. where 1 is the unit matrix. This must hol for arbitrary ρ. We have seen how to hanle such an expression. By successively putting in the four ensity basis matrices, we obtain the trace constraint 4 λ E E = N-imensional space: trace constraint The N-imensional case works just like the twoimensional case. It follows from Eq.(5) that A can be viewe as an hermitian matrix. It has real eigenvalues. Its complex eigenvectors Eir satisfy the orthonormality conitions IV. COMPLETE POSITIVITY The final constraint is positivity. This says, given an arbitrary N-imensional vector v, that the expectation value of the ensity matrix ρ is non-negative. This constraint, applie to Eq.(8), is 0 v ρ v = λ v E ρe v = λ v ρ v, (10) where we have efine E v v. Positivity of ρ ensures v ρ v 0. Thus, we see from Eq.(10), if all the λ s are non-negative, then ρ will be positive too. However λ 0, while just shown to be sufficient for ρ to be positive, is not necessary. In the next section, we shall give an example where an eigenvalue is negative, yet ρ is positive! Therefore, a stronger conition than positivity is necessary to ensure that λ 0. This conition, presente after the example, is complete positivity. i=1 r=1 E ire β ir = T re E β = δ β. (7) With A written in terms of its eigenvalues an eigenvectors, Eq.(4) becomes the evolution equation ρ ij = ρ = λ N E ire js ρ rs or λ E ρe. (8) λ an E epen upon t t, but we shall not write that epenence until it is neee. Next, imposition of the trace constraint on Eq.(8), with T rρ = 1 = T r1ρ, gives T r λ E E 1 ρ = 0. Using the ensity matrix basis as in Eq.(6) et seq., we obtain the trace constraint: λ E E = 1. (9) By taking the trace of Eq.(9) an using Eq.(7), we fin the interesting relation λ = N. 5. Two-imensional space: example of a positive ensity matrix with negative eigenvalue This example uses the Pauli+1 eigenvectors E = 1/ σ an 1/ 1. (Note that the trace constraint (9) is satisfie, provie 4 λ =, since the square of each of the Pauli+1 matrices is 1/ 1.) Choose λ 1 = λ = λ 3 = λ 4 = 1: ρ = 1 σ1 ρσ 1 + σ ρσ σ 3 ρσ 3 + 1ρ1 = ρ ρ 1. ρ 1 ρ 11 ρ is just ρ with its iagonal elements exchange. Thus, because ρ is positive, then ρ is positive. This is a particularly simple example of a more general case iscusse in Appenix A. 6. N-imensional space: efinition of complete positivity It is not positivity but, rather, complete positivity that makes the non-negative eigenvalue conition necessary. Here is what it means. A to our system a non-interacting an non-evolving aitional system in its own N-imensional hilbert space. The enlarge hilbert space is of imension. The simplest state vector in the enlarge space is a irect prouct φ i χ j : φ i is a vector from the original hilbert space, χ j is a vector from the ae system. The general state vector in the joint space is the sum of such proucts with c-number coefficients. Form an arbitrary ensity matrix R for the enlarge system. Suppose it evolves accoring to Eq.(8), where

6 6 E is replace by E 1 (i.e., the evolution has no effect on the vectors of the ae system.) Complete positivity says that the resulting ensity matrix R must be positive. 7. Two-imensional space: complete positivity Complete positivity says, given the evolution equation (8), that w R w 0 for an arbitrary imensional vector w an for any initial ensity matrix R in the enlarge hilbert space. We wish to prove that complete positivity implies the eigenvalues are non-negative. What we shall o is juiciously choose a single vector w an four pure ensity matrices R so that the expressions w R w = λ w E RE w, are λ β, with a positive constant of proportionality. Therefore, for complete positivity to hol, λ β must be non-zero. Here are choices that will o the job. We shall choose the maximally entangle vector w 4 φ r χ r. r=1 ( w w = 4, but it nee not be normalize to 1). construct the state vectors 4 ψ β E β ij φ i χ j, i,j=1 We an use them to make four pure ensity matrices ψ β β ψ. (Note that T r ψ β β ψ = 1 because of the orthogonality relation Eq.(7)). Then, for one β, R = 4 i,j,i,j =1 E β ij Eβ j i φ i χ j φ i χ j. Putting this into Eq.(8), the complete positivity conition is 4 0 w R w = λ E β ij Eβ j i,r,r,i,j,i,j =1 φ r χ r E φ i χ j φ i χ j E φ r χ r 4 = λ E β ij Eβ j i E jie ij = =,i,j,i,j =1 4 λ T r(e β E )T r(e β E ) 4 λ (δ β ) = λ β (using the orthogonality relation (7)). Thus, complete positivity implies λ β N-imensional space: complete positivity We follow the same proceure in the N-imensional case. However, to be a bit more general, we shall use an arbitrary vector w, an an arbitrary pure ensity matrix R: w D mn φ m χ n (11a) R m,n=1 k,l,k,l =1 C kl C k l φ k χ l φ k χ l (11b) where C kl, D mn are yet to be specifie complex constants. The unit trace of ρ in Eq.(11b) requires T rc C = 1. Then, the complete positivity conition is 0 w R w = λ w E RE w = = λ N 1 D m n D mnc kl C k l E m ke mk δ n lδ l n λ T rcd E T re DC. (1) Now, choose DC = E β, for any particular β. This choice can be mae in many ways. Two are C = E β, D = 1 (the choice mae in the two-imensional example just iscusse) or D = E β, C = N 1/ 1 (note, both choices respect T rc C = 1). With this choice in Eq.(1), an with use of the orthonormality conitions Eq.(7), we obtain as the consequence of complete positivity: 0 λ (δ β ) = λ β for 1 β. (13) V. KRAUS REPRESENTATION We have now applie all the constraints neee to obtain a vali ensity matrix ρ at a later time t from an earlier ensity matrix ρ at time t. This relation is Eq.(8), supplemente by the orthonormality conitions (7), the trace constraint (9) an the conition of non-negative eigenvalues (13). It is customary to efine M λ E, so that Eqs.(8, 9) can be written in terms of M alone: ρ = M ρm, M M = 1. (14a) (14b) (However, the orthonormality conitions, written in terms of M, now epen upon λ ). Eq.(14a) is calle

7 7 the Kraus representation an {M } are calle Kraus operators5. We have prove the necessity of the Kraus representation, but it is also sufficient. That is, for any {M } satisfying Eqs.(14a,14b), even for more than operators, also with no orthonormality conitions impose, all the constraints on ρ are satisfie. It is easy to see that hermiticity, trace 1 an positivity are satisfie. Complete positivity requires a bit more work, an that is given in Appenix B. This general statement of the Kraus representation might seem to imply a larger class than we have erive as necessary, but that is not so. Since the Kraus representation is hermitian, trace 1 an completely positive, it may be written in the form Eq.(8), as we have shown. For β, the eigenvalues vanish. Note that the conition T re β 1 = 0 says that these eigenvectors are orthogonal to E 1. An, inee, in this case, Eq.(8) becomes the ientity ρ(t) = 1ρ(t)1 = N 1 N 1ρ(t) 1 N 1. (18) B. Eigenvectors an eigenvalues when t = t + t When t = t + t, the eigenvalues an eigenvectors change infinitesimally. Accoringly we write VI. LINDBLAD EQUATION Now that we have satisfie all the constraints on the ensity matrix ρ ρ(t ), we can let t = t + t, an obtain the ifferential equation satisfie by ρ(t). For the rest of this paper we shall only treat the N-imensional case since the argument is precisely ientical for the twoimensional case, except that N =. A. Eigenvectors an eigenvalues when t = t First, lets see what we can say about the eigenvectors an eigenvalues when t = t. Then, Eq.(8) says ρ ij = 0 = λ N N E ire sj ρ rs or λ E ire sj δ riδ js ρ rs. (15) As we have one before, successive replacement of ρ by the members of the ensity matrix basis results in δ ri δ js = λ E ire sj. (16) Multiply Eq.(16) by E β js an sum over j, s. Use of the orthonormality relation (7) gives δ ri T re β = λ β E β ir. (17) If T re β 0 an λ β 0, Eq.(17) says that all the eigenvectors are 1. But only one of a set of orthogonal eigenvectors can be proportional to the ientity. Therefore, for the rest of the eigenvectors, λ β = 0 an T re β = 0. Call one eigenvector E N 1/ 1. From Eq.(17), we fin the associate eigenvalue λ = N. λ (t) = N1 c t, λ (t) = c t ( ), E (t) = 1 N 1 + Bt, E (t) = K ( ), (19) where the c are constants. We o not inclue a term t in the expression for E (t) since, because λ (t) t, it woul contribute a negligible term (t) to Eqs.(8,9). Because the eigenvalues must be positive, an because the eigenvalues sum to N (equation following Eq.(9)), we see that c 0 (all ). B an K are restricte by the orthonormality conitions, which we shall look at later. C. The evolution equation Putting Eqs.(19) into the evolution equation (8) gives ρ(t + t) = 1 c t1 + Btρ(t)1 + B t 1 +t c K ρ(t)k, or in the limit t 0, N 1 ρ(t) = cn ρ(t) + Bρ(t) + ρ(t)b + c K ρ(t)k. t Putting Eqs.(19) into the trace constraint (9) gives N1 c t1 + B t1 + Bt 1 +t c K K = 1, or, 1 c 1 = B + B + c K K. (1) (0)

8 8 Using (1) to replace c in (0) (specifically, c ρ = (1/)c 1ρ + ρc 1) results in 1 t ρ(t) = (B B ), ρ(t) 1 3 c K K ρ(t) + ρ(t)k K K ρ(t)k. D. The Linbla equation () If we efine ih (1/)(B B) an L = c K, the evolution equation () becomes the Linbla equation (3): t ρ(t) = i H, ρ(t) 1 L L ρ(t) + ρ(t)l L L ρ(t)l. 1 E. Orthonormality conitions (3) It is a consequence of this erivation that the Linbla operators L in Eq.(3) are not arbitrary operators, because they are restricte by the orthonormality conitions (7). Putting Eqs.(19) into Eq.(7) constrains B, K. For what follows, we recall from the iscussion in section II that {K } can be regare in two ways. In one way, they are regare as 1 operators acting on vectors in an N imensional space, with matrix elements Kij (i, j = 1...N). In the other way, they are regare as 1 vectors in an imensional space, each with components (K11, K1,...KNN ). In particular, the trace of two operators is the same as the scalar prouct of two vectors, as in Eq.(7) The orthonormality relation (7), applie successively to ( = β = ), (, β = ), (, β ), with use of Eqs.(19), are T rb + B = 0, T rk = 0,, ( = 1,... 1) T rk K β = δ β (, β = 1,... 1). (4a) (4b) (4c) Eq.(4a) says that the hermitian part of B vanishes. This provies no restriction at all on H, which is the anti-hermitian part of B. Eq.(4c) says that the vectors K are orthonormal. Eq.(4b) says that T rk 1 = 0, which implies that N 1/ 1 completes the orthonormal set. F. The general Linbla form We shall now show that the Linbla equation (3) with arbitrary Linbla operators (no constraints whatsoever) can be transforme to new, constraine, Linbla operators of Eq.() by aing a constant (achieving the vanishing trace constraint (4b)) followe by a unitary transformation(achieving the orthogonality constraint (4c)). First, we see that we can transform the arbitrary Linbla operators L to Linbla operators L which are traceless. Define L L + k 1, where the k are 1 constants, an substitute that into the Linbla equation, obtaining : t ρ(t) = i H + i(k L k L ), ρ(t) 1 3 L L ρ(t) + ρ(t)l L L ρ(t)l. With a reefinition of H, this is again the Linbla equation, expresse in terms of L. By choosing k = N 1 T rl, the new Linbla operators satisfy T rl = 0. Now write the 1 Linbla operators L (hereafter assume traceless) in terms of 1 new operators L β (, β = 1,... 1) using the linear transformation L = 1 β=1 U,β Lβ. (5) We ask what the matrix U,β must be in orer that the Linbla form be unchange ( L replacing L in Eq.(3)). With arbitrary operators A 1, A, A 3, 1 A 1 L A L A 3 = 1,β,β =1 U,β U,β A 1 Lβ A Lβ A 3. (The three terms in the Linbla equation have two of the A = 1 while the thir A = ρ). This equals 1 β=1 A Lβ 1 A Lβ A 3, leaving the Linbla equation unchange in form, if an only if the matrix U,β is unitary, 1 U β, U,β = δ β,β. Inverting Eq.(5), we see that L β is traceless. Now, consier the matrix T rl L. It is hermitian, so its eigenvalues are real an its eigenvectors are orthogonal. It can be brought to iagonal form by properly choosing our unitary transformation, so we obtain T r L β Lβ = 1 U β, U,β T rl L = c β δ ββ. (6) This is almost the orthogonality constraint (4c). The eigenvalues c β are non-negative, since it follows from

9 9 Eq.(6) that c β = T r L β Lβ = i,j=1 L β ij L β ij 0. Thus, accoring to Eq.(6), the L β are orthogonal vectors, with square norm c β. We can a one more vector 1 to complete the set, orthogonal to the rest since T r L β 1 = T r L β = 0. We can efine new operators K β which are orthonormal an traceless, by L β c β K β. In terms of these operators, the Linbla equation (3) written in terms of L β becomes t ρ(t) = i H, ρ(t) 1 1 c K K ρ(t) + ρ(t) K K K ρ(t) K. (7) This is precisely Eq.() with K, c, replacing K, c. Moreover, the orthonormality constraints Eqs.(4b,4c) on K are satisfie by K. G. Concluing remarks We have shown that the constraints on the ensity matrix ρ manate its evolution equation (), where c 0 an B, K satisfy the constraints (4). We then showe that this is completely equivalent to the Linbla Eq.(3), with no constraints at all on the 1 Linbla operators L. However, there is no nee to restrict the Linbla equation to no more than 1 operators. We conclue our presentation by showing that the Linbla equation with any number of operators has all the require properties. (Of course, from what we have shown, such an equation may be reuce to one with no more than 1 operators.) Looking at the Linbla equation in that case, t ρ(t) = i H, ρ(t) 1 L L ρ(t) + ρ(t)l L L ρ(t)l,(8) it is easy to see that hermiticity, an trace 1 (in the form T rρ(t)/t = 0, with T rρ(t 0 ) = 1) are satisfie. Complete positivity requires a bit more work. We have shown in Appenix B that complete positivity hols for the Kraus form (14a) subject to the trace constraint (14b), with an arbitrary number of operators. So, if Eq.(8) can be written in the Kraus form with the trace constraint, we have shown it is completely positive. Accoringly, we write Eq.(8) as ρ(t + t) = 1 t(ih + 1 L L ) ρ(t) 1 t( ih + 1 L L ) + t L ρ(t)l. (9) Ientifying the Kraus operators as M 0 = 1 t(ih gives the Kraus form Eq.(14a), ρ(t + t) = L L ), M 0 = tl, (30) M ρ(t)m. (31) Now, take the trace of Eq.(31). Since Eq.(8) implies T rρ(t + t) = 1 = T r1ρ(t), the result is 1 = T r M M ρ(t) or T r M M 1 ρ(t). (3) As we have one before, by successively replacing ρ(t) by the members of the ensity matrix basis, we obtain the Kraus trace constraint (14b): T r M M = 1. (33) Therefore, the Linbla form with an arbitrary number of Linbla operators, is completely positive. Appenix A: A Class of Positive Density Matrices With Negative Superoperator Eigenvalues Consier the Pauli+1 matrices (multiplie by 1/ ) as the eigenvectors E β. As note in Section IV, they satisfy the trace constraint (9) if 4 λ =. Eq.(8) becomes ρ = 4 λ E ρe = 1 λ1 σ 1 ρσ 1 + λ σ ρσ + λ 3 σ 3 ρσ 3 + λ 4 1ρ1 = 1 (λ 1 + λ )ρ + (λ 3 + λ 4 )ρ 11, (λ 1 λ )ρ 1 (λ 3 λ 4 )ρ 1 (λ 1 λ )ρ 1 (λ 3 λ 4 )ρ 1, (λ 1 + λ )ρ 11 + (λ 3 + λ 4. )ρ (A1) Now, we can successively replace ρ by the four ensity basis matrices, an eman that the λ be chosen so ρ is positive for all. Since the sum of these four ρ s with positive coefficients (aing up to 1) is the most general two-imensional ensity matrix, then the most general ρ will be positive.

10 10 The first ensity basis matrix has ρ 11 = 1 an the rest of the matrix elements vanishing. Then, ρ = 1 (λ 3 + λ 4 ) 0 0 (λ 1 + λ. (A) ) Since the eigenvalues of the ensity matrix must lie between 1 an 0, we obtain the two conitions: λ 1 + λ 0 (A3a) λ 3 + λ 4 0. (A3b) The secon ensity basis matrix, with ρ = 1 an the rest of the matrix elements vanishing, gives the same results. The thir ensity basis matrix is ρ = (1/)1 + σ 1. Using 4 λ = to simplify the result, we obtain: ρ = 1 1 (λ 1 + λ 4 1) (λ 1 + λ 4. (A4) 1) 1 The eigenvalues of ρ here are (λ 1 +λ 4 )/, 1 (λ 1 +λ 4 )/, so the conition that they lie between 0 an 1 is λ 1 + λ 4 0. (A5) The fourth ensity basis matrix is ρ = (1/)1 + σ. Using 4 λ = to simplify the result, we obtain: ρ = 1 1 i(λ + λ 4 1) i(λ + λ 4. (A6) 1) 1 The eigenvalues of ρ here are (λ +λ 4 )/, 1 (λ +λ 4 )/, so the conition that they lie between 0 an 1 is λ + λ 4 0. (A7) So, we have obtaine the result that ρ will be positive if Eqs.(A3a, A3b,A5,A7) an 4 λ = (A8) are satisfie. The sum of Eqs.(A3b, A5, A7))minus (A8) tells us that λ 4 1. Eq.(A8) an the constraint bounaries are threeimensional hyperplanes in the four imensional λ-space. Their intersections elineate the allowe areas for the eigenvalues. We shall be content here to set λ = λ 1 (in which case Eqs.(A3a) simplifies to 1 λ 1 0). Then, Eq.(A8) escribes a plane in {λ 1, λ 3, λ 4 } space, an its intersection with the constraint bounary planes can be rawn. This is shown in Fig. 1. There are two regions where one of the eigenvalues is negative an the other two are positive: the points in the heavily outline upper left triangle have λ 3 0, an the points in the heavily outline lower right triangle have λ 4 0. FIG. 1: Allowe regions of eigenvalues for a positive ensity matrix ρ. The two ark-outline triangular regions are where an eigenvalue is negative. They abut an isosceles triangle, the restricte region of complete positivity, where all the eigenvalues are positive. Appenix B: Complete Positivity of the Kraus representation The Kraus form Eq.(14a) an the Kraus constraint Eq.(14b), generalize to any number of arbitrary operators M, are respectively ρ = M ρm, M M = 1. (B1a) (B1b) We want to show complete positivity. Call any one of the M M. If we can show complete positivity for MρM for an arbitrary M, then Eq.(B1a), which involves a sum of such terms, will be completely positive. An, it is only necessary to prove complete positivity for R = ρ ρ µ, where ρ µ is any possible basis ensity matrix (escribe in the paragraph following Eq.(6)) in the ae hilbert space, since the most general ensity matrix in the irect prouct hilbert space is the linear sum of such terms. We now calculate v R v = v MρM ρ µ v. (B) for arbitrary v = N m,n=1 D mn φ m χ n. There is no loss of generality if we pick the basis vectors χ n in the ae hilbert space any way we like. We shall pick them to be the eigenstates of ρ µ. Now, we note that each ρ µ has one eigenvalue 1 an the remaining N 1 eigenvalues are 0. Call the eigenvector χ 1 which correspons to the

11 11 eigenvalue 1. Then, v MρM ρ µ v = N = m =1 m,m =1 N Dm 1 φ m M ρ m =1 D m 1D m1 φ m MρM φ m D m1 M φ m 0, (B3) Therefore, the Kraus form with arbitrary M s is completely positive. 1 B. Espagnat, Conceptual Founations of Quantum Mechanics (W. A. Benjamin, California 1971). Superoperator, was first use by I. Prigogine an coworkers, e.g., I. Prigogine, C. George, F. Henin, P. Manel, J. W. Turner, Proc. Nat l. Aca. Sci. USA66, 709 (1970). 3 G. Linbla, Commun. Math. Phys (1976). 4 V. Gorini, A. Kossakowski an E. C. G. Suarshan, Journ. of Math. Phys. 17, 81 (1976). 5 K. Kraus, Ann. Phys. (N.Y.) 64, 311 (1971), States, Effects an Operations (Spring-Verlag, Berlin, 1983).