Random Organization Criticality with Long-Range Hydrodynamic Interactions

Basic Information

• Paper ID: 2506.14330
• Title: Random organization criticality with long-range hydrodynamic interactions
• Authors: Tristan Jocteur, Cesare Nardini, Eric Bertin, Romain Mari
• Classification: cond-mat.soft cond-mat.stat-mech
• Publication Date: June 17, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2506.14330

Abstract

This paper investigates the reversible-irreversible transition (RIT) in athermal soft matter systems, a phase transition from absorbing state to active phase. By extending the random organization model (ROM) to include long-range hydrodynamic interactions, the authors find that when the interaction decay exponent is smaller than the spatial dimension, critical properties continuously depend on the decay exponent. As the interaction range increases, the transition changes from concave to convex (order parameter exponent β > 1), fluctuations transition from divergent to vanishing, and hyperuniformity at the transition disappears. The authors rationalize this critical behavior through a local mean-field model.

Research Background and Motivation

Problem Background

1. Kinetic arrest in soft matter systems: Soft materials such as foams, granular materials, and emulsions exhibit kinetic arrest under external driving, forming "absorbing states"
2. RIT in periodically sheared suspensions: Under small-amplitude oscillatory shear, the system reaches a limit cycle where dynamics freeze under stroboscopic observation; beyond critical driving amplitude, stroboscopic dynamics become diffusive
3. Limitations of ROM model: Existing random organization models neglect hydrodynamic interactions between suspended particles

Research Motivation

1. Theoretical challenges: Experiments and numerical simulations show RIT as a convex transition (β > 1), inconsistent with the conservative directed percolation (CDP) universality class predicted by ROM
2. Hydrodynamic effects: Long-range hydrodynamic interactions generated by particle contacts affect all other particles, creating new activity generation mechanisms
3. Universality questions: Understanding how different boundary conditions (bulk vs. confined systems) affect RIT properties

Core Contributions

1. Extended ROM model: Proposes a mediated ROM incorporating long-range hydrodynamic interactions with power-law decay of interaction strength
2. Continuous variation of critical exponents: Demonstrates that when decay exponent α < d, critical exponents continuously depend on α
3. Discovery of convex transitions: Proves that long-range interactions cause transitions to change from concave (β < 1) to convex (β > 1)
4. Disappearance of hyperuniformity: Reveals that long-range interactions destroy the characteristic hyperuniformity of CDP-class transitions
5. Theoretical framework: Establishes mean-field theory based on Hébraud-Lequeux-type models, explaining diffusion-induced activity generation mechanisms

Methodology Details

Task Definition

Study critical behavior of reversible-irreversible transitions in periodically sheared suspensions under long-range hydrodynamic interactions, particularly how critical exponents depend on the interaction decay exponent α.

Model Architecture

Original ROM Model

• System setup: Cubic domain with periodic boundary conditions in d dimensions, containing N spherical particles of diameter D
• Activity definition: Particles i and j are defined as active when their separation r_ < D
• Dynamics rules: Active particles receive random kicks δ_{a,i} at each time step, with typical magnitude Δ_a = D

Mediated ROM Model

Core innovation lies in introducing long-range induced motion of passive particles:

1. Spatial coarse-graining: Partition system into cubic grid of size ℓ, defining activity field: $A_b^t = \frac{1}{\ell^d}\sum_{i \in b} n_i$
2. Long-range interaction kernel: $G(r) = \frac{c}{(1 + r^2)^α}$ where α controls interaction decay rate
3. Passive particle dynamics: Passive particle i receives random kicks δ_{p,i} with amplitude: $\Delta_{p,i} = \sqrt{\sum_{b'} G(r_{b'b})A_{b'}}$
4. Boundary condition dependence:
   • Bulk systems: α = d-1
   • Confined systems: α = d

Technical Innovations

1. Spatial treatment of long-range noise: Unlike previous homogenization approaches, this work considers spatially distributed long-range interactions
2. Systematic α-dependence study: First systematic investigation of critical property dependence on interaction range
3. Multi-scale analysis: Combines particle-level dynamics with continuum field theory analysis

Experimental Setup

Numerical Simulation Parameters

• Two-dimensional systems: Maximum system size L = 8192, up to N ≈ 1×10^7 particles
• Three-dimensional systems: Maximum system size L = 768, up to N ≈ 4×10^7 particles
• Interaction range: α varies from 0.5 to 3.5
• Initial conditions: Particle positions uniformly randomly distributed

Evaluation Metrics

1. Order parameter exponent β: Mean activity ⟨A⟩ ∼ δφ^β, where δφ = φ - φ_c
2. Fluctuation exponent γ': Activity fluctuations ⟨δA²⟩ ∼ δφ^{-γ'}
3. Relaxation exponent δ: Activity decay at critical point A(t) ∼ t^{-δ}
4. Structure factor: S(q) ∼ q^ζ (q → 0), characterizing hyperuniformity

Comparison Methods

• Standard ROM: Short-range limit α → ∞
• CDP universality class: Theoretically expected critical exponents
• MROM∞ model: Mean-field version with infinite jumps of active particles
• Lévy-SHL model: Theoretical model based on anomalous diffusion

Experimental Results

Main Results

Evolution of Critical Exponents

Key Findings

1. Convexity transition threshold:
   • d=2: α ≈ 1.5
   • d=3: α ≈ 2.5
2. CDP stability boundary: CDP behavior maintained when α ≳ d
3. Fluctuation property reversal: γ' transitions from positive values (divergent fluctuations) to negative values (vanishing fluctuations), synchronized with convexity emergence

Hyperuniformity Analysis

• CDP region: Exhibits typical q^{0.45} scaling behavior
• Long-range interaction region: Hyperuniformity gradually disappears, completely vanishing at α ≲ 1.5
• Structure factor: Plateau appears at small q, indicating recovery of normal fluctuations

Theoretical Verification

Verification through MROM∞ model confirms mean-field picture:

• Measured pseudo-gap exponent θ and order parameter exponent β satisfy relation β = 2θ (when 1 < 2θ < 2)
• Confirms 2P → 2A channel dominates transition at small α

Related Work

Absorbing Phase Transitions in Soft Matter

1. Yield transitions: Similar convex behavior in elastoplastic models
2. Active matter: Freezing transitions in conditionally active systems
3. Granular materials: Reversible-irreversible transitions in dry granular systems

Critical Phenomena with Long-Range Interactions

1. Equilibrium systems: Fisher-Ma-Nickel theoretical framework
2. Non-equilibrium systems: Lévy flights and anomalous transport
3. Directed percolation: Corrections to CDP from long-range transport

Hydrodynamic Interactions

1. Stokes flow: Long-range velocity fields at low Reynolds numbers
2. Boundary effects: Hydrodynamic corrections in confined geometries
3. Colloidal suspensions: Non-CDP behavior observed in experiments

Conclusions and Discussion

Main Conclusions

1. Breaking of universality: Long-range hydrodynamic interactions break CDP universality class, leading to new critical behavior
2. Geometric dependence: Bulk (α = d-1) and confined (α = d) systems belong to different universality classes
3. Mechanism switching: Transition from A+P → 2A (CDP mechanism) to 2P → 2A (diffusion-induced mechanism)
4. Experimental predictions:
   • Bulk systems: Convex transitions, no hyperuniformity
   • Confined systems: CDP behavior, preserved hyperuniformity

Limitations

1. Model simplifications: Neglects temporal evolution of contacts during shear cycles
2. Interaction form: Employs simplified power-law decay; actual hydrodynamics is more complex
3. Finite-size effects: Some regimes may be limited by system size
4. Missing continuum theory: Complete continuum field theory description remains undeveloped
5. Experimental validation: Primarily based on numerical simulations; more experimental verification needed

Future Directions

1. Continuum theory development: Establish field theory incorporating long-range induced diffusion
2. Experimental verification: Validate theoretical predictions under different boundary conditions
3. Unified description: Seek unified theoretical framework for RIT and yield transitions
4. Application extensions: Explore applications in hyperuniform structure preparation

In-Depth Evaluation

Strengths

1. Theoretical innovation: First systematic study of hydrodynamic interaction effects on RIT, revealing new physical mechanisms
2. Numerical precision: Large-scale simulations provide high-precision critical exponent measurements
3. Multi-scale analysis: Successfully connects microscopic particle dynamics with macroscopic critical phenomena
4. Experimental relevance: Resolves long-standing theory-experiment discrepancies

Weaknesses

1. Model assumptions: Certain simplifications (e.g., uncorrelated random displacements) may affect quantitative results
2. Theoretical completeness: Lacks complete continuum field theory; relies primarily on mean-field analysis
3. Insufficient experimental validation: Primarily based on numerical simulations; requires more experimental verification

Impact

1. Academic value: Provides important supplement to soft matter phase transition theory, challenging existing universality assumptions
2. Practical significance: Provides predictive framework for suspension behavior under different geometric conditions
3. Cross-disciplinary connections: Establishes deep links between RIT and yield transitions

Applicable Scenarios

1. Colloidal suspensions: Predicting shear response under different boundary conditions
2. Microfluidic devices: Designing rheological behavior in confined geometries
3. Material preparation: Guidance for preparing hyperuniform structures via phase transitions
4. Theoretical research: Provides paradigm for studying other long-range interaction systems

References

This paper cites 79 important references spanning soft matter physics, phase transition theory, hydrodynamics, and other fields, providing solid theoretical foundation for understanding critical phenomena under long-range interactions.

Overall Assessment: This is a high-quality theoretical physics paper that, through clever model design and precise numerical computation, reveals profound effects of long-range hydrodynamic interactions on soft matter phase transitions. The paper not only resolves long-standing theory-experiment discrepancies but also opens new research directions in related fields. Despite room for improvement in theoretical completeness, its innovation and significance make it an important contribution to soft matter physics.