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:2105.01657 · Quantum

QuantumCumulants.jl: A Julia framework for generalized mean-field equations in open quantum systems

David Plankensteiner, Christoph Hotter, Helmut Ritsch

Published Jan 4, 2022• 84 citations
arXivHTMLPDF

Abstract

A full quantum mechanical treatment of open quantum systems via a Master equation is often limited by the size of the underlying Hilbert space. As an alternative, the dynamics can also be formulated in terms of systems of coupled differential equations for operators in the Heisenberg picture. This typically leads to an infinite hierarchy of equations for products of operators. A well-established approach to truncate this infinite set at the level of expectation values is to neglect quantum correlations of high order. This is systematically realized with a so-called cumulant expansion, which decomposes expectation values of operator products into products of a given lower order, leading to a closed set of equations. Here we present an open-source framework that fully automizes this approach: first, the equations of motion of operators up to a desired order are derived symbolically using predefined canonical commutation relations. Next, the resulting equations for the expectation values are expanded employing the cumulant expansion approach, where moments up to a chosen order specified by the user are included. Finally, a numerical solution can be directly obtained from the symbolic equations. After reviewing the theory we present the framework and showcase its usefulness in a few example problems.

DOI
10.22331/q-2022-01-04-617

Related papers

Don't truncate, decompose: mean-field dynamics of long-range quantum systems from strongly correlated states

Federico Carollo

We challenge the widespread consensus that mean-field theory fails to describe long-range open quantum systems in the presence of symmetry breaking and/or when starting from strongly correlated states (e.g., macroscopic superpositions). While recent literature relies on cumulant expansions to capture such systems, this approach rests on truncations with no clear justification. Here, we show that it is, at best, conceptually redundant in the strong long-range regime. We show that the evolution can be decomposed into, and fully reconstructed from, independent mean-field dynamics. This decomposition generates the entire hierarchy of cumulants and, as a byproduct, identifies---to our knowledge, for the first time---a regime in which cumulant expansions exactly predict low-order cumulants. We illustrate the power of our findings with two applications: we compute the moment generating function for nonequilibrium $\mathcal{Z}_2$ symmetry breaking, and construct states restoring time-translation symmetry in time crystals. In both cases, our method fully reproduces the exact many-body dynamics, which is out of reach of cumulant expansions. Our results reclaim the exactness of mean-field theory, offering a transparent framework for large-scale open quantum systems.

cond-mat.stat-mechCross-listedarXiv:2607.25434•Jul 29, 2026•0 citations

Determination of Quantum Expectation Values Via Fluctuation Expansion

Metin Demіralp

The goal of this work is, the utilization of a method recently developed by the author, in quantum dynamical problems. Main idea is based on the relation between quantum and classical dynamics such that quantum dynamical problems become their classical dynamical counterparts when they are formulated via expectation values of certain operators and the probability density tends to be sharply localized. Expectation values can be expanded into a series in ascending appearence multiplicity of the complement of the projection operator which projects the space spanned by the wave function. This expansion contains fluctuation terms as unknowns besides the expectation values. We construct an infinite set of differential equations. The finite truncations of this set enables us to approximate the quantum dynamics of the system under consideration.

CRC Press eBooks•May 4, 2022•2 citations

Universal structure of computing moments for exact quantum dynamics: Application to arbitrary system–bath couplings

Ruihao Bi, Wei Liu, Wenjie Dou

We introduce a general procedure for computing higher-order moments of correlation functions in open quantum systems, extending the scope of our recent work on Memory Kernel Coupling Theory (MKCT) [Liu et al., arXiv:2407.01923 (2024)]. This approach is demonstrated for arbitrary system-bath coupling that can be expressed as a polynomial, HSB=V̂(α0+α1q̂+α2q̂2+⋯) , where we show that the recursive commutators of a system operator obey a universal hierarchy. Exploiting this structure, the higher-order moments are obtained by evaluating the expectation values of the system and bath operators separately, with bath expectation values derived from the derivatives of a generating function. We further apply MKCT to compute the dipole autocorrelation function for the spin-boson model with both linear and quadratic coupling, achieving agreement with the hierarchical equations of motion approach. Our findings suggest a promising path toward accurate dynamics for complex open quantum systems.

arXiv:2504.00649•The Journal of Chemical Physics•Jun 11, 2025•6 citations

Cumulant expansion for studying damped quantum solitons

Eduard Schmidt, L. Knöll, Dirk‐Gunnar Welsch

The quantum statistics of damped optical solitons is studied using cumulant-expansion techniques. The effect of absorption is described in terms of ordinary Markovian relaxation theory, by coupling the optical field to a continuum of reservoir modes. After introduction of local bosonic field operators and spatial discretization, pseudo-Fokker-Planck equations for multidimensional s-parametrized phase-space functions are derived. These partial differential equations are equivalent to an infinite set of ordinary differential equations for the cumulants of the phase-space functions. Introducing an appropriate truncation condition, the resulting finite set of cumulant evolution equations can be solved numerically. Solutions are presented in a Gaussian approximation and the quantum noise is calculated, with special emphasis on squeezing and the recently measured spectral photon-number correlations [Sp\"alter et al., Phys. Rev. Lett. 81, 786 (1998)].

arXiv:quant-ph/9810018•Physical Review A•Mar 1, 1999•16 citations

A perturbative probabilistic approach to quantum many-body systems

Andrea Di Stefano, Massimo Ostilli, Carlo Presilla A

In the probabilistic approach to quantum many-body systems, the ground-state energy is the solution of a nonlinear scalar equation written either as a cumulant expansion or as an expectation with respect to a probability distribution of the potential and hopping (amplitude and phase) values recorded during an infinitely lengthy evolution. We introduce a perturbative expansion of this probability distribution which conserves, at any order, a multinomial-like structure, typical of uncorrelated systems, but includes, order by order, the statistical correlations provided by the cumulant expansion. The proposed perturbative scheme is successfully tested in the case of pseudo spin 1/2 hard-core boson Hubbard models also when affected by a phase problem due to an applied magnetic field.

arXiv:1210.3976•Journal of Statistical Mechanics Theory and Experiment•Apr 2, 2013•1 citation

The Cumulant Expansion Approach: The Good, The Bad and The Ugly

Johannes Kerber, Helmut Ritsch, Laurin Ostermann

The configuration space, i.e. the Hilbert space, of compound quantum systems grows exponentially with the number of its subsystems: its dimensionality is given by the product of the dimensions of its constituents. Therefore a full quantum treatment is rarely possible analytically and can be carried out numerically for fairly small systems only. Fortunately, in order to obtain interesting physics, approximations often very well suffice. One of these approximations is given by the cumulants expansion, where expectation values of products of operators are approximated by products of expectation values of said operators, neglecting higher-order correlations. The lowest order of this approximation is widely known as the mean field approximation and used routinely throughout quantum physics. Despite its ubiquitous presence, a general criterion for applicability and convergence properties of higher order cumulants expansions remains to be found. In this paper, we discuss two problems in quantum electrodynamics and quantum information, namely the collective radiative dissipation of a dipole-dipole interacting chain of atoms and the factorization of a bi-prime by annealing in an adiabatic quantum simulator. In the first case we find smooth, convergence behavior, where the approximation performs increasingly better with higher orders, while in the latter going beyond mean field turns out useless and, even for small system sizes, we are puzzled by numerically challenging and partly non-physical solutions.

quant-pharXiv:2511.20115•Nov 26, 2025•0 citations
A full quantum mechanical treatment of open quantum systems via a Master equation is often limited by the size of the underlying Hilbert space. As an alternative, the dynamics can also be formulated in terms of systems of coupled differential equations for operators in the Heisenberg picture. This typically leads to an infinite hierarchy of equations for products of operators. A well-established approach to truncate this infinite set at the level of expectation values is to neglect quantum correlations of high order. This is systematically realized with a so-called cumulant expansion, which decomposes expectation values of operator products into products of a given lower order, leading to a closed set of equations. Here we present an open-source framework that fully automizes this approach: first, the equations of motion of operators up to a desired order are derived symbolically using predefined canonical commutation relations. Next, the resulting equations for the expectation values are expanded employing the cumulant expansion approach, where moments up to a chosen order specified by the user are included. Finally, a numerical solution can be directly obtained from the symbolic equations. After reviewing the theory we present the framework and showcase its usefulness in a few example problems.
We challenge the widespread consensus that mean-field theory fails to describe long-range open quantum systems in the presence of symmetry breaking and/or when starting from strongly correlated states (e.g., macroscopic superpositions). While recent literature relies on cumulant expansions to capture such systems, this approach rests on truncations with no clear justification. Here, we show that it is, at best, conceptually redundant in the strong long-range regime. We show that the evolution can be decomposed into, and fully reconstructed from, independent mean-field dynamics. This decomposition generates the entire hierarchy of cumulants and, as a byproduct, identifies---to our knowledge, for the first time---a regime in which cumulant expansions exactly predict low-order cumulants. We illustrate the power of our findings with two applications: we compute the moment generating function for nonequilibrium Z2\mathcal{Z}_2Z2​ symmetry breaking, and construct states restoring time-translation symmetry in time crystals. In both cases, our method fully reproduces the exact many-body dynamics, which is out of reach of cumulant expansions. Our results reclaim the exactness of mean-field theory, offering a transparent framework for large-scale open quantum systems.
The goal of this work is, the utilization of a method recently developed by the author, in quantum dynamical problems. Main idea is based on the relation between quantum and classical dynamics such that quantum dynamical problems become their classical dynamical counterparts when they are formulated via expectation values of certain operators and the probability density tends to be sharply localized. Expectation values can be expanded into a series in ascending appearence multiplicity of the complement of the projection operator which projects the space spanned by the wave function. This expansion contains fluctuation terms as unknowns besides the expectation values. We construct an infinite set of differential equations. The finite truncations of this set enables us to approximate the quantum dynamics of the system under consideration.
We introduce a general procedure for computing higher-order moments of correlation functions in open quantum systems, extending the scope of our recent work on Memory Kernel Coupling Theory (MKCT) [Liu et al., arXiv:2407.01923 (2024)]. This approach is demonstrated for arbitrary system-bath coupling that can be expressed as a polynomial, HSB=V̂(α0+α1q̂+α2q̂2+⋯) , where we show that the recursive commutators of a system operator obey a universal hierarchy. Exploiting this structure, the higher-order moments are obtained by evaluating the expectation values of the system and bath operators separately, with bath expectation values derived from the derivatives of a generating function. We further apply MKCT to compute the dipole autocorrelation function for the spin-boson model with both linear and quadratic coupling, achieving agreement with the hierarchical equations of motion approach. Our findings suggest a promising path toward accurate dynamics for complex open quantum systems.
The quantum statistics of damped optical solitons is studied using cumulant-expansion techniques. The effect of absorption is described in terms of ordinary Markovian relaxation theory, by coupling the optical field to a continuum of reservoir modes. After introduction of local bosonic field operators and spatial discretization, pseudo-Fokker-Planck equations for multidimensional s-parametrized phase-space functions are derived. These partial differential equations are equivalent to an infinite set of ordinary differential equations for the cumulants of the phase-space functions. Introducing an appropriate truncation condition, the resulting finite set of cumulant evolution equations can be solved numerically. Solutions are presented in a Gaussian approximation and the quantum noise is calculated, with special emphasis on squeezing and the recently measured spectral photon-number correlations [Sp\"alter et al., Phys. Rev. Lett. 81, 786 (1998)].
In the probabilistic approach to quantum many-body systems, the ground-state energy is the solution of a nonlinear scalar equation written either as a cumulant expansion or as an expectation with respect to a probability distribution of the potential and hopping (amplitude and phase) values recorded during an infinitely lengthy evolution. We introduce a perturbative expansion of this probability distribution which conserves, at any order, a multinomial-like structure, typical of uncorrelated systems, but includes, order by order, the statistical correlations provided by the cumulant expansion. The proposed perturbative scheme is successfully tested in the case of pseudo spin 1/2 hard-core boson Hubbard models also when affected by a phase problem due to an applied magnetic field.
The configuration space, i.e. the Hilbert space, of compound quantum systems grows exponentially with the number of its subsystems: its dimensionality is given by the product of the dimensions of its constituents. Therefore a full quantum treatment is rarely possible analytically and can be carried out numerically for fairly small systems only. Fortunately, in order to obtain interesting physics, approximations often very well suffice. One of these approximations is given by the cumulants expansion, where expectation values of products of operators are approximated by products of expectation values of said operators, neglecting higher-order correlations. The lowest order of this approximation is widely known as the mean field approximation and used routinely throughout quantum physics. Despite its ubiquitous presence, a general criterion for applicability and convergence properties of higher order cumulants expansions remains to be found. In this paper, we discuss two problems in quantum electrodynamics and quantum information, namely the collective radiative dissipation of a dipole-dipole interacting chain of atoms and the factorization of a bi-prime by annealing in an adiabatic quantum simulator. In the first case we find smooth, convergence behavior, where the approximation performs increasingly better with higher orders, while in the latter going beyond mean field turns out useless and, even for small system sizes, we are puzzled by numerically challenging and partly non-physical solutions.