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.20115

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

Johannes Kerber, Helmut Ritsch, Laurin Ostermann

quant-phphysics.opticsSubmitted Nov 25, 2025• Announced Nov 26, 2025• 0 citations• Revised Sep 4, 2026 (v2)
v1 · Nov 25, 2025v2 · Sep 4, 2026
arXivHTMLPDF

Abstract

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.

Comments
10 pages, 10 figures and 3 tables. Minor revisions and formatting updates throughout

Related papers

Light-scattering properties beyond weak-field excitation in atomic ensembles

Chung-Hsien Wang, Nai-Yu Tsai, Yi-Cheng Wang +1 more

In the study of optical properties of large atomic systems, a weak laser driving is often assumed to simplify the system dynamics by linearly coupled equations. Here we investigate the light-scattering properties of atomic ensembles beyond weak-field excitation through the cumulant expansion method. By progressively incorporating higher-order correlations into the steady-state equations, an enhanced accuracy can be achieved in comparison to the exact solutions obtained by solving a full density matrix. Our analysis reveals that, in the regime of weak dipole-dipole interaction, the first-order expansion yields satisfactory predictions for optical depth, while denser atomic configurations necessitate consideration of higher-order correlations. As the intensity of incident light increases, atom saturation effects become noticeable, giving rise to significant changes in light transparency, energy shift, and decay rate. This saturation phenomenon extends to subradiant atom arrays even under weak driving conditions, leading to substantial deviations from the linear model. Our findings demonstrate that the mean-field model is a good extension to linear models as it balances both accuracy and computational complexity. However, the crucial role of higher-order cumulants in large and dense atom systems remains unclear, since it is challenging theoretically owing to the exponentially increasing Hilbert space in such light-matter interacting systems.

arXiv:2310.17106•Physical Review A•Jul 22, 2024•3 citations

Cumulant expansion approach to nonlinear dynamics of inhomogeneous ensembles in cavity QED

Matthias Zens

In the last couple of years much effort has been dedicated to the development and the investigation of hybrid quantum devices, which combine the advantages of very dissimilar quantum systems for the realization of efficient quantum computation and communication technologies. Particularly interesting in this regard are solid-state quantum memories based on ensembles of spins (e.g. nitrogen vacancy defects in diamond or rare-earth doped crystals) coupled to superconducting microwave cavities. Here the spin-ensemble acts as a robust memory where the collective coupling to the cavity mode allows for the coherent transfer of quantum information. The major downside of ensembles inside solid-state systems is their natural tendency to exhibit inhomogeneous broadening of the transition frequencies, which makes experiments as well as their theoretical description very challenging. On the theoretical side it is primarily the computational complexity due to the exponential growth of the associated Hilbert space with the system size that inevitably calls for approximation schemes. To date, the theoretical description of large spin-cavity systems is mainly limited to weak driving fields, where the number of spin-excitations in the ensemble is negligible, corresponding to a Holstein-Primakoff approximation, or to mean-field approaches, where correlations within the system are neglected. In this thesis we aim to develop a model that accurately accounts for the dynamics of large inhomogeneously broadened spin ensembles, coupled to a single cavity mode, even in the case of strong driving fields. In particular we employ the generalized cumulant expansion method to account for correlations within the system, which allows us to go beyond the Holstein-Primakoff and mean-field approximations. As a first step, we demonstrate the applicability of the cumulant expansion technique to the well-studied Jaynes-Cummings model and investigate its main limitations. We then move on to the Tavis-Cummings model with multiple spins inside the cavity and generalize our approach to ensembles containing a very large number of spins (ca. 10 12). Our model aims not only to provide a proper description of ensemble-based quantum memories, but also to serve as a novel tool for the investigation of interesting new physics arising from cooperative effects in inhomogeneous ensembles.

•Sep 1, 2016•1 citation

Determining the validity of cumulant expansions for central spin models

Piper Fowler-Wright, Kristín B. Arnardóttir, Peter Kirton +2 more

For a model with many-to-one connectivity it is widely expected that mean-field theory captures the exact many-particle $N\ensuremath{\rightarrow}\ensuremath{\infty}$ limit, and that higher-order cumulant expansions of the Heisenberg equations converge to this same limit whilst providing improved approximations at finite $N$. Here we show that this is in fact not always the case. Instead, whether mean-field theory correctly describes the large-$N$ limit depends on how the model parameters scale with $N$, and the convergence of cumulant expansions may be nonuniform across even and odd orders. Further, even when a higher-order cumulant expansion does recover the correct limit, the error is not monotonic with $N$ and may exceed that of mean-field theory.

arXiv:2303.04410•Physical Review Research•Sep 1, 2023•21 citations

Approximating Time-Dependent Quantum Statistical Properties

Sara Bonella, Giovanni Ciccotti

Computing quantum dynamics in condensed matter systems is an open challenge due to the exponential scaling of exact algorithms with the number of degrees of freedom. Current methods try to reduce the cost of the calculation using classical dynamics as the key ingredient of approximations of the quantum time evolution. Two main approaches exist, quantum classical and semi-classical, but they suffer from various difficulties, in particular when trying to go beyond the classical approximation. It may then be useful to reconsider the problem focusing on statistical time-dependent averages rather than directly on the dynamics. In this paper, we discuss a recently developed scheme for calculating symmetrized correlation functions. In this scheme, the full (complex time) evolution is broken into segments alternating thermal and real-time propagation, and the latter is reduced to classical dynamics via a linearization approximation. Increasing the number of segments systematically improves the result with respect to full classical dynamics, but at a cost which is still prohibitive. If only one segment is considered, a cumulant expansion can be used to obtain a computationally efficient algorithm, which has proven accurate for condensed phase systems in moderately quantum regimes. This scheme is summarized in the second part of the paper. We conclude by outlining how the cumulant expansion formally provides a way to improve convergence also for more than one segment. Future work will focus on testing the numerical performance of this extension and, more importantly, on investigating the limit for the number of segments that goes to infinity of the approximate expression for the symmetrized correlation function to assess formally its convergence to the exact result.

Entropy•Dec 27, 2013•3 citations

Benchmarking the Cumulant Expansion Method Using Dicke Superradiance

Martin Fasser, Claudiu Genes, Helmut Ritsch +1 more

Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task requiring approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a newly found exact solution. The application of two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The long-term population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show unphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, on longer time scales the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when used for long-time evolution or reliably finding steady states.

Photonics•Oct 10, 2025•0 citations

Quantum-mechanical cumulant expansions and their application to phase-space and to phase distributions

Alfred Wünsche

Starting from the characteristic function of an operator, we investigate cumulant expansions in quantum optics and apply them to two-dimensional distributions for the canonical variables of the phase space in the case of one degree of freedom (Wigner quasiprobability and its Fourier transform, uncertainty matrix) and to one-dimensional distributions (phase operator, time evolution operator to Hamiltonian). In the relations between cumulants and moments, we make emphasis on the central moments of an operator. It is shown that the determinant of the uncertainty matrix (modified uncertainty product) is invariant with respect to rotation and squeezing of the state in the phase space, whereas the uncertainty sum is only invariant with respect to rotations. We examine some problems for exponentials of the phase operator and show how mean values and variances are connected with the cumulants. The Hilbert–Schmidt distance of a state during time evolution to an initial state is discussed by cumulants.

Physica Scripta•Jun 1, 2015•8 citations
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.
In the study of optical properties of large atomic systems, a weak laser driving is often assumed to simplify the system dynamics by linearly coupled equations. Here we investigate the light-scattering properties of atomic ensembles beyond weak-field excitation through the cumulant expansion method. By progressively incorporating higher-order correlations into the steady-state equations, an enhanced accuracy can be achieved in comparison to the exact solutions obtained by solving a full density matrix. Our analysis reveals that, in the regime of weak dipole-dipole interaction, the first-order expansion yields satisfactory predictions for optical depth, while denser atomic configurations necessitate consideration of higher-order correlations. As the intensity of incident light increases, atom saturation effects become noticeable, giving rise to significant changes in light transparency, energy shift, and decay rate. This saturation phenomenon extends to subradiant atom arrays even under weak driving conditions, leading to substantial deviations from the linear model. Our findings demonstrate that the mean-field model is a good extension to linear models as it balances both accuracy and computational complexity. However, the crucial role of higher-order cumulants in large and dense atom systems remains unclear, since it is challenging theoretically owing to the exponentially increasing Hilbert space in such light-matter interacting systems.
In the last couple of years much effort has been dedicated to the development and the investigation of hybrid quantum devices, which combine the advantages of very dissimilar quantum systems for the realization of efficient quantum computation and communication technologies. Particularly interesting in this regard are solid-state quantum memories based on ensembles of spins (e.g. nitrogen vacancy defects in diamond or rare-earth doped crystals) coupled to superconducting microwave cavities. Here the spin-ensemble acts as a robust memory where the collective coupling to the cavity mode allows for the coherent transfer of quantum information. The major downside of ensembles inside solid-state systems is their natural tendency to exhibit inhomogeneous broadening of the transition frequencies, which makes experiments as well as their theoretical description very challenging. On the theoretical side it is primarily the computational complexity due to the exponential growth of the associated Hilbert space with the system size that inevitably calls for approximation schemes. To date, the theoretical description of large spin-cavity systems is mainly limited to weak driving fields, where the number of spin-excitations in the ensemble is negligible, corresponding to a Holstein-Primakoff approximation, or to mean-field approaches, where correlations within the system are neglected. In this thesis we aim to develop a model that accurately accounts for the dynamics of large inhomogeneously broadened spin ensembles, coupled to a single cavity mode, even in the case of strong driving fields. In particular we employ the generalized cumulant expansion method to account for correlations within the system, which allows us to go beyond the Holstein-Primakoff and mean-field approximations. As a first step, we demonstrate the applicability of the cumulant expansion technique to the well-studied Jaynes-Cummings model and investigate its main limitations. We then move on to the Tavis-Cummings model with multiple spins inside the cavity and generalize our approach to ensembles containing a very large number of spins (ca. 10 12). Our model aims not only to provide a proper description of ensemble-based quantum memories, but also to serve as a novel tool for the investigation of interesting new physics arising from cooperative effects in inhomogeneous ensembles.
For a model with many-to-one connectivity it is widely expected that mean-field theory captures the exact many-particle $N\ensuremath{\rightarrow}\ensuremath{\infty}$ limit, and that higher-order cumulant expansions of the Heisenberg equations converge to this same limit whilst providing improved approximations at finite NNN. Here we show that this is in fact not always the case. Instead, whether mean-field theory correctly describes the large-NNN limit depends on how the model parameters scale with NNN, and the convergence of cumulant expansions may be nonuniform across even and odd orders. Further, even when a higher-order cumulant expansion does recover the correct limit, the error is not monotonic with NNN and may exceed that of mean-field theory.
Computing quantum dynamics in condensed matter systems is an open challenge due to the exponential scaling of exact algorithms with the number of degrees of freedom. Current methods try to reduce the cost of the calculation using classical dynamics as the key ingredient of approximations of the quantum time evolution. Two main approaches exist, quantum classical and semi-classical, but they suffer from various difficulties, in particular when trying to go beyond the classical approximation. It may then be useful to reconsider the problem focusing on statistical time-dependent averages rather than directly on the dynamics. In this paper, we discuss a recently developed scheme for calculating symmetrized correlation functions. In this scheme, the full (complex time) evolution is broken into segments alternating thermal and real-time propagation, and the latter is reduced to classical dynamics via a linearization approximation. Increasing the number of segments systematically improves the result with respect to full classical dynamics, but at a cost which is still prohibitive. If only one segment is considered, a cumulant expansion can be used to obtain a computationally efficient algorithm, which has proven accurate for condensed phase systems in moderately quantum regimes. This scheme is summarized in the second part of the paper. We conclude by outlining how the cumulant expansion formally provides a way to improve convergence also for more than one segment. Future work will focus on testing the numerical performance of this extension and, more importantly, on investigating the limit for the number of segments that goes to infinity of the approximate expression for the symmetrized correlation function to assess formally its convergence to the exact result.
Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task requiring approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a newly found exact solution. The application of two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The long-term population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show unphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, on longer time scales the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when used for long-time evolution or reliably finding steady states.
Starting from the characteristic function of an operator, we investigate cumulant expansions in quantum optics and apply them to two-dimensional distributions for the canonical variables of the phase space in the case of one degree of freedom (Wigner quasiprobability and its Fourier transform, uncertainty matrix) and to one-dimensional distributions (phase operator, time evolution operator to Hamiltonian). In the relations between cumulants and moments, we make emphasis on the central moments of an operator. It is shown that the determinant of the uncertainty matrix (modified uncertainty product) is invariant with respect to rotation and squeezing of the state in the phase space, whereas the uncertainty sum is only invariant with respect to rotations. We examine some problems for exponentials of the phase operator and show how mean values and variances are connected with the cumulants. The Hilbert–Schmidt distance of a state during time evolution to an initial state is discussed by cumulants.