---
title: Chiral Lennard-Jones Fluid Model
url: https://www.emergentmind.com/topics/chiral-lennard-jones-fluid-model
type: topic
---

# Chiral Lennard-Jones Fluid Model

The chiral Lennard-Jones (LJ) fluid model is a two-dimensional statistical mechanical framework for assemblies of Brownian disks that interact via both standard LJ potentials and a non-conservative, pairwise transverse (chiral) force, designed to mimic the collective behavior of colloidal particles rotating at a prescribed angular speed. This model exhibits a rich phase diagram, including gas-liquid phase separation, interface-driven edge currents, and pronounced chirality-induced melting of solid phases, thus providing a systematic theoretical and computational foundation for understanding chiral particle ensembles [2307.03528].

## 1. Model Definition

The system is composed of $N$ two-dimensional Brownian disks, each of diameter $\sigma_d$ and mass $m$. The particles' interactions are governed by two core contributions:

- **Lennard-Jones Potential:**  
  The interaction is truncated at $r_c = 5\sigma_d$ and shifted by $\epsilon$ for continuity:
  $$
  U(r) = \epsilon \left[ \left( \frac{\sigma_d}{r} \right)^{12} - \left( \frac{\sigma_d}{r} \right)^6 \right] + \epsilon, \quad r < r_c
  $$
- **Transverse (“Chiral”) Force:**  
  Each pair within $r_c$ experiences a short-range, non-conservative force:
  $$
  \mathbf{f}_{ij} = \zeta \left( \frac{\omega}{\omega_0} \right) \hat{z} \times \frac{(\mathbf{r}_i - \mathbf{r}_j)} {|\mathbf{r}_i - \mathbf{r}_j|^3}
  $$
  where $\omega$ is the spin rate, $\hat{z}$ is normal to the plane, $\omega_0 = 1/\tau$ sets the unit of time with $\tau = \sqrt{m \sigma_d^2 / \epsilon}$, and $\zeta = m \sigma_d^2 / \tau^3$.

- **Overdamped Langevin Dynamics:**  
  In the overdamped regime ($\Gamma \gg m/\tau$), the equations of motion are:
  $$
  m\ddot{\mathbf{r}}_i + \Gamma \dot{\mathbf{r}}_i = -\sum_{j \ne i} \nabla_i U(|\mathbf{r}_i - \mathbf{r}_j|) + \sum_{j \ne i} \mathbf{f}_{ij} + \sqrt{2\Gamma k_B T} \, \nu_i(t)
  $$
  where $\Gamma$ is the friction coefficient and $\nu_i(t)$ is Gaussian white noise with standard correlations.

**Key Control Parameters:**

| Parameter                | Definition                                                  | Meaning                                         |
|--------------------------|-------------------------------------------------------------|-------------------------------------------------|
| $T^* = \frac{k_BT}{\epsilon}$   | Reduced temperature                               | Thermal energy scale                            |
| $\phi = \frac{\pi N}{4V}$      | Surface fraction                                   | Packing density                                 |
| $\Omega = \omega \tau$         | Dimensionless chirality                            | Chirality strength                              |
| $Pe_R = \frac{\omega \sigma_d^2}{D_t}$ | Rotational Péclet number ($D_t = \frac{k_BT}{\Gamma}$) | Chirality vs. Brownian diffusion         |
| $\chi = \zeta \Omega$          | Transverse-to-conservative force ratio             | Relative non-conservative forcing               |

## 2. Thermodynamic Framework and Phase Coexistence

Despite explicit violation of time-reversal and parity by the transverse force, the system admits a well-defined mechanical equation of state. The pressure is computed via the Irving–Kirkwood (IK) stress tensor:
$$
\sigma^{ab} = -\frac{1}{V} \left[ \sum_i m v_i^a v_i^b + \frac{1}{2} \sum_{i \ne j} r_{ij}^a F_{ij}^b \right], \quad F_{ij} = -\nabla_i U(r_{ij}) + f_{ij}
$$
The scalar pressure is $P = -\text{Tr}\, \sigma/2$.

Simulation yields $P(\phi)$ curves with pronounced van der Waals–type (“Mayer–Wood”) loops. Coexistence (binodal) densities and pressures can be quantitatively located via the Maxwell equal-area construction, with the rule:
$$
\int_{\phi_\ell}^{\phi_g} [P(\phi) - P_\text{coex}]\, d(1/\phi) = 0
$$
This construction, although equilibrium-based, remains applicable to the out-of-equilibrium chiral system and matches direct histogram-based density measurements.

For finite-size systems, the interface free-energy excess $\Delta F = \gamma L_y$ shifts the area of the pressure loop:
$$
A = \oint P(\phi)\, d(1/\phi) = 2\gamma L_y / V
$$
where $\gamma$ denotes the surface tension and $L_y$ the interface length. Increasing chirality ($\Omega$) yields a larger $A$ and thus a higher $\gamma$.

## 3. Surface Tension and Interface Structure

Surface tension is characterized using spatially resolved stress tensor components for a slab with gas–liquid interface along $y$:
$$
\gamma = -\frac{1}{2}\int dx\, [\sigma^{xx}(x) - \sigma^{yy}(x)]
$$
Numerical evaluation shows that $\gamma(\Omega)$ increases monotonically with chirality. Edge currents, generated by the transverse force, enhance tangential momentum transfer at the interface, amplifying $\sigma^{yy}$ and thus the normal–tangential stress imbalance.

These features indicate that the non-conservative chiral drive not only shifts thermodynamic coexistence but also fundamentally alters interfacial mechanics, resulting in stiffer interfaces and novel flow behavior.

## 4. Edge Currents and Rotational Viscosity

At liquid-gas coexistence, the interface sustains persistent edge currents: droplets or strips of chiral liquid display a unidirectional flow along their perimeters. The radial profile of the azimuthal velocity $v_t(r)$ near the edge follows
$$
v_t(r) \simeq v_e e^{-(R - r)/\delta}, \quad \text{for}\ r \lesssim R
$$
where $R$ is the droplet radius and $\delta = \sqrt{\frac{\eta_0 + \eta_R}{\mu}}$ is the edge current penetration depth; here, $\eta_0$ is shear viscosity, $\eta_R$ is rotational viscosity, and $\mu = \Gamma$ is substrate friction. The edge velocity magnitude is
$$
v_e \simeq \frac{2\omega \delta \eta_R}{\eta_0 + \eta_R}
$$

Rotational viscosity $\eta_R$ is central to emergent chiral hydrodynamics. It is determined microscopically either by inverting the steady-state velocity profile,
$$
\eta_R = \frac{\mu \delta v_e}{2\omega}
$$
or directly from the stress tensor in a homogeneous chiral liquid,
$$
\eta_R = \frac{\sigma^{xy}}{2\omega}
$$
Both estimates yield consistent values within approximately $30\%$.

## 5. Chirality-Induced Melting and Hexatic Patch Formation

At reduced temperatures $T^* \lesssim 0.5$, the equilibrium 2D LJ system forms a solid with quasi–long-range orientational (hexatic) order. Introduction of the chiral force induces internal stresses that disrupt global crystalline order above a critical chirality $\Omega_c(\phi, T)$. The local hexatic order parameter is
$$
\psi_{6,i} = \frac{1}{n_i} \sum_{k \in \partial_i} e^{i6\theta_{ik}}
$$
where the sum is over nearest neighbors. The global and spatial correlations are:
$$
|\Psi_6| = \frac{1}{N} \sum_{i=1}^N |\psi_{6,i}|, \qquad G_6(r) = \langle \psi_6(r) \psi_6(0) \rangle
$$
In the chiral solid, $G_6(r)$ decays exponentially, $G_6(r) \sim e^{-r/\xi}$, with hexatic correlation length $\xi \sim \Omega^{-1}$. This $\Omega$-dependence indicates that stronger spinning fragments the solid into finer “hexatic patches,” each retaining local sixfold order but exhibiting swirling dynamics.

Defects localize at the boundaries of these rotating patches, and particle displacements predominantly occur along domain walls, generating a mosaic of dynamically rotating hexatic domains that fluidize the dense phase without loss of local orientational order.

## 6. Theoretical Significance and Outlook

The chiral Lennard-Jones fluid model provides a statistically robust framework for non-equilibrium chiral fluids, revealing that:

- Phase coexistence can be captured by equilibrium-inspired thermodynamic constructs augmented by interface stress contributions.
- Surface tension increases monotonically with applied chirality due to persistent edge currents driven by non-conservative transverse forces.
- The emergence of rotational viscosity $\eta_R$ is intrinsic and quantifiable both at the hydrodynamic and microscopic levels.
- Chirality melts otherwise stable 2D solids into dynamic hexatic patchworks, with the characteristic domain size set by the inverse of chirality.

These insights lay the foundation for advanced theoretical treatment and simulation of chiral particle assemblies, bridging equilibrium and non-equilibrium statistical mechanics in active matter and soft condensed phases [2307.03528].

Source: https://www.emergentmind.com/topics/chiral-lennard-jones-fluid-model