advanced27 min studyLesson 360 in Physics

Let conserved momentum remember a molecular impulse

Connect collective fluctuations, correlation tails, and limits of simple relaxation models.

On this page

What you’ll understand

  • Explain why conserved momentum creates slow collective memory.
  • Integrate a model correlation tail to assess transport convergence.
  • Distinguish long-time tails from instabilities and numerical artifacts.

An exponential collision model predicts a single finite memory time. A momentum-conserving fluid can retain a much slower memory through its surrounding flow: a moving particle excites a collective disturbance that spreads and partly returns correlated motion to the particle. This produces long-time tails in equilibrium correlations. The effect connects kinetic and hydrodynamic reasoning, and it warns that a simple relaxation fit may describe early data while missing the asymptotic transport contribution.

Local relaxation leaves conserved fields behind

Collisions relax nonconserved velocity distortions rapidly, but they do not destroy total momentum. At long wavelengths a transverse momentum disturbance spreads diffusively through viscosity. Its characteristic radius grows as the square root of elapsed time, so the affected volume grows as time to half the spatial dimension. A fixed total impulse distributed through that increasing volume gives a local amplitude with an algebraic, rather than exponential, decay. Coupling a tagged velocity to this field creates persistent correlation.

Hydrodynamic fluctuations add fluctuating stress and heat flux to deterministic balance laws, with strengths linked to dissipation at thermal equilibrium. These terms represent unresolved microscopic fluctuations after coarse graining; they are not arbitrary external forcing. Averaging the noise to zero does not erase its effects on correlations, because products of fluctuations can survive averaging. Nonlinear couplings among slow modes are central to the tail mechanism and are absent from a purely independent single-particle relaxation ansatz.

∂tu⊥=ν∇2u⊥,R(t)∼νt,Cv(t)∼A t−d/2\partial_t\mathbf u_\perp=\nu\nabla^2\mathbf u_\perp,\qquad R(t)\sim\sqrt{\nu t},\qquad C_v(t)\sim A\,t^{-d/2}
ν is kinematic viscosity in m²/s, R a spreading length in m, and d spatial dimension. Cv is one-component velocity correlation in m²/s²; the scaling is a hydrodynamic asymptotic model, not a universal exact exponent in every confined or two-dimensional fluid.↔ Scroll wide equations horizontally.

A small tail can matter to an integral

The diffusion coefficient is the time integral of a velocity correlation. An algebraic tail can therefore contribute noticeably even after the correlation becomes visually small. In three dimensions a t⁻³ᐟ² tail is integrable and leaves a slowly decaying cutoff error. In a simple two-dimensional t⁻¹ approximation, the integral grows logarithmically with the upper cutoff. Refined self-consistent theories can introduce logarithmic changes, so the simple power law should not be presented as a universal final theorem for every two-dimensional fluid.

Finite simulations introduce their own cutoffs: periodic images return sound and momentum disturbances, and small systems restrict the longest wavelengths. A late-time feature that moves when box size changes can be a finite-size artifact. The original Alder–Wainwright calculations explicitly compared molecular dynamics with collective hydrodynamic flow and considered image interference. Reporting a diffusion coefficient from a finite trace therefore requires specifying the correlation integration window, boundary conditions, and evidence that the relevant tail or cutoff is controlled.

D(T)=∫0TCv(t) dt,Cv(t)=At−3/2 ⟹ D(∞)−D(T)=2ATD(T)=\int_0^TC_v(t)\,dt,\qquad C_v(t)=At^{-3/2}\ \Longrightarrow\ D(\infty)-D(T)=\frac{2A}{\sqrt T}
T is a late-time cutoff in s, distinct from thermodynamic temperature in this display. For the assumed three-dimensional tail, A has units m²/s¹ᐟ² and the omitted area has diffusivity units m²/s.↔ Scroll wide equations horizontally.

Slow decay is not spontaneous growth

A positive long-time tail is a decaying equilibrium memory, not proof of a kinetic instability. Instabilities require a specified nonequilibrium state and a growing perturbation under its dynamics. Similarly, an unstable discretization can create artificial growth even when every physical mode is damped. A convincing calculation separates physical growth rates, conservation errors, numerical resolution, and finite-size changes. Basic BGK or deterministic diffusion alone cannot reproduce all correlated fluctuation spectra merely by fitting one coefficient.

The sampling lab illustrates another limit of independent models. Reset and set success probability to 0.5 and sample size to 25; predict sample-mean variance 0.01. Reset and set probability to 0.5 and size to 100; predict 0.0025. Its exact binomial law assumes independent samples. Correlated fluid observations can have a larger variance and a different effective sample size. Playback only highlights exact outcomes, so it cannot generate correlated hydrodynamic time traces or discover a long-time tail.

Var⁡ ⁣(1N∑i=1NXi)=σ2N[1+2∑k=1N−1(1−kN)ρk]\operatorname{Var}\!\left(\frac1N\sum_{i=1}^NX_i\right)=\frac{\sigma^2}{N}\left[1+2\sum_{k=1}^{N-1}\left(1-\frac{k}{N}\right)\rho_k\right]
Xi is a stationary sampled observable with variance σ² and normalized lag correlation ρk. This identity assumes finite second moments; independent samples have ρk=0 for k>0. N here is sample count, not gas particle number.↔ Scroll wide equations horizontally.

Estimate an unmeasured tail

For times t≥1 s, an illustrative three-dimensional correlation tail is Cv(t)=0.020(t/(1 s))⁻³ᐟ² m²/s². A diffusion estimate integrates measured data only to 4 s. Find the missing tail area and how it changes if the cutoff is extended to 16 s.

  1. Write the amplitude with its units: A=0.020 m²/s¹ᐟ² in Cv=A t⁻³ᐟ². The tail is integrable at infinity, but its contribution from a finite cutoff must still be counted.

  2. The missing area after 4 s is 2A/√(4 s)=0.020 m²/s. This is a correction to the integrated diffusion coefficient, not an additional velocity or a correlation value at the cutoff.

  3. At 16 s the missing area is 2A/√(16 s)=0.010 m²/s. Extending observation time by a factor of four halves the remaining tail error; it does not reduce it by four as a naive inverse-time estimate would.

The result

The omitted diffusion contribution is 0.020 m²/s at a 4 s cutoff and 0.010 m²/s at 16 s under the specified tail model.

Explore the model

Interactive lab

Interactive lab · sampling

Locate the independence assumption

Reset and set p=0.5, n=25; predict sample-mean variance 0.01. Reset and set p=0.5, n=100; predict 0.0025.

Locate the independence assumptionExact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Current readings: Mean sample proportion: 0.5; Sample-mean variance: 0.0125; Uniform CDF error bound: ≤ 0.1062. All parameters are adjustable with the labeled controls below.Exact distribution of the mean of 20 independent Bernoulli variables with success probability 0.5; playback highlights outcome 0.00.250.50.75100.0580.1150.1730.23Sample mean · k / nProbability massNormal bin masses are an approximationk = 0 · exact probability = 9.54e-7
Exact binomial massNormal interval massHighlighted outcome / mean
Mean sample proportion
0.5
Sample-mean variance
0.0125
Uniform CDF error bound
≤ 0.1062
t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe exact binomial instrument has independent Bernoulli observations and highlights outcomes without resampling. It does not represent correlated fluid trajectories, fluctuating stress, or a long-time tail.

Check your understanding

Which observation would most directly warn that a late-time correlation feature is affected by finite box size?

Try it yourself

For an illustrative two-dimensional model Cv(t)=0.030(t/(1 s))⁻¹ m²/s² beyond 1 s, calculate the extra diffusion integral from 2 s to 8 s. Separately, two unit-variance samples have correlation 0.5; find their mean’s variance.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Integrate 1/t over a finite interval. For two samples, include both cross-covariance terms.

Compare with the solution

The tail contribution is 0.030 ln(8/2)=0.030 ln4≈0.04159 m²/s. Continued logarithmic growth follows this model, not a claim of a universal uncorrected two-dimensional exponent. The two-sample mean has variance (1+1+2×0.5)/4=0.75, compared with 0.5 for independent samples, showing why correlated observations do not obey the naive independent variance reduction.

Take it with you

  • Conserved collective fields can retain algebraic memory after local collision relaxation.
  • Transport estimates must control tail integrals and finite-size cutoffs.
  • Equilibrium tails, kinetic instabilities, and numerical growth are different phenomena.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.