Anomalous Hall Crystals in Rhombohedral Multilayer Graphene II: General Mechanism and a Minimal Model

Tomohiro Soejima,Junkai Dong,Taige Wang,Tianle Wang,Michael P. Zaletel,Ashvin Vishwanath,Daniel E. Parker

2024-08-21

Abstract:We propose a minimal "three-patch model" for the anomalous Hall crystal (AHC), a topological electronic state that spontaneously breaks both time-reversal symmetry and continuous translation symmetry. The proposal for this state is inspired by the recently observed integer and fractional quantum Hall states in rhombohedral multilayer graphene at zero magnetic field. There, interaction effects appear to amplify the effects of a weak moiré potential, leading to the formation of stable, isolated Chern bands. It has been further shown that Chern bands are stabilized in mean field calculations even without a moiré potential, enabling a realization of the AHC state. Our model is built upon the dissection of the Brillouin zone into patches centered around high symmetry points. Within this model, the wavefunctions at high symmetry points fully determine the topology and energetics of the state. We extract two quantum geometrical phases of the non-interacting wavefunctions that control the stability of the topologically nontrivial AHC state. The model predicts that the AHC state wins over the topological trivial Wigner crystal in a wide range of parameters, and agrees very well with the results of full self-consistent Hartree-Fock calculations of the rhombohedral multilayer graphene Hamiltonian.

Strongly Correlated Electrons,Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the Quantum Anomalous Hall Effect (QAHE) in Rhombohedral Multilayer Graphene (RMG) under zero - magnetic - field conditions. Specifically, the researchers observed that in the absence of an external periodic potential, electron - electron interactions can spontaneously break continuous translational symmetry and form stable, isolated Chern Insulators. However, how this phenomenon can stably exist under the influence of a weak moiré potential remains an unsolved mystery.

The paper proposes a minimal "three - patch model" to explain this phenomenon by analyzing the high - symmetry points in the Brillouin zone. The core of this model is that the behavior of the wave function at these high - symmetry points completely determines the topological properties and energy states of the system. It was found that even in the absence of a moiré potential, the Anomalous Hall Crystal (AHC) state can be achieved by spontaneously breaking translational symmetry. In addition, this model predicts that, over a wide range of parameters, the AHC state has an advantage over the topologically trivial Wigner crystal and is in very good agreement with the results of Self - Consistent Hartree - Fock (SCHF) calculations.

### Key problem summary:
1. **Spontaneous translational symmetry breaking**: How do electron - electron interactions lead to the spontaneous breaking of continuous translational symmetry in the absence of an external periodic potential?
2. **Stability of Chern insulators**: How do Chern insulators remain stable under the influence of a weak moiré potential?
3. **Generation mechanism of topologically non - trivial states**: How is the topologically non - trivial AHC state generated through spontaneous translational symmetry breaking?
4. **Verification of the theoretical model**: Verify the validity of the three - patch model through Self - Consistent Hartree - Fock calculations.

### Summary of mathematical formulas:
- **C3 symmetry action**:
\[
\hat{C}_3|q, a\rangle=\exp\left(\frac{2\pi i\ell_a}{3}\right)|C_3q, a\rangle, \quad \ell_a\in\mathbb{Z}/3\mathbb{Z}
\]
- **Single - particle wave function**:
\[
\phi_a^q(r)=\langle r, a|\phi_q\rangle = \chi_a(q)e^{i q\cdot r}
\]
- **Hamiltonian**:
\[
\hat{H}_0=\hat{h}+\hat{H}_{\text{int}}=\hat{h}+\frac{1}{2}\int d^2r d^2r' V(r - r'):\hat{\rho}(r)\hat{\rho}(r'):
\]
- **Hartree and Fock energies**:
\[
E_H(\xi, \xi')=\int d^2r d^2r'\psi_a^\xi(r)\psi_a^\xi(r)V(r - r')\psi_b^{\xi'}(r)\psi_b^{\xi'}(r)
\]
\[
E_F(\xi, \xi')=-\int d^2r d^2r'\psi_a^{\xi'}(r)\psi_a^\xi(r)V(r - r')\psi_b^\xi(r')\psi_b^{\xi'}(r')
\]
- **Energy difference**:
\[
\Delta E:=\epsilon_H(1)\sqrt{3}+\epsilon_F(1)-\min_{C = 0, - 1}\left[\cdots\right]