High-harmonic generation in systems with chiral Bloch states: application to rhombohedral graphene

This paper demonstrates that rhombohedral graphene serves as a promising platform for exploring rich nonlinear optical phenomena, specifically showing how the linear scaling of chiral Bloch state winding with layer number $n$ and interaction-induced valley splitting govern the system's high-harmonic generation and circular dichroism.

The Gist

Imagine you have a piece of graphene, which is essentially a single layer of carbon atoms arranged in a honeycomb pattern (like a chicken wire fence). Now, imagine stacking several of these layers on top of each other, but not just any way. You stack them in a specific "rhombic" pattern, like a spiral staircase where each step is slightly rotated. This creates a material called Rhombohedral Graphene.

The scientists in this paper asked a simple question: What happens when we blast this material with a super-powerful laser?

Specifically, they wanted to see if the material would spit out new colors of light (harmonics) and, more importantly, if the shape of the electrons inside the material would leave a fingerprint on that light.

Here is the breakdown of their discovery using everyday analogies:

1. The "Spiral Staircase" of Electrons

In normal materials, electrons move in straight lines or simple loops. But in this stacked graphene, the electrons live in a special kind of "quantum landscape." Because of the way the layers are stacked, the electrons have a property called chirality (handedness).

Think of the electrons as dancers on a stage. In a normal crowd, they might just shuffle around. But in this rhombohedral graphene, the dancers are forced to spin in a specific direction as they move.

• The "Winding": The paper found that the more layers you stack (let's say $n$ layers), the more times the dancers have to spin to get from one side of the stage to the other. If you have 4 layers, they spin 4 times. If you have 6, they spin 6 times. This "winding" is the key to the whole story.

2. The Laser as a DJ

The researchers hit this material with an intense laser pulse. Imagine the laser is a DJ playing a beat (the fundamental frequency).

• The Reaction: Usually, when you play a beat, the crowd just bounces to that beat. But because the electrons in this graphene are spinning so wildly (due to the "winding" mentioned above), they start shouting out new beats that are much faster than the original song.
• High-Harmonic Generation (HHG): This is the scientific term for the material spitting out these new, higher-pitched "notes" (light frequencies).

3. The Magic Number: $2n \pm 1$

Here is the coolest part of the discovery. The scientists found a direct link between the number of layers and the "notes" the material sings.

• If you have $n$ layers, the material doesn't just sing random notes. It sings its loudest, strongest notes at frequencies that are $2n + 1$ and $2n - 1$ times the original laser frequency.
• Analogy: It's like if you have a 3-string guitar, and you strum it, it magically produces a chord that is exactly the 7th and 5th harmonics of the note you played. The number of layers dictates the pitch of the loudest sound.

4. The "Left-Handed" vs. "Right-Handed" Light

The laser used in the experiment wasn't just a straight beam; it was circularly polarized, meaning the light itself was spinning like a corkscrew (either left-handed or right-handed).

• Because the electrons in the graphene are also "chiral" (they have a handedness), they react differently to left-spinning light vs. right-spinning light.
• The paper measured something called Circular Dichroism. Think of it as a "chirality detector."
• The Surprise: When they looked at the light coming out, they found that the "handedness" of the output light changed depending on how many layers were in the stack. For some layer counts, the material preferred left-handed light; for others, it flipped and preferred right-handed light. It's as if the material's "personality" changed based on its thickness.

5. The "Ring" of Activity

The scientists also discovered where this magic happens. The electrons don't react everywhere in the material.

• Imagine a bullseye target. The center is boring, and the very edge is boring. But there is a specific ring in the middle where all the action happens.
• When the laser hits, the electrons on this specific ring start dancing the hardest. The paper showed that the pattern of this "dance floor" is directly imprinted on the light that comes out.

Why Does This Matter?

This isn't just about making cool light shows.

1. New Tech: This proves that we can use light to "read" the internal quantum structure of materials without breaking them. It's like using X-ray vision to see the spin of electrons.
2. Information Processing: Since the material reacts differently to left and right spinning light, it could be used to build ultra-fast optical computers that process information using the "handedness" of light instead of electricity.
3. Predictability: The fact that the output frequency scales perfectly with the number of layers ($2n \pm 1$) means engineers can design materials to produce specific colors of light just by stacking a precise number of graphene sheets.

In a nutshell: The researchers stacked graphene layers like a spiral staircase, hit them with a spinning laser, and found that the material sings specific, predictable notes that reveal the "handedness" and "spin" of the electrons inside. It's a new way to turn quantum physics into a tool for future technology.

Technical Summary

1. Problem Statement

The paper investigates the nonlinear optical response, specifically High-Harmonic Generation (HHG), in rhombohedral multilayer graphene (RMG). While HHG has been extensively studied in single-layer graphene and twisted bilayer graphene, the unique properties of RMG remain under-explored in this context.

Key motivations include:

• Chiral Bloch States: RMG lacks $C_{2z}$ rotational symmetry and inversion symmetry, leading to chiral Bloch states with non-zero Berry curvature ($\Omega \neq 0$).
• Layer-Dependent Topology: The "winding" of the Bloch states scales linearly with the number of layers ($n$), and the Berry curvature is concentrated on a ring in momentum space.
• Valley Imbalance: RMG exhibits strongly correlated physics, including spontaneous valley polarization (imbalance between the $K$ and $K'$ valleys), which breaks time-reversal symmetry.
• Gap: The authors aim to understand how the interplay between the layer-dependent quantum geometry, valley chirality, and interaction-induced valley splitting affects HHG spectra and circular dichroism (CD).

2. Methodology

The authors employ a combination of analytical modeling and numerical simulations:

• Effective Two-Band Model: They utilize a minimal continuum model for the low-energy physics of $n$-layer RMG. The Bloch Hamiltonian for a valley $\xi = \pm$ is:
  $$h^{(n)}_\xi(\mathbf{k}) = \begin{pmatrix} -\mu + w & -\gamma_1(-k^*_\xi/k_c)^n \\ -\gamma_1(-k_\xi/k_c)^n & -\mu - w \end{pmatrix}$$
  where $k_\xi = k_x + i\xi k_y$. This model captures the $n$-th order momentum dependence characteristic of RMG.
• Light-Matter Interaction: The system is coupled to a strong pulsed laser field using the length gauge. The dynamics are governed by the Semiconductor Bloch Equations (SBEs), which track the population of bands ($n_{p,\xi}$) and interband coherence ($\pi_\xi$).
• Current Calculation: The total current $\mathbf{J}(t)$ is decomposed into intraband (anomalous velocity + group velocity) and interband (polarization derivative) contributions.
• Spectral Analysis: The HHG spectrum is obtained via the Fourier transform of the current (Larmor formula).
• Circular Dichroism (CD): To probe chirality, the authors calculate the dimensionless CD parameter for the $l$-th harmonic:
  $$CD_l = \frac{I^+_{RCP}(l) - I^-_{LCP}(l)}{I^+_{RCP}(l) + I^-_{LCP}(l)}$$
  where $I^\pm$ represents the intensity for Right/Left Circularly Polarized light.
• Valley Imbalance: To simulate interaction effects, a mean-field valley imbalance order parameter $\Phi_V$ is introduced, creating an energy splitting between the two valleys ($\xi = +$ and $\xi = -$).

3. Key Contributions and Results

A. Single-Valley Physics and Quantum Geometry

• Momentum-Space Distribution: At the onset of the laser pulse, the population in the conduction band is not determined by band energy but by the dipole matrix elements. These are maximized on a ring of finite radius in momentum space, mirroring the distribution of the Berry curvature.
• Harmonic Scaling: The dominant harmonic order ($l$) above the bandgap scales linearly with the number of layers:
  $$l \approx 2n \pm 1$$
  This scaling is directly linked to the winding number of the Bloch states. The interband current dominates these high-order harmonics.
• Intraband vs. Interband:
  • Interband: Dominates for $l \geq 2n-1$.
  • Intraband: Shows distinct behavior based on the parity of $n$. For odd $n$, it dominates up to the $n$-th harmonic; for even $n$, it is significant only for the first harmonic.
• Circular Dichroism (CD):
  • The CD for the dominant harmonic ($l=2n+1$) is nearly saturated ($|CD| \approx 1$), indicating strong chirality.
  • The ratio $c_n = CD_{2n+1} / CD_1$ provides a sensitive measure of the winding degree. The authors find that $\ln|c_n|$ scales linearly with $n$, offering a quantitative probe of the topological winding.
• Doping Effects: While the net Berry curvature ($\Phi_B$) of occupied states varies significantly with chemical potential ($\mu$), the dominant CD ($CD_{2n+1}$) remains saturated over a wide range of $\mu$. However, the ratio $c_n$ tracks the Berry curvature more closely, peaking near charge neutrality.

B. Two-Valley Interplay and Valley Imbalance

• Valley Splitting: When both valleys are active with an interaction-induced gap splitting ($\Delta_{gap} \propto \Phi_V$), the system exhibits a complex interplay.
• Chirality Competition: The two valleys have opposite chiralities.
  • For small $n$ ($n \leq 3$), the valley with the smaller gap dominates the CD signal.
  • For larger $n$ ($n \geq 4$), the valley with the larger gap determines the sign of the CD.
• Sign Reversal: A critical finding is that the sign of the dominant circular dichroism ($CD_{2n+1}$) changes as a function of $n$. At $n=3$, there is a near-cancellation of the opposite chiralities. At $n=6$, the CD approaches the maximum value of $+1$, determined by the larger-gap valley.

4. Significance and Implications

• Probing Quantum Geometry: The paper establishes HHG as a powerful tool to probe the quantum geometry (specifically the winding of Bloch states) in multilayer systems. The linear scaling of the dominant harmonic order with layer count ($n$) is a unique fingerprint of RMG.
• Valleytronics: The study demonstrates that HHG and Circular Dichroism can detect and quantify valley imbalance and interaction-induced symmetry breaking, which are crucial for valleytronic applications.
• Topological Quantification: The proposed metric $c_n = CD_{2n+1}/CD_1$ offers a robust, experimentally accessible way to quantify the degree of chirality and topological winding in chiral multiband systems.
• Material Platform: Rhombohedral graphene is identified as a promising platform for exploring rich nonlinear optical phenomena, bridging the gap between correlated electron physics and ultrafast optics.

In summary, the work provides a comprehensive theoretical framework linking the layer-dependent topological properties of rhombohedral graphene to its nonlinear optical response, predicting distinct signatures (harmonic scaling and CD sign reversal) that can be used to experimentally verify the chiral nature of these materials.