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

Modelling many-body quantum dynamics with stochastic trajectories: a critical test on the Tavis-Cummings model

A. Leonau, S. Chuchurka, V. Sukharnikov, A. Benediktovitch, N. Rohringer

quant-phSubmitted Nov 11, 2025• Announced Nov 13, 2025• 0 citations
arXivHTMLPDF

Abstract

We critically explore the applicability of a recently proposed framework to sample the quantum dynamics of a many-body quantum system interacting with light by stochastic trajectories, applying it to the closed and open Tavis-Cummings model (TCM). The stochastic differential equations (SDEs) sample the positive P phase-space representation by analog complex-valued dynamical variables that are linked to the quantum operators. Statistical average over the stochastic trajectories yields the evolution of the quantum mechanical expectation values. However, numerical implementation of these SDEs for the TCM indicates divergent solutions, also known from other phase-space methods. This limits the applicability of the framework to finite propagation times, that are strongly dependent on the physical parameters and initial conditions of the system. We outline the underlying mathematical reason for these divergences and show that their contribution to the averages are, however, essential. To attempt to regularize the divergences, we transform the SDEs to an equivalent set of SDEs with different noise realisations, thereby pushing the valid time boundary. Quantum collapse and revival of the TCM, however, cannot be recovered by the stochastic trajectory approach, pointing to the general difficulty of the applicability of stochastic phase-space sampling methods to systems with strong quantum features.

Related papers

Path-integral formulation of bosonic Markovian open quantum dynamics with Monte Carlo stochastic trajectories using the Glauber–Sudarshan P, Wigner, and Husimi Q functions and hybrids

Toma Yoneya, Kazuya Fujimoto, Yuki Kawaguchi

The Monte Carlo trajectory sampling of stochastic differential equations based on the quasiprobability distribution functions, such as the Glauber–Sudarshan P, Wigner, and Husimi Q functions, enables us to investigate bosonic open quantum many-body dynamics described by the Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) equation. In this method, the Monte Carlo samplings for the initial distribution and stochastic noises incorporate quantum fluctuations, and thus, we can go beyond the mean-field approximation. However, description using stochastic differential equations is possible only when the corresponding Fokker–Planck equation has a positive-semidefinite diffusion matrix. In this work, we analytically derive the stochastic differential equations for arbitrary Hamiltonian and jump operators based on the path-integral formula, independently of the derivation of the Fokker–Planck equation. In the course of the derivation, we formulate the path-integral representation of the GKSL equation by using the -ordered quasiprobability distribution function, which systematically describes the aforementioned quasiprobability distribution functions by changing the real parameter . The essential point of this derivation is that we employ the Hubbard–Stratonovich transformation in the path integral, and its application is not always feasible. We find that the feasible condition of the Hubbard–Stratonovich transformation is identical to the positive-semidefiniteness condition of the diffusion matrix in the Fokker–Planck equation. In the benchmark calculations, we confirm that the Monte Carlo simulations of the obtained stochastic differential equations well reproduce the exact dynamics of physical quantities and non-equal time correlation functions of numerically solvable models, including the Bose–Hubbard model. This work clarifies the applicability of the approximation and gives systematic and simplified procedures to obtain the stochastic differential equations to be numerically solved.

arXiv:2508.01991•Annals of Physics•Dec 16, 2025•0 citations

Universal Dilation of Linear Itô SDEs: Quantum Trajectories and Lindblad Simulation of Second Moments

Hsuan-Cheng Wu, Xiantao Li

We present a universal framework for simulating $N$-dimensional linear Itô stochastic differential equations (SDEs) on quantum computers with additive or multiplicative noises. Building on a unitary dilation technique, we establish a rigorous mapping from the general linear SDEs \[ dX_t = A(t) X_t\,dt + \sum_{j=1}^J B_j(t)X_t\,dW_t^j \] to stochastic Schrödinger equations (SSE) on a dilated Hilbert space. Crucially, this embedding is pathwise exact in that the classical solution is recovered as a projection of the dilated quantum state for each fixed noise realization. We demonstrate that the resulting SSEs are {naturally implementable} on digital quantum processors, where the stochastic Wiener increments are encoded directly by preparing the ancillary qubits. Exploiting this physical mapping, we develop two algorithmic strategies: (1) a trajectory-based approach that uses sequential weak measurements to realize efficient stochastic integrators, including a second-order scheme, and (2) an ensemble-based approach that maps moment evolution to a deterministic Lindblad quantum master equation, enabling simulation without Monte Carlo sampling. We provide error bounds based on a stochastic light-cone analysis and validate the framework with numerical experiments.

quant-pharXiv:2601.05928•Quantum•Jan 12, 2026•0 citations

The Dynamics of the Hubbard Model Through Stochastic Calculus and Girsanov Transformation

Detlef Lehmann

Abstract As a typical quantum many body problem, we consider the time evolution of density matrix elements in the Bose-Hubbard model. For an arbitrary initial state, these quantities can be obtained from an SDE or stochastic differential equation system. To this SDE system, a Girsanov transformation can be applied. This has the effect that all the information from the initial state moves into the drift part, into the mean field part, of the transformed system. In the large N limit with $$g=UN$$ g = U N fixed, the diffusive part of the transformed system vanishes and as a result, the exact quantum dynamics is given by an ODE system which turns out to be the time dependent discrete Gross Pitaevskii equation. For the two site Bose-Hubbard model, the GP equation reduces to the mathematical pendulum and the difference of expected number of particles at the two lattice sites is equal to the velocity of that pendulum which is either oscillatory or it can have rollovers which then corresponds to the self trapping or insulating phase. As a by-product, we also find an equivalence of the mathematical pendulum with a quartic double well potential. Collapse and revivals are a more subtle phenomenom, in order to see these the diffusive part of the SDE system or quantum corrections have to be taken into account. This can be done with an approximation and collapse and revivals can be reproduced, numerically and also through an analytic calculation. Since expectation values of Fresnel or Wiener diffusion processes, we write the density matrix elements exactly in this way, can be obtained from parabolic second order PDEs, we also obtain various exact PDE representations. The paper has been written with the goal to come up with an efficient calculation scheme for quantum many body systems and as such the formalism is generic and applies to arbitrary dimension, arbitrary hopping matrices and, with suitable adjustments, to fermionic models.

arXiv:2205.02010•International Journal of Theoretical Physics•May 28, 2024•1 citation

Criteria for Feasible Monte Carlo Stochastic Simulations of Bosonic Markovian Open Quantum Dynamics

Toma Yoneya, Kazuya Fujimoto, Yuki Kawaguchi

The Monte Carlo sampling of the stochastic differential equations (SDEs) based on the quasiprobability distribution function, such as the Glauber--Sudarshan P, Wigner, and Husimi Q functions provides a powerful framework for investigating bosonic open quantum many-body dynamics described by the Gorini--Kossakowski--Sudarshan--Lindblad (GKSL) equation, while considering the effects of quantum fluctuations beyond the mean-field approximation. However, the stochastic Monte Carlo simulation is possible only when the corresponding Fokker--Planck equation has a positive-semidefinite diffusion matrix, and the general conditions for the diffusion matrix to be positive semidefinite have remained unclear. In this work, starting from the path integral formulation, we first derive the sufficient conditions under which the diffusion matrix is positive semidefinite for an arbitrary Hamiltonian, jump operators, and choice of quasiprobability distribution functions. We also analytically derive the corresponding SDEs to be solved. We then investigate the dynamics of the GKSL equation in the thermodynamic limit and show that, depending on the form of the jump operators, the mean-field approximation may fail to describe the dynamics accurately, making stochastic Monte Carlo simulations indispensable. Furthermore, we derive the sufficient conditions under which the higher-order quantum fluctuation terms beyond the Fokker--Planck description vanish identically, even when the jump operators contain quadratic terms. Under these conditions, whenever the corresponding SDEs can be derived, the stochastic Monte Carlo simulation reproduces the exact dynamics. These results clarify the conditions under which the stochastic Monte Carlo simulations are both feasible and necessary for accurately describing the dynamics governed by the GKSL equation in phase space.

cond-mat.quant-gasCross-listedarXiv:2608.06056•Annals of Physics•Aug 7, 2026•0 citations

Can Quantum Field Theory be Recovered from Time-Symmetric Stochastic Mechanics? Part II: Prospects for a Trajectory Interpretation

Simon Friederich, Mritunjay Tyagi

In a companion paper we derived a unique time-reversal-invariant stochastic generalization of the Liouville equation and showed that it coincides with the evolution equation for the Husimi $Q$-function in a broad class of bosonic quantum field theories. Here we investigate the prospects for interpreting that evolution equation in terms of underlying stochastic trajectories. Drawing on Drummond's time-symmetric stochastic action formalism, we show that the traceless diffusion Fokker-Planck equation defines a natural measure over stochastic trajectories conditional on mixed-time boundary conditions. However, we identify a significant gap: it has not been established that every $Q$-function can be represented as a weighted average of these conditional probabilities over boundary values. The trajectory interpretation holds for ensembles with fixed boundary conditions but does not straightforwardly extend to arbitrary quantum states. Despite this limitation, we show that Drummond's trajectory dynamics are fundamentally non-Markovian -- a natural consequence of combining stochasticity with time-reversal invariance. This non-Markovianity places the dynamics outside the scope of the ontological models framework and thereby explains why the major no-go theorems for hidden-variable theories do not rule out the approach. These results clarify both the achievements and the remaining challenges in the project of understanding quantum field theory as the statistical mechanics of time-symmetric stochastic processes.

quant-pharXiv:2603.28983•Apr 1, 2026•0 citations

Quantum stochastic trajectories for particles and fields based on positive P-representation

Stasis Chuchurka, Andrei Benediktovitch, Nina Rohringer

In this work we introduce a phase-space description based on the positive P representation for bosonic fields interacting with a system of quantum emitters. The formalism is applicable to collective light-matter interactions and open quantum systems with decoherence. Conservation of particle numbers is considered, and a Jordan-Schwinger transformation enables the representation of multi-level quantum emitters. The evolution of the phase-space description of the combined system of emitters and field is formulated in terms of stochastic trajectories and we derive the rules of mapping from traditional quantum mechanics to this stochastic formalism. The resulting equations of motion encode deterministic, classical evolution with quantum effects incorporated by stochastic noise terms. The framework's equations and properties are provided without specifying the Hamiltonian, aiming for broad applicability in diverse research domains. A potential future application is the quantum mechanical description of collective spontaneous emission of an incoherently pumped ensemble of atoms.

arXiv:2306.17474•Jun 30, 2023•0 citations
We critically explore the applicability of a recently proposed framework to sample the quantum dynamics of a many-body quantum system interacting with light by stochastic trajectories, applying it to the closed and open Tavis-Cummings model (TCM). The stochastic differential equations (SDEs) sample the positive P phase-space representation by analog complex-valued dynamical variables that are linked to the quantum operators. Statistical average over the stochastic trajectories yields the evolution of the quantum mechanical expectation values. However, numerical implementation of these SDEs for the TCM indicates divergent solutions, also known from other phase-space methods. This limits the applicability of the framework to finite propagation times, that are strongly dependent on the physical parameters and initial conditions of the system. We outline the underlying mathematical reason for these divergences and show that their contribution to the averages are, however, essential. To attempt to regularize the divergences, we transform the SDEs to an equivalent set of SDEs with different noise realisations, thereby pushing the valid time boundary. Quantum collapse and revival of the TCM, however, cannot be recovered by the stochastic trajectory approach, pointing to the general difficulty of the applicability of stochastic phase-space sampling methods to systems with strong quantum features.
The Monte Carlo trajectory sampling of stochastic differential equations based on the quasiprobability distribution functions, such as the Glauber–Sudarshan P, Wigner, and Husimi Q functions, enables us to investigate bosonic open quantum many-body dynamics described by the Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) equation. In this method, the Monte Carlo samplings for the initial distribution and stochastic noises incorporate quantum fluctuations, and thus, we can go beyond the mean-field approximation. However, description using stochastic differential equations is possible only when the corresponding Fokker–Planck equation has a positive-semidefinite diffusion matrix. In this work, we analytically derive the stochastic differential equations for arbitrary Hamiltonian and jump operators based on the path-integral formula, independently of the derivation of the Fokker–Planck equation. In the course of the derivation, we formulate the path-integral representation of the GKSL equation by using the -ordered quasiprobability distribution function, which systematically describes the aforementioned quasiprobability distribution functions by changing the real parameter . The essential point of this derivation is that we employ the Hubbard–Stratonovich transformation in the path integral, and its application is not always feasible. We find that the feasible condition of the Hubbard–Stratonovich transformation is identical to the positive-semidefiniteness condition of the diffusion matrix in the Fokker–Planck equation. In the benchmark calculations, we confirm that the Monte Carlo simulations of the obtained stochastic differential equations well reproduce the exact dynamics of physical quantities and non-equal time correlation functions of numerically solvable models, including the Bose–Hubbard model. This work clarifies the applicability of the approximation and gives systematic and simplified procedures to obtain the stochastic differential equations to be numerically solved.
We present a universal framework for simulating NNN-dimensional linear Itô stochastic differential equations (SDEs) on quantum computers with additive or multiplicative noises. Building on a unitary dilation technique, we establish a rigorous mapping from the general linear SDEs dXt=A(t)Xt dt+∑j=1JBj(t)Xt dWtj dX_t = A(t) X_t\,dt + \sum_{j=1}^J B_j(t)X_t\,dW_t^j dXt​=A(t)Xt​dt+∑j=1J​Bj​(t)Xt​dWtj​ to stochastic Schrödinger equations (SSE) on a dilated Hilbert space. Crucially, this embedding is pathwise exact in that the classical solution is recovered as a projection of the dilated quantum state for each fixed noise realization. We demonstrate that the resulting SSEs are {naturally implementable} on digital quantum processors, where the stochastic Wiener increments are encoded directly by preparing the ancillary qubits. Exploiting this physical mapping, we develop two algorithmic strategies: (1) a trajectory-based approach that uses sequential weak measurements to realize efficient stochastic integrators, including a second-order scheme, and (2) an ensemble-based approach that maps moment evolution to a deterministic Lindblad quantum master equation, enabling simulation without Monte Carlo sampling. We provide error bounds based on a stochastic light-cone analysis and validate the framework with numerical experiments.
Abstract As a typical quantum many body problem, we consider the time evolution of density matrix elements in the Bose-Hubbard model. For an arbitrary initial state, these quantities can be obtained from an SDE or stochastic differential equation system. To this SDE system, a Girsanov transformation can be applied. This has the effect that all the information from the initial state moves into the drift part, into the mean field part, of the transformed system. In the large N limit with g=UNg=UNg=UN g = U N fixed, the diffusive part of the transformed system vanishes and as a result, the exact quantum dynamics is given by an ODE system which turns out to be the time dependent discrete Gross Pitaevskii equation. For the two site Bose-Hubbard model, the GP equation reduces to the mathematical pendulum and the difference of expected number of particles at the two lattice sites is equal to the velocity of that pendulum which is either oscillatory or it can have rollovers which then corresponds to the self trapping or insulating phase. As a by-product, we also find an equivalence of the mathematical pendulum with a quartic double well potential. Collapse and revivals are a more subtle phenomenom, in order to see these the diffusive part of the SDE system or quantum corrections have to be taken into account. This can be done with an approximation and collapse and revivals can be reproduced, numerically and also through an analytic calculation. Since expectation values of Fresnel or Wiener diffusion processes, we write the density matrix elements exactly in this way, can be obtained from parabolic second order PDEs, we also obtain various exact PDE representations. The paper has been written with the goal to come up with an efficient calculation scheme for quantum many body systems and as such the formalism is generic and applies to arbitrary dimension, arbitrary hopping matrices and, with suitable adjustments, to fermionic models.
The Monte Carlo sampling of the stochastic differential equations (SDEs) based on the quasiprobability distribution function, such as the Glauber--Sudarshan P, Wigner, and Husimi Q functions provides a powerful framework for investigating bosonic open quantum many-body dynamics described by the Gorini--Kossakowski--Sudarshan--Lindblad (GKSL) equation, while considering the effects of quantum fluctuations beyond the mean-field approximation. However, the stochastic Monte Carlo simulation is possible only when the corresponding Fokker--Planck equation has a positive-semidefinite diffusion matrix, and the general conditions for the diffusion matrix to be positive semidefinite have remained unclear. In this work, starting from the path integral formulation, we first derive the sufficient conditions under which the diffusion matrix is positive semidefinite for an arbitrary Hamiltonian, jump operators, and choice of quasiprobability distribution functions. We also analytically derive the corresponding SDEs to be solved. We then investigate the dynamics of the GKSL equation in the thermodynamic limit and show that, depending on the form of the jump operators, the mean-field approximation may fail to describe the dynamics accurately, making stochastic Monte Carlo simulations indispensable. Furthermore, we derive the sufficient conditions under which the higher-order quantum fluctuation terms beyond the Fokker--Planck description vanish identically, even when the jump operators contain quadratic terms. Under these conditions, whenever the corresponding SDEs can be derived, the stochastic Monte Carlo simulation reproduces the exact dynamics. These results clarify the conditions under which the stochastic Monte Carlo simulations are both feasible and necessary for accurately describing the dynamics governed by the GKSL equation in phase space.
In a companion paper we derived a unique time-reversal-invariant stochastic generalization of the Liouville equation and showed that it coincides with the evolution equation for the Husimi QQQ-function in a broad class of bosonic quantum field theories. Here we investigate the prospects for interpreting that evolution equation in terms of underlying stochastic trajectories. Drawing on Drummond's time-symmetric stochastic action formalism, we show that the traceless diffusion Fokker-Planck equation defines a natural measure over stochastic trajectories conditional on mixed-time boundary conditions. However, we identify a significant gap: it has not been established that every QQQ-function can be represented as a weighted average of these conditional probabilities over boundary values. The trajectory interpretation holds for ensembles with fixed boundary conditions but does not straightforwardly extend to arbitrary quantum states. Despite this limitation, we show that Drummond's trajectory dynamics are fundamentally non-Markovian -- a natural consequence of combining stochasticity with time-reversal invariance. This non-Markovianity places the dynamics outside the scope of the ontological models framework and thereby explains why the major no-go theorems for hidden-variable theories do not rule out the approach. These results clarify both the achievements and the remaining challenges in the project of understanding quantum field theory as the statistical mechanics of time-symmetric stochastic processes.
In this work we introduce a phase-space description based on the positive P representation for bosonic fields interacting with a system of quantum emitters. The formalism is applicable to collective light-matter interactions and open quantum systems with decoherence. Conservation of particle numbers is considered, and a Jordan-Schwinger transformation enables the representation of multi-level quantum emitters. The evolution of the phase-space description of the combined system of emitters and field is formulated in terms of stochastic trajectories and we derive the rules of mapping from traditional quantum mechanics to this stochastic formalism. The resulting equations of motion encode deterministic, classical evolution with quantum effects incorporated by stochastic noise terms. The framework's equations and properties are provided without specifying the Hamiltonian, aiming for broad applicability in diverse research domains. A potential future application is the quantum mechanical description of collective spontaneous emission of an incoherently pumped ensemble of atoms.