Chapman-Enskog expansion for chirally colliding disks

This paper demonstrates that a two-dimensional fluid of hard disks with chiral, energy- and momentum-conserving collisions obeys the H-theorem despite breaking time-reversal symmetry, and derives analytical expressions for its transport coefficients via Chapman-Enskog expansion that are validated by molecular dynamics simulations.

The Gist

Imagine a crowded dance floor where everyone is dancing to the same beat, bumping into each other, and moving around. In the world of physics, this is a gas. Usually, we think of these gas particles as perfect, round billiard balls. When two billiard balls hit, they bounce off each other in a way that is perfectly symmetrical: if you played the movie backward, it would look exactly the same. This is called "time-reversal symmetry."

But what if the dancers weren't perfect spheres? What if they were slightly lopsided, or if they had a "handedness" (like a left hand vs. a right hand)?

This paper explores a world where gas particles are like chiral dancers. They are still round disks, but they have a secret rule: they treat left-handed collisions differently than right-handed ones.

Here is the story of their discovery, broken down into simple concepts:

1. The "Lopsided" Dance Floor

Usually, when two round balls collide, it doesn't matter if they hit from the left or the right; the result is the same. In this study, the scientists imagined a gas where the particles have a tiny "bias."

Think of it like a turnstile at a subway station.

• If you approach the turnstile from the left, it lets you through easily.
• If you approach from the right, it jams a little bit, or maybe it lets you through at a slightly different angle.

The particles in this gas don't actually change shape. They are still perfect circles. But their collision rule is biased. A collision that happens "clockwise" feels slightly different than one that happens "counter-clockwise." This breaks the symmetry of the dance floor.

2. The Big Surprise: Order from Chaos

In physics, there is a famous rule called the H-theorem. It basically says: "If you have a messy, chaotic system, it will eventually settle down into a calm, orderly state (equilibrium)."

Usually, this rule relies on the idea that the laws of physics look the same whether you play the movie forward or backward. Since our "chiral dancers" break this rule (the movie backward looks different), physicists thought, "Uh oh, maybe these particles will never settle down. Maybe they'll just spin in chaos forever."

The paper's big discovery: Even though the particles are biased and "chiral," the H-theorem still works!
They proved mathematically that even with this weird bias, the gas still finds a calm, happy equilibrium. It's like a dance floor where everyone has a favorite spin direction, but eventually, everyone still settles into a steady rhythm.

3. The "Odd" Viscosity (The Mysterious Glue)

When you stir honey, it resists your spoon. That's viscosity (thickness). In normal fluids, if you push the fluid to the right, it flows to the right.

But in this chiral gas, something magical happens. Because the particles are biased, they create a new kind of resistance called "Odd Viscosity."

The Analogy:
Imagine you are trying to push a block of jelly across a table.

• Normal Viscosity: The jelly just resists your push. It's like friction.
• Odd Viscosity: If you push the jelly to the right, it doesn't just resist; it also tries to slide up or down the table, perpendicular to your push. It's like the fluid has a "twist" built into it.

This "odd" behavior is usually only seen in exotic quantum systems or under strong magnetic fields. This paper is special because they showed you can get this "twist" just by having particles that prefer one collision direction over the other, without needing any magnets or quantum magic.

4. The Math and the Simulation

The scientists used a famous mathematical tool called the Chapman-Enskog expansion.

• The Metaphor: Imagine trying to predict how a crowd moves. You start by assuming everyone is standing still (equilibrium). Then, you add small corrections for people walking slowly, then people running, then people bumping into each other. This "expansion" allows them to calculate exactly how thick the fluid is (viscosity) and how well it conducts heat.

They derived formulas for these "odd" properties. But math is great, and simulations are better. They built a computer model of these chiral disks and ran thousands of collisions.

• The Result: The computer simulation matched their math perfectly. The "odd viscosity" and "odd thermal conductivity" (heat flowing sideways) appeared exactly as predicted.

Why Does This Matter?

This is a "first principles" discovery. Before this, we knew odd viscosity existed in weird, complex systems (like electrons in a magnetic field or active bacteria). This paper says: "You don't need complex physics to get this. You just need particles that are slightly 'handed' in how they bump into each other."

It opens the door to designing new materials—like special fluids or active gels—that can flow in twisted, unexpected ways, simply by engineering the shape or collision rules of their tiny building blocks.

In a nutshell:
They took a gas of round balls, gave them a "lefty-loosey, righty-tighty" collision rule, proved they still calm down, and discovered that this simple trick creates a fluid that flows sideways when pushed, a phenomenon previously thought to require much more complex physics.

Technical Summary

1. Problem Statement

The paper addresses the kinetic theory of chiral fluids in two dimensions. Specifically, it investigates how microscopic chirality (handedness) in particle collisions affects macroscopic transport coefficients, particularly odd viscosity and odd thermal conductivity.

• Context: In standard kinetic theory, gases are modeled as symmetric hard spheres. Transport coefficients like shear viscosity and thermal conductivity arise from the breaking of time-reversal symmetry at the macroscopic scale, while microscopic interactions remain reversible.
• The Gap: While odd viscosity (a non-dissipative, parity-odd transport coefficient) is known to exist in systems with external magnetic fields or active particles driven by external torques, there was no first-principles derivation for a system where chirality arises intrinsically from the collision geometry itself, without external fields or active driving.
• The Challenge: Introducing chirality often breaks time-reversal and parity symmetries at the microscopic level. A major theoretical hurdle is ensuring that such a system still obeys the H-theorem (guaranteeing relaxation to a well-defined thermal equilibrium) to allow for a hydrodynamic expansion.

2. Methodology

The authors employ a multi-faceted approach combining a novel microscopic model, analytical kinetic theory, and numerical simulations.

A. The Microscopic Model: Chiral Hard Disks

The authors propose a "toy model" of hard disks where chirality is introduced via effective collision radii that depend on the collision's handedness (chirality $c$).

• Collision Rule: Two particles with relative velocity $\mathbf{g}$ collide. The chirality is defined as $c = \text{sign}[\hat{z} \cdot (\mathbf{n} \times \mathbf{g})]$, where $\mathbf{n}$ is the vector connecting the centers.
• Effective Radius: Depending on $c$, the effective radius becomes $r(1 + c\epsilon)$.
  • If $c=+1$, radius is $r(1+\epsilon)$.
  • If $c=-1$, radius is $r(1-\epsilon)$.
• Conservation Laws: Crucially, the collision rule (elastic scattering) remains unchanged; only the rate or probability of collision (determined by the impact parameter) is biased by chirality. This ensures strict conservation of energy, linear momentum, and angular momentum.

B. Analytical Approach: Chapman-Enskog Expansion

The authors perform a Chapman-Enskog expansion on the Boltzmann equation for this system.

1. H-Theorem Verification: They rigorously prove (Appendix B) that despite the breaking of microscopic time-reversal and parity symmetries, the H-theorem holds ($\partial H / \partial t \leq 0$). This guarantees the existence of a Maxwell-Boltzmann equilibrium distribution.
2. Linearized Boltzmann Equation: They expand the distribution function $f = f_0(1 + \Phi)$ around the local equilibrium $f_0$.
3. Sonine Polynomial Expansion: The correction term $\Phi$ is expanded in Sonine polynomials to solve the linearized collision integral. This allows for the analytical derivation of transport coefficients up to a specific order of truncation ($N$).

C. Numerical Validation: Non-Equilibrium Molecular Dynamics (NEMD)

To verify the analytical results, the authors perform NEMD simulations:

• System: A 2D periodic box of hard disks.
• Shear Flow: Shear is induced using the SLLOD algorithm.
• Extrapolation: To avoid non-linear effects (shear thinning), simulations are run at various shear rates $\gamma$, and the results are extrapolated to the limit $\gamma \to 0$.
• Cross-Section Validation: They also simulate "ratchet-shaped" particles with intrinsic rotation to confirm that the toy model's effective radii accurately reproduce the scattering cross-sections of more complex chiral geometries.

3. Key Contributions

1. First-Principles Derivation of Odd Viscosity: This is the first work to calculate odd viscosity and odd thermal conductivity for a passive gas of hard disks where chirality is an intrinsic property of the collision geometry, rather than an external field or active drive.
2. Proof of H-Theorem for Chiral Collisions: The authors demonstrate that microscopic breaking of time-reversal symmetry does not preclude the existence of a thermal equilibrium state, provided energy and momentum are conserved.
3. Analytical Expressions: They derive closed-form analytical expressions for the transport coefficients in the dilute limit.
4. Validation: The analytical predictions are quantitatively confirmed by NEMD simulations, bridging the gap between kinetic theory and many-body simulations for chiral fluids.

4. Results

The study yields explicit formulas for the shear viscosity ($\eta_e$), odd viscosity ($\eta_o$), thermal conductivity ($\kappa_e$), and odd thermal conductivity ($\kappa_o$) as functions of the chirality parameter $\epsilon$, temperature $T$, and particle mass $m$.

Zeroth-order transport coefficients (Eq. 14):

• Shear Viscosity: $\eta_e^{(0)} \propto \frac{1}{d(16+\epsilon^2)}\sqrt{m k_B T}$
• Odd Viscosity: $\eta_o^{(0)} \propto -\frac{2\epsilon}{d(16+\epsilon^2)}\sqrt{m k_B T}$
  • Key Finding: Odd viscosity is linear in the chirality parameter $\epsilon$ and vanishes if $\epsilon=0$. It is significantly smaller than shear viscosity (by an order of magnitude).
• Thermal Conductivity: $\kappa_e^{(0)} \propto \frac{1}{d(16+\epsilon^2)}\sqrt{\frac{k_B T}{m}}$
• Odd Thermal Conductivity (Righi-Leduc effect): $\kappa_o^{(0)} \propto -\frac{8\epsilon}{d(16+\epsilon^2)}\sqrt{\frac{k_B T}{m}}$

Numerical Findings:

• The NEMD results show strong agreement with the analytical predictions derived from the Chapman-Enskog expansion (truncated at $N=5$).
• Systematic Deviations:
  • At low temperatures, numerical viscosity is slightly underestimated due to non-linear shear rate dependence.
  • At high temperatures, odd viscosity shows a slight upward bias, attributed to finite resolution effects (the length scale of odd viscosity is $d\epsilon$, which becomes comparable to the simulation grid resolution at high $T$).
• Signal-to-Noise: The statistical error for odd viscosity is high because its magnitude is much smaller than that of shear viscosity.

5. Significance

• Theoretical Foundation: The paper establishes a rigorous kinetic theory framework for chiral passive fluids. It resolves the question of whether a system with intrinsic chiral collisions can reach equilibrium, confirming that the H-theorem remains valid.
• New Mechanism for Odd Transport: It identifies a new mechanism for generating odd viscosity that does not require magnetic fields (Hall effect) or active energy injection (active matter). This suggests that odd transport could be a generic feature of any system with chiral scattering cross-sections.
• Experimental Relevance: The model provides a theoretical basis for understanding odd transport in colloidal suspensions, granular materials, or biological systems where particles have chiral shapes or effective chiral interactions, even in the absence of external fields.
• Methodological Benchmark: The agreement between the exact Chapman-Enskog expansion and NEMD simulations serves as a benchmark for future studies of non-equilibrium transport in complex fluids.