Brownian Inextensible Elastic Filaments

• Brownian inextensible elastic filaments are slender, flexible structures governed by elasticity, hydrodynamic drag, inextensibility, and thermal fluctuations, modeling both biological and synthetic polymers.
• The underlying mathematical framework incorporates elastic bending energy, nonlocal hydrodynamics, and stochastic forces, with dimensionless parameters delineating transitions among tumbling, buckling, and snaking regimes.
• Advanced numerical methods and experimental validations enhance our understanding of cytoskeletal mechanics, polymer rheology, and soft-matter systems in microfluidic environments.

A Brownian inextensible elastic filament is a slender, flexible object whose centerline dynamics are governed by elasticity, hydrodynamic drag in a viscous solvent, geometrically exact inextensibility, and thermal (Brownian) fluctuations. This model applies to a broad class of biological filaments, such as actin and microtubules, as well as synthetic semiflexible polymers in solution. The interplay of stochastic thermal forcing, flow-induced stresses, elastic bending, and the nonlinear inextensibility constraint yields a rich spectrum of nonequilibrium behaviors, central to cytoskeletal mechanics, polymer rheology, and soft-matter flows.

1. Mathematical Formulation and Physical Principles

The centerline of a filament is described by a time-dependent curve $\mathbf{r}(s,t)$, parameterized by arclength $s \in [0, L]$ and constrained to be inextensible, $|\partial_s \mathbf{r}| = 1$. The elastic energy functional is given by the Euler–Bernoulli form,

$$E_{\mathrm{bend}}[\mathbf{r}] = \frac{B}{2} \int_0^L |\partial_s^2 \mathbf{r}(s)|^2 \, ds,$$

where $B = k_B T\,\ell_p$ is the bending rigidity and $\ell_p$ is the persistence length. The resulting force per unit length is

$$\mathbf{f}(s, t) = B\,\partial_s^4 \mathbf{r} - \partial_s[\sigma(s)\partial_s \mathbf{r}] + \mathbf{f}^{\text{br}}(s, t),$$

where $\sigma(s)$ is a Lagrange multiplier enforcing inextensibility and $\mathbf{f}^{\text{br}}$ is a Gaussian white noise with covariance set by the fluctuation–dissipation theorem: $$\langle \mathbf{f}^{\text{br}}(s, t) \mathbf{f}^{\text{br}}(s', t')^T \rangle = 2 k_B T \Lambda^{-1} \delta(s - s') \delta(t - t').$$ The local, zero-Reynolds-number hydrodynamics are often represented by resistive force or slender-body theory, leading to a velocity–force relation that is nonlocal in $s$ and encodes anisotropic drag.

2. Dimensionless Parameters and Dynamical Regimes

The behavior of Brownian inextensible elastic filaments is controlled by three principal dimensionless numbers:

• Elasto-viscous number: $\bar{\mu} = 8\pi \mu \dot{\gamma} L^4 / B$ (steady shear) or $\bar{\mu}_m = 8\pi \mu \dot{\gamma}_m L^4 / (B c)$ (oscillatory shear), comparing viscous loading to elastic bending resistance. Here, $\mu$ is solvent viscosity, $\dot{\gamma}$ is shear rate, $L$ is filament length, and $c = -\ln(e \varepsilon^2)$ is the logarithmic slenderness parameter.
• Persistence ratio: $\ell_p/L$, quantifying filament semiflexibility versus length.
• Dimensionless period: $\rho = \dot{\gamma}_m T$, the ratio of imposed oscillation period to characteristic time scales.

Critical thresholds for morphological transitions are sharply defined in the high-stiffness limit: the global buckling threshold is $\bar{\mu}_c^{(1)}/c \simeq 306.4$; the “snaking” (U-turn) threshold is $\bar{\mu}_c^{(2)}/c \simeq 1700$ (Liu et al., 2018, Bonacci et al., 2022, Bonacci et al., 2 Jan 2026).

3. Morphological Transitions under Shear and Flow

In steady shear, filaments exhibit a sequence of morphologies as the elasto-viscous number increases (Liu et al., 2018):

1. Quasi-periodic tumbling (“Jeffery” regime): For $\bar{\mu} \lesssim \bar{\mu}_c^{(1)}$, filaments remain nearly straight and rotate periodically.
2. Buckled (“C-shaped”) regime: For $\bar{\mu}$ exceeding the global buckling threshold, a linear instability develops, and filaments buckle, transiently forming C-shaped conformations.
3. Snaking (“U-turn”) regime: Above the higher threshold $\bar{\mu}_c^{(2)}$, strongly deformed conformations with localized, propagating high-curvature regions (“snaking” or U-turns) emerge. This transition is captured by solvability of a minimal J-shape model (Liu et al., 2018).

In time-dependent (oscillatory) shear, the phase space becomes richer: reversibility, chaos, and “attractor-hopping” are observed, with stochastic switching between symmetry-related limit cycles (Bonacci et al., 2 Jan 2026). Flexibility and thermal fluctuations amplify orientation noise, resulting in irreversibility for sufficiently large $\rho$ and $\bar{\mu}_m$.

4. Role of Thermal Fluctuations and Persistence Length

Thermal fluctuations (Brownian noise) are introduced as additive Gaussian white noise forces, consistent with fluctuation–dissipation and hydrodynamics (Nedelec et al., 2009, Moreau et al., 2017, Maxian et al., 2024). The stochastic amplitude scales as $\sqrt{L/\ell_p}$, making noise effects most pronounced in filaments close to or below their persistence length. However, the sharp dynamical thresholds $\bar{\mu}_c^{(1,2)}$ are essentially independent of persistence ratio for $\ell_p/L \gtrsim 1$; noise simply broadens transitions into finite-width regimes and enhances rare events such as attractor switching (Liu et al., 2018, Bonacci et al., 2 Jan 2026).

5. Numerical Methods and Computational Schemes

Computational approaches for simulating Brownian inextensible elastic filaments must handle stiff elastic forces, nonlinear inextensibility constraints, hydrodynamic coupling, and detailed-balance–consistent stochasticity (Nedelec et al., 2009, Moreau et al., 2017, Maxian et al., 2024). Key algorithms include:

Method (Ref)	Enforcement of Inextensibility	Hydrodynamics
Projection/Lagrange Multiplier (Nedelec et al., 2009)	Explicit constraints via projection or Lagrange multipliers at each time step	Local drag, optionally nonlocal
Asymptotic Coarse-Graining (Moreau et al., 2017)	Built-in via tangent-angle discretization	Resistive Force Theory (RFT)
Spectral-Chebyshev, SDP Mobility (Maxian et al., 2024)	Saddle-point system for tangents and multipliers	Nonlocal Rotne-Prager regularized, SPD-split

Implicit/IMEX time-stepping is necessary to deal with elastic stiffness. Efficient projection methods (e.g., orthogonal or saddle-point projections) maintain inextensibility to numerical precision.

Special quadrature and fattened regularization allow accurate and scalable nonlocal hydrodynamics for large filament ensembles (Maxian et al., 2024). Dynamic steric (excluded volume) and cross-linker interactions are added for network and bundle simulations.

6. Experimental and Theoretical Validation

Experiments employing fluorescent actin filaments in microfluidic shear flows demonstrate quantitative agreement between observed morphologies, tumbling dynamics, and numerically simulated transitions (Liu et al., 2018, Bonacci et al., 2022). Parameter regimes for actin typically span $L=4$–$40\,\mu$m, $\ell_p\approx17~\mu$m, $\mu\approx5.6~$mPa·s, and shear rates $\dot{\gamma}\approx0.5$–$10~$s⁻¹, resulting in elasto-viscous numbers spanning $O(10^2)$–$O(10^7)$.

Phase diagrams in $(\bar{\mu}/c, \ell_p/L)$ space identify sharp boundaries for the three canonical regimes—Jeffery tumbling, global buckling, and snaking U-turns—their topology robust to moderate Brownian noise (Liu et al., 2018). Advanced simulations capture equilibrium shape statistics, end-to-end distributions, and transition rates between dynamic attractors (Maxian et al., 2024, Bonacci et al., 2 Jan 2026).

7. Applications, Implications, and Extensions

Brownian inextensible elastic filament models underlie predictive theories for cytoskeletal rheology, flagellar propulsion, single-filament manipulation, and microfluidic sorting/alignment strategies (Liu et al., 2018, Bonacci et al., 2022). The same formalism generalizes to cross-linked actin networks, bundles (enhanced power transduction via flexibility (Perilli et al., 2018)), active matter, and collective Langevin systems (Nedelec et al., 2009, Maxian et al., 2024).

A plausible implication is that the universal structure of elastic instabilities and fluctuating constraints informs improved constitutive models for both biological and synthetic complex fluids, and guides optimal control of filament-driven soft-matter systems. The robustness of key transitions and the ability to simulate large numbers of filaments with full nonlocal hydrodynamics and thermodynamic consistency provides a foundation for future studies of emergent, nonequilibrium behaviors in disordered, confined, and interacting filamentous assemblies.