Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2604.01165

Variational Dynamics of Open Quantum Spin Systems in Phase Space

Jacopo Tosca, Zejian Li, Francesco Carnazza, Cristiano Ciuti

quant-phcond-mat.otherSubmitted Apr 1, 2026• Announced Apr 2, 2026• 0 citations
arXivHTMLPDF

Abstract

We introduce a variational method for simulating the dynamics of interacting open quantum spin systems. The method is based on the spin phase-space representation and variationally targets the Husimi-$Q$ function with an ansatz based on a multi-dimensional mixture of spin-coherent states. Crucially, the mixture coefficients are allowed to take negative values, enabling the faithful capture of quantum correlations beyond semiclassical descriptions. The resulting equations of motion are derived from the Dirac-Frenkel variational principle and can be evaluated efficiently without resorting to Monte Carlo sampling by exploiting the analytical structure of the ansatz. As a first application, we demonstrate that this approach accurately captures both the full quantum dynamics and the non-equilibrium steady states of the transverse-field quantum Ising model, in excellent agreement with exact diagonalization. Furthermore, we show that the method scales efficiently to large two-dimensional lattices, a regime that remains challenging for other techniques.

Comments
8 pages, 4 figures

Related papers

Variational Multi-Gaussian Phase-Space Bosonic Dynamics via Automatic Differentiation

Jacopo Tosca, Francesco Carnazza, Luca Giacomelli +1 more

We introduce a variational method for simulating the dynamics of interacting open quantum bosonic systems deep in the quantum regime. The method is based on a multidimensional Wigner phase-space representation and employs a variational multi-Gaussian , whose accuracy is systematically controlled by the number of Gaussian components. The variational equations of motion are derived from the Dirac-Frenkel principle and evaluated efficiently by combining the analytical structure of Gaussian functions with automatic differentiation. As a key physical application, we study a driven-dissipative two-dimensional Bose-Hubbard lattice with two-boson coherent driving and two-body losses. Using our dynamical approach, we compute the finite-size scaling of the Liouvillian spectral gap—extracted from the relaxation dynamics—which vanishes in the thermodynamic limit. Our results reveal critical slowing down with dynamical exponents of the 2D quantum Ising universality class, demonstrating the power of our method to capture complex quantum dynamics in large open systems.

arXiv:2507.14076•Physical Review X•Jul 9, 2026•0 citations

Variational dynamics of open quantum systems in phase space

Debbie Eeltink, Filippo Vicentini, Vincenzo Savona

We present a method to simulate the dynamics of large driven-dissipative many-body open quantum systems using a variational encoding of the Wigner or Husimi-Q quasi-probability distributions. The method relies on Monte-Carlo sampling to maintain a polynomial computational complexity while allowing for several quantities to be estimated efficiently. As a first application, we present a proof of principle investigation into the physics of the driven-dissipative Bose-Hubbard model with weak nonlinearity, providing evidence for the high efficiency of the phase space variational approach.

arXiv:2307.07429•Jul 14, 2023•0 citations

Variational Multi-Gaussian Quantum Trajectories

Zejian Li, Jacopo Tosca

We propose a variational framework for the efficient simulation of stochastic quantum dynamics in interacting bosonic systems. Our method generalizes the time-dependent variational principle to accommodate quantum state diffusion processes and is tailored for describing quantum trajectories unraveled from a master equation. We adopt a wavefunction ansatz composed of a coherent superposition of Gaussian wavepackets, which allows the fully analytical assembly of the variational equations of motion while capturing both the quantum trajectory dynamics and the Lindblad evolution beyond semiclassical approximations. The method is carefully benchmarked on a Bose-Hubbard dimer, and applied to extended Bose-Hubbard lattices in both 1D and 2D. In particular, we show the emergence of a symmetry-breaking phase transition in 2D lattices, which is absent in 1D chains.

quant-pharXiv:2609.30225•Sep 25, 2026•0 citations

Variational Quantum Monte Carlo Method with a Neural-Network Ansatz for Open Quantum Systems

Alexandra Nagy, Vincenzo Savona

The possibility to simulate the properties of many-body open quantum systems with a large number of degrees of freedom (d.o.f.) is the premise to the solution of several outstanding problems in quantum science and quantum information. The challenge posed by this task lies in the complexity of the density matrix increasing exponentially with the system size. Here, we develop a variational method to efficiently simulate the nonequilibrium steady state of Markovian open quantum systems based on variational Monte Carlo methods and on a neural network representation of the density matrix. Thanks to the stochastic reconfiguration scheme, the application of the variational principle is translated into the actual integration of the quantum master equation. We test the effectiveness of the method by modeling the two-dimensional dissipative XYZ spin model on a lattice.

arXiv:1902.09483•Physical Review Letters•Jun 28, 2019•227 citations

Simulation of the dynamics of many-body quantum spin systems using phase-space techniques

Ray Jia Hong Ng, Erik Schwartz Sørensen, Piotr Deuar

We reformulate the full quantum dynamics of spin systems using a phase-space representation based on SU(2) coherent states which generates an exact mapping of the dynamics of any spin system onto a set of stochastic differential equations. This representation is superior in practice to an earlier phase-space approach based on Schwinger bosons, with the numerical effort scaling only linearly with system size. By also implementing extrapolation techniques from quasiclassical equations to the full quantum limit, we are able to extend useful simulation times severalfold. This approach is applicable in any dimension including cases where frustration is present in the spin system. The method is demonstrated by simulating quenches in the transverse-field Ising model in one and two dimensions.

arXiv:1307.3786•Physical Review B•Oct 4, 2013•18 citations

Accurate Two-Time Correlation from a Pseudomode-Based Variational Approach with the Multiple Davydov D1 Density Operator

Yiying Yan, Zhiguo LÜ, Yuguo Su +2 more

The pseudomode framework provides an efficient description of open quantum systems in non-Markovian and nonperturbative coupling regimes by utilizing simple Markovian Lindblad master equations that incorporate a set of damped bosonic pseudomodes. Developing an efficient numerical scheme to solve these equations accurately is highly desirable. In this work, we propose a modified multiple Davydov D1 density operator ansatz and combine it with the Dirac-Frenkel time-dependent variational principle to simulate the dynamics governed by these master equations. We apply this variational method to the spin-boson model, both with and without a periodic external driving field, and compute both single- and two-time expectation values in non-Markovian, strong-coupling regimes. The accuracy of our variational approach is validated by comparing its predictions with results obtained from the small matrix path integral method. We demonstrate that the variational and path-integral approaches are in excellent agreement for both single- and two-time expectations across nondriven and driven cases, spanning both low- and high-temperature regimes. Our results establish the reliability and broad applicability of the pseudomode-based variational approach for calculating two-time correlation functions in open quantum systems.

Journal of Chemical Theory and Computation•Oct 16, 2026•0 citations
We introduce a variational method for simulating the dynamics of interacting open quantum spin systems. The method is based on the spin phase-space representation and variationally targets the Husimi-QQQ function with an ansatz based on a multi-dimensional mixture of spin-coherent states. Crucially, the mixture coefficients are allowed to take negative values, enabling the faithful capture of quantum correlations beyond semiclassical descriptions. The resulting equations of motion are derived from the Dirac-Frenkel variational principle and can be evaluated efficiently without resorting to Monte Carlo sampling by exploiting the analytical structure of the ansatz. As a first application, we demonstrate that this approach accurately captures both the full quantum dynamics and the non-equilibrium steady states of the transverse-field quantum Ising model, in excellent agreement with exact diagonalization. Furthermore, we show that the method scales efficiently to large two-dimensional lattices, a regime that remains challenging for other techniques.
We introduce a variational method for simulating the dynamics of interacting open quantum bosonic systems deep in the quantum regime. The method is based on a multidimensional Wigner phase-space representation and employs a variational multi-Gaussian , whose accuracy is systematically controlled by the number of Gaussian components. The variational equations of motion are derived from the Dirac-Frenkel principle and evaluated efficiently by combining the analytical structure of Gaussian functions with automatic differentiation. As a key physical application, we study a driven-dissipative two-dimensional Bose-Hubbard lattice with two-boson coherent driving and two-body losses. Using our dynamical approach, we compute the finite-size scaling of the Liouvillian spectral gap—extracted from the relaxation dynamics—which vanishes in the thermodynamic limit. Our results reveal critical slowing down with dynamical exponents of the 2D quantum Ising universality class, demonstrating the power of our method to capture complex quantum dynamics in large open systems.
We present a method to simulate the dynamics of large driven-dissipative many-body open quantum systems using a variational encoding of the Wigner or Husimi-Q quasi-probability distributions. The method relies on Monte-Carlo sampling to maintain a polynomial computational complexity while allowing for several quantities to be estimated efficiently. As a first application, we present a proof of principle investigation into the physics of the driven-dissipative Bose-Hubbard model with weak nonlinearity, providing evidence for the high efficiency of the phase space variational approach.
We propose a variational framework for the efficient simulation of stochastic quantum dynamics in interacting bosonic systems. Our method generalizes the time-dependent variational principle to accommodate quantum state diffusion processes and is tailored for describing quantum trajectories unraveled from a master equation. We adopt a wavefunction ansatz composed of a coherent superposition of Gaussian wavepackets, which allows the fully analytical assembly of the variational equations of motion while capturing both the quantum trajectory dynamics and the Lindblad evolution beyond semiclassical approximations. The method is carefully benchmarked on a Bose-Hubbard dimer, and applied to extended Bose-Hubbard lattices in both 1D and 2D. In particular, we show the emergence of a symmetry-breaking phase transition in 2D lattices, which is absent in 1D chains.
The possibility to simulate the properties of many-body open quantum systems with a large number of degrees of freedom (d.o.f.) is the premise to the solution of several outstanding problems in quantum science and quantum information. The challenge posed by this task lies in the complexity of the density matrix increasing exponentially with the system size. Here, we develop a variational method to efficiently simulate the nonequilibrium steady state of Markovian open quantum systems based on variational Monte Carlo methods and on a neural network representation of the density matrix. Thanks to the stochastic reconfiguration scheme, the application of the variational principle is translated into the actual integration of the quantum master equation. We test the effectiveness of the method by modeling the two-dimensional dissipative XYZ spin model on a lattice.
We reformulate the full quantum dynamics of spin systems using a phase-space representation based on SU(2) coherent states which generates an exact mapping of the dynamics of any spin system onto a set of stochastic differential equations. This representation is superior in practice to an earlier phase-space approach based on Schwinger bosons, with the numerical effort scaling only linearly with system size. By also implementing extrapolation techniques from quasiclassical equations to the full quantum limit, we are able to extend useful simulation times severalfold. This approach is applicable in any dimension including cases where frustration is present in the spin system. The method is demonstrated by simulating quenches in the transverse-field Ising model in one and two dimensions.
The pseudomode framework provides an efficient description of open quantum systems in non-Markovian and nonperturbative coupling regimes by utilizing simple Markovian Lindblad master equations that incorporate a set of damped bosonic pseudomodes. Developing an efficient numerical scheme to solve these equations accurately is highly desirable. In this work, we propose a modified multiple Davydov D1 density operator ansatz and combine it with the Dirac-Frenkel time-dependent variational principle to simulate the dynamics governed by these master equations. We apply this variational method to the spin-boson model, both with and without a periodic external driving field, and compute both single- and two-time expectation values in non-Markovian, strong-coupling regimes. The accuracy of our variational approach is validated by comparing its predictions with results obtained from the small matrix path integral method. We demonstrate that the variational and path-integral approaches are in excellent agreement for both single- and two-time expectations across nondriven and driven cases, spanning both low- and high-temperature regimes. Our results establish the reliability and broad applicability of the pseudomode-based variational approach for calculating two-time correlation functions in open quantum systems.