Linear Viscoelasticity of Semiflexible Polymers with Hydrodynamic Interactions

This study utilizes Brownian dynamics simulations of coarse-grained bead-spring chains with hydrodynamic interactions to demonstrate that semiflexible polymers exhibit power-law stress relaxation and frequency-dependent moduli that align with theoretical predictions for rigid rods and experimental data across a broad range of stiffness and frequencies.

The Gist

Imagine a world inside a drop of water where tiny, microscopic strings are floating around. These aren't just any strings; they are biopolymers—the building blocks of life like DNA, collagen (which gives your skin its bounce), and the scaffolding inside your cells.

Some of these strings are floppy and loose, like a cooked spaghetti noodle. Others are stiff and rigid, like a dry uncooked spaghetti noodle or a fishing rod. The scientists in this paper wanted to understand exactly how these "stiff strings" behave when you wiggle them, stretch them, or let them relax. This behavior is called viscoelasticity—a fancy word that means they act like a mix of a spring (elastic) and honey (viscous).

Here is the story of what they did, explained simply:

1. The Problem: The "Too Hard" vs. "Too Soft" Dilemma

To study these strings on a computer, scientists usually build a model. They imagine the string is made of little beads connected by springs.

• The "Too Soft" Model: If you use normal springs, the string stretches too easily. It's like a slinky. It doesn't act like a stiff rod, so the computer simulation gives the wrong answer for stiff biological strings.
• The "Too Hard" Model: If you make the string completely rigid (like a steel rod), the computer gets stuck. It's like trying to calculate the movement of a rigid stick in water; the math gets incredibly messy and slow because the stick can't bend at all, making the simulation crash or take forever.

The Goal: The authors wanted a "Goldilocks" model—a simulation that is stiff enough to act like a real biological rod, but flexible enough for the computer to run quickly and accurately.

2. The Solution: The "Smart Spring"

The team invented a new way to connect their beads using a special type of spring called a FENE-Fraenkel spring.

Think of this spring like a bungee cord with a hard stop.

• If you pull it gently, it stretches like a normal rubber band.
• But if you try to pull it too hard, it hits a "hard stop" and refuses to stretch any further, acting exactly like a rigid rod.

By tuning the settings of this "smart spring," they could make the computer model behave exactly like a stiff, unbreakable rod without the computer crashing. They also added a "bending penalty" (like a stiffener in a fishing rod) to control how much the string wanted to curl up.

3. The "Invisible Water" (Hydrodynamics)

One of the biggest challenges in studying these strings is the water they swim in. When a string moves, it pushes the water, which pushes back on the string. This is called hydrodynamic interaction.

• The Old Way: Many scientists ignored this "push-back" to save time, assuming the water didn't matter much.
• The New Discovery: The authors found that for stiff rods (like a fishing pole), the water doesn't matter much; the rod is too heavy and stiff to be bothered by the water's push.
• The Twist: But as the strings get more flexible (like a cooked noodle), the water becomes a huge factor. The water drags on the noodle, changing how it relaxes. The authors showed that if you ignore the water for flexible strings, your predictions are wrong.

4. The "Relaxation Dance"

The team simulated what happens when you stretch these strings and then let go. They watched how the "stress" (the tension) relaxed over time. They found a beautiful pattern:

• The Stiff Rods: When you let go of a stiff rod, it relaxes in a very specific, predictable way. It's like a stiff ruler snapping back. The math showed a specific "power law" (a mathematical rule) that matched previous theories perfectly.
• The Flexible Strings: When you let go of a floppy string, it relaxes differently, more like a tangled ball of yarn uncoiling.
• The Middle Ground: The most interesting part was the strings that were somewhat stiff. They showed a mix of behaviors. The authors mapped out exactly how the relaxation speed changes as the string goes from "stiff" to "floppy."

5. Does it Match Reality?

Finally, they compared their computer results with real-world experiments done on actual proteins (like collagen and PBLG).

• The Result: Their "Smart Spring" model matched the real experimental data almost perfectly across a wide range of speeds (frequencies).
• The Winner: The old theories (which assumed the strings were perfectly rigid) started to fail when the strings got a bit more flexible. The new model, however, worked for everything from stiff rods to floppy chains.

The Big Takeaway

This paper is like building a better simulator for a video game.

• Before, you could only simulate stiff sticks or floppy noodles, but not the "in-between" stuff that real life is made of.
• Or, you could simulate the stiff sticks, but the game ran so slowly it was unusable.
• Now, the authors have created a tool that runs fast, handles the "invisible water" correctly, and accurately predicts how real biological strings behave.

This helps scientists understand how cells move, how tissues heal, and how to design better synthetic materials that mimic nature's strength and flexibility. It's a bridge between the rigid world of math and the messy, fluid world of biology.

Technical Summary

1. Problem Statement

Semiflexible biopolymers (e.g., F-actin, collagen, hyaluronic acid) are critical to cellular mechanics, yet a comprehensive quantitative theory for their linear viscoelastic response across a wide range of bending stiffness remains elusive.

• The Gap: Existing models face a trade-off. Bead-rod models accurately capture inextensibility and rod-like behavior but are computationally expensive, especially when incorporating Hydrodynamic Interactions (HI). Bead-spring models are computationally efficient but typically fail to reproduce the correct high-frequency viscoelastic response because they lack the inextensibility constraint required for stiff chains.
• The Challenge: Developing a mesoscopic model that can replicate the linear viscoelastic response of both rigid rods and flexible chains in dilute solutions while explicitly including fluctuating hydrodynamic interactions, without the prohibitive computational cost of bead-rod simulations.

2. Methodology

The authors employed Brownian Dynamics (BD) simulations of coarse-grained bead-spring chains in the limit of infinite dilution.

• Model Formulation:

  • Chain Representation: A single chain consisting of $N_b$ beads connected by springs.
  • Spring Force Law: The study utilizes the FENE-Fraenkel spring force law. This is a hybrid potential that combines the Finite Extensible Nonlinear Elastic (FENE) and Fraenkel (non-zero rest length) laws. By tuning the spring stiffness ($H_R$) and extensibility parameters, this force law can mimic an inextensible rod (bead-rod behavior) while retaining the computational simplicity of a bead-spring model.
  • Bending Potential: A bending potential ($U_B \propto 1 - \cos\theta$) is applied to control chain stiffness, defined by the ratio of contour length to persistence length ($L/l_p$).
  • Hydrodynamic Interactions (HI): HI are incorporated using the Rotne-Prager-Yamakawa (RPY) tensor, which accounts for long-range fluid dynamic coupling between beads.
  • Excluded Volume: Modeled using the Soddemann-Duenweg-Kramer (SDK) potential (validated in athermal solvents), though excluded volume effects were largely neglected in the main study to match theta-solvent theoretical comparisons.

• Simulation Protocol:

  • Simulations were performed using a first-order Euler integration scheme.
  • Parameters were varied across a wide range of $L/l_p$ (from flexible $L/l_p \to \infty$ to stiff $L/l_p = 0.125$) and spring stiffness ($H_R$ from $10^2$ to $10^7$).
  • Validation: The code was validated against Multi-Particle Collision Dynamics (MPCD) simulations by Nikoubashman et al. regarding dynamic properties (mean square displacement and end-to-end vector autocorrelation).

• Property Calculation:

  • Relaxation Modulus $G(t)$: Calculated via the Green-Kubo relation from the stress tensor autocorrelation function.
  • Dynamic Moduli ($G', G''$): Obtained by Fourier transforming $G(t)$.
  • Zero-Shear Viscosity: Calculated by integrating $G(t)$.

3. Key Contributions

1. Development of a Unified Mesoscopic Model: The paper demonstrates that the FENE-Fraenkel bead-spring model can successfully replicate the linear viscoelastic response of inextensible bead-rod chains when appropriate spring stiffness and extensibility parameters are chosen. This bridges the gap between computationally expensive bead-rod models and standard bead-spring models.
2. Systematic Study of Hydrodynamic Interactions: This is the first study to comprehensively evaluate the effect of fluctuating hydrodynamic interactions on the linear viscoelasticity of semiflexible chains across the entire spectrum of stiffness (from rigid rods to flexible coils).
3. Validation Against Theory and Experiment: The model is rigorously validated against:
   • Semiflexible Rod Theory (SRT): Analytical predictions for stiff, free-draining rods.
   • Existing Simulations: Bead-rod simulation data.
   • Experimental Data: Dynamic moduli of Poly-γ-benzyl-L-glutamate (PBLG) and Collagen in infinitely dilute solutions.

4. Key Results

A. Free-Draining Limit (No HI)

• Spring Law Selection: FENE-Fraenkel springs were found to be superior to Hookean or standard FENE springs. They reproduce the distinct two-stage relaxation (fast spring relaxation followed by slow rod-like relaxation) and the terminal plateau predicted by SRT.
• Stiffness and Discretization: Increasing spring stiffness ($H_R$) and the number of beads ($N_b$) extends the time window where the chain behaves as a rigid rod.
• Power Law Scaling: The stress relaxation modulus $G(t)$ exhibits distinct power-law behaviors at intermediate times:
  • Stiff Chains ($L/l_p \ll 1$): Exponent $\alpha \approx -5/4$ (consistent with SRT).
  • Flexible Chains ($L/l_p \gg 1$): Exponent $\alpha \approx -1/2$ (consistent with Rouse theory).
  • Transition: As flexibility increases, the $-5/4$ regime vanishes, and the initial $-3/4$ regime (predicted for very short times in stiff rods) becomes accessible for more flexible chains.
• Long-Time Behavior: Rigid chains exhibit orientational relaxation (rigid rod behavior), while flexible chains exhibit Rouse-like relaxation.

B. Effect of Hydrodynamic Interactions (HI)

• Stiff Chains: For stiff semiflexible chains ($L/l_p \lesssim 10$), hydrodynamic interactions have a negligible effect on the intermediate power-law scaling. The behavior remains dominated by bending modes.
• Flexible Chains: As chains become more flexible ($L/l_p > 10$), HI significantly alter the dynamics.
  • Free-Draining: Intermediate slope $\alpha = -1/2$ (Rouse).
  • With HI: Intermediate slope shifts to $\alpha = -2/3$ (Zimm).
• Crossover: The study identifies a crossover regime around $L/l_p \approx 10$ where the influence of HI becomes significant, even though the chains are still technically "semiflexible."

C. Comparison with Experiments

• The model shows excellent quantitative agreement with experimental data for PBLG and Collagen across a broad frequency range.
• Unlike the Semiflexible Rod Theory (SRT), which deviates from experimental data at higher frequencies (especially for chains with $L/l_p > 1$), the FENE-Fraenkel bead-spring model accurately captures both low- and high-frequency behaviors.

5. Significance

• Computational Efficiency: The proposed model offers a computationally tractable alternative to bead-rod simulations, enabling the study of semiflexible polymers with hydrodynamic interactions over long time scales and large system sizes.
• Unified Framework: It provides a single theoretical framework that seamlessly transitions from rigid rod dynamics to flexible coil dynamics, capturing the correct power-law regimes for both limits.
• Biological Relevance: By accurately modeling the viscoelasticity of biopolymers like actin and collagen in dilute solutions, the model aids in understanding the mechanical properties of the cytoskeleton and extracellular matrices.
• Foundation for Future Work: The methodology establishes a robust basis for future studies on concentrated solutions, entangled networks, and crosslinked systems, where hydrodynamic screening and inter-chain correlations become critical.

In conclusion, the paper successfully validates the FENE-Fraenkel bead-spring model as a powerful tool for simulating the linear viscoelasticity of semiflexible polymers, effectively resolving the long-standing challenge of incorporating hydrodynamic interactions into inextensible chain models.