In closed quantum dynamics, energy eigenstates represent physical excitations which, for systems with a classical limit, localize on quantized orbits in phase space. Such a picture is lacking in open systems: Liouvillians have distinct right and left eigenstates that decay with time and are difficult to interpret. We recover a physical correspondence by introducing a gauge-invariant quasiprobability measure combining right and left eigenstates that interprets a Liouvillian mode as a set of coherences, allowing one to visualize the modes as quasiprobability distributions in a ``doubled'' phase space. Using this measure we show that, remarkably, a damped oscillator retains its quantized ``orbits'' for a large class of linear and nonlinear loss, with measurable signatures. A thermal bath broadens these orbits, similar to energy levels, leading to a characteristic multistage relaxation. This provides an experimentally accessible open-system analog of the eigenstate-orbit correspondence underlying semiclassical quantization. For quadratic systems, our results yield dissipative analogs of invariant tori - a central concept in Hamiltonian systems.
In closed quantum dynamics, energy eigenstates represent physical excitations which, for systems with a classical limit, localize on quantized orbits in phase space. Such a picture is lacking in open systems: Liouvillians have distinct right and left eigenstates that decay with time and are difficult to interpret. We recover a physical correspondence by introducing a gauge-invariant quasiprobability measure combining right and left eigenstates that interprets a Liouvillian mode as a set of coherences, allowing one to visualize the modes as quasiprobability distributions in a doubled'' phase space. Using this measure we show that, remarkably, a damped oscillator retains its quantized orbits'' for a large class of linear and nonlinear loss, with measurable signatures. A thermal bath broadens these orbits, similar to energy levels, leading to a characteristic multistage relaxation. This provides an experimentally accessible open-system analog of the eigenstate-orbit correspondence underlying semiclassical quantization. For quadratic systems, our results yield dissipative analogs of invariant tori - a central concept in Hamiltonian systems.
We develop an approach for understanding the dynamics of open quantum systems by analyzing individual quantum trajectories in the eigenbasis of the Liouvillian superoperator. From trajectory-eigenstate overlaps, we construct a quasiprobability distribution that characterizes the degree of localization of the trajectories in the Liouvillian eigenbasis. Contrary to the common wisdom that late-time dynamics are governed solely by the steady state and the slowest-decaying modes, we show that trajectories can remain well spread over transient eigenstates deep within the bulk of the Liouvillian spectrum even at late times. We demonstrate this explicitly using numerical simulations of interacting spin chains and bosonic systems. Moreover, we find that the delocalization of the trajectory strongly correlates with the purity of the trajectory-averaged steady state, establishing a further link between the trajectory and ensemble pictures of open quantum dynamics.
The classical-quantum correspondence is developed in its most natural setting, i.e., that of phase space distributions. Specifically, the Wigner–Weyl representation of the two fundamental quantum densities ‖n≳<n‖ and ‖n≳<m‖, which are eigenstates of the quantum Liouville operator, are shown, for integrable systems, to yield stationary and nonstationary classical Liouville eigenstates in the ℏ→0 limit. The former are uniform distributions on semiclassically quantized tori, whereas the latter are nonuniform distributions on similar tori at intermediate actions. Applications to semiclassical mechanics in action-angle variable, the Heisenberg correspondence for matrix elements and the computation of vibrational intensities are provided.
Open quantum systems lose energy and phase coherence through different dissipative processes, but these processes can produce overlapping dynamical signatures. The Liouvillian spectrum summarizes how such a system relaxes, yet it is not obvious how much of that spectrum is needed to distinguish the underlying dissipation rates. We study this question for amplitude damping and dephasing in a six-qubit Lindblad model whose spectrum can be derived analytically. We retain only the slowest non-steady spectral modes and ask how many are required before each dissipative rate becomes recoverable. We show that population modes contain no dephasing information, which creates a lower bound of D = 2^n retained modes for uniform dephasing identifiability in the relevant rate regime. The measured recovery threshold reaches this bound at n = 4,5,6, while n = 3 remains above it. At n = 6, least squares achieves a mean joint absolute error of order 10^-9, compared with 4.355 x 10^-4 for four tabular learning methods. Robustness tests show that this advantage weakens when the spectra are perturbed and when a transverse field breaks the commuting structure. These results show that the amount and structure of retained spectral information can determine whether dissipative parameters are recoverable, independently of the estimator used. The present conclusions apply to noise-free simulator spectra rather than measurement-derived spectra.
Open quantum systems with nearly degenerate energy levels have been shown to exhibit long-lived metastable states in the approach to equilibrium, even when modelled with certain Lindblad-form quantum master equations. This is a result of dramatic separation of timescales due to differences between Liouvillian eigenvalues. These metastable states often have nonzero coherences that die off only in the long-time limit once the system reaches thermal equilibrium. We examine two distinct situations that give rise to this effect: one in which dissipative dynamics couple together states only within a nearly degenerate subspace, and one in which they give rise to jumps over finite-energy splittings, between nearly degenerate subspaces. We find, in each case, that a change of basis can often lead to a representation that more naturally captures the impact of the system-bath interaction than does the energy eigenbasis, revealing that separate timescales are associated with separate processes (e.g., decoherence into a nonenergy eigenbasis, decay of population correlations to the initial state). This approach is paired with the inspection of quantum trajectories, which further provide intuition as to how open system evolution is characterized when coherent oscillations, thermal relaxation, and decoherence all occur simultaneously. Published by the American Physical Society 2024
We consider an open quantum system with dissipation, described by a Lindblad Master equation (LME). For dissipation locally acting and sufficiently strong, a separation of the relaxation timescales occurs, which, in terms of the eigenvalues of the Liouvillian, implies a grouping of the latter in distinct vertical stripes in the complex plane at positions determined by the eigenvalues of the dissipator. We derive effective LME equations describing the modes within each stripe separately, and solve them perturbatively, obtaining for the full set of eigenvalues and eigenstates of the Liouvillian explicit expressions correct at order 1/Γ included, where Γ is the strength of the dissipation. As an example, we apply our general results to quantum XYZ spin chains coupled, at one boundary, to a dissipative bath of polarization.
This record contains the manuscript “Canonical Invariants, Quantum Noise, and Memory in Gaussian Open Systems.” We develop an invariant-based formulation of Gaussian open-system dynamics for a time-dependent quantum harmonic oscillator subject to dissipation and noise. The construction combines canonical invariants with stochastic and retarded-mode descriptions, allowing Markovian thermal relaxation and finite-memory non-Markovian dynamics to be treated within a common framework. A central result is that, away from equilibrium, the instantaneous energy is not determined by the entropy alone. For local relaxation, the energy admits a differential state representation of the form[E=E!\left(S,\dot S;\frac{\omega(t)}{\nu_B(t)},\gamma_{\rm eff}(t)\right),]so that the entropy rate (\dot S) carries essential nonequilibrium information. Finite-memory dynamics enlarges the required state data while preserving this basic structure. The formulation also makes explicit the consistency relation between dissipation and the asymmetric quantum-noise spectrum, connects the Markovian limit with thermal GKLS dynamics, and provides a finite-dimensional localization of colored-noise memory kernels. The framework is intended as a basis for studying driven Gaussian open systems and for possible extensions to quantum-field and electromagnetic settings. This version has been submitted to Physical Review Research.