Next Article in Journal
A Comprehensive Fault Diagnosis Method for Rolling Bearings Based on Refined Composite Multiscale Dispersion Entropy and Fast Ensemble Empirical Mode Decomposition
Next Article in Special Issue
Maximum Entropy Production Theorem for Transitions between Enzyme Functional States and Its Applications
Previous Article in Journal
Coexisting Attractors and Multistability in a Simple Memristive Wien-Bridge Chaotic Circuit
Previous Article in Special Issue
The Thermodynamics of Network Coding, and an Algorithmic Refinement of the Principle of Maximum Entropy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time–Energy and Time–Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations

by
Gian Paolo Beretta
Department of Mechanical and Industrial Engineering, Università di Brescia, via Branze 38, 25123 Brescia, Italy
Entropy 2019, 21(7), 679; https://doi.org/10.3390/e21070679
Submission received: 13 June 2019 / Revised: 7 July 2019 / Accepted: 8 July 2019 / Published: 11 July 2019
(This article belongs to the Special Issue Entropy Production and Its Applications: From Cosmology to Biology)

Abstract

:
In the domain of nondissipative unitary Hamiltonian dynamics, the well-known Mandelstam–Tamm–Messiah time–energy uncertainty relation τ F Δ H / 2 provides a general lower bound to the characteristic time τ F = Δ F / | d F / d t | with which the mean value of a generic quantum observable F can change with respect to the width Δ F of its uncertainty distribution (square root of F fluctuations). A useful practical consequence is that in unitary dynamics the states with longer lifetimes are those with smaller energy uncertainty Δ H (square root of energy fluctuations). Here we show that when unitary evolution is complemented with a steepest-entropy-ascent model of dissipation, the resulting nonlinear master equation entails that these lower bounds get modified and depend also on the entropy uncertainty Δ S (square root of entropy fluctuations). For example, we obtain the time–energy-and–time–entropy uncertainty relation ( 2 τ F Δ H / ) 2 + ( τ F Δ S / k B τ ) 2 1 where τ is a characteristic dissipation time functional that for each given state defines the strength of the nonunitary, steepest-entropy-ascent part of the assumed master equation. For purely dissipative dynamics this reduces to the time–entropy uncertainty relation τ F Δ S k B τ , meaning that the nonequilibrium dissipative states with longer lifetime are those with smaller entropy uncertainty Δ S .

1. Introduction

Recent advances in quantum information and quantum thermodynamics (QT) have increased the importance of estimating the lifetime of a given quantum state, for example to engineer decoherence correction protocols aimed at entanglement preservation. In the same spirit as fluctuation theorems that allow to estimate some statistical features of the dynamics from suitable state properties, also the Mandelstam–Tamm–Messiah time–energy uncertainty relations (MTM-TEURs) have been long known to provide bounds on lifetimes of quantum decaying states under Hamiltonian (non-dissipative) evolution. For practical applications, however, such bounds are insufficient when Hamiltonian dynamics must be complemented by models of dissipation and decoherence.
The time–energy uncertainty relation has remained an open and at times controversial issue throughout the history of quantum theory. Several reviews are available on the pioneering discussions and the subsequent developments [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]. In the present paper, we are motivated by the past two decades of important advancements in our understanding of the general structure of dynamical models for non-equilibrium thermodynamics, including non-equilibrium quantum thermodynamic models. Such revival has been prompted and paralleled by a steady advancement of experimental techniques dealing with single ion traps [16,17], qubits [18,19], neutron interferometry [20,21], and a countless and growing number of other quantum-information developments since then, e.g., nonlinear quantum metrology [22,23]. Within these applications, TEURs can provide useful information and practical bounds for parameter estimation. But since dissipation and decoherence are often the limiting factors, there is a need to generalize the MTM-TEURs to frameworks where microscopic few-particle quantum setups exhibit non-unitary dissipative dynamical behavior.
A recent review paper on the physical significance of TEURs provides 300 references and the following conclusion [24]: “We have shown that the area of energy–time uncertainty relations continues to attract attention of many researchers until now, and it remains alive almost 90 years after its birth. It received a new breath in the past quarter of century due to the actual problems of quantum information theory and impressive progress of the experimental technique in quantum optics and atomic physics. It is impossible to describe various applications of the TEURs to numerous different physical phenomena in this minireview.”
The main objective in the present paper is to extend the time–energy uncertainty relations to the framework of dissipative quantum dynamical systems. But differently from the most popular and traditional model of dissipation in open quantum systems, which is based on the well-known Kossakowski–Lindblad–Gorini–Sudarshan (KLGS) master equations [25,26,27,28,29,30,31,32,33], we assume the less known locally steepest-entropy-ascent (LSEA) model of dissipation. We make this choice not only to avoid some drawbacks (outlined in more details in Appendix A) of the KLGS master equation from the point of view of full and strong consistency with the general principles of thermodynamics, causality, and far non-equilibrium, but more importantly because we have shown in References [34,35] that the LSEA principle—by providing the minimal but essential elements of thermodynamic consistency, near as well as far from stable (maximal entropy) equilibrium states—has the potential to unify all the successful frameworks of non-equilibrium modeling, from kinetic theory to chemical kinetics, from stochastic to mesoscopic to extended irreversible thermodynamics, as well as the metriplectic structure or, in more recent terms, the General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) structure. In addition, it is noteworthy that a particular but broad class of KLGS master equations has been recently shown to fall into a LSEA (entropic gradient flow) structure [36,37], and hence some of the TEURs we derive here hold also for such class of models.
Steepest-entropy-ascent (SEA) nonlinear master equations have proved to be effective tools to model dissipative dynamics, thermalization, transport processes and, in general, entropy production in a wide range of frameworks of non-equilibrium thermodynamics. In essence, SEA models are explicit implementations of the general principle of maximum local entropy production. In recent mathematical terms, SEA models are entropic gradient flows. From the fundamental point of view, the general structure and nonlinearity of the SEA master equations are instrumental to providing strong compatibility with the second law of thermodynamics by guaranteeing, within the model, the Hatsopoulos–Keenan statement of existence and uniqueness of the stable equilibrium (maximum entropy) states. Here we focus on the quantum thermodynamic modeling framework of application and show how the entropy production modifies the usual TEURs.
The usual time–energy uncertainty relation—as interpreted according to the Mandelstam–Tamm– Messiah intrinsic-time approach [38,39] based on unitary Hamiltonian dynamics—is modified by the presence of a maximally dissipative term in the dynamical law, which models at the single- or few-particle quantum level the so-called maximum entropy production principle (MEPP) [40,41,42,43,44,45,46,47]. TEURs obtained in other frameworks [48,49,50,51,52,53,54,55,56] such as attempts to define time or “tempus” operators, entropic uncertainties, and measurement times are beyond our scope here.
The class of MEPP master equations we designed in References [57,58,59,60,61,62] is suitable to model dissipation phenomenologically not only in open quantum systems in contact with macroscopic baths, but also in closed isolated systems, as well as strongly coupled and entangled composite systems (references below). These master equations are capable to describe the natural tendency of any initial nonequilibrium state (read: density operator) to relax towards canonical or partially-canonical thermodynamic equilibrium (Gibbs state), i.e., capable of describing the irreversible tendency to evolve towards the highest entropy state compatible with the instantaneous mean values of the energy (and possibly other constants of the motion and other constraints). They do so by preserving exactly the conserved properties while pulling every nonequilibrium state in the SEA direction with respect to the local dissipation metric that is part of the nonequilibrium description of the system [34]. This dissipative tendency is simultaneous and in competition with the usual non-dissipative Hamiltonian unitary evolution.
Our original approach—when understood as an attempt to develop thermodynamically consistent modeling approaches that merge mechanics and thermodynamics following Hatsopoulos and Gyftopoulos [63,64,65,66]—can perhaps be considered a first pioneering “resource theory” of quantum thermodynamics equipped with a nonlinear dissipative dynamical structure capable to describe relaxation even from arbitrarily far from equilibrium and to entail the second law as a theorem of the dynamical law. Several other pioneering aspects of QT resource theories were present in References [63,64,65,66]. For example, the energy versus entropy diagram to represent nonequilibrium states in the QT framework, first introduced in Reference [65], has recently found interesting applications in [67]. Again, it provides definitions and expressions for adiabatic availability and available energy with respect to a heat bath, work element, heat interaction, etc. which are currently discussed intensely in the QT community. It must also be mentioned that this first QT resource theory was proposed in years when talking of quantum thermodynamics was considered heresy by the orthodox physical community. Considering that it is little cited and still not well known, we give more details below and in Appendix B.
We provide in the two appendices a brief review of some practical and conceptual issues of the prevailing model of irreversibility, and a discussion of the original motivation that lead us to develop a quantum maximal entropy production formalism. We do not repeat here the geometrical derivations of our nonlinear MEPP dynamical law, nor the discussions of its many intriguing mathematical–physics implications, because they are available in many previous papers. Here, we simply adopt that master equation without derivation, and focus on its consequences related to TEURs, illustrated also by some numerical simulations. Quantum statistical mechanics and quantum thermodynamics practitioners have so far essentially dismissed and ignored our class of SEA master equations on the basis that they do not belong to the standard class of KLGS master equations and hence cannot be the correct description of the reduced dynamics of a system in interaction with one or more thermal baths. However, at least when used as phenomenological modeling tools, SEA master equations have recently proved [68,69,70,71,72,73,74,75,76] to offer in a variety of fields important advantages of broader or complementary applicability for the description and correlation of near- and far-non-equilibrium behavior.
In the quantum framework, the local state of a subsystem is represented by the local density operator ρ and its lifetime may be characterized by the intrinsic characteristic times τ F of the dynamical variables associated with the linear functionals Tr ( ρ F ) . If the local dynamics is non-dissipative and described by the usual unitary evolution, we show below that the Heisenberg–Robertson inequality entails the usual MTM-TEURs τ F Δ H / 2 , while the Schroedinger inequality entails sharper and more general exact TEURs [Equation (23)].
For simultaneous unitary+dissipative dynamics, the usual TEUR is expectedly replaced by less restrictive relations and additional characteristic times acquire physical significance. In particular, we focus our attention to the characteristic time associated with the rate of change of the von Neumann entropy functional k B Tr ( ρ ln ρ ) . For unitary+LSEA evolution, in Section 7 [Equation (83)] we obtain an interesting time–energy and time–entropy uncertainty relation ( 2 τ F Δ H / ) 2 + ( τ F Δ S / k B τ ) 2 1 where τ is the main dissipation time that defines the strength of the dissipative component of the assumed dynamical law. With the help of numerical simulations, we illustrate this relation and several other even more precise uncertainty relations, that in the framework of QT resource theories may have a useful application in quantifying the lifetime of quantum states.
The structure of the paper is outlined at the end of the next section, where we first introduce the particular class of nonlinear dissipative quantum master equations on which we restrict our attention in the first part of the paper.

2. Assumed Structure of the Nonlinear Dissipative Quantum Master Equation

Let H ( dim H ) be the Hilbert space and H the Hamiltonian operator that in standard Quantum Mechanics we associate with a given isolated (or adiabatic, see below) and uncorrelated system. We assume that the quantum states are one-to-one with the linear hermitian operators ρ on H with Tr ( ρ ) = 1 and ρ ρ 2 , and we assume a dynamical equation of the form
d ρ d t = ρ E ( ρ ) + E ( ρ ) ρ ,
where E ( ρ ) is an operator-valued function of ρ that we may call the “evolution generator” which may in general be non-hermitian and nonlinear in ρ , but must be such as to preserve ρ unit trace and non-negative definite. Without loss of generality, we write E = E + + i E where E + = ( E + E ) / 2 and E = ( E E ) / 2 i are hermitian operators. Then, the dynamical law takes the form
d ρ d t = i [ E , ρ ] + { E + , ρ } ,
where [ · , · ] and { · , · } are the usual commutator and anticommutator, respectively. In Appendix A we discuss the reasons why we adopt this form, and exclude terms like V ( ρ ) ρ V ( ρ ) which appear instead in the the celebrated KSGL class of (linear) quantum master equations.
In preparation for our SEA construction in Section 7, we assume E = H / (independent of ρ ), where H is the Hamiltonian operator of the system and the reduced Planck constant, and rewrite E + as E + = Δ M ( ρ ) / 2 k B τ where k B is the Boltzmann constant, τ a positive constant (or state functional) that in the SEA framework we will interpret as an intrinsic dissipation time of the system, because it essentially fixes the rate at which the state evolves along the path of SEA in state space, and Δ M ( ρ ) a hermitian operator-valued nonlinear function of ρ that we call the “nonequilibrium Massieu operator” and until Section 7 we do not define explicitly, except for the assumption that it satisfies the condition
Tr [ ρ Δ M ( ρ ) ] = 0
as well as the condition that it preserves the nonnegativity of ρ (both forward and backwards in time!). As a result, Equation (1) takes the form
d ρ d t = i [ H , ρ ] + 1 2 k B τ { Δ M ( ρ ) , ρ } .
Let us note that as in standard unitary dynamics, we say that the particle is either isolated or adiabatic, respectively, if the Hamiltonian operator H is either time independent or time dependent, for example, through one or more external control parameters.
In Section 7, we will consider for Δ M ( ρ ) the explicit SEA form for the simplest case, first proposed in [57,58,61]. Reference [59] proposed also a general LSEA form for a composite quantum system, which will not be considered here, but has clear applications in the description of decoherence and lifetime of entanglement (see Reference [68]).
In the present paper we limit the discussion to the derivation of general inequalities, and to illustrative considerations and a numerical example valid within the simplest framework of steepest-entropy-ascent conservative dynamics. The application to structured composite systems based on our LSEA version [59,68] of operator Δ M ( ρ ) will be discussed elsewhere.
The specific physical interpretations of the uncertainty relations that follow from dynamical law (4) will depend on the theoretical or modeling context in which such time evolution is assumed. For example, the problem of designing well-behaved nonlinear extensions of the standard unitary dynamical law of quantum mechanics has been faced in the past few decades with a variety of motivations, and is recently seeing a vigorous revival in connection with questions about the foundations of quantum mechanics and the need for thermodynamically sound phenomenological models (recently referred to as “resource theories” [67]) that arise from the current developments of quantum information technologies and related single-particle and single-photon experiments to test quantum computing components and devices and fundamental questions about entanglement, decoherence, nonlocality, and measurement theory.
In our original development [57,58,59,60,61], Equation (4) was designed as part of an ad-hoc fundamental dynamical postulate needed [77,78,79,80,81,82] to complete the Hatsopoulos-Gyftopoulos attempt [63,64,65,66,83,84] to unify mechanics and thermodynamics into a generalized quantum theory by building the Hatsopoulos–Keenan statement of the second law [85,86] directly into the microscopic level of description. In particular, the key ansatz in References [63,64,65,66] is the assumption that even for a strictly isolated system, there exists a broad class of genuine states (homogeneous preparations, in von Neumann language [57,82,87,88]) that require non-idempotent density operators, i.e., such that ρ 2 ρ . Two decades later, this ansatz has been re-proposed in Reference [89], and our nonlinear dynamical Equation (4) has been re-discovered and studied in References [90,91,92,93], where it is shown to be well-behaved from various perspectives including a relativistic point of view. An important feature is that it entails non-unitary evolution only for non-idempotent ( ρ 2 ρ ) density operators, whereas for idempotent ( ρ 2 = ρ ) density operators it entails the standard unitary evolution (see, e.g., References [58,94]).
However, the present results are valid also in any other framework, theoretical discussion, modeling context, or resource theory whereby—for example to study decoherence, dissipation, quantum thermal engines, quantum refrigerators, and so on—the usual Liouville-von Neumann equation for the density operator is modified, linearly or nonlinearly, into form (4).
Since many of the relations we derive here are valid and nontrivial in all these contexts, in Section 3, Section 4, Section 5 and Section 6 we begin by presenting the results that do not depend on assuming a particular form of operator Δ M ( ρ ) . Thus, independently of the interpretation, the context of application, and the specific form of master Equation (4), the uncertainty relations derived in the first part of the paper extend the usual relations to the far non-equilibrium domain and in general to all non-zero-entropy states.
In Section 7 and Section 8, to fix ideas and be able to present numerical results and qualitative considerations, we specialize the analysis to the simplest nontrivial form of Equation (4) that implements our conservative steepest-entropy-ascent dynamical ansatz, namely, a model for irreversible relaxation of a four-level qudit.
Appendix A discusses our reasons for not considering, in the present context, the extension of our results to a full Kossakowski-Lindblad form of the evolution equation.
Appendix B gives a brief review of the original motivations that lead us to develop the SEA and LSEA formalism in the early quantum thermodynamics scenario, and of the subsequent developments that in recent years have shown how the locally steepest-entropy-ascent principle not only gives a clear, explicit, and unambiguous meaning to the MEPP but it also constitutes the heart of (and essentially unifies) all successful theories of nonequilibrium.

3. General Uncertainty Relations

We consider the space L ( H ) of linear operators on H equipped with the real scalar product
( F | G ) = Tr ( F G + G F ) / 2 = ( G | F ) ,
and the real antisymmetric bilinear form
( F \ G ) = i Tr ( F G G F ) / 2 = ( G \ F ) = ( F | i G ) ,
so that for any hermitian F in L ( H ) the corresponding mean-value state functional can be written as F = Tr ( ρ F ) = Tr ( ρ F ρ ) = ( ρ | ρ F ) , and can therefore be viewed as a functional of ρ , the square-root density operator, obtained from the spectral expansion of ρ by substituting its eigenvalues with their positive square roots. When ρ evolves according to Equation (1) and F is time-independent, the rate of change of Tr ( ρ F ) can be written as
d Tr ( ρ F ) / d t = Tr ( F d ρ / d t ) = 2 ρ F ρ E ( ρ ) .
In particular, for the evolution Equation (1) to be well defined, the functional Tr ( ρ I ) where I is the identity on H must remain equal to unity at all times; therefore, d Tr ( ρ I ) / d t = 2 ρ I ρ E ( ρ ) = 0 or, equivalently, in view of Equation (4), Equation (3) rewrites as
ρ ρ Δ M ( ρ ) = 0 .
For F and G hermitian in L ( H ) , we introduce the following shorthand notation
Δ F = F Tr ( ρ F ) I ,
σ F G = Δ F Δ G = ( ρ Δ F | ρ Δ G ) = 1 2 Tr ( ρ { Δ F , Δ G } ) = σ G F ,
Δ F = σ F F = Δ F Δ F ,
η F G = [ F , G ] / 2 i = ( ρ Δ F \ ρ Δ G ) = 1 2 i Tr ( ρ [ F , G ] ) = η F G * = η G F ,
For example, we may write the rate of change of the mean value of a time-independent observable F as
d Tr ( ρ F ) d t = [ F , H ] / 2 i / 2 + Δ F Δ M k B τ = η F H / 2 + σ F M k B τ ,
from which we see that not all operators F that commute with H correspond to constants of the motion, but only those for which Δ F Δ M = 0 , i.e., such that ρ Δ F is orthogonal to both i ρ Δ H and ρ Δ M , in the sense of scalar product (5). For an isolated system, conservation of the mean energy functional Tr ( ρ H ) requires an operator function Δ M ( ρ ) that maintains ρ Δ M always orthogonal to ρ Δ H , so that Δ H Δ M = 0 for every ρ .
From the Schwarz inequality, we readily verify the following generalized Schrödinger uncertainty relation
Δ F Δ F Δ G Δ G Δ F Δ G 2 + [ F , G ] / 2 i 2 ,
usually written in the form σ F F σ G G σ F G 2 | η F G | . It follows from the Cauchy–Schwarz inequality ( f , f ) ( g , g ) | ( f , g ) | 2 and the identity | ( f , g ) | 2 = ( f | g ) 2 + ( f \ g ) 2 where ( f | g ) = [ ( f , g ) + ( g , f ) ] / 2 , ( f \ g ) = i [ ( f , g ) ( g , f ) ] / 2 , and f, g are vectors in some complex Hilbert space (strict equality iff f = λ g for some scalar λ ). In the space L ( H ) of linear operators on H equipped with the complex scalar product ( f , g ) = Tr ( f g ) , we note that ( f , f ) = ( f | f ) and obtain the inequality ( f | f ) ( g | g ) ( f | g ) 2 + ( f \ g ) 2 and hence inequality (14) by setting f = ρ Δ F and g = ρ Δ G . Note that the strict equality in (14) holds iff ρ Δ F = λ ρ Δ G for some scalar λ (in which case we have [ F , G ] / 2 i = 0 iff either λ * = λ or ρ Δ F = 0 or both). This proof was given as footnote 7 of Reference [95]. For Schroedinger’s original proof and an alternative one see Reference [96]. Relation (14) is a generalization of the inequality first appeared in [97,98] and later generalized in [99] to the form det σ = σ F F σ G G σ F G 2 det η = η F G 2 = η G F 2 , suitable for generalizations to more than two observables. Early proofs of relation (14) were restricted to pure state operators ( ρ 2 = ρ ). To our knowledge, the earliest proof valid for general (mixed and pure) states ρ is that in [6]. For further inequalities in the case of position and momentum operators see [14] and references therein. Notice also that by using our proof of the Schrodinger inequality (14), just given above, Relation (22) of the main theorem in the review paper [15] can be made sharper and read | Tr ( R [ A , B ] ) | 2 + | Tr ( R { A , B } ) 2 n λ n Tr ( P n A P n B ) | 2 4 f 2 ( R , A ) f 2 ( R , B ) .
Relation (14) obviously entails the less precise and less symmetric Heisenberg-Robertson uncertainty relation
Δ F Δ F Δ G Δ G [ F , G ] / 2 i 2 ,
usually written in the form Δ F Δ G | η F G | .
For further compactness, we introduce the notation
r F G = σ F G / σ F F σ G G , c F G = η F G / σ F F σ G G ,
where clearly, r F G represents the cosine of the angle between the ‘vectors’ ρ Δ F and ρ Δ G in L ( H ) , and r F G 2 1 . Inequality (14) may thus be rewritten as
r F G 2 + c F G 2 1
and clearly implies
c F G 2 1 1 + ( r F G 2 / c F G 2 ) 1 r F G 2 1 .
Next, for any hermitian F we define the characteristic time of change of the corresponding property defined by the mean value of the linear functional F = Tr ( ρ F ) as follows
τ F ( ρ ) = Δ F / | d F / d t | .
As is well known [1,2,3,4,5,7,8,10,11,14,15,38,39,48], τ F represents the time required for the statistical distribution of measurements of observable F to be appreciably modified, i.e., for the mean value F to change by an amount equal to the width Δ F of the distribution.
Now, defining the nonnegative, dimensionless functional
a τ = Δ M / 2 k B τ Δ H ,
we rewrite (13) in the form
d F / d t = 2 Δ F Δ H ( c F H + a τ r F M ) /
and, substituting into (19), we obtain the general exact uncertainty relation
/ 2 τ F Δ H = | c F H + a τ r F M | .
For non-dissipative dynamics [ Δ M ( ρ ) / τ = 0 ], a τ = 0 , Equation (22) yields the time–energy uncertainty relations
2 / 4 τ F 2 σ H H = c F H 2 1 1 + ( r F H 2 / c F H 2 ) 1 r F H 2 1 ,
which entail but are more precise than the usual time–energy uncertainty relation, in the same sense as Schrödinger’s relation (14) entails but is more precise than Heisenberg’s relation (15). According to (19), the last inequality in (23) implies that property F cannot change at rates faster than 2 Δ F Δ H / .
For dissipative dynamics let us first consider an observable A that commutes with H, so that [ A , H ] / 2 i = 0 while Δ A Δ H 0 ; in other words, an observable conserved by the Hamiltonian term in the dynamical law (4), but not conserved by the dissipative term. Then Equation (22) yields the equivalent time–energy uncertainty relations
/ 2 τ A Δ H = a τ | r A M | a τ ,
k B τ τ A Δ M = | r A M | 1 .
We note that while r A M 2 1 , the value of a τ depends on how Δ M ( ρ ) / τ is defined and, a priori, could well be larger than unity, in which case there could be some observables A for which τ A Δ H / 2 . If instead we impose that the operator function Δ M ( ρ ) / τ is defined in such a way that a τ 1 , i.e.,
τ Δ M / 2 k B Δ H ,
then we obtain that even in dissipative dynamics the usual time–energy uncertainty relations are never violated by observables A commuting with H. In Section 8 we will consider a numerical example for a case with non-constant τ given by Equation (26) with strict equality, for a qualitative comparison with the same case with constant τ .
However, in general, if the dynamics is dissipative [ Δ M ( ρ ) / τ 0 ] there are density operators for which | c F H + a τ r F M | > 1 so that τ F Δ H takes a value less than / 2 and thus the usual time–energy uncertainty relation is violated. The sharpest general time–energy uncertainty relation that is always satisfied when both Hamiltonian and dissipative dynamics are active is (proof in Section 5)
2 / 4 τ F 2 σ H H 1 + a τ 2 + 2 a τ c M H ,
which may also take the equivalent form
τ F 2 σ H H 2 / 4 + τ F 2 σ M M k B 2 τ 2 ( ρ ) + τ F 2 Δ M Δ H c M H k B τ / 4 1 .
The upper bound in the rate of change of property F becomes
Δ F σ H H 2 / 4 + σ M M k B 2 τ + Δ M Δ H c M H k B τ / 4 .
As anticipated, because the dissipative term in Equation (4) implies an additional dynamical mechanism, this bound (29), valid for the particular nonunitary dynamics we are considering, is higher than the standard bound valid in unitary hamiltonian dynamics, given by 2 Δ F Δ H / . For observables commuting with H, however, (25) provides the sharper general bound Δ F Δ M / k B τ , solely due to dissipative dynamics, which is lower than (29).
Because in general | c M H | < 1 , (28) obviously implies the less precise relation
2 / 4 τ F 2 σ H H ( 1 + a τ ) 2 .
However, as for the dynamics we discuss in Section 7, if the Massieu operator Δ M ( ρ ) is a linear combination (with coefficients that may depend nonlinearly on ρ ) of operators that commute with either ρ or H, then it is easy to show that c M H = 0 . Therefore, in such important case, (28) becomes
2 / 4 τ F 2 σ H H 1 + a τ 2 ,
clearly sharper than (30). If in addition Δ M ( ρ ) / τ is such that (26) is satisfied, then (31) implies τ F Δ H / 2 2 .

4. Characteristic Time of the Rate of Entropy Change

We now consider the entropy functional S = Tr ( ρ S ) = k B Tr ( ρ ln ρ ) = k B ρ ρ ln ( ρ ) 2 and its rate of change, which using Equations (4) and (8) may be written as
d Tr ( ρ S ) / d t = 2 ρ S ρ E ( ρ ) = Δ S Δ M / k B τ = Δ S Δ M r S M / k B τ ,
where S is the entropy operator defined as follows
S = k B P Ran ρ ln ρ = k B ln ( ρ + P Ker ρ ) ,
where P Ran ρ and P Ker ρ are the projection operators onto the range and kernel of ρ . Operator S, introduced in [58,62], is always well defined for any ρ ρ 2 , even if some eigenvalues of ρ are zero. It is the null operator when ρ 2 = ρ . In models where S is always multiplied by ρ or ρ , the operators P Ran ρ (or P Ker ρ ) in Equation (33) could be omitted, because in general ρ S = k B ρ ln ρ and ρ S = k B ρ ln ρ . But, in models of decoherence and composite systems based on the LSEA equation of motion proposed in [59], further discussed in [100], and applied for example in [68,69], their role is important because the LSEA master evolution equation involves the operators
( H ) J = Tr J ¯ [ ( I J ρ J ¯ ) H ] ,
( S ) J = Tr J ¯ [ ( I J ρ J ¯ ) S ] ,
that we call “locally perceived overall-system energy operator” and “locally perceived overall-system entropy operator,” respectively, associated with a mean-field-like measure of how the overall-system energy and entropy operators, H and S, are “perceived” locally within the J-th constituent subsystem. The symbol J ¯ denotes the composite of all subsystems except the J -th one. As discussed in full details in [59,100] the dissipative term in our LSEA master equation points in the direction of the local constrained gradient of the “locally perceived overall-system entropy” TrJ[ρJ(S)J], constrained by the condition of orthogonality with respect to the local gradient of the “locally perceived overall-system energy” TrJ[ρJ(H)J]. Operators (S)J and, hence, the LSEA models just mentioned, would not be well defined without P Ran ρ (or P Ker ρ ) in Equation (33).
Interestingly, the rate of entropy change, being proportional to the correlation coefficient between entropy measurements and M measurements, under the assumptions made so far, may be positive or negative, depending on how Δ M ( ρ ) is defined, i.e., depending on the specifics of the physical model in which Equation (4) is adopted.
The characteristic time of change of the entropy functional, defined as
τ S = Δ S / | d S / d t | ,
gives rise to the following equivalent exact time–energy uncertainty relations
/ 2 τ S Δ H = a τ | r S M | a τ ,
k B τ τ S Δ M = | r S M | 1 ,
where r S M is defined as in (16) using operators Δ M ( ρ ) and Δ S = S S . The physical interpretation of (38) is that the entropy cannot change in time at a rate faster than Δ S Δ M / k B τ , as immediately obvious also from (32).
We notice from (37) that if the nonequilibrium Massieu operator satisfies condition (26) then a τ 1 and, therefore, the characteristic time of entropy change, τ S , satisfies the usual uncertainty relation τ S Δ H / 2 and the rate of entropy change cannot exceed 2 Δ S Δ H / .
We conclude this Section by noting that, in general, the equality in (37) may be used to rewrite Relation (27) in the form
a τ 1 + a τ | r S M | τ S τ F 1 + a τ 2 + 2 a τ c M H 1 + a τ τ F ,
where the last inequality follows from | c M H | 1 . This relation shows, on one hand, that the entropy change characteristic time τ S is not necessarily the shortest among the characteristic times τ F associated with observables of the type F = Tr ( ρ F ) according to the Mandelstam–Tamm definition (19). On the other hand, it also shows that the left-hand side defines a characteristic-time functional
τ U D = a τ 1 + a τ | r S M | τ S τ F ,
which constitutes a general lower bound for all τ F ’s, and may therefore be considered the shortest characteristic time of simultaneous unitary+dissipative dynamics as described by Equation (4). This observation prompts the discussion in the next section.

5. Shortest Characteristic Times for Purely-Unitary and Purely-Dissipative Dynamics

The Mandelstam–Tamm definition (19) of characteristic times has been criticized for various reasons (see for example References [101,102,103]) mainly related to the fact that depending on which observable F is investigated, as seen by inspecting (23), the bound τ F / 2 Δ H may be very poor whenever c F H 2 is much smaller than 1.
Therefore, different attempts have been made to define characteristic times that (1) refer to the quantum system as a whole rather than to some particular observable, and (2) bound all the particular τ F ’s from below. Notable examples are the characteristic times τ E S and τ L K , respectively defined by Eberly and Singh [101] and Leubner and Kiener [102].
Here, however, we consider the shortest characteristic times that emerge from the following geometrical observations. The functional Δ F may be interpreted as the norm of ρ Δ F (viewed as a vector in L ( H ) ) in the sense that it equals ( ρ Δ F | ρ Δ F ) , therefore, we may use it to define the (generally non hermitian) unit norm vector in L ( H )
F ˜ ρ = ρ Δ F / Δ F .
As a result, Equation (13) may be rewritten in the form
1 Δ F d F d t = Δ H / 2 ( F ˜ ρ | i H ˜ ρ ) + Δ M k B τ ( F ˜ ρ | M ˜ ρ ) = ( F ˜ ρ | C ) ,
where for shorthand we define the operator
C = i Δ H H ˜ ρ / 2 + Δ M M ˜ ρ k B τ = 2 ρ E ( ρ ) ,
directly related [see Equation (7)] with the evolution operator function E ( ρ ) defined in Section 2, which determines the rates of change of all linear functionals of the state operator ρ , i.e., all observables of the linear type Tr ( ρ F ) , by its projection onto the respective directions F ˜ ρ .
Each characteristic time τ F can now be written as
τ F = Δ F / | d F / d t | = 1 / | ( F ˜ ρ | C ) | .
Because F ˜ ρ is unit norm, | ( F ˜ ρ | C ) | is bounded by the value attained for an operator F ˜ ρ that has the same ‘direction’ in L ( H ) as operator C, i.e., for
F ˜ ρ = ± C / ( C | C ) ,
in which case | ( F ˜ ρ | C ) | = ( C | C ) = Tr ( C C ) . Thus we conclude that, for any, F,
1 / ( C | C ) τ F ,
and, therefore, we introduce the shortest characteristic time for the combined unitary+dissipative dynamics described by Equation (4),
τ U D = 1 / ( C | C ) ,
which binds from below all τ F ’s. From (43) and (46), and the identities ( i H ˜ ρ | i H ˜ ρ ) = ( M ˜ ρ | M ˜ ρ ) = 1 and ( i H ˜ ρ | M ˜ ρ ) = ( M ˜ ρ | i H ˜ ρ ) = c M H we obtain
1 τ F 2 1 τ U D 2 = ( C | C ) = σ H H 2 / 4 + σ M M k B 2 τ 2 ( ρ ) + Δ M Δ H c M H k B τ / 4 = σ H H 2 / 4 ( 1 + a τ 2 + 2 a τ c M H ) ,
which proves relations (27) and (28).
For nondissipative (purely Hamiltonian, unitary) dynamics the same reasoning (or substitution of τ = , a τ = 0 in the above relations) leads to the definition of the shortest characteristic time of unitary dynamics
τ U = / 2 Δ H ,
with which the usual time–energy relation reduces to
τ F τ U .
Its physical meaning is that when the energy dispersion (or uncertainty or spread) Δ H is small, τ U is large and τ F must be larger for all observables F, therefore, the mean values of all properties change slowly [15,104,105], i.e., the state ρ has a long lifetime. In other words, states with a small energy spread cannot change rapidly with time. Conversely, states that change rapidly due to unitary dynamics, necessarily have a large energy spread.
Another interesting extreme case obtains from Equation (4) when Δ M ( ρ ) is such that the condition [ ρ , H ] = 0 implies [ Δ M ( ρ ) , H ] = 0 for any ρ , as for the steepest-entropy-ascent dynamics discussed in Section 7 and Section 8. In this case, it is easy to see that if the state operator ρ commutes with H at one instant of time then it commutes with H at all times and, therefore, the entire time evolution is purely dissipative. Then, the reasoning above leads to the definition of the shortest characteristic time of purely dissipative evolution
τ D = k B τ / Δ M .
It is noteworthy that τ D can be viewed as the characteristic time associated not with the (generally nonlinear) Massieu functional M = Tr ( ρ M ( ρ ) ) but with the linear functional A = Tr ( ρ A ) corresponding to the time-independent operator A which at time t happens to coincide with M ( ρ ( t ) ) .
For purely dissipative dynamics, the bound τ F τ D = k B τ / Δ M implies that when Δ M / k B τ , i.e., the ratio between the uncertainty in our generalized nonequilibrium Massieu observable represented by operator M and the intrinsic dissipation time τ , is small, then τ D is large and τ F must be larger for all observables F, therefore, the state ρ has a long lifetime. This may be a desirable feature in quantum computing applications where the interest is in engineering states ρ that preserve the entanglement of component subsystems. Conversely, if some observable changes rapidly, τ F is small and since τ D must be smaller, we conclude that the spread Δ M (more precisely, the ratio Δ M / k B τ ) must be large.
In terms of τ U and τ D we can rewrite (20), (38) and (48) as
a τ = τ U / τ D ,
1 τ S = | r S M | τ D 1 τ D ,
1 τ F 2 = c F H τ U + r F M τ D 2 1 τ U D 2 = 1 τ U 2 + 1 τ D 2 + 2 c M H τ U τ D 1 τ U + 1 τ D 2 .
Equation (53) implies that the entropy cannot change rapidly with time if the ratio Δ M ( ρ ) / k B τ is not large. The first equality in (54) follows from ( F ˜ ρ | i H ˜ ρ ) = c F H and ( F ˜ ρ | M ˜ ρ ) = r F M , which also imply that Equation (42) may take the form
d F d t = Δ F c F H τ U + r F M τ D ,
and operator C defined in (43) takes also the forms
C = i H ˜ ρ τ U + M ˜ ρ τ D = i ρ Δ H Δ H τ U + ρ Δ M Δ M τ D ,
and its norm is 1 / τ U 2 + 1 / τ D 2 + 2 c M H / τ U τ D .
Similarly, the rate of entropy change (32) takes the form
d S d t = Δ S τ D ( S ˜ ρ | M ˜ ρ ) = Δ S r S M τ D
which, because | r S M | 1 , implies the bounds [equivalent to (38) and (53)],
Δ S τ D d S d t Δ S τ D .

6. Occupation Probabilities

An important class of observables for a quantum system are those associated with the projection operators. For example, for pure states evolving unitarily, the mean value P = Tr ( ρ ( t ) P ) where P = | ϕ 0 ϕ 0 | = ρ ( 0 ) represents the survival probability of the initial state, and is related to several notions of lifetimes [15,104,105].
We do not restrict our attention to pure states, and we discuss first results that hold for any projector P associated with a yes/no type of measurement. Let P = P = P 2 be an orthogonal projector onto the g-dimensional subspace P H of H . Clearly, g = Tr ( P ) , the variance Δ P Δ P = p ( 1 p ) where p = P = Tr ( ρ P ) denotes the mean value and represents the probability in state ρ of obtaining a ‘yes’ result upon measuring the associated observable. The characteristic time of the rate of change of this occupation probability is defined according to (19) by
1 τ P = | d p / d t | p ( 1 p ) = 2 d d t arccos ( p ) = 2 d d t arcsin ( p ) 1 τ U D ,
where the inequality follows from (48). Therefore,
1 2 τ U D d d t arccos ( p ) 1 2 τ U D ,
or, over any finite time interval of any time history p ( t ) ,
arccos ( p ( t 2 ) ) arccos ( p ( t 1 ) ) t 1 t 2 d t 2 τ U D ( t ) .
This result generalizes the results on lifetimes obtained in [103] where the focus is restricted to full quantum decay [ p ( ) 0 ] of an initially fully populated state [ p ( 0 ) 1 ] and τ U (here τ U D ) is assumed constant during the time interval. It is also directly related to some of the results in [15,104,105], where a number of additional inequalities and bounds on lifetimes are obtained for unitary dynamics, and may be straightforwardly generalized to the class of simultaneous unitary/dissipative dynamics described by our Equation (4).
Because p ( 1 p ) attains its maximum value when p = 1 / 2 , we also have the inequality
d p d t 1 2 τ U D .
which, analogously to what noted in [103], implies that no full decay nor full population can occur within a time 2 τ U D , so that this time may be interpreted as a limit to the degree of instability of a quantum state.
Next, we focus on the projectors onto the eigenspaces of the Hamiltonian operator H, assumed time-independent. Let us write its spectral expansion as H = n e n P e n where e n is the n-th eigenvalue and P e n the projector onto the corresponding eigenspace. Clearly, H P e n = e n P e n , P e n P e m = δ n m P e n , g n = Tr ( P e n ) is the degeneracy of eigenvalue e n , p n = P e n = Tr ( ρ P e n ) the occupation probability of energy level e n , Δ P e n Δ P e m = p n ( δ n m p m ) the covariance of pairs of occupations, and Δ P e n Δ P e n = p n ( 1 p n ) the variance or fluctuation of the n-th occupation. Because [ P e n , H ] = 0 , c P e n H = 0 and by (55) we have
d p n d t = Δ P e n r P e n M τ D ,
and the corresponding characteristic time is
1 τ P e n = | r P e n M | τ D 1 τ D .
Energy level occupation probabilities p n are used in Section 8 for numerical illustration/validation of inequalities (64) within the steepest-entropy-ascent dynamical model outlined in the next Section.

7. Example. Steepest-Entropy-Ascent Master Equation for Conservative Dissipative Dynamics

So far we have not assumed an explicit form of the operator Δ M ( ρ ) except for the condition that it maintains ρ unit trace ((3) or (8)) and nonnnegative definite. In this section, we illustrate the above results by further assuming a particular form of steepest-entropy-ascent, conservative dissipative dynamics. For our generalized nonequilibrium Massieu operator we assume the expression
Δ M ( ρ ) = Δ S Δ H ( ρ ) / θ ( ρ ) ,
where S is the entropy operator defined in Equation (33),
Δ H ( ρ ) = Δ H ν ( ρ ) · Δ N ,
H is the Hamiltonian operator, N = { N 1 , , N r } a (possibly empty) set of operators commuting with H that we call non-Hamiltonian generators of the motion (for example, the number-of-particles operators or a subset of them, or the momentum component operators for a free particle) and that must be such that operators ρ Δ H and ρ Δ N are linearly independent, and—most importantly— θ ( ρ ) and ν ( ρ ) = { ν 1 ( ρ ) , , ν r ( ρ ) } are a set of real functionals defined for each ρ by the solution of the following system of linear equations
Δ S Δ H θ + i = 1 r Δ N i Δ H ν i = Δ H Δ H ,
Δ S Δ N j θ + i = 1 r Δ N i Δ N j ν i = Δ H Δ N j ,
which warrant the conditions that Δ H Δ M = 0 and Δ N j Δ M = 0 , and hence that the mean values Tr ( ρ H ) and Tr ( ρ N ) are maintained time invariant by the dissipative term of the resulting SEA master equation [Equation (4) together with Equations (65)–(68)].
As a result, our assumption may be rewritten as follows
Δ M ( ρ ) = M ( ρ ) I Tr [ ρ M ( ρ ) ]
where I is the identity and the nonequilibrium Massieu operator M ( ρ ) is the following nonlinear function of ρ
M ( ρ ) = S ( ρ ) H θ ( ρ ) + ν ( ρ ) · N θ ( ρ ) ,
and we note that at a thermodynamic equilibrium (Gibbs) state,
ρ e = 1 Z exp H μ e · N T e ,
its mean value belongs to the family of entropic characteristic functions introduced by Massieu [106], i.e.,
M e = S e H e T e + μ e · N e T e ,
where S e , H e , N e , T e = θ ( ρ e ) and μ e = ν ( ρ e ) are the (grand canonical) equilibrium entropy, energy, amounts of constituents, temperature and chemical potentials, respectively.
Notice that operator M, its eigenvalues and its mean value Tr ( ρ M ) for a given state ρ , that we first termed “nonequilibrium Massieu operator” in References [62,94,107], differ substantially from the “nonequilibrium Massieu potentials” defined recently in References [108,109]. Their nonequilibrium Massieu construct is defined by the difference between the entropy and a linear combination of the conserved properties, with coefficients that are weighted averages of the fixed temperatures and other entropic potentials of the reservoirs interacting with the system. In our nonequilibrium Massieu construct, instead, the coefficients θ and ν of the linear combination are truly nonequilibrium functionals of the state ρ , that evolve in time with ρ , and that only when the system has relaxed to equilibrium can be identified with the inverse temperature 1 / T of the system and the entropic potentials μ / T of the other conserved properties.
The non-Hamiltonian generators of the motion represent the other conserved properties of the system, however, this condition may be relaxed in the framework of a resource theory of a quantum thermodynamic subsystem that, via the Hamiltonian part of the master equation, exchanges with other systems or a thermal bath some non-commuting quantities or “charges”, as recently envisioned in Reference [110].
Operators ρ Δ H and ρ Δ M are always orthogonal to each other, in the sense that Δ M Δ H = 0 for every ρ . It follows that, in general, Δ S Δ H = Δ H Δ H / θ ,
Δ S Δ M = Δ M Δ M = Δ S Δ S Δ H Δ H θ 2 ( ρ ) 0 ,
and hence the rate of entropy generation (32) is always strictly positive except for Δ M Δ M = 0 (which occurs iff ρ Δ M = 0 ), i.e., for ρ nd Δ S nd = ( ρ nd Δ H μ nd · ρ nd Δ N ) / T nd , for some real scalars T nd and μ nd , that is, for density operators (that we call non-dissipative [58,62,94,107]) of the following Gibbs (or partially Gibbs, if B I ) form
ρ nd = B exp [ ( H μ nd · N ) / k B T nd ] B Tr B exp [ ( H μ nd · N ) / k B T nd ] ,
where B is any projection operator on H ( B 2 = B ).
The nonlinear functional
θ ( ρ ) = σ H H σ S H = Δ H Δ S r S H
may be interpreted in this framework as a natural generalization to nonequilibrium of the temperature, at least insofar as for t + , while the state operator ρ ( t ) approaches a non-dissipative operator of form (74), θ ( ρ ( t ) ) approaches smoothly the temperature T nd of the non-dissipative thermodynamic equilibrium (stable, if B = I , or unstable, if B I ) or of the unstable limit cycle (if [ B , H ] 0 ), and ν ( ρ ( t ) ) / θ ( ρ ( t ) ) approach smoothly the corresponding entropic potentials μ nd / T nd .
Because here we assumed that H always commutes with M, c M H = 0 and ( M ˜ | i H ˜ ) = 0 , which means that ρ Δ M ( ρ ) is always orthogonal to i ρ Δ H . This reflects the fact that the direction of steepest-entropy-ascent is orthogonal to the (constant entropy) orbits that characterize purely Hamiltonian (unitary) motion (which maintains the entropy constant by keeping invariant each eigenvalue of ρ ).
Here, for simplicity, we have assumed that dissipation pulls the state in the direction of steepest-entropy-ascent with respect to the uniform Fisher–Rao metric (see [62]). However, we have discussed elsewhere (see [34,35]) that, in general, a most important and characterizing feature of the nonequilibrium states of a system is the metric with respect to which the system identifies the direction of steepest-entropy-ascent. In most cases, it is a non-uniform metric, such as for a material with a nonisotropic thermal conductivity or, in the quantum framework, for a spin system in a magnetic field that near equilibrium obeys the Bloch equations [111] with different relaxation times along the field and normal to the field.
Inequality (73), which follows from r S M 2 1 , implies that σ M M σ S S and 0 r S M = Δ M / Δ S 1 or, equivalently,
τ K = k B τ / Δ S τ D ,
where for convenience we define the characteristic time τ K , which is simply related to the entropy uncertainty, but cannot be attained by any rate of change, being shorter than τ D . In addition, we have the identities
r S M 2 = σ M M σ S S = τ K 2 τ D 2 = τ K τ S = 1 σ H H θ 2 σ S S = 1 r S H 2 ,
and, from r S H 2 1 , the bounds
| θ | Δ H Δ S or Δ S Δ H 1 θ Δ S Δ H ,
where the equality | θ | = Δ H / Δ S holds when and only when the state is non-dissipative [Equation (74)]. Additional bounds on our generalized nonequilibrium temperature θ obtain by combining (77) with the inequality 4 r S M 2 ( 1 r S M 2 ) 1 (which clearly holds because r S M 2 1 ), to obtain 4 r S M 2 r S H 2 1 and, therefore,
2 Δ M Δ H | θ | σ S S 1 or σ S S 2 Δ M Δ H 1 θ σ S S 2 Δ M Δ H .
At equilibrium, Δ M = 0 and (79) implies no actual bound on θ , but in nonequilibrium states bounds (79) may be tighter than (78), as illustrated by the numerical example in Section 8.
Notice that whereas in steepest-entropy-ascent dynamics τ K is always shorter than τ D and obeys the identity
τ S τ K = τ D 2 ,
in general it is not necessarily shorter than τ D and obeys the identity
Δ M Δ S τ D 2 τ S τ K = | r S M | .
In summary, we conclude that within steepest-entropy-ascent, conservative dissipative quantum dynamics, the general uncertainty relations (28), (37) and (38) that constitute the main results of this paper, yield the time–energy/time-Massieu uncertainty relation
τ F Δ H / 2 2 + τ F Δ M k B τ 2 1 or τ F 2 τ U 2 + τ F 2 τ D 2 1 ,
which implies the interesting time–energy and time–entropy uncertainty relation
τ F Δ H / 2 2 + τ F Δ S k B τ 2 1 or τ F 2 τ U 2 + τ F 2 τ K 2 1 ,
and the time–entropy uncertainty relation
τ K τ S = k B τ τ S Δ S = r S M 2 r S M 1 ,
which implies that the rate of entropy generation never exceeds σ S S / k B τ , i.e.,
d S d t = k B d d t Tr ( ρ ln ρ ) = σ M M k B τ Δ S Δ M k B τ σ S S k B τ .
If in addition the dynamics is purely dissipative, such as along a trajectory ρ ( t ) that commutes with H for every t, then (83) may be replaced by the time–entropy uncertainty relation
τ K τ F = k B τ τ F Δ S 1 .
As shown in References [58,62], the dissipative dynamics generated by Equation (4) with Δ M ( ρ ) as just defined and a time-independent Hamiltonian H: (i) maintains ρ ( t ) ρ 2 ( t ) at all times, both forward and backward in time for any initial density operator ρ ( 0 ) (see also [90,91]); (ii) maintains the cardinality of ρ ( t ) invariant; (iii) entails that the entropy functional is an S-function in the sense defined in [112] and therefore that maximal entropy density operators (Gibbs states) obtained from (74) with B = I are the only equilibrium states of the dynamics that are stable with respect to perturbations that do not alter the mean values of the energy and the other time invariants (if any): this theorem of the dynamics coincides with the Hatsopoulos-Keenan statement of the second law of thermodynamics [86]; (iv) entails Onsager reciprocity in the sense defined in [113]; (v) can be derived from a variational principle [90,91], equivalent to our steepest-entropy-ascent geometrical construction, by maximizing the entropy generation rate subject to the Tr ( ρ ) , Tr ( ρ H ) , and Tr ( ρ N ) conservation constraints and the additional constraint ( ρ E | ρ E ) = c ( ρ ) .
Operator ρ E is a ‘vector’ in L ( H ) and determines through its scalar product with ρ F and ρ S [Equations (7) and (32)] the rates of change of Tr ( ρ F ) and Tr ( ρ S ) , respectively. From (32) and the Schwarz inequality ( ρ S | ρ E ) 2 ( ρ S | ρ S ) ( ρ E | ρ E ) , we see that for a given ρ , among all vectors ρ E with given norm ( ρ E | ρ E ) = c ( ρ ) , the one maximizing ( ρ S | ρ E ) has the same direction as ρ S . In general, along such direction Tr ( ρ H ) and Tr ( ρ N ) are not conserved because ρ S is not always orthogonal to ρ H and ρ N . Instead, dynamics along the direction of steepest-entropy-ascent compatible with such conservation requirements, as first postulated and formulated in [57,58,62], obtains when ρ E has the direction of the component of ρ S orthogonal to ρ H and ρ N . This is precisely how Δ M ( ρ ) is defined through Equations (65)–(68). See also Reference [100].
We finally note that assuming in Equation (4) a Δ M ( ρ ) / τ that satisfied Equation (26) with strict equality, we obtain the most dissipative (maximal entropy generation rate) dynamics in which the entropic characteristic time τ S (Equation (36)) is always compatible with the time–energy uncertainty relation τ S Δ H / 2 and the rate of entropy generation is always given by 2 Δ M Δ H / .
The physical meaning of relations (28), (37), (38), (83), (84) are worth further investigations and experimental validation in specific contexts in which the dissipative behavior is correctly modeled by a dynamical law of form (4), possibly with Δ M ( ρ ) / τ of form (65). One such context may be the currently debated so-called “fluctuation theorems” [114,115,116,117] whereby fluctuations and, hence, uncertainties are measured on a microscopic system (optically trapped colloidal particle [118,119], electrical resistor [120]) driven at steady state (off thermodynamic equilibrium) by means of a work interaction, while a heat interaction (with a bath) removes the entropy being generated by irreversibility. Another such context may be that of pion-nucleus scattering, where available experimental data have recently allowed partial validation [121] of “entropic” uncertainty relations [122,123,124]. Yet another is within the model we propose in Reference [94] for the description of the irreversible time evolution of a perturbed, isolated, physical system during relaxation toward thermodynamic equilibrium by spontaneous internal rearrangement of the occupation probabilities. We pursue this example in the next section.

8. Numerical Results for Relaxation within a Single N-Level Qudit or a One-Particle Model of a Dilute Boltzmann Gas of N-Level Particles

To illustrate the time dependence of the uncertainty relations derived in this paper, we consider an isolated, closed system composed of noninteracting identical particles with single-particle eigenstates with energies e i for i = 1 , 2,..., N, where N is assumed finite for simplicity and the e i ’s are repeated in case of degeneracy, and we restrict our attention to the class of dilute-Boltzmann-gas states in which the particles are independently distributed among the N (possibly degenerate) one-particle energy eigenstates. This model is introduced in Reference [94], where we assume an equation of form (4) with Δ M ( ρ ) given by (65) with the further simplification that Δ H ( ρ ) = Δ H so that our generalized nonequilibrium Massieu operator is simply
M ( ρ ) = S H / θ ( ρ ) ,
and, therefore,
Δ M ( ρ ) = Δ S Δ H / θ ( ρ ) .
For simplicity and illustrative purposes, we focus on purely dissipative dynamics by considering a particular trajectory ρ ( t ) that commutes with H at all times t, assuming that H is time independent and has a nondegenerate spectrum. As a result, the energy-level occupation probabilities p n coincide with the eigenvalues of ρ , and the dynamical equation reduces to the simple form [94]
d p n d t = 1 τ p n ln p n + p n S k B + p n e n H k B θ ,
where
S = k B n p n ln p n ,
H = n p n e n ,
θ = σ H H / σ H S ,
σ H H = n p n e n 2 H 2 ,
σ H S = k B n p n e n ln p n H S .
The same model describes relaxation to the Gibbs state of an N-level qudit with time independent Hamiltonian H from arbitrary initial states ρ ( 0 ) that commute with H.
To obtain the plots in Figure 1 and Figure 2, that illustrate the main inequalities derived in this paper for a sample trajectory, we consider an initial state with cardinality equal to 4, with nonzero occupation probabilities only for the four energy levels e 1 = 0 , e 2 = u / 3 , e 3 = 2 u / 3 , and e 4 = u , and with mean energy H = 2 u / 5 (u is arbitrary, with units of energy). Moreover, as done in [94], we select an initial state ρ ( 0 ) at time t = 0 such that the resulting trajectory ρ ( t ) passes in the neighborhood of the partially canonical nondissipative state ρ nd ft that has nonzero occupation probabilities only for the three energy levels e 1 , e 2 , and e 4 , and mean energy H = 2 u / 5 ( p nd 1 ft = 0 . 3725 , p nd 2 ft = 0 . 3412 , p nd 3 ft = 0 , p nd 4 ft = 0 . 2863 , θ nd ft = 3 . 796 u / k B ). As shown in Figure 1, during the first part of the trajectory, this nondissipative state appears as an attractor, an approximate or ‘false target’ equilibrium state; when the trajectory gets close to this state, the evolution slows down, the entropy generation drops almost to zero and the value of θ gets very close ( 3 . 767 u / k B ) to that of θ nd ft ; however eventually the small, but nonzero initial occupation of level e 3 builds up and a new rapid rearrangement of the occupation probabilities takes place, and finally drives the system toward the maximal entropy state ρ nd pe with energy H = 2 u / 5 and all four active levels occupied, with canonical (Gibbs) distribution p nd 1 e = 0 . 3474 , p nd 2 e = 0 . 2722 , p nd 3 e = 0 . 2133 , p nd 4 e = 0 . 1671 , and characterized by the equilibrium temperature T e = 1 . 366 u / k B .
The trajectory is computed by integrating Equation (89) numerically, both forward and backward in time, starting from the chosen initial state ρ ( 0 ) , and assuming for Figure 1a and Figure 2a that the dissipation time τ is a constant, and for Figure 1b and Figure 2b that it is given by (26) with strict equality ( a τ = 1 , τ D = τ U ), i.e., assuming
τ = / 2 k B Δ M Δ H = / 2 k B σ S S σ H H 1 θ 2 ,
σ S S = k B 2 n p n ( ln p n ) 2 S 2 .
The system of ordinary differential Equation (89) is highly nonlinear, especially when τ is assumed according to (95), nevertheless it is sufficiently well behaved to allow simple integration by means of a standard Runge–Kutta numerical scheme. Of course, we check that at all times < t < each p n remains nonnegative, n p n remains equal to unity, n p n e n remains constant at the value 2 u / 5 fixed by the selected initial state, and the rate of change of S is always nonnegative.
In each Figure, the top subfigure shows for ease of comparison the plots of the four nonzero occupation probabilities as functions of dimensionless time: t / τ , in Figure 1a and Figure 2a; u t / , in Figure 1b and Figure 2b. The dots on the right represent the maximal entropy distribution, p n ( + ) = p n e ; the dots at the left represent the lowest-entropy or ‘primordial’ distribution, p n ( ) = p nd ls , which for the particular trajectory selected here, corresponds to a nondissipative state ρ nd ls that has only two occupied energy levels, e 1 and e 4 , with probabilities p nd 1 ls = 0 . 6 and p nd 4 ls = 0 . 4 (and temperature T nd ls = 2 . 466 u / k B ); in fact the four-level system has no lower entropy states ρ that commute with H, have energy 2 u / 5 , and have zero occupation probabilities [94]. The dots in the middle represent the nondissipative state ρ nd ft which appears as the false target state during the first part of the trajectory, plotted at the instant in time when the entropy of the time-varying trajectory is equal to the entropy of this distribution.
It is interesting to observe from Figure 1 (bottom subfigures) that during the early part of the trajectory, τ D almost exactly coincides with τ P e 2 while in the late part it almost exactly coincides with τ P e 3 , and the switch occurs when the trajectory slows down in the neighborhood of the ‘false target’ nondissipative state.
In Figure 1, the second subfigures show the time dependence of the dimensionless entropy S / k B ; the third subfigures show its rate of change (proportional to σ M M ) and compares it with σ S S and σ H H / θ 2 , to illustrate relation (73); the fourth show the time dependence of our generalized ‘nonequilibrium temperature’ θ (properly nondimensionalized) and compares it with Δ H / Δ S and 2 Δ M Δ H / σ S S to illustrate relations (78) and (79); the fifth subfigures show the time dependence of 1 / τ D (which here is proportional to the square root of the rate of entropy generation, third subfigures) and compares it with 1 / τ S and 1 / τ K to illustrate relations (53) and (76); the sixth subfigures show 1 / τ P e n for each of the four occupation probabilities and compares them with 1 / τ D to illustrate relation (64), which for this particular trajectory has the feature we just discussed.
In Figure 2, the second subfigures illustrate again relation (64) for each of the four observables p n = P e n ; the third subfigures illustrate the time–entropy uncertainty relation (86) for the same observables; the fourth illustrate inequality (62); the fifth illustrate relations (53) and (84).
By comparing subfigures (a) and (b) in both Figure 1 and Figure 2, it is noted that most qualitative features remain the same when τ is changed from constant to the state-dependent functional defined by Equation (95), except for the almost singular behavior near the false target partially canonical nondissipative state, where Δ M approaches zero and so does the dissipative time τ [Equation (95)]. The approach to final equilibrium in this case is not exponential in time as for τ = const. This puzzling behavior suggests that assumption (95) may hardly be physically sensible. However, as already noted after (24), it represents an interesting extreme behavior, i.e., the minimum dissipative time functional τ by which observables that commute with H, like the occupations P e n , never violate the usual time–energy uncertainty relations τ P e n Δ H / 2 , even though their time dependence is not determined here by unitary dynamics but by purely dissipative dynamics. These usual time–energy uncertainty relations, τ P e n τ U , are illustrated by the second row subfigures of Figure 2, because in this case τ U = τ D .

9. Conclusions

The Mandelstam–Tamm–Messiah time–energy uncertainty relation τ F Δ H / 2 provides a general lower bound to the characteristic times of change of all observables of a quantum system that can be expressed as linear functionals of the density operator ρ . This has been used to obtain estimates of rates of change and lifetimes of unstable states, without explicitly solving the time dependent evolution equation of the system. It may also be used as a general consistency check in measurements of time dependent phenomena. In this respect, the exact relation and inequalities (22) [that we derive for standard unitary dynamics based on the generalized Schrödinger inequality (14)] provide, for unitary evolution, a more general and sharper chain of consistency checks than the usual time–energy uncertainty relation.
The growing interest during the last three or four decades in quantum dynamical models of systems undergoing irreversible processes has been motivated by impressive technological advances in the manipulation of smaller and smaller systems, from the micrometer scale to the nanometer scale, and down to the single atom scale. The laws of thermodynamics, that fifty years ago were invariably understood as pertaining only to macroscopic phenomena, have gradually earned more attention and a central role in studies of mesoscopic phenomena first, and of microscopic and quantum phenomena more recently. In this paper we do not address the controversial issues currently under discussion about interpretational matters, nor do we attempt a reconstruction and review of the different views, detailed models and pioneering contributions that propelled during the past two decades this fascinating advance of thermodynamics towards the realm of few particle and single particle systems.
Motivated by this context and background, we derive various extensions of the usual time–energy uncertainty relations that may become useful in phenomenological studies of dissipative phenomena. We do so by focusing on a special but broad class of model evolution equations, that has been designed for the description of dissipative quantum phenomena and for satisfying a set of strict compatibility conditions with general thermodynamic principles. In this framework, we derive various forms of considerably precise time–energy and time–entropy uncertainty relations, and other interesting general inequalities, that should turn out to be useful at least as additional consistency checks in measurements of nonequilibrium states and time-dependent dissipative phenomena. To illustrate the qualitative features and the sharpness of the bounds provided by this set of inequalities, we show and discuss a numerical example obtained by integration (forward and backward in time) of the nonlinear evolution equation in the specific form introduced by this author for the description of steepest-entropy-ascent dynamics of an isolated system far from thermodynamic equilibrium.

Funding

This research received no external funding.

Acknowledgments

The author is grateful to Victor Dodonov for an interesting discussion during his visit in Brescia on 11 July 2007 on the viewgraphs I had just presented at the Conference on “Quantum Theory: Reconsideration of Foundations–4” in Växjö, Sweden, 11–16 June 2007, based on the early version of the present paper I had uploaded in ArXiv in 2005 [95], also presented at the 10th International Conference on “Squeezed States and Uncertainty Relations,” Bradford, UK, 2 April 2007 with the title “Time-Energy and Time-Entropy Uncertainty Relations in Steepest-Entropy-Ascent Dissipative Dynamics.” For the sake of historical record, early versions of this paper were submitted to Physical Review Letters (LL10220, December 2005) and to Physical Review A (LL10220A, April 2006) but were rejected after five mixed peer reviews.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Reasons for Not Assuming a Kossakowski–Lindblad form of the Master Equation

With various motivations, fundamental or phenomenological, dissipative quantum dynamical models, i.e., evolution equations for the density operator ρ that do not conserve the functional Tr ( ρ ln ρ ) , are almost invariably based on the KLGS master equations. For example, in theories of open systems in contact with a heat bath, or subsystems of a composite system which as a whole evolves unitarily, a variety of successful model evolution equations for the reduced density operator of the system have the KLGS form [25,26,27,28,29,30,31,32,33]
d ρ d t = i [ H , ρ ] + 1 2 j 2 V j ρ V j { V j V j , ρ } ,
where the V j ’s are operators on H (each term within the summation, often written in the alternative form [ V j , ρ V j ] + [ V j ρ , V j ] , is obviously traceless). Evolution equations of this form are linear in the density operator ρ and preserve its hermiticity, nonnegativity and trace.
For example, in a number of successful models of dissipative quantum dynamics of open subsystems, operators V j are in general interpreted as creation and annihilation, or transition operators. For example, by choosing V j = c r s | r s | , where c r s are complex scalars and | s eigenvectors of the Hamiltonian operator H, and defining the transition probabilities w r s = c r s c r s * , Equation (A1) becomes
d ρ d t = i [ H , ρ ] + r s w r s | s ρ s | 1 2 { | s s | , ρ } ,
or, equivalently, for the n m -th matrix element of ρ in the H representation,
d ρ n m d t = i ρ n m ( E n E m ) + δ n m r w n r ρ r r ρ n m 1 2 r ( w r n + w r m ) ,
which, for the occupation probabilities p n = ρ n n , is the Pauli master equation
d p n d t = r w n r p r p n r w r n .
Equation (A1) has also the intriguing feature of generating a completely positive dynamical map. However, Reference [125] argues quite clearly that the requirement of complete positivity of the reduced dynamics is too restrictive, as it is physically unnecessary to assure preservation of positivity of the density operator of the composite of any two noninteracting, uncorrelated systems.
Our objective here, instead, is to consider a class of model evolution equations applicable not only to open systems but also to closed isolated systems, capable of describing, simultaneously with the usual Hamiltonian unitary evolution, the natural tendency of any initial nonequilibrium state to relax towards canonical or partially-canonical thermodynamic equilibrium, i.e., capable of describing the irreversible tendency to evolve towards the highest entropy state compatible with the instantaneous mean values of the energy, the other constants of the motion, and possibly other constraints. To avoid the severe restrictions imposed by the linearity of the evolution equation, we open our attention to nonlinearity in the density operator ρ [84]. Therefore, it may at first appear natural to maintain the Kossakowski–Lindblad form (A1) and simply assume operators V j that are functions of ρ . This is true only in part for the evolution Equation (4) that we assume. Indeed, our hermitian operator Δ M ( ρ ) / k B τ can always be written as j V j ( ρ ) V j ( ρ ) and therefore our anticommutator term may be viewed as a generalization of the corresponding term in (A1).
However, in our Equations (1) and (4) we suppress the term corresponding to j V j ρ V j in (A1). The reason for this suppression is the following. Due to the terms V j ρ V j , whenever the state operator ρ is singular, i.e., it has one or more zero eigenvalues, Equation (A1) implies that these zero eigenvalues may change at a finite rate. This can be seen clearly from (A4) by which d p n / d t is finite whenever there is a nonzero transition probability w n r from some other populated level ( p r 0 ), regardless of whether p n is zero or not. When this occurs, for one instant in time the rate of entropy change is infinite, as seen clearly from the expression of the rate of entropy change implied by (A1),
d S d t = k B j Tr ( V j V j ρ ln ρ V j ρ V j ln ρ ) = k B j r n ( V j ) n r * ( V j ) n r ( ρ r ρ n ) ln ρ r ,
where ρ r denotes the r-th eigenvalue of ρ and ( V j ) n r the matrix elements of V j in the ρ representation.
We may argue that an infinite rate of entropy change can be tolerated, because it would last only for one instant in time. But the fact that zero eigenvalues of ρ in general could not survive, i.e., would not remain zero (or close to zero) for longer than one instant in time, is an unphysical feature, at least because it is in contrast with a wealth of successful models of physical systems in which great simplification is achieved by limiting our attention to a restricted subset of relevant eigenstates (forming a subspace of H that we call the effective Hilbert space of the system [81]). Such common practice N-level models yield extremely good results, that being reproducible, ought to be relatively robust with respect to including in the model other less relevant eigenstates. In fact, such added eigenstates, when initially unpopulated, are irrelevant if they remain unpopulated (or very little populated) for long times, so that neglecting their existence introduces very little error. The terms V j ρ V j , instead, would rapidly populate such irrelevant unpopulated eigenstates and void the validity of our so successful simple N-level models, unless we deliberately overlook this instability problem by highly ad-hoc assumption, e.g., by forcing the V j ’s to be such that ( V j ) n r = 0 whenever either ρ n = 0 or ρ r = 0 , in which case, however, we can no longer claim true linearity with respect to ρ .
To avoid the unphysical implications of this seldom recognized [81,94] problem of linear evolution equations of form (A1), we consider in this paper only equations of form (4). We do not exclude that it may be interesting to investigate also the behavior of equations that include nonlinear terms of the form V j ( ρ ) ρ V j ( ρ ) . However, at least when the system is strictly isolated, the operator-functions V j ( ρ ) should be such that ( V j ( ρ ) ) n r = 0 whenever either ρ n = 0 or ρ r = 0 .
Another important general physical reason why we exclude terms that generate nonzero rates of change of zero eigenvalues of ρ , is that if such terms are construed so as to conserve positivity in forward time, in general they cannot maintain positivity in backward time. The view implicitly assumed when Equation (A1) is adopted, is that the model is “mathematically irreversible” (a distinguishing feature if not a starting point of the theory of completely positive linear dynamical semigroups on which it is based), in the sense that neither uniqueness of solutions in forward time nor existence in backward time are required (and granted). Such mathematical irreversibility of the initial value problem, is often accepted, presented and justified as a natural counterpart of physical irreversibility. However, it is more related to the principle of causality than to physical irreversibility. The strongest form of the non-relativistic principle of causality—a keystone of traditional physical thought—requires that future states of a system should unfold deterministically from initial states along smooth unique trajectories in state domain defined for all times (future as well as past). Accepting mathematical irreversibility of the model dynamics implies giving up such causality requirement. The point is that such requirement is not strictly necessary to describe physical irreversibility, at least not if we are willing to give up linearity instead. The proof of this statement is our Equation (4) which, together with the additional assumptions made in Section 7 to describe relaxation within an isolated system, is mathematically reversible, in the sense that it features existence and uniqueness of well-defined solutions both in forward and backward time, and yet it does describe physically irreversible time evolutions, in the sense that the physical property described by the entropy functional k B Tr ( ρ ln ρ ) is a strictly increasing function of time for all states except the very restricted subset defined by Equation (74), where it is time invariant.

Appendix B. How Did Locally Steepest Entropy Ascent Come About?

The KLGS master equation emerges from a bottom-up phenomenological approach, whereby one considers a weakly interacting system+bath isolated composite evolving under the phenomenological assumption that the large number of degrees of freedom of the bath dilutes and destroys the correlations that build up under the standard unitary evolution due to the interaction term in the Hamiltonian, so that for the purpose of computing the evolution of the reduced density operator ρ S of the system, the overall state can be assumed to evolve through uncorrelated states, i.e., the system+bath density operator can at all times be written as ρ = ρ S ρ B . So, we may say that the derivation is “bottom-up” because it starts from the fundamental unitary evolution of the composite system, but it is also “phenomenological” because the assumption of loss of correlations is an approximation that depends on the bath details and requires neglecting some terms during the partial tracing over the bath subspace.
By contrast, the locally steepest-entropy-ascent master equation was originally constructed (not derived) from a “top-down” approach, meant to see what modifications of standard quantum mechanics would be required if one wants to embed the second law of thermodynamics directly into the fundamental law of description. This heretic, but certainly intriguing and thought-provoking theoretical exercise, belongs to the early history of quantum thermodynamics and should not be forgotten. In the 70’s, the need for a quantum thermodynamics had been addressed boldly and explicitly only by Hatsopoulos and Gyftopoulos [63,64,65,66] and—with a very different approach that we do not review here—by Prigogine and the Brussels school of thermodynamics [126,127,128]. As an additional note pertaining to the history of thermodynamics, the course 2.47 J/22.58 J, listed in the MIT Bulletin of the academic year 1970-71 and taught jointly by George N. Hatsopoulos and Elias P. Gyftopoulos in the Spring of 1971, is the first official course entitled “Quantum Thermodynamics” that we are aware of.
An emphatic way to explain the (philosophical?) intuition of these pioneers of quantum thermodynamics, is the set of “what if” questions posed by the present author in a famous conference on the frontiers of nonequilibrium statistical physics held in Santa Fe in 1984 [129]: “what if entropy, rather than a statistical, information theoretic, macroscopic or phenomenological concept, were an intrinsic property of matter in the same sense as energy is universally understood to be an intrinsic property of matter? What if irreversibility were an intrinsic feature of the fundamental dynamical laws obeyed by all physical objects, macroscopic and microscopic, complex and simple, large and small? What if the second law of thermodynamics, in the hierarchy of physical laws, were at the same level as the fundamental laws of mechanics, such as the great conservation principles? Is it inevitable that the gap between mechanics and thermodynamics be bridged by resorting to the usual statistical, phenomenological, or information-theoretic reasoning, and by hinging on the hardly definable distinction between microscopic and macroscopic reality? Is it inevitable that irreversibility be explained by designing ad hoc mechanisms of coupling with some heat bath, reservoir or environment, and ad hoc mechanisms of loss of correlation? What if, instead, mechanics and thermodynamics were both special cases of a more general unified fundamental physical theory valid for all systems, including a single strictly isolated particle, such as a single isolated harmonic oscillator or a single isolated two-level spin system?”
In References [63,64,65,66] Hatsopoulos and Gyftopoulos showed that the only price we have to pay [130,131] to gain a possible positive answer to these questions is the reinterpretation of the physical meaning of the density operator, abandoning the standard (statistical mechanics and information theoretic) interpretation whereby it represents the epistemic ignorance of which particular pure state the system is ‘really’ in. Instead, the density operator acquires an ‘ontic’ status and represents the individual state of the isolated and uncorrelated atom (or particle or indivisible entity of the system’s model). Equivalently, in terms of ensembles, the density operator represents the measurement statistics from a homogeneous ensemble—homogeneous in the sense defined by von Neumann and discussed in References [57,82,87,88], i.e., such that no subsensemble can be identified which gives rise to different measurement statistics.
As a result of such ansatz, the severe restrictions imposed by linearity on the evolution equation become unnecessary, and we must open up our attention to evolution equations nonlinear in the density operator ρ . This is what prompted the search for a generalization of the time-dependent Schrödinger equation for pure density operators, to a broader fundamental kinematics in which not only every pure density operator, but also every non-pure density operator represents a real ontological object, the ‘true’ state of the system, which can be mixed even if the system is isolated and uncorrelated (no entanglement) from the rest of the universe. Among the desiderata [81] that drove the search for an extension, the most important was that the second law should emerge as a theorem of the equation of motion, which we considered the strongest way to enforce strong compatibility of a dynamical model (or law of motion or resource theory...) with thermodynamics.
For such purpose, the Hatsopoulos–Keenan statement of the second law is particularly suited, because it is directly linked with stability features of the equilibrium states of the dynamics. This somewhat still overlooked statement of the second law asserts ([86], p. 62) that for any well-defined (i.e., separable and uncorrelated) system, among the set of states that share the same values of the parameters of the Hamiltonian and the (mean) values of the energy and the amounts of constituents, there exists one and only one (conditionally [112,132]) stable equilibrium state, which turns out to be the one with the maximal entropy, often called the Gibbs state in the recent QT literature.
The Hatsopoulos–Keenan statement of the second law not only can be proved to entail the better known statements (Kelvin–Planck [86], p. 64; Clausius [86], p. 134; Caratheodory [86], p. 121), but—quite importantly for the current developments of quantum thermodynamics—it supports a rigorous operational definition of entropy as a general property of any uncorrelated (and unentangled) state of any well-separated system, valid not only for the stable equilibrium states of macroscopic systems but also for their nonequilibrium states (see Reference [133] and references therein) and providing a possible basis for its extension to systems with only few particles and quantum systems. Its extendability to correlated states of interacting or non-interacting systems is instead still the subject of intense debate, because the correlation entropy (often called mutual information), like the mean energy of interaction between the subsystems, is a well defined feature of the overall state of a composite system, but there is no unique way nor fundamental reason to allocate it among the subsystems and assign it to their local (reduced, marginal) states.
Using a similar operational definition, Hatsopoulos and Gyftopoulos in References [63,64,65,66] showed that the von Neumann entropy functional k B Tr ( ρ ln ρ ) fulfills the definition of entropy. It was for such pioneering QT framework that the present author designed the nonlinear master equation [58,59,60] and shortly thereafter [61] proved it admits a steepest-entropy-ascent (SEA) variational formulation that embodies at the local microscopic level the principle of maximal entropy production [43,44,62,90,94,100,134].
Starting in the mid eighties, during times when quantum thermodynamics was considered interesting only by a handful of pioneers, the present author has addressed both the quantum-foundations and mathematical–physics communities [61,77,78,80,100,107,129,130,131,135,136,137,138] and the engineering-thermodynamics and non-equilibrium-thermodynamics communities [79,139,140,141,142,143,144] to raise awareness about the requirement that fundamental and phenomenological models of dissipative and irreversible processes must incorporate, i.e., must not violate, general thermodynamic principles.
More recently, the SEA principle has been shown to encompass all the major levels of description of nonequilibrium dynamics and irreversible processes [94], including the general modeling structure known as metriplectic dynamics [145,146,147] or GENERIC (see [35] and references therein), which even more recently has been shown to bear deep connections also with the mathematical theories of gradient flows [148,149,150] and large fluctuations [151,152].
In other words, the locally steepest-entropy-ascent model of far-non-equilibrium dissipative evolution in QT can be considered the most general precursor of all more recent and successful theories of nonequilibrium dynamical systems.
In several instances and different fields of application the LSEA approach has shown the ability to provide new nontrivial modeling capabilities not only in the realm of QT and related phenomenological resource theories (see, e.g., References [68,70,76]) but also in materials science [71,72,73,75] and transport theory [74].
Also the unified theory presented in References [63,64,65,66], if one puts aside the epistemic interpretation and considers it as an effective resource theory, represents in our view the pioneering precursor of many quantum thermodynamics results that have been re-derived in recent years (free energy versus available energy, energy versus entropy diagram, work element, adiabatic availability, etc.).

References

  1. Aharonov, Y.; Bohm, D. Time in the quantum theory and the uncertainty relation for time and energy. Phys. Rev. 1961, 122, 1649–1658. [Google Scholar] [CrossRef]
  2. Allcock, G.R. The time of arrival in quantum mechanics I. Formal considerations. Ann. Phys. 1969, 53, 253–285. [Google Scholar] [CrossRef]
  3. Allcock, G.R. The time of arrival in quantum mechanics II. The individual measurement. Ann. Phys. 1969, 53, 286–310. [Google Scholar] [CrossRef]
  4. Allcock, G.R. The time of arrival in quantum mechanics III. The measurement ensemble. Ann. Phys. 1969, 53, 311–348. [Google Scholar] [CrossRef]
  5. Bauer, M.; Mello, P.A. The time-energy uncertainty relation. Ann. Phys. 1978, 111, 38–60. [Google Scholar] [CrossRef]
  6. Dodonov, V.V.; Kurmyshev, E.V.; Man’ko, V.I. Generalized uncertainty relation and correlated coherent states. Phys. Lett. A 1980, 79, 150–152. [Google Scholar] [CrossRef]
  7. Busch, P. On the energy-time uncertainty relation. Part I: Dynamical time and time indeterminacy. Found. Phys. 1990, 20, 1–32. [Google Scholar] [CrossRef]
  8. Busch, P. On the energy-time uncertainty relation. Part II: Pragmatic time versus energy indeterminacy. Found. Phys. 1990, 20, 33–43. [Google Scholar] [CrossRef]
  9. Landauer, R. Information is physical. Phys. Today 1991, 44, 23–29. [Google Scholar] [CrossRef]
  10. Hilgevoord, J. The uncertainty principle for energy and time I. Am. J. Phys. 1996, 64, 1451–1456. [Google Scholar] [CrossRef]
  11. Hilgevoord, J. The uncertainty principle for energy and time. II. Am. J. Phys. 1998, 66, 396–402. [Google Scholar] [CrossRef]
  12. Hu, B.L.; Zhang, Y. Uncertainty Principle for Quantum Open Systems. Int. J. Mod. Phys. A 1995, 10, 4537–4561. [Google Scholar] [CrossRef]
  13. Trifonov, D.A. Generalizations of Heisenberg uncertainty relation. Eur. Phys. J. B 2002, 29, 349–353. [Google Scholar] [CrossRef]
  14. Dodonov, V.V. Purity- and entropy-bounded uncertainty relations for mixed quantum states. J. Opt. B 2002, 4, S98–S108. [Google Scholar] [CrossRef]
  15. Pfeifer, P.; Fröhlich, J. Generalized time-energy uncertainty relations and bounds on lifetimes of resonances. Rev. Mod. Phys. 1995, 67, 759–779. [Google Scholar] [CrossRef] [Green Version]
  16. Sackett, C.A.; Kielpinski, D.; King, B.E.; Langer, C.; Meyer, V.V.; Myatt, C.J.; Rowe, M.; Turchette, Q.A.; Itano, W.M.; Wineland, D.J.; et al. Experimental entanglement of four particles. Nature 2000, 404, 256–259. [Google Scholar] [CrossRef]
  17. Kielpinski, D.; Monroe, C.; Wineland, D.J. Architecture for a large-scale ion-trap quantum computer. Nature 2002, 417, 709–711. [Google Scholar] [CrossRef] [Green Version]
  18. Ekert, A.K.; Alves, C.M.; Oi, D.K.L.; Horodecki, M.; Horodecki, P.; Kwek, L.C. Direct estimations of linear and nonlinear functionals of a quantum state. Phys. Rev. Lett. 2002, 88, 217901. [Google Scholar] [CrossRef]
  19. Bovino, F.A.; Castagnoli, G.; Ekert, A.; Horodecki, P.; Alves, C.M.; Sergienko, A.V. Direct measurement of nonlinear properties of bipartite quantum states. Phys. Rev. Lett. 2005, 95, 240407. [Google Scholar] [CrossRef]
  20. Shull, C.G.; Atwood, D.K.; Arthur, J.; Horne, M.A. Search for a Nonlinear variant of the Schrödinger equation by neutron interferometry. Phys. Lett. Lett. 1980, 44, 765. [Google Scholar] [CrossRef]
  21. Martinez, J.C.; Polatdemir, E. Measurement of tunneling time via electron interferometry. Appl. Phys. Lett. 2004, 84, 1320–1322. [Google Scholar] [CrossRef]
  22. Demkowicz-Dobrzański, R.; Czajkowski, J.; Sekatski, P. Adaptive quantum metrology under general markovian noise. Phys. Rev. X 2017, 7, 041009. [Google Scholar] [CrossRef]
  23. Beau, M.; Del Campo, A. Nonlinear quantum metrology of many-body open systems. Phys. Rev. Lett. 2017, 119, 010403. [Google Scholar] [CrossRef] [PubMed]
  24. Dodonov, V.V.; Dodonov, A.V. Energy–time and frequency–time uncertainty relations: Exact inequalities. Phys. Scr. 2015, 90, 074049. [Google Scholar] [CrossRef]
  25. Kraus, K. General state changes in quantum theory. Ann. Phys. 1971, 64, 311–335. [Google Scholar] [CrossRef]
  26. Kossakowski, A. Necessary and sufficient conditions for a generator of a quantum dynamical semigroup. Bull. Acad. Pol. Sci. Ser. Math. Astr. Phys. 1972, 20, 1021–1025. [Google Scholar]
  27. Kossakowski, A. On quantum statistical mechanics of non-Hamiltonian systems. Rep. Math. Phys. 1972, 3, 247–274. [Google Scholar] [CrossRef]
  28. Ingarden, R.S.; Kossakowski, A. On the connection of nonequilibrium information thermodynamics with non-Hamiltonian quantum mechanics of open systems. Ann. Phys. 1975, 89, 451–485. [Google Scholar] [CrossRef]
  29. Lindblad, G. On the generators of quantum dynamical semigroups. Commun. Math. Phys. 1976, 119, 119–130. [Google Scholar] [CrossRef]
  30. Gorini, V.; Kossakowski, A.; Sudarshan, E.C.G. Completely positive dynamical semigroups of N-level systems. J. Math. Phys. 1976, 17, 821–825. [Google Scholar] [CrossRef]
  31. Spohn, H. Approach to equilibrium for completely positive dynamical semigroups of N-level systems. Rep. Math. Phys. 1976, 10, 189–194. [Google Scholar] [CrossRef]
  32. Spohn, H. Entropy production for quantum dynamical semigroups. J. Math. Phys. 1978, 19, 1227–1230. [Google Scholar] [CrossRef]
  33. Asorey, M.; Kossakowski, A.; Marmo, G.; Sudarshan, E.C.G. Dynamical maps and density matrices. J. Phys. Conf. Ser. 2009, 196, 012023. [Google Scholar] [CrossRef]
  34. Beretta, G.P. Steepest Entropy Ascent Model for Far-Non-Equilibrium Thermodynamics. Unified Implementation of the Maximum Entropy Production Principle. Phys. Rev. E 2014, 90, 042113. [Google Scholar] [CrossRef] [PubMed]
  35. Montefusco, A.; Consonni, F.; Beretta, G.P. Essential equivalence of the general equation for the nonequilibrium reversible-irreversible coupling (GENERIC) and steepest-entropy-ascent models of dissipation for nonequilibrium thermodynamics. Phys. Rev. E 2015, 91, 042138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Mittnenzweig, M.; Mielke, A. An entropic gradient structure for Lindblad equations and couplings of quantum systems to macroscopic models. J. Stat. Phys. 2017, 167, 205–233. [Google Scholar] [CrossRef]
  37. Kantner, M.; Mittnenzweig, M.; Koprucki, T. Hybrid quantum-classical modeling of quantum dot devices. Phys. Rev. B 2017, 96, 205301. [Google Scholar] [CrossRef] [Green Version]
  38. Messiah, A. Quantum Mechanics; Wiley: New York, NY, USA, 1976; p. 320. [Google Scholar]
  39. Mandelstam, L.; Tamm, I.G. The uncertainty relation between energy and time in non-relativistic quantum mechanics. In Selected Papers; Bolotovskii, B.M., Frenkel, V.Y., Peierls, R., Eds.; Springer: Berlin/Heidelberg, Germany, 1991; pp. 115–123. [Google Scholar] [CrossRef]
  40. Ziegler, H. An attempt to generalize Onsager’s principle, and its significance for rheological problems. J. Appl. Math. Phys. 1958, 9, 748–763. [Google Scholar] [CrossRef]
  41. Swenson, R. Emergent attractors and the law of maximum entropy production: Foundations to a theory of general evolution. Syst. Res. 1989, 6, 187–197. [Google Scholar] [CrossRef]
  42. Struchtrup, H.; Weiss, W. Maximum of the Local Entropy Production Becomes Minimal in Stationary Processes. Phys. Rev. Lett. 1998, 80, 5048–5051. [Google Scholar] [CrossRef] [Green Version]
  43. Dewar, R.C. Maximum entropy production and the fluctuation theorem. J. Phys. A 2005, 38, L371–L381. [Google Scholar] [CrossRef] [Green Version]
  44. Martyushev, L.M.; Seleznev, V.D. Maximum entropy production principle in physics, chemistry and biology. Phys. Rep. 2006, 426, 1–45. [Google Scholar] [CrossRef]
  45. Martyushev, L.M. The maximum entropy production principle: Two basic questions. Philos. Trans. R. Soc. B 2010, 365, 1333–1334. [Google Scholar] [CrossRef] [PubMed]
  46. Martyushev, L.M. Entropy and entropy production: Old misconceptions and new breakthroughs. Entropy 2013, 15, 1152–1170. [Google Scholar] [CrossRef]
  47. Martyushev, L.M.; Seleznev, V.D. The restrictions of the maximum entropy production principle. Physica A 2014, 410, 17–21. [Google Scholar] [CrossRef] [Green Version]
  48. Gislason, E.A.; Sabelli, N.H.; Wood, J.W. New form of the time-energy uncertainty relation. Phys. Rev. A 1985, 31, 2078–2081. [Google Scholar] [CrossRef]
  49. Bialynicki-Birula, I.; Mycielski, J. Uncertainty relations for information entropy in wave mechanics. Commun. Math. Phys. 1975, 44, 123–132. [Google Scholar] [CrossRef]
  50. Kobe, D.H.; Aguilera-Navarro, V.C. Derivation of the time-energy uncertainty relation. Phys. Rev. A 1994, 50, 933–938. [Google Scholar] [CrossRef]
  51. Majernik, V.; Richterek, L. Entropic uncertainty relations. Eur. J. Phys. 1997, 18, 79–89. [Google Scholar] [CrossRef]
  52. Aharonov, Y.; Reznik, B. “Weighing” a closed system and the time-energy uncertainty principle. Phys. Rev. Lett. 2000, 84, 1368–1370. [Google Scholar] [CrossRef]
  53. Aharanov, Y.; Massar, S.; Popescu, S. Measuring energy, estimating Hamiltonians, and the time-energy uncertainty relation. Phys. Rev. A 2002, 66, 052107. [Google Scholar] [CrossRef] [Green Version]
  54. Busch, P. The time-energy uncertainty relation. In Time in Quantum Mechanics; Muga, J., Mayato, R.S., Egusquiza, I., Eds.; Lecture Notes in Physics; Springer: Berlin/Heidelberg, Germany, 2008; Volume 734. [Google Scholar] [CrossRef]
  55. Brunetti, R.; Fredenhagen, K. Remarks on time-energy uncertainty relations. Rev. Math. Phys. 2002, 14, 897–903. [Google Scholar] [CrossRef]
  56. Gillies, G.T.; Allison, S.W. Experimental test of a time-temperature formulation of the uncertainty principle via nanoparticle fluorescence. Found. Phys. 2005, 18, 65–74. [Google Scholar] [CrossRef]
  57. Beretta, G.P. On the General Equation of Motion of Quantum Thermodynamics and the Distinction between Quantal and Nonquantal Uncertainties. Ph.D. Thesis, MIT, Cambridge, MA, USA, 1981. [Google Scholar]
  58. Beretta, G.P.; Gyftopoulos, E.P.; Park, J.L.; Hatsopoulos, G.N. Quantum thermodynamics. A new equation of motion for a single constituent of matter. Nuovo Cimento B 1984, 82, 169–191. [Google Scholar] [CrossRef]
  59. Beretta, G.P.; Gyftopoulos, E.P.; Park, J.L. Quantum thermodynamics. A new equation of motion for a general quantum system. Nuovo Cimento B 1985, 87, 77–97. [Google Scholar] [CrossRef] [Green Version]
  60. Maddox, J. Uniting mechanics and statistics. Nature 1985, 316, 11. [Google Scholar] [CrossRef]
  61. Beretta, G.P. Steepest entropy ascent in quantum thermodynamics. In The Physics of Phase Space, Nonlinear Dynamics and Chaos, Geometric Quantization, and Wigner Function, Proceedings of the First International Conference on the Physics of Phase Space, University of Maryland, College Park, MD, USA, 20–23 May 1986; Kim, Y.S., Zachary, W.W., Eds.; Springer: Berlin/Heidelberg, Germany, 1987; pp. 441–443. [Google Scholar] [CrossRef]
  62. Beretta, G.P. Nonlinear quantum evolution equations to model irreversible adiabatic relaxation with maximal entropy production and other nonunitary processes. Rep. Math. Phys. 2009, 64, 139–168. [Google Scholar] [CrossRef] [Green Version]
  63. Hatsopoulos, G.N.; Gyftopoulos, E.P. A unified quantum theory of mechanics and thermodynamics. Part I. Postulates. Found. Phys. 1976, 6, 15–31. [Google Scholar] [CrossRef]
  64. Hatsopoulos, G.N.; Gyftopoulos, E.P. A unified quantum theory of mechanics and thermodynamics. Part IIa. Available energy. Found. Phys. 1976, 6, 127–141. [Google Scholar] [CrossRef]
  65. Hatsopoulos, G.N.; Gyftopoulos, E.P. A unified quantum theory of mechanics and thermodynamics. Part IIb. Stable equilibrium states. Found. Phys. 1976, 6, 439–455. [Google Scholar] [CrossRef]
  66. Hatsopoulos, G.N.; Gyftopoulos, E.P. A unified quantum theory of mechanics and thermodynamics. Part III. Irreducible quantal dispersions. Found. Phys. 1976, 6, 561–570. [Google Scholar] [CrossRef]
  67. Brandão, F.G.L.S.; Horodecki, M.; Oppenheim, J.; Renes, J.M.; Spekkens, R.W. Resource theory of quantum states out of thermal equilibrium. Phys. Rev. Lett. 2013, 111, 250404. [Google Scholar] [CrossRef] [PubMed]
  68. Cano-Andrade, S.; Beretta, G.P.; von Spakovsky, M.R. Steepest-entropy-ascent quantum thermodynamic modeling of decoherence in two different microscopic composite systems. Phys. Rev. A 2015, 91, 013848. [Google Scholar] [CrossRef]
  69. Smith, C. Comparing the models of steepest entropy ascent quantum thermodynamics, master equation and the difference equation for a simple quantum system interacting with reservoirs. Entropy 2016, 18, 176. [Google Scholar] [CrossRef]
  70. Tabakin, F. Model dynamics for quantum computing. Ann. Phys. 2017, 383, 33–78. [Google Scholar] [CrossRef] [Green Version]
  71. Yamada, R.; von Spakovsky, M.R.; Reynolds, W.T., Jr. A method for predicting non-equilibrium thermal expansion using steepest-entropy-ascent quantum thermodynamics. J. Phys. 2018, 30, 325901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Yamada, R.; von Spakovsky, M.R.; Reynolds, W.T., Jr. Methodology of an application of the steepest-entropy- ascent quantum thermodynamic framework to physical phenomena in materials science. Comput. Mater. Sci. 2019, 166, 215–264. [Google Scholar] [CrossRef]
  73. Yamada, R.; von Spakovsky, M.R.; Reynolds, W.T., Jr. Predicting continuous and discontinuous phase decompositions using steepest-entropy-ascent quantum thermodynamics. Phys. Rev. E 2019, 99, 052121. [Google Scholar] [CrossRef] [Green Version]
  74. Li, G.; von Spakovsky, M.R.; Hin, C. Steepest entropy ascent quantum thermodynamic model of electron and phonon transport. Phys. Rev. B 2018, 97, 024308. [Google Scholar] [CrossRef] [Green Version]
  75. Kusaba, A.; Li, G.; von Spakovsky, M.R.; Kangawa, Y.; Kakimoto, K. Modeling the non-equilibrium process of the chemical adsorption of ammonia on GaN(0001) reconstructed surfaces based on steepest-entropy-ascent quantum thermodynamics. Materials 2017, 10, 948. [Google Scholar] [CrossRef]
  76. Militello, B. Steepest entropy ascent for two-state systems with slowly varying Hamiltonians. Phys. Rev. E 2018, 97, 052113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Beretta, G.P. A general nonlinear evolution equation for irreversible conservative approach to stable equilibrium. In Frontiers of Nonequilibrium Statistical Physics, Proceedings of the NATO Advanced Study Institute on the Frontiers of Nonequilibrium Statistical Physics, Santa Fe, New Mexico, 3–16 June 1984; Moore, G.T., Scully, M.O., Eds.; NATO ASI Series B; Springer: Boston, MA, USA, 1986; Volume 135, pp. 193–204. [Google Scholar] [CrossRef]
  78. Beretta, G.P. The role of stability in the unification of mechanics and thermodynamics. In Proceedings of the ANAE-ONERA-CERT International Symposium on Stability, Paris, France, 3–5 November 1987; pp. 87–104. [Google Scholar]
  79. Beretta, G.P. The Schroedinger-Park paradox about the concept of “state” in quantum statistical mechanics and quantum information theory is still open. One more reason to go beyond? In Beyond the Quantum, Proceedings of the Workshop “Beyond the Quantum”, Lorentz Center of Leiden University, Leiden, The Netherlands, 28 May–3 June 2006; Nieuwenhuizen, T.M., Spicka, V., Mehmani, B., Jafar-Aghdami, M., Khrennikov, A.Y., Eds.; World Scientific: Singapore, 2007; pp. 355–365. ISBN 981-277-117-4. [Google Scholar] [CrossRef]
  80. Beretta, G.P. Quantum thermodynamics: New light on the physical meaning of entropy and the origin of irreversibility. In Proceedings of the Taormina Conference on Thermodynamics, Taormina, Italy, 18–22 February 1991; Giaquinta, P.V., Vicentini Missoni, M., Wanderlingh, F., Eds.; Atti della Accademia Peloritana dei Pericolanti-Classe di Scienze Fisiche Matematiche e Naturali: Messina, Italy, 1992; Volume 70, pp. 61–99. [Google Scholar]
  81. Beretta, G.P. Nonlinear extensions of Schroedinger-von Neumann quantum dynamics: A set of necessary conditions for compatibility with thermodynamics. Mod. Phys. Lett. A 2005, 20, 977–984. [Google Scholar] [CrossRef]
  82. Beretta, G.P. The Hatsopoulos-Gyftopoulos resolution of the Schroedinger-Park paradox about the concept of “state” in quantum statistical mechanics. Mod. Phys. Lett. A 2006, 21, 2799–2811. [Google Scholar] [CrossRef]
  83. Park, J.L.; Band, W. Generalized two-level quantum dynamics. III. Irreversible conservative motion. Found. Phys. 1978, 8, 239–254. [Google Scholar] [CrossRef]
  84. Simmons, R.F., Jr.; Park, J.L. The essential nonlinearity of N-level quantum thermodynamics. Found. Phys. 1981, 11, 297–305. [Google Scholar] [CrossRef]
  85. Hatsopoulos, G.N.; Keenan, J.H. Principles of General Thermodynamics; Wiley: New York, NY, USA, 1965. [Google Scholar]
  86. Gyftopoulos, E.P.; Beretta, G.P. Thermodynamics: Foundations and Applications; Dover Publications: Mineola, NY, USA, 2005. [Google Scholar]
  87. Von Neumann, J. Mathematical Foundations of Quantum Mechanics; Beyer, R.T., Ed.; Princeton University Press: Princeton, NJ, USA, 1955; pp. 295–346. [Google Scholar]
  88. Park, J.L. Nature of quantum states. Am. J. Phys. 1968, 36, 211–226. [Google Scholar] [CrossRef]
  89. Bóna, P. Geometric formulation of nonlinear quantum mechanics for density matrices. arXiv 1999, arXiv:quant-ph/9910011. [Google Scholar]
  90. Gheorghiu-Svirschevski, S. Nonlinear quantum evolution with maximal entropy production. Phys. Rev. A 2001, 63, 022105. [Google Scholar] [CrossRef] [Green Version]
  91. Gheorghiu-Svirschevski, S. Addendum to “Nonlinear quantum evolution with maximal entropy production.”. Phys. Rev. A 2001, 63, 054102. [Google Scholar] [CrossRef]
  92. Gheorghiu-Svirschevski, S. Quantum nonlocality and quantum dynamics. arXiv 2002, arXiv:quant-ph/0203153. [Google Scholar]
  93. Gheorghiu-Svirschevski, S. A general framework for nonlinear quantum dynamics. arXiv 2002, arXiv:quant-ph/0207042. [Google Scholar]
  94. Beretta, G.P. Nonlinear model dynamics for closed-system, constrained, maximal-entropy-generation relaxation by energy redistribution. Phys. Rev. E 2006, 73, 026113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Beretta, G.P. Time-energy and time-entropy uncertainty relations in dissipative quantum dynamics. arXiv 2006, arXiv:quant-ph/0511091. [Google Scholar]
  96. Rigolin, G. A simple derivation of the Schrödinger uncertainty relation. Eur. J. Phys. 2015, 36, 065007. [Google Scholar] [CrossRef]
  97. Schrödinger, E. About Heisenberg Uncertainty Relation. Proc. Prussian Acad. Sci. Phys. Math. 1930, 19, 296–303. [Google Scholar]
  98. Robertson, H.P. A general formulation of the uncertainty principle and its classical interpretation. Phys. Rev. 1930, 35, 667. [Google Scholar] [CrossRef]
  99. Robertson, H.P. An indeterminacy relation for several observables and its classical interpretation. Phys. Rev. 1934, 46, 794. [Google Scholar] [CrossRef]
  100. Beretta, G.P. Maximum entropy production rate in quantum thermodynamics. J. Phys. Conf. Ser. 2010, 237, 012004. [Google Scholar] [CrossRef]
  101. Eberly, J.H.; Singh, L.P.S. Time operators, partial stationarity, and the energy-time uncertainty relation. Phys. Rev. D 1973, 7, 359. [Google Scholar] [CrossRef]
  102. Leubner, C.; Kiener, C. Improvement of the Eberly-Singh time-energy inequality by combination with the Mandelstam-Tamm approach. Phys. Rev. A 1985, 31, 483. [Google Scholar] [CrossRef]
  103. Bhattacharyya, K. Quantum decay and the Mandelstam-Tamm-energy inequality. J. Phys. A Math. Gen. 1983, 16, 2993. [Google Scholar] [CrossRef]
  104. Pfeifer, P. How fast can a quantum state change with time? Phys. Rev. Lett. 1993, 70, 3365. [Google Scholar] [CrossRef] [PubMed]
  105. Pfeifer, P. Erratum: How fast can a quantum state change with time? Phys. Rev. Lett. 1993, 71, 306. [Google Scholar] [CrossRef]
  106. Massieu, F. Sur les fonctions caractéristiques des divers fluides et sur la théorie des vapeurs. C. R. Acad. Sci. 1859, 69, 858–862. [Google Scholar]
  107. Beretta, G.P. Well-behaved nonlinear evolution equation for steepest-entropy-ascent dissipative quantum dynamics. Int. J. Quantum Inf. 2007, 5, 249–255. [Google Scholar] [CrossRef]
  108. Rao, R.; Esposito, M. Conservation laws shape dissipation. New J. Phys. 2018, 20, 023007. [Google Scholar] [CrossRef] [Green Version]
  109. Rao, R.; Esposito, M. Detailed fluctuation theorems: A unifying perspective. Entropy 2018, 20, 635. [Google Scholar] [CrossRef]
  110. Yunger Halpern, N.; Faist, P.; Oppenheim, J.; Winter, A. Microcanonical and resource-theoretic derivations of the thermal state of a quantum system with noncommuting charges. Nat. Commun. 2016, 7, 12051. [Google Scholar] [CrossRef]
  111. Bedeauxa, D.; Mazur, P. Mesoscopic non-equilibrium thermodynamics for quantum systems. Physica A 2001, 298, 81–100. [Google Scholar] [CrossRef]
  112. Beretta, G.P. A theorem on Lyapunov stability for dynamical systems and a conjecture on a property of entropy. J. Math. Phys. 1986, 27, 305308. [Google Scholar] [CrossRef]
  113. Beretta, G.P. Quantum thermodynamics of nonequilibrium. Onsager reciprocity and dispersion-dissipation relations. Found. Phys. 1987, 17, 365–381. [Google Scholar] [CrossRef] [Green Version]
  114. Evans, D.J.; Searles, D.J. Equilibrium microstates which generate second law violating steady states. Phys. Rev. E 1994, 50, 1645. [Google Scholar] [CrossRef]
  115. Crooks, G.E. Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences. Phys. Rev. E 1999, 60, 2721. [Google Scholar] [CrossRef]
  116. Dewar, R. Information theory explanation of the fluctuation theorem, maximum entropy production and self-organized criticality in non-equilibrium stationary states. J. Phys. A Math. Gen. 2003, 36, 631. [Google Scholar] [CrossRef]
  117. Van Zon, R.; Cohen, E.G.D. Extension of the fluctuation theorem. Phys. Rev. Lett. 2003, 91, 110601. [Google Scholar] [CrossRef] [PubMed]
  118. Wang, G.M.; Sevick, E.M.; Mittag, E.; Searles, D.J.; Evans, D.J. Experimental demonstration of violations of the second law of thermodynamics for small systems and short time scales. Phys. Rev. Lett. 2002, 89, 050601. [Google Scholar] [CrossRef] [PubMed]
  119. Nieuwenhuizen, T.M.; Allahverdyan, A.E. Comment on “Experimental violations of the second law of thermodynamics for small systems and short timescales.”. arXiv 2002, arXiv:cond-mat/0207587. [Google Scholar]
  120. Garnier, N.; Ciliberto, S. Nonequilibrium fluctuations in a resistor. Phys. Rev. E 2005, 71, 060101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Ion, D.B.; Ion, M.L.D. Entropic lower bound for the quantum scattering of spinless particles. Phys. Rev. Lett. 1998, 81, 5714. [Google Scholar] [CrossRef]
  122. Deutsch, D. Uncertainty in quantum measurements. Phys. Rev. Lett. 1983, 50, 631. [Google Scholar] [CrossRef]
  123. Partovi, M.H. Entropic formulation of uncertainty for quantum measurements. Phys. Rev. Lett. 1983, 50, 1883. [Google Scholar] [CrossRef]
  124. Maassen, H.; Uffink, J.B.M. Generalized entropic uncertainty relations. Phys. Rev. Lett. 1988, 60, 1103. [Google Scholar] [CrossRef] [PubMed]
  125. Simmons, R.F.; Park, J.L. On completely positive maps in generalized quantum dynamics. Found. Phys. 1981, 11, 47–55. [Google Scholar] [CrossRef]
  126. Prigogine, I.; Mayne, F.; George, C.; De Haan, M. Microscopic theory of irreversible processes. Proc. Natl. Acad. Sci. USA 1977, 74, 4152–4156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  127. Grecos, A.; Prigogine, I. On the derivation of linear irreversible thermodynamics for classical fluids. Proc. Natl. Acad. Sci. USA 1978, 75, 1632–1636. [Google Scholar] [CrossRef] [Green Version]
  128. Misra, B.; Prigogine, I.; Courbage, M. Lyapounov variable: Entropy and measurement in quantum mechanics. Proc. Natl. Acad. Sci. USA 1979, 76, 4768–4772. [Google Scholar] [CrossRef] [Green Version]
  129. Beretta, G.P. Intrinsic entropy and intrinsic irreversibility for a single isolated constituent of matter: Broader kinematics and generalized nonlinear dynamics. In Frontiers of Nonequilibrium Statistical Physics, Proceedings of the NATO Advanced Study Institute on the Frontiers of Nonequilibrium Statistical Physics, Santa Fe, New Mexico, 3–16 June 1984; Moore, G.T., Scully, M.O., Eds.; NATO ASI Series B; Springer: Boston, MA, USA, 1986; Volume 135, pp. 205–212. [Google Scholar] [CrossRef]
  130. Beretta, G.P. What if quantum thermodynamics were a fundamental extension of quantum mechanics? In Perimeter Institute Recorded Seminar Archive; PIRSA Quantum Foundations Series; Perimeter Institute of Theoretical Physics: Waterloo, ON, Canada, 2007. [Google Scholar]
  131. Beretta, G.P. Mechanics and Thermodynamics Can be Fundamentally United by Density Operators with an Ontic Status Obeying a Locally Maximum Entropy Production Dynamics. But at What Price? In Perimeter Institute Recorded Seminar Archive; PIRSA Collection: PIAF 09, New Perspectives on the Quantum State; Perimeter Institute of Theoretical Physics: Waterloo, ON, Canada, 2009. [Google Scholar]
  132. Hiai, F.; Ohya, M.; Tsukada, M. Sufficiency, KMS condition and relative entropy in von Neumann algebras. Pac. J. Math. 1981, 96, 99–109. [Google Scholar] [CrossRef] [Green Version]
  133. Beretta, G.P.; Zanchini, E. New definitions of thermodynamic temperature and entropy not based on the concepts of heat and thermal reservoir. Atti della Accademia Peloritana dei Pericolanti-Classe di Scienze Fisiche Matematiche e Naturali 2019, 97 (Suppl. 1), A1. [Google Scholar] [CrossRef]
  134. Sieniutycz, S. From a least action principle to mass action law and extended affinity. Chem. Eng. Sci. 1987, 42, 2697–2711. [Google Scholar] [CrossRef]
  135. Beretta, G.P. Quantum thermodynamics: Microscopic foundations of entropy and of entropy generation by Iirreversibility. Atti della Accademia Peloritana dei Pericolanti Classe di Scienze Fisiche Matematiche e Naturali 2008, 86, 1–22. [Google Scholar] [CrossRef]
  136. Beretta, G.P. Nonlinear dynamical equation for irreversible, steepest-entropy-ascent relaxation to stable equilibrium. In Quantum Theory: Reconsideration of Foundations—4, Proceedings of the IV International Conference on Quantum Theory: Reconsideration of Foundations, Vaxjo, Sweden, 11–16 June 2007; Adenier, G., Khrennikov, A.Y., Lahti, P., Man’ko, V.I., Nieuwenhuizen, T.M., Eds.; American Institute of Physics: Melville, NY, USA, 2007; Volume 962, pp. 233–237. [Google Scholar] [CrossRef]
  137. Hatsopoulos, G.N.; Beretta, G.P. Where is the entropy challenge? In Meeting the Entropy Challenge, Proceedings of the International Thermodynamics Symposium in Honor and Memory of Professor Joseph H. Keenan, MIT, Cambridge, MA, USA, 4–5 October 2007; Beretta, G.P., Ghoniem, A.F., Hatsopoulos, G.N., Eds.; American Institute of Physics: Melville, NY, USA, 2008; Volume 1033, pp. 34–54. [Google Scholar] [CrossRef]
  138. Beretta, G.P. Steepest entropy ascent model for far-non-equilibrium dissipative evolution in quantum thermodynamics: A 1984 precursor of GENERIC (1997), gradient flows (1998), maximal entropy production (2001), and SEAQT (2014). In Proceedings of the Workshop on Mathematical Aspects of Non-Equilibrium Thermodynamics, RWTH, Aachen University, Aachen, Germany, 5–7 March 2018. [Google Scholar]
  139. Beretta, G.P. A new approach to constrained-maximization nonequilibrium problems. In Computer-Aided Engineering of Energy Systems: Second Law Analysis and Modeling, Proceedings of the Winter Annual Meeting of the ASME, Anaheim, CA, USA, 7–12 December 1986; Gaggioli, R.A., Ed.; Book H0341C-AES; ASME: New York, NY, USA, 1986; Volume 3, pp. 129–134. [Google Scholar]
  140. Beretta, G.P. Steepest-ascent constrained approach to maximum entropy. In Second Law Analysis of Heat Transfer in Energy Systems, Proceedings of the Winter Annual Meeting of the ASME, Boston, MA, USA, 13–18 December 1987; Boehm, R.F., Lior, N., Eds.; Book G00390; ASME: New York, NY, USA, 1987; Volume HTD, pp. 31–38. [Google Scholar]
  141. Beretta, G.P. Dynamics of smooth constrained approach to maximum entropy. In Second Law Analysis of Thermal Systems, Proceedings of the 4th International Symposium on Second Law Analysis of Thermal Systems, Roma, Italy, 25–29 May 1987; Moran, M.J., Sciubba, E., Eds.; Book I00236; ASME: New York, NY, USA, 1987; pp. 17–24. [Google Scholar]
  142. Beretta, G.P. Axiomatic definition of entropy for nonequilibrium states. In Proceedings of the 9th Joint European Thermodynamics Conference JETC07, Saint Etienne, France, 12–15 June 2007. [Google Scholar]
  143. Beretta, G.P. The second law from locally maximal entropy generation quantum dynamics. In Meeting the Entropy Challenge, Proceedings of the International Thermodynamics Symposium in Honor and Memory of Professor Joseph H. Keenan, MIT, Cambridge, MA, USA, 4–5 October 2007; Beretta, G.P., Ghoniem, A.F., Hatsopoulos, G.N., Eds.; American Institute of Physics: Melville, NY, USA, 2008; Volume 1033, pp. 180–187. [Google Scholar] [CrossRef]
  144. Beretta, G.P.; Hadjiconstantinou, N.G. Steepest entropy ascent models of the Boltzmann equation. Comparisons with hard-sphere dynamics and relaxation-time models for homogeneous relaxation from highly non-equilibrium states. In Proceedings of the ASME 2013 International Mechanical Engineering Congress and Exposition IMECE2013, San Diego, CA, USA, 15–21 November 2013. [Google Scholar] [CrossRef]
  145. Kaufman, A.N. Dissipative hamiltonian systems: A unifying principle. Phys. Lett. A 1984, 100, 419–422. [Google Scholar] [CrossRef] [Green Version]
  146. Morrison, P.J. Bracket formulation for irreversible classical fields. Phys. Lett. A 1984, 100, 423–427. [Google Scholar] [CrossRef]
  147. Grmela, M. Bracket formulation of dissipative fluid mechanics equations. Phys. Lett. A 1984, 102, 355–358. [Google Scholar] [CrossRef]
  148. Jordan, R.; Kinderlehrer, D.; Otto, F. The variational formulation of the Fokker–Planck equation. SIAM J. Math. Anal. 1998, 29, 1–17. [Google Scholar] [CrossRef]
  149. Otto, F. The geometry of dissipative evolution equations: The porous medium equation. Commun. Partial Differ. Equ. 2001, 26, 101–174. [Google Scholar] [CrossRef]
  150. Mielke, F. A gradient structure for reaction–diffusion systems and for energy-drift-diffusion systems. Nonlinearity 2011, 24, 1329–1346. [Google Scholar] [CrossRef]
  151. Mielke, A.; Peletier, M.A.; Renger, D.R.M. On the relation between gradient flows and the large-deviation principle, with applications to Markov chains and diffusion. Potential Anal. 2014, 41, 1293–1327. [Google Scholar] [CrossRef]
  152. Montefusco, A.; Peletier, M.A.; Öttinger, H.C. Coarse-graining via the fluctuation-dissipation theorem and large-deviation theory. arXiv 2018, arXiv:1809.07253. [Google Scholar]
Figure 1. (color online) Time-dependent relaxation results obtained by integrating the steepest-entropy- ascent master Equation (89) for the four-level qudit with equally spaced energy levels, for two different choices of τ : (a) τ = const and (b) τ state-dependent according to Equation (95) and time non-dimensionalized by / u where u is the energy difference between the highest and lowest energy levels of the system. First row subfigures: Time evolution of the four occupation probabilities p n . Second row: dimensionless entropy S / k B . Third row: rate of entropy change (proportional to σ M M ) compared with σ S S and σ H H / θ 2 , to illustrate relation (73). Fourth row: generalized ‘nonequilibrium temperature’ θ (nondimensionalized by u / k B ) compared with Δ H / Δ S and 2 Δ M Δ H / σ S S (also nondimensionalized) to illustrate relations (78) and (79). Fifth row: characteristic time of purely dissipative evolution τ D (here proportional to the inverse of the square root of the rate of entropy generation, shown in the third row subplots) compared with τ S and τ K to illustrate relations (53) and (76). Sixth row: characteristic times of the four occupation probabilities τ P e n compared with τ D to illustrate relation (64).
Figure 1. (color online) Time-dependent relaxation results obtained by integrating the steepest-entropy- ascent master Equation (89) for the four-level qudit with equally spaced energy levels, for two different choices of τ : (a) τ = const and (b) τ state-dependent according to Equation (95) and time non-dimensionalized by / u where u is the energy difference between the highest and lowest energy levels of the system. First row subfigures: Time evolution of the four occupation probabilities p n . Second row: dimensionless entropy S / k B . Third row: rate of entropy change (proportional to σ M M ) compared with σ S S and σ H H / θ 2 , to illustrate relation (73). Fourth row: generalized ‘nonequilibrium temperature’ θ (nondimensionalized by u / k B ) compared with Δ H / Δ S and 2 Δ M Δ H / σ S S (also nondimensionalized) to illustrate relations (78) and (79). Fifth row: characteristic time of purely dissipative evolution τ D (here proportional to the inverse of the square root of the rate of entropy generation, shown in the third row subplots) compared with τ S and τ K to illustrate relations (53) and (76). Sixth row: characteristic times of the four occupation probabilities τ P e n compared with τ D to illustrate relation (64).
Entropy 21 00679 g001
Figure 2. Time evolution of various other ratios of characteristic times for the same cases of Figure 1. First row subfigures: Time evolution of occupation probabilities p n (same as first row of Figure 1, repeated here for ease of comparison). Second row: ratios τ D / τ P e n for each of the four occupation probabilities to illustrate again relation (64). Third row: τ K / τ P e n to illustrate relation (86). Fourth row: 2 τ D | p ˙ n | to illustrate relation (62). Fifth row: τ D / τ S and τ K / τ S to illustrate relations (53) and (84).
Figure 2. Time evolution of various other ratios of characteristic times for the same cases of Figure 1. First row subfigures: Time evolution of occupation probabilities p n (same as first row of Figure 1, repeated here for ease of comparison). Second row: ratios τ D / τ P e n for each of the four occupation probabilities to illustrate again relation (64). Third row: τ K / τ P e n to illustrate relation (86). Fourth row: 2 τ D | p ˙ n | to illustrate relation (62). Fifth row: τ D / τ S and τ K / τ S to illustrate relations (53) and (84).
Entropy 21 00679 g002

Share and Cite

MDPI and ACS Style

Beretta, G.P. Time–Energy and Time–Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations. Entropy 2019, 21, 679. https://doi.org/10.3390/e21070679

AMA Style

Beretta GP. Time–Energy and Time–Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations. Entropy. 2019; 21(7):679. https://doi.org/10.3390/e21070679

Chicago/Turabian Style

Beretta, Gian Paolo. 2019. "Time–Energy and Time–Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations" Entropy 21, no. 7: 679. https://doi.org/10.3390/e21070679

APA Style

Beretta, G. P. (2019). Time–Energy and Time–Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations. Entropy, 21(7), 679. https://doi.org/10.3390/e21070679

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop