Convergence Rates for Quantum Evolution and Entropic Continuity Bounds in Infinite Dimensions

By extending the concept of energy-constrained diamond norms, we obtain continuity bounds on the dynamics of both closed and open quantum systems in infinite dimensions, which are stronger than previously known bounds. We extensively discuss applications of our theory to quantum speed limits, attenuator and amplifier channels, the quantum Boltzmann equation, and quantum Brownian motion. Next, we obtain explicit log-Lipschitz continuity bounds for entropies of infinite-dimensional quantum systems, and classical capacities of infinite-dimensional quantum channels under energy-constraints. These bounds are determined by the high energy spectrum of the underlying Hamiltonian and can be evaluated using Weyl’s law.


Introduction
Infinite-dimensional quantum systems play an important role in quantum theory. The quantum harmonic oscillator, which is the simplest example of such a system, has various physical realizations, e.g. in vibrational modes of molecules, lattice vibrations of crystals, electric and magnetic fields of electromagnetic waves, etc. Even though much of quantum information science focusses on finite-dimensional quantum systems, the relevance of infinite-dimensional (or continuous variable) quantum systems in quantum thermodynamics, quantum computing, and various other quantum technologies, has become increasingly apparent (see e.g. [SL,E06] and references therein).
In this paper we make a detailed analysis of the time evolution of time-independent, infinite-dimensional quantum systems. The dynamics of such a system is described by a quantum dynamical semigroup (QDS) (T t ) t≥0 under the Markovian approximation, which is valid under the assumption of weak coupling between the system and its environment. In the Schrödinger picture, this is a one-parameter family of linear, completely positive, trace-preserving maps (i.e. quantum channels) acting on states of the quantum system. In the Heisenberg picture, the dynamics of observables is given by the adjoint semigroup (T * t ) t≥0 where ∀ t ≥ 0, T * t is a linear, completely positive, unital map on the space of bounded operators acting on the system. 1 There are different notions of continuity of QDSs. The case of uniformly continuous QDSs is the simplest, and is easy to characterize (see Sect. 2.1 for a compendium on semigroup theory). A semigroup is uniformly continuous if and only if the generator is bounded. In this paper, we focus on the analytically richer case of strongly continuous semigroups, which appear naturally when the generator is unbounded.
QDSs are used to describe the dynamics of both closed and open quantum systems. 2 Open quantum systems are of particular importance in quantum information theory since systems which are of relevance in quantum information-processing tasks undergo unavoidable interactions with their environments, and hence are inherently open. In fact, any realistic quantum-mechanical system is influenced by its interactions with its environment, which typically has a large number of degrees of freedom. A prototypical example of such a system is an atom interacting with its surrounding radiation field. In quantum information-processing tasks, interactions between a system and its environment leads to loss of information (encoded in the system) due to processes such as decoherence and dissipation. QDSs are useful in describing these processes. The theory of open quantum systems has also found applications in various other fields including condensed matter theory and quantum optics.
Infinite-dimensional closed quantum systems to which our results apply are e.g. described by time-independent Schrödinger operators H = − + V , which are ubiquitous in the literature. Examples of infinite-dimensional open quantum systems, to which our results apply, include, among others, amplifier and attenuator channels, the Jaynes-Cummings model of quantum optics, quantum Brownian motion, and the quantum Boltzmann equation (which describes how the motion of a single test particle is affected by collisions with an ideal background gas). These will be discussed in detail in Sect. 5.
1.1. Rates of convergence for quantum evolution. Let us focus on the defining property of a strongly continuous semigroup (T t ) t≥0 on a Banach space X , namely, the convergence property for all x ∈ X lim t→0 + T t x = x.
In this paper, we are interested in a refined analysis of this convergence, i.e., our aim is to determine the rate at which T t converges to the identity map I as t goes to zero, and study applications of it.
The rate of convergence lim t→0 + T t − I of a semigroup (T t ) t on a Banach space X is linear in time, if and only if the generator A of the semigroup is a bounded operator. To see this, observe that by the fundamental theorem of calculus and d ds T (s) = T (s)A For general strongly continuous semigroups with unbounded generators, however, one merely knows that lim t→0 T t x − x = 0 by strong continuity, and there is no information on the rate of convergence. If the generator, A, of the semigroup is unbounded, all elements x ∈ X that are also in the domain, D(A), of the generator still satisfy a linear time asymptotics by (1.1). This is because Ax is well-defined for x ∈ D(A), and thus (1.1) holds. However, if the generator A is unbounded, then the bound (1.1) is not uniform for normalized x ∈ D(A), since ||Ax|| can become arbitrarily large.
To obtain more refined information on the rate of convergence, we study spaces that interpolate between the convergence with linear rate t 1 [that holds for elements in the domain D(A) ⊆ X of the generator, by (1.1)] and the convergence without an a priori rate, which we might formally interpret as t 0 , for general elements of the space X . More precisely, we consider interpolation spaces, known as Favard spaces F α = F α ((T t ) t ) in semigroup theory [T78], of elements x ∈ X such that for some C x > 0 T t x − x ≤ C x t α with α ∈ (0, 1], for all t > 0.
In order to study convergence rates and analyze continuity properties of QDSs we need to choose a suitable metric on the set of quantum channels. 3 A natural metric which is frequently used is the one induced by the so-called completely bounded trace norm or diamond norm, denoted as • . However, it has been observed in [W17] that if the underlying Hilbert space H is infinite-dimensional, then the convergence generated by the diamond norm is, in general, too strong to capture the empirical observation that channels whose parameters differ only by a small amount, should be close to each other. Examples of Gaussian channels for which convergence in the diamond norm does hold are, for example, studied in [Wi18].
In this case, a weaker norm, namely the energy-constrained diamond norm, (or ECD norm, in short), introduced independently by Shirokov [Shi18,(2)] and Winter [W17, Definition 3], proves more useful for studying convergence properties of QDSs in the Schrödinger picture (see Example 1). It is denoted as • E , where E characterizes the energy constraint.
In this paper, we introduce a one-parameter family of ECD norms, • S,E 2α ; here S denotes a positive semi-definite operator, E is a scalar taking values above the bottom of the spectrum of S, and α ∈ (0, 1] is a parameter (see Definition 2.3). We refer to these norms as α-ECD norms. They reduce to the usual ECD norm for the choice α = 1/2, when S is chosen to be the Hamiltonian of the system. A version of the α = 1/2-ECD norm, for S being the number operator, was first introduced in the context of bosonic channels by Pirandola et al. [PLOB17,(98)].
To illustrate the power of the α-ECD norms over the standard diamond norm, and even over the usual ECD norm, we discuss the example of the (single mode bosonic quantumlimited) attenuator channel with time-dependent attenuation parameter η(t) := e −t (see Example 5 for details): Example 1 (Attenuator channel). Let N := a * a be the number operator, with a * , a being the standard creation and annihilation operators. Consider the attenuator channel att t , with time-dependent attenuation parameter η(t) := e −t . This one is uniquely defined by its action on coherent states |α = e −|α| 2 /2 ∞ (1. 2) The family ( att t ) is then a QDS. As pointed out in [W17], the diamond norm is too strong in many situations. In fact, for any times t = s it is shown in [W17, Proposition 1] that att t − att s = 2.
Thus, no matter whether t and s are close or far apart, the diamond norm is always equal to 2. The ECD norm serves to overcome this problem, since it follows from [ However, as we will show in Example 5, considering the entire family of α-ECD norms provides further improvement, since it allows us to capture the rate of convergence of the channels as t converges to s: for some constant C α > 0 that is explicitly given in Example 5.

Quantum speed limits.
The bounds which we obtain on the dynamics of closed and open quantum systems, immediately lead to lower bounds on the minimal time needed for a quantum system to evolve from one quantum state to another. Such bounds are known as quantum speed limits. Mandelstam and Tamm [MT91] were the first to derive a bound on the minimal time, t min , needed for a given pure state to evolve to a pure state orthogonal to it. It is given by 4 where E is the variance of the energy of the initial state. From the work of [ML98,LT09] it followed then that the minimal time needed to reach any state of expected energy E, which is orthogonal to the initial state, satisfies Moreover, this bound was shown to be tight. If one includes physical constants and formally studies the semiclassical limit → 0, one discovers that the lower bound in t min vanishes. However, it was shown in [SCMC18] that speed limits also exist in the classical regime. The study of speed limits was generalized in [P93] to the case of initial and target pure states which are not necessarily orthogonal, but are instead separated by arbitrary angles. It has also been generalized to mixed states and open systems with bounded generators. Although the quantum speed limit for closed quantum systems that we obtain from the ECD norm (i.e. for α = 1/2), stated in (3.11), is smaller than (1.3), we obtain better estimates on the quantum speed limit for many states using different α-ECD norm. In particular, the approach pursued in this article allows us to deal with: • open quantum systems with unbounded generators, • states with infinite expected energy, and • systems whose dynamics is generated by an operator which is different from that which penalizes the energy.

Explicit convergence rates for entropies and capacities.
It is well-known that on infinite-dimensional spaces, the von Neumann entropy is discontinuous [We78]. Hence, in order to obtain explicit bounds on the difference of the von Neumann entropies of two states, it is necessary to impose further restrictions on the set of admissible states. In [W15], continuity bounds for the von Neumann entropy of states of infinite-dimensional quantum systems were obtained by imposing an additional energy-constraint condition on the states, and imposing further assumptions on the class of admissible Hamiltonians. The latter are assumed to satisfy the so-called (Gibbs hypothesis). Under the energy-constraint condition and the Gibbs hypothesis it is true that for any energy E above the bottom of the spectrum of the Hamiltonian H , the Gibbs state γ (E) = e −β(E)H /Z H (β(E)) 5 is the maximum entropy state of expected energy E [GS11,p. 196]. Bounds on the difference of von Neumann entropies stated in [W15] are fully explicit up to the occurrence of the entropy of a Gibbs state of the form γ (E/ε), where ε is an upper bound on the trace distance of the two states.
Since entropic continuity bounds are tight in the limit ε ↓ 0, we study (in Sect. 7) the entropy of such a Gibbs state in this limit. Note that for the Gibbs state γ (E/ε), the limit ε ↓ 0 translates into a high energy limit. By employing the so-called Weyl law [I16], which states that certain classes of time-independent Schrödinger operators H = − + V have asymptotically the same high energy spectrum, we show that the asymptotic behaviour of the entropy of the Gibbs state is universal for such classes of operators. This in turn yields fully explicit convergence rates both for the von Neumann entropy and for the conditional entropy (see Proposition [Entropy convergence]).
In finite dimensions, continuity bounds on conditional entropies have found various applications, e.g. in establishing continuity properties of capacities of quantum channels [LS09] and entanglement measures [CW03,YHW08], and in the study of so-called approximately degradable quantum channels [SSRW15]. Analogously, in infinite dimensions, continuity bounds on the conditional entropy for states satisfying an energy constraint [W15], were used by Shirokov [Shi18] to derive continuity bounds for various constrained classical capacities of quantum channels. 6 These bounds were once again given in terms of the entropy of a Gibbs state of the form γ (E/ε). Here, ε denotes the upper bound on the ECD norm distance between the pair of channels considered, and E denotes the energy threshold appearing in the energy constraint. Our result on the high energy asymptotics of Gibbs states yields a refinement of Shirokov's results, by providing the explicit behaviour of these bounds for small ε.
The bounds that we obtain on the dynamics of closed and open quantum systems (see Proposition 3.2 and Theorem 1) also allow us to identify explicit time intervals over which the evolved state is close to the initial state. Since entropic continuity bounds require such a smallness condition for the trace distance between pairs of states, we can then bound the entropy difference between the initial state and the time-evolved state (see Example 12).
We start the rest of the paper with some mathematical preliminaries in Sect. 2. These include a discussion of QDSs, definition and properties of the α-ECD norms, and some basic results from functional analysis that we use. In Sect. 3 we state our main results. These consist of (i) rates of convergence for quantum evolution in both closed and open quantum systems, and (ii) explicit convergence rates for entropies and certain constrained classical capacities of quantum channels. The results concerning (i) are proved in Sects. 4 and 5, while those on (ii) are proved in Sect. 7. In Sect. 6 we discuss some interesting applications of our results, in particular to generalized relative entropies and quantum speed limits. We end the paper with some open problems in Sect. 8. Certain auxiliary results and technical proofs are relegated to the appendices.

Mathematical Preliminaries
Notation In the sequel, all Hilbert spaces H are infinite-dimensional, separable and complex. We denote the space of trace class operators on a Hilbert space H by T 1 (H), that of Hilbert-Schmidt operators by T 2 (H), and the q-th Schatten norm by • q , see e.g. [ We denote the space of bounded linear operators between normed spaces X, Y as B(X, Y ), and as B( If there is a constant C > 0 such that x ≤ C y we use the notation x = O( y ). For closable operators A, B the tensor product A ⊗ B is also closable on D(A) ⊗ D(B) and we denote the closure by A ⊗ B as well. For Banach spaces X, Y one has the projective cross norm on the algebraic tensor product X ⊗ Y The completion of the tensor product space with respect to the projective cross norm is denoted by X ⊗ π Y. In particular, H ⊗ π H is naturally identified with the space of trace class operators on H.
Let A, B be positive operators, Bx . Furthermore, we say B is relatively A-bounded with A-bound a and bound b, if D(A) ⊆ D(B) and for all ϕ ∈ D(A): Bϕ ≤ a Aϕ + b ϕ . Strongly continuous semigroups (T t ) that are defined on Hilbert spaces H, can be extended to act on states We employ a version of Baire's theorem [RS1, Theorem 3.8] in our proofs: Theorem (Baire). Let X = ∅ be a complete metric space and (A n ) n∈N a family of closed sets covering X , then there is k 0 ∈ N for which A k 0 has a non-empty interior.

Quantum dynamical semigroups (QDS).
A quantum dynamical semigroup (QDS) (T t ) t≥0 in the Schrödinger picture is a one-parameter family of bounded linear operators T t : T 1 (H) → T 1 (H) on some Hilbert space H with the property that T 0 = id (where id denotes the identity operator between operator spaces and I the identity acting on the underlying Hilbert space), and T t T s = T t+s for all t, s ≥ 0 (the semigroup property). 7 In addition, they are completely positive (CP) and trace-preserving (TP). The adjoint semigroup is denoted as (T * t ), where for each t ≥ 0, T * t is a bounded linear operator on B(H), which is CP and unital, i.e. T * t (I ) = I for all t ≥ 0. Moreover, T * t is the adjoint of T t with respect to the Hilbert Schmidt inner product. Due to unitality, the QDS (T * t ) is said to be a quantum Markov semigroup (QMS).
For our purposes we consider the following notions of continuity for semigroups (S t ) defined on a Banach space X : • strong continuity if for all x ∈ X : lim t↓0 S t x = x, and • weak * continuity if for all y ∈ X * , where X * is the predual Banach space of X , and x ∈ X the map t → (S t x)(y) is continuous.
Uniformly continuous semigroups describe the quantum dynamics of autonomous systems with bounded generators (see e.g. [EN00, Theorem 3.7]). More precisely, every uniformly continuous semigroup (T t ) is of the form T t = e t A for some bounded linear operator A. Such an operator A is called the generator of the QDS. Strongly continuous semigroups describe the quantum dynamics of closed and open quantum systems with unbounded generators in the Schrödinger picture, and will be the main object of interest in this paper. The QMS in the Heisenberg picture on infinite-dimensional spaces is, in general, only weak * continuous: Denoting this QMS by ( * t ) for an open quantum system, we have that for all y ∈ B(H) * ≡ T 1 (H) and x ∈ B(H) the map t → ( * t x)(y) is continuous. The predual of a weak * continuous semigroup is known to be strongly continuous [EN06, Theorem 1.6].
The generator of a strongly continuous semigroup (T t ) on a Banach space X is the operator A on X such that In this case, d dt T t x = AT t x = T t Ax and by integrating we obtain for all x ∈ D(A) is called a contraction semigroup if T t ≤ 1 for all t ≥ 0 and for any λ > 0 the generator A of such a semigroup satisfies the dissipativity condition For λ > 0 the resolvent of the generator of a contraction semigroup can then be expressed by Its complement is the resolvent set r (A), i.e. the set of λ for which (λI − A) −1 exists as a bounded operator. Let A, B be two operators defined on the same space and λ ∈ r (A) ∩ r (B) then the following resolvent identity holds For any self-adjoint operator S on some Hilbert space H there is, by the spectral theorem, a spectral measure E S mapping Borel sets to orthogonal projections such that the self-adjoint operator S can be decomposed as [RS1, Sect. VII] In particular, this representation allows us to define a functional calculus for S, i.e. we can define operators f (S), by setting for any Borel measurable function f : In particular, if f is bounded, then f (S) is a bounded operator as well.
The dynamics of a closed quantum system is described by strongly continuous oneparameter QDSs 8 according to the following definition: (2.6) Since the self-adjoint time-independent Hamiltonian H fully describes the above QDSs, we will refer to both T S t and T vN t as H -associated QDSs. 8 As mentioned earlier, since a QDS for a closed system consists of unitary operators, it extends to a group.

A generalized family of energy-constrained diamond norms.
Motivated by the ECD norm introduced in [Shi18,W17] we introduce a generalized family of such energyconstrained norms labelled by a parameter α ∈ (0, 1], which coincides with the ECD norm for α = 1/2. We refer to these norms as α-energy-constrained diamond norms, or α-ECD norms in short. The notion of a regularized trace is employed in the definition of these norms.
Definition 2.3 (α-Energy-constrained diamond (α-ECD) norms). Let S be a positive semi-definite operator and E > inf(σ (S)) (where σ (S) denotes the spectrum of S) then we define for quantum channels T , acting between spaces of trace class operators, the α-energy constrained diamond norms induced by S for α ∈ (0, 1] as follows: where ρ H = tr C n ρ. Moreover, any α-ECD norm can be expressed as a standard ECD norm by rescaling both the operator and parameter E as The diamond norm is obtained by setting E = ∞ in the above definition. The maximum distance of the α-ECD norm between two quantum channels is two. Of particular interest to us will be (i) the 1/2-ECD norm • S,E 1 , which reduces to the ECD norm • E considered in [Shi18,W17] when S is chosen to be the underlying Hamiltonian, as well as (ii) the 1-ECD norm • S,E 2 , since they penalize the first and second moments of the operator S, respectively. Although the operator S in the ECD norm is not necessarily an energy observable (i.e. Hamiltonian), we will refer to the condition E 2α ≥ tr(S 2α ρ H ) as an energy-constraint.
We show that by studying the entire family of norms, we obtain a more refined analysis for convergence rates of QDSs. Moreover, we allow the generator of the dynamics of the QDS to be different from the operator penalizing the states in the condition E 2α ≥ tr(S 2α ρ H ). This does not only allow greater flexibility but also enables us to study open quantum systems since the generator of the dynamics of an open quantum system is not self-adjoint in general and therefore also not positive.
By extending the properties for the ECD norm with α = 1/2 stated in [W17, Lemma 4], we conclude that: • The α-ECD norm • S,E 2α defines a norm on the space of hermitian preserving superoperators.
• The α-ECD norm • S,E 2α is increasing in the energy parameter E and satisfies for E ≥ E > inf(σ (S)) • In the limit E → ∞ we recover the actual diamond norm sup E>inf(σ (S)) • The following calculation shows that the topology, for α ≤ β, induced by the β norm is not stronger than the topology induced by α , i.e.

Rates of convergence for quantum evolution.
Our first set of results concerns bounds on the dynamics of both closed and open quantum systems. The following quantities arise in the bounds for α ∈ (0, 1]: (3.1) When α = 1/2, the above two expressions reduce to ζ 1/2 = 2 √ 2 and g 1/2 = 2. Our first Proposition provides a bound on the dynamics of the Schrödinger equation (2.5), both in the time-independent and time-dependent setting: Proposition 3.1 (Closed systems 1). Consider a closed quantum system whose dynamics is governed by an unbounded self-adjoint time-independent Hamiltonian H according to (2.5). Let |ϕ 0 ∈ D(|H | α ) with α ∈ (0, 1]. Then the one-parameter group (T S t ) (c.f. (2.5) of Definition 2.1) satisfies, with g α as in (3.1) and t, s ≥ 0 For the non-autonomous Schrödinger equation where H 0 and V (t) are self-adjoint and The bound (3.2) shows that the dynamics governed by the Schrödinger equation is α-Hölder continuous in time on sets of |ϕ ∈ H with uniformly bounded |H | α |ϕ . The bound is also tight, at least for α = 1, as the prefactor becomes exactly one as α → 1 which is illustrated in Fig. 1. From the bound on the dynamics of the Schrödinger equation in Proposition 3.1, we obtain an analogous result for the dynamics of the von Neumann equation (2.6). The latter result generalizes and improves the bound in [W17, Theorem 6], by providing a bound with rate t 1/2 rather than t 1/3 for the ECD norm, which implies faster convergence to zero [see (3.6) of the following Proposition and (3.5) In particular, when α = 1/2 we find for the ECD norm Moreover, for times |t − s| α ≤ 1/( √ 2g α ), any n ∈ N, and pure states |ϕ ϕ| ∈ D(H ⊗ C n ) satisfying the energy constraint condition tr |ϕ ϕ| (|H | 2α ⊗ I C n ) ≤ E 2α one can slightly ameliorate (3.5) such that In Fig. 2 we see that estimate (3.5) globally improves the estimate stated in [W17, Theorem 6]. For times larger than the time interval [0, 1/4] that is shown in Fig. 2   of two quantum channels and therefore only provide trivial bounds. The bound on the pure states (3.7) however, is especially an improvement over the other two (3.5) for large times.
The above results which are proved in Sect. 4 provide estimates on the dynamics of closed quantum systems. In Sect. 5 we develop perturbative methods to obtain bounds on the evolution of open quantum systems which have the same time-dependence, i.e. α-Hölder continuity in time, as the estimates on the dynamics of closed quantum systems stated in Proposition 3.2.
We focus on open quantum systems governed by a QDS ( t ) with a generator which is unbounded but still has a GKLS-type form. The latter is obtained by a direct extension of Theorem GKLS under some straightforward assumptions, which are discussed in detail in Sect. 5. To state our results on open systems, we define (3.8) In the sequel, we write ω • to denote either one of them.

Theorem 1 (Open systems). Let H be a self-adjoint operator on a Hilbert space H and (L l ) l∈N a family of Lindblad-type operators, generalizing the Lindblad operators L l of
. Then, let α ∈ (0, 1] and let either of the following conditions be satisfied:

Assume that K is relatively H -bounded with H -bound a and bound b. If G := K − i H on D(H ) is the generator of a contraction semigroup, then for energies
of the open system in the Schrödinger picture, generated by L as in (5.6), satisfies, for any c > 0 the α-Hölder continuity estimate For α = 1/2 the above inequality reduces to is the generator of a contraction semigroup, then for energies E > inf(σ (|K |)) the QDS ( t ) of the open system in the Schrödinger picture satisfies, for any c > 0, the α-Hölder continuity estimate In (3.10) For α = 1/2 this expression reduces to • Consider an initial state ρ(0) = ρ 0 to the von Neumann equation (2.6) with E 2α ≥ tr(|H | 2α ρ 0 ). The minimal time for it to evolve to a state ρ(t) which is at a Bures angle relative to ρ(0), satisfies Consider an open quantum system governed by a QDS ( t ) satisfying the conditions of Theorem 1. Let ρ 0 denote an initial state, with purity p start = tr(ρ 2 0 ), for which E 2α ≥ tr(|H | 2α ρ 0 ) (or E 2α ≥ tr(|K | 2α ρ 0 )). Then the minimal time needed for this state to evolve to a state with Bures angle θ , satisfies either for ω H or ω K as in (3.8), where the choice of ω • depends on whether one considers the situation (1) or (2) in Theorem 1, (3.14) Moreover, the minimal time to reach a state with purity p fin satisfies

Explicit convergence rates for entropies and capacities.
Our next set of results comprises explicit convergence rates for entropies of infinite-dimensional quantum states and several classical capacities of infinite-dimensional quantum channels, under energy constraints. See Sect. 7 for definitions, details and proofs. The Hamiltonian arising in the energy constraint is assumed to satisfy the Gibbs hypothesis. Continuity bounds on these entropies and capacities rely essentially on the behaviour of the entropy of the is the partition function, for some positive semi-definite Hamiltonian H ) in the limit E → ∞. This asymptotic behaviour is studied in Theorem 3, and discussed for standard classes of Schrödinger operators in Example 11.
Assumption 1 (Gibbs hypothesis). A self-adjoint operator H satisfies the Gibbs hypothesis, if for all β > 0 the operator e −β H is of trace class such that the partition function The asymptotic behaviour of the entropy of the Gibbs states allows us then to obtain explicit convergence rates for entropies of quantum states and capacities of quantum channels. Consider the following auxiliary functions λλ which depend only on the spectrum of H. We obtain the following explicit convergence rates for the von Neumann entropy S(ρ) of a state ρ, and the conditional entropy t 2 1−εt a function on (0, 1 2ε ], with ε ∈ (0, 1). Proposition (Entropy convergence). Let H be a positive semi-definite operator, with E H := inf(σ (H )) ≥ 0, on a quantum system A satisfying the Gibbs hypothesis and assume that the limit ξ := lim λ→∞ and ε and δ as in (2), the conditional entropy (7.2) satisfies as ε ↓ 0 For the constrained product-state classical capacity C (1) , whose expression is given by (7.17), and the constrained classical capacity C, defined through (7.18), we obtain the following convergence results:

Proposition (Capacity convergence). Consider positive semi-definite operators H A on a Hilbert space H A and H B on a Hilbert space H B , where H B satisfies the Gibbs hypothesis. Moreover, let E H
. We also assume that the limit ξ := Then for t ∈ (0, 1 2ε ] the capacities satisfy +2g(εr ε (t)) + 4h(εt), as ε ↓ 0. (3.18)

Closed Quantum Systems
In this section we study the dynamics of closed quantum systems in α-ECD norms.
From Proposition A.1 in the appendix it follows that if a state ρ = ∞ i=1 λ i |ϕ i ϕ i | satisfies the energy constraint tr(S 2α ρ) < ∞ for some positive operator S, then all |ϕ i for which λ i = 0, are contained in the domain of S α . However, the expectation value tr(Sρ) of an operator S in a state ρ can be infinite even if all the eigenvectors of ρ are in the domain of S. This is shown in the following example.
Proposition 3.2 implies that any group T vN t (ρ) = e −it H ρe it H , with self-adjoint operator H , is continuous with respect to the ECD norm induced by |H | without any further assumptions on H besides self-adjointness. Before proving this result, we start with the definition of the Favard spaces [EN00, Chap. 2., Sect.5.5.10] or [BF15,Sect. 4] and an auxiliary lemma: Definition 4.1 (Favard spaces). Let (T t ) be a contraction semigroup, i.e. for all x ∈ X : T t x ≤ x , on some Banach space X , then for each α ∈ (0, 1] we introduce Favard spaces of the semigroup: In order to link Favard spaces to QDSs, we require a characterization of these spaces in terms of the resolvent of the associated generator. (4.2) We rewrote A = λI + (A − λI ) to get (1) and we used the representation of the resolvent as in (2.3) for (2). Hence, it follows that by taking the norm of (4.2) where we used the definition of the Favard spaces F α in (1), and computed the integral to obtain (2). Conversely, let where now x λ ∈ D(A). Then, using identity (2.1) we get (1) where (2) follows from contractivity of the semigroup, and used the definition of x λ and K to obtain (3) and (4), respectively. For y λ , the triangle inequality and contractivity of the semigroup imply that Combining both estimates (4.3) and (4.4) shows by the triangle inequality Optimizing the right-hand side over λ > 0 proves that x ∈ F α , since the right-hand side is finite, Here, we used the functional calculus, see Sect. 2.2, in (1), optimized over λ to show (2), and used again the functional calculus in (3) which implies the claim.
It is known that if the generator A is defined on a Hilbert space H, then the Favard space F 1 coincides with the operator domain D(A) [EN00, Corollary 5.21]. As all QDSs are contractive, it suffices to establish a bound at t = 0, since by contractivity of the semigroup for t ≥ t 0 ≥ 0 : (4.6) The above lemma then implies Proposition 3.1, which provides a bound on the dynamics of the Schrödinger equation (T S t ) as shown below (Fig. 3). Proof of Proposition 3.1. The result on the autonomous dynamics follows directly by rearranging the estimate and using (4.6) to transfer the result to arbitrary times t, s. The non-autonomous result follows from the variation of constant formula such that by using the result for the autonomous semigroup we obtain where g α is given by (3.1). The general result follows by considering the initial state U s |ϕ 0 at initial time t 0 = s.
Before extending the above result to the dynamics of the von Neumann equation (2.6) for states on the product space H ⊗ C n , we need another auxiliary Lemma on the action of the Schrödinger dynamics on states: Lemma 4.3. The tensor product of the strongly continuous one-parameter group (4.7) Proof. By Proposition A.1 we can assume that all eigenvectors (ϕ i ) of ρ are in D(|H | α ⊗ I C n ) as the right-hand side in (4.7) is infinite otherwise. The generator of (T S t ⊗ id B(C n ) ) acting on trace class operators is the operator −i H ⊗ I C n acting on some set of trace class operators [NS86, Section A-I 3.7]. Using the results from Lemma 4.2 it suffices to bound for λ > 0 accordingly. From the spectral decomposition ρ = ∞ i=1 λ i |ϕ i ϕ i | of a state, the claim then follows immediately from the following bound where we applied Hölder's inequality in (1), used the spectral decomposition of the state and the functional calculus, as in Sect. 2.2, in (2), optimized over λ and applied the functional calculus again in (3), and used in (4) again the spectral decomposition of the state, as well as From estimate (4.7) we can then state the proof of Proposition 3.2: Proof of Proposition 3.2. From a simple application of the triangle inequality and the unitary quantum evolution we conclude that such that by applying Lemma 4.3 in (1) and the energy constraint in (2), we obtain the result for the ECD norms The estimate on pure states follows immediately from Proposition 3.1 after expressing the trace distance in terms of the Hilbert space norm.
The preceding Propositions 3.1 and 3.2 show that the quantum dynamics of closed quantum systems generated by some self-adjoint operator H is always continuous with respect to the α-ECD norm induced by the absolute value of the same operator H .
We now do a perturbation analysis for the convergence in α-ECD norm: such that In particular, if S α is also relatively |H | α -bounded, then the α-ECD norms • S,E 2α and • |H |,E 2α are equivalent in the sense there are constants c 1 , c 2 > 0 such that for all quantum channels T Proof. Consider a density matrix with spectral decomposition ρ = ∞ i=1 λ i |ϕ i ϕ i |. If any of the |ϕ i / ∈ D(S α ) then tr(S 2α ρ) = ∞ as in Proposition A.1. Thus, we may assume that all |ϕ i ∈ D(S α ). Therefore, if tr(S 2α ρ) ≤ E 2α then also tr(|H | 2α ρ) ≤ a tr(S 2α ρ) + b ≤ a E 2α + b which proves the Proposition, since the estimate follows from Proposition 3.2.
The previous result allows us to study QDSs generated by complicated Hamiltonians using more accessible operators penalizing the states in the ECD norms. We illustrate this in the following example where we see that it suffices to penalize the kinetic energy of a state and still obtain convergence for the semigroup of the Coulomb Hamiltonian.
Example 3 (Coulomb potential). If H is relatively S bounded and both H and S are positive, then it follows from [RS2,Theorem X.18] that H is also S form-bounded. This is to say that √ H is also √ S bounded. Iterating this idea, we find that H 2 −n is relatively S 2 −n bounded for all n ∈ N 0 . Let H := − + 1 |x| and S := − on L 2 (R 3 ), then H is relatively S-bounded, see for example [RS2,Theorem X.15]. Thus, the semigroup T vN t (ρ) := e −it H ρe it H is 2 −n -Hölder continuous in time with respect to • S,E 2 1−n . We provide a simple example showing that it is impossible to select arbitrary unbounded self-adjoint operators to penalize the energy in the diamond norm and still have the same convergence rates in time: Example 4 (Harmonic oscillator). Let H osc := − +|x| 2 be the dimensionless Hamiltonian of the harmonic oscillator on D(H osc ) := ϕ ∈ L 2 (R d ); ϕ, |x| 2 ϕ ∈ L 2 (R d ) .
The perturbation result, Proposition 4.4, essentially relies on operator boundedness and provides explicit bounds to compare the two different α-ECD norms induced by the perturbed and unperturbed operator. This result is a special case of a more abstract result, stated as Proposition B.1 in Appendix B, that relies on the special geometry of the space of trace class operators. It yields the same rate t α for the convergence with respect to the perturbed and unperturbed norms. However, it does not provide an explicit prefactor.

Open Quantum Systems
We start with an auxiliary Lemma that provides sufficient conditions under which a perturbation of the generator of a contraction semigroup leaves its Favard spaces invariant: Lemma 5.1 (Perturbation of Favard spaces). Let A 0 and A = A 0 +B be two generators of contraction semigroups on some Banach space X . Furthermore, we fix some α ∈ (0, 1]. Let λ > 0 and B be relatively A 0 -bounded with A 0 -bound a ≥ 0 and bound b ≥ 0. Then, for any k ≥ 0 such that In particular, the Favard space F α of the semigroup generated by A 0 is contained in the Favard space F α of the semigroup generated by A. Proof. Fix c > 0, then for λ ∈ (0, c] it follows that λ α A(λI − A) −1 ≤ 2λ α ≤ 2c α where we used that by (2.2) and the triangle inequality, For λ > c we obtain from the resolvent identity (2.4) and the triangle inequality By relative A 0 -boundedness of B we obtain for the first term on the right-hand side of (5.2) by splitting up For the second term on the right-hand side of (5.2), we can use (5.1) and submultiplicativity to bound Again, using the relative A 0 -boundedness of B we can estimate the last term Thus, since A 0 generates a contraction semigroup, it follows by (2.2) that λ(λI − A 0 ) −1 ≤ 1, and since λ > c such that we finally obtain the claim of the lemma by putting all estimates together and using Lemma 4.2 The most general form of the generator of a uniformly continuous QMS is the socalled GKLS representation, named after Lindblad [Lin76] and Gorini, Kossakowski and Sudarshan [GKS76]. This construction has been generalized by Davies [Da77] to unbounded generators which is discussed below: [Da77]. Let G : D(G) ⊆ H → H be the generator of a contractive strongly continuous semigroup, that we denote by (P t ) t≥0 in the sequel, and consider Lindblad-type operators (L l ) l∈N  where the series converges in trace norm. To keep the notation short, we write X = X ⊗ I C n for operators X on H and X = X ⊗ id B(C n ) for super-operators. Then, by inserting (5.4) into (5.5) it follows that for all S ∈ B(H ⊗ C n ) and x, y ∈ D(G) ⊗ C n x, * t (S)y = x, Sy + x, * t (S)y = P t x, S P t y + ∞ l=1 t 0 L l P t−s x, * s (S) L l P t−s y ds.

Extension of GKLS theorem to unbounded generators
(5.8) By the representation of the QMS in (5.8), bounds on the dynamics of the full, possibly intricate, QDS ( t ) can be found using the simpler semigroup P t as the subsequent Lemma shows: Lemma 5.2. For arbitrary n ∈ N and states ρ ∈ D(H ⊗ C n ) we have Proof. Consider an approximation of ρ ∈ D(H ⊗ C n ) in trace norm by finite-rank operators ρ m := m i=1 λ i |u i u i | with |u i ∈ D(G) ⊗ C n and λ i ≥ 0 such that ρ m − −−− → m→∞ ρ in trace norm. This one exists by a two-step argument. First, we record that we can always approximate density operators by their finite-rank approximations using the spectral decomposition. Thus, it suffices to approximate operators N i=1 λ i |ϕ i ϕ i | with |u i ∈ H ⊗ C n and arbitrary N . Since all norms on finite-dimensional spaces are equivalent, it suffices to observe that since which shows the claim. Then we estimate, using that ρ m : where we expressed the norm in a weak formulation in (1), applied the spectral decomposition of ρ m in (2), and used (5.8) to obtain (3). The two terms in the second-to-last line of (5.9) satisfy, again by the spectral decomposition of ρ m , Here, we used Hölder's inequality and contractivity of the semigroup ( P t ) to get (1) and then used that the trace norm is the same for any operator and its adjoint to conclude (2). For the last term in (5.9) we obtain by contractivity of the QMS where we used (5.3) in (1), that G is the generator of (P t ) to obtain (2), and finally the fundamental theorem of calculus to obtain (3). We can then rewrite this term by decomposing it as follows where we used cyclicity of the trace in (1). To obtain (2) we used Hölder's inequality together with the contractivity of the semigroup P * t and the fact that the trace norm for operators and their adjoints coincide. Estimating (5.9) by (5.10) and (5.12), we can let m tend to infinity and obtain the bound stated in the lemma.
We are now able to prove Theorem 1 which shows that the uniform continuity for the α-ECD norm which we obtained for closed quantum systems in Proposition 3.2 applies to open quantum systems as well: Proof of Theorem 1. We start by proving the first part of the theorem: That G is the generator of a contraction semigroup if a < 1 follows from [EN00, Theorem 2.7].
First, we observe that K ⊗ I C n is still relatively H ⊗ I C n -bounded with the same bound a [Si15, Theorem 7.1.20].
According to Lemmas 4.2 and 5.2 it suffices to obtain bounds on the rate of convergence for the semigroups ( P t ) on density operators ρ ∈ D(H ⊗ C n ) with spectral where we used Hölder's inequality to get (1), the spectral decomposition of ρ in (2), and Lemma 5.1 to get (3). Then, by expanding the expression above (1) and using the Cauchy-Schwarz inequality (2) we find Applying (4.8) yields the desired estimate on the semigroup ( P t ) and Lemma 5.2 the one on ( t ). By (4.6), we then conclude that t − s

Corollary 5.3. For open quantum systems satisfying the assumptions of Theorem 1 the change in purity is bounded for states
Proof. Applying Theorem 1 to the following estimate yields the claim We continue with a discussion of applications of Theorem 1. Let us start by continuing our study of the quantum-limited attenuator and amplifier channels that we started in Example 1: The domain of the generator consists of all such trace class operators for which the time-derivative at zero, of the semigroup ( att t ) in trace norm exists. The space V on which the above identity holds is a core for L att by [C15, Proposition 6.3.6]. The QDS generated by is denoted as ( Hence, by Theorem 1 with H = 0, parameters a = b = 0, and K = N it follows that (5.13) At least for α = 1/2, we can compare the above asymptotics with the explicit bound that was obtained in [N18]: Consider attenuation parameters η = 1, for the initial state, and η = e −t , for the time evolved state, as in [N18]. If we assume for simplicity that the energy E is integer-valued, then the energy-constrained minimum fidelity, that is the infimum of the fidelity over all pure states of expected energy less or equal to E evolved under the attenuator channel with parameters η, η respectively, defined in [N18, (11) . By the Fuchs-van de Graaf inequality as in (6.3) this yields the short-time asymptotics which has the same scaling both in time and energy as the above estimate (5.13). In analogy to (5.13), we find for the amplifier channel, since inf(σ (N )) = 0 and inf(σ Finally, since M = N + I it follows that Mϕ 2 ≤ 2( N ϕ 2 + ϕ 2 ) and thus by Proposition 4.4 Example 6 (Linear quantum Boltzmann equation [A02,HV09]). Since this example describes scattering effects, that depend on the ratio of mass parameters, we exceptionally include physical constants in this example. Consider a particle of mass M whose motion without an environment is described by the self-adjoint Schrödinger operator H 0 = − 2 2M + V. The linear quantum Boltzmann equation describes the motion of the particle in the presence of an additional ideal gas of particles with mass m distributed according to the Maxwell-Boltzmann distribution μ β ( p) = 1 π 3/2 p 3 β e −| p| 2 / p 2 β where p β = √ 2m/β. Here, we discuss for simplicity the linear quantum Boltzmann equation under the Born approximation of scattering theory [HV09]: Let m red = m M/(m + M) be the reduced mass and n gas the density of gas particles. We assume that the scattering potential between the gas particles and the single particle is of short-range and smooth such that is the Schwartz space [RS1]. In the Born approximation the scattering amplitude becomes The presence of the ideal gas leads then to a constant energy shift H per = −2π 2 n gas m red Re( f (0)) in the Hamiltonian H = H 0 + H per and also to an additional dissipative part [HV08]: Let P = −i ∇ x be the momentum operator, then we introduce operators By combining the attenuator channel with the amplifier channel, and using an operator proportional to the number operator N as the Hamiltonian part, we obtain the example of a damped and pumped harmonic oscillator which found, for example, applications in quantum optics, to describe a single mode of radiation in a cavity [A02]: Example 7 (Harmonic oscillator [A02a]). We consider a scaled number operator as the Hamiltonian H = ζ a * a for some ζ > 0 and damping V (ρ) := γ ↓ aρa * and pumping W (ρ) := γ ↑ a * ρa operators and transition rates γ ↓ , γ ↑ ≥ 0. The damping and pumping processes are described by Lindblad operators L ↓ := √ γ ↓ a and L ↑ := √ γ ↑ a * .
The operator K = − 1 2 L * ↓ L ↓ + L * ↑ L ↑ is then dissipative and self-adjoint, such that Theorem 1 applies, and implies that the QDS ( t ) satisfies for any E > 0 Next, we study the evolution of quantum particles under Brownian motion which is obtained as the diffusive limit of the quantum Boltzmann equation that we discussed in Example 6 [HV09, Sect. 5].
Example 8 (Quantum Brownian motion [AS04,V04]). Consider the Hamiltonian of a harmonic oscillator H = − d 2 dx 2 + x 2 and Lindblad operators for j ∈ {1, 2} given by In particular, choosing γ j = β j turns L j into the annihilation operator L j = γ j d dx + x and L * into the creation operator L * j = γ j − d dx + x which have been considered in the previous example. The Lindblad equation for quantum Brownian motion reads and λ = Im γ * 1 β 1 + γ * 2 β 2 . The operator K = − 1 2 2 j=1 L * j L j is then relatively H -bounded and G = i H − K is the generator of a contraction semigroup on D(H ). By Theorem 1, the QDS ( t ) of quantum Brownian motion satisfies for E > inf(σ (H )) and α ∈ (0, 1] The field of quantum optics is a rich source of open quantum systems to which the convergence Theorem 1 applies and we discuss a few of them in the following example: Example 9 (Quantum optics/Jaynes-Cummings model [CGQ03]). Systems that consist of a harmonic oscillator coupled to two-level systems are among the common illustrative examples considered in quantum optics and within this theory are called Jaynes-Cummings models. A particular example of a Jaynes-Cummings model is a two-level ion coupled to a harmonic trap of strength ν > 0 located at the node of a standing light wave. For a detuning parameter and Rabi frequency , a Master equation with Hamiltonian where η is the Lamb-Dicke parameter, and with Lindblad operators L = √ σ − , L * = √ σ + has been proposed in [CBPZ92] for this model. Here, is the decay rate of the excited state of the two-state ion. The Hilbert space is therefore 2 (N) ⊗ C 2 and as the Lindblad operators are just bounded operators, all conditions of Theorem 1 are trivially satisfied. Thus, it follows that for E > 0 the QDS ( t ) satisfies More generally, plenty of models in quantum optics are special cases of the following form [CGQ03]: acting on a Hilbert space H = 2 (N) ⊗N ⊗ C M . The Lindblad operators are of the form where a k is the annihilation operator of the k-th factor of the tensor product 2 (N) ⊗N and λ k ≥ 0 a positive semi-definite matrix acting on C M . Hence, the operators − 1 2 L * k L k are self-adjoint and dissipative and for a large class of Hamiltonians H the asymptotics of Theorem 1 can be applied.

Generalized Relative Entropies and Quantum Speed Limits
We start with some immediate consequences of Propositions 3.1, 3.2, and Theorem 1 on certain generalized relative entropies and distance measures which are dominated by the trace norm: Of particular interest to us are the α = 1/2-divergences: The α = 1/2-Tsallis divergence is, up to a prefactor, the square of the Hellinger distance [RSI] and satisfies The Consider also the fidelity of two states ρ, σ that we denote by It is related to the trace distance via the Fuchs-van de Graaf inequalities [FG99]: Bures angle θ and Bures distance d B are respectively defined as θ(ρ, σ ) := arccos (F(ρ, σ )) and d B (ρ, σ ) := 2 (1 − F(ρ, σ )).
As a Corollary of Proposition 3.2 for closed quantum systems and Theorem 1 for open quantum systems, we obtain: Corollary 6.2. For closed quantum systems and states ρ ∈ D(H) such that tr ρ |H | 2α ≤ E 2α it follows with the notation introduced in Definition 2.1 that: • The Bures distance and Bures angle satisfy where (a) + := max{a, 0}. For open quantum systems satisfying the conditions of Theorem 1 and states ρ satisfying tr ρ |H | 2α ≤ E 2α (or tr ρ |K | 2α ≤ E 2α ) we obtain for ω • as in (3.8) • For the 1/2-divergences it follows that • For the Bures distance and Bures angle, we obtain implies that D Tsallis 1/2 (ρ||σ ) ≤ ρ − σ 1 .
The study of quantum speed limits, see also the review article [DC17], is concerned with the minimal time for the system needed to evolve from one state of expected energy E to another state that is a certain distance away from the initial state. It has been shown in [ML98,LT09] that the minimal time of a closed quantum system to evolve from an initial state |ϕ 0 to another state that is orthogonal to it, under the evolution given by the Schrödinger equation (2.5), with positive semi-definite Hamiltonian H , satisfies and showed that this bound can be saturated. For arbitrary (Bures) angles θ ∈ [0, π/2] this bound was (partially numerically) extended by Giovannetti, Lloyd, and Maccone [GLM03,GLM03a,GLM04] to While the quantum speed limits for closed quantum system still yield non-trivial statements for dynamics generated by unbounded operators, non-trivial estimates for open quantum systems with unbounded operators do not seem to exist. Let us begin by mentioning some results that hold for open quantum systems with bounded generators. In [CEPH13,UK16] a bound on the purity has been stated saying that to reach a purity p fin := tr(ρ(t) 2 ) from a purity p start := tr(ρ(0) 2 ) the minimal time needed is bounded from below by where L k are the Lindblad operators, • 2 the Hilbert-Schmidt norm, and L is the generator of the associated QDS. Furthermore, a bound on the quantum speed limit in terms of the operator norm of the generator has been derived in [DL13]. In the following remark we see that all these bounds have a pathological behaviour for certain infinitedimensional systems and cannot be sharp in general: Remark 1. Consider a closed system with Hamiltonian S = − d 2 dx 2 on R. The state ψ ∈ L 2 (R) with Fourier transform F(ψ)(x) = c (1+x 2 ) 1/2 where c > 0 is such that ψ is of unit norm. Then, Sψ, ψ = ∞ whereas tr(S α ρ) < ∞ for α < 1/4. Thus, the above bounds (6.5) and (6.6) reduce to the trivial bound t min ≥ 0.
For infinite-dimensional open quantum systems, the first term in the bound on the purity (6.7) reduces to zero if the Lindblad operators are not Hilbert-Schmidt, which is the case for all examples presented in Sect. 5. In particular, if the Lindblad operators are unbounded, then the bound simplifies to t min ≥ 0.
We can now state the proof of Theorem 2: Proof of Theorem 2. The first estimate on the minimal time of the Schrödinger dynamics, follows from the polarization identity of the Hilbert space inner product and Proposition 3.1, which after rearranging yields the claim. For the estimates on the Bures angle we rearrange the estimates in Corollary 6.2, and for the estimate on the purity we rearrange the estimate in Corollary 5.3.

Entropy and Capacity Bounds
In this section, we obtain explicit continuity bounds for different families of entropies of quantum states, and various constrained classical capacities of quantum channels in infinite dimensions. The capacity of a channel is the maximal rate at which information can be transmitted through it reliably. Unlike a classical channel, a quantum channel has various different capacities. These depend, for example, on the nature of the information transmitted (classical or quantum), the nature of the input states (product or entangled), the nature of the allowed measurements at the output of the channel (individual or collective), the availability of any auxiliary resource (e.g. prior shared entanglement between the sender and the receiver), the presence or absence of feedback from the receiver to the sender, etc. From a practical point of view, it is natural, and in fact necessary, to impose energy constraints on the input states to the quantum channels in computing their capacities, since one cannot physically prepare states of infinite energy. Due to the energy constraint, the resulting capacity is called the constrained capacity of the channel. Here we consider three different constrained capacities for transmission of classical information through an infinite-dimensional quantum channel: (i) the constrained product-state capacity, which is the capacity evaluated under the additional constraint that the inputs are product states, (ii) the constrained classical capacity, for which the only constraint is the energy constraint, and (iii) the constrained entanglement-assisted classical capacity, which corresponds to the case in which the sender and the receiver have prior shared entanglement. 12 If : denotes an infinite-dimensional quantum channel, then the energy constraint on an input state ρ to the channel is given by tr(H A ρ) ≤ E, where H A is the Hamiltonian of the input system A. 13 For n identical copies of the channel, the energy constraint is tr(H A n ρ A n ) ≤ n E, where ρ A n ∈ D(H ⊗n A ) and The capacities are evaluated in the asymptotic limit (n → ∞). For their operational definitions see [H03]. Obviously these capacities depend not only on the channel, , but also on H A and E. We denote the three different classical capacities introduced above as follows: (i) C (1) ( , H A , E), (ii) C( , H A , E), and (iii) C ea ( , H A , E). Expressions for these capacities have been evaluated [H03] and are given by equations (7.17), (7.18) and (7.15), respectively.
Besides classical capacities, we also study convergence of entropies in this section. It has been shown by Lindblad in [Lin74,Lin73] that the quantum relative entropy for states ρ, σ ∈ D(H) is well-defined for a complete orthonormal basis of eigenvectors ( ϕ n ) n of either ρ or σ as If the underlying Hilbert space is infinite-dimensional, the von Neumann entropy depends discontinuously on the states and is even unbounded in every neighbourhood: More precisely, let ε > 0, then in the ε-neighbourhood (in trace distance) of any state ρ, there is another state ρ (say) for which S(ρ ) = ∞ [We78]. In general the von Neumann entropy is only lower semicontinuous i.e. given a state ρ, if (ρ n ) n∈N denotes a sequence of states such that ρ n − ρ 1 − −−→ n→∞ 0, then S(ρ) ≤ lim inf n S(ρ n ) [We78]. Although, this explains why there are no continuity bounds for the entropy of states in infinite dimensions, the following observation shows that under additional assumptions, such continuity estimates can indeed be derived: Let H be a self-adjoint operator such that a Gibbs state γ (β) := e −β H tr(e −β H ) ∈ D(H) is well-defined for all β > 0, 14 the sequence of states (ρ n ) converge in trace norm ρ n − ρ 1 → 0, and the energies tr(ρ n H ), tr(ρ H ) are uniformly bounded, then the entropies converge S(ρ) = lim n→∞ S(ρ n ) as well [We78]. Thus, continuity bounds on the von Neumann entropy can be expected to hold for energy-constrained states when the underlying Hamiltonian defines a Gibbs state for all inverse temperatures. Indeed, in [W15] for entropies and [Shi18] for capacities, such continuity estimates have been established which are fully explicit up to the asymptotic behaviour of the Gibbs state for high energies. It is precisely this asymptotic behaviour that we discuss in this section.
We now want to compare the delicate continuity properties of the von Neumann entropy with the properties of the Tsallis-(T q ) and Rényi-(S q ) entropies: Definition 7.1. The q-Tsallis entropy is for q > 1, using the q-Schatten norm, defined by The q-Rényi entropy is for q > 1, using the q-Schatten norm, defined by Unlike the von Neumann entropy, our next Proposition shows that the Tsallis and Rényi entropies are Lipschitz continuous, without any assumptions on the expected energy of the state or the Hamiltonian: Proposition 7.2. Let ρ, σ ∈ D(H) be two states and α ∈ (0, 1]. Then the q-Tsallis entropy satisfies the global Lipschitz estimates Assume now that there is additionally some δ > 0 such that ρ q ≥ δ > 0 and ρ − σ q ≤ ε < δ. Then the q-Rényi entropy satisfies the local Lipschitz condition In particular, under the assumptions of Theorem 1, it follows that for states ρ with tr |H | 2α ρ ≤ E 2α (or tr |K | 2α ρ ≤ E 2α ) we obtain for any c > 0 and t, s > 0 for the QDS ( t ) of an open quantum system with ω • as in (3.8) If the initial state satisfies additionally ρ q ≥ δ > 0 then up to sufficiently short times Proof. The statement on the Tsallis entropy follows directly from where we used the mean-value theorem for the function f (ρ) = ρ q q on states for which ρ q ≤ 1 and the inverse triangle inequality in (1), and ρ q ≤ ρ 1 in (2).
The additional assumptions for the Rényi entropy imply that σ q ≥ δ− ρ − σ q ≥ δ − ε > 0 which we need for the local Lipschitz condition on the logarithm. Proceeding as for the Tsallis entropy this shows .
For any finite energy, the term (1) is a finite sum, while (2) is an infinite sum (since the operator H is unbounded). Thus, if the energy would remain finite, as β(E) ↓ 0, then (1) is finite whereas (2) becomes infinite. Conversely, if the temperature would remain finite (β > 0) as E → ∞, then (1) tends to negative infinity while (2) vanishes by the dominated convergence theorem.
A straightforward calculation shows that the entropy of the Gibbs state satisfies [W17, p. 7] S(γ (β(E))) = log tr e −β(E)H + β(E)E. (7.4) In the proof of [Shi06, Proposition 1] it is shown that lim ε↓0 εS(γ (β(E/ε))) = 0. In the following, we want to derive precise asymptotics of (7.4) in the high energy limit and discuss applications of it. Before entering the general theory, let us study the fully explicit case of the harmonic oscillator first: We stress that this shows that for the special case when the Hamiltonian is the harmonic oscillator, then S(γ (β(E))) behaves like log(E) as E → ∞.
Our aim in this section is to show that, in some sense, the logarithmic divergence of the entropy of the Gibbs state, as E → ∞, is not a special feature of the harmonic oscillator but universal for many classes of Hamiltonians. This result allows us then to state explicitly a rate of convergence in continuity bounds on entropies and capacities. We start with some preliminary related ideas: Let H be a self-adjoint operator with compact resolvent on L 2 (X, dν(x)). The spectral function e H of H , is defined as [Hö07,(17.5 where ϕ j are the eigenfunctions corresponding to the eigenvalue λ j of the operator H . The number of eigenvalues of H that are at most of energy E, counted with multiplicities, is then given by The famous Weyl law [I16] gives an asymptotic description of N H for certain classes of operators in the limit of high energies, and shows that this distribution is universal. In many cases, even the precise asymptotics of eigenvalues is known. We will show that to estimate the entropy of the Gibbs state at high energies for arbitrary Hamiltonians, it suffices to estimate the ratio of the following two auxiliary functions for high energies λλ . (7.5) We also observe that the simple estimate 2λλ ≤ λ 2 + λ 2 implies that N where Weyl's law ensures that these two functions have a universal asymptotic behaviour as E → ∞ for large classes of operators. The next theorem shows that the high energy asymptotics for the entropy of the Gibbs state is uniquely determined by the high energy spectrum of the Hamiltonian expressed in terms of functions defined in (7.5). is well-defined. Let the inverse temperature β(E) be given as the solution of (7.3). For high energies, the inverse temperature satisfies the asymptotic law In the same high energy limit the partition function satisfies Finally, the entropy of the Gibbs state satisfies Proof. The derivative of the inverse temperature as a function of the inverse energy satisfies where we used (7.3) in the second equality. We obtain then for the two-sided Laplace transform of the auxiliary function N where we used the definition of the two-sided Laplace transform in (1), Fubini's theorem to get (2), and by computing the integral we obtained (3). By an analogous calculation, we find that for G(β) := λ,λ ∈σ (H ) λλ e −β(λ+λ ) , (7.10) The quotient of (7.9) and (7.10) allows us to recover the factor appearing in (7.8) (7.11) From the existence of the limit ξ = lim λ→∞ in the assumption of the theorem, we conclude that for any ε > 0 there is λ 0 > 0 large enough such that for all λ ∈ R Hence, by applying the two-sided Laplace transform to this inequality we infer that for Example 11. The entropy of the Gibbs state for the quantum harmonic oscillator as in Example 10 satisfies S(γ (E)) = log(E)(1 + o(1)), as E → ∞.
The entropy of the Gibbs state for regular Sturm-Liouville operators defined through The entropy of the Gibbs state for multi-dimensional second order differential operators [Hö07,Sect. 17.5] on bounded open subsets of R n with smooth boundary, Dirichlet boundary condition, and positive semi-definite matrix (g jk ) on such that H is self-adjoint on L 2 (X, dμ(x)) satisfies Calculation. Instead of just referring to Example 10 for the harmonic oscillator, we apply Theorem 3: Harmonic oscillator: By applying the Cauchy product formula, we find from the Harmonic oscillator spectrum {n + 1/2; n ∈ N 0 } N ↑ H (n + 1) = n k=0 k i=0 i + 1 2 2 = (2n 2 + 6n + 3)(2 + n)(1 + n) 24 and By Theorem 3 it follows that S(γ (E)) = log(E)(1 + o(1)), as E → ∞.   1 + o(1)), as E → ∞. (7.14) In Fig. 4a we compare the true inverse temperature β(E) of the Gibbs state for the quantum harmonic oscillator as in Example 10 with the asymptotic law β(E) ≈ 1 E obtained from Theorem 3. In Fig. 4b we compare the inverse temperature of the Gibbs state for the Hamiltonian describing a particle in a box of length 1 √ 8 with the asymptotic law β(E) ≈ 1 2E we obtained in Example 11. The following Proposition, which relies on Theorem 3, shows that for large generic classes of Schrödinger operators with compact resolvent, the entropy of the Gibbs states obeys a universal high energy asymptotic behaviour.
The Proposition [Entropy convergence] then follows as an application of Theorem 3, which provides an explicit rate of convergence for entropies on infinite-dimensional Hilbert spaces: Proof. (Prop. ent. conv.) Under the assumptions stated in the Proposition and if E H = 0, Lemmas 15 and 16 in [W17] show that the von Neumann entropy satisfies The conditional entropy satisfies by [W17, Lemma 17] However, to apply [W17], H has to satisfy the assumption E H = 0.
This implies that .
Combining these results with Theorem 3 yields the claim of the Proposition.
Another correlation measure for a bipartite state ρ AB ∈ D(H A ⊗H B ) is the quantum mutual information (QMI) and is defined in terms of the relative entropy (7.1). Let : T 1 (H A ) → T 1 (H B ) be a quantum channel and H A a positive semi-definite operator on H A . The energyconstrained entanglement assisted capacity, C ea ( , H A , E) was proved to be given by the following expression (7.15) by Holevo and Shirokov [HS13], where ρ is a pure state in D( The two Corollaries C.1 and C.2 of Theorem 3 that are stated in Appendix C provide convergence rates on QMI and hence on C ea .
We continue our discussion of attenuator and amplifier channels, that were defined in Example 5 by studying their convergence of entropies.
Example 12 (Entropy bounds for attenuator and amplifier channels). We start by discussing how the expected energy of output states of these channels with time-dependent attenuation and amplification parameters behave as a function of time.
Hence, it follows that tr(Nρ att (t)) = tr(Nρ att 0 )e −t and tr(Mρ amp (t)) = tr(Mρ amp 0 )e t . Let ε > 0 and t 0 be sufficiently small such that t 0 ≤ 1 Remark 3. If the ensemble is of the form μ = ∞ i=1 p i δ ρ i for probabilities p i ≥ 0 summing up to one ∞ i=1 p i = 1 and delta distributions associated with states ρ i ∈ D(H) then the ensemble is also called discrete. In this case the barycenter state is just Let be a quantum channel, then the pushforward ensemble of such a discrete ensemble becomes just * (μ) = ∞ i=1 p i δ (ρ i ) . Discrete ensembles play a particularly important role in the study of capacities. Let D E be the set of discrete ensembles with barycenter state ρ of energy less than E under a positive semi-definite Hamiltonian. Let be a channel, H a positive semidefinite Hamiltonian, and μ a discrete ensemble. The constrained product-state classical capacity is known to be given by the Holevo capacity χ * ( ) := sup μ∈D E χ( * (μ)), defined in terms of the Holevo quantity, by C (1) ( , H, E) = χ * ( ). (7.17) The full classical capacity is given in terms of C (1) as follows

Open Problems
Concerning the first part of the paper, it would be desirable to study extensions of our work to non-autonomous systems, such as systems described by a Schrödinger operator with time-dependent potentials. For the Schrödinger equation, an application of the variation of constants formula yields a bound for such systems as well (see Proposition 3.1). This should also work, under suitable assumptions, for non-autonomous open quantum systems. However, more mathematical care may be needed for the latter.
To answer the important questions: (i) "How fast can entropy increase?"-for any infinite-dimensional open quantum system whose dynamics is governed by a QDS, and (ii) "How fast can information be transmitted?"-through any quantum channel (obtained by freezing the time parameter in the QDS), it seems necessary to find bounds on the evolution of the expected energy for the state of the underlying open quantum system over time (as has been done for the case of the attenuator and amplifier channels in Example 12. 15 ) To our knowledge, such bounds have not been obtained in full generality yet. See also [BN88,DKSW18,OCA] for related results on question (i).
The first step to answer these two questions was provided by Winter [W15] and Shirokov [Shi18], who derived continuity bounds on entropies and capacities, respectively. Our paper provides, as a second step, a time-dependent bound on the evolution of the expected energy of the state of the open quantum system, which enters these continuity bounds through the energy constraint. Understanding the behaviour of this expected energy as a function of time is needed in order to infer, from the continuity bounds, how fast entropies and capacities can change.
It would be furthermore desirable to extend Theorem 3 to higher-order terms. In Fig. 4a we see that the leading-order approximation for the inverse temperature provided by Theorem 3 is almost indistinguishable from the true solution for the harmonic oscillator whereas the leading-order approximation in Fig. 4b for the particle in a box seems to converge somewhat slower than the true solution. A better understanding of higher order terms should be able to capture these behaviours more precisely. Note that by considering SE S [0,n] instead of just S we work with a bounded operator that commutes with the series and also the application of SE S [0,n] to elements ϕ i is welldefined.

Appendix B. Dynamics of QDS in ECD Norms
The following Proposition is an adaptation of the uniform boundedness principle to the α-ECD norm and can be applied as a perturbation theorem for convergence in α-ECD norms.