Fundamental Speed Limits to the Generation of Quantumness.

Quantum physics dictates fundamental speed limits during time evolution. We present a quantum speed limit governing the generation of nonclassicality and the mutual incompatibility of two states connected by time evolution. This result is used to characterize the timescale required to generate a given amount of quantumness under an arbitrary physical process. The bound is found to be tight under pure dephasing dynamics. More generally, our analysis reveals the dependence on the initial and final states and non-Markovian effects.

Quantum speed limits (QSLs) provide a upper bound to the rate at which a physical system can evolve. Due to their fundamental nature, QSLs have found applications in a wide variety of fields including quantum information processing123, quantum metrology45, quantum simulation6, quantum thermodynamics7, quantum critical dynamics89, quantum control10111213 and other quantum technologies.

The first rigorous QSL was derived as a time-energy uncertainty relation providing a lower bound to the required passage time τ for a system to evolve from an initial state |ψ0〉 to a final state |ψ〉 = U(τ, 0)|ψ0〉, where U(τ, 0) is the time-evolution operator associated with the driving Hamiltonian H. It was shown that τ ≥ arccos (|〈ψ0|ψ〉|)/ΔE, where ΔE is the energy dispersion of the initial state141516171819. The modern formulation of QSL for unitary processes takes into account an alternative expression as an upper bound for the speed of evolution, the mean energy of the system, that can replace the role of energy dispersion ΔE1232021. A geometric interpretation provides an intuitive understanding of the QSL bound as a brachistochrone22 where the geodesic2324 set by the Fubini-Study metric in (projective) Hilbert space is travelled at the maximum speed of evolution achievable under a given Hamiltonian dynamics25262728. Time-optimal evolutions are often explored in the context of quantum control theory, where the existence of a QSL has been shown to limit the performance of algorithms aimed at identifying optimal driving protocols1112. More recently, QSLs have been extended to open quantum dynamics where the system of interest is embedded in an environment2930313233. The evolution need not be restricted to a master equation and can be alternatively described by general quantum channels2930. These new QSLs to non-unitary evolution have been formulated in terms of a variety of norms of the generator of the dynamics. Similar bounds can be expected to apply to classical processes as well34.

While for certain applications it might suffice to characterize QSL exclusively through the properties of the generator of the dynamics435, a reference to the initial and time-evolving states generally becomes unavoidable. This is particularly the case for externally driven systems or open quantum systems exhibiting non-Markovian effects resulting from the finite-memory of the environment. We further notice that when the dynamics of a system is registered by monitoring a given observable, the standard QSL governing the fidelity decay can become too conservative, and even fail to capture the right scaling of the time scale of interest with the system parameters. A prominent example is provided by thermalization, where the identification of the relevant time scale remains an open problem36.

Identifying the minimal passage time for arbitrary physical processes is as well crucial to understand the quantum-to-classical transition37. This transition is of particular relevance in composite quantum systems exhibiting non-classical correlations, with applications to a variety of fields38. To characterize the crossover between the quantum and classical worlds of a single physical system, one can define the notion of quantumness on the non-commutativity of the algebra of observables39 in a way that it is experimentally measurable40. In this framework a system is found to be classical if all its accessible states commute with each other.

In this work, we exploit the definition of quantumness involving the non-commutativity of the initial and final states of the system of interest. We derive lower bound for the timescale required to generate a given amount of quantumness, that quantifies the degree to which the time-evolving state is mutually incompatible with the initial state under arbitrary dynamics. The new bound allows one to classify different dynamics according to their power to generate nonclassicality and is shown to be saturated under pure dephasing dynamics, whether it is induced by a Markovian or a non-Markovian environment.

Quantum speed limit to the dynamics of quantumness

The nonclassicality of quantum systems can be conveniently quantified using the Hilbert-Schmidt norm of the commutator of two states, which is proposed to witness the “state incompatibility” between any two admissible states ρ and ρ3940. The “quantumness” is then defined as

where the pre-factor is required for normalization and is the Hilbert-Schmidt norm of A. As a quantumness witness,

and Q(ρ, ρ) = 0 iff [ρ, ρ] = 03940. Choosing ρ = ρ0 and ρ = ρ, Q (ρ0, ρ) allows one to quantify the capacity of an arbitrary physical process to generate or sustain quantumness in case of [ρ0, ρ] ≠ 0. Clearly, if ρ0 is a diagonal density matrix in a given basis and time evolution just alters the weight distribution without generating coherences, the quantumness between the initial and time-evolving states Q(ρ0, ρ) remains zero. Consequently we can generally expect a QSL different in nature from those previously derived for the fidelity decay, which would remain valid as weaker lower bounds. As shown in the section Methods, we obtain the lower-bound of evolution time associated with quantumness,

where the time-average is denoted by , and is the Liouville super-operator describing the time derivative of density matrix: .

Results

The lower bound obtained in Eq. (3) constitutes the main result of this work. In what follows, this bound is analyzed in a series of relevant scenarios, that will be used to identify the salient physical principles governing the generation of quantumness. After a discussion of its dynamics in isolated systems we consider a system embedded in an environment, exhibiting possible non-Markovian effects, and discuss the limits of pure dephasing and dissipative processes.

Unitary quantum dynamics

Consider a general two-parameter unitary transformation for a two-level system (setting ħ ≡ 1 from now on)

where θ and α are arbitrary real functions of time and σ is the Pauli operator41 (We apply the standard conventions that σ = |1〉〈0| + |0〉〈1| and σ = i|1〉〈0| − i|0〉〈1|). When the system is prepared in an initially pure states ρ0 = |0〉〈0|, it evolves into ρ = |ψ〉〈ψ| with . In this case, the system Hamiltonian is found to be

For the creation of quantumness, we require [ρ0, ρ] ≠ 0, i.e., sin 2θ ≠ 0. It follows from Eq. (3), that

where

Specially in the case that both and α are constant numbers, . According to Eq. (5), the dynamics is generated by . Then due to Eq. (6), the exact evolution time saturates the lower bound τ = τ in the regime that 0 < θ < π/4. It is worth pointing out that in this case, the independence of the result on the angle α is coincident with the result in the Bloch-vector formalism set up by the angles θ and α. When the final state reaches θ = π/4, the lower bound τ attains its maximal value and starts to decline when θ goes over this point. After that, Eq. (6) remains valid while loosing the tightness in the regime 0 < θ < π/4. In another special case with both θ and constant during the evolution, one can find that by Eq. (5). Then for a target state characterized by a nonvanishing α, τ(θ = π/4) = τ/|α| by Eq. (6). Therefore, the QSL ruling the evolution of quantumness exhibits a pronounced dependence on the initial and final states.

A similar analysis can be extended to higher-dimensional systems. Consider the stimulated Raman adiabatic passage (STIRAP)42 in a three-level atomic system, under the Hamiltonian as

The system can be transferred from ρ0 to ρ = |ψ〉〈ψ|, where without disturbing the quasistable state |1〉. The QSL bound also becomes tight and matches the exact time of evolution τ = τ when θ < π/4 and and α are time-independent.

Nonunitary process

In this part, we consider the scenario of open quantum systems. We will use the quantum-state-diffusion (QSD) equation4344 as a general framework to derive the exact master equation before discussing the relevant QSL. In doing so, we treat both Markovian and non-Markovian environments in a unified way. In particular, we consider an Ornstein-Uhlenbeck process for the environmental noise. The correlation function reads

where Γ implies the system-environment coupling strength and γ is inversely proportional to the memory time of the environment. Here, 0 < γ < ∞ and the lower and upper limits of γ correspond to the strongly non-Markovian and Markov environments, respectively. For a single two-level system, the system-environment Hamiltonian is

where ω and ω are the frequencies of the system and the λ-th environmental mode, respectively, and g is their coupling strength. L is the coupling operator and a () is the annihilation (creation) operator for the λ-th environmental mode.

When L = σ, QSD equation describes a pure dephasing process. In the rotating frame with respect to the system bare Hamiltonian, the exact super-operator is found to be

where . If ρ0 = |ψ0〉〈ψ0| where , then 〈0|ρ|0〉 = 〈0|ρ0|0〉, 〈1|ρ|1〉 = 〈1|ρ0|1〉, and , where . In the Markov limit, and then . After ρ0, ρ and are substituted into Eq. (3), it is found that

where β ≡ 2Γτ. Remarkably, the bound is tight and reachable under pure-dephasing dynamics, when τ = τ as shown in Fig. 1. Equation (12) also applies to the non-Markovian case as long as β in Eq. (11) is modified into . Qualitatively, QSL timescale depends on the choice of initial state parameter θ, specifically, the initial population distribution determined by . QSL is therefore symmetric as a function of θ with respect to θ = π/4.

Figure 1

The quantum speed limit timescales τ (based on quantumness) and τ (based on the fidelity) as a function of quantumness Q in the Markovian pure-dephasing processes with different initial states: , where.

Under pure dephasing the bound τ is shown to be identical to the exact time τ in which quantumness is generated.

In Fig. 1, for different initial states, we compare the new QSL timescale τ obtained in Eq. (3) and that τ based on the fidelity evolving with time (see ref. 30) in the presence of a Markovian dephasing environment. Specifically, it was then shown that for an initially pure state, the minimum time for the (squared) fidelity or relative purity F(t) = Tr[ρ0, ρ] to decay to a given value F(τ) is lower bounded by whenever the dynamics is governed by a master equation of the form . Figure 1 illustrates that τ not only provides a tighter bound than τ for the generation of quantumness, but also it actually captures the real evolution pattern. Furthermore, it is known that as the system progressively goes to a steady state, which depends on the initial coherence between the up and down states, the dephasing rate should be asymptotically slowed down. This pattern has not been captured by τ. When the quantumness approaches a final value determined by θ, τ increases rapidly while the rate of τ is nearly invariant.

Next we consider the effect of the environmental memory, which is parameterized by γ, on the QSL timescale. In Fig. 2, τ is evaluated for a fixed initial state (with θ = π/5) and the other parameters except γ and Q. The dependence of τ on the quantumness Q of system and environment is monotonic. The environmental memory timescale is inversely proportional to γ. As an upshot, in the presence of a strongly non-Markovian environment the evolution speed is greatly suppressed, resulting in larger values of τ. Yet it is found that at the end of the dephasing process, the quantum speed limit timescale quickly approaches the same asymptotical value. The difference between the QSL timescale of the system in the extremely non-Markovian environment [τ(γ/Γ = 0.1)] and that in a nearly Markov environment [τ(γ/Γ = 2.0)] is increased with increasing Q before the system goes to the steady state.

Figure 2

Dependence of the quantum speed limit timescale τ on the memory parameter γ in the non-Markovian pure-dephasing dynamics as a function of quantumness Q (θ = π/5).

The inset shows the bound τ derived from the fidelity decay, which fails to capture the correct behavior.

For an n-qubit system in a common dephasing environment, we can rigorously discuss the scaling behavior of QSL for certain states. By a treatment in the Kraus representation45, a general GHZ state evolves into ρ = C(t)○ρ0. Here ○ denotes the entry-wise product and effectively C(t) (as well as ρ) can be expressed in a 2 × 2 matrix expanded by |0⊗〉 and |1⊗〉, where the off-diagonal terms are with and the diagonal terms are unity. By Eqs (11) and (12), we can find that when is sufficiently small (e.g., with a Markov environment, r = e−2Γ is sufficiently close to unity in the short time limit), both the quantumness Q and QSL time τ scale with the number of qubits n as n2.

When L = σ−, the total Hamiltonian describes a dissipation (energy relaxation) model, whose exact super-operation is found to be

where and ∂p(t, s) = P(t)p(t, s) with p(s, s) = 1. Starting from the same pure state as that of the pure-dephasing model, here the time-evolving density matrix satisfies and , where is a complex function of time. In the Markov limit, P(t) = Γ/2 and then ξ = Γτ/2. In the non-Markovian situation, P(t) satisfies ∂P(t) = Γγ/2 − γP(t) + P2(t) with P(0) = 0. Consequently, according to Eq. (3), it is found that

where , , and d ≡ d(t) = Im[P(t)e−]. Note here is not allowed to be zero, otherwise, Q(ρ0, ρ) will vanish according to its definition in Eq. (1). Equations (14) and (15) indicate that in the dissipation model, it is hard to find a closed analytical expression for τ, and one has to resort to the numerical evaluation.

In Fig. 3, we demonstrate the dependence of the QSL timescale on the environmental memory parameter γ, measured in units of the system-environment coupling strength Γ, for a fixed initial state. From the numerically exact dynamics, we find that τ monotonically decreases with increasing γ. With a nearly Markovian environment (see e.g., the dot-dashed line for γ/Γ = 2.0), τ approaches a steady value. As expected in an environment with short memory time, the energy dissipated into the environment from the system has nearly no chance to come back to the system. The dissipation process becomes therefore irreversible. This favors the evolution of the system towards a final incompatible state. As a result, two different regimes are observed. For nearly memoryless dynamics, γ/Γ ≥ 1, the QSL timescale is found to rapidly increase as the system approaches the steady state through a roughly exponential decay. Regarding the spectral function G(t, s), a smaller γ then yields a lesser damping rate of the system. In the strong non-Markovian regime 0.1 ≤ γ/Γ < 1, the pattern becomes complex and the QSL timescale appears to be greatly enhanced by decreasing γ. In this regime, it is difficult for the time-evolving state to become classically incompatible with the initial state.

Figure 3

Dependence of the quantum speed limit timescale on the memory parameter γ of the non-Markovian dissipative process as a function of the quantumness Q, for θ = π/4.

Discussion

We have studied the generation of nonclassicality via the quantumness witness defined as the Hilbert-Schmidt norm of the commutator of the initial and the final quantum states, resulting from time evolution. For arbitrary physical processes we have derived a quantum speed limit that sets the minimum timescale τ for the generation of a given amount of quantumness. This novel QSL has been computed and analyzed in a variety of relevant scenarios including unitary evolution, pure dephasing (of both single- and multiple-qubit system), and energy dissipation. In addition, we have discussed the generation of quantumness in non-unitary evolutions, by employing the exact quantum-state-diffusion equations.

While standard quantum speed limits characterizing the fidelity decay become too conservative and even fail to capture the correct dependence of this timescale on the parameters of the system, the new bound is tight and is saturated under pure dephasing dynamics, whether induced by a Markovian or non-Markovian environment.

Method

We consider the time-evolution of the initial density matrix to be described by a master equation of the form

where is the Liouville super-operator. The rate at which quantumness can vary is then exactly given by

As an example, for unitary dynamics, i.e., in a closed system. Using the Cauchy-Schwarz inequality, i.e., and by virtue of , it follows from the definition of quantumness in Eq. (1) that

To derive a quantum speed limit we integrate from t = 0 to t = τ. Note that Q(ρ0, ρ0) = 0, and . As an upshot, the time in which quantumness can emerge is lower-bounded by

where the time-average is denoted by . We note that even if is explicitly time-independent, i.e., the parameters in the equation of motion are constants, then can not be reduced to since ρ is a function of time.

It is worth pointing out that Eq. (19) suggests

as an upper bound for the speed of evolution of quantumness. Clearly, this quantity can be further upper bounded using the triangular and Cauchy-Schwarz inequalities by . The resulting bound closely resembles the QSL derived by studying the reduced dynamics of an open quantum system in terms of the fidelity decay303132. We note however that this bound is more conservative than that given by τ in Eq. (19). Weaker bounds could be derived as well exploiting the fact that , or conversely using the adjoint of the generator 30. We shall not pursue this goal here.

Additional Information

How to cite this article: Jing, J. et al. Fundamental Speed Limits to the Generation of Quantumness. Sci. Rep. 6, 38149; doi: 10.1038/srep38149 (2016).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.