🔬 condensed matter

How Phase Coexistence affects the mechanical properties of heterogeneous 2D suspensions

This study demonstrates that phase coexistence in two-dimensional Lennard-Jones suspensions significantly alters their rheological properties, establishing a comprehensive map where the interplay of density, temperature, and coexisting phases dictates the transition from viscous to elastic-dominated behavior under shear.

Original authors: Raul Molina-Prados Lallena, Jose Martin-Roca, Kristian Thijssen, Tyler Shendruk, Angelo Cacciuto, Chantal Valeriani

Published 2026-02-24

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Raul Molina-Prados Lallena, Jose Martin-Roca, Kristian Thijssen, Tyler Shendruk, Angelo Cacciuto, Chantal Valeriani

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are making a giant batch of soup. Usually, when scientists study how liquids flow (rheology), they look at smooth, uniform soups—like a perfectly blended tomato soup where every spoonful is exactly the same.

But in the real world, things are rarely that perfect. Sometimes you have a soup with big chunks of vegetables floating in clear broth, or a mixture where ice crystals are forming in a slushy drink. This paper explores exactly that: what happens to the "flow" of a material when it's a messy mix of different states at the same time.

The researchers used a computer to simulate a 2D world filled with tiny, sticky balls (particles). They watched what happened when these balls were in a state of "phase coexistence"—meaning some were acting like a gas (flying around loosely), some like a liquid (huddling together), and some like a solid (stuck in a rigid grid).

Here is the breakdown of their findings using everyday analogies:

1. The "Perfect Soup" vs. The "Chunky Stew"

The Homogeneous State (The Perfect Soup): When the temperature is high, the particles move freely and uniformly. If you stir this soup, it flows easily. It behaves like a standard fluid (Newtonian). It has one "relaxation time," meaning if you stop stirring, it settles down quickly and predictably.
The Coexistence State (The Chunky Stew): When they cooled the soup down, it didn't just get thicker; it split. You had pockets of gas, pockets of liquid, and eventually, pockets of solid.
- The Discovery: Even though it was a mix, at first, it still acted like a simple fluid. But as it got colder and the "solid chunks" grew, the behavior changed drastically.

2. The "Traffic Jam" Effect

The researchers found that when the system enters a Gas/Solid coexistence (like a crowd where some people are running free, but a large group has formed a rigid, unmoving wall):

Stress gets stuck: Imagine trying to push a crowd. If everyone is moving, you can push through. But if a large group has formed a solid wall, your push gets "stuck" against them.
The Result: The material stops acting like a liquid and starts acting like a spring. Instead of flowing away, it stores the energy you put into it (elasticity). The "solid" parts act like a skeleton that holds the whole system together, making it incredibly hard to move.

3. The "Shear Thinning" Surprise

One of the most interesting findings is about Shear Thinning.

The Analogy: Think of ketchup. When it sits in the bottle, it's thick and won't move. But if you shake the bottle hard (apply force), it suddenly becomes runny and flows out.
The Finding: This "sticky" mixture of gas and solid particles behaves exactly like ketchup. No matter how cold or "solid" it gets, if you push it hard enough (apply high shear), it suddenly becomes easier to flow.
Why? The force of your push breaks up the rigid "solid" structures, forcing the particles to line up and slide past each other, temporarily turning the solid back into a fluid.

4. Why This Matters

Scientists often study "active matter" (like bacteria or self-driving robots) to understand why complex mixtures behave strangely. They thought these weird behaviors were caused by the particles being "alive" or moving on their own.

The Big Twist: This paper shows that you don't need living, moving particles to get these complex behaviors. You just need a mixture of states (coexistence). Even a passive, "dead" system will act weirdly and become elastic or shear-thinning just because it's a mix of gas, liquid, and solid all at once.

Summary

Think of this research as a map for navigating a landscape of materials.

Old Map: "If it's hot, it's a liquid. If it's cold, it's a solid."
New Map: "If it's a mix of hot and cold states, it becomes a super-material that acts like a spring when you push it gently, but turns into a runny liquid when you push it hard."

This helps engineers and scientists design better materials (like new paints, gels, or biological tissues) by understanding that the mixture of states is just as important as the temperature or density itself.

1. Problem Statement

While the rheology of homogeneous colloidal suspensions (liquid, crystalline, or glassy states) is well-studied, the mechanical properties of heterogeneous systems where distinct phases coexist (e.g., gas/liquid or gas/solid) remain poorly understood.

The Gap: Most rheological investigations assume spatial uniformity. However, in systems with attractive interactions (like Lennard-Jones fluids), thermodynamic conditions can lead to phase separation. It is unclear how the coexistence of domains with different densities and mobilities alters stress transmission, relaxation pathways, and the transition from viscous to elastic behavior.
The Specific Question: How does the presence of thermodynamic phase coexistence (specifically gas/liquid and gas/solid) modify the viscoelastic response and viscosity of a 2D passive system compared to its homogeneous counterparts?

2. Methodology

The authors employed two-dimensional Brownian dynamics simulations to model a passive, thermalized system of attractive particles.

System Model:
- Particles: $N = 5000$ disks of diameter $\sigma$ interacting via a truncated Lennard-Jones (LJ) potential ( $U_{ij}$ ) with a cutoff at $r_c = 2.5\sigma$ .
- Dynamics: Governed by the Langevin equation (overdamped limit) including friction ( $\gamma$ ) and thermal noise.
- Ensemble: NVT (constant number of particles, volume, and temperature) with square periodic boundary conditions.
Phase Space Exploration:
- Simulations spanned the gas/liquid coexistence region and the gas/solid coexistence region (below the triple point).
- Parameters varied: Density ( $\rho$ ) from 0.3 to 0.7 and Temperature ( $T$ ) from 0.25 to 0.55 (in LJ reduced units).
Rheological Analysis Tools:
- Linear Response (Equilibrium): Used the Green-Kubo relation to compute the shear relaxation modulus $G(t)$ from the autocorrelation of the stress tensor ( $\sigma_{xy}$ ).
- Frequency Domain: Applied Fourier transforms to $G(t)$ to obtain the complex shear modulus $G^*(\omega) = G'(\omega) + iG''(\omega)$ , where $G'$ is the storage modulus (elastic) and $G''$ is the loss modulus (viscous).
- Effective Viscosity: Calculated $\eta(\omega)$ from the moduli.
- Nonlinear Response (Steady Shear): Performed constant strain-rate simulations to measure steady-state shear stress $\langle \sigma_{xy} \rangle$ and fit the data to the Herschel-Bulkley model ( $\sigma = \sigma_0 + k\dot{\gamma}^n$ ) to determine the shear-thinning exponent $n$ .

3. Key Results

A. Relaxation Dynamics and Moduli

Homogeneous/High-T Regime: At high temperatures (above the binodal) and low densities, the system exhibits a single-exponential decay in $G(t)$ . The storage ( $G'$ ) and loss ( $G''$ ) moduli cross over at a specific frequency $\omega^\dagger$ , characteristic of a Maxwell fluid (viscous at low frequencies, elastic at high frequencies).
Gas/Liquid Coexistence: Even within the gas/liquid coexistence region, the system largely retains Maxwell-like behavior with a single dominant relaxation timescale, though the crossover frequency shifts to lower values as temperature decreases (slowing dynamics).
Gas/Solid Coexistence (Below Triple Point):
- Multi-scale Relaxation: $G(t)$ deviates from single-exponential decay, exhibiting two distinct timescales (fast and slow).
- Elastic Dominance: The crossover between $G'$ and $G''$ disappears within the accessible frequency window; $G'$ consistently dominates $G''$ , indicating a transition to elastic-dominated behavior.
- Non-monotonic Loss Modulus: $G''(\omega)$ develops a local maximum and minimum at low temperatures, unlike the single peak seen in pure fluids. This suggests complex viscous dissipation mechanisms due to structural heterogeneity.

B. Viscosity and Shear Thinning

Effective Viscosity ( $\eta(\omega)$ ):
- High T: A clear low-frequency Newtonian plateau exists.
- Low T (Coexistence): The Newtonian plateau vanishes. Viscosity continues to rise as frequency decreases, indicating long-lived stress correlations caused by extended solid domains that hinder structural relaxation.
Steady Shear Behavior:
- All systems exhibit shear thinning (exponent $n < 1$ ) at high shear rates.
- Below the Triple Point: The Herschel-Bulkley exponent $n$ remains less than unity for all shear rates, confirming persistent shear thinning in gas/solid coexistence phases.
- Consistency: The frequency-dependent viscosity at high frequencies matches the steady-state viscosity at high shear rates, confirming that shear thinning is a robust feature driven by the competition between imposed deformation and internal relaxation.

4. Key Contributions

Mapping Phase Coexistence to Rheology: The study establishes an extended map showing that phase coexistence (not just density or temperature) is a primary driver of rheological transitions. It defines how the system moves from viscous to elastic-dominated behavior as it enters gas/solid coexistence.
Distinction from Active Matter: The authors compare their passive results to recent findings on Motility-Induced Phase Separation (MIPS) in active matter. They demonstrate that complex rheological features (like non-Newtonian behavior and elastic dominance) can arise purely from thermodynamic phase coexistence in passive systems, without the need for self-propulsion.
Identification of Relaxation Mechanisms: The work identifies that the emergence of multiple relaxation timescales and the suppression of the Newtonian plateau are direct signatures of the percolation of solid domains in 2D, which fundamentally alters stress transmission.

5. Significance

Benchmark for Heterogeneous Fluids: This paper provides a critical benchmark for understanding fluids where complex dynamic states coexist, a scenario common in soft matter, biological tissues, and industrial suspensions.
Re-evaluation of "Active" Features: It suggests that some rheological complexities previously attributed solely to "activity" in biological or synthetic active matter may actually stem from general mechanisms of structural heterogeneity and phase coexistence.
Material Design: The findings offer a framework for predicting the mechanical stability and flow properties of 2D materials (e.g., colloidal monolayers, lipid membranes) under varying thermodynamic conditions, aiding in the design of materials with tunable viscoelastic properties.

In summary, the paper demonstrates that phase coexistence fundamentally alters the mechanical landscape of 2D suspensions, inducing a transition from simple viscous flow to complex, elastic-dominated, shear-thinning behavior driven by the formation of long-lived solid-like structures.