Simulating Open Quantum Systems with Hamiltonian Ensembles and the Nonclassicality of the Dynamics
Abstract
The incoherent dynamical properties of open quantum systems are generically attributed to an ongoing correlation between the system and its environment. Here, we propose a novel way to assess the nature of these systemenvironment correlations by examining the system dynamics alone. Our approach is based on the possibility or impossibility to simulate opensystem dynamics with Hamiltonian ensembles. As we show, such (im)possibility to simulate is closely linked to the systemenvironment correlations. We thus define the nonclassicality of opensystem dynamics in terms of the nonexistence of a Hamiltonianensemble simulation. This classifies any nonunital opensystem dynamics as nonclassical. We give examples for opensystem dynamics that are unital and classical, as well as unital and nonclassical.
Introduction.—When a quantum system interacts with its environment, its dynamical behavior will, in general, deviate from the dynamics of a strictly isolated one Leggett et al. (1987); Breuer and Petruccione (2002); Weiss (2012); Zhang et al. (2012); Xiong et al. (2015); Chen et al. (2015b). As a result of an ongoing bipartite correlation arising from the systemenvironment interaction, the system dynamics may display incoherent characteristics, such as dephasing or damping processes. Formally, such processes are captured by quantum master equations, replacing the von Neumann equation for isolated systems.
However, incoherent dynamics can also arise as a consequence of a purely classical averaging procedure over distinct autonomous evolutions. For example, the double slit experiment can, when exposed to a disordered potential and after averaging, encounter similar decoherence as if whichslit information had leaked into an environment Gneiting et al. (2016). In this sense, disordered quantum systems described by Hamiltonian ensembles can behave in an analogous manner as open quantum systems, even if individual realizations are strictly isolated (Fig. 1) Gneiting et al. (2016); Kropf et al. (2016); Gneiting and Nori (2017, 2017b).
Here, we exploit this dynamical correspondence to assess the nature of the systemenvironment correlations in terms of the system properties alone. As we show, the impossibility to simulate is necessarily linked to nonclassical systemenvironment correlations. On the other hand, if such a simulation is possible, then there always exists a systemenvironment model which reproduces the system dynamics by relying only on classical correlations. This leads us to defining the nonclassicality of opensystem dynamics in terms of the nonexistence of a Hamiltonianensemble simulation.
Alternative definitions for the nonclassicality of system dynamics have been proposed RahimiKeshari et al. (2013); Sabapathy (2016). In these definitions, the dynamics is considered classical if the state preserves classicality during the temporal evolution. Typically, the classicality of states in these approaches is formulated in terms of the Wigner function or the GlauberSudarshan representation Wigner (1932); Glauber (1963); Sudarshan (1963); Miranowicz et al. (2010); Bartkowiak et al. (2011); Miranowicz et al. (2015). While these definitions also rely on system properties alone, their applicability is limited to systems amenable for such a phase space description, excluding other cases of interest. We, instead, propose to discuss the nonclassicality of opensystem dynamics separately, based on the systemenvironment correlations and independent of the nature of the system.
As an immediate consequence of our definition, all nonunital dynamics, e.g., dissipative processes, are nonclassical. The spinboson model, in contrast, which displays unital dynamics and—on the level of the model—quantum correlations, can, as we find, be simulated and thus exhibits classical opensystem dynamics according to our definition, since these correlations cannot be certified by considering the system dynamics alone. In the case of an extended spinboson model, however, where the environment is complemented by a second qubit and the system dynamics remains unital, we prove the nonexistence of a simulating Hamiltonian ensemble for certain spectral densities; i.e., the system’s evolution is in these cases witnessed to be manifestly nonclassical.
Dynamics of Hamiltonian ensembles.—An isolated quantum system is described by a Hamiltonian ensemble (HE) , if the autonomous Hamiltonian of the system is drawn from a probability distribution [see Fig. 1(a)]. Such HEs are applicable to describing disordered quantum systems. Here, we relate HEs to open quantum systems.
The dynamics of the ensemble averaged state exhibits features distinct from the dynamics of any single realization. The latter is governed by the unitary evolution , with the initial state and , whereas the dynamics of the averaged state is given by the unital (i.e., identity invariant) map
(1) 
Note that an evolution equation for cannot be reduced to some effective Hamiltonian alone but must, in general, take the form of a quantum master equation Gneiting et al. (2016); Kropf et al. (2016) (see Xu et al. (2013) for an experimental implementation of this).
A seminal and instructive example considers a single qubit subject to spectral disorder; i.e., the Hamiltonians in the ensemble differ only in their eigenvalues, while they share a common basis of eigenstates Kropf et al. (2016). The HE may be given by , with the probability distribution kept general. The resulting master equation reads
(2) 
where the effective energy and the decoherence rate follow from the dephasing factor
(3) 
Depending on the underlying probability distribution , the master equation (2) can range from timeconstant dephasing to a strongly oscillating incoherent behavior, the latter even giving rise to purity revivals Kropf et al. (2016).
It is worthwhile to recall that the occurrence of incoherent dynamics in the case of HEs is a consequence of the averaging procedure. Nevertheless, it is reminiscent of open quantum systems, where, in contrast, an ongoing correlation between the system and environment gives rise to the incoherent dynamics. This Letter explores the possibility to simulate open quantum systems with HEs, and vice versa, and the implications on the systemenvironment correlations.
Simulating open quantum systems with Hamiltonian ensembles.—We now show that nonclassical systemenvironment correlations are necessarily linked to the impossibility to simulate the opensystem dynamics with a HE. To this end, we show—conversely—that, if system and environment are persistently classically correlated, then the reduced system state is described by a HE.
A systemenvironment arrangement is characterized by an autonomous total Hamiltonian , with the system , the environment , and the interaction Hamiltonian [see Fig. 1(b)]. The total system evolves unitarily as , with . We say that an open system is described by a HE if the reduced system state allows a decomposition of the form (1), where the probabilities and the Hamiltonians of the ensemble are determined by and the initial state .
Instead of further specifying the total Hamiltonian , we now assume that the total state remains at all times classically correlated, displaying neither quantum discord Ollivier and Zurek (2001); Dakić et al. (2010) nor entanglement. Under this condition, we argue that, for every initial state of the form (with a timeindependent probability distribution and a basis of the environment), the reduced system state can be described by a HE.
The detailed proof is presented in the Supplemental Material sup (). Here, we outline the central steps. First, as a direct consequence of the classical correlations, there exists an environmental basis (in general different from ), such that
(4) 
where the operators act on the system and satisfy for each .
To demonstrate that the are unitary, we again use the zerodiscord assumption, which implies that each environmental offdiagonal term vanishes, i.e., for . This, in turn, implies that there exists a bijection between and such that is nonzero only when its two indices match the bijection, i.e., . Unitarity of the then follows directly.
Finally, we address the time dependence of the . Expressing the bijection as a unitary operator and safely neglecting the index , we can recast in a separable form:
(5) 
The group properties of are thus inherited by the operators and ; i.e., due to the time independence of the total Hamiltonian, we can write and , with timeindependent Hermitian operators and realvalued constants. Consequently, Eq. (1) corresponds to a timeindependent HE when tracing over the environment.
Simulating Hamiltonian ensembles with open quantum systems.—The impossibility to simulate an open system with a HE certifies the quantum nature of the systemenvironment correlations. Notably, this is achieved by considering system properties alone, i.e., without explicit reference to the environment. We now show that, on the other hand, the existence of a simulating HE always admits the possibility of classical systemenvironment correlations; i.e., the latter cannot be excluded by considering only the system.
We explicitly construct a systemenvironment arrangement which reproduces an arbitrary HE relying only on classical correlations. To this end, we write (with the average ) and choose the interaction to be of the form ; i.e., it associates with each of the ensemble a distinct state of an (arbitrary) basis of the environment. Note that the index is generic and may be continuous and/or a multiindex. The environment must then be chosen appropriately to accommodate the complexity of the HE. Moreover, we take the system Hamiltonian to be the average and the bath Hamiltonian to be diagonal in the same basis as , i.e., .
With a separable initial state, , and (i.e., , and the probabilities of the Hamiltonian ensemble are assigned to the environmental populations), the timeevolved total state reads , with . Rewriting , with , we obtain
(6) 
If we now trace over the environment, , we recover the desired decomposition (1) in terms of the HE . Moreover, it is easy to see that the total state (6) is exclusively classically correlated, as desired.
As an instructive example, we consider a pair of qubits coupled to each other via a controlledNOT gate, where a control (C) qubit determines the operation on a target (T) qubit. If the state of the C qubit resides in the classical mixture Chen et al. (2015a), the reduced dynamics of the T qubit will be described by the mixture of evolutions , with and the coupling strength. We thus recover the HE . In this example, the C qubit plays the role of environment, and the qubit pair is at most classically correlated.
Nonclassicality of the dynamics.—It appears natural to classify opensystem dynamics according to their correlation with the environment; i.e., if the system and environment are merely classically correlated, the dynamics may be considered classical; if they are quantum correlated, one may call the dynamics nonclassical. In most cases, however, one does not have (full) access to the environment, rendering such an immediate definition problematic.
We now suggest to classify opensystem dynamics by the (im)possibility to describe the system dynamics by a HE. On the one hand, this definition relies only on system properties, as desired from a practical point of view. On the other hand, as we have shown, it directly links to the systemenvironment correlations, as desirable from a conceptual perspective. Whenever such a simulation exists, it is impossible to exclude classical systemenvironment correlations by knowledge of the system dynamics alone, and we call the latter classical. If the simulation does not exist, quantum correlations must be involved; hence, the dynamics is nonclassical.
As a direct consequence of our definition, any nonunital dynamics is classified nonclassical—a simulating HE is manifestly excluded, certifying the presence of quantum correlations. This includes, e.g., dissipative processes such as the spontaneous decay of an atom. On the other hand, according to our operational definition, we may even call an opensystem dynamics classical if the actual systemenvironment correlations are quantum. This is because our approach is deliberately ignorant of the actual environment and relies only on the possibility to explain the system dynamics with classical correlations. Next, we give an example for this.
Simulating the spinboson model.—We now show that the system dynamics of the spinboson model can be simulated by a HE, even though the actual model displays quantum correlations Bera et al. (2014); IlesSmith et al. (2014). The spinboson model
(7) 
has been extensively studied and is analytically solvable Breuer and Petruccione (2002). Tracing over the environment, the qubit system exhibits pure dephasing dynamics characterized by the dephasing factor
(8) 
In contrast to Eq. (3), which results from averaging over a HE, the dephasing factor (8) incorporates the information of the interaction and the environment into , where is the environmental spectral density. The above solution assumes that the initial state is a direct product, and that the environment is initially thermalized at temperature .
We now construct a simulating HE. In view of Eq. (2), we deduce that individual member Hamiltonians in the ensemble must be of the form , which leaves us with determining the corresponding probabilities. Given a probability distribution , the averaged dynamics can be determined by Eq. (3). Conversely, the underlying distribution function leading to a specific dephasing factor (8) is obtained via the inverse Fourier transform
(9) 
It is clear that the effect of is merely to shift .
To be a legitimate probability distribution function, the resulting in Eq. (9) must be normalized [], real [], and positive []. Normalization is easily seen, since follows from the fact that the pure dephasing dynamics, characterized by Eq. (8), should be completely positive and trace preserving. We therefore have . Moreover, since one is generically interested in the dynamical properties only for , we can deliberately extend the time domain to the full real axis such that is even and . This guarantees that is real: .
The positivity of is less obvious, due to the sinusoidal factors of the integrand in Eq. (9). In the following, we invoke Bochner’s theory Bochner (1933) to prove the general positivity of . To this end, we first introduce the notion of positive definiteness. A function is called positive definite if it satisfies for any finite number of pairs . Note that positive definiteness of a function is different from a positive function, since the latter may not necessarily be positive definite and vice versa. Rather, it corresponds to the positive semidefiniteness of a Hermitian matrix , formed by the function values in accordance with a certain set of indices . As one can show, in Eq. (8) is indeed positive definite. The proof is given in Ref. sup ().
Bochner’s theorem states that a function , defined on , is the Fourier transform of unique positive measure with density function , if and only if is continuous and positive definite Loomis (1953); Rudin (1990). We can thus conclude that in Eq. (8) is the Fourier transform of a certain valid probability distribution [Eq. (9)], i.e., an analog to Eq. (3).
In summary, we have proven that there exists a unique HE, , which simulates the system dynamics exactly, irrespective of the spectral density and the associated, possibly intricate systemenvironment entanglement. We thus call this dynamics classical.
Extended spinboson model.—Unless the system dynamics is nonunital, proving the nonexistence of a simulating HE is, in general, a nontrivial task. We now accomplish this for an extended spinboson model, at the same time deducing the presence of quantum correlations from system properties alone.
Our model consists of two qubits coupled to a common boson environment. The system and the interaction Hamiltonian are replaced by (, ) and , respectively, while the environment Hamiltonian is kept as in Eq. (Simulating Open Quantum Systems with Hamiltonian Ensembles and the Nonclassicality of the Dynamics). Note that the two qubits do not interact directly. The coupling constants are, in general, complex numbers. In order to reveal the nonclassical effects caused by their relative phase, we assume, for simplicity, that they have the same amplitude, i.e., .
In the interaction picture, the total system evolves according to , with the timeordering operator, , and , respectively. In contrast to the conventional spinboson model, time ordering plays a nontrivial role here Reina et al. (2002). (For details, see Supplemental Material sup ().)
In the following, we regard one qubit as the system and the other as part of the environment. The reduced dynamics of the system qubit is then pure dephasing with the dephasing factor [cf. Eq. (8)]
(10) 
where
(11)  
In the second line, we have manually inserted . This ensures that and . The presence of , however, will, in general, result in the violation of positivity. Note that, similar to the conventional spinboson model, individual member Hamiltonians in the HE must be of the form , which allows us to follow the same line of argument.
To demonstrate this violation, we consider the Ohmic spectral density in the zerotemperature limit and a degenerate system Hamiltonian, i.e., . In Fig. 2(a), we depict the legitimate probability distribution for the conventional spinboson model at (blue curve) and (red curve), while in Fig. 2(b), we show for our extended model with (solid curves) and (dashed curves). The latter display a manifest violation of positivity. In Fig. 2(c), we show the landscape of negative contributions to against and for . The gray dashed lines highlight and , chosen in Fig. 2(b).
Conclusions.—We propose a way to classify opensystem dynamics according to their systemenvironment correlations, i.e., if the latter are classical or quantum. As we showed, this can be tested by knowledge of the system evolution alone, based on the (im)possibility to simulate the opensystem dynamics with a Hamiltonian ensemble. According to our definition, any nonunital dynamics is nonclassical. Some unital system evolutions, however, such as in the spinboson model, are classified as classical, even though the model displays quantum correlations. This highlights the operational nature of our definition.
With the extended spinboson model, we provide an example for unital dynamics which is nonclassical according to our definition. Let us note that one may be able to simulate a larger class of unital dynamics with timedependent Hamiltonian ensembles. It is, for example, known that, in the case of qubits, any unital dynamics can be simulated with an ensemble of timedependent Hamiltonians, albeit only if also the probabilities are allowed to be time dependent Landau and Streater (1993); Audenaert and Scheel (2008). However, in the case of autonomous systemenvironment arrangements, which we consider here, such generalization appears unjustified. Finally, let us remark that demonstrating the nonexistence of a Hamiltonianensemble simulation is, in the case of unital evolutions, in general, a nontrivial task. An equivalent but simpler test appears desirable.
Acknowledgements.—This work is supported partially by the National Center for Theoretical Sciences and Ministry of Science and Technology, Taiwan, Grants No. MOST 1032112M006017MY4 and No. MOST 1052811M006059, the MURI Center for Dynamic MagnetoOptics via AFOSR Award No. FA95501410040, the Japan Society for the Promotion of Science (KAKENHI), the IMPACT program of JST, JSPSRFBR Grant No. 175250023, CREST Grant No. JPMJCR1676, RIKENAIST Challenge Research Fund, and the Sir John Templeton Foundation.
Supplemental material: Simulating Open Quantum Systems with Hamiltonian Ensembles and the Nonclassicality of the Dynamics
Appendix A TimeIndependent Hamiltonian Ensemble
In the following, we elaborate in detail the proof that classical bipartite correlations allow for a timeindependent Hamiltonian ensemble decomposition of the reduced system dynamics. In this proof, we do not assume a specific form of the total Hamiltonian . To be precise, we now show that, if the total Hamiltonian is timeindependent, and if, for every initial state of the form (with being any timeindependent probability distribution), the timeevolved total state is always classically correlated between system and environment, displaying neither quantum discord nor entanglement at any time, then the reduced system dynamics admits a timeindependent Hamiltonian ensemble decomposition.
Proof.
Due to the zerodiscord assumption, there exists an environmental basis (in general timedependent and different from ), such that
(1) 
where are operators acting on the system Hilbert space satisfying for each . At this point, we are not yet clear about the timedependence of and nor the unitarity of .
Crucially, the condition
(2) 
should hold for any , due to the zerodiscord assumption. Therefore each term in the above equation vanishes individually. The only possibility to reconcile Eqs. (1) and (2) is the existence of a specific bijection between and , such that for each , i.e., is nonzero only when its two indices match the bijection. Then the unitarity of can then be confirmed according to
(3) 
The bijection between and can be expressed in terms of a unitary operator , such that . The unitary evolution operator can then be recast in a separable form,
(4) 
In the following discussion, we can, in order to keep the notation simple, safely neglect the index .
Since forms a group isomorphism on , we have the oneparameter group property
(5) 
for and infinitesimal . Due to the unitarity of , it can be expressed in terms of an Hermitian generator in the Lie algebra on the environmental Hilbert space such that . Together with Eq. (4), the left hand side of Eq. (5) can be written as
On the right hand side of Eq. (A), we expand around to first order in . Notably, since we do not know the timedependence and commutativity of at this point, we can only achieve a formal expansion in Eq. (A).
Meanwhile, the right hand side of Eq. (5) reads
We again expand around . However, unlike the formal expansion in Eq. (A), we now obtain an explicit expansion in Eq. (A), since commutes with any operator.
Comparing Eqs. (A) and (A), we conclude from their first terms that the group property holds and, combined with the unitarity inferred in Eq. (3), that timeindependent Hermitian operators exist, such that , as well.
To reconcile the second terms of Eqs. (A) and (A), should be diagonalized in the basis , such that , with real parameters . Moreover, should satisfy
(8) 
To guarantee its validity, should commute with , since the latter is the generator of . Consequently, the timedependence of each can be of first order, such that , with real constants . ∎
Consequently, the total state in Eq. (1) can be rewritten as
(9) 
with , which corresponds to a timeindependent Hamiltonian ensemble when tracing over the environment.
Finally, let us remark that, while we restrict ourselves to a timeindependent total Hamiltonian, some of our conclusions can be easily generalized to the timedependent case. This is because Eqs. (14) are consequences of the zerodiscord assumption alone, regardless of the timedependence of the total Hamiltonian. Therefore, we can also achieve the ensemble form with timevarying member Hamiltonians for a timedependent total Hamiltonian. However, as discussed in the main article, in the case of autonomous systemenvironment arrangements, i.e., in the absence of external control, such generalization appears unjustified.
Additionally, we note that, for the case of timeindependent total Hamiltonians, the separable form (4) not only guarantees a persistently classically correlated total state, but also keeps the environmental basis intact without rotation, up to a phase angle .
Appendix B Positive Definiteness
Here, we present the proof of the positive definiteness of the dephasing factor . For completeness, we recall the definition of positive definiteness.
Positive definiteness: A function defined on is called positive definite if it satisfies
(10) 
for any finite number of pairs .
We now show that, if defines a CPTP pure dephasing dynamics, and if is even and , then defined on is positive definite.
Note that describing a CPTP pure dephasing dynamics implies that , , and for any . This means that the coherence of the system can never exceed its initial value. These properties will be frequently used in the following proof.
Proof.
To simplify the problem, we first observe that the positive definiteness of is equivalent to that of , since
(11) 
Correspondingly, we can assume that without loss of generality.
Since Eq. (10) must be valid for any number of pairs, we give the proof in an inductive manner.
In the case of only one pair , Eq. (10) is trivially satisfied. We therefore start with the case of two pairs. As stated in the main article, Eq. (10) is equivalent to the positive semidefiniteness of the Hermitian matrix:
(12) 
It is automatically satisfied according to the CPTP dynamics defined by .
We proceed to show the positive semidefiniteness of the Hermitian matrix
(13) 
for the case of three pairs. In the above matrix, and hereafter, the abbreviation has been adopted. Since is threedimensional, it is generically hard to write down an analytic expression for its three eigenvalues . Nevertheless, analyzing its characteristic polynomial gives us substantial knowledge on the eigenvalues:

follows from the invariance of the trace.

equals to the sum of all principal minors of of order and is consequently nonnegative, since each principal minor is nonnegative, following the positive semidefiniteness of .

. The positivity of the product of eigenvalues is verified with the help of a simple geometric visualization shown in Fig. 1. Explicitly expanding the determinant leads to
(14) with the notation , , and . This can be interpreted in terms of a triangle with circumcircle (dashed circle) of diameter less than . With the help of , the positivity of Eq. (14) and, consequently, of the product of eigenvalues is then inferred.
Combining (i)(iii), we can conclude that the three eigenvalues are nonnegative each and, therefore, that is positive semidefinite.
Before proceeding to the case of four pairs, it is worthwhile to discuss how the minimum of Eq. (10) is achieved. For the case of three pairs, the LHS of Eq. (10) is equivalent to the summation over entries in the following array:
(15) 
If we first determine the amplitudes and adjust their arguments and , it is clear that the diagonal elements in the array (15) are all positive and, to reduce the resulting summation, the possible negative contributions are given by the offdiagonal elements. If we choose three pairs such that the angles between any two within them is less than , as show in Fig. 2(a), the summation of all offdiagonal elements is positive. Therefore, we must choose appropriate pairs such that all their relative arguments strictly exceed , as show in Fig. 2(b). To maximize the negative contributions, we assume and . We therefore draw the conclusion that for the case of maximized relative arguments between three .
However, in the case of four pairs, it is impossible to insert the fourth such that all relative arguments are strictly larger than , as shown in Fig. 2(c). According to the above discussion, to deal with the array of four pairs,
(16) 
we can at most group three with all three relative arguments strictly larger than by setting their corresponding equal. Then the array (16) reduces to a simpler one:
(17) 
Again, in accordance with the positive semidefiniteness of , we can guarantee the validity of Eq. (10) for the case of four pairs.
For the case of five or more pairs, a similar procedure can be applied to continuously reduce the problem to an equivalent or case. This implies the validity of Eq. (10) for the general case. ∎
Let us remark that the above proof already indicates the general impossibility of a Hamiltonian ensemble description for arbitrary pure dephasing dynamics. Many conclusions in the above proof hold since the phase angle of is directly proportional to time . This is particularly manifest in Eq. (11). However, this is in general not the case, e.g., in the extended spinboson model below. We consequently may obtain invalid (or quasi) distributions in the extended spinboson model.
Appendix C Extended SpinBoson Model
We proceed with the details of the extended spinboson model, which consists of two qubits coupled to a common boson environment. The system and the interaction Hamiltonian of the conventional spinboson model are thus replaced by
(18) 
Note that the two qubits do not interact with each other directly. Let us remark that, while we consider two qubits here, our treatment can straightforwardly be generalized to more than two qubits.
Transforming to the interaction picture with respect to , the total system evolves according to the unitary evolution operator
(19) 
where is the timeordering operator, , and , respectively. In the conventional spinboson model with a single qubit, timeordering plays no significant role, since it merely introduces a global phase to the unitary evolution operator. However, this is not the case for extended models with more than one qubit, where one must carefully deal with the effect of timeordering . We therefore have
(20) 
with
By using the prescription , the operator can be recast into