\textit{Ab initio} \textit{GW}-BSE theory of optical activity in $\alpha$-quartz

This paper presents an *ab initio* many-body theory of optical activity in solids using the GW-BSE framework, demonstrating that a combination of exciton envelope modulation and sum-over-exciton-states formulations is necessary to accurately capture the full frequency-dependent spectral dispersion of $\alpha$-quartz, a feat unachievable by independent-particle or simple local-field approximations.

The Gist

The Big Picture: Twisting Light with a Crystal

Imagine you have a beam of white light. Now, imagine that light is made of two types of "spirals": one twisting clockwise (right-handed) and one twisting counter-clockwise (left-handed).

In most materials, these two spirals travel at the exact same speed. But in optically active materials (like the crystal α-quartz found in quartz watches), something magical happens: the crystal treats the two spirals differently. One slows down more than the other.

When these two spirals recombine after passing through the crystal, they form a new beam of light, but its "twist" (polarization) has rotated. This is called optical rotation. It's the same effect that makes sugar water or certain crystals look different when viewed through a polarized lens.

For over 200 years, scientists have known this happens in quartz. But they have struggled to predict exactly how much the light will twist at different colors (frequencies) using only the laws of physics, without just guessing based on experiments.

The Problem: The "Solo" vs. The "Duet"

To predict this, scientists use computer simulations.

• The Old Way (Independent Particles): Imagine trying to predict the behavior of a crowded dance floor by watching each dancer move alone, ignoring everyone else. This is the "Independent Particle Approximation." In the world of crystals, this means looking at electrons as if they are solo dancers.

  • The Result: This method failed miserably for quartz. It predicted the wrong amount of rotation, sometimes even the wrong direction! It was like predicting a dance routine by ignoring the fact that dancers hold hands.

• The Missing Link (Excitons): In reality, when light hits a crystal, it doesn't just bump into a single electron. It creates a "duet" between an electron (which gets excited) and a "hole" (the empty space it left behind). They are attracted to each other by electricity and dance together as a pair. This pair is called an Exciton.

  • To get the right answer, you have to simulate the duet, not just the solo.

The New Theory: The "GW-BSE" Framework

The authors of this paper, Xiaoming Wang and Yanfa Yan, developed a super-advanced simulation method called GW-BSE.

• GW: This part calculates the energy levels of the solo dancers very accurately.
• BSE (Bethe-Salpeter Equation): This is the "duet" engine. It forces the computer to calculate how the electron and hole dance together.

They applied this to α-quartz, the classic crystal used to test these theories.

The Two Formulas: Two Ways to Describe the Dance

The paper presents two different mathematical ways to describe how this "exciton duet" interacts with the twisting light. Think of them as two different ways to describe a complex dance move:

1. The "Envelope Modulation" Method:

   • Analogy: Imagine the dancers are holding a giant, stretchy blanket (the "envelope"). This method looks at how the shape of the blanket changes slightly as the light wave passes through.
   • Performance: It works great for low-energy light (like deep red or infrared). It's simple and captures the basic shape of the dance.

2. The "Sum-Over-Exciton-States" (SOXS) Method:

   • Analogy: This method is like listing every single possible move the dancers could make, adding them all up to see the final result. It's much more detailed and accounts for the fact that the dancers might change their steps depending on the music's pitch.
   • Performance: This is the heavy lifter. It captures the full frequency range. It correctly predicts how the rotation changes as the light goes from red to blue to ultraviolet.

The Verdict: The "Envelope" method is a good shortcut for simple cases, but the "Sum-Over-States" method is the only one that gets the entire song right.

Why This Matters

Before this paper, scientists could only guess the optical rotation of quartz accurately at a single, static point (like a snapshot). They couldn't predict how it would behave across the whole spectrum of light.

By using this new GW-BSE theory with the SOXS method, the authors finally achieved a "perfect match" between their computer simulation and real-world experiments.

• They fixed the sign errors (making sure the light twists the right way).
• They got the numbers right across all colors of light.

The Takeaway

This paper is like solving a 200-year-old puzzle. It proves that to understand how light twists in crystals, you can't just look at individual electrons. You have to look at the electron-hole duets and how they dance together.

This is a huge step forward for materials science. If we can predict exactly how a material twists light, we can design new, better materials for:

• 3D glasses and VR headsets (which rely on polarized light).
• Chemical sensors (to detect if a drug molecule is the "right" or "left" handed version).
• Quantum computers (which use light to carry information).

In short: They finally figured out the choreography of the crystal, and now we can design new crystals to dance exactly how we want them to.

Technical Summary

1. Problem Statement

Optical activity (the rotation of the polarization plane of light) is a fundamental property of chiral materials, yet a predictive ab initio theory for the full frequency dependence of optical activity in periodic solids has remained elusive.

• Theoretical Gap: While molecular optical activity is well-described by multipole expansions, applying this to solids is hindered by the ill-defined position operator ($\mathbf{r}$) in periodic boundary conditions.
• Limitations of Previous Methods:
  • Independent-Particle Approximation (IPA): Previous studies using IPA significantly underestimated optical rotation in $\alpha$-quartz.
  • Local-Field Corrections (LFCs): While LFCs improved static-limit predictions, they relied on empirically tuned band gaps and failed to capture the full frequency dependence (optical rotatory dispersion).
  • Missing Many-Body Effects: Existing approaches largely ignored electron-hole interactions (excitonic effects), which are crucial for wide-gap insulators like $\alpha$-quartz.
• Goal: To develop a rigorous, predictive ab initio framework that incorporates many-body effects (specifically within the GW-BSE formalism) to accurately describe optical activity across the entire frequency spectrum.

2. Methodology

The authors formulated optical activity within the GW-BSE (Bethe-Salpeter Equation) framework, expanding the dielectric tensor $\epsilon_{ij}(\omega, \mathbf{q})$ to the first order in the wavevector $\mathbf{q}$:
$$\epsilon_{ij}(\omega, \mathbf{q}) = \epsilon_{ij}(\omega, 0) + i q_l \gamma_{ijl}(\omega) + O(q^2)$$
where $\gamma_{ijl}$ is the optical activity tensor.

They developed two complementary formulations to handle the $\mathbf{q}$-dependence of the exciton oscillator strength $\rho_\lambda(\mathbf{q})$:

A. Exciton Envelope Modulation

• Concept: Expands the exciton envelope function $A^\lambda_{cv\mathbf{k}}(\mathbf{q})$ and the oscillator strength directly.
• Approximation: Assumes the linear-in-$\mathbf{q}$ contribution from the BSE Hamiltonian vanishes (valid for 3D bulk materials near the zone center).
• Implementation: Uses a multipole expansion (magnetic dipole and electric quadrupole) derived from the envelope-modulated Bloch states.
• Velocity Operator: Tested with different velocity operators: $v_{DFT}$ (from DFT Hamiltonian), $v_{GW}$ (from GW Hamiltonian), and $v_{opt}$ (scissors-shifted to match the BSE optical gap).

B. Sum-Over-Exciton-States (SOXS) Expansion

• Concept: A more rigorous approach that inserts a complete set of exciton states between the position and velocity operators in the oscillator strength definition.
• Derivation:
  $$\rho_\lambda(\mathbf{q}) = \frac{1}{2} \sum_\mu \left[ \langle \Psi_\lambda | e^{i\mathbf{q}\cdot\mathbf{r}} | \Psi_\mu \rangle \langle \Psi_\mu | \mathbf{v} | 0 \rangle + \langle \Psi_\lambda | \mathbf{v} | \Psi_\mu \rangle \langle \Psi_\mu | e^{i\mathbf{q}\cdot\mathbf{r}} | 0 \rangle \right]$$
• Key Innovation: This formulation explicitly uses the velocity operator consistent with the excitonic Hamiltonian ($\mathbf{v} = (i/\hbar)[H_{BSE}, \mathbf{r}]$), avoiding gauge ambiguities.
• Result: Defines exciton magnetic dipole ($M$) and electric quadrupole ($Q$) contributions via inter-exciton transition dipole moments ($R_{\lambda\mu}$).

Computational Details:

• Material: $\alpha$-quartz (Space group $P3_221$).
• Codes: VASP (DFT, GW, BSE).
• Parameters: PBE functional for DFT; $GW_0$ scheme for quasiparticle corrections; BSE with 18 valence and 24 conduction bands; $10\times10\times10$ k-mesh.

3. Key Contributions

1. First Full-Frequency GW-BSE Theory: The paper presents the first ab initio many-body theory capable of predicting the full frequency dependence of optical activity in solids, moving beyond the static limit.
2. Dual Formulations: The derivation of two distinct but complementary approaches (Envelope Modulation and SOXS) to handle spatial dispersion in the presence of strong excitonic effects.
3. Resolution of Velocity Operator Inconsistency: The authors identified that using a single-particle velocity operator ($v_{DFT}$ or $v_{GW}$) within the excitonic framework leads to discrepancies. The SOXS formulation resolves this by ensuring the velocity operator is consistent with the excitonic Hamiltonian.
4. Quantitative Agreement: The work successfully reproduces experimental optical rotatory dispersion for $\alpha$-quartz, a system where previous theories failed.

4. Results

• Static Limit ($\omega \to 0$):
  • IPA: Predicted a negative rotation ($\bar{\rho} \approx -0.7$), contradicting the experimental sign and magnitude.
  • LFC (with DFT): Improved magnitude ($\bar{\rho} \approx 6.5$) but still relied on empirical tuning.
  • GW-BSE (Envelope Modulation): With $v_{GW}$, the result was $5.8$; with $v_{opt}$ (scissors shift), it was $5.1$. These are close to the experimental value of $4.6$, but the method is sensitive to the band gap alignment.
• Frequency Dependence (Optical Rotatory Dispersion):
  • Envelope Modulation: Failed to reproduce the experimental frequency dependence, particularly at higher energies, because the simple scissors shift only aligns the lowest exciton, neglecting dynamical electron-hole effects for higher states.
  • SOXS Formulation: Successfully reproduced the experimental frequency dependence across the entire spectrum. The curve matched experimental data points from Ref. [36] with high fidelity.
• Role of Many-Body Effects: The comparison confirmed that excitonic effects (electron-hole attraction) are decisive. Without them (IPA), the theory fails. With them (GW-BSE), the correct spectral dispersion is captured.

5. Significance

• Fundamental Physics: This work resolves a long-standing challenge in condensed matter physics by providing a rigorous theoretical foundation for optical activity in periodic solids, bridging the gap between molecular multipole theories and solid-state many-body physics.
• Predictive Power: The ability to predict optical activity ab initio without empirical tuning allows for the screening and design of new chiral materials.
• Chiroptoelectronics: The methodology opens a pathway for the rational design of chiroptoelectronic materials (e.g., for circularly polarized light detection or generation) by accurately modeling light-matter interactions in chiral crystals.
• Generalizability: The SOXS formulation is applicable to other chiral systems, including topological materials and halide perovskites, where excitonic effects and spatial dispersion are critical.

In conclusion, Wang and Yan demonstrate that a consistent many-body treatment of excitons, specifically using the Sum-Over-Exciton-States expansion within the GW-BSE framework, is essential for accurately describing the optical activity of solids, successfully resolving discrepancies that plagued previous independent-particle and local-field correction approaches.