Particle Displacement Fields

• Particle displacement fields are quantitative measures of particle translational deviations from reference states, driven by stresses, thermal fluctuations, and interactions.
• They employ tensor decompositions, Eshelby-like defect superpositions, and stochastic increments to characterize correlations across amorphous, granular, and Brownian systems.
• These fields underpin improved continuum mechanics, enable efficient simulations, and refine measurement techniques in image velocimetry for complex flow scenarios.

Particle displacement fields quantitatively describe the collective or individual translational deviations of particles from their initial or reference positions, arising due to applied stresses, thermal fluctuations, interactions, or flow fields. These fields play a central role in continuum and statistical descriptions of amorphous and crystalline solids, granular materials under load, suspensions, and microfluidic environments. Contemporary research formulates and characterizes such fields at multiple scales: via tensorial correlation functions in amorphous solids, as superpositions of Eshelby-type defects in networks, through stochastic increments in Brownian dynamics, or as deterministic drift fields in hydrodynamics and image-based velocimetry.

1. Correlation Functions and Tensor Decomposition in Amorphous Solids

In isotropic amorphous solids, particle-displacement correlations are fully characterized by two-point tensors that encode both transverse and longitudinal fluctuations. The real-space correlation tensor,

$$C_{ij}(\mathbf{x}) = \langle u_i(\mathbf{x}) u_j(0) \rangle,$$

admits the standard decomposition in $d$ dimensions,

$$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$

where $A(r)$ and $B(r)$ quantify the amplitudes of transverse and longitudinal correlations, respectively. In three-dimensional isotropic elasticity, both decay as $1/r$, with their relative magnitudes fixed by the Poisson ratio via the Lamé moduli $(\lambda,\mu)$ (Wittmer, 2024).

Fourier transforms yield

$$C_{ij}(\mathbf{q}) = k_1(q) \delta_{ij} + k_2(q) q_i q_j,$$

from which the longitudinal and transverse spectra,

$$C_L(q) = k_1(q) + k_2(q) q^2, \qquad C_T(q) = k_1(q),$$

relate to $q$-dependent longitudinal, $d$0, and shear, $d$1, moduli: $d$2 For $d$3 (large length scales), Gaussianity dominates, justifying continuum elasticity, but non-Gaussian corrections and anisotropy become prominent for $d$4, which marks the breakdown of the continuum description at roughly $d$5 particle diameters. Strain correlations display strong non-monotonic $d$6-dependence, with minima near $d$7 (Wittmer, 2024).

2. Non-Affine Displacement Fields and Defect Superposition in Granular Materials

Under shear, the non-affine displacement field $d$8 quantifies the deviation of each particle position from its affine prediction. In jammed disk packings subjected to athermal, quasistatic (AQS) shear, these non-affine fields are well-approximated by linear superpositions of quadrupolar Eshelby-like inclusions. Each inclusion corresponds to a Delaunay triangle in the packing that exhibits a local stiffness mismatch (Willmarth et al., 7 Jul 2025).

The far-field of a single quadrupole is

$d$9

where $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$0 is proportional to the magnitude of the eigenstrain. Multiple inclusions sum linearly to reconstruct the global non-affine field.

The likelihood that a single quadrupole dominates the response increases with pressure and is maximized when a single contact—nearly aligned with low-frequency vibrational modes—breaks. The core size and strength of the quadrupolar field are set by the local geometry and missing contact cluster. This microscopic-to-macroscopic derivation provides a mechanical basis for observed non-affine, highly inhomogeneous displacement fields in amorphous solids under deformation (Willmarth et al., 7 Jul 2025).

3. Stochastic Displacement Fields in Brownian Dynamics

In suspensions, stochastic particle displacements arise from the interplay of hydrodynamic interactions and thermal fluctuations. The Rotne-Prager-Yamakawa (RPY) tensor, split into real- and wave-space contributions, provides the mobility mapping

$$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$1

where $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$2 is the many-body mobility. By a "positively-split Ewald" (PSE) decomposition, displacements are generated as

$$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$3

with $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$4 and $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$5 independent Wiener processes. Real-space samples treat short-range correlations via a Krylov method; wave-space samples exploit FFT-based convolutions. This approach guarantees symmetry, positive-definiteness, and linear (or near-linear) computational complexity, supporting simulations with $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$6 particles (Fiore et al., 2016). The PSE construction maintains exact fluctuation-dissipation balance for arbitrary configurations.

4. Particle Displacement Fields in Hydrodynamic and Microfluidic Contexts

In microfluidic systems, particle displacement fields encode the deterministic and irreversible motion of suspended particles under steady flows at zero Reynolds number. The canonical equation of motion for a particle near boundaries is

$$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$7

where $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$8 denotes particle radius, $$C_{ij}(\mathbf{x}) = A(r) \delta_{ij} + B(r) \frac{x_i x_j}{r^2},$$9 the imposed flow, and $A(r)$0 wall-induced corrections. The net drift field,

$A(r)$1

emerges from Faxén and wall effects (Liu et al., 5 Aug 2025).

Key to sustained net displacement is the breaking of flow symmetries. In fully symmetric vortices, closed particle orbits preclude net motion. Introducing asymmetry—either in the flow or via obstacles—permits controlled systematic particle displacement, accumulation at fixed points, or capture at wall locations. The scaling of the net displacement in "deterministic lateral displacement" arrays follows

$A(r)$2

with $A(r)$3 and $A(r)$4 the obstacle size and aspect ratio, and $A(r)$5 the flow angle. This produces strongly size-selective separation, even in purely hydrodynamic regimes (Das et al., 15 Nov 2025).

5. Displacement Fields and Fluctuation Bounds in Gibbsian Systems

In two-dimensional Gibbsian particle systems with reasonable repulsive interactions, particle displacements relative to ideal lattice locations are at least of order $A(r)$6 for a system of size $A(r)$7. This lower bound, established by a measure-preserving shift transformation $A(r)$8, demonstrates that large-scale positional order cannot be maintained, a manifestation of the Mermin-Wagner theorem. The displacement field thus exhibits slow, typically unbounded, fluctuations dictated by system size and interaction range (Fiedler et al., 2019).

Such universal lower bounds apply across a broad class of models—including continuum Potts, Widom-Rowlinson, and Lennard-Jones systems—and imply the absence of true Bragg peaks and the divergence of mean-square displacement in two-dimensional equilibrium ensembles.

6. Particle Displacement Field Measurement in Imaging Velocimetry

Particle image velocimetry (PIV) maps particle displacement fields by registering sequential images, with conventional algorithms assuming straight-line motion. Recent advances—"diffeomorphic PIV"—redefine the displacement field $A(r)$9 as the streamline-integral of the estimated velocity field: $B(r)$0 where $B(r)$1 solves $B(r)$2 over $B(r)$3 (Lee et al., 2021). This construction accurately captures curved particle trajectories, reducing bias in high-curvature flows. The practical PIV workflow iteratively refines $B(r)$4 by minimizing measures of image misalignment,

$B(r)$5

where $B(r)$6 is a cross-correlation, L2, or neural network similarity metric. Quantitative tests confirm substantial reductions in velocity estimation error, especially in high-curvature regions.

7. Summary Table of Principal Contexts for Particle Displacement Fields

Context	Governing Mechanism	Mathematical Description
Amorphous solids/correlations	Elastic/tensor correlators	$B(r)$7
Granular packings/non-affine fields	Eshelby inclusion superposition	$B(r)$8
Brownian dynamics	Hydrodynamic + random increments	$B(r)$9 (PSE split)
Microfluidic hydrodynamics	Faxén+wall corrections, symmetry	$1/r$0, limit cycles
Gibbsian statistical systems	Entropic/interaction constraints	Displacement $1/r$1
Imaging velocimetry	Data-driven, streamline-integral	$1/r$2

This cross-disciplinary framework provides a unified basis for characterizing how and why particles displace relative to reference states across disordered, driven, thermal, or engineered environments, with precise mathematical models and measurement methodologies now available in each domain (Wittmer, 2024, Willmarth et al., 7 Jul 2025, Fiore et al., 2016, Liu et al., 5 Aug 2025, Das et al., 15 Nov 2025, Fiedler et al., 2019, Lee et al., 2021).