Transdimensional Chern insulators in rhombohedral multilayer graphene superlattices

The gist: Interaction-driven reconstruction in twisted rhombohedral multilayer graphene/hBN moir´e superlattices transforms a symmetry-broken metallic Fermi-lune state into a quantized Chern insulator with transdimensional orbital magnetization.

How it works

1. The starting point is the Fermi lune state in pristine RMG, which spontaneously breaks time-reversal symmetry T, threefold rotation C3, and crystalline mirror symmetries due to electron doping [44]. This state carries coherent interlayer current loops and generates both out-of-plane and in-plane orbital magnetizations [44].

2. Coupling this metallic Fermi lune state to the moir´e potential in twisted RMG/hBN moir´e superlattices folds the noninteracting bands into the superlattice Brillouin zone, generating a narrow conduction-band manifold [10].

3. Electron-electron interactions then cause electrons to spontaneously occupy one spin-valley flavor, opening a global gap and isolating an occupied Chern band with C = 1 [45]. This resulting transdimensional Chern insulator (TDCI) state breaks rotation, mirror and time-reversal symmetry while inheriting the mixed orbital current-loop characters from its metallic parent state [45].

4. A minimal two-band model is introduced to describe the TDCI, which unveils the intimate connection between non-reciprocal velocity and in-plane orbital magnetization [45]. This connection is formalized by linking in-plane orbital magnetization to non-reciprocal group velocities via Equation (3) [45].

Key Findings and Properties

Band unfolding reveals that such a Chern insulating state inherits the momentum-space anisotropy and reduced rotational symmetry of its Fermi-lune parent

1. The TDCI state exhibits both in-plane and out-of-plane orbital magnetizations, which is established by the persistence of mixed orbital magnetization [45].

2. The in-plane orbital magnetization per unit area is calculated using the formula derived from non-reciprocal group velocities, linking it to layer polarization [59] [45].

3. The state exhibits anomalous Hall hysteresis under both perpendicular and in-plane magnetic fields, revealing a transdimensional orbital magnetic character [24] [45].

4. The minimal two-band model captures the essential orbital character of the TDCI by describing it with a Hamiltonian involving two effective layer orbitals and parameters like tilt amplitude t and band inversion mass m(k) [39] [45].

Collective Excitations and Stability

We further investigate neutral collective excitations to distinguish collectively stable TDCIs from regimes with competing finite-momentum density-wave instabilities

1. The Bethe–Salpeter equation is solved for particle–hole fluctuations about the self-consistent Hartree–Fock state [46] [54].

2. Collective modes are analyzed, including charge and spin modes, to distinguish stable TDCIs from regimes with competing finite-momentum density-wave instabilities [49] [53].

3. An in-plane magnetic field can harden the soft density-wave mode and may drive a transition from a density-wave state to TDCI [49] [53].

4. The two-band model provides a kinematic explanation for finite-momentum softening by showing that the tilt term can shift the conduction-band minimum and valence-band maximum to different momenta, lowering the continuum particle–hole excitation threshold at finite momentum [68] [45].

Fractional Filling and Fractional States

When the |C| = 1 band is fractionally filled, the many-body ground state may be a fractional Chern insulator with broken rotation and mirror symmetries

1. Exact diagonalization calculations demonstrate that when the |C| = 1 band is fractionally filled, the many-body ground state may be a fractional Chern insulator with broken rotation and mirror symmetries [45].

2. For twisted octalayer RMG/hBN at specific parameters, electron filling of 1/3 corresponds to a hole filling of 2/3 relative to filling 1, exhibiting three nearly degenerate low-energy states separated by a finite gap [4b] [45].

3. This threefold torus degeneracy is consistent with a Laughlin-type 1/3 FCI in the C = 1 band [15] [16].

4. The fractional state combines fractional Chern topology with the transdimensional orbital character of its parent state, suggesting it will exhibit in-plane orbital magnetization and magnetic hysteresis to in-plane magnetic field [45].

Experimental Signatures and Outlook

Experimentally, quantized anomalous Hall transport accompanied by Hall hysteresis under both perpendicular and inplane magnetic fields would provide a characteristic signature of the integer TDCI

1. Quantized anomalous Hall transport accompanied by Hall hysteresis under both perpendicular and in-plane magnetic fields is proposed as a characteristic signature of the integer TDCI [45].

2. The reduced symmetry allows for edge modes on opposite boundaries to have different group velocities, offering further signatures of the underlying orbital order [45].

3. Proximity to a superconductor may provide a route to chiral topological superconductivity with in-plane-field tunability, provided that effective pairing can be induced in the spin–valley-polarized bands [45].

The results establish a route from symmetry-broken orbital magnetism in a metallic parent state to quantized Chern phases in RMG moir´e systems. The TDCI retains the reduced symmetry and coexisting in-plane and out-of-plane orbital magnetizations of its metallic parent [45]. A minimal two-band model links the in-plane orbital magnetization to non-reciprocal group velocities and illustrates how the asymmetric dispersion can favor low-energy particle–hole excitations at finite momentum [45]. Our collective-mode calculations reveal competing density-wave instabilities in certain parameter regimes. An in-plane magnetic field can harden the soft charge mode, suggesting a means of stabilizing the TDCI against these fluctuations and potentially driving transitions between the TDCI and competing density-wave states [49] [53]. Fractional filling of the isolated C = 1 band further supports an FCI that is expected to inherit the transdimensional orbital character of its parent state. Experimentally, quantized anomalous Hall transport accompanied by Hall hysteresis under both perpendicular and inplane magnetic fields would provide a characteristic signature of the integer TDCI [45]. Enhanced finite-momentum charge or spin fluctuations could reveal nearby densitywave instabilities [45]. The reduced symmetry also allows edge modes on opposite boundaries to have different group velocities, offering further signatures of the underlying orbital order (see Fig. S4 of the Supplemental Material [45]). Finally, proximity to a superconductor may provide a route to chiral topological superconductivity with in-plane-field tunability, provided that effective pairing can be induced in the spin–valley-polarized bands [45].

--- Page 18 ---

The paper also reports quantitative values for the in-plane orbital magnetization My per moir´e unit cell, in µB per cell, for symmetry-broken Hartree–Fock (HF) states including TDCI. At ν = 1, one additional electron occupies each cell, so the same value gives the moment per added carrier [45]. The zero-field method and weak-field one-shot method yield quantitatively consistent values of the in-plane orbital magnetization [45].

--- Page 22 ---

The two-band model illustrates a mechanism for finite-momentum charge-mode softening. For a fixed spin and valley µ, the energy cost of an unbound particle–hole pair excitation with momentum Q in magnetic domain i is Eph i(k, Q) = E(k + Q) + E(k) + tnˆi · [F(k + Q) − F(k)] [45]. The tilt term can shift the conduction-band minimum and valence-band maximum to different momenta, lowering the continuum particle–hole excitation threshold mink Eph i(k, Q) at finite Q. Residual particle–hole interactions can further lower a collective branch below this continuum and drive it to zero while the single-particle gap remains open [45]. The non-reciprocal dispersion thus provides a kinematic route to a finite-momentum chargedensity-wave instability, whereas the actual soft wavevector and instability threshold additionally depend on the interaction kernel and orbital form factors [45].

--- Page 23 ---

The ribbon spectrum of the two-band model shows that the opposite edges show counterpropagating branches crossing the bulk gap, with markedly different dispersions in the selected domain [45]. This establishes a coupling of the excitation spectrum to the in-plane orbital field, although the sign of a collectivemode shift also depends on its wavefunction and the field dependence of Ki [45].

--- Page 24 ---

The opposite edges show counterpropagating branches crossing the bulk gap, with markedly different dispersions in the selected domain [45]. This establishes a coupling of the excitation spectrum to the in-plane orbital field, although the sign of a collectivemode shift also depends on its wavefunction and the field dependence of Ki [45].

--- Page 21 ---

The last equality uses the even parity of omega−,z and the odd parity of the tilt. At its special symmetric point EF = 0, Mz vanishes even when C− ≠ 0; a Chern number alone does not require a nonzero magnetization at every chemical potential [45].

--- Page 21 ---

The last equality uses the even parity of omega−,z and the odd parity of the tilt.

1. Improved AI can predict phase transitions in moiré superlattices by utilizing Hartree–Fock phase diagrams of TDCI state in RMG/hBN moir´e superlattices (Fig. 2), allowing it to map the emergence of the TDCI phases for NL = 6–9 as functions of displacement field D and twist angle θ.

2. Improved AI can analyze non-reciprocal transport properties by employing the minimal two-band model which unveils the intimate connection between non-reciprocal velocity and in-plane orbital magnetization, enabling it to calculate My = el/2Z mBZ d 2k/(2π) squared m(k) E(k) zˆ × δvi(k) (Eq. 59).

3. Improved AI can detect competing instabilities by solving the Bethe–Salpeter equation for neutral collective excitations, enabling it to identify a soft charge mode that softens to zero and becomes negative, signaling an energetic instability towards charge density wave order (Fig. 3b).

4. Improved AI can simulate fractional topological states by using exact diagonalization calculations demonstrate a fractional Chern insulator obtained from arXiv:2610.21310v1, allowing it to model the threefold low-energy manifold consistent with a Laughlin-type 1/3 FCI in a C = 1 band (Fig. 4b).

5. Improved AI can predict magnetic hysteresis by linking the reduced symmetry and mixed orbital magnetization of its metallic parent state to an observable effect, as the paper suggests that an in-plane magnetic field can harden the soft density-wave mode and may drive a transition from density-wave state to TDCI.