A non-linear evolution equation for the density operator is proposed which models the irreversible dissipative time evolution of a system in contact with its surroundings. Various properties of the proposed equation are discussed. As an illustration the three-state system is studied in detail.
We consider the temporal evolution of an isolated system of particles from a non-equilibrium state of entropy S = S′ to the equilibrium state of maximum entropy S = S max, S′ ≤ S max. The application of usual density matrix theory to the temporal development of S leads us to the unexpected result P = dS/dt = 0. The entropy does not change in time. This result is valid in case of the equilibrium state. However it is wrong for a non-equilibrium state, which could not irreversibly change from S′ to S max by entropy production P > 0. This contradicts the second law of thermodynamics. The paradoxial result can be traced back to the Liouville theorem. To overcome this difficulty, we introduce the t–P space. The time t and the entropy production P are observables in quantum thermodynamics. They are defined by the eigenvalue equations t|t〉 = t|t〉 and P|P〉 = P|P〉. The operators t and P do not commute. We have a non-vanishing commutator, [t,P] = ik, and consequently a t–P uncertainty relation, ΔtΔP ≥ k/2. The Boltzmann constant k is characteristic of quantum thermodynamics and represents the atomic entropy unit. — The time evolution operator U(τ) = exp[−(i/k)τP] shifts the time t by τ, U(τ)|t〉 = |t + τ〉. The generator of τ is the operator P of entropy production. Then the eigenvector |P〉 is also an eigenket of U(t), and it follows with progressing time τ ≥ 0 that τP ≥ 0 and hence P = dS/dt ≥ 0. This is the second law of thermodynamics, which is correctly reproduced in quantum thermodynamics. — Exchanging −τ for +τ replaces the equations of forward time +τ by the equations of backward time −τ. Independent of the direction of the arrow of time, the entropy of the system always increases or remains constant in the limit of the equilibrium state, dS ≥ 0. This is characteristic of irreversibility, for the system does not return to its initial state on time reversal. — The t–P uncertainty relation allows to discuss the temporal evolution of the expectation value 〈P〉 of entropy production and the lifetime Δt of the transient states, the system passes through on approaching the equilibrium state. In case of a non-equilibrium state 〈P〉, ΔP, and Δt are finite. On approaching the equilibrium state, 〈P〉 and ΔP decrease, consequently Δt increases. Finally, in the case of equilibrium (eq), we have 〈P〉eq = 0, ΔP eq = 0, and Δt eq → ∞. The lifetime of the equilibrium state is unlimited, because the entropy production 〈P〉 and the uncertainty of entropy production ΔP exactly vanish. The temporal evolution of the system can be described by means of a phenomenological equation.
In this paper, we outline the main features of the nonlinear quantum evolution equation proposed by the present author. Such an equation may be used as a model of reduced subsystem dynamics to complement various historical and contemporary efforts to extend linear Markovian theories of dissipative phenomena and relaxation based on master equations, Lindblad and Langevin equations, to the nonlinear and far nonequilibrium domain. It may also be used as the fundamental dynamical principle in theories that attempt to unite mechanics and thermodynamics, such as the Hatsopoulos–Gyftopoulos unified theory which motivated the original development of this well-behaved general nonlinear equation for the evolution of the density operator capable of generating irreversible deterministic relaxation to thermodynamic equilibrium from any far nonequilibrium state even for an isolated system.
In this contribution we prove that the entropy of an N-body isolated system can not decrease and the entropy production should be non-negative provided the system possesses an equilibrium state. We define the entropy as a functional of the set of n-particle reduced density operators (n=0,....,N) generalizing the von Neumann fine-grained entropy formula. Additionally, as a consequence of our analysis we find the expression of the equilibrium n-particle reduced density operators which enter the definition of the entropy as well as the dissipated energy in an irreversible process.
Abstract We consider the temporal evolution of an isolated system of N particles from a non-equilibrium state of entropy S = S′ to the equilibrium state of maximum entropy, S = S max, S′ ≤ S max. The application of usual density matrix theory to the temporal development of S leads us to dS/dt = 0: The entropy S does not change in time t. Thereby it is irrelevant, whether we consider a non-equilibrium state or an equilibrium state. Consequently, the system cannot irreversibly change by entropy production dS/dt > 0 from S′ to S max. This is a paradoxial result, which contradicts the experience. It can be traced back to the von Neumann equation, which in principle describes reversible processes and hence is unsuitable for calculating the irreversible evolution of the entropy S in time t. Each irreversible process is accompanied by a positive entropy production P = dS/dt ≥ 0 inside the system, which only vanishes in case of the equilibrium state. In order to overcome the above mentioned difficulties, we assign an operator P to the entropy production P, which is defined by the eigenvalue equation P |u〉 = P|u〉 with the state vector |u〉 of the N-particle system. There was an extensive discussion about the relation between the production of entropy P on one hand and the progression of time t on the other. Making use of this concept, we combine the operator of entropy production P with the time-development operator U (t, t 0) of the system and finally deduce the infinitesimal unitary operator U (t+τ,t) = 1 + (i/k)τ P by means of very general assumptions. Here P means the generator of the infinitesimal progression of time τ = t − t 0 and k the Boltzmann constant, representing the atomic entropy unit. Similarly to P we also treat the time t as an observable, defined by t |u〉 = t|u〉. We apply the infinitesimal time-evolution operator U (t+τ,t) to the operator of time t and finally obtain the P – t commutation relation i[ P , t ] = k, which is independent of τ. It shows that the operators P and t do not commute, and hence P and t are not sharply defined simultaneously. Instead we have uncertainties ΔP and Δt on measuring P and t, which are given by the P–t uncertainty relation ΔPΔt ≥ k/2. It readily allows a discussion of the evolution of the entropy S of the isolated system in time t from S′ to S max. Now, the irreversible steps are correctly described by the entropy production P = dS/dt > 0, while the thermal equilibrium is given by P = 0, ΔP = 0, and thus the lifetime of the equilibrium state Δt = ∞. According to the P–t uncertainty relation, the Boltzmann constant k is similarly important to the quantum thermodynamics of irreversible processes like Planck′s constant h to usual quantum mechanics.
We consider the temporal evolution of the entropy S of an isolated system from an initial non-equilibrium state 2· of entropy S = S i to the equilibrium state of maximum entropy S = S max > S i. The application of usual density matrix theory of quantum mechanics leads us to the surprising result P = dS/dt = 0. The entropy does not change in time. This result is valid in case of the equilibrium state. However it is wrong for a non-equilibrium state, because the entropy of the isolated system could not spontaneously increase from S i to S max by entropy production P > 0. This contradicts the second law of thermodynamics. We avoid this discrepancy by changing from quantum mechanics to quantum thermodynamics. Now the entropy production P = dS/dt and the time t are the relevant observables. In order to treat the temporal development of S and P, we at first prepare an initial state i of entropy production P i at time t < 0. Subsequently at ϑ t ≥ 0 the initial state i interacts with a continuum of final states f, which transfers the entropy production from i to f. The temporal evolution of the entropy production P(t) of our system is treated by means of the time-dependent perturbation theory. For that reason we split the operator of entropy production P = P 0 + W into two parts. The operator P 0 describes the unperturbed basic vectors of entropy production | i > and | f >, while the perturbation operator W determines the interaction between | i > and | f >. We expand the time-dependent state function < t | ψ > of the isolated system into the complete set of stationary eigenvectors | i > and | f > and calculate the time-dependent expansion coefficients c i(t) and c(f, t) by means of the dynamic equation of quantum thermodynamics − i k (d/dt) < t | ψ > = < t | P | ψ >. Here k means the Boltzmann constant, which is characteristic of quantum thermodynamics and represents the atomic entropy unit. The expansion coefficients are given by a system of coupled differential equations, from which we deduce an integral-differential equation for the coefficient c i(t) of the initial state i. Unfortunately it has been impossible to obtain an analytic solution of the complicated equation. However approximations are available. Accordingly, in the vicinity of the equilibrium state, the probability of still finding the initial state i of entropy production P i at time t is given by ω i(t) = | c i(t) |2 = exp (−Lt). The rate constant L ≥ σ depends on the density ρ(P f) of final states f of entropy production P f near by P i, and on the squared modulus of the perturbation matrix elements W if = < i | W | f >, which determines the interaction between | i > and | f >. The entropy production P(t) = P i ω i(t) = P i exp (−ϑLt) of the isolated system decreases exponentially and irreversibly with proceeding time t. Its rate-determining step is the transfer of entropy production from i to f. The integral over the entropy production of the final states f vanishes. Therefore the continuum states f do not contribute to the total entropy production P(t). The temporal evolution of the entropy S(t) can easily be calculated by integration of P(t). In conclusion, the entropy production of a non-equilibrium state is finite P(t) > 0 and only vanishes in the case of the equilibrium state P(∞) = 0. Thus, quantum thermodynamics agrees with the second law of thermodynamics.
We discuss a new analytical approach to real-time evolution in quantum many-body systems. Our approach extends the framework of continuous unitary transformations such that it amounts to a novel solution method for the Heisenberg equations of motion for an operator. It is our purpose to illustrate the accuracy of this approach by studying dissipative quantum systems on all timescales. In particular, we obtain results for non-equilibrium correlation functions for general initial conditions. We illustrate our ideas for the exactly solvable dissipative oscillator and, as a non-trivial model, for the dissipative two-state system.