Flow-history-dependent orientational relaxation in dilute polydisperse colloidal rod suspensions

This study demonstrates that in dilute polydisperse colloidal rod suspensions, the orientation relaxation time following flow cessation is not a fixed material constant but is flow-history-dependent, systematically decreasing with higher pre-shear rates as the dominant relaxing sub-population shifts from longer to shorter rods.

The Gist

Imagine a crowded dance floor filled with thousands of tiny, rigid sticks (like toothpicks) floating in a thick, sugary syrup. This is essentially what scientists call a colloidal suspension of cellulose nanocrystals (CNCs). These "sticks" are actually microscopic wood fibers, and they are the building blocks for many advanced materials, from 3D printing inks to smart food packaging.

The paper you shared investigates a fascinating question: How do these sticks behave when you spin them around and then suddenly stop?

Here is the breakdown of the study using simple analogies:

1. The Setup: The "Sticky" Dance Floor

The researchers put these microscopic sticks into a special machine called a Taylor-Couette cell. Think of this as two large cylinders, one inside the other, with a tiny gap between them. The inner cylinder spins, dragging the syrup and the sticks along with it.

• The Problem: In the real world, these sticks aren't all the same size. Some are long (like a ruler), and some are short (like a matchstick). This is called polydispersity.
• The Physics: Long sticks are heavy and sluggish; they take a long time to turn around. Short sticks are light and zippy; they spin and reorient very quickly.

2. The Experiment: The "Spin and Stop"

The researchers did two things:

1. Spun it: They rotated the cylinder at different speeds (shear rates). This forced the sticks to line up, like soldiers marching in formation.
2. Stopped it: They abruptly halted the cylinder.

When the spinning stops, the sticks don't just snap back to being random instantly. They have to "relax" back to a messy, random state. The researchers used a special high-speed camera with polarized light (like 3D glasses) to watch this relaxation happen in real-time. They measured birefringence, which is basically how much the light gets "twisted" by the aligned sticks. The more aligned the sticks, the more the light twists.

3. The Big Discovery: It's Not Just About Size, It's About History

In a perfect world where all sticks are the same size, the time it takes for them to relax back to random would be a fixed constant. It would be like a pendulum swinging back to rest: it always takes the same amount of time.

But this isn't a perfect world. The researchers found something surprising: The relaxation time changes depending on how fast you spun them before stopping.

• Slow Spin: If you spin the machine slowly, only the long, lazy sticks get a chance to line up. The short, zippy ones are too busy jiggling around to get organized. When you stop, the system relaxes slowly because it's dominated by those sluggish long sticks.
• Fast Spin: If you spin the machine very fast, you force both the long and the short sticks to line up. When you stop, the short sticks snap back to random almost instantly. Because the short sticks are now part of the "aligned group," the whole system relaxes much faster.

The Analogy:
Imagine a group of people trying to stand in a straight line.

• If you give them a gentle nudge (low speed), only the tall, slow-moving people manage to get in line. When you tell them to scatter, they take a long time to leave the line.
• If you give them a hard shove (high speed), even the short, fast people get pushed into the line. When you tell them to scatter, the short people run away immediately, making the whole group disperse much faster.

4. The "Magic Number" (The Peclet Number)

The researchers created a mathematical model (a Fokker-Planck equation) to predict this behavior. They found that if you look at the data through the lens of a specific number called the Péclet number (which compares the speed of the spin to the natural wiggling of the sticks), everything lines up perfectly.

It's like realizing that whether you are running a 100-meter dash or a marathon, your performance depends on your speed relative to your own energy level. Once you account for that, the data from all different stick sizes and spinning speeds collapses into a single, predictable curve.

5. Why Does This Matter?

This study proves that for materials made of mixed-size particles, you cannot just say "this material relaxes in 5 seconds." The relaxation time depends on how the material was treated before.

• For Engineers: If you are designing a 3D printer or a new type of plastic film, you can't just assume the material will behave the same way every time. You have to consider the "flow history." Did the material get squeezed fast or slow before it set? That history dictates how the final product will look and behave.
• For Scientists: It solves a long-standing puzzle about why relaxation times in messy, mixed suspensions seem to change. It's not that the material is changing; it's that the mix of particles contributing to the signal changes based on the flow.

Summary

In short, the paper shows that in a crowd of mixed-size sticks, the speed at which you spin them decides which sticks get to join the "alignment party."

• Slow spin? Only the long, slow sticks join. Relaxation is slow.
• Fast spin? Everyone joins, including the short, fast sticks. Relaxation is fast.

The "relaxation time" isn't a fixed property of the material; it's a memory of the flow. This helps scientists better predict and control the behavior of complex fluids in everything from industrial manufacturing to biological systems.

Technical Summary

1. Problem Statement

Rod-like colloidal particles, such as cellulose nanocrystals (CNCs), are fundamental to the optical and mechanical properties of soft materials. Their behavior under flow is governed by orientational dynamics. While the dynamics of monodisperse (uniform size) suspensions are well-understood, polydisperse suspensions (containing rods of varying lengths) present a significant challenge.

• The Core Issue: In polydisperse systems, rods of different lengths possess distinct rotational diffusion timescales ($D_r \propto l^{-3}$). It remains unclear how the pre-shear rate (flow history) prior to flow cessation determines which sub-populations of rods are aligned and, consequently, how this selection influences the subsequent orientational relaxation time.
• The Gap: Previous studies often treated relaxation time as an intrinsic material constant or struggled to isolate flow history effects in non-uniform flows where steady alignment was not achieved. The relationship between the applied shear rate, the specific rod sub-population contributing to the signal, and the resulting relaxation dynamics in a dilute polydisperse system was not quantitatively established.

2. Methodology

The authors employed a combined experimental and theoretical approach to isolate the effects of flow history on relaxation.

Experimental Setup

• Material: Dilute suspensions of Cellulose Nanocrystals (CNCs) in glycerol-water mixtures. CNCs were characterized via Atomic Force Microscopy (AFM) to determine their length ($l$) and diameter ($d$) distributions.
• Flow Geometry: A narrow-gap Taylor-Couette cell (inner radius 10 mm, gap 1 mm) was used to impose a well-controlled, steady simple shear flow. This geometry approximates planar Couette flow, ensuring a uniform shear rate across the gap.
• Measurement Technique: High-speed polarization imaging was used to measure flow-induced birefringence ($\Delta n$) and the orientation angle ($\phi$).
  • Steady State: Samples were sheared at various rates ($\dot{\gamma}$) until steady alignment was reached.
  • Relaxation: The flow was abruptly stopped ($t=0$), and the decay of birefringence was recorded at high frame rates (up to 250 fps) to capture the relaxation dynamics.
• Control: The study utilized five different concentrations and solvent viscosities to span a wide range of Péclet numbers ($Pe = \dot{\gamma}/D_r$).

Theoretical Framework

• Model: A polydisperse Fokker-Planck equation was formulated to describe the orientational probability density function ($\psi$) of rigid Brownian rods.
• Polydispersity Handling: The model incorporated the experimentally measured lognormal length distribution ($f_0(l)$). To account for how different lengths contribute to the optical signal, the authors introduced a moment-weighted distribution $f_n(l) \propto l^n f_0(l)$, where the parameter $n$ controls the emphasis on longer rods.
• Key Parameters:
  • Effective Rotational Diffusion ($D_r^{eff}$): Defined based on a weighted mean length $\langle l \rangle_n$.
  • Averaged Relaxation Time ($\bar{\tau}_b$): Defined as the integral of the normalized birefringence decay curve, providing a single scalar metric for the multi-modal relaxation process.

3. Key Contributions

1. Quantitative Link Between Flow History and Relaxation: The study demonstrates that the apparent relaxation time in polydisperse suspensions is not an intrinsic material constant. Instead, it is a flow-history-dependent property determined by which rod sub-populations are aligned by the pre-shear.
2. Validation of Moment-Weighted Scaling: The authors identified that a specific moment-weighted distribution ($n=2$) provides the most consistent scaling parameter to collapse both steady-state and transient relaxation data across different concentrations and shear rates.
3. Mechanism of Sub-population Selection: The work elucidates the mechanism by which increasing the pre-shear rate shifts the dominant contribution to the birefringence signal from longer rods (at low $Pe$) to shorter rods (at high $Pe$). Since shorter rods diffuse faster, this shift leads to a systematic decrease in the observed relaxation time.

4. Key Results

Steady-State Alignment

• Order Parameter ($S$) and Angle ($\phi$): Both experimental data and theoretical predictions showed that as the shear rate increases, rods align more strongly with the flow direction ($S \to 1$, $\phi \to 0^\circ$).
• Data Collapse: When plotted against the Péclet number based on the effective diffusion coefficient ($Pe = \dot{\gamma}/D_r^{eff}$ with $n=2$), the experimental data for different concentrations collapsed onto a single master curve. This confirmed that a single effective length scale can describe the competition between shear alignment and Brownian diffusion in polydisperse systems.

Relaxation Dynamics (Post-Cessation)

• Non-Exponential Decay: Unlike monodisperse systems which exhibit single-exponential decay, polydisperse suspensions showed multi-modal relaxation.
• Flow-Dependent Relaxation Time: The averaged relaxation time ($\bar{\tau}_b$) decreased systematically as the pre-shear rate increased.
  • Low Shear: Only long rods align; they have slow diffusion, leading to slow relaxation.
  • High Shear: Shorter rods also align; they have fast diffusion, leading to rapid relaxation.
• Theoretical Agreement: The polydisperse Fokker-Planck model, parameterized with the $n=2$ weighting, accurately predicted the shape of the decay curves and the trend of $\bar{\tau}_b$ versus $Pe$. The dimensionless relaxation time ($\bar{\tau}_b D_r^{eff}$) collapsed onto a universal curve distinct from the monodisperse prediction (which is constant at $1/6$).

5. Significance and Implications

• Rheo-Optical Interpretation: This study provides a quantitative framework for interpreting flow birefringence in unsteady flows. It warns that assuming a constant relaxation time or applying the Stress-Optic Law (SOL) directly to transient birefringence in polydisperse systems can lead to significant errors due to hysteresis and flow-history dependence.
• Material Characterization: The ability to link the relaxation dynamics directly to the rod length distribution and flow conditions offers a potential method for characterizing polydispersity or monitoring flow-induced microstructure changes in real-time.
• Theoretical Advancement: By moving beyond the "effective medium" approximation to a rigorous polydisperse Fokker-Planck treatment, the authors show that the "apparent" relaxation time is an emergent property of the flow-selected ensemble, fundamentally changing how relaxation in complex fluids should be modeled.

In conclusion, the paper establishes that in dilute polydisperse rod suspensions, relaxation is a flow-history-dependent ensemble property, not a fixed material constant. The pre-shear rate acts as a filter, selecting specific rod sub-populations that dictate the subsequent relaxation speed.