Cumulant dynamics in finite-memory diffusion

This paper extends the standard Fickian diffusion model to Maxwell–Cattaneo diffusion to account for finite current-relaxation time, deriving closed evolution equations that reveal how this memory effect suppresses, shifts, and reshapes the non-monotonic behavior of conserved charge cumulants in the quark-gluon plasma.

The Gist

The Big Picture: The "Heavy" Current

Imagine a crowd of people (representing particles in a Quark-Gluon Plasma, or QGP) trying to move from a crowded room to an empty one.

In the standard, old-school way of thinking (called Fickian diffusion), the crowd moves instantly. As soon as someone sees a gap, they step into it immediately. The flow of people is perfectly synchronized with the empty space available. It's like a light switch: flip it, and the light turns on instantly.

However, the authors of this paper argue that in the extreme conditions of a heavy-ion collision (where a tiny fireball of super-hot matter is created), the "people" (the flow of charge) are actually a bit sluggish. They have inertia. When the crowd sees a gap, they don't step into it instantly; it takes a tiny fraction of a second for them to react, accelerate, and get moving.

This paper studies what happens when we account for that delay. They call this Maxwell-Cattaneo diffusion. It's like saying the current (the flow) has a "memory" of where it was a moment ago, rather than just reacting to where it is right now.

The Problem: The "Freeze-Out" Snapshot

In these high-energy physics experiments, the fireball expands and cools down incredibly fast. Eventually, it "freezes out"—the particles stop interacting and fly off to detectors. Scientists take a snapshot of this moment to count how many particles are in a specific window.

They aren't just counting the average number of particles; they are looking at the fluctuations (the randomness).

• Cumulants: Think of these as different ways to measure the "shape" of the crowd's randomness.
  • 2nd Cumulant (Variance): How much does the crowd size vary? (Is it always 100 people, or sometimes 90, sometimes 110?)
  • 3rd & 4th Cumulants (Skewness & Kurtosis): These measure if the crowd is lopsided or if there are extreme outliers. These are the "sensitive" detectors for finding a Critical Point (a special state of matter where the rules of physics change dramatically).

The Experiment: Running the Simulation

The authors built a mathematical model to simulate how these fluctuations evolve over the short life of the fireball. They compared two scenarios:

1. The Old Way (Fickian): The crowd reacts instantly.
2. The New Way (Maxwell-Cattaneo): The crowd has a "reaction time" (memory).

They ran this simulation along different paths through the "phase diagram" (a map of temperature and density), including paths that go right near the mysterious Critical Point.

The Findings: Why the Delay Matters

1. The "Lag" Effect
In the standard model, the crowd tries to keep up with the changing environment but falls slightly behind (a "diffusive lag").
In the new model, because the flow has inertia, it falls even further behind. It's like a heavy truck trying to turn a corner; it doesn't just turn slowly; it overshoots or undershoots because it can't stop or start instantly.

2. The Critical Point is a Bumpy Road
When the system is far from the Critical Point, the "bumpy road" (the changing environment) is smooth. The delay just makes the truck arrive a few seconds later. The results look mostly the same as the old model.

But when the system passes near the Critical Point, the road becomes very bumpy and erratic. The environment changes rapidly.

• The Result: The "heavy truck" (the current with memory) reacts very differently here. Instead of just lagging, it starts to oscillate (wobble) and reshape the fluctuations.
• The Analogy: Imagine trying to walk through a crowd that is suddenly pushing and pulling you. If you are light and fast (instant reaction), you adjust instantly. If you are heavy and slow (memory), you might stumble, wobble, or get pushed in a different direction than expected.

3. Higher-Order Numbers Tell the Story
The most important finding is that this "memory effect" is barely noticeable in simple counts (2nd cumulant). However, it dramatically changes the complex shapes (3rd and 4th cumulants).

• The paper shows that the "wobble" caused by the delay can shift the peaks and valleys of these complex measurements.
• It can even flip the sign (positive to negative) of the measurements in certain areas.

The Conclusion: Don't Ignore the "Heavy" Flow

The authors conclude that if scientists want to find the QCD Critical Point using these fluctuation measurements, they cannot assume the flow of particles is instant.

If they ignore the finite memory (the delay in the current), they might misinterpret the data. They might think a signal is coming from the Critical Point when it's actually just the "inertia" of the flow, or they might miss the Critical Point entirely because the signal looks different than the "instant" models predicted.

In short: The paper says that in the chaotic, fast-moving world of particle collisions, the flow of matter has a "reaction time." Ignoring this reaction time leads to a distorted picture of the most interesting physics happening near the Critical Point. To get the right answer, you have to treat the flow like a heavy truck, not a light switch.

Technical Summary

Technical Summary: Cumulant Dynamics in Finite-Memory Diffusion

Problem Statement
The interpretation of conserved charge fluctuations (specifically net-proton number fluctuations) as signatures of the Quantum Chromodynamics (QCD) critical point in heavy-ion collisions requires a dynamical description of how fluctuation correlators evolve during the finite lifetime of the quark-gluon plasma (QGP). The standard theoretical baseline for this evolution is Fickian diffusion, where the diffusive current responds instantaneously to the local density gradient. However, this instantaneous limit may fail to capture delayed-response effects when the current-relaxation time becomes comparable to the relaxation time of relevant fluctuation modes. While critical slowing down is known to delay the evolution of fluctuation correlators relative to their local equilibrium values, the specific impact of a finite relaxation time in the constitutive relation of the current itself (memory effects) on multi-point correlators and experimental cumulants remains to be fully quantified.

Methodology
The authors extend the standard Fickian diffusion framework to Maxwell–Cattaneo (MC) diffusion, where the diffusive current is treated as an independent dynamical variable that relaxes to its Fickian value over a finite time scale $\tau$. This introduces a "memory" effect into the charge transport dynamics.

1. Stochastic Formulation: The dynamics are formulated as a first-order Markovian stochastic system for the charge density $n$ and an auxiliary "generalized momentum" variable $g \equiv \partial_t n$ (related to the current). The system includes Gaussian white noise satisfying the fluctuation-dissipation relation (KMS condition).
2. Evolution Equations: By applying Itô calculus to the stochastic equations, the authors derive closed, deterministic evolution equations for the equal-time multi-point Wigner functions ($W_N$) for $N=2, 3, 4$. These equations describe the time evolution of connected correlators in an enlarged variable space.
3. Phenomenological Setup: The equations are solved numerically along representative trajectories in the QCD phase diagram (fixed baryon chemical potential $\mu$, time-dependent temperature $T(t)$). The thermodynamic inputs (susceptibilities $\chi_k$) and transport coefficients are parameterized using a mapping to the 3D Ising universality class to model the critical region.
4. Observables: The evolved correlators at the freezeout surface are converted into acceptance-dependent cumulants ($C_N$) and their ratios ($S\sigma = C_3/C_2$, $\kappa\sigma^2 = C_4/C_2$) using a specific acceptance kernel that maps the experimental momentum-rapidity window to a spatial interval in the effective 1D description.

Key Contributions

• Derivation of Closed Hierarchies: The paper provides the first derivation of closed evolution equations for three- and four-point Wigner functions within the Maxwell–Cattaneo framework, reducing to the known Fickian hierarchy in the limit $\tau \to 0$.
• Memory vs. Diffusive Lag: The work distinguishes between the "diffusive lag" (where correlators fail to track instantaneous equilibrium due to finite diffusion rates) and the additional "memory effect" introduced by finite current relaxation ($\tau > 0$).
• Impact on Higher-Order Cumulants: The study demonstrates that while finite current relaxation acts as a small correction to the variance ($C_2$) away from the critical point, it induces significant qualitative changes in higher-order cumulants ($C_3, C_4$) and their ratios, particularly near the critical region.

Results

• Two-Point Functions ($W_2$): For modes within the physical momentum window, finite $\tau$ introduces an additional delay relative to the Fickian baseline. If the relaxation time is large enough that the characteristic MC scale $q_\star \sim 1/\sqrt{\tau\gamma}$ falls within the window, the response can become weakly underdamped. Near the critical point, where the equilibrium target varies non-monotonically, this delay reshapes the time profile of the correlator.
• Higher-Point Functions ($W_3, W_4$): Higher-order correlators exhibit greater sensitivity to finite $\tau$ than $W_2$. The delayed response leads to shifts in the extrema (peaks/valleys), changes in magnitude, and modifications to the sign structure of the relaxation profiles. These effects are more pronounced for the larger relaxation time ($\tau_{large}$) and near the critical trajectory.
• Freezeout Cumulants:
  • $C_2$ (Variance): The behavior of $C_2$ is dominated by low-momentum modes. Finite memory shifts and broadens the cumulant curve relative to both the instantaneous equilibrium estimate and the Fickian baseline. The specific ordering depends on the freezeout condition and acceptance window.
  • $C_3$ and $C_4$: The higher-order cumulants show robust qualitative signatures of memory. The skewness ($S\sigma$) and kurtosis ($\kappa\sigma^2$) ratios exhibit shifted peaks and modified sign structures. Notably, $\kappa\sigma^2$ can undergo a qualitative change in its sign structure relative to the Fickian baseline over a substantial range of $\mu$.
  • Smoothing: While the fixed-$q$ correlators show underdamped oscillations, the acceptance integration smooths these oscillations. However, the net finite-$\tau$ effects remain visible in the height, position, and sign structure of the final cumulant curves.

Significance and Claims
The paper claims that finite current relaxation is not merely a quantitative correction to diffusion but introduces a genuine memory effect that can systematically distort the interpretation of critical signatures.

• Systematic Distortion: The authors argue that current memory can reshape the observable cumulants (location, height, and even sign of peaks) even if the underlying equation of state is fixed. This suggests that neglecting finite relaxation times could lead to misinterpretations of experimental data regarding the location and nature of the QCD critical point.
• Diagnostic Power: Higher-order cumulants and their ratios, particularly the kurtosis ratio $\kappa\sigma^2$, are identified as the sharpest diagnostics for these memory effects.
• Minimal Model: The work positions Maxwell–Cattaneo diffusion as a "minimal laboratory" to isolate how memory in constitutive relations affects fluctuation evolution, distinct from the more complex Hydro+ frameworks where memory arises from integrating out slow critical modes.

The authors acknowledge limitations, including the use of prescribed (rather than self-consistently solved) background trajectories, the restriction to a 1D effective description, and the lack of a full particlization map to hadronic observables. They suggest that future work should embed this mechanism in expanding backgrounds with energy-momentum fluctuations and connect more directly to experimental hadronic cumulants.