Two-field theory for phase coexistence of active Brownian particles

Basic Information

• Paper ID: 2504.13327
• Title: Two-field theory for phase coexistence of active Brownian particles
• Authors: Pablo Pérez-Bastías, Rodrigo Soto (Department of Physics, University of Chile)
• Classification: cond-mat.soft (Soft Condensed Matter Physics)
• Publication Date: October 15, 2025 (arXiv v2)
• Paper Link: https://arxiv.org/abs/2504.13327v2

Abstract

Active Brownian particles (ABPs) serve as a minimal model of active matter systems. When ABPs are sufficiently persistent, they undergo a liquid-gas phase separation and, in the presence of obstacles, accumulate around them, forming a wetting layer. Here, we perform simulations of ABPs in a quasi-one-dimensional domain in the presence of a wall, studying the dynamics of the polarization field. Over the course of time, we observe a transition from a homogeneous (where all particles are aligned) to a heterogeneous (where particles align only at the interface) polarization regime. We propose coarse-grained equations for the density and polarization fields based on microscopic and phenomenological arguments that correctly account for the observed phenomena.

Research Background and Motivation

Research Problem

The core problem addressed in this study is how to accurately describe the dynamics of the density field and polarization field of highly persistent active Brownian particles (ABPs) during phase separation, particularly the wetting layer formation mechanism in the presence of a wall.

Problem Significance

1. Theoretical Importance: ABPs as a minimal model of active matter systems are crucial for understanding collective behavior in biological systems (such as bacteria and microalgae) and artificial active colloids
2. Practical Value: Related mechanisms have potential applications in drug delivery and self-assembling materials
3. Physical Insight: Reveals phase transition mechanisms in non-equilibrium systems where detailed balance is broken

Limitations of Existing Methods

1. Adiabatic Approximation Failure: In the high persistence limit (Dr → 0), traditional local density gradient approximations fail to capture system dynamics
2. Divergence Issues: At interfaces, density gradient approximations lead to divergence of the polarization field
3. Time Scale Separation Failure: At high persistence, the relaxation time scale of the polarization field is long and cannot be simply enslaved to the density field

Research Motivation

Based on the Dean method combined with phenomenological parameters, establish a two-field theory of density and polarization to accurately describe the polarization dynamics transition across different time scales.

Core Contributions

1. Proposed Two-field Continuum Theory: Established coupled dynamical equations for density ρ and polarization q, overcoming limitations of single-field theory at high persistence
2. Discovered Polarization Dynamics Transition: Identified temporal evolution from homogeneous polarization (all particles aligned) to heterogeneous polarization (alignment only at interfaces)
3. Established Phenomenological Closure: Obtained phenomenological expressions for tensor T^αβ through numerical simulations, achieving closure of the equation system
4. Verified Theoretical Predictions: Numerical solutions of continuum equations show excellent agreement with ABP simulations across different time scales

Methodology Details

Task Definition

Study phase separation dynamics of ABPs in a quasi-one-dimensional domain in the presence of a wall, with particular focus on spatiotemporal evolution of the polarization field.

Input: Microscopic ABP dynamical equations Output: Continuum dynamical equations for macroscopic density field ρ(x,t) and polarization field q^α(x,t) Constraints: High persistence condition (Dr ≪ v₀/σ)

Model Architecture

1. Microscopic Dynamical Equations

The equations of motion for ABPs are:

ṙᵢ = v₀nᵢ + Fᵢ
θ̇ᵢ = √(2Dr)ηᵢ(t)

where nᵢ = (cos θᵢ, sin θᵢ) is the self-propulsion direction, and Fᵢ is the inter-particle interaction force.

2. Dean Method Derivation

Define macroscopic fields:

• Density field: ρ(r,t) = Σᵢ δ(rᵢ - r)
• Polarization field: q^α(r,t) = Σᵢ nᵢ^α δ(rᵢ - r)
• Nematic tensor: Q^αβ(r,t) = Σᵢ δ(rᵢ - r)nᵢ^α nᵢ^β - δ^αβ/2

Through the Dean method, exact field equations are obtained:

∂ₜρ = -v₀∇α[q^α + G^α]                    (5)
∂ₜq^α = -Drq^α - v₀∇βT^αβ + √(2Dr)ρμ^α   (6)

3. Phenomenological Closure

Approximation for G^α:

G^α = -ρ/ρₗ q^α - ∇α δF/δρ               (7)

The first term represents effective drag, the second term is the thermodynamic contribution.

Approximation for T^αβ:

T^αβ ≈ h(ρ)δ^αβ + g(ρ)[q²/2 δ^αβ - q^αq^β]  (12)

where h(ρ) = aρ(1-ρ/ρ₀)², g(ρ) = b1-(ρ/ρ₀)²/ρ³.

4. Final Two-field Equations

∂ₜρ = -v₀∇ₓ[(1-ρ/ρₗ)q] + v₀∇²ₓ δF/δρ                    (13)

∂ₜq = -Drq - v₀∇ₓ[aρ(1-ρ/ρ₀)² - b(1-[ρ/ρ₀]²)q²/(2ρ³)] - v₀δF/δq  (14)

Technical Innovations

1. Two-field Coupling Mechanism: Unlike traditional density gradient approximations, explicitly considers independent dynamics of the polarization field
2. Phenomenological Tensor: Captures polarization-induced stress through anisotropic form of T^αβ
3. Multi-timescale Description: Unified description of both short-time homogeneous polarization and long-time heterogeneous polarization states
4. Finite Persistence Correction: Provides finite polarization solution q^x = ±√(2h(ρ)/g(ρ)) in the Dr → 0 limit

Experimental Setup

Numerical Simulation Parameters

• Number of Particles: N = 1200
• System Size: Lₓ = 410σ, Lᵧ = 10σ
• Self-propulsion Velocity: v₀ = 1
• Rotational Diffusion Coefficient: Dr = 10⁻³v₀/σ
• Persistence Length: ℓ = v₀/Dr = 1000σ
• Size Ratio: K = Lₓ/ℓ = 0.4 < 1
• Simulation Time: T = 5×10³v₀/σ

Interaction Potential

WCA-type repulsive potential is employed to prevent particle overlap, with particle diameters of σ and 1.4σ respectively to prevent crystallization.

Data Processing

Interface position x₀ is determined by fitting with hyperbolic tangent function, with all simulation results averaged after interface alignment.

Experimental Results

Main Results

1. Polarization Dynamics Transition

• t < 1/Dr: Homogeneous polarization state, all particles in wetting film oriented toward wall
• t > 1/Dr: Heterogeneous polarization state, non-zero polarization exists only at interfaces

2. Numerical Verification

Continuum equation numerical solutions show excellent agreement with ABP simulations in:

• Growth law of wetting layer width over time
• Exponential decay of polarization field ∼ exp(-t/τq)
• Steady-state density and polarization distribution profiles

3. Parameter Fitting Results

• a = 0.4, b = 0.1, ρ₀ = ρₗ = 0.77
• Free energy density: f(ρ) = α(ρ - ρₗ/2)⁴, α = 1
• Gradient coefficient: γ = 1, λ = 1

Behavior at Different K Values

• K < 1: Transition between two polarization states observed
• K > 1: Only heterogeneous polarization state appears
• Validates the theoretical prediction of K = Lₓ/ℓ classification criterion

Stress Tensor Analysis

The existence of non-variational stress tensor σ^αβ_a is key to maintaining wetting film stability, embodying the non-equilibrium nature of active systems.

Related Work

Traditional MIPS Theory

• Cahn-Hilliard-type equations consider only density field
• Effective velocity model v(ρ) = v₀(1-ρ/ρₗ)
• Adiabatic approximation of polarization field q ≈ -1/(2Dr)∇v(ρ)ρ

Existing Two-field Theories

This work extends early theories by Bialké et al. and Speck et al., particularly:

• Corrects divergence problems at high persistence
• Introduces anisotropic structure of phenomenological tensor T^αβ
• Provides unified multi-timescale description

Experimental Relevance

Resonates with active colloid experiments by Fins Carreira et al. and inertial active particle research by Caprini et al.

Conclusions and Discussion

Main Conclusions

1. Established a two-field theoretical framework describing phase separation of highly persistent ABPs
2. Successfully predicted and explained polarization field dynamics transition phenomena
3. Provided mechanism for describing phase separation without requiring double-well potential
4. Verified theory validity across different time scales

Limitations

1. Phenomenological Parameters: T^αβ expressions based on steady-state data fitting; transient behavior may require additional corrections
2. One-dimensional Restriction: Phenomenological closure verified only in quasi-one-dimensional geometry
3. Nematic Contributions: Tangential contributions of nematic tensor not fully considered
4. Noise Effects: Multiplicative noise in field equations not thoroughly investigated

Future Directions

1. Extension to more complex two-dimensional/three-dimensional geometries
2. Three-field theory including nematic field
3. Study effects of field equation noise on wetting layer detachment
4. Explore applicability under different boundary conditions

In-depth Evaluation

Strengths

1. Theoretical Innovation: First systematic establishment of two-field continuum theory for ABPs, solving fundamental defects of single-field theory
2. Phenomenological Insight: Identifies universal classification of polarization dynamics through K = Lₓ/ℓ parameter
3. Numerical Verification: Quantitative consistency between theoretical predictions and simulation results is convincing
4. Physical Intuition: Clearly elucidates microscopic mechanism of interface polarization caused by bias effects

Weaknesses

1. Parameter Dependence: Phenomenological functions h(ρ) and g(ρ) require refitting for specific systems
2. Dimensional Limitation: Theory verification limited to quasi-one-dimensional case; high-dimensional generalization requires further work
3. Transient Accuracy: Certain deviations exist in quantitative prediction of polarization decay
4. Free Energy Choice: Specific form of f(ρ) significantly affects numerical stability

Impact

1. Academic Value: Provides new theoretical framework for active matter theory, expected to inspire subsequent research
2. Practical Prospects: Provides guidance for designing functional materials based on active particles
3. Methodological Contribution: Combination of Dean method and phenomenological closure provides paradigm for modeling complex active systems

Applicable Scenarios

• Phase separation studies of highly persistent active particle systems
• Behavior prediction of active matter under complex boundary conditions
• Theoretical modeling of bacterial aggregation and biofilm formation in biological systems
• Design and optimization of active colloid materials

References

The paper cites 44 related references, primarily including:

• Foundational active matter theory (Ramaswamy 2017, Cates & Tailleur 2015)
• MIPS theory development (Bialké et al. 2013, Speck et al. 2014-2022)
• Experimental verification work (Fins Carreira et al. 2024)
• Numerical methods (Dean 1996, Burns et al. 2020)

Overall Assessment: This is a high-quality paper with significant contributions to active matter theory. By cleverly combining microscopic derivation and phenomenological modeling, it successfully establishes a two-field theory describing phase separation of highly persistent ABPs. The excellent agreement between theoretical predictions and numerical verification demonstrates method effectiveness, laying a solid foundation for subsequent research in this field.