[Paper Review] Anomalous diffusion of distinguished particles in bead-spring networks

Scott A. McKinley|ArXiv.org|Nov 22, 2009

Granular flow and fluidized beds20 references3 citations

TL;DR

This paper establishes a rigorous mathematical framework for anomalous sub-diffusion in bead-spring networks by modeling distinguished particle dynamics as sums of Ornstein-Uhlenbeck processes (ΣOU processes). It proves that the anomalous exponent ν in the mean-squared displacement ⟨x²(t)⟩ ∼ t^ν arises from the spectral structure of the network's spring connectivity, demonstrating universality of the ν = 1/2 scaling in Rouse chains and showing that altering the graph dimension or introducing repulsive forces enables tunable exponents ν ∈ (0,1).

ABSTRACT

We consider the anomalous sub-diffusion of a class of Gaussian processes that can be expressed in terms of sums of Ornstein-Uhlenbeck processes. As a generic class of processes, we introduce a single parameter such that for any $ν\in (0,1)$ the process can be tuned to produce a mean-squared displacement with $\E{x^2(t)} \sim t^ν$ for large $t$. The motivation for the specific structure of these sums of OU processes comes from the Rouse chain model from polymer kinetic theory. We generalize the model by studying the general dynamics of individual particles in networks of thermally fluctuating beads connected by Hookean springs. Such a set-up is similar to the study of Kac-Zwanzig heat bath models. Whereas the existing heat bath literature places its assumptions on the spectrum of the Laplacian matrix associated to the spring connection graph, we study explicit graph structures. In this setting we prove a notion of universality for the Rouse chain's well-known $\E{x^2(t)} \sim t^{1/2}$ scaling behavior. Subsequently we demonstrate the existence of other anomalous behavior by changing the dimension of the connection graph or by allowing repulsive forces among the beads.

Motivation & Objective

• To provide a rigorous mathematical interpretation of the intermediate-time anomalous diffusion profile ⟨x²(t)⟩ ∼ t^ν observed in polymer Rouse chains and related networks.
• To establish that the anomalous exponent ν = 1/2 in the Rouse model is a universal consequence of the network's geometric structure, not specific to polymer dynamics.
• To generalize the Rouse model by allowing non-local spring interactions and repulsive forces, enabling a broader class of anomalous exponents.
• To demonstrate that the long-time behavior of distinguished particles in such networks can be described by ΣOU processes with tunable spectral parameters.
• To explore the scaling limits of these systems, showing convergence to stochastic PDEs such as the stochastic heat equation under appropriate rescaling.

Proposed method

• Model the dynamics of a distinguished particle in a bead-spring network as a ΣOU process: x(t) = c₀B₀(t) + Σₖ₌₁^{N−1} cₖ zₖ(t), where each zₖ satisfies dzₖ = −λₖ zₖ dt + dBₖ.
• Relate the diffusive spectrum {λₖ} to the spring connectivity graph’s Laplacian eigenvalues, using spectral theory to derive the mean-squared displacement (MSD) scaling.
• Introduce a generalized Rouse spectrum λₖ,ₙ = (k/N)^ρ τ₁⁻¹ to tune the anomalous exponent to ν = 1 − 1/ρ for any ρ > 1.
• Use Szegö’s theorem on eigenvalue distribution to analyze the asymptotic behavior of the MSD in the large-N limit.
• Incorporate repulsive forces by modifying the spring potential, leading to higher-order Fourier modes in the shape function φ(x) = sin²ⁿ(x), enabling ν = 1 − 1/(2n).
• Analyze the continuum limit by rescaling spring constants as N² and noise as √N, showing convergence to the stochastic heat equation ∂ₜu = Δu + W(dy,dt).

Experimental results

Research questions

• RQ1 What determines the anomalous sub-diffusive exponent ν in the mean-squared displacement of a distinguished particle in a bead-spring network?
• RQ2 Is the ν = 1/2 scaling in the Rouse chain model universal across different network geometries, or specific to linear polymers?
• RQ3 Can the anomalous exponent ν be tuned to arbitrary values in (0,1) by modifying the network’s spring connectivity or introducing repulsive forces?
• RQ4 What is the limiting behavior of the distinguished particle process when the network size N → ∞ under appropriate rescaling of spring constants and noise?
• RQ5 How does the inclusion of higher-order spatial derivatives (e.g., via repulsive forces) affect the roughness and scaling of the particle trajectories?

Key findings

• The mean-squared displacement of a distinguished particle in a Rouse chain exhibits the profile ⟨x²(t)⟩ ∼ t^ν with ν = 1/2 for intermediate times, and this behavior is universal across all networks with the same spectral structure.
• By generalizing the Rouse spectrum to λₖ,ₙ = (k/N)^ρ τ₁⁻¹, the model achieves any anomalous exponent ν = 1 − 1/ρ for ρ > 1, with the MSD scaling as ⟨x²(t)⟩ ∼ t^ν in the intermediate regime.
• Introducing repulsive forces via a shape function φ(x) = sin²ⁿ(x) leads to a spectral parameter ρ = 2n, resulting in an anomalous exponent ν = 1 − 1/(2n) for any n ∈ ℕ.
• The coefficient family {κⱼ} for such systems is given by κⱼ = (−1)^{j+1} ₙCⱼ, ensuring symmetry and enabling the construction of systems with desired ν.
• Under a rescaling where spring constants scale as N² and noise as √N, the discrete bead-spring system converges to the stochastic heat equation ∂ₜu = Δu + W(dy,dt), with the limiting process exhibiting ν = 1/2 at short times.
• The limiting process u(y,t) is locally rougher than Brownian motion, with finite quartic variation, and its sample paths are not of bounded variation.