Microscopic derivation of a class of non-linear dissipative schrödinger-like equations

Microscopic derivation of a class of non-linear dissipative schrödinger-like equations

Physicu Publishing IllA (1982) 364-370 North-Hollnnd MICROSCOPIC NON-LINEAR DERIVATION DISSIPATIVE Co. OF A CLASS SCHRGDINGER-LIKE OF EQUATIO...

392KB Sizes 0 Downloads 34 Views

Physicu

Publishing

IllA (1982) 364-370 North-Hollnnd

MICROSCOPIC NON-LINEAR

DERIVATION

DISSIPATIVE

Co.

OF A CLASS

SCHRGDINGER-LIKE

OF EQUATIONS

N. GISIN* Dipartemrnf

de Physique

Theoriqrre.

Unicrrsite

Received

The general

formalism

of master equations

the pure states to derive a class of non-linear spin-i coupled

to a bath of two-level

make a connection

with Sz-Nagy’s

t/e GenPre.

6 August

I?!

I GvnPw 4.

1981

is used together Schrodinger-like

with a projection equations.

systems and the Rloch equations theorem

Swifrerland

on the dilations

which preserves

The result is applied to a are recovered.

of contracting

Finally

we

semi-groups.

1. Introduction

The master equation formalism has been widely studies and applied with success to numerous problems involving an irreversible dynamical evolution of an open quantum system C, coupled to its surrounding (bath) X2. The idea is to project the unitary evolution of the whole, system + bath, on the state space of the system alone. The bath, on which the thermodynamic limit has been taken, is assumed to be in a fixed state, not influenced by the interaction master equation thus obtained, which is Hamiltonian Hlz. The generalized formally exact, limits’,‘,3).

is then

simplified

by

some

Markovian

approximations

or

The Pauli master equation is deduced from the partial-trace projection. Note that it does not preserve the pure states. We present another example of the above scheme, in which the pure states are preserved. Such an example is interesting on several accounts. From a fundamental point of view we believe it is important to see by an explicit example that a dissipative quantum evolution is not necessarily stochastic. but may be deterministic in the sense of preserving the pure states4). Moreover, although we use exactly the same approximations and limits that lead to the Pauli master equation, we expect that far from equilibrium large differences exist between the different approaches (section 2). Next our proposal puts in a new perspective Sz-Nagy’s well-known theorem on dilations of contracting semi-groups (section 3). Last, but not least, our approach is a simple example of how a non-linear evolution can be deduced from the linear Schrodinger equation. *Partially

0378437

supported

by the Swiss National

I /82/0000-0000/$02.7S

@

Science

Foundation.

1982 North-Holland

MASTER

2. Formal deduction

EQUATION

FOR STATE VECTORS

365

and exact results

Let %‘, (resp. 2,) be the Hilbert space corresponding to the system (resp. to the bath). We suppose that to the whole, system and bath, is associated the tensor usual

product

X1 @ X2’). The Hamiltonian

of the whole

is assumed

of the

form: H = H,@1++@H~z+hH,z.

(I)

Let XOE 2’2 be the state of the bath which, according of master equations, we assume to be time-independent, E2,yo. As projection

to the general in particular

theory H2x0 =

we choose

P = I @ PO, where

(2)

PO is the projection

The theory of general equation of

onto the state x0. We identify master

equations

gives

X, with P - Xl @I X2.

the formally

exact

evolution

(pt = P e-‘“‘cPo@ x0.

(3)

It reads

+t = - t(H, + APH,,P)cp,

- A*

I

PHn(I

- P)

0 x

exp[ - i(H, + HZ+ h(PHnP

x (I - P)H,zPq,--s

- P)))s]

ds.

(4)

Going over to the interaction functions of the bath decrease approximation Cp,= - i(H, + APHr2P)~, where

+ (I - P)Hn(l

picture, and assuming sufficiently fast, one

that gets

the the

correlation Markovian

- A’(B + iA)cp,,

A and B are two linear

operators,

formally

(5) self-adjoint,

defined

by

D

B + iA =

I

Pei(H1+Hz’sH12(Z- P)e-i(H~+H~s(l - P)HIZP

ds.

(6)

0

The relationship between the eqs. (4) and (5) can be studied by the same techniques as for the usual Pauli master equation. In particular the singular coupling and the van Hove limit can be applied to the eq. (4). In the latter one lets the coupling constant A go to zero, and simultaneously one lets the time t go to infinity in such a way that tA* = T remains fixed; r is called the resealed time. A general theorem due to E.B. Davie8) gives the exact conditions under

N. GISIN

366

which the solution of the equation (4) tends, solution of (5) with the resealed time. With our choice decreasing function

of P, defined of time. But,

when

A goes

to zero,

to the

by (2), the norm of the vector cp, is a as long as cp,f 0, C,C+ defines a ray of the

Hilbert space X,, i.e. ‘p, defines a pure state of the system. In order to avoid possible difficulties with the limit t + 3~ (cpf could tend to zero), and to be consistent with the usual at each time:

The evolution

equation

+,t = -iH&

choice

of normalization,

of $, corresponding

+ AI((

we normalize

to (5) is easily

the vector

cp,

deduced:

- B)&,

(7)

where Ho = H, + APHlzP + A’A and (B),, = (C/JIB/J/~). The equation (7) defines a class of non-linear dissipative Schrodinger-like equations. We shall call them the Master equations for state vectors, since our formal deduction shows a close connection with the usual Pauli master equation. operators,

Using’) one easily proves that if Ho and then the eq. (7) has a unique solution:

B are

both

bounded

-,H”f-~h~RI rlrr =

(‘eiH,t-AtBr

(8)

. e-iHc,t?~Bf),4Jl/2’

Let us analyze this solution in the simple case discrete spectra. First, if Ho and B commute, then i(B),,

Ho and

this projection

is non-vanishing.

B have

(9)

= -2A2((B2)+,-(B);,)~0

and $cIItends as t + 130 to the (normalized) an eigensubspace of B. The corresponding which

where

projection of the initial state & on eigenvalue is the lowest one for Next,

if

Ho and B do not commute,

but if A is small enough, then the asymptotic state I,!J~ is close to an eigenvector C#Iof Ho for which (B), is minimal. As an example let us consider the evolution of one spin-i coupled to a bath of spins-$, i.e. a bath of two-level systems: X,=C?, The whole matrices:

X*=&C’. 1-I is supposed

to be in a magnetic

field,

so that,

with

u the Pauli

MASTER

EQUATION

FOR STATE

VECTORS

where al acts on the ith spin. The coupling Hamiltonian HI2 =

9 Ai(UlU++ da_),

367

is assumed to be

where (++= ux ‘2 ic+y

i=l

and the state of the bath is chosen as

Consequently Pe

PHIZP = 0 and

i(H,+Hz)sH,2,-i(HI+Hz)SH12p

In the thermodynamic limit (N + 00, n + 00, n/N = 8 = constant) the sums become integrals. If the distributions of the A: as functions of the oi are Lorentzian, then, as is well known, the correlation functions of the bath decrease exponentially and one gets m B+iA=h*

II

Oe- i(oo-0)s

e-s/r

c 2l -

(+z +

(1

_

f)),+i(w-~)s,-s/T

++,.p

0

-2 =

(coo-

:2+

l/7: rl

1 i + (i-e);+$wo-o) [

&+ (i- e)i(w,c I

0)

I

.I

err P, I

(10)

where W (resp. ,i) is the mean Larmor frequency (resp. mean coupling constant), and T, is the correlation time. According to (7) and (10) the Markovian master equation for the state vector of our one spin-f system reads:

(11) This equation has been studied in detail in*). The right-hand side of (11) introduces in the corresponding m = (u)$, a term, known as the Landau-Lifshitz friction used by several authors to describe the damping of a

second term on the evolution equation of term, which has been classical spin (see for

N. GISIN

368

example’)).

Note

that

(11) is stable

which are macroscopically

under

variations

undistinguishable.

From

of the state

of the bath

(11) one deduces:

(u;)ti, = th where

(12)

Co is an integration

ingly (a,)$, tends

constant

asymptotically

depending

on the initial

state

$0. Accord-

to 1 (resp.

- l), whenever 0 > 4 (resp. 0
eq. (11) gives a microscopic model of the wellBloch equationn). that if one uses an equation of the same type as (11) of a spin toward the instantaneous direction of a field B(t) = B. + B,(t), where B,(t) is a rotating field, is the so-called quasiperiodic state (or permanent the phenomenological theory of Redfield’<).

3. Inverse problem and dilations In this conditions

section we briefly study the following question: under can an evolution equation of the form of our eq. (7). i.e.

4, = -i HO+, + (WI),,

- HI)&.

with Ho and H, two self-adjoint operators, be associated evolution on a larger space in the sense of section 2?

which

(13) to an

unitary

A simple necessary and sufficient condition is that HI is bounded below. Indeed this condition is necessary since, as mentioned in the previous section, the norm of the solution of eq. (5) decreases with time. Consequently the operator B defined by (6) is necessarily positive. But it may contain a multiple of the identity, and this term cancels in eq. (7) due to the normalization. Thus the operator HI of (13) must be of the form HI = B - M - I, where B is a positive operator, and M a real number. This condition is also sufficient. Indeed if HI is bounded below, then the operator B defined by B = HI + M - I, where M is the lower bound of HI, is positive, and ((B)* - B)$ =

MASTER

EQUATION

FOR STATE VECTORS

369

((HI), - HI)+ for all +. Now exp{( -iHo- B)t} is a contraction semi-group and, according to Sz-Nagy’s theorem”), it can be dilated to a unitary evolution. Accordingly the dissipative evolution equation (13) with H, = k . Ho (k > 0) can be considered as describing the reduced dynamics of a system coupled to a heat bath if and only if Ho is bounded below. This is interesting since with this choice the equation (13) implies d(Ho)/dt = -2k(AHo)*, i.e. the mean value of Ho is a decreasing function of time for all states except the eigenstates. The eigenstates turn out to be asymptotic states, all semi-stable except for the ground-state which is stable (if it exist)“). One can apply this last result to the Hamiltonian of the harmonic oscillator and we then find an evolution which is analogous to the classical damped oscillator: If the initial state is the coherent state defined by the point (~0of the classical phase space, then the solution is at each time a coherent state and the corresponding function (Y,is the solution of the classical equations of the damped oscillator”).

4. Conclusion

We have shown that, by using the same approximations which are necessary to deduce the Pauli master equation, it is possible to construct a master equation for state vectors. The theorems which give the exact relationship between the dynamics of the reduced density matrix and the dynamics of the whole2v6), also apply to our approach. The non-linear term in our eq. (7) comes from the fact that one associates a ray (pure state) of the whole to a ray (pure state) of the system and that representating rays by normalized vectors (defined up to a global phase), this correspondence is not linear. The results of the simple examples of sections 2 and 3 are physically sound. Finally let us make two remarks. First we wish to emphasize that the master equation (4) is reversible, in spite of the appearance of an integral for which one often gives the antropomorphic interpretation of a memory2). Only the Markovian approximation describes an irreversible evolution. The latter approximation is usually applied to the Pauli master equation which governs the evolution of the density matrix of the system. But, and this leads to our second remark, one could a priori also think of making some approximation on the Schrodinger equation which governs the evolution of the state vector of the system. This could be done as presented in this paper, or, at least in a first step, by simply adding non-linear phenomenological terms to the Schrddinger equation, in the same way as one adds friction terms to the the different classical canonical equations. Due to these non-linearities

N. GISIN

370

microscopic towards

distributions

distributions

accordingly, equilibrium suggested

corresponding

which

no evolution however,

one

by the example

to one given

may correspond equation

for

expects

these

of section

2.

density

to different the

density

matrix

density matrix

non-linearities

evolve

matrices exist.

Close

to be irrelevant,

and, to as

Acknowledgement I wish to thank Profs E.B. Davies, for several enlightening discussions.

G.G. Emch,

Ph.A.

Martin

and C. Piron

References I) L. Prigogine and P. Resibois, Physica 27 (1961) 629. F. Haake, Springer Tracts in Modern Physics, Vol. 66, Springer. (Berlin, Heidelberg. New York, 1973). 2) V. Gorini et al., Rep. Math. Phys. 13(1978) 149. 3) Ph.A. Martin, Lecture Notes in Physics 103 (Springer, Berlin, 1979). 4) C. Piron, to be published in Essays in Honor of Wolfgang Yourgrau, Ed. Alwyn van der Merwe (Plenum, New York). 5) D. Aerts, The one and the many, these de doctorat. Universite Libre de Bruxelles. 1981. 6) E.B. Davies, Math. Ann. 219 (1976) 147. 7) M.S. Berger, Nonlinearity and Functional Analysis, Acad. Press. 1977 (Theorem 3.1.23). 8) N. &sin, Spin relaxation and dissipative Schrddinger-like evolution equations, to appear in Helv. Phys. Acta. 9) R. Kubo and N. Hashitsume, Prog. Th. Phys. Suppl. 46 (1970) 210. IO) B. Sz-Nagy and C. Foias. Harmonic Analysis of Operators on Hilbert Space (North-Holland. Amsterdam, 1970). 11) N. Gisin, J. Phys. Al4 (1981) 2259. 12) E.B. Davies, Quantum theory of an open systems (Acad. Press, Oxford, 1976) p. 154. 13) A. Abragam, The Principles of Nuclear Magnetism (Clarendon, Oxford, 1979). 14) R.H. Young and W.J. Deal, J. Math. Phys. ll(l970) 3298. See also G. Lochak and A. Alaoui. Ann. Foundation Louis de Broglie 2 (1977) 87. 15) A.G. Redfield, Phys. Rev. 96 (1955) 1787.