Superconductivity and competing orders in honeycomb $t$-$J$ model: interplay of lattice geometry and next-nearest-neighbor hopping

この論文では、ハニカム格子の拡張$t$-$J$モデルに対して大規模DMRG計算とスレーブボソン平均場理論を適用し、次 nearest neighbor ホッピング$t'$と境界幾何学が$d$波超伝導とストライプ秩序の競合をどのように制御するかを解明し、$t'$の導入により$d$波超伝導相が安定化される可能性を示唆しています。

要約

この論文は、**「電子が踊るハチの巣（六角形）のダンスフロア」**で、どんな条件下で「超電導（電気抵抗ゼロの魔法）」が起きるのか、そしてそのダンスの邪魔をする「整列した列（ストライプ秩序）」とどう戦っているのかを解明した研究です。

研究者たちは、巨大なコンピューターシミュレーション（DMRG）と、理論的な近似計算（SBMFT）という 2 つの異なる「レンズ」を使って、この複雑な現象を覗き見しました。

以下に、専門用語を排して、わかりやすい比喩で説明します。

1. 舞台設定：ハチの巣と電子のダンス

まず、ハチの巣（六角形）の格子を想像してください。これは、グラフェンや「ツイストされた二層グラフェン」といった、最近注目されている物質の構造です。

• 電子：このダンスフロアを動き回る「ダンサー」です。
• 半導体状態：ダンサーたちが固まって動けない状態（絶縁体）。
• ドープ（不純物添加）：ダンサーの一部を抜いて、少し隙間を作ること。これでダンサーたちは動き出せます。

通常、この隙間を作ると、ダンサーたちは「超電導（みんなで手を取り合って滑らかに動く状態）」になりたがりますが、同時に「ストライプ秩序（壁を作って整列して並ぶ状態）」という別の動きも競い合います。どちらが勝つかは、ダンサーたちの「隣の人との距離感（ hopping）」によって決まります。

2. 重要な要素：「次の隣人」への飛び方（$t'$）

この研究で注目したのは、「次の隣人（Next-Nearest Neighbor）」への飛び方です。

• 通常の飛び方（$t$）：すぐ隣の席へ移動する。
• 次の隣人への飛び方（$t'$）：斜め向かいの席へジャンプする。

この「斜めへのジャンプ（$t'$）」の強さを調整することで、ダンスフロアの雰囲気がどう変わるかを探りました。

3. 発見その 1：最適な「ジャンプ力」がある

シミュレーションの結果、驚くべきことがわかりました。

• ジャンプ力が弱いとき：ダンサーたちは少し固まって動き、超電導はあまり強くない。
• ジャンプ力が強すぎるとき：ダンサーたちはバラバラになり、整列（ストライプ）が強すぎて、超電導が邪魔される。
• 絶妙な「ジャンプ力」（$t' \approx 0.4$）：ここが**「黄金点」**です。この強さのジャンプ力があるとき、超電導が最も活発になります。

まるで、ダンスのテンポが速すぎず遅すぎず、**「ちょうどいいリズム」**で音楽が流れていると、ダンサーたちが最高のパフォーマンス（超電導）を見せるようなものです。

4. 発見その 2：「部屋の形」で勝敗が決まる（境界の重要性）

ここがこの論文の最も面白い点です。ダンサーたちのいる「部屋（シミュレーションに使われる円筒形のモデル）」の形を変えると、結果がガラリと変わりました。

• タイプ A の部屋（YC4-0）：壁の形が「アームチェア（肘掛け椅子）」の方向に沿っている。
  • 結果：超電導が強く、少しだけ整列（ストライプ）も混ざっているが、超電導が勝利。
• タイプ B の部屋（XC8-0）：壁の形が「ジグザグ」の方向に沿っている。
  • 結果：超電導は消えてしまい、整列（ストライプ）が完全勝利。

これは、**「ダンスフロアの壁の形が、ダンサーたちの動きを制限してしまう」**ことを意味します。同じルール（同じ物質）でも、容器の形（境界条件）によって、超電導が起きるか、整列が起きかが変わってしまうのです。

5. 2 次元の世界（本当の物質）ではどうなる？

シミュレーションは「細長い部屋（円筒）」で行いましたが、実際の物質は「広い平らな床（2 次元）」です。
研究者たちは、この「部屋の形による影響」を理論計算で補正し、広い床での真の姿を推測しました。

• 結論：広い床では、**「アームチェア型の整列（a-stripe）」が基本的な状態ですが、ジャンプ力（$t'$）が強くなると、「整列が溶けて、均一な超電導」**に変わる可能性があります。
• つまり、**「ジャンプ力（$t'$）を調整すれば、超電導を安定して実現できる」**という希望が持てます。

まとめ：この研究が意味すること

この研究は、**「ハチの巣構造を持つ物質」**において、超電導を発生させるための「魔法のレシピ」を見つけました。

1. 斜めへのジャンプ（$t'$）が重要：ただ隣に移動するだけでなく、斜めへの移動を少し加えることで、超電導が強化されます。
2. 最適な強さがある：ジャンプ力は「強すぎず、弱すぎず」がベストです。
3. 形状の罠：実験やシミュレーションでは、物質の形（境界）によって見かけの現象が変わるため、注意が必要ですが、それを逆手に取って異なる現象を調べることができます。

これは、**「ツイストグラフェン」**のような新しい物質で、より高い温度で超電導を実現するための道しるべとなる可能性があります。まるで、ダンサーたちの「ジャンプのタイミング」を微調整することで、魔法のような滑らかな動き（超電導）を引き出せるようになったようなものです。

技術概要

以下は、提示された論文「Superconductivity and competing orders in honeycomb t-J model: interplay of lattice geometry and next-nearest-neighbor hopping（ハニカム格子 t-J モデルにおける超伝導と競合秩序：格子幾何学と次近接ホッピングの相互作用）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

強相関電子系におけるドープされたモット絶縁体の物理学は、現代凝縮系物理学の中心的な課題の一つです。特に、高温超伝導体（銅酸化物）の理解を深めるために、正方格子におけるハバードモデルや t-J モデルは広く研究されてきました。しかし、最近のツイストド・グラフェンや遷移金属ダイカルコゲナイド（TMD）などのモアレ超格子系の実験的進展により、ハニカム格子における電子相関の理解が急務となっています。

ハニカム格子の半充填状態では、反強磁性（AFM）モット相、量子スピン液体（QSL）、価電子結合結晶（VBS）など、多様な基底状態が予言されています。ドープを導入すると、ホールの運動エネルギーと電子間相互作用の競合が生じ、超伝導（SC）と電荷密度波（CDW）秩序（ストライプ秩序など）のどちらが基底状態となるかが議論の的となっています。
特に、**次近接ホッピング項（$t'$）**および対応する超交換相互作用（$J'$）が、超伝導の安定性や CDW 秩序との競合にどのような役割を果たすかは、ハニカム格子 t-J モデルにおいて未解明な部分が多く残されています。また、従来の研究では有限サイズのシミュレーション結果を 2 次元熱力学的極限へ外挿する際に、境界条件や格子幾何学の影響を十分に考慮できていないという課題がありました。

2. 手法 (Methodology)

本研究では、以下の 2 つの相補的な手法を組み合わせることで、ハニカム格子上の拡張 t-t'-J-J' モデルを系統的に調査しました。

1. 大規模密度行列再正規化群（DMRG）シミュレーション:

   • 有限幅の円筒（cylinder）幾何学を用いて、基底状態の性質を高精度に計算しました。
   • 境界条件の影響を評価するために、3 種類の異なる円筒幾何学（YC4-0, YC4-4, XC8-0）を使用しました。これらは、周期的境界条件を適用する方向と、格子結合の向き（armchair 方向または zigzag 方向）が異なります。
   • システムサイズは $N=384$ サイト（ドープ濃度 $\delta = 1/8, 1/12$）、ブロック状態数 $m=24000$ まで保持し、収束性を確保しました。
   • 超伝導相関関数（ペア - ペア相関）と電荷密度分布を解析し、ルッティング指数（$K_{sc}, K_c$）を抽出して秩序の長距離性を評価しました。

2. スレーブボソン平均場理論（SBMFT）:

   • DMRG で観測された競合する秩序状態（a-ストライプ、z-ストライプ、双方向 CDW など）の 2 次元熱力学的極限における安定性を評価するために採用しました。
   • 電子をフェルミオン的スピンオンとボソニックホールンに分解し、粒子 - 粒子チャネルおよび粒子 - hole チャネルで平均場近似を行いました。
   • 異なる初期電荷分布パターンから出発し、エネルギーが最小となる安定な相を特定しました。

3. 主要な結果 (Key Results)

A. YC4-0 円筒におけるロバストな超伝導相

YC4-0 円筒（armchair 方向に結合を持つ）上での DMRG 結果は、広範な $t'$ の範囲で**準長距離的な nematic d 波超伝導（SC）**が支配的であることを示しました。

• SC の最適化: 超伝導ルッティング指数 $K_{sc}$ は $t'$ に対して非単調な「谷型」依存性を示します。$t'_{op} \approx 0.4$（$t=2$ の場合）で $K_{sc}$ が最小（$\approx 0.82$）となり、超伝導が最も強化されます。これは、未ドープ相の AFM 領域（$J' \sim 0.02J$）に位置しており、量子臨界点や QSL 領域ではなく、ドープされた AFM 状態内のスピン揺らぎが超伝導を駆動している可能性を示唆しています。
• CDW との共存: 低 $t'$ 領域では、SC とともに「armchair 方向のストライプ（a-stripe）」が副次的な秩序として共存します。しかし、$t' > 0.5$ になると、$\delta=1/8$ の場合、CDW 相関が著しく抑制され、一様な電荷密度を持つ nematic d 波 SC 相への移行が観測されました。

B. 幾何学に依存する秩序の競合

異なる円筒幾何学（XC8-0, YC4-4）では、同じモデルパラメータでも異なる基底状態が観測されました。

• XC8-0 円筒: $t' > 0.5$ の領域で、**長距離的な zigzag 方向のストライプ（z-stripe）**が安定し、超伝導は消失します。
• YC4-4 円筒: $t'=0$ 付近では SC と CDW が共存する z-ストライプが観測されますが、$t'$ が増加すると位相分離領域に入ります。
• この結果は、有限幅シミュレーションにおける秩序状態の安定性が、境界幾何学に強く依存することを示しており、2 次元極限への外挿の難しさを浮き彫りにしました。

C. 2 次元極限における SBMFT の予測

SBMFT による 2 次元極限の相図は、DMRG の幾何学依存性を補完する重要な知見を提供しました。

• $\delta = 1/12$: 広範な $t'$ 領域において、a-stripe 秩序が他の競合状態（z-stripe など）よりもエネルギー的に安定であることが示されました。
• $\delta = 1/8$: 小 $t'$ 領域では a-stripe が支配的ですが、$t' > 0.5$ 付近で一様な電荷密度を持つ nematic d 波 SC 相への転移が起こります。これは、YC4-0 円筒で観測された CDW の抑制傾向と一致しています。

4. 貢献と意義 (Contributions and Significance)

1. $t'$ 誘起超伝導の確立:
   ハニカム格子のドープされた拡張 t-J モデルにおいて、次近接ホッピング $t'$ が超伝導を強化し、CDW 秩序を抑制するメカニズムを初めて包括的に実証しました。特に、最適 $t'$ 付近で AFM 揺らぎが超伝導 pairing の主要な駆動力となる可能性を指摘しました。

2. 幾何学効果の解明と手法の統合:
   有限サイズシミュレーション（DMRG）における境界条件の重要性を明確にし、それを「秩序の競合を制御する新しいパラメータ」として捉える戦略を提案しました。DMRG（準 1 次元系での高精度計算）と SBMFT（2 次元極限での相図評価）を組み合わせることで、複雑な競合秩序を持つ 2 次元強相関系を解明する有効なアプローチを示しました。

3. 実験への示唆:
   ツイストド・モアレ超格子（例：twisted bilayer graphene や TMD 類）では、ホッピングパラメータ $t'$ が実験的に制御可能です。本研究は、$t' \approx 0.4t$ 付近の領域で非従来型超伝導が実現しやすいことを示唆しており、これらの物質系における超伝導探索の指針となります。

結論

本論文は、ハニカム格子におけるドープされたモット絶縁体の基底状態が、ドープ濃度、次近接ホッピング、そして格子幾何学の微妙なバランスによって決定されることを示しました。特に、$t'$ の導入が nematic d 波超伝導を安定化させるという発見は、高温超伝導のメカニズム理解と、新しいモアレ超格子材料における超伝導制御の両面で重要な意義を持ちます。