Two-dimensional active polar semiflexible polymer under shear flow

This study employs numerical simulations to demonstrate that shear flow induces distinct conformational changes, tumbling dynamics, and rheological behaviors—including negative viscosity—in two-dimensional active polar semiflexible polymers, with activity significantly influencing their response until high shear rates restore passive-like characteristics.

The Gist

Imagine a long, wiggly noodle made of a semi-stiff material, like a piece of cooked spaghetti or a human hair. Now, imagine this noodle is alive. It has tiny motors running along its length, pushing it forward, making it swim on its own. This is what scientists call an "active semiflexible polymer."

In this paper, researchers put these "living noodles" into a fluid and started stirring the fluid (creating shear flow). They wanted to see how the noodle's shape, movement, and the fluid's thickness changed when the noodle was both trying to swim on its own and being dragged by the current.

Here is the story of what they found, broken down into simple concepts:

1. The Setup: The Tangled Dance Floor

Think of the fluid as a crowded dance floor.

• The Noodle (Polymer): It's not a floppy piece of string (flexible); it's a bit stiff, like a garden hose.
• The Motors (Activity): The noodle has little engines on it that push it forward. If the dance floor is still, the noodle might curl up into a tight spiral or a ball because it's so active and energetic.
• The Current (Shear Flow): The researchers started moving the dance floor in one direction (like a conveyor belt).

2. The Three Acts of the Noodle's Life

Act I: The "Unfolding" (Low Speed)
When the conveyor belt moves slowly, the noodle is still curling up into tight spirals because of its own energy. But as the belt speeds up just a little, the current acts like a hand gently pulling the noodle straight.

• The Result: The tight spirals unravel. The noodle stretches out long and thin, aligning perfectly with the direction of the flow. It's like a dancer who was spinning in a tight circle suddenly being pulled into a long, straight line by a partner.

Act II: The "Tumbling" (Medium Speed)
As the conveyor belt gets faster, things get chaotic. The noodle can't stay straight. It starts to fold over itself, forming shapes like a "U" or an "S." It flips and flops, tumbling end-over-end.

• The Surprise: In this middle zone, the noodle shrinks much faster in the direction perpendicular to the flow than a normal, non-living noodle would. It's as if the noodle is actively fighting the current by curling up tighter than it naturally would.

Act III: The "Passive Surrender" (Very High Speed)
If the conveyor belt goes super fast, the noodle's own motors can't keep up. The current is too strong. The noodle stops fighting, stops tumbling, and just gets dragged along like a dead stick. It behaves exactly like a normal, non-living noodle. The "life" of the noodle is drowned out by the speed of the water.

3. The Weird Physics: "Negative Viscosity"

This is the most mind-bending part. Usually, when you stir a fluid, it gets thicker (more viscous) or thinner (shear-thinning) in a predictable way.

• The Discovery: For these active noodles, at certain speeds, the fluid actually acts like it has negative viscosity.
• The Analogy: Imagine you are pushing a heavy box across the floor. Usually, friction resists you. But with these noodles, it's as if the box suddenly starts pushing you forward as you try to drag it. The active motors inside the noodles are pushing against the flow so hard that they actually help the fluid move, rather than resisting it. It's like the noodles are "swimming" in a way that reduces the drag on the whole system.

4. Why Does This Matter?

You might ask, "Who cares about wiggly noodles?"

• Real Life: This isn't just about plastic chains. This models real biological things like DNA, microtubules (the skeleton of cells), and even tiny worms like C. elegans.
• The Takeaway: Cells are full of these active, semi-stiff structures. They are constantly being pushed and pulled by fluid flows inside the body. Understanding how they stretch, curl, and tumble helps us understand how cells organize their DNA, how they move, and how diseases might disrupt these delicate mechanical dances.

Summary

The researchers found that when you mix self-propulsion (swimming) with shear flow (stirring), you get a unique dance that is different from anything seen in non-living matter.

1. Slow flow: Unravels the active curls.
2. Medium flow: Creates a unique, fast-tumbling motion and a "negative viscosity" where the noodles help the flow instead of hindering it.
3. Fast flow: The current wins, and the noodles act like dead sticks.

It's a reminder that in the microscopic world, "alive" things don't just get pushed around; they push back, creating entirely new rules for how fluids behave.

Technical Summary

1. Problem Statement

The study investigates the nonequilibrium structural and dynamical properties of active polar semiflexible polymers subjected to linear shear flow. While the behavior of passive polymers and flexible active polymers under shear is well-documented, the specific interplay between stiffness (semiflexibility), self-propulsion (activity), and shear flow remains less understood.

Key questions addressed include:

• How does semiflexibility alter the conformational response (stretching, folding, tumbling) compared to flexible active polymers?
• What are the scaling laws governing the polymer's dimensions and tumbling dynamics in the presence of both activity and shear?
• How does activity influence the rheological properties, specifically the shear viscosity, of semiflexible filaments?

The study is motivated by biological systems such as actin filaments, microtubules, and collective worm movements, which exhibit semiflexible characteristics and active dynamics.

2. Methodology

The authors employed numerical simulations using a coarse-grained model in two dimensions.

• Polymer Model: The polymer is modeled as a wormlike chain consisting of $N=51$ beads connected by harmonic bonds.
  • Interactions: Includes bond stretching ($U_{bond}$), bending rigidity ($U_{bend}$), and excluded volume interactions (truncated Lennard-Jones potential, $U_{ex}$).
  • Activity: Active forces are applied tangentially along the bond vectors, mimicking self-propulsion. The total active force on a bead depends on the orientation of adjacent bonds.
• Hydrodynamics: The system is immersed in a fluid modeled by Brownian Multiparticle Collision Dynamics (B-MPC). This approach introduces stochastic collisions with virtual solvent particles to generate thermal noise and friction but excludes hydrodynamic interactions (momentum is not conserved globally). This setup mimics "dry" active systems like filaments on a motility assay or gliding microorganisms.
• Simulation Parameters:
  • Dimensionless Numbers:
    • Persistence Length Ratio ($L_p/L$): 0.4 (flexible/semiflexible) and 2.0 (stiffer).
    • Péclet Number ($Pe$): Ratio of active to thermal energy ($0 \le Pe \le 2 \times 10^4$).
    • Weissenberg Number ($Wi$): Ratio of shear rate to polymer relaxation time ($1 \le Wi \le 10^4$).
  • Integration: Newton's equations of motion are integrated using the velocity-Verlet algorithm.

3. Key Contributions and Results

A. Conformational Dynamics

• Zero Shear ($Wi=0$):
  • Flexible ($L_p/L=0.4$): At high activity ($Pe \gtrsim 10^3$), polymers form long-lived spiral and folded states due to excluded volume interactions.
  • Stiff ($L_p/L=2$): Activity has a minimal effect on conformation; the polymer remains extended.
• Under Shear Flow:
  • Low Shear ($Wi \lesssim 10$): Shear disentangles the active spiral structures, leading to significant stretching along the flow direction. This effect is most pronounced for flexible polymers with high activity.
  • Intermediate Shear: The polymer undergoes tumbling motion, folding into U- or S-shaped structures (hairpins). This reduces the end-to-end distance along the flow direction.
  • High Shear ($Wi \to \infty$): Shear dominates over activity. The polymer behaves like a passive polymer, fully stretched along the flow direction.

B. Scaling Laws and Critical Exponents

The study identifies distinct scaling regimes that differ from both passive polymers and flexible active polymers:

1. Gradient Direction Shrinkage ($\langle R_{ey}^2 \rangle$):

   • Passive: Decays as $Wi^{-1/2}$.
   • Flexible Active (3D): Decays as $Wi^{-4/3}$.
   • Semiflexible Active (2D - This Study): In an intermediate, activity-dependent regime, the decay follows $\langle R_{ey}^2 \rangle \sim Wi^{-1}$. This exponent $\nu \approx 1$ is attributed to the specific conformational constraints of semiflexible polymers.

2. Tumbling Time ($\tau_\phi$):

   • Passive: Scales as $\tau_\phi \sim Wi^{-2/3}$.
   • Flexible Active (3D): Scales as $\tau_\phi \sim Wi^{-1/3}$.
   • Semiflexible Active (2D - This Study): In the intermediate regime, the tumbling time scales as $\tau_\phi \sim Wi^{-3/4}$. This steeper dependence indicates a faster tumbling dynamics driven by the interplay of stiffness and activity.

3. Alignment:

   • The alignment angle $\phi_m$ relative to the flow direction scales as $\tan(2\phi_m) \sim Wi^{-1}$ in the intermediate active regime, extending the alignment regime to higher Weissenberg numbers compared to passive polymers (which scale as $Wi^{-1/3}$).

C. Rheological Properties (Viscosity)

• Negative Viscosity: A significant finding is the emergence of negative viscosity ($\eta < 0$) for semiflexible polymers at low to intermediate shear rates and high activity ($Pe \gtrsim 10^2$).
  • In this regime, the active contribution to shear stress dominates over intramolecular forces.
  • The negative viscosity scales roughly as $\eta \sim -Wi^{-2}$.
  • This contrasts with flexible active polymers (which show positive viscosity) and "wet" active systems (where negative viscosity is usually associated with pusher-type swimmers interacting with flow fields). Here, it arises in a "dry" system due to the specific mechanics of semiflexible filaments.
• Transition: As $Wi$ increases, shear dominates, the viscosity becomes positive, and the system reverts to standard shear-thinning behavior ($\eta \sim Wi^{-1/2}$) typical of passive polymers.

4. Significance and Conclusion

This paper provides a comprehensive characterization of how stiffness modulates the response of active matter to shear flow.

• Distinct Physics of Semiflexibility: The results demonstrate that semiflexible active polymers do not simply interpolate between passive and flexible active behaviors. They exhibit unique scaling exponents ($\nu \approx 1$ for shrinkage, $\tau_\phi \sim Wi^{-3/4}$ for tumbling) driven by their ability to form compact spiral states that are subsequently unfolded by shear.
• Negative Viscosity in Dry Systems: The observation of negative viscosity in a dry, semiflexible system challenges previous assumptions that such phenomena are exclusive to hydrodynamic "wet" systems with pusher-type swimmers.
• Biological Relevance: The findings offer theoretical insights into the rheology and dynamics of cytoskeletal networks (actin/microtubules) and collective worm behaviors, suggesting that activity can significantly enhance shear-thinning and induce non-standard flow responses in biological filaments.

The study concludes that while the system asymptotically approaches passive polymer behavior at very high shear rates, the intermediate regime is dominated by a complex interplay of activity and stiffness, resulting in unique dynamical and rheological signatures.