Double-Bracket Master Equations: Phase-Space Representation and Classical Limit

TL;DR

The paper addresses how double-bracket quantum master equations—including energy dephasing from double commutators and nonlinear, trace-preserving dynamics from double anticommutators—behave in the classical limit. By reformulating the dynamics in phase space with the Wigner function and performing an $\hbar$-expansion, it identifies a classical generator that combines a Poisson bracket flow with dissipative terms, and reveals a gradient-flow structure for these evolutions. It provides explicit analyses for a simple harmonic oscillator and a driven anharmonic oscillator, illustrating decoherence, Wigner negativity, and the quantum-classical transition, and demonstrates how higher-order nested brackets model spectral filtering techniques in quantum many-body contexts. The results offer a geometric interpretation of nonunitary dynamics, enable semiclassical analyses of decoherence and chaos, and suggest extensions to interacting systems and advanced filtering methods with practical relevance to quantum control and open-system dynamics.

Abstract

We investigate the classical limit of quantum master equations featuring double-bracket dissipators. Specifically, we consider dissipators defined by double commutators, which describe dephasing dynamics, as well as dissipators involving double anticommutators, associated with fluctuating anti-Hermitian Hamiltonians. The classical limit is obtained by formulating the open quantum dynamics in phase space using the Wigner function and Moyal products, followed by a systematic $\hbar$-expansion. We begin with the well-known model of energy dephasing, associated with energy diffusion. We then turn to master equations containing a double anticommutator with the system Hamiltonian, recently derived in the context of noisy non-Hermitian systems. For both classes of double-bracket equations, we provide a gradient-flow representation of the dynamics. We analyze the classical limit of the resulting evolutions for harmonic and driven anharmonic quantum oscillators, considering both classical and nonclassical initial states. The dynamics is characterized through the evolution of several observables as well as the Wigner logarithmic negativity. We conclude by extending our analysis to generalized master equations involving higher-order nested brackets, which provide a time-continuous description of spectral filtering techniques commonly used in the numerical analysis of quantum systems.

Figures (7)

• Figure 1: Gradient flow description of the time evolution, showing the orbit associated with the trajectory $\hat{\mathcal{O}}$ in the manifold $\mathcal{M}$. The tangent space, shown as the plane surface with grids, is denoted by $T_{\hat{\sigma}_t}$, at the point $\hat{\sigma}_t$ on $\hat{\mathcal{O}}$, for some parameter $t$ parametrizing the trajectory. An element on $T_{\hat{\sigma}_t}$ is marked as $[\hat{X},\hat{\sigma}_t]$.
• Figure 2: Evolution of the Wigner function of an initial Gaussian wavepacket with mean at $(x,p) = (1,1)$ in the semiclassical limit of double-commutator (top) and double-anticommutator (bottom) master equations ($\Gamma = 0.3$).
• Figure 3: Average position, momentum, and energy of the Wigner function evolved using the semiclassical limit of the master equation with double commutator (left) and double anticommutator (middle). We can also observe the monotonic decrease of energy moments $\mu_2$ and $\mu_4$ (right) for the case of double anticommutator.
• Figure 4: Wigner logarithmic negativity $\mathcal{W}$ and time of classical emergence $t_c$ for different parameter values. Physical parameters in the Hamiltonian are set to $m = 1.0$, $A = 1.0$, $B = 0.1$, $\kappa = 0.2$, and $\omega_d = 1.0$. The initial state is taken to be a Gaussian wavepacket with minimum uncertainty, centered near the right minimum at $x_0 = 2.19$.
• Figure 5: Wigner logarithmic negativity $\mathcal{W}$ for different strengths of the nonlinearity parameter $\kappa \in [0,1]$ plotted for different decoherence strengths $\gamma$ and $\Gamma$ using the Gaussian initial state in the driven anharmonic oscillator.