Frozen Gaussian sampling algorithms for simulating Markovian open quantum systems in the semiclassical regime

Simulating Markovian open quantum systems in the semiclassical regime poses a grand challenge for computational physics, as the highly oscillatory nature of the dynamics imposes prohibitive resolution requirements on traditional grid-based methods. To overcome this barrier, this paper introduces an efficient Frozen Gaussian Sampling (FGS) algorithm based on the Wigner-Fokker-Planck phase-space formulation. The proposed algorithm exhibits two transformative advantages. First, for the computation of physical observables, its sampling error is independent of the semiclassical parameter $\varepsilon$, thus fundamentally breaking the prohibitive computational scaling faced by grid methods in the semiclassical limit. Second, its mesh-free nature entirely eliminates the boundary-induced instabilities that constrain long-time grid-based simulations. Leveraging these capabilities, the FGS algorithm serves as a powerful investigatory tool for exploring the long-time behavior of open quantum systems. Specifically, we provide compelling numerical evidence for the existence of steady states in strongly non-harmonic potentials-a regime where rigorous analytical results are currently lacking.

💡 Research Summary

This paper addresses the formidable computational challenge of simulating Markovian open quantum systems in the semiclassical regime (where the scaled Planck constant ε ≪ 1). In this limit, the highly oscillatory nature of quantum dynamics imposes prohibitive spatial and temporal resolution requirements on traditional grid-based numerical methods, leading to an exponential scaling of cost with 1/ε. To overcome this barrier, the authors introduce a novel and efficient Frozen Gaussian Sampling (FGS) algorithm based on a phase-space formulation using the Wigner-Fokker-Planck (WFP) equation.

The core innovation lies in approximating the system’s Wigner function as a superposition of evolving Gaussian wave packets. Instead of discretizing space on a fixed grid, the method uses a Monte Carlo approach to sample these wave packets. Each sampled packet’s parameters—its phase-space center (q, p), its covariance matrix Σ (describing its shape), and its amplitude A—are propagated in time according to a set of ordinary differential equations (ODEs) rigorously derived from the WFP equation via asymptotic analysis. This derivation, a key theoretical contribution of the work, fully incorporates the effects of dissipation, diffusion, and drift from the environment into the wave packet dynamics.

The proposed FGS algorithm offers two transformative advantages over conventional methods. First, its sampling error for computing physical observables is essentially independent of the semiclassical parameter ε. This fundamentally breaks the curse of dimensionality faced by grid-based techniques, allowing for efficient simulation deep into the semiclassical limit where ε → 0. Second, as a mesh-free method that tracks trajectories in an unbounded phase space, it entirely eliminates the boundary-induced instabilities that plague long-time simulations on finite grids. This makes the algorithm exceptionally robust for studying long-time and steady-state behaviors.

The authors provide a crucial mathematical proof ensuring the stability of their framework: they demonstrate that the positive definiteness of the covariance matrix Σ(t) is preserved throughout the time evolution, guaranteeing that the Gaussian ansatz remains physically well-defined.

Leveraging these capabilities, the paper employs the FGS algorithm not only as a numerical solver but also as a powerful investigative tool. It presents compelling numerical evidence for the existence of steady states in open quantum systems subject to strongly anharmonic potentials (e.g., a double-well potential)—a regime where rigorous analytical results are currently lacking. The simulations clearly show convergence to the same steady-state Wigner distribution from different initial conditions.

In summary, this work develops a comprehensive theoretical and computational framework for semiclassical open quantum dynamics. The FGS algorithm achieves a breakthrough in computational efficiency and long-time stability, opening the door to accurate simulations of complex dissipative quantum phenomena that were previously intractable. It bridges advanced semiclassical analysis with practical numerical computation, offering a new paradigm for exploring the rich behavior of open quantum systems.