---
title: Multiparticle Cumulant Method in Flow Analysis
url: https://www.emergentmind.com/topics/multiparticle-cumulant-method
type: topic
---

# Multiparticle Cumulant Method in Flow Analysis

The multiparticle cumulant method is a connected-correlation framework in which lower-order disconnected contributions are subtracted from multiparticle observables so that the remaining quantity isolates genuine collective or correlated structure. In relativistic collision studies, it is most often formulated for azimuthal anisotropy through cumulants such as $c_n\{2\}$, $c_n\{4\}$, symmetric cumulants, and asymmetric cumulants, implemented with Q-vectors and, in small systems, with pseudorapidity-separated subevents to reduce non-flow. Closely related cumulant constructions are also used for identified-particle multiplicity fluctuations, multibin long-range correlations, multi-particle balance functions, and femtoscopic many-body correlation functions [1802.00374] [1101.1907] [2310.07618] [2107.10227].

## 1. Statistical basis and connected correlations

At its most general level, the method is organized by generating functions and their logarithms. For multiplicity fluctuations of a species with event-by-event multiplicity $N$, the ordinary cumulants $C_n$ are defined from the moment generating function $M(t)=\langle e^{tN}\rangle$ and the cumulant generating function $K(t)=\ln M(t)$ through
\[
C_n=\left.\frac{\partial^n K(t)}{\partial t^n}\right|_{t=0}.
\]
The same formalism admits factorial moments and factorial cumulants through
\[
F(s)=\langle(1+s)^N\rangle,\qquad G(s)=\ln F(s),\qquad \kappa_n=\left.\frac{\partial^n G(s)}{\partial s^n}\right|_{s=0},
\]
with
\[
G(s)=K(\ln(1+s)),\qquad K(t)=G(e^t-1).
\]
The ordinary and factorial cumulants are related by Stirling-number expansions such as
\[
C_n=\sum_{k=1}^n S(n,k)\kappa_k,\qquad \kappa_n=\sum_{k=1}^n s(n,k)C_k
\]
[2510.15353].

For long-range multiplicity correlations in several bins, the basic observables are the joint factorial moments
\[
F_{n_1,\dots,n_m}=\left\langle \prod_{i=1}^m N_i^{(n_i)}\right\rangle,
\]
where $N_i^{(n)}=N_i(N_i-1)\cdots(N_i-n+1)$, together with the multivariate factorial moment generating function
\[
H(\{x_i\})=\left\langle \prod_{i=1}^m (1+x_i)^{N_i}\right\rangle.
\]
The joint factorial cumulants are obtained from $K(\{x_i\})=\ln H(\{x_i\})$ by differentiation at the origin. This construction makes explicit that mixed cumulants vanish if the bins are statistically independent, whereas nonzero mixed cumulants signal connected long-range structure [1101.1907].

The same connected-correlation logic appears in differential correlation functions. Ordinary cumulant densities $C_k(x_1,\dots,x_k)$ are obtained from inclusive densities $\rho_k$ by cluster expansion, and integral factorial cumulants are their phase-space integrals. In the balance-function framework, this distinction between differential cumulants and their acceptance-integrated counterparts is essential because differential objects retain the full $(p_T,y,\phi)$ dependence while integral cumulants compress that information into acceptance-dependent fluctuation measures [2310.07618].

## 2. Azimuthal cumulants and the Q-cumulant formalism

For azimuthal anisotropy, the starting point is the Fourier expansion
\[
\frac{dN}{d\phi}\propto 1+2\sum_{n=1}^{\infty} v_n\cos[n(\phi-\Psi_n)].
\]
Multiparticle azimuthal correlators are built from particle angles and averaged first over distinct tuples within an event and then over events. For harmonic $n$, the fundamental two- and four-particle correlators are
\[
\langle\!\langle e^{in(\phi_1-\phi_2)}\rangle\!\rangle=\langle\cos[n(\phi_1-\phi_2)]\rangle,
\]
\[
\langle\!\langle e^{in(\phi_1+\phi_2-\phi_3-\phi_4)}\rangle\!\rangle=\langle\cos[n(\phi_1+\phi_2-\phi_3-\phi_4)]\rangle.
\]
The corresponding standard cumulants are
\[
c_n\{2\}=\langle\!\langle e^{in(\phi_1-\phi_2)}\rangle\!\rangle,
\]
\[
c_n\{4\}=\langle\!\langle e^{in(\phi_1+\phi_2-\phi_3-\phi_4)}\rangle\!\rangle-2\langle\!\langle e^{in(\phi_1-\phi_2)}\rangle\!\rangle^2.
\]
In the pure-flow limit,
\[
\langle\!\langle e^{in(\phi_1-\phi_2)}\rangle\!\rangle\approx \langle v_n^2\rangle,\qquad
\langle\!\langle e^{in(\phi_1+\phi_2-\phi_3-\phi_4)}\rangle\!\rangle\approx \langle v_n^4\rangle,
\]
so that
\[
v_n\{2\}=\sqrt{c_n\{2\}},\qquad v_n\{4\}=(-c_n\{4\})^{1/4}\quad\text{for }c_n\{4\}<0
\]
[1802.00374].

The standard implementation uses event-wise flow vectors
\[
Q_n=\sum_{j=1}^{M} e^{in\phi_j}
\]
or, in weighted analyses,
\[
Q_n=\sum_{i=1}^{M} w_i e^{in\phi_i}.
\]
For unit weights, the two-particle correlator is
\[
\langle 2\rangle_n=\frac{|Q_n|^2-M}{M(M-1)}.
\]
For weighted reference particles, one convenient form is
\[
c_n\{2\}=\frac{|Q_n|^2-\sum_i w_i^2}{(\sum_i w_i)^2-\sum_i w_i^2}.
\]
Higher-order correlators are expressed through combinations of $Q_n$, $Q_{2n}$, and higher harmonics, and can be generated recursively in a single loop over particles and events [2112.12236] [2212.12282].

This formalism extends beyond the basic two- and four-particle case. The recursive cumulant relations
\[
c_n\{2\}=\langle\!\langle 2\rangle\!\rangle_n,\quad
c_n\{4\}=\langle\!\langle 4\rangle\!\rangle_n-2\langle\!\langle 2\rangle\!\rangle_n^2,
\]
\[
c_n\{6\}=\langle\!\langle 6\rangle\!\rangle_n-9\langle\!\langle 4\rangle\!\rangle_n\langle\!\langle 2\rangle\!\rangle_n+12\langle\!\langle 2\rangle\!\rangle_n^3,
\]
\[
c_n\{8\}=\langle\!\langle 8\rangle\!\rangle_n-16\langle\!\langle 6\rangle\!\rangle_n\langle\!\langle 2\rangle\!\rangle_n-18\langle\!\langle 4\rangle\!\rangle_n^2+144\langle\!\langle 4\rangle\!\rangle_n\langle\!\langle 2\rangle\!\rangle_n^2-144\langle\!\langle 2\rangle\!\rangle_n^4
\]
generalize systematically to higher order [2212.12282] [2005.07974].

## 3. Subevent cumulants and the suppression of non-flow

In small systems, low multiplicity makes non-flow especially problematic. Jets, dijets, resonance decays, Coulomb/HBT correlations, and global momentum conservation can dominate standard cumulants. The subevent cumulant method addresses this by splitting the event into pseudorapidity-separated regions and requiring particles in a correlator to come from different subevents [1802.00374].

For two subevents $A$ and $B$,
\[
c_n\{2\}^{\text{sub}}=\langle \cos[n(\phi_a-\phi_b)]\rangle,\qquad
\langle 2\rangle_n^{A|B}=\frac{\mathrm{Re}(Q_n^A(Q_n^B)^*)}{M_A M_B}.
\]
For three subevents, a particularly clean four-particle construction is
\[
\langle 4\rangle_n^{A,A|B|C}=\langle\!\langle e^{in(\phi_1^A+\phi_2^A-\phi_3^B-\phi_4^C)}\rangle\!\rangle,
\]
\[
c_n\{4\}^{\text{three-sub}}=\langle 4\rangle_n^{A,A|B|C}-2\langle 2\rangle_n^{A|B}\langle 2\rangle_n^{A|C}.
\]
The same logic extends to symmetric cumulants,
\[
SC(m,n)=\langle v_m^2 v_n^2\rangle-\langle v_m^2\rangle\langle v_n^2\rangle,
\]
and to their two- and three-subevent variants [1802.00374].

The physical rationale is explicit. Requiring particles to come from different $\eta$-separated subevents removes same-jet and resonance pairs, three-subevent constructions break typical dijet topologies, and $\eta$ gaps reduce Coulomb/HBT contributions. In AMPT for $p$+Pb at $\sqrt{s_{NN}}=5.02$ TeV, the standard method gave $c_2\{4\}$ that changed sign, becoming positive at smaller $N_{\rm ch}$ where non-flow dominates, whereas both two- and three-subevent $c_2\{4\}$ remained negative across the full $N_{\rm ch}$ range. In the same study, standard-method $SC(2,3)$ and $SC(2,4)$ were comparable to data, while the subevent results behaved differently, indicating contamination of the standard method, especially when the number of produced particles is small [1802.00374].

Small-system measurements at RHIC reinforce the same point in a different way. In $d$+Au collisions, PHENIX observed real-valued $v_2\{4\}$, that is $c_2\{4\}<0$, at $\sqrt{s_{NN}}=200$, $62.4$, $39$, and $19.6$ GeV, while in $p$+Au at $200$ GeV the four-particle cumulant had the opposite sign. At $200$ GeV in $d$+Au, $c_2\{6\}$ agreed with $c_2\{4\}$, so that $v_2\{6\}\approx v_2\{4\}$, which PHENIX interpreted as evidence that nonflow effects are subdominant in the four-particle signal [1707.06108].

## 4. Mixed-harmonic, differential, and high-order generalizations

The method extends naturally from single-harmonic cumulants to correlations among different harmonics. The symmetric cumulant
\[
SC(m,n)=\langle\!\langle e^{im(\phi_1-\phi_2)+in(\phi_3-\phi_4)}\rangle\!\rangle
-\langle\!\langle e^{im(\phi_1-\phi_2)}\rangle\!\rangle\langle\!\langle e^{in(\phi_1-\phi_2)}\rangle\!\rangle
\]
measures whether events with larger $v_m$ tend to have larger or smaller $v_n$ [1802.00374]. The asymmetric cumulant
\[
ac_2\{3\}=\left\langle e^{i\,2(\phi_1+\phi_2-2\phi_3)}\right\rangle
\]
projects the nonlinear coupling between $v_2$ and $v_4$, while the four-particle symmetric cumulant
\[
sc_{2,4}\{4\}=\left\langle e^{i\,2(\phi_1-\phi_2)+i\,4(\phi_3-\phi_4)}\right\rangle
-\left\langle e^{i\,2(\phi_1-\phi_2)}\right\rangle\left\langle e^{i\,4(\phi_3-\phi_4)}\right\rangle
\]
probes correlations between $v_2$ and $v_4$. Their normalized forms,
\[
nsc_{2,4}\{4\}=\frac{sc_{2,4}\{4\}}{v_2\{2\}^2\,v_4\{2\}^2},\qquad
nac_2\{3\}=\frac{ac_2\{3\}}{v_2\{2\}^2\sqrt{v_4\{2\}^2}},
\]
were calculated from transverse momentum conservation and flow and found to be in good agreement with recent ATLAS subevent measurements in $pp$ collisions [2403.05782].

For rare probes, the cumulant construction becomes differential. In PbPb collisions at $\sqrt{s_{NN}}=5.02$ TeV, prompt $D^0$ mesons were analyzed with a two-subevent reference built from HF$-$ and HF$+$ towers. The differential four-particle cumulant is
\[
d_n\{4\}=\langle 4'\rangle-2\langle 2'\rangle\langle 2\rangle,
\]
and the corresponding four-particle flow coefficient is
\[
v_n\{4\}(D^0)= -\,\frac{d_n\{4\}}{(-c_n\{4\})^{3/4}}.
\]
The reported $v_2\{4\}(D^0)$ rose with $p_T$ to a maximum near $p_T\approx 3.5$ GeV and remained systematically below $v_2\{2\}$, providing direct access to event-by-event fluctuations of charm-quark azimuthal anisotropy [2112.12236].

The same logic generalizes to arbitrary mixtures of particles of interest and inclusive reference particles. Mixed correlators of the form
\[
\langle m_P,m_R\rangle_n=\Big\langle e^{in\big(\sum_{a=1}^{m_P}\psi_{p,a}-\sum_{b=1}^{m_R}\phi_{r,b}\big)}\Big\rangle
\]
lead to mixed cumulants $d_n^{(P|R)}\{2\}$, $d_n^{(P|R)}\{4\}$, $d_n^{(P|R)}\{6\}$, and $d_n^{(P|R)}\{8\}$, together with rare-probe estimators such as
\[
v_n^P\{2\}=\frac{d_n^{(P|R)}\{2\}}{\sqrt{c_n\{2\}}},\qquad
v_n^P\{4\}=-\,\frac{d_n^{(P|R)}\{4\}}{[-c_n\{4\}]^{3/4}}.
\]
This construction was proposed specifically to study fluctuations and correlations in the azimuthal anisotropies of rare probes [2307.16796].

At higher order, several general frameworks now exist. A generic recursive algorithm produced explicit new 10-, 12-, 14-, and 16-particle single-harmonic cumulants and the corresponding $v_n$ coefficients [2005.07974]. A decomposition based on mathematical induction generated exact event-level expressions for high-order Q-cumulants and tabulated the needed integer coefficients and basis elements up to 14-particle order [2212.12282]. A later partition-based method made the fast calculation of arbitrary-order $v_n\{2k\}$ practical and was validated up to 40th order [2506.00535]. In parallel, a multidimensional generating-function analysis classified 33 distinct cumulants with orders $2,3,4,5$ and harmonics $2,3,4,5$, including the genuine three-particle correlation function that contains information of all third-order cumulants [2005.04742].

## 5. Multiplicity cumulants, balance functions, and femtoscopic connected correlations

Although azimuthal flow dominates present usage, the same cumulant logic is also central in fluctuation analyses of identified-particle multiplicities. In Ag+Ag collisions at $\sqrt{s_{NN}}=2.55$ GeV, HADES used ordinary cumulants $C_n$, factorial cumulants $\kappa_n$, and normalized factorial cumulants together with a probabilistic particle-identification method in which each track carries a species weight
\[
w_i(x_k)=\frac{f_i(x_k)}{\sum_j f_j(x_k)}.
\]
Event-wise proxies such as
\[
W_a=\sum_{k\in\text{event}} w_a(x_k)
\]
permit the reconstruction of factorial moments and cumulants without hard PID cuts. Because MC closure tests showed that binomial efficiency corrections fail in the studied phase space, HADES used unfolding and full detector-response corrections, and removed centrality fluctuations with a data-driven event-mixing procedure [2510.15353].

The differential-cumulant framework also yields multi-particle balance functions. For the two-particle case,
\[
B_2^{+-}(x_1,x_2)=\frac{1}{\langle N^{-}\rangle}\left[C_2^{+-}(x_1,x_2)-C_2^{--}(x_1,x_2)\right],
\]
while higher-order balance functions are built from connected $C_n$ combinations so that two- and three-prong resonance decays are suppressed. In full acceptance, these balance functions satisfy simple sum rules determined by quantum number conservation, and the paper argues that they arise as an intrinsic component of high-order net-charge cumulants [2310.07618].

A distinct but structurally related use of cumulants appears in femtoscopy. There, the measured three-body correlation function is decomposed as
\[
C_3(Q_3)=C_3^{2b}(Q_3)+c_3(Q_3),
\]
with
\[
c_3(Q_3)=C_3(Q_3)-C_3^{2b}(Q_3),
\qquad
C_3^{2b}(Q_3)=C_3^{12}(Q_3)+C_3^{23}(Q_3)+C_3^{31}(Q_3)-2.
\]
The lower-order contribution is obtained by projecting measured pair correlations into the triplet variable $Q_3$ through an analytic kernel
\[
W(k_1,Q_3)=\frac{16(\alpha\gamma-\beta^2)^{3/2}k_1^2}{\pi Q_3^4\gamma^2}
\sqrt{\gamma Q_3^2-(\alpha\gamma-\beta^2)k_1^2},
\]
which isolates the genuine three-body connected correlation $c_3$ after subtraction [2107.10227].

## 6. Interpretation, model discrimination, and open issues

A central virtue of the method is its sensitivity to fluctuations. Under widely used assumptions for the $v_2$ distribution,
\[
v_2\{2\}^2\approx \langle v_2\rangle^2+\sigma_{v_2}^2,\qquad
v_2\{4\}^2\approx \langle v_2\rangle^2-\sigma_{v_2}^2,
\]
so ratios such as $v_2\{4\}/v_2\{2\}$ directly probe fluctuation strength [2112.12236]. In $p$+Pb, AMPT reproduced $v_2\{2\}$ but yielded a $c_2\{4\}$ magnitude smaller than ATLAS data, which the study interpreted as either underestimated collectivity or missing non-Gaussian $v_n$ fluctuations [1802.00374].

The sign structure of cumulants is physically informative but not unique to a single mechanism. A parton-model initial-state calculation in $pA$ collisions reproduced negative $c_2\{4\}$, the ordering
\[
v_2\{2\}>v_2\{4\}\approx v_2\{6\}\approx v_2\{8\},
\]
and the qualitative behavior of symmetric cumulants, and argued that coherent multiple scattering is essential to obtain a positive definite $v_2\{4\}$ [1706.06260]. In a related CGC analysis, the 3-dipole correlator up to order $1/N_c^4$ led to the conjecture
\[
c_n\{m\}\sim \frac{1}{N_c^{2(m-1)}},
\]
with the dilute-limit calculation of $c_2\{4\}$ supporting that scaling [1901.01778]. These results imply that negative $c_2\{4\}$ and higher-order hierarchy do not by themselves uniquely identify hydrodynamic flow.

Global constraints can also bias cumulants. For transverse momentum conservation, the leading contribution to the $2k$-particle cumulant was derived as
\[
c_n\{2k\}=( -1)^{nk}\frac{(nk)!}{(n!)^{2k}}
\frac{p_1^n p_2^n\cdots p_{2k}^n}{(N-2k)^{nk}\langle p^2\rangle_F^{nk}},
\]
so the sign is $( -1)^{nk}$ and the magnitude falls as $(N-2k)^{-nk}$. The interplay of TMC and flow can change the signs of $c_2\{4\}$, $c_3\{2\}$, and $c_3\{4\}$ at some multiplicities [2204.01038]. The same interplay was shown to describe $sc_{2,4}\{4\}$ and $ac_2\{3\}$ measured with the subevent cumulant method, indicating that subevent cumulants suppress short-range non-flow but do not remove global TMC [2403.05782].

Taken together, these developments indicate that the multiparticle cumulant method is not a single observable but a hierarchy of connected-correlation constructions. In large systems it constrains the full joint probability density of flow magnitudes and symmetry planes; in small systems it separates collectivity from jets, resonances, and global constraints as far as the chosen subevent topology allows; in multiplicity, balance-function, and femtoscopic applications it isolates connected many-body structure from acceptance, PID overlap, and lower-order backgrounds. A plausible implication is that future progress will depend less on introducing entirely new observables than on combining high-order, mixed-harmonic, differential, and subevent cumulants with response-based modeling and detector-accurate inference frameworks [2005.07974] [2212.12282] [2506.00535].

Source: https://www.emergentmind.com/topics/multiparticle-cumulant-method