Quantifying chirality of phonons

✨

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 a crystal not as a rigid, static block of stone, but as a bustling city where every atom is a dancer. Even when the crystal looks perfectly still, these atomic dancers are constantly vibrating, wiggling, and spinning in place. In physics, these vibrations are called phonons.

For a long time, scientists knew that some of these dancers could spin in a specific direction, like a figure skater doing a pirouette. This spinning gives them "handedness" (chirality)—they can be left-handed or right-handed, just like your hands are mirror images of each other.

However, measuring exactly how chiral a crystal is, or how much "spin" is happening overall, has been like trying to count the number of people spinning in a crowded dance hall without a clear way to distinguish the left-spinners from the right-spinners. It was a bit of a mystery.

This paper introduces a new, clever way to quantify (measure) this spin. Here is the breakdown of their discovery using simple analogies:

1. The Problem: The "Spinning Crowd"

Imagine a massive dance floor (the crystal).

In some crystals (like Centrosymmetric ones, e.g., Silicon), for every dancer spinning left, there is an identical partner spinning right. They cancel each other out perfectly. The net spin is zero.
In other crystals (like Chiral ones, e.g., Quartz), the whole dance floor is biased. Maybe more dancers are spinning left, or the way they spin is inherently twisted.
In a third group (like Gallium Arsenide), the dancers can spin individually, but the crowd is so mixed that the total spin still averages to zero.

The challenge was: How do we measure the "twistiness" of the whole crowd, not just individual dancers?

2. The Solution: Two New "Spin Meters"

The authors created two new tools to measure this:

Tool A: The "Spotlight" (Momentum-Resolved Chirality)

Think of this as a spotlight that scans the dance floor one specific dancer at a time.

It asks: "Are you spinning with the flow of the music (right-handed) or against it (left-handed)?"
This allows scientists to see a map of the entire crystal, showing exactly where the left-spinners and right-spinners are located.
The Result: In chiral crystals (like Quartz), the spotlight shows a clear, consistent bias across the whole floor. In non-chiral crystals, the spotlight sees spins, but they are scattered and random.

Tool B: The "Crowd Count" (Bulk Dynamical Chirality)

This is the big innovation. Instead of looking at one dancer, this tool counts the entire crowd at once, taking into account how hot the room is (temperature).

Imagine a thermometer that doesn't measure heat, but measures the net imbalance of the spinning.
If the room is full of equal numbers of left and right spinners, the meter reads Zero.
If the room is full of mostly left spinners (because the crystal structure forces them to be), the meter reads a Positive Number.
If the room is full of mostly right spinners, the meter reads a Negative Number.

3. The "Handedness" Test

The authors tested this on different materials:

The "Perfectly Balanced" (Silicon): The meter reads 0. Even though atoms wiggle, the left and right spins cancel out perfectly.
The "Twisted" (Quartz, Selenium, Tellurium): The meter reads a strong number. Crucially, if you take a "Left-Handed" crystal of Quartz, the meter reads a negative number. If you take its mirror-image "Right-Handed" twin, the meter reads a positive number.
- Analogy: It's like a scale that can tell the difference between a left glove and a right glove just by weighing them, even though they look identical from a distance.
The "Confused" (Gallium Arsenide): The individual dancers spin, but the total count is 0. The crystal looks twisted locally, but globally, it's balanced.

Why Does This Matter?

Previously, scientists had to look at the crystal's structure (its shape) to guess if it was chiral. This new method looks at the dynamics (the movement).

It's a fingerprint: It can distinguish between two crystals that look identical in shape but have opposite "handedness" in their vibrations.
It's a new language: It gives us a way to talk about "chirality" not just as a static shape, but as a living, breathing property of how atoms move and transfer energy.

The Bottom Line

This paper gives us a mathematical ruler to measure the "spin" of atoms in a crystal. It proves that we can now detect the "handedness" of a material simply by listening to how its atoms vibrate, distinguishing between materials that are truly twisted and those that are just pretending to be. It turns the invisible dance of atoms into a measurable quantity.

1. Problem Statement

While "chiral phonons" (lattice vibrations carrying angular momentum and exhibiting handedness) have attracted significant attention due to their interactions with circularly polarized light and role in angular momentum transfer, a rigorous quantitative framework for characterizing their chirality has been lacking.

The Gap: Existing approaches typically rely on the helicity of single modes or specific experimental probes. There is no intrinsic, dynamical measure that defines the degree of chirality independent of specific experimental setups.
The Challenge: It remains an open question how to define and evaluate the "degree of chirality" for phonon modes in a way that is intrinsic to lattice dynamics and can distinguish between enantiomers (left- and right-handed crystal structures) and different types of achiral materials (centrosymmetric vs. noncentrosymmetric).

2. Methodology

The authors propose a theoretical framework to quantify phonon chirality using first-principles calculations based on Density Functional Perturbation Theory (DFPT) implemented in the abinit code. They introduce two distinct quantitative measures:

A. Momentum-Resolved Dynamical Chirality

Definition: Based on Barron's definition of "true chirality," the authors define the chirality of a phonon mode with branch $j$ and wave vector $\mathbf{k}$ as the inner product of the phonon angular momentum $\mathbf{L}_j(\mathbf{k})$ and the wave vector $\mathbf{k}$ .
Formula: $\chi(\mathbf{k}) = \mathbf{L}_j(\mathbf{k}) \cdot \hat{\mathbf{k}}$ , where $\hat{\mathbf{k}}$ is the normalized wave vector.
Properties:
- This quantity is odd under mirror reflection ( $M$ ) and even under time reversal ( $T$ ).
- Positive and negative values correspond to right- and left-handed chiral phonon modes, respectively.
- Normalization constrains values to the interval $[-\hbar, \hbar]$ .
Purpose: Provides a mode- and wave-vector-resolved visualization of chirality across the Brillouin zone.

B. Bulk Dynamical Chirality

Motivation: The momentum-resolved measure does not yield a single material-level quantity suitable for comparing different crystal systems.
Definition: A bulk quantity obtained by summing the momentum-resolved chirality over the entire Brillouin zone, weighted by the Bose-Einstein distribution $f(\omega)$ (representing thermal population at $T=300$ K) and a structure factor $\mathbf{F}_1(\mathbf{k})$ derived from the crystal's point-group symmetry.
Key Indicators: The authors project the phonon chirality onto the lowest-order basis functions compatible with the lattice symmetry, specifically using:
- Electric Toroidal Monopole ( $G_0$ ): Captures the isotropic component of bulk true chirality.
- Electric Toroidal Quadrupole ( $G_u$ ): Describes uniaxial anisotropy (or $G_4$ for cubic systems).
Formulation:
$G_0(T) = \frac{1}{N_0\hbar} \sum_{j,\mathbf{k}} f(\omega_j(\mathbf{k})) [\mathbf{L}_j(\mathbf{k}) \cdot \mathbf{F}_1(\mathbf{k})]$
(A similar expression defines $G_u$ ).
Structure Factors: Explicitly derived for point groups $D_3$ (chiral), $O_h$ (centrosymmetric achiral), and $T_d$ (noncentrosymmetric achiral) using group theory projection operators.

3. Key Contributions

Theoretical Framework: Established the first quantitative indicators ( $G_0$ and $G_u$ ) for the collective dynamical chirality of thermally populated phonons.
Symmetry-Based Classification: Demonstrated how these measures distinguish between three classes of materials:
- Chiral Crystals: Possess structural handedness.
- Centrosymmetric Achiral Crystals: Possess inversion symmetry.
- Noncentrosymmetric Achiral Crystals: Lack inversion symmetry but possess other symmetries preventing net chirality.
Validation: Validated the framework using first-principles calculations on eight materials and cross-verified results using an alternative finite-displacement method (VASP/Phonopy).

4. Results

The authors computed these quantities for representative materials: $\alpha$ -quartz ( $D_3$ ), Se, Te, $\alpha$ -HgS (all chiral), and Si ( $O_h$ ), GaAs, GaP, ZnTe (achiral).

Momentum-Resolved Chirality:
- Chiral Materials ( $\alpha$ -quartz, Se, Te, HgS): Show non-zero chirality values across the Brillouin zone. Crucially, the sign of the chirality flips completely between left- and right-handed enantiomers (e.g., L-SiO $_2$ vs. R-SiO $_2$ ).
- Centrosymmetric Achiral (Si): The chirality vanishes ( $\mathbf{L} \cdot \mathbf{k} = 0$ ) at all wave vectors due to the combined presence of spatial inversion and time-reversal symmetries.
- Noncentrosymmetric Achiral (GaAs, GaP, ZnTe): Finite, non-zero chirality values appear at specific wave vectors (e.g., W point in GaAs) where local symmetry is low. However, these values are not uniform across the zone.
Bulk Dynamical Chirality ( $G_0, G_u$ ):
- Chiral Materials: Both $G_0$ and $G_u$ take finite, non-zero values. The signs reverse between opposite enantiomers (e.g., $G_0 \approx -0.44$ for L-SiO $_2$ and $+0.42$ for R-SiO $_2$ ). This indicates a net population imbalance of left- vs. right-handed phonons.
- Achiral Materials: Despite having finite local chirality in noncentrosymmetric cases (like GaAs), the bulk dynamical chirality vanishes ( $G_0 \approx 0, G_u \approx 0$ ) within numerical accuracy. This is because the contributions from left- and right-handed modes cancel out perfectly over the entire Brillouin zone due to the crystal's symmetry.
Temperature Dependence: The bulk quantities $G_0$ and $G_u$ scale approximately linearly with temperature ( $T$ ) in the high-temperature limit, consistent with the Bose-Einstein distribution behavior.

5. Significance and Discussion

Distinction of Enantiomers: The bulk dynamical chirality serves as a robust theoretical tool to distinguish enantiomers of chiral crystals without direct structural measurement, reflecting the underlying structural handedness through dynamical properties.
Order Parameter Status: The authors argue that while $G_0$ and $G_u$ fulfill some criteria of an order parameter (vanishing in high-symmetry phases, changing sign between enantiomers), they are not true Landau order parameters because they are not directly measurable observables and do not couple to a well-defined conjugate field. Instead, they are best viewed as dynamical quantities quantifying the net population imbalance of chiral phonons.
Implications: This framework bridges the gap between structural chirality and dynamical phenomena. It suggests that bulk dynamical chirality could be linked to experimentally accessible quantities, such as chiral optical responses and angular momentum transfer processes, providing a new route for experimental validation.
Future Work: The paper identifies the need to find external fields or non-equilibrium driving schemes that can selectively couple to these chiral phonons to control chirality dynamically.

In summary, the paper successfully moves the field from qualitative descriptions of chiral phonons to a rigorous, quantitative, and symmetry-resolved characterization, providing a new metric ( $G_0, G_u$ ) to understand and distinguish chiral materials based on their lattice dynamics.