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:2307.07429

Variational dynamics of open quantum systems in phase space

Debbie Eeltink, Filippo Vicentini, Vincenzo Savona

Published Jul 14, 2023• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2307.07429

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 Spin Systems in Phase Space

Jacopo Tosca, Zejian Li, Francesco Carnazza +1 more

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.

quant-pharXiv:2604.01165•Apr 2, 2026•0 citations

Stochastic and Neural Network methods for many-body open quantum systems

Filippo Vicentini

This thesis is devoted to the theoretical study of driven-dissipative many-body quantum systems with a particular focus on advanced computational methods. After reviewing several state-of-the-art numerical methods to determine the steady-state and dynamics of open quantum systems described by a Lindblad master equation for the density matrix, this thesis manuscript presents two original methods: (i) an optimal unraveling of stochastic quantum trajectories for disordered systems; (ii) a variational neural network method. The stochastic methods have been applied to investigate the driven-dissipative Bose-Hubbard lattice model that can for instance describe lattices of coupled nonlinear electromagnetic resonators such those obtained with semiconductor microcavities or superconducting quantum circuits. Our study has shown that a critical slowing down of the dynamics occurs for some parameters in dimensions higher than one. The phenomena is linked to a first-order dissipative phase transition and to lattice-level collective bistability. Furthermore, we reveal the role of disorder on the spatial correlation properties of these phenomena. In the last part of the thesis, we present the strong connection between fundamental concepts of machine learning and variational approaches in quantum physics. We present a novel computational method for open quantum systems based on neural network ans\"{a}tze that is compatible with a stochastic minimisation of a cost function through Monte Carlo sampling. As a first application, we have applied this method to the dissipative transverse quantum Ising model. Exciting perspectives and possible developments are also drawn in the concluding remarks.

•Dec 11, 2019•1 citation

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

Unraveling Open Quantum Dynamics with Time-Dependent Variational Monte Carlo

Christian Apostoli, Jacopo D'Alberto, Marco G. Genoni +2 more

We introduce a method to simulate open quantum many-body dynamics by combining time-dependent variational Monte Carlo (tVMC) with quantum trajectory techniques. Our approach unravels the Lindblad master equation into an ensemble of stochastic Schrödinger equations for a variational ansatz, avoiding the exponential cost of density matrix evolution. The method is compatible with generic ansatzë, including expressive neural-network wavefunctions. We derive the nonlinear stochastic equations of motion for the variational parameters and employ suitable Stratonovich numerical solvers. To validate our approach, we simulate quenches in the locally dissipative long-range Ising model in a transverse field, accurately capturing non-equilibrium magnetization and spin squeezing dynamics relevant to trapped-ion and Rydberg atom experiments. The framework is computationally efficient, scalable on high-performance computing platforms, and can be readily integrated into existing tVMC implementations. This work enables large-scale simulations of complex, dissipative quantum systems in higher dimensions, with broad implications for quantum technology and fundamental science.

arXiv:2506.23928•Jun 30, 2025•0 citations

Quantum Monte Carlo approach to the non-equilibrium steady state of open quantum systems

Alexandra Nagy

Many-body open quantum systems are exposed to an essentially uncontrollable environment that acts as a source of decoherence and dissipation. As the exact treatment of such models is generally unfeasible, it is favourable to formulate an approximate description by means of the dynamics of the reduced density operator of the system. When studying models with a weak coupling to a memoryless environment, the existence of such solution is granted and the reduced dynamics is governed by the Lindblad quantum master equation. In recent years, open quantum systems have evolved into a major field of studies. Focus of these studies are the characterization of emergent phenomena and dissipative phase transitions, as well as the ongoing debate about whether quantum computing schemes are still hard to simulate classically --- and thus achieve quantum supremacy --- when in presence of some degree of noise-induced decoherence. However, the challenge posed by this task lies in the complexity of the density matrix that increases exponentially with the system size, and the quest for efficient numerical methods is a research field that is still in its infancy. Here, we first develop a real-time full configuration interaction quantum Monte Carlo technique that stems from a class of methods generally known as projector Monte Carlo. The approach enables the stochastic sampling of the Lindblad time evolution of the density matrix thanks to a massively parallel algorithm, thus providing estimates of observables on the non-equilibrium steady state. We present the underlying theory and introduce an initiator technique and importance sampling to reduce statistical error. We demonstrate the efficiency of the approach by applying it to the dissipative two-dimensional XYZ spin-1/2 model on a lattice. As the importance sampling of projector approaches is often combined with variational results, we then introduce a novel method that is based on the 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 Lindblad quantum master equation. We test the effectiveness of the method by modeling the steady state of the dissipative two-dimensional XYZ spin-1/2 model through a dissipative phase transition, and also the real-time dynamics of the dissipative Ising model. In addition, we discuss the application of the developed methods and the open questions of the field.

•Jan 1, 2020•0 citations
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 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 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.
This thesis is devoted to the theoretical study of driven-dissipative many-body quantum systems with a particular focus on advanced computational methods. After reviewing several state-of-the-art numerical methods to determine the steady-state and dynamics of open quantum systems described by a Lindblad master equation for the density matrix, this thesis manuscript presents two original methods: (i) an optimal unraveling of stochastic quantum trajectories for disordered systems; (ii) a variational neural network method. The stochastic methods have been applied to investigate the driven-dissipative Bose-Hubbard lattice model that can for instance describe lattices of coupled nonlinear electromagnetic resonators such those obtained with semiconductor microcavities or superconducting quantum circuits. Our study has shown that a critical slowing down of the dynamics occurs for some parameters in dimensions higher than one. The phenomena is linked to a first-order dissipative phase transition and to lattice-level collective bistability. Furthermore, we reveal the role of disorder on the spatial correlation properties of these phenomena. In the last part of the thesis, we present the strong connection between fundamental concepts of machine learning and variational approaches in quantum physics. We present a novel computational method for open quantum systems based on neural network ans\"{a}tze that is compatible with a stochastic minimisation of a cost function through Monte Carlo sampling. As a first application, we have applied this method to the dissipative transverse quantum Ising model. Exciting perspectives and possible developments are also drawn in the concluding remarks.
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 introduce a method to simulate open quantum many-body dynamics by combining time-dependent variational Monte Carlo (tVMC) with quantum trajectory techniques. Our approach unravels the Lindblad master equation into an ensemble of stochastic Schrödinger equations for a variational ansatz, avoiding the exponential cost of density matrix evolution. The method is compatible with generic ansatzë, including expressive neural-network wavefunctions. We derive the nonlinear stochastic equations of motion for the variational parameters and employ suitable Stratonovich numerical solvers. To validate our approach, we simulate quenches in the locally dissipative long-range Ising model in a transverse field, accurately capturing non-equilibrium magnetization and spin squeezing dynamics relevant to trapped-ion and Rydberg atom experiments. The framework is computationally efficient, scalable on high-performance computing platforms, and can be readily integrated into existing tVMC implementations. This work enables large-scale simulations of complex, dissipative quantum systems in higher dimensions, with broad implications for quantum technology and fundamental science.
Many-body open quantum systems are exposed to an essentially uncontrollable environment that acts as a source of decoherence and dissipation. As the exact treatment of such models is generally unfeasible, it is favourable to formulate an approximate description by means of the dynamics of the reduced density operator of the system. When studying models with a weak coupling to a memoryless environment, the existence of such solution is granted and the reduced dynamics is governed by the Lindblad quantum master equation. In recent years, open quantum systems have evolved into a major field of studies. Focus of these studies are the characterization of emergent phenomena and dissipative phase transitions, as well as the ongoing debate about whether quantum computing schemes are still hard to simulate classically --- and thus achieve quantum supremacy --- when in presence of some degree of noise-induced decoherence. However, the challenge posed by this task lies in the complexity of the density matrix that increases exponentially with the system size, and the quest for efficient numerical methods is a research field that is still in its infancy. Here, we first develop a real-time full configuration interaction quantum Monte Carlo technique that stems from a class of methods generally known as projector Monte Carlo. The approach enables the stochastic sampling of the Lindblad time evolution of the density matrix thanks to a massively parallel algorithm, thus providing estimates of observables on the non-equilibrium steady state. We present the underlying theory and introduce an initiator technique and importance sampling to reduce statistical error. We demonstrate the efficiency of the approach by applying it to the dissipative two-dimensional XYZ spin-1/2 model on a lattice. As the importance sampling of projector approaches is often combined with variational results, we then introduce a novel method that is based on the 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 Lindblad quantum master equation. We test the effectiveness of the method by modeling the steady state of the dissipative two-dimensional XYZ spin-1/2 model through a dissipative phase transition, and also the real-time dynamics of the dissipative Ising model. In addition, we discuss the application of the developed methods and the open questions of the field.