Fluctuations of driven probes reveal nonequilibrium transitions in complex fluids

Complex fluids subjected to localized microscopic energy inputs, typical of active microrheology setups, exhibit poorly understood nonequilibrium behaviors because of the intricate self-organization of their mesoscopic constituents. In this work we show how to identify changes in the microstructural conformation of the fluid by monitoring the variance of the probe position, based on a general method grounded in the breakdown of the equipartition theorem. To illustrate our method, we perform large-scale Brownian dynamics simulations of an effective model of micellar solution, and we link the different scaling regimes in the variance of the probe’s position to the transitions from diffusive to jump dynamics, where the fluid intermittently relaxes the accumulated stress. This suggests stored elastic stress may be the physical mechanism behind the nonlinear friction curves recently measured in micellar solutions, pointing at a mechanism for the observed multi-step rheology. Our approach overcomes the limitations of continuum macroscopic descriptions and introduces an empirical method, applicable in experiments, to detect nonequilibrium transitions in the structure of complex fluids.

💡 Research Summary

The paper addresses the challenge of probing nonequilibrium structural changes in complex fluids subjected to localized energy input, a situation typical of active microrheology. The authors propose a simple yet powerful diagnostic: monitor the variance of a driven probe’s position. In equilibrium, equipartition dictates that the variance of the probe’s displacement in a harmonic trap equals (k_BT/\kappa). When a constant‑velocity force is applied, this relation breaks down, and the excess variance (\Delta x(v)=\kappa,\mathrm{Var}(r_{0,x})/(k_BT)-1) becomes a quantitative marker of nonequilibrium activity.

To validate the idea, they construct a coarse‑grained two‑dimensional Brownian dynamics model of a micellar solution. The probe is a particle of radius (\sigma_0) pulled by a harmonic potential moving at speed (v). The surrounding fluid is represented by flexible chains of length (L) (“multi‑blob” model). Each chain consists of Gaussian‑repelling beads linked by harmonic springs, mimicking the elasticity and excluded volume of real polymers or micelles. The chains can effectively “break” when the probe’s repulsion forces two neighboring beads apart, thereby modeling scission events that release stored elastic stress.

Simulation results reproduce the experimentally observed multi‑step friction curves: the effective friction (\gamma_{\text{eff}}(v)=\kappa\langle r_{0,x}\rangle/v) shows a low‑velocity plateau, an intermediate plateau, and for longer chains a shear‑thickening peak. Crucially, the probe’s positional variance displays a strikingly similar pattern. At very low speeds the variance is indistinguishable from the equilibrium value; as speed increases, (\Delta x) grows as (v^2), a universal quadratic scaling dictated by symmetry (the linear term vanishes when switching to a frame moving with the opposite velocity).

A characteristic speed (v^(L)) separates the quadratic regime from a higher‑speed regime with a different power‑law. Empirically, (v^(L)\propto L^{-3/2}). This scaling matches the Peclet number (Pe=vR_g/D_L) where the chain diffusion coefficient (D_L\sim L^{-1}) (Rouse scaling) and the equilibrium gyration radius (R_g\sim L^{1/2}). The crossover at (Pe\approx1) marks a transition from diffusion‑dominated to advection‑dominated dynamics: chains no longer have time to diffuse away from the approaching probe, leading to pronounced deformation and stress buildup.

Beyond the quadratic regime, the variance’s growth slows, and simultaneously the friction curve exhibits a steep rise. Visualizing the local gyration tensor (G(r_{\text{cm}})) reveals that at intermediate (U=v/v^*) the chains stretch ahead of the probe and align with the flow, while at high (U) the chains in front become shorter and the local bead density drops, indicating stress release via “jump dynamics”. The authors argue that this jump dynamics—elastic energy stored in the deformed chains being released when it exceeds the Gaussian repulsion barrier—underlies the nonlinear friction peaks observed experimentally.

The discussion emphasizes that two distinct mechanisms govern the observed crossovers: (1) a diffusion‑advection crossover captured by the Peclet number, and (2) an elastic‑stress‑release crossover associated with chain scission‑like events. The proposed variance‑based indicator is experimentally accessible because microrheology setups already track probe trajectories with high precision. Thus, without any additional structural probes, one can detect nonequilibrium transitions, identify the onset of stress accumulation, and infer microscopic rearrangements in complex fluids.

Overall, the work provides a clear theoretical framework linking probe fluctuations to microstructural dynamics, validates it with large‑scale simulations, and offers a practical tool for experimentalists to uncover hidden nonequilibrium transitions in polymeric, micellar, or other soft matter systems.