Charge-Charge Correlation (QQC) Overview

• Charge-charge correlation (QQC) is a family of observables quantifying how charge-carrying degrees of freedom fluctuate and co-vary in diverse physical systems.
• QQC plays a crucial role in QCD thermodynamics by linking generalized susceptibilities and mixed cumulants to critical dynamics and hadronic interactions.
• QQC techniques also bridge theory and experiment in studies of charge order in cuprate superconductors and photon correlations in mesoscopic and optical systems.

Charge-charge correlation denotes a class of observables that quantify how charge-carrying degrees of freedom fluctuate, co-vary, or organize collectively. The precise object is strongly context dependent. In QCD thermodynamics it usually means generalized susceptibilities or mixed cumulants of conserved baryon number, electric charge, and strangeness; in correlated-electron systems it can denote charge-order correlations tied to the reconstructed electronic structure; in mesoscopic and optical systems it can denote correlations induced by a shared fluctuating charge environment; and in collider theory it can denote detector-level charge-flow correlators or related event shapes (Wen et al., 2019, Zhao et al., 2017, Purohit et al., 2015, Monni et al., 1 Aug 2025). This suggests that QQC is best understood as a family of charge-sensitive correlators rather than a single canonical observable.

1. Formal definitions and scope

In equilibrium QCD, the basic QQC objects are generalized susceptibilities of the pressure with respect to dimensionless chemical potentials. A standard definition is

$$\chi_{ijk}^{BQS}(T,\mu_B,\mu_Q,\mu_S)= \frac{\partial^{i+j+k}(p/T^4)}{\partial \hat\mu_B^i\,\partial \hat\mu_Q^j\,\partial \hat\mu_S^k}, \qquad \hat\mu_X=\mu_X/T,$$

with mixed cumulants such as

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$

and analogous higher-order cumulants $\chi_{21}^{XY}$, $\chi_{12}^{XY}$, $\chi_{31}^{XY}$, $\chi_{22}^{XY}$, and $\chi_{13}^{XY}$ for $XY\in\{BQ,BS,QS\}$ (Wen et al., 2019). A related formulation uses the conserved-flavor covariance matrix

$\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$

with $a,b\in\{u,d,s\}$, so that diagonal elements measure same-flavor fluctuations and off-diagonal elements measure whether different charges are carried together by the same effective degrees of freedom (Pratt, 2012).

Outside equilibrium thermodynamics, the same terminology refers to different but structurally related observables. In deep-inelastic scattering, the one-point charge correlator in the Breit frame is defined by

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$0

and its azimuthally differential version by

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$1

with an operator realization through a charge-flow detector built from the electromagnetic current $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$2 (Cao et al., 25 Jun 2026). In $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$3 annihilation, the pairwise charge-charge correlation is defined as

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$4

and can be rewritten as a correlator of light-ray charge operators $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$5 (Monni et al., 1 Aug 2025).

These definitions are not interchangeable. In some settings QQC is a thermodynamic derivative, in others a two-point detector observable, a pair distribution, or a fluctuation-induced proxy. What they share is a common task: isolating how charge information is encoded in the microscopic state and transmitted to measurable correlators.

2. Conserved-charge QQC in equilibrium QCD

The QCD literature uses QQC most systematically for correlations among baryon number $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$6, electric charge $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$7, and strangeness $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$8. In a $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$9 flavor Polyakov-loop improved quark-meson model solved with the functional renormalization group, mixed cumulants are computed up to fourth order under heavy-ion conditions, with the central conclusion that higher-order mixed charge correlations are much more sensitive to critical dynamics than second-order ones (Wen et al., 2019). In a $\chi_{21}^{XY}$0 flavor PNJL analysis, fourth-order mixed correlations such as $\chi_{21}^{XY}$1, $\chi_{21}^{XY}$2, $\chi_{21}^{XY}$3, together with $\chi_{21}^{XY}$4, $\chi_{21}^{XY}$5, and $\chi_{21}^{XY}$6, change rapidly through the chiral crossover, become discontinuous at the first-order transition, and diverge at the critical point; all studied fourth-order correlations except $\chi_{21}^{XY}$7 and $\chi_{21}^{XY}$8 are described as excellent probes of the QCD critical point (Fu et al., 2010).

The physical interpretation is degree-of-freedom sensitive. In an independent parton gas,

$\chi_{21}^{XY}$9

so off-diagonal correlators vanish. In a hadron gas,

$\chi_{12}^{XY}$0

and off-diagonal structure is generically nonzero because hadrons bundle several conserved charges into the same excitation (Pratt, 2012). Lattice gauge calculations cited in this context show that above the crossover region the off-diagonal elements largely disappear, while the diagonal elements approach about $\chi_{12}^{XY}$1 of the parton-gas expectation, which supports a picture in which up, down, and strange charges move largely independently in the quark-gluon plasma (Pratt, 2012).

Particular channels isolate specific physics. The ratio

$\chi_{12}^{XY}$2

roughly measures how much of strangeness is tied to baryonic degrees of freedom, since strange baryons carry both $\chi_{12}^{XY}$3 and $\chi_{12}^{XY}$4 whereas strange mesons carry strangeness but no baryon number (Bellwied et al., 2019). The baryon-electric charge correlator

$\chi_{12}^{XY}$5

is unusually sensitive to hadronic interactions: in an isospin-symmetric system the $\chi_{12}^{XY}$6 strange baryons do not contribute, because the $\chi_{12}^{XY}$7 is neutral and the $\chi_{12}^{XY}$8 triplet cancels, so the observable becomes a particularly clean probe of the pion-nucleon sector (Lo et al., 2017). In an S-matrix treatment based on empirical phase shifts, the interaction contribution to $\chi_{12}^{XY}$9 at $\chi_{31}^{XY}$0 MeV is suppressed by about $\chi_{31}^{XY}$1 relative to the hadron resonance gas model, and this improves agreement with lattice QCD in the hadronic and crossover region up to about $\chi_{31}^{XY}$2 MeV (Lo et al., 2017).

3. Heavy-ion phenomenology, magnetic fields, and experimental proxies

For heavy-ion applications, QQC is rarely evaluated in a completely unconstrained ensemble. A common choice imposes strangeness neutrality,

$\chi_{31}^{XY}$3

and a fixed electric-charge to baryon-density ratio,

$\chi_{31}^{XY}$4

which determines $\chi_{31}^{XY}$5 and $\chi_{31}^{XY}$6 self-consistently (Wen et al., 2019). In hybrid descriptions of heavy-ion collisions, the microscopic bookkeeping is written as

$\chi_{31}^{XY}$7

with the balancing part constrained by

$\chi_{31}^{XY}$8

The nonlocal piece is evolved through a hydrodynamic QGP stage and then through a hadronic Boltzmann cascade, after which it is projected onto charge balance functions in $\chi_{31}^{XY}$9, $\chi_{22}^{XY}$0, and $\chi_{22}^{XY}$1 (Pratt et al., 2018). A complementary formulation relates the final hadronic correlators

$\chi_{22}^{XY}$2

to the pre-hadronization charge matrix $\chi_{22}^{XY}$3, the hadronic susceptibility matrix $\chi_{22}^{XY}$4, and the QGP diffusion width $\chi_{22}^{XY}$5, thereby providing an explicit bridge from lattice-like susceptibilities to measurable hadronic balance functions (Pratt, 2012).

This equilibrium-to-experiment bridge motivates hadronic proxies. A lattice-plus-HRG study proposed

$\chi_{22}^{XY}$6

as a proxy for $\chi_{22}^{XY}$7, because $\chi_{22}^{XY}$8 baryons dominate the measured $\chi_{22}^{XY}$9 sector while kaons dominate the strangeness variance among measured species; at $\chi_{13}^{XY}$0 the proxy tracks the full ratio well over the temperature range around the crossover (Bellwied et al., 2019). Experimental measurements by ALICE use second-order cumulant ratios

$\chi_{13}^{XY}$1

built from net-proton, net-kaon, and net-charge multiplicities. These observables deviate from the Poissonian baseline, show a significant impact of resonance decays, and are described better by a canonical Thermal-FIST treatment with finite correlation volume than by a grand-canonical treatment, which underscores the role of local charge conservation (Collaboration, 24 Mar 2025).

Magnetic-field sensitivity has become a distinct QQC subtopic. Lattice QCD with physical-mass HISQ quarks in external magnetic fields finds that the baryon-electric charge correlator $\chi_{13}^{XY}$2 along the transition line starts to increase rapidly for

$\chi_{13}^{XY}$3

and rises by about a factor of $\chi_{13}^{XY}$4 at

$\chi_{13}^{XY}$5

while at $\chi_{13}^{XY}$6 MeV the continuum-estimated enhancement is about a factor of $\chi_{13}^{XY}$7 at the same field strength (Ding et al., 2023). The same study finds significant magnetic-field dependence in $\chi_{13}^{XY}$8, with the normalized leading-order double ratio reaching about $\chi_{13}^{XY}$9 in Au+Au and Pb+Pb and about $XY\in\{BQ,BS,QS\}$0 in Ru+Ru at $XY\in\{BQ,BS,QS\}$1 (Ding et al., 2023). A later three-flavor PNJL analysis at $XY\in\{BQ,BS,QS\}$2 extended this logic to fourth order and concluded that $XY\in\{BQ,BS,QS\}$3 is more magnetic-field sensitive than the second-order $XY\in\{BQ,BS,QS\}$4 and more sensitive than the other fourth-order mixed correlations considered, making it a better magnetometer of QCD in that model (Mao et al., 14 May 2026).

Charge-sensitive observables also appear in searches for anomalous transport such as the chiral magnetic effect. In the Color Glass Condensate, initial-state $XY\in\{BQ,BS,QS\}$5 correlations generate a rapidity-independent negative pedestal and a rapidity-dependent Pauli-blocking contribution that is positive at small $XY\in\{BQ,BS,QS\}$6 and negative at larger $XY\in\{BQ,BS,QS\}$7, providing a same-charge background to CME-sensitive observables in small systems (Kovner et al., 2017). STAR’s charge multiplicity asymmetry correlations further show that the observed charge separation signal could not be explained by CME alone, appears proportional to event-by-event $XY\in\{BQ,BS,QS\}$8, and is consistent with zero in events with $XY\in\{BQ,BS,QS\}$9, which argues for a substantial background from event anisotropy and correlated particle production (Wang, 2012).

4. Charge order and reconstructed electronic structure in cuprates

In correlated-electron materials, QQC can refer not to conserved-charge susceptibilities but to charge-order correlations emerging from low-energy electronic structure. In a $\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$0 analysis of cuprate superconductors, charge-charge correlations are treated as a manifestation of a charge-order instability of the reconstructed electron Fermi surface, with the central claim that the characteristic charge-order wave vector is controlled by the second-neighbor hopping $\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$1 (Zhao et al., 2017).

Within the charge-spin separation fermion-spin formulation, the normal-state electron Green’s function is written as

$\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$2

with bare dispersion

$\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$3

and spectral function

$\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$4

The low-energy spectral weight at $\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$5 forms two contours, $\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$6 and $\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$7, associated with the front and back sides of a Fermi pocket created by the pseudogap. Antinodal weight is strongly suppressed, leaving disconnected segments whose tips meet at hot spots. The charge-order wave vector is identified with the hot-spot-to-hot-spot vector,

$\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$8

so charge order is not introduced phenomenologically as an independent density modulation but emerges from an instability of the reconstructed Fermi surface (Zhao et al., 2017).

A useful reformulation of the self-energy,

$\chi_{ab}\equiv \frac{\langle Q_a Q_b\rangle}{V},$9

shows that the pseudogap controls the quasiparticle scattering rate through

$a,b\in\{u,d,s\}$0

with the pseudogap largest near the antinode and smallest near the hot spots. Increasing $a,b\in\{u,d,s\}$1 from $a,b\in\{u,d,s\}$2 to $a,b\in\{u,d,s\}$3 at fixed doping $a,b\in\{u,d,s\}$4 shifts the hot spots toward the nodal region and changes the Fermi-pocket geometry so that $a,b\in\{u,d,s\}$5 increases approximately linearly with $a,b\in\{u,d,s\}$6. For $a,b\in\{u,d,s\}$7, $a,b\in\{u,d,s\}$8, and $a,b\in\{u,d,s\}$9, the calculated $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$00, close to the observed $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$01 in Bi$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$02Sr$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$03La$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$04CuO$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$05. The proposed family-to-family variation of $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$06 is therefore traced to differences in $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$07, rather than to nearly universal $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$08 and $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$09 (Zhao et al., 2017).

5. Mesoscopic, optical, and classical many-body realizations

In solid-state open quantum systems, QQC can mean correlations induced by a shared fluctuating charge environment rather than direct charge transport. For a pair of driven quantum-dot exciton qubits interacting with nearby charge traps, the total Hamiltonian includes the bare quantum-dot energies, photon bath, laser driving, QD-photon coupling, and QD-charge coupling, with the charge-induced Stark shift represented by couplings $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$10 between fluctuator $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$11 and dot $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$12 (Purohit et al., 2015). The central observable is the photon-photon cross-correlation

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$13

which, using input-output theory, becomes a correlator of the dot lowering and raising operators. A common fluctuator modulates both transition energies via the DC Stark effect and therefore induces measurable photon bunching or antibunching. Same-sign laser detunings yield positive correlation and bunching, opposite-sign detunings yield antibunching for sufficiently large detuning, and the decay of $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$14 occurs on the charge-switching timescale rather than the optical lifetime when the fluctuator is slow enough. With two fluctuators of different switching rates, the correlation function can develop a two-plateau structure, which allows several charges fluctuating at different rates to be distinguished (Purohit et al., 2015).

In classical Coulomb matter, charge-charge correlations are studied through pair distributions rather than cumulants. For $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$15 classical charges $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$16 in a harmonic trap,

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$17

the fluid-phase global pair-separation distribution is found to be remarkably close to the bulk three-dimensional one-component plasma pair correlation function for sufficiently large $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$18, despite strong shell structure in the density profile (Wrighton et al., 2011). Shell-resolved angular correlations map naturally onto chord-length or arc-length interpretations of bulk OCP correlations, and for particles constrained to a single spherical shell the agreement with the two-dimensional OCP on a sphere is described as remarkable (Wrighton et al., 2011). In the ordered phase, shell correlations are well reproduced by thermally broadened Thomson-problem point sets, so the relevant structure is no longer liquid-like but a disordered or broadened spherical Coulomb crystal (Wrighton et al., 2011).

These realizations broaden the meaning of QQC. In one case it is inferred indirectly from emitted photons and the structure of environmental noise; in the other it is a pair-separation statistic in an inhomogeneous Coulomb system. The common element is not a shared formula but a shared emphasis on how charge interactions imprint spatial or temporal correlations.

6. Charge correlators in conformal and collider field theory

Large-charge field theory introduces another use of charge correlators. In non-relativistic conformal field theory, a large-$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$19 operator insertion is represented semiclassically by a Goldstone master field,

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$20

which generates large-charge $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$21-point functions directly from the path integral (Beane et al., 2024). The corresponding $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$22 field takes the form of an emergent harmonic trap,

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$23

and the operator insertion is interpreted as a superfluid droplet whose spatial size scales with Euclidean time separation (Beane et al., 2024). In conformal collider physics at large charge, the two-point charge detector correlator has leading behavior

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$24

while the first connected correction is a nontrivial EFT prediction with collinear enhancement

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$25

which is attributed to phonon propagation and described as a sound-jet effect (Cuomo et al., 27 Mar 2025).

In deep-inelastic scattering, the Breit-frame one-point charge correlator is designed to measure angular charge flow using charged tracks only. It is infrared and collinear safe, factorizes in the forward limit into ordinary DIS hard coefficients times a new nonperturbative nucleon charge correlator $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$26, and in the back-to-back limit it reduces to standard TMD factorization with a charge-weighted jet function $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$27 (Cao et al., 25 Jun 2026). The collinear logarithms of the forward correlator are resummed to NLL accuracy, while the back-to-back TMD logarithms are resummed to $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$28LL for the unpolarized distribution and $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$29LL for the Sivers asymmetry; the singular terms agree with full QCD through $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$30 (Cao et al., 25 Jun 2026).

The $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$31 charge-charge correlator provides a still more direct charge analogue of the energy-energy correlator. Although generally divergent beyond leading order, it becomes infrared and collinear safe at leading power in the back-to-back limit $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$32, a property termed leading-power safety (Monni et al., 1 Aug 2025). In that regime the observable admits an SCET factorization theorem with hard functions and a QQC jet function,

$$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$33

which replaces the energy-weighted jet function of the EEC (Monni et al., 1 Aug 2025). The logarithmic behavior is determined analytically up to four loops, and the back-to-back QQC is resummed to $$\chi_{11}^{XY}=\frac{1}{VT^3}\langle(\delta N_X)(\delta N_Y)\rangle,$$34LL accuracy, with Event2 validation of the predicted singular terms (Monni et al., 1 Aug 2025). This establishes a precise sense in which charge flow, like energy flow, can become a high-precision event-shape observable once the relevant kinematic limit is identified.

Across these field-theoretic settings, QQC no longer refers to a thermodynamic susceptibility. Instead it refers to correlation functions generated by heavy charged operators, light-ray charge detectors, or angular charge-flow measurements. The shared structure is the extraction of charge transport or charge organization from symmetry-constrained correlation functions, often with factorization and EFT control.