arxiv: v1 [math.fa] 1 Oct 2018

Transcription

1 GRADIENT FLOW STRUCTURE AND EXPONENTIAL DECAY OF THE SANDWICHED RÉNYI DIVERGENCE FOR PRIMITIVE LINDBLAD EQUATIONS WITH GNS-DETAILED BALANCE arxiv:18196v1 [mathfa] 1 Oct 18 YU CAO, JIANFENG LU, AND YULONG LU Abstract We study the entropy production of the sandwiched Rényi divergence under a primitive Lindblad equation with GNS-detailed balance We prove that the Lindblad equation can be identified as the gradient flow of the sandwiched Rényi divergence for any order, This extends a previous result by Carlen and Maas [Journal of Functional Analysis, 735, ] for the quantum relative entropy ie, = 1 Moreover, we show that the sandwiched Rényi divergence with order, decays exponentially fast under the time-evolution of such a Lindblad equation Keywords: gradient flow dynamics; exponential convergence; sandwiched Rényi divergence; Lindblad equation 1 Introduction and main result The characterization of dynamics for open classical and quantum systems is an essential and fundamental topic One approach to characterize the dynamics for open systems is entropy production [7, 1, 36, 46], from the perspective of information theory For Markovian dynamics, the entropy production is positive, which relates to the data processing inequality [, 31, 35] and complies with the second law of thermodynamics [8, 17] The entropy production is also related with the gradient flow structure on well-chosen Riemannian manifolds for those Markov semigroups [14, 19,, 8, 3, 41] Furthermore, the validity of log-sobolev inequalities [6, 4, 43] helps to prove the exponential decay of those entropies under respective Markov semigroups In this paper, we focus on the entropy production of sandwiched Rényi divergences, under primitive Lindblad equations with GNS-detailed balance, for finitedimensional quantum systems Let n be the dimension of the quantum system The Lindblad equation is an ordinary differential equation ODE ρ t = L ρ t, where ρ t are n-by-n density matrices, and the Lindblad super-operator L is a linear operator acting on the space of density matrices The Lindblad equation with GNS-detailed balance see Theorem 1 below has the form 1 ρ t = L ρ t = J j=1 e ωj/ [ V j ρ t, V j ] [ + Vj, ρ t Vj ], Date: October 3, 18 1

2 YU CAO, JIANFENG LU, AND YULONG LU where positive integer J n 1; ω j R; V j are n-by-n matrices; the constraints for ω j and V j will be given in Theorem 1 below; the commutator of matrices A and B is defined as [A, B] := AB BA The notion GNS-detailed balance will be explained in Sec below The Lindblad equation is called primitive also sometimes referred as ergodic, eg, see [1, 14] if KerL is spanned by I, which is the n-by-n identity matrix Before delving more into Lindblad equations for open quantum systems, we would like to review some existing results of Fokker-Planck equations for open classical systems, for the purpose of comparison In fact, almost all the conclusions in this paper are motivated by the classical results for Fokker-Planck equations Let us consider, for instance, a n-dimensional Brownian particle under an external potential Ψx with x R n Then the time evolution of the probability density function p t x of the Brownian particle is described by the Fokker-Planck equation t p t x = p t x Ψx + p t x It has been shown by Jordan, Kinderlehrer, and Otto [8] that such types of Fokker- Planck equations can be viewed as the Wasserstein gradient flow dynamics of the relative entropy or the free energy functional Moreover, based on the method of entropy production and the log-sobolev inequality [4], one can prove that the solution of the Fokker-Planck equation converges to the Gibbs stationary distribution exponentially fast in the relative entropy In our recent work [11], we generalized these results to the classical Rényi divergence for any order, : 1 the same Fokker-Planck equation can be identified as a gradient flow dynamics of the classical Rényi divergence of any order with respect to a transportation distance defined by a certain metric tensor; the solution of the Fokker-Planck equation converges exponentially fast to the Gibbs distribution in the classical Rényi divergence Apart from Fokker-Planck equations, several other dissipative dynamics can also be viewed as the gradient flow dynamics of various transport metrics, for instance, the heat equation [19], finite Markov Chains [3] and porous medium equations [, 41] As for open quantum systems, Lindblad equations [5, 3], generally regarded as the quantum analog of Fokker-Planck equations, could be derived from the weakcoupling limit of the system and baths [9,18] Due to the similar role that Lindblad and Fokker-Planck equations play for open quantum and classical systems respectively, it is interesting to study Lindblad equations from the perspectives of gradient flows In a previous work by Carlen and Maas, they showed that fermionic Fokker- Planck equation can be identified as the gradient flow dynamics of the quantum relative entropy [13]; in their more recent work, they extended results from [13] and showed that for finite dimensional quantum systems, primitive Lindblad equations with GNS-detailed balance can indeed be identified as the gradient flow dynamics of the quantum relative entropy [14] The underlying quantum Wasserstein distance for density matrices is defined via Benamou-Brenier formalism [6] The generalization of their results from finite dimensional matrices to infinite dimensional noncommutative algebras has been recently considered in [7, 5] In the present paper, we focus on the finite dimensional setting and aim to characterize the gradient flow structure of Lindblad equations in terms of a family of quantum divergences sandwiched Rényi divergences [38, 49], including the quantum relative entropy as a specific instance Moreover, we are also interested in quantifying

3 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 3 the decaying property of the solution ρ t of the Lindblad equation 1 in terms of the sandwiched Rényi divergence The sandwiched Rényi divergence is defined as follows Definition 11 Sandwiched Rényi divergence [38] Given two positive semidefinite matrices ρ, such that ρ and ρ, the sandwiched Rényi divergence D ρ of order, is defined by 3 D ρ := 1 1 log Tr [ 1 1 ] ρ,, 1 1, ; Tr ρ logρ ρ log, = 1 The notation ρ means the kernel of is a subspace of the kernel of ρ The term Trρ in [38, Definition ] is simply set as one, since we only consider ρ to be a density matrix in this paper is assumed to be a full-rank density matrix, so ρ is always satisfied for this paper For fixed ρ and full-rank, the sandwiched Rényi divergence is continuous with respect to order, in particular, D ρ = lim 1 D ρ The sandwiched Rényi divergence is a quantum analog of the classical Rényi divergence, and has recently received much attention due to its use in quantum information theory; see eg, [33, 49] The sandwiched Rényi divergence unifies many entropy measures: when = 1, it reduces to the quantum relative entropy; when, it converges to the quantum relative max-entropy; when = 1/ and =, it is closely connected to the quantum fidelity and χ divergence [47], respectively The sandwiched Rényi divergence satisfies the data processing inequality for order 1 for any completely positive trace-preserving CPTP map [4,,38] For fixed density matrices ρ and, it is monotonically increasing with respect to the order, [38, Theorem 7] The sandwiched Rényi divergence has also been used to study the quantum second laws of thermodynamics [8], which is a major motivation of our work Its global convergence rate and mixing time under quantum Markov semigroups has been studied in [34] We remark that there is another well-studied family of quantum Rényi divergences known as Petz-Rényi divergence [4], which is defined by D ρ := 1 1 log Tr ρ 1 When ρ and commute, the sandwiched Rényi divergence is identical to the Petz-Rényi divergence We also remark that [] has studied more generalized, z-quantum Rényi divergences, which include sandwiched Rényi divergences and Petz-Rényi divergences as special instances; the discussion on the data processing inequality of, z-quantum Rényi divergences can be found in a recent review paper [1] and the references therein The first main result of this paper is the following theorem Theorem 1 Consider a primitive Lindblad equation with GNS-detailed balance 1 and suppose its unique stationary state is full-rank For any order,, consider the energy functional D ρ with respect to ρ In the space of strictly positive density matrices D + equipped with the metric tensor defined in the Definition 36, the gradient flow dynamics of the sandwiched Rényi divergence D ρ is exactly the Lindblad equation 1

4 4 YU CAO, JIANFENG LU, AND YULONG LU This theorem immediately implies the monotonicity of the sandwiched Rényi divergence under the evolution of such Lindblad equations, summarized in the following Corollary 13 For 1, this corollary can be proved simply by applying the data processing inequality [] In our case, the monotonicity also holds for, 1 for primitive Lindblad equations with GNS-detailed balance Corollary 13 Under the same conditions as in Theorem 1, suppose ρ : [, D + is the solution of that Lindblad equation For any,, the sandwiched Rényi divergence D ρ t, is non-increasing with respect to time t, ie t D ρ t Next, we would like to study the exponential decay of the sandwiched Rényi divergence For Fokker-Planck equations in classical systems, the exponential decay of the relative entropy could be proved via the classical log-sobolev inequality, namely, the relative entropy is bounded above by the relative Fisher information This approach has been generalized in our previous work [11] to obtain the exponential decay of the classical Rényi divergence by introducing the relative -Fisher information For quantum systems, in the same flavor, let us first define the quantum relative -Fisher information as follows Definition 14 Suppose L is the generator of a primitive Lindblad equation with GNS-detailed balance and the unique stationary state is full-rank The quantum relative -Fisher information between ρ D + and is defined as { δdρ δρ, L ρ, if 1; 4 I ρ := logρ log, L ρ, if = 1 When = 1, we denote I ρ I 1 ρ and call it quantum relative Fisher information for simplicity The inner product, herein is the Hilbert-Schmidt inner product, defined by A, B := TrA B for any matrices A, B Notice that in the definition above, the quantum relative -Fisher information depends on the Lindblad super-operator L, which is different from the classical setting where the relative -Fisher information only depends on the Gibbs stationary distribution see [11] Also, observe that if ρ t solves the Lindblad equation 1, then I ρ t = t D ρ t Since by Corollary 13 t D ρ t, the quantum relative - Fisher information I ρ t is non-negative Definition 15 Suppose density matrix is the unique stationary state of the Lindblad super-operator L For a fixed order,, density matrix is said to satisfy the quantum -log Sobolev inequality or quantum -LSI in short if there exists K >, such that 5 D ρ 1 I ρ, ρ D, K where D is the space of density matrices K is called -log Sobolev constant When = 1, we call 5 quantum log-sobolev inequality and denote K K 1 for simplicity The -log Sobolev constant K has been considered in [34, Definition 31] When = 1, the above form reduces to the quantum log-sobolev inequality in [14, Eq 817]

5 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 5 Remark 16 Another definition of quantum -LSI was given previously under the framework of noncommutative L spaces see eg, [9, Definition 11] and [5, Sec 5] We will briefly revisit this concept in Sec 5 and we shall denote the log- Sobolev constants in this setting up by κ The essential difference is that for the definition of κ, the entropy function uses the quantum relative entropy for a density matrix raising to the power of, instead of using sandwiched Rényi divergence see 3 for more details For Lindblad equations with GNS-detailed balance, the primitivity of Lindblad equations is equivalent to the validity of the quantum LSI, summarized in the following Proposition 17 It is important to note that by Lemma 4 below, if the Lindblad equation satisfies GNS-detailed balance, then L is a positive semidefinite operator, with respect to the inner product, 1/ which is defined by A, B 1/ := Tr A 1/ B 1/ for any matrices A, B As a consequence, primitivity of L implies that L has a positive spectral gap, denoted by λ L The proposition below provides two-side estimates of K in terms of λ L Proposition 17 Consider the Lindblad equation with GNS-detailed balance 1 and suppose a full-rank density matrix is its stationary state satisfies the quantum LSI iff the Lindblad equation is primitive More quantitatively, if the Lindblad equation is primitive, we have 1 λ min log λ min 1 log logλmin 1 log1+e 1 where λ min is the smallest eigenvalue of λ L K λ L, Proposition 17 follows mainly from [34, Theorem 5], [34, Theorem 53], [9, Theorem 16], and the quantum Stroock-Varopoulos inequality [5, Corollary 16] The proof is postponed to Sec 43 It is not hard to see that if the quantum -LSI 5 is valid, then it follows immediately from Grönwall inequality that see also [34, Theorem 3] 6 D ρ t D ρ e Kt, t The validity of 5 with a positive K is in general difficult to establish For depolarizing semigroup, K has been computed in [34, Theorem 33] and K has been studied in [37] For depolarizing semigroup with stationary state being the maximally mixed state, both upper bound and lower bound of K with have been provided in [34, Theorem 34] Moreover, it has been proved in [34, Theorem 41] that the -log Sobolev constant can not exceed the spectral gap, namely for any 1, K λ L In the present paper, instead of pursuing the global decay of the sandwiched Rényi divergence D ρ t as in 6, we focus on the exponential decay of D ρ t in the large time regime More specifically, our second main result Theorem 18 shows that the sandwiched Rényi divergence D ρ t decays exponentially with the sharp rate λ L when t is large Theorem 18 Under the same conditions as in Theorem 1, for any given solution ρ : [, D + of the Lindblad equation and any ɛ, there, λmin

6 6 YU CAO, JIANFENG LU, AND YULONG LU exist C,ɛ and τ,ɛ such that 7 D ρ t C,ɛ D ρ e λ Lt, t τ,ɛ, where and C,ɛ = edρ 1 D ρ exp Θ λ L T,ɛ,, ; τ,ɛ = 1 T,ɛ + max, K logd ρ /ɛ, > In the above, λ min is the smallest eigenvalue of ; T,ɛ = 1 Kη log 1 is a function depending on and ɛ; Θ is the Heaviside function; η = min 1, J min j=1 Λ = λ max λ min exp e ωj 1 Λ, 1 + e ωj Λ λ min ɛ ɛ λ min λ min, ɛ where ω j are some constants determined by the stationary state see Sec The proof of Theorem 18 can be found in Sec 4 Here we sketch the strategies used in the proof First, we prove this theorem for the case = by using a uniform lower bound of I ρ see 5 Next, we prove that the quantum LSI 5 holds see Proposition 17 Based on the quantum LSI see 5 for = 1, we prove a quantum comparison theorem see Proposition 46, which implies the exponential decay of the sandwiched Rényi divergence for order 1, Finally, for the case <, the above result follows from the monotonicity of the sandwiched Rényi divergence with respect to order ; for the case >, the above result follows from the quantum comparison theorem Some remarks of this theorem are in order: Remark 19 1 The upper bounds for τ,ɛ and C,ɛ may not be sharp; a better bound might be important if one wants to study the mixing time for quantum decoherence, which however will not be pursued here For fixed >, when ɛ becomes smaller, T,ɛ decreases and we have a smaller prefactor C,ɛ, at the expense of waiting longer time τ,ɛ since the term D ρ /ɛ blows up as ɛ + The artificial ɛ and the complicated expression of T,ɛ come from Proposition 46 There is a different way to quantify η under the assumption of the quantum -LSI in the sense of noncommutative L spaces; please see Sec 5 for background and the discussion in Sec 46 for more details We prefer this version because the log- Sobolev inequality for = 1 is also used for the classical result in [11, Theorem 1] One could remove the dependence of the parameter ɛ by simply choosing, eg, ɛ = λ min /8 in this theorem, then Λ = λ max λ min e3, η = min 1, J min j=1 e 3/ ωj λmin e λ max ωj λmax 1 + e λ min e3

7 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 7 If we further assume that is the maximally mixed state, ie, = 1 ni, then Λ = e 3, are independent of the dimension n η = e 3/ 1 + e 3 Our contribution We summarize here the contributions of the current work, which consists of two parts: 1 We extend the work of Carlen and Maas in [14] by proving that primitive Lindblad equations with GNS-detailed balance can be identified as the gradient flow dynamics of the sandwiched Rényi divergence for any order, ; see Theorem 1 We prove that the sandwiched Rényi divergence decays exponentially fast under the evolution of primitive Lindblad equations with GNS-detailed balance; see Theorem 18 Our result 7 shows that the sandwiched Rényi divergence decays with an asymptotic rate λ L, which is sharp according to [34, Theorem 41] To prove Theorem 18, we first prove a quantum comparison theorem see Proposition 46, which is of interest in its own right It is esssentially a hypercontractivity-type estimate of the solution of the Lindblad equation, enabling us to bound D 1 ρ T by D ρ with some T > when 1 > As a direct consequence of the comparison theorem, Proposition 48 shows that the sandwiched Rényi divergence exponentially decays for all orders, Such a comparison result for classical Rényi divergences was proved for the Fokker-Planck equation in [11, Theorem 1] For quantum dynamics, a similar comparison result was also obtained in [5, Corollary 17] However, a major difference between our result and theirs is that we use K in the estimation of the waiting time, while they use κ see 8 The rest of the paper is organized as follows In Sec, we will recall several preliminary results on, eg, gradient flows, Lindblad equations with detailed balance, and noncommutative L spaces Sec 3 will be devoted to the proof of Theorem 1 and Corollary 13 Proposition 17 and Theorem 18 will be proved in Sec 4 Finally, we will discuss in Sec 5 some potential directions in the future, in particular, the generalization of quantum Wasserstein distance involving the sandwiched Rényi divergence Preliminaries In this section, we recall some technical notions: gradient flow dynamics on Riemannian manifold, Lindblad equations with GNS-detailed balance, quantum analog of gradient and divergence operator, chain rule identity, and quantum noncommutative L spaces Most results from Sec to Sec 4 follow from [14]; Sec 5 is mostly based on [5] Notations Throughout this section and the rest of this paper, we will use the following notations

8 8 YU CAO, JIANFENG LU, AND YULONG LU Denote M as the space of all n n complex matrices The space of Hermitian matrices is A; the space of strictly positive matrices is P + ; the space of traceless Hermitian matrices is A The space of density matrices is D and the space of invertible density matrices is D + Unity I is the n-by-n identity matrix; I is the identity super-operator on M Hilbert-Schmidt inner product is A, B := Tr A B for A, B M We extend the inner product, to vector fields over M by A, B := J j=1 A j, B j, where A = A 1 A A J, B = B1 B B J and Aj, B j M for all 1 j J; here J is a fixed positive integer The space of such vector fields is denoted by M J Such a convention is also applied analogously to all inner products on M Denote λ min ρ as the smallest eigenvalue of a density matrix ρ and λ max ρ as the largest one Asterisk in the superscript means Hermitian-conjugate; dagger means the adjoint operator with respect to the Hilbert-Schmidt inner product To simplify the notation, a weight operator Γ : M M is defined by 8 Γ A := 1 A 1, for any matrix A M Hence, for any γ R, Γ γ A = γ/ A γ/ 1 Gradient flow dynamics Consider a Riemannian manifold M, g ρ : M is a smooth manifold and g ρ, is an inner product on the tangent space T ρ M for any ρ M Given any energy functional E on M, the gradient, denoted by grade ρ at any state ρ M, is defined as the element in the tangent space T ρ M such that 9 g ρ grade, ν = d ρ dɛ ɛ= Eρ + ɛν, ν T ρ M Then, the gradient flow dynamics refers to following autonomous differential equation d 1 dt ρ t = grade ρt This gradient flow dynamics is essentially a generalization of steepest descend dynamics in the Euclidean space In this paper, we will take M to be the space of strictly positive density matrices D + The choice of g ρ is more subtle and technical; thus we will explain it below in Sec 3 Quantum Markov semigroup with GNS-detailed balance For a continuous time-parameterized semigroup P t t, acting on the space of linear operators on a finite dimensional Hilbert space C n, it is called a quantum Markov semigroup QMS if P t is completely positive and P t I = I, for any t The generator of the QMS is denoted by L, ie, P t = e tl Physically, this semigroup P t usually refers to the time evolution super-operator for observables in the Heisenberg picture

9 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 9 To explain quantum detailed balance, let us define inner products, s parameterized by s [, 1], for a given full-rank density matrix, in the following way: for any A, B M, 11 A, B s := Tr s A 1 s B Essentially,, s is a -weighted Hilbert-Schmidt inner product The QMS is said to satisfy GNS-detailed balance if P t is self-adjoint with respect to the inner product, 1 for all t, that is to say, P t A, B 1 = A, P t B 1 for all A, B M and t In some literature, the case s = 1/ is considered instead and the quantum Markov semigroup is said to satisfy KMS-detailed balance if P t is self-adjoint with respect to the inner product, 1/ It is not hard to verify that the QMS satisfies KMS-detailed balance iff 1 Γ 1 L Γ = L An important relation is that GNS-detailed balance is a more restrictive situation and it implies KMS-detailed balance; reversely, it is not true More discussion on the relation between GNS and KMS detailed balance could be found in [14], in particular, Theorem 9 therein Another important consequence is that if P t satisfies detailed balance either KMS or GNS, then is invariant under dynamical evolution P t, that is to say, P t = We add the prefix GNS and KMS to explicitly distinguish these two conditions The general form of quantum Markov semigroups with GNS-detailed balance is given by the theorem below Theorem 1 [14, Theorem 31] Suppose that the QMS P t satisfies GNS-detailed balance, then 13 LA = = J j=1 J j=1 e ωj/ Vj [A, V j ] + e ωj/ [V j, A] Vj e ωj/ V j [A, V j ] + [ V j, A ] V j, where ω j R and J n 1 Moreover {V j } J j=1 satisfy the following: 1 TrV j = and V j, V k = δ j,k c j for all 1 j, k J; constants c j > for all 1 j J come from normalization; for each j, there exists 1 j J such that V j = V j ; 3 V j = e ωj V j, where the modular operator : M M is defined as A := A 1, for A M; 4 if V j = V j then c j = c j and ω j = ω j

10 1 YU CAO, JIANFENG LU, AND YULONG LU Then it could be easily derived that see also Eq 37 in [14] 14 L A = = J j=1 J j=1 e ωj/ [ V j A, V j e ωj/ [ V j A, V j which is well known as the Lindblad super-operator ] + e ω j/ [ V j, AV j ] ] [ + Vj, AVj ], Example Consider the depolarizing semigroup with generator L depol A = γ TrA A, and γ > It is not hard to verify that L depol A = γ TrAI A and that L depol is self-adjoint with respect to the inner product, 1 Thus this depolarizing semigroup satisfies the GNS-detailed balance Lemma 3 Suppose QMS P t satisfies the GNS-detailed balance, then 15 [P t, ] =, and equivalently [L, ] = 3 Noncommutative analog of gradient and divergence Given those V j from 13, define noncommutative partial derivative j by 16 j A := [V j, A], for matrix A M; as a result, its adjoint operator j has the form j A = [ V j, A ] Let M J be the vector fields over M, ie, each element A M J has the form A = A 1 A A J where Aj M for all 1 j J The noncommutative gradient : M M J is defined by 17 A := 1 A A J A M J The noncommutative divergence, acting on M J, is defined by 18 div A = J j A j = j=1 J j=1 [ Aj, Vj ], for A = A 1 A A J M J The divergence is defined as this so that it is the adjoint operator of gradient, with a multiplicative factor 1, analogous to the classical case Lemma 4 [14, Eq 53] For a quantum Markov semigroup with GNS-detailed balance, 19 A, A 1/ = A, LA 1/ Thus L is a positive semi-definite operator on the Hilbert space M,, 1/ Lemma 5 [14, Theorem 53] Suppose the Lindblad equation satisfies GNSdetailed balance Then the kernel of L is equal to the commutant of {V j } J j=1 In other words, KerL = Ker

11 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 11 The lemma below gives a sufficient condition for Lindblad equations with GNSdetailed balance to be primitive, though the analysis in the rest of this manuscript does not rely on this result However, the condition J = n 1 is not necessary for Lindblad equations with GNS-detailed balance to be primitive For example, the case n =, = 1 I, V 1 = X and V = Z with J = still leads to a primitive Lindblad equation, though J < n 1; X and Z here are just Pauli matrices Lemma 6 If the Lindblad equation satisfies the GNS-detailed balance and J = n 1, then the Lindblad equation is primitive Proof Note that {I, V 1, V, V J } form a non-normalized basis of M Hence for any matrix A M, one can decompose A = TrA n I + à where Trà = By the definition of primitive Lindblad equation and Lemma 5, it is equivalent to show that the commutant of {V j } J j=1 is spanned by I, that is, if A is in the commutant of {V j } J j=1, then à = ] If A is in the commutant, then [Ã, Vj = for all j Since {V j } J j=1 spans the space ] of traceless matrices, = [Ã, j k for all j k, that is, à j k = j k à Considering the matrix element a b, we have Ãa,jδ k,b = Ãk,bδ a,j Case 1: take a = b = k, then we know that Ãk,j =, hence off-diagonal terms of à are all zero Case : take a = j and b = k, then we know that Ãj,j = Ãk,k; since = Trà = j Ãj,j, then Ãj,j = for all j In summary, à =, thus completes the proof 4 Chain rule identity One of the key steps in connecting the Lindblad equation and the gradient flow dynamics of the quantum relative entropy is the following chain rule identity In Sec 3 below, this chain rule identity is also the key in proving Lemma 31 Lemma 7 Chain rule identity [14, Lemma 55] For any V M, X P + and ω R, [X] ω V loge ω/ X loge ω/ XV = e ω/ V X e ω/ XV, where the operator [X] ω : M M is defined by 1 [X] ω A := 1 e ωs 1/ X s AX 1 s ds This operator [X] ω is a noncommutative multiplication of the operator X It can be extended and defined for vector fields M J For A = A 1 A A J M J and ω = ω 1 ω ω J R J, define the operator [X] ω : M J M J by [X] ω A := [X] ω1 A 1 [X] ω A [X] ωj A J As a reminder, these ω j with 1 j J come from the spectrum of the modular operator see Theorem 1 Lemma 8 If X P +, then for any ω R, [X] ω is a strictly positive operator on M As a consequence, [X] 1 ω is a well-defined positive operator

12 1 YU CAO, JIANFENG LU, AND YULONG LU Proof Notice that A, [X] ω A = = 1 1 e ωs 1/ TrA X s AX 1 s ds e ωs 1/ X s/ AX 1 s/, X s/ AX 1 s/ ds Moreover, A, [X] ω A = iff A = Thus [X] ω is a strictly positive operator on M for any ω R This completes the proof More explicitly, the inverse of [X] ω is see [14, Lemma 58] 3 [X] 1 ω A = 1 t + e ω/ X A 1 t + e ω/ X dt By the expression 1 and 3, it is straightforward to verify the following Lemma Lemma 9 [14, Lemma 58] For any ω R and A M, 4 [X] ω A = [X] ω A, [X] 1 ω A = [X] 1 ω A 5 Non-commutative L spaces In this subsection, we briefly recall some concepts from non-commutative L spaces and the log-sobolev constants in that setting [5, 9, 39], We also recall an implication of the quantum Stroock-Varopoulos inequality [5], which compares the magnitudes of various log-sobolev constants Let us first introduce the weighted L norm, defined as Γ 1/ A, := Tr A 1/ 5 Then we need to introduce the power operator, for any, β R, A M, Γ I β, A := Γ 1 1 β β A When A is positive and commutes with, I β, A = A /β Two most important ingredients in non-commutative L spaces are entropy function and Dirichlet form -Entropy function for X P + is defined as [ ] Ent, X := Tr Γ 1/ X log Γ 1/ X 6 [ ] Tr X log X, log, Γ 1/ and the -Dirichlet form, for, 1 1,, is defined as X, 7 E,L X := 4 I, X, LX 1/, where 1/ + 1/ = 1 When = 1, the Dirichlet form is defined by taking the limit 1, ie, E 1,L X := lim 1 E,L X = 1 4 logγ X log, LX 1/ The -log Sobolev constant in this setting is defined as E,L X 8 κ κ L := inf X> Ent, X,

13 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 13 and thus the quantum -LSI refers to 9 κ Ent, X E,L X, X > In order to compare 5 and 9, for any matrix X >, let us introduce ρ := Γ X/ Tr Γ X, which is a strictly positive density matrix By inverting the operator Γ, X = Tr Γ X Γ 1 ρ By the fact that both E,L X and Ent, X are homogeneous with respect to X see, eg, Proposition 4 ii and Proposition 8 ii in [5], we have κ = inf ρ D + = inf ρ D + E,L Γ 1 ρ Ent, Γ 1 ρ 4 1 Γ 1 Γ 1 ρ ZD Γ 1 / 1 ρ, LΓ 1 ρ 1/, /Z [ where Z = Tr Γ 1 / ] ρ, which turns out to be the exponent in the sandwiched Rényi divergence 3 Provided that Lindblad equations satisfy KMSdetailed balance, by 1 and Lemma 31 see below, Z Γ 1 Γ 1 1 ρ, L ρ κ = inf ρ D + D Γ 1 / ρ /Z 3 δdρ 4 δρ, L ρ 4 = inf = inf I ρ ρ D + D ρ /Z ρ D + D ρ /Z Γ 1 / Γ 1 / As one might observe, the -Dirichlet form becomes the quantum relative -Fisher information up to a multiplicative constant, after the above transformation However, the denominator in the above equation is different from the sandwiched Rényi divergence, which makes two definitions of the quantum -LSI in 5 and 9 different When = 1, this log-sobolev constant κ 1 becomes 31 κ 1 = 1 4 inf I ρ ρ D + D ρ = 1 4 K = K Therefore, both definitions 5 and 9 coincide when = 1 The following implication of the quantum Stroock-Varopoulos inequality see [5, Theorem 14] characterizes the relationship between various κ Theorem 1 [5, Corollary 16] Consider primitive Lindblad equations with GNS-detailed balance Then κ is non-increasing with respect to [, ] In particular, 3 κ 1 κ

14 14 YU CAO, JIANFENG LU, AND YULONG LU 3 Lindblad equations as the gradient flow dynamics of the sandwiched Rényi divergence This section devotes to the proof of Theorem We will focus on the case that 1 since the case = 1 has been treated previously in [14] Alternatively, one may easily adapt the proof below to the case = 1 The proof follows from a sequence of lemmas whose proofs are postponed to the end of this section 31 Proof of Theorem 1 To simplify notations, for any density matrix ρ D, define 33 ρ := Γ 1 ρ 1 1 ρ Then the sandwiched Rényi divergence for 1 can then be rewritten as D ρ = 1 1 log Trρ The following lemma, connecting the functional derivative of the sandwiched Rényi divergence and Lindblad equations, plays an essential role in the proof of Theorem 1 Lemma 31 Functional derivative The functional derivative of the sandwiched Rényi divergence in the space D is 34 and moreover for ρ D + δd ρ δρ 35 [ρ],ωj j δd ρ δρ Tr [X],ω := = 1 Γ 1 ρ 1 Tr ρ = e ωj/ V j ρ e ωj/ ρv j, where the operator [X],ω : M M, for positive X P + and ω R, is defined by 36 Γ 1 X [ [Γ 1 Γ 1 Γ 1 ] X ω/, X 1 ] 1 1ω/ Γ 1 Remark 3 1 When = 1, Γ 1 = I and ρ = ρ Then [ρ],ω reduces to [ρ] ω Hence, [ρ],ω is the generalization of [ρ] ω to the case of sandwiched Rényi divergence When =, the above operator [ρ],ω has a simple form [ρ],ω = Trρ Γ = Tr 1/ ρ 1/ ρ 37 Γ Thus the mapping ρ, 1/ ρ 1/ A, 1/ A 1/ 38 A, ρ A, [ρ],ω A = is the multiplication of quadratic terms with respect to both A and ρ The properties of the operator [X] ω stated in Lemma 8 and Lemma 9 generalize to the operator [X],ω, as we summarize in the next lemma

15 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 15 Lemma 33 If X P + and ω R, then 1 [X],ω : M M is strictly positive; For any matrix A, [X],ω A = [X], ω A The operator [X],ω can be similarly extended to the space M J For vector fields A = A 1 A A J M J and ω = ω 1 ω ω J R J, we define 39 [X], ω A := [X],ω1 A 1 [X],ω A [X],ωJ A J, and analogously, 4 [X] 1, ω A := [X] 1,ω 1 A 1 [X] 1,ω A [X] 1,ω J A J The following lemma will be useful to characterize the tangent space of density matrices Lemma 34 Assume that the Lindblad equation 1 is primitive and X P + Denote the space of traceless matrices by M and denote the space of traceless Hermitian matrices by A, that is, A = M A 1 Given ω R J, define the operator D,X, ω : M M by D,X, ω A := div [X], ω A This operator D,X, ω is strictly positive, and thus invertible For any A M, D,X, ω A = D,X, ω A, and the restriction of D,X, ω on A is still strictly positive and invertible In the case = 1, this lemma has been proved in [45, Lemma 3] Before defining the metric tensor, we need to first introduce an inner product, weighted by [ρ], ω Definition 35 Inner product Given,, ρ D +, for vector fields A = A 1 A A J and B = B1 B B J, an inner product on vector fields is defined as J 41 A, B := A j, [ρ],ωj B j,ρ,l j=1 Given ρ D +, the tangent space is A For ν 1 and ν A, by Lemma 34 part, there is a unique U k A, such that 4 ν k = D,ρ, ω U k = div [ρ], ω U k, k = 1, Finally, we are ready to define the metric tensor we need

16 16 YU CAO, JIANFENG LU, AND YULONG LU Definition 36 Metric tensor Given ρ D +, for ν k T ρ D A, define metric tensor g,ρ,l by 43 g,ρ,l ν 1, ν := U 1, U,ρ,L where U k and ν k are linked via 4 J j=1 j U 1, [ρ],ωj j U, Proof of Theorem 1 Up to now, the Riemannian structure D +, g,ρ,l has been specified, as well as the energy functional Eρ D ρ Then, we will verify that the gradient flow dynamics defined in 1 is exactly the Lindblad equation in 1 To achieve that, we need to compute the gradient at any ρ D + Recall the notion of gradient in the Riemannian manifold 9, grade ρ is defined as the traceless Hermitian matrix such that for any differentiable trajectory ρ t D + with ρ at time t =, g,ρ,l grade, ρ t ρ t= d dt Eρ t d t= dt D ρ t t= By direct computation, d δd ρ dt D ρ t = t= δρ, ρ t ρ=ρ t= δd ρ = δρ, div [ρ ], ω U ρ=ρ = δd ρ δρ, [ρ ], ω U ρ=ρ = δd ρ δρ, U, ρ=ρ,ρ,l where we used the following fact that follows from Lemma 34: there exists a unique family of traceless Hermitian matrices {U t } t such that ρ t = div[ρ t ], ω U t After combing the last two equations, we have g,ρ,l grade, ρ t = ρ t= δd ρ δρ By the definition of metric tensor g,ρ,l, 43, grade δd ρ = D,ρ, ω ρ δρ ρ=ρ = div [ρ ], ω δd ρ δρ, U ρ=ρ ρ=ρ Consequently, gradient flow dynamics 1 is ρ t = grade = div [ρ t ], ω δd ρ δρ ρt ρ=ρt,ρ,l

17 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 17 Let us rewrite the term on the right hand side 44 div [ρ t ], ω δd ρ δρ ρ=ρt J ] 18 δd ρ = [[ρ t ],ωj j, Vj δρ 35 = j=1 j=1 ρ=ρt J [ e ωj/ V j ρ t e ωj/ ρ t V j, Vj 14 = L ρ t ] Thus the gradient flow dynamics is exactly the Lindblad equation in 1 Proof of Corollary 13 If ρ t is the solution of the Lindblad equation, then U = δdρ δρ Hence, ρ=ρ d dt D ρ t = δd ρ t= δρ, U ρ=ρ,ρ,l = δd ρ δρ, δd ρ ρ=ρ δρ ρ=ρ,ρ,l 3 Proof of Lemmas Here we present the proofs of lemmas we used in the previous subsection Proof of Lemma 31 For any differentiable curve ρ t D passing through ρ at time t =, that is, ρ = ρ, we have δd ρ, ρ d δρ dt D ρ t t= = Tr ρ 1 1 = Tr 1 Then it is straightforward to obtain Tr ρ ρ 1 ρ 1 1 ρ Tr ρ

18 18 YU CAO, JIANFENG LU, AND YULONG LU Next, we compute the partial derivative of the functional derivative of the sandwiched Rényi divergence In fact, 45 δd ρ j δρ = 1 Tr ρ = 1 Tr ρ 1 [V j, 1 = 1 Tr ρ Γ 1 1 Tr ρ Γ 1 = = Tr ρ Γ 1 1 ρ 1 1 ] Vj 1 ρ 1 e 1 ωj V j ρ 1 [ρ ] 1 1 ωj ρ 1 1 e 1 ωj ρ 1 V j V j loge 1 ωj ρ 1 Vj 1 1 loge 1 ωj ρ 1 V j [ ρ 1 ] 1 V j loge ω j ρ loge ω j ρ V j ωj To get the fifth line, we have used with X = ρ 1 Note that since ρ, D +, ρ, ρ 1 P + and ω = 1 ω j and V = V j Furthermore, applying again with X = ρ, ω = ω j / and V = V j, [ρ ] ωj/ V j log e ωj/ ρ log e ω j/ ρ Vj = e ωj/ V j ρ e ωj/ ρ V j = e ωj/ 1 1 Vj 1 = Γ 1 e ωj/ V j ρ e ωj/ ρv j ρ 1 e ω j/ 1 ρ 1 Vj 1 1 Therefore, 46 e ωj/ V j ρ e ωj/ ρv j = Γ 1 [ρ ] ωj/ V j log e ωj/ ρ log e ω j/ ρ Vj After plugging it back into 45, one could straightforwardly obtain 35 after arranging a few terms and 47 [ρ] 1,ω j := Tr ρ Γ 1 [ ρ 1 ] 1 ωj Γ 1 [ρ ] ωj/ 1 This becomes 36 after we invert this operator

19 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 19 Proof of Lemma 33 1 Since X P +, D +, then X := Γ 1 X is strictly 1 positive For any A M, let à := Γ A, then A, [X],ω A = TrX = TrX 1 = TrX = TrX 3 = TrX = TrX A, Γ 1 Ã, [X ] [ ω/ X 1 ] [X ] ω/ [ X 1 ] 1 e ω s 1 Tr à X s e ω s 1 Tr e ω s 1 Tr ] 1 1ω/ X 1 s X 1 s e ω s 1 X s ÃX 1 s 1ω/ Γ 1 1ω/ à [X ] 1 1 1ω/ à A X 1 s ds à X s X s [X ] 1 1 1ω/ à à X s [X ] 1 1, [ X 1 ] 1 X 1 s ds X s 1ω/ ÃX 1 s ds X s 1ω/ ÃX 1 s ds, since [ X 1 is a strictly positive operator by Lemma 8 Additionally, A, [X],ω A = if and only if = X s ÃX 1 s X s Γ 1 A X 1 s, which implies that A = For any matrix A, by Lemma 9, Tr X [X],ω A = = Tr X = Tr X = Tr X Γ 1 Γ 1 Γ 1 Γ 1 = [X], ω A [X ] ω/ [ X 1 Γ 1 1ω/ A ] 1 [X ] 1 1 [X ] ω/ Γ 1 1ω/ [X ] ω/ [ X 1 ] 1 1ω/ [X ] ω/ [ X 1 ] 1 A Γ 1 A Γ 1 1ω/ A Proof of Lemma 34 1 Consider a traceless matrix A M It is easy to verify that Tr D,X, ω A = From the definition of the divergence operator 18, we have J A, D,X, ω A = j A, [X],ωj j A j=1 Moreover, since [X],ωj is strictly positive by Lemma 33, we know that A, D,X, ω A = iff j A = for all 1 j J Thanks to Lemma 5 and the assumption that the Lindblad equation 1 is primitive, we have that A = since A is assumed to be traceless Therefore, D,X, ω is a strictly positive operator on the space of traceless matrices M

20 YU CAO, JIANFENG LU, AND YULONG LU We only need to prove that D,X, ω A = D,X, ω A The rest simply follows from the first part of the lemma Note that D,X, ω A = div [X], ω A = [ ] Vj, [X],ωj j A j = Vj [X],ωj V j A AV j [X],ωj V j A AV j Vj j Hence by Theorem 1 and Lemma 33, we have D,X, ω A = [X], ωj A Vj Vj A V j V j [X], ωj A Vj Vj A j = [X] A,ωj V j V j A Vj V j [X] A,ωj V j V j A j = D,X, ω A 4 Exponential decay of sandwiched Rényi divergences This section is devoted to proving Theorem 18, ie, the exponential decay of sandwiched Rényi divergences for primitive Lindblad equations with GNS-detailed balance We start by recalling the definition of a spectral gap and then we prove a Poincaré inequality see Proposition 41 With this Poincaré inequality, we derive a uniform lower bound of the quantum relative -Fisher information, which immediately implies Theorem 18 in the case = Next we prove the Proposition 17, which shows that for primitive Lindblad equations with GNS-detailed balance, the quantum LSI 5 holds for = 1 ie, there exists K > Then we prove a quantum comparison theorem see Proposition 46, which implies that the exponential decay rates of sandwiched Rényi divergences along the Lindblad equation are the same for all > 1 This together with the monotonicity of sandwiched Rényi divergences with respect to order concludes the proof of Theorem Spectral gap and Poincaré inequality Let us recall from Sec that if the Lindblad equation satisfies GNS-detailed balance, then it also satisfies KMSdetailed balance ie, L is self-adjoint with respect to the inner product, 1/ Moreover, by Lemma 4, L is a positive semi-definite operator with respect to the inner product, 1/ By the additional assumption that the Lindblad equation is primitive, we know I is the only eigenvector of L with respect to the eigenvalue zero Hence, spectral theory shows that there exists an orthonormal basis {L 1, L,, L n 1, I} such that 48 LL j = θ j L j, θ 1 θ θ n 1 > ; LI = The spectral gap of the operator L is λ L := θ n 1 It is worth remarking that {L 1, L,, L n 1, I} is an orthonormal basis in the Hilbert space M equipped with the inner product, 1/, but not necessarily with other inner products The following Poincaré inequality follows directly from the definition of the spectral gap

21 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 1 Proposition 41 Poincaré inequality Assume that the primitive Lindblad equation satisfies GNS-detailed balance 1 For any A M such that I, A 1/ = or equivalently, tra =, 49 A, A 1/ = A, LA 1/ λ L A, A 1/ Moreover, the equality can be achieved when A = L n 1 Proof For any A M such that I, A 1/ =, we know that A span{l j } n 1 j=1 Thus by the definition of the spectral gap, A, LA 1/ λ L A, A 1/ Moreover, when A = V n 1, the equality holds Then 49 follows immediately from 19 Remark 4 In fact, the above Poincaré inequality is a special case of a whole family of Poincaré inequalities see [44, p 19]: given any function ϕ :,,, we could consider C ϕ A, ϕ A 1/ C ϕ A, ϕ A 1 A, ϕ LA 1 A, ϕ LA 1/, for all A M such that TrA =, where ϕx := x 1/ ϕx The inequality in 49 simply refers to the case ϕx = x 1/ Because [ L, ] = see Lemma 3, L and are simultaneously diagonalizable Also note that both L and are positive semi-definite operators, so they have non-negative spectrum By expanding A = n 1 j=1 a jl j where L j are given in 48, it is not hard to show that C ϕ = λ L in our situation Therefore, all such Poincaré inequalities have the same prefactor C ϕ for primitive Lindblad equations with GNS-detailed balance 4 Proof of Theorem 18 for the case = Proposition 43 Lower bound of I ρ Assume that the primitive Lindblad equation 1 satisfies GNS-detailed balance Then 5 I ρ λ L 1 e Dρ, ρ D Proof The case = is easier to deal with as certain quantities become more explicit: δd ρ δρ = 1/ ρ 1/ ρ, 1/ ρ 1/ = Γ 1 ρ Γ 1 ρ, Γ 1 ρ ; 1/ 1/ ρ 1/ ρ = log D ρ = log Tr Γ 1 ρ, Γ 1 ρ 1/

22 YU CAO, JIANFENG LU, AND YULONG LU Then, by the definition of the quantum relative -Fisher information 4 Γ 1 ρ I ρ = Γ 1 ρ, Γ 1 ρ, L ρ 1/ = 1 = = 49 Γ 1 ρ, Γ 1 ρ 1/ Γ 1 ρ, Γ 1 ρ 1/ Γ 1 ρ, Γ 1 ρ 1/ Γ 1 ρ, Γ 1 ρ 1/ Γ 1 ρ, Γ 1 L Γ Γ 1 ρ 1/ Γ 1 ρ, LΓ 1 ρ 1/ Γ 1 ρ, LΓ 1 ρ 1/ λ L Γ 1 ρ, Γ 1 ρ 1/ It is not hard to verify that 51 Γ 1 ρ, Γ 1 ρ = 1 + Γ 1 1/ ρ, Γ 1 ρ 1/ It is interesting to note that the term Γ 1 quantum χ -divergence χ 1/ I ρ which yields 5 ρ, Γ 1 ρ turns out to be 1/ ρ, studied in [47] Then Γ 1 ρ, Γ 1 ρ 1 1/ λ L Γ 1 ρ, Γ 1 ρ 1/ = λ L e Dρ e Dρ 1 With Proposition 43, we can immediately show the exponential decay of the sandwiched Rényi divergence with order = Proof of Theorem 18 for the case = From the definition of the quantum relative -Fisher information in 4 and by Proposition 43, d dt D ρ t = I ρ t λ L 1 e Dρt Then d dt log e Dρt 1 λ L After integrating it from time to t, one could find after some straightforward simplification that D ρ t log 1 + e Dρ 1e λ Lt e Dρ 1e λlt Thus D ρ t decays exponentially fast with rate λ L Apparently, prefactor C,ɛ and waiting time τ,ɛ could be chosen as 1 C,ɛ = e D ρ 1, τ,ɛ = D ρ The lower bound of I ρ in Proposition 43 immediately leads into a lower bound of the quantum -log Sobolev constant K in 5,

23 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 3 Corollary 44 [34, Theorem 4] For primitive Lindblad equations with GNSdetailed balance, the quantum -log Sobolev constant K is bounded below by the spectral gap; specifically, 5 K 1 λ min log λ min 1 λ L Proof From Proposition 43, we know that I ρ D ρ λ 1 e Dρ L D ρ 1 λ min λ L logλ min 1, because x 1 e x x is monotonically decreasing on, and D ρ logλ min 1 [34, Lemma 1] The above result follows immediately by taking the minimization over ρ for both sides 43 Proof of the quantum LSI In this subsection, we prove the Proposition 17, which shows the equivalence of primitivity and the validity of quantum LSI, for Lindblad equations with GNS-detailed balance Proof of Proposition 17 First, recall that [9, Theorem 16] has proved that K λ L by linearizing the quantum LSI, thus K > implies λ L > As for the other direction, assume that the Lindblad equation is primitive Suppose ρ t = e tl ρ is the solution of the Lindblad equation and define the mixing time t ɛ := inf {t : D ρ t log1 + ɛ, ρ D} We first recall that see [34, Theorem 5] 53 t ɛ 1 K log logλmin 1 log1 + ɛ Indeed, from 6 for =, D ρ t D ρ e Kt logλ min 1 e Kt By solving logλ min 1 e Kt log1+ɛ, one obtains that t 1 K log logλmin 1 log1+ɛ and the above upper bound for t ɛ follows immediately from its definition Also recall from [34, Theorem 53] that 54 t e 1 κ 1,

24 4 YU CAO, JIANFENG LU, AND YULONG LU Combining the estimates, we have K 31 = κ 1 3 κ t e 1 log K logλmin 1 log1+e 1 1 λ min log λ min 1 log λ L logλmin 1 log1+e 1 > Remark 45 In the proof above, we use the mixing time t e 1 as a bridge to connect κ and K, and use the quantum Stroock-Varopoulos inequality to connect K and K To the best of our knowledge, a direct proof of the quantum LSI 5 has not appeared in the literature For classical systems, Bakry-Émery [3] condition is an important criterion for classical LSI to hold for Fokker-Planck equations It is an interesting open question to see whether there is a quantum analog 44 Proof of Theorem 18 for, via a quantum comparison theorem As a quantum analog of [11, Lemma 36] for Fokker-Planck equations and [43, Theorem 33] for the Ornstein-Uhlenbeck process, in Proposition 46, we will show that under primitive Lindblad equations with GNS-detailed balance 1, the sandwiched Rényi divergence with higher order 1 can be bounded above by the sandwiched Rényi divergence with lower order 1 < < 1, at the expense of some time delay This result helps to prove the exponential decay for sandwiched Rényi divergences under the Lindblad equations, summarized in Proposition 48 Proposition 46 Quantum comparison theorem Let ρ t be the solution of the Lindblad equation with GNS-detailed balance 1 Assume that the quantum LSI 5 holds with constant K Assume also that 55 D ρ ɛ for some fixed ɛ, λ min Then for any two orders 1 < 1, we have 56 D 1 ρ T D ρ, with 57 T = 1 Kη log 1 1, 1 where η = min 1, J min j=1 Λ = λ max λ min exp e ωj 1 Λ, 1 + e ωj Λ λ min ɛ ɛ λ min λ min ɛ

25 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 5 The proof is postponed to Sec 45 Later in Sec 46, we will discuss this quantum comparison theorem in details As a reminder, we have, in fact, implicitly assumed that Lindblad equations under consideration are primitive see Proposition 17 Remark 47 1 The assumption 55 on the initial condition ρ is merely a technical assumption allowing us to obtain a neat expression for the delay time T This assumption 55 is not essential for the inequality 56 to hold; one could remove the assumption at the expense of longer delay time In fact, the validity of the quantum LSI 5 leads to the exponential decay of the quantum relative entropy, so that Dρ t ɛ holds for some t Taking ρ t in place of ρ, one sees that this would imply 56 with a larger decay time T Another reason for imposing the assumption 55 is to avoid technicalities Thanks to the quantum Pinsker s inequality 58 D ρ 1 Tr ρ, the assumption that D ρ ɛ with < ɛ < λmin implies that ρ is a strictly positive density matrix; in fact, we can even show that ρ t D +, for any t This assumption can help avoid many technical issues arising from degenerate density matrices and it validates the usage of results in the previous two sections, eg, properties of operator [X] ω which require X to be strictly positive Proposition 48 Let ρ t be the solution of the Lindblad equation with GNS-detailed balance 1 Assume that the quantum LSI 5 holds with constant K If the sandwiched Rényi divergence decays exponentially fast with rate R for some order > 1, then it decays exponentially fast for any, with rate at least R Proof of Proposition 48 The exponential decay of the sandwiched Rényi divergence means there exists some time t and constant C such that D ρ t C D ρ e Rt, t t Case I: By the monotonicity of sandwiched Rényi divergences with respect to order [38, Theorem 7], D ρ t D ρ t C D ρ e Rt = C D ρ e Rt, for all t t, where C = CD ρ D ρ Case II: > Pick an arbitrary < ɛ < λmin Notice that the quantum LSI implies the exponential decay of the quantum relative entropy, namely, D ρ t D ρ e Kt, one has D ρ t ɛ Moreover, if t T + Hence when t 1 K log Dρ ɛ 1 max t, K log Dρ ɛ with the choice 1 =, then one obtains from Proposition 46 that, where T is the time delay given in the Proposition 46 D ρ t D ρ t T C D ρ e Rt T = C D ρ e Rt,

26 6 YU CAO, JIANFENG LU, AND YULONG LU where C = CD ρ ert D ρ Combining both cases above finishes the proof of Proposition 48 Proof of Theorem 18 Set = and R = λ L in the Proposition 48 Theorem 18 follows immediately from the Proposition 48 and the case = proved in Sec 4; the validity of the quantum LSI has been shown in Proposition 17 The expression of T,ɛ follows immediately from the Proposition Proof of quantum comparison theorem We now turn to the proof of Proposition 46 Consider the density matrix Γ 1 ρ ϱ = ρ Z Z, where normalization constant Z = Tr ρ and ρ D + Then apparently, ϱ D + The log-sobolev inequality 5, with ρ in 5 chosen as ϱ above, becomes Tr ρ log ρ Z logz Tr ρ log 59 1 K Tr L ρ logρ log This inequality can be regarded as a variant of the log-sobolev inequality and it will be used later Let us define the operator G X,ω s : M M by 6 G X,ω sa := e ωs X s AX s + e ω1 s X 1 s AX 1 s The lemma belows collects some properties of G X,ω s which will be useful in the proof of Proposition 46 Lemma 49 Assume that X P + and ω R Then the operator G X,ω s satisfies: i For s 1 < s 1, G X,ωs 1 G X,ωs Moreover, the equality is obtained iff ω = and X = ci for some constant c > ii For s [, 1 ], a bound for the spectrum of G X,ωs is n e ω λ minx λ max X I G X,ωs n As a consequence, ηg X,ω s 1 G X,ω s, 1 + e ω λ maxx λ min X s 1, s [, 1 ] I with η = e ω λminx λ maxx 1 + e ω λmaxx λ minx

27 SANDWICHED RÉNYI DIVERGENCE UNDER LINDBLAD EQUATION 7 Proof of Lemma 49 1 For any positive number b, notice that b s + b 1 s is a non-increasing convex function with respect to s on the interval [, 1 ] Moreover, the range of function b s + b 1 s on this interval is [ b, 1 + b]; b s1 + b 1 s1 = b s + b 1 s for some s 1 < s 1 iff b = 1 Let us write the spectral decomposition of X as X = n k=1 λ k k k where λ k > by assumption Then for any matrix A, n A, G X,ω sa = e ω λ s k + e ω λ 1 s k A k,l λ l λ l k,l=1 Since b s +b 1 s is non-increasing with respect to s, so is G X,ω s To achieve equality, we need e ω λ k λ l = 1 for all k, l Hence ω = and there exists some c = λ k for all 1 k n ie, X = ci From the last equation, by using the range of b s + b 1 s, we immediately have A, G X,ω sa 1 + e ω λ k A k,l n 1 + e ω λ maxx A, A λ l λ min X k,l A, G X,ω sa e ω λ k A k,l n e λ ω λ minx A, A, l λ max X k,l which proves part Proof of Proposition 46 Let us first introduce some useful notations Suppose ρ t is evolving according to the Lindblad equation 1 and = 1 + 1e ηkt is also changing with respect to time for some positive constant η, independent of any specific initial condition ρ The parameter η, 1] is yet to be determined later We are interested in the time interval [, T ] such that t= =, t=t = 1 More specifically, T = 1 Kη log Let us define an important quantity ρ,t by 1 61 β ρ,t := Γ t ρ t We also define Z t := Tr Γ 1 ρ t Tr ρ βt,t and F t := 1 log Z t After taking the derivative of F t with respect to time t and arranging terms, βt Z t F t = Ż t Z t logz t We claim that to prove the inequality 56 it suffices to show that F t is monotonically decreasing, ie, F t on the interval [, T ] In fact, from the definition of F and the quantum Rényi divergence 3, we immediately obtain from F T F that 6 D 1 ρ T 1 1 D ρ D ρ 1 1