[論文レビュー] Kekulé spiral order in magic-angle graphene: a density matrix renormalization group study

Tianle Wang, Daniel E. Parker+5 more|arXiv (Cornell University)|Nov 4, 2022

Graphene research and applications参考文献 75被引用数 6

ひとこと要約

本研究では、結合次元χ=24576を用いた大規模で偏りのない密度行列逐次縮約法（DMRG）計算を用い、魔性角ねじれビレイヤヤグラフェン（TBG）におけるヘテロストレイン効果を、ν=−3で調査した。0.05%のヘテロストレインで、量子化された異常ホール絶縁体から、モアレ周期と不整合な位相モードを持つ修正された並進対称性を保つが、谷および並進対称性を破る非周期的ケクレースピRAL（IKS）相へのストレイン駆動的転移が特定された。さらに高応力下では通常金属への転移が起こり、非整列系における理論と実験の不一致が解消された。

ABSTRACT

When the two layers of a twisted moiré system are subject to different degrees of strain, the effect is amplified by the inverse twist angle, e.g., by a factor of 50 in magic angle twisted bilayer graphene (TBG). Samples of TBG typically have heterostrains of 0.1-0.7%, increasing the bandwidth of the "flat'' bands by as much as tenfold, placing TBG in an intermediate coupling regime. Here we study the phase diagram of TBG in the presence of heterostrain with unbiased, large-scale density matrix renormalization group calculations (bond dimension $χ=24576$), including all spin and valley degrees of freedom. Working at filling $ν= -3$, we find a strain of $0.05\%$ drives a transition from a quantized anomalous Hall insulator into an incommensurate-Kekulé spiral (IKS) phase. This peculiar order, proposed and studied at mean-field level by Kwan et al (PRX 11, 041063), breaks both valley conservation and translation symmetry $\hat{T}$, but preserves a modified translation symmetry $\hat{T}'$ with moiré-incommensurate phase modulation. Even higher strains drive the system to a fully symmetric metal.

研究の動機と目的

• 強相関理論が予測するν=−3における量子化された異常ホール（QAH）絶縁体と、非整列TBG試料における実験的観測（C=0の絶縁体または金属）との不一致を解消すること。
• グラフェンスケールでは微小なものの、逆ねじれ角度によって増幅される実際のヘテロストレインが、TBGの低エネルギー物理学にどのように影響を与えるかを調査すること。
• すべてのスピンおよび谷自由度を含む偏りのない大規模DMRGシミュレーションを用いて、ν=−3におけるTBGの基底状態の性質を特定すること。
• バンド幅がストレインによって顕著に拡大される中間結合領域におけるTBGの相図を明確にすること。強相関近似では到達できない領域である。

提案手法

• ν=−3におけるねじれビレイヤヤグラフェンのハーバードハミルトニアンを、結合次元χ=24576の大規模DMRG計算でシミュレートした。
• 多体相関を偏りなく取り扱うために、DMRGシミュレーションにすべてのスピンおよび谷自由度を組み込んだ。
• 2つのグラフェン層にヘテロストレインεGr=0.1–0.7%を適用し、逆ねじれ角度によって増幅されることで、モアレスーパーラティスが顕著に歪み、フラットバンドのバンド幅が増加した。
• DMRG基底状態エネルギーおよびオーダーパラメータ（例：谷およびスピン密度波、量子もつれエントロピー）を用いて、異なるストレイン値における相の識別を行った。
• 量子数（τz, sz）の分析により、相の分類と対称性の破れパターン（特に非周期的ケクレースピRAL（IKS）秩序）を特定した。
• 異なる量子数セクター間の基底状態エネルギーを比較し、グローバル基底状態および相境界を同定した。

実験結果

リサーチクエスチョン

• RQ1 グラフェンスケールでは微小なにもかかわらず、ヘテロストレインが魔性角ねじれビレイヤヤグラフェンの多体基底状態にどのように影響を与えるか？
• RQ2 非整列TBG試料における実験的に観測されたν=−3の絶縁体状態の性質は何か？また、現実的なストレイン効果で説明可能か？
• RQ3 もともと平均場論で予想された非周期的ケクレースピRAL（IKS）秩序が、偏りのないDMRGシミュレーションでも安定な相として出現するか？
• RQ4 ストレインが増加するに従い、TBGの相図は強相関領域から中間結合領域へどのように変化するか？
• RQ5 IKS相における対称性の破れの役割は何か？また、従来の量子ホールフェロ磁体とどのように異なるか？

主な発見

• ヘテロストレインが0.05%に達する段階で、量子化された異常ホール（QAH）絶縁体から非周期的ケクレースピRAL（IKS）相への転移が発生した。
• IKS相は谷保存則および並進対称性（T）を破るが、モアレ周期と不整合な位相モードを持つ修正された並進対称性T′を保つ。
• より高いストレイン（εGr ≥ 0.2%）では、系は完全に対称な通常金属に転移し、(0,0)量子数セクターがグローバル基底状態となった。
• IKS相はスピン非対称な間接谷相関および一時的なスピン密度波（SDW）秩序を示すが、ネットスピン極化は観測されなかった。
• 相間のエネルギー差は小さいが、数値精度を上回っており、QAH（低ストレイン）→ IKS（中間ストレイン）→ 通常金属（高ストレイン）という明確な三相構造が安定であることを示唆している。
• 相境界は系サイズおよび誘電率に敏感であるが、全体的な三相構造は予想される安定性を示している。