Kubo response and fluctuation–dissipation
Relate a weak force response to equilibrium correlations, derive the overdamped harmonic example, and state the thermal and linear-response assumptions needed for the relation.
On this page
The question for today
What you’ll understand
- Interpret a causal response kernel and its units.
- Verify an equilibrium fluctuation–dissipation identity for a harmonic coordinate.
- Distinguish linear response, thermal consistency, and arbitrary driven behavior.
Equilibrium fluctuations contain information about how a system responds to a small external disturbance. This is the content of fluctuation–dissipation relations, with Kubo theory providing a general response framework. The relation is powerful because it connects measurements performed without driving to predictions about controlled forcing. It also has precise assumptions: a response around equilibrium, the correct observable and conjugate force, and a bath whose noise and dissipation are thermally consistent. We will verify the statement in a solvable overdamped trap.
A response kernel measures a controlled change
Let a weak force f couple through an energy perturbation . The change in mean position is a convolution of f with a retarded kernel R. Retarded means that the response vanishes before the perturbation: future forcing cannot affect an earlier observation. Since integrating R times force over time gives a length, R has units of length divided by force and time. Its time integral is the static susceptibility, with units of displacement per force. Confusing the kernel with this integrated quantity loses a timescale.
Linear response retains only the part proportional to the perturbation amplitude. It does not mean the unperturbed system has no fluctuations, nor that every driven problem has an exactly linear solution. A nonlinear system may have a valid weak-field response around one equilibrium while behaving differently under large forcing. In the harmonic example below, the mean equation happens to be exactly linear for any applied force within the model. That special solvability should not be promoted into a general claim about strongly driven matter.
Verify the equilibrium relation in a thermal trap
For and constant mobility mu, the Itô equation has drift and noise amplitude . Averaging gives a deterministic linear equation for the mean. An impulse force produces a response for positive t. Separately, equilibrium Ornstein–Uhlenbeck fluctuations have covariance . Differentiating this covariance for positive time reproduces the response after multiplication by minus inverse thermal energy.
The equality ties a dissipative mobility measured by forcing to a correlation measured while no force is applied. Its minus sign is essential because the covariance decays. The equilibrium variance fixes the static susceptibility through . These expressions concern a time-even classical coordinate with the stated coupling convention. Quantum Kubo relations involve operator commutators and frequency-dependent thermal factors, while velocity observables and magnetic fields require careful time-reversal conventions. One must match the formula to the observable and regime.
A mismatch is evidence to interpret
If an experiment violates the equilibrium relation, possibilities include active forcing, a nonstationary preparation, unresolved slow variables, measurement filtering, or an inappropriate response protocol. The mismatch alone does not identify a unique mechanism. One must first compare correlations and responses at matched timescales, units, and baseline conditions. Replacing temperature by a fitted number can sometimes match one observable or frequency range, but does not prove that the entire driven system has a thermodynamic temperature or obeys all equilibrium identities.
A system can also be far from equilibrium while admitting useful linear response about its driven steady state. That response is linear in an additional small perturbation, but need not obey the equilibrium correlation formula written above. Generalized response theories can involve extra dynamical correlations and information about steady currents. Conversely, fluctuation theorems in later lessons can remain valid for large driving when their own microscopic-reversibility assumptions hold. Linear-response restrictions and far-from-equilibrium identities are different statements with different domains of validity.
Work it through
Predict a step response from trap fluctuations
A thermal trap has , , and . A constant force 0.2 pN is switched on from equilibrium. Find the mean shift and unforced equilibrium covariance at .
The relaxation rate is . At the stated time the surviving exponential is one-half, so the step response has reached half its final displacement.
The final mean shift is . Integrating the response kernel gives the present shift .
The unforced equilibrium variance is . Its covariance across the stated lag is . This is an equilibrium reference correlation, not an extra displacement caused by the applied force.
The mean moves by 0.05 micrometres, while the unperturbed equilibrium covariance at the same time lag is 0.001 square micrometres.
The living notebook
Explore the model
Interactive lab · function
Stiffer quadratic confinement gives a larger restoring slope
For each fresh comparison, Reset and set BOTH Coefficient a = 0.5 and Position x = 1, then Reset and set a = 1 and x = 1. Predict tangent slopes 1 and 2 for y = ax². Inspect the paused tangent or Step and compare how doubling quadratic curvature changes the slope at fixed displacement.
- y at chosen x
- 0.608
- Tangent slope · 2ax
- 1.35
- Signed area · ax³/3
- 0.182
This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. 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 forThis deterministic quadratic shows stiffness geometry only. It has no thermal noise, correlation function, temperature, or response-kernel readout, so it cannot test fluctuation–dissipation numerically.
Pause & reason
Check your understanding
Your turn
Try it yourself
A thermal coordinate has measured equilibrium covariance , where t is in seconds, and . Assuming the harmonic overdamped model, find k, mu, D, and the initial response-kernel value.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The covariance amplitude gives kBT/k; its decay rate gives mu times k. The initial response kernel equals mobility.
Compare with the solution
The stiffness is . The decay rate gives . Thus and . Differentiating C gives in the same units, independently verifying the relation.
Keep the idea
Take it with you
- A retarded kernel describes a causal response to a specified conjugate perturbation.
- Equilibrium fluctuations predict response only with the correct thermal and symmetry assumptions.
- Linear response around a driven state need not satisfy an equilibrium fluctuation–dissipation formula.
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.