Magic-angle twisted bilayer graphene (MATBG) hosts normal-state nearly-flat bands with nonzero Euler numbers and shows superconductivity. In this work, we study the effects of the nontrivial normal-state band topology on the intervalley $C_{2z}\mathcal{T}$-invariant mean-field Cooper pairing order parameter in MATBG. We show that the pairing order parameter can always be split into a trivial channel and an Euler obstructed channel in all gauges for the normal-state basis, generalizing the previously-studied channel splitting in the Chern gauge. The nonzero normal-state Euler numbers require the pairing gap function of the Euler obstructed channel to have zeros, while the trivial channel can have a nonvanishing pairing gap function. When the pairing is spontaneously nematic, we find that a sufficiently-dominant Euler obstructed channel with two zeros typically leads to nodal superconductivity. Under the approximation of exactly-flat bands, we find that the mean-field zero-temperature superfluid weight is generally bounded from below, no matter whether the Euler obstructed channel is dominant or not, generalizing the previously-derived bound for the uniform s-wave pairing. We numerically verify these statements for pairings derived from a local attractive interaction. Our work suggests that Euler obstructed Cooper pairing may play an essential role in the superconducting MATBG.
論文ID : 2202.02353タイトル : Euler Obstructed Cooper Pairing in Twisted Bilayer Graphene: Nematic Nodal Superconductivity and Bounded Superfluid Weight著者 : Jiabin Yu, Ming Xie, Fengcheng Wu, Sankar Das Sarma分類 : cond-mat.supr-con(超伝導)、cond-mat.mes-hall(メゾスコピック物理)発表時期 : 2022年2月(arXiv プレプリント)論文リンク : https://arxiv.org/abs/2202.02353 本論文は、魔角捻れ二層グラフェン(MATBG)における非自明な正常状態バンド位相相がクーパー対形成超伝導に及ぼす影響を研究している。著者らは、対形成秩序パラメータが自明なチャネルとオイラー阻害チャネルに分解でき、後者は非ゼロのオイラー数のため必ずゼロ点を持つことを発見した。対形成が自発的にネマティック化する場合、支配的なオイラー阻害チャネルは節点超伝導をもたらす。厳密な平坦バンド近似下では、ゼロ温度超流動重みは普遍的な下界を有する。
MATBGの特異性 : 魔角捻れ二層グラフェンは電荷中性点付近でほぼ平坦なバンドを有し、これらのバンドは非ゼロのオイラー数を持ち、非自明な位相的性質を特徴付ける。超伝導現象の観測 : MATBGは部分充填時に超伝導性を示し、実験的にネマティック(向列)および節点超伝導の特性が観測されている。理論的課題 : 既存理論は、正常状態の位相的性質がいかに超伝導対形成に影響するか、特にネマティック節点超伝導の形成機構を完全には理解していない。正常状態のオイラー数と超伝導対形成の間の直接的な関連性を確立する ネマティック超伝導がいかに節点超伝導をもたらすかを理解する 超流動重みの理論的下界を提供する ゲージ不変なオイラー阻害クーパー対形成理論を確立 し、対形成秩序パラメータを自明なチャネル(Δ⊥)とオイラー阻害チャネル(Δ∥)に分解したオイラー阻害チャネルは必ずゼロ点を持つことを証明 し、ゼロ点の総巻数は正常状態のオイラー数により決定される:
∑ i W b , i = N + − ( − 1 ) b N − \sum_i W_{b,i} = N_+ - (-1)^b N_- ∑ i W b , i = N + − ( − 1 ) b N − ネマティック超伝導が節点超伝導をもたらすメカニズムを解明 :オイラー阻害チャネルが支配的で2つのゼロ点を有する場合、自発的ネマティック対形成は安定な節点超伝導を生成する超流動重みの普遍的下界を導出 :
Tr [ D S F ] ≥ ⟨ [ ∣ Δ ⊥ ( k ) ∣ − ∣ Δ ∥ ( k ) ∣ ] 2 [ ∣ Δ ⊥ ( k ) ∣ − ∣ Δ ∥ ( k ) ∣ ] 2 + μ 2 ⟩ g 4 e 2 π N + \text{Tr}[D_{SF}] \geq \left\langle\frac{[|\Delta_\perp(k)| - |\Delta_\|(k)|]^2}{\sqrt{[|\Delta_\perp(k)| - |\Delta_\|(k)|]^2 + \mu^2}}\right\rangle_g \frac{4e^2}{\pi}N_+ Tr [ D SF ] ≥ ⟨ [ ∣ Δ ⊥ ( k ) ∣ − ∣ Δ ∥ ( k ) ∣ ] 2 + μ 2 [ ∣ Δ ⊥ ( k ) ∣ − ∣ Δ ∥ ( k ) ∣ ] 2 ⟩ g π 4 e 2 N + ウィルソン線とゲージ不変演算子に基づき、著者らは以下を定義した:
P Δ ( k ) = ∣ u + , k ⟩ Δ ( k ) ⟨ u − , − k C 2 z T ∣ P_{\Delta}(k) = |u_{+,k}\rangle\Delta(k)\langle u^{C_{2z}T}_{-,-k}| P Δ ( k ) = ∣ u + , k ⟩ Δ ( k ) ⟨ u − , − k C 2 z T ∣
ゲージ不変な量子幾何テンソルを導入することで:
Q α ( k ) = N α , k 0 N α Q α , k 0 ( k ) Q_\alpha(k) = \frac{N_{\alpha,k_0}}{N_\alpha}Q_{\alpha,k_0}(k) Q α ( k ) = N α N α , k 0 Q α , k 0 ( k )
対形成秩序パラメータのチャネル分解を実現した:
Δ ( k ) = Δ ⊥ ( k ) + Δ ∥ ( k ) \Delta(k) = \Delta_\perp(k) + \Delta_\|(k) Δ ( k ) = Δ ⊥ ( k ) + Δ ∥ ( k )
オイラー阻害チャネルに対して、ゼロ点の巻数は以下を満たす:
N + − ( − 1 ) b N − ≠ 0 N_+ - (-1)^b N_- \neq 0 N + − ( − 1 ) b N − = 0 の場合、Δ b ( k ) \Delta_b(k) Δ b ( k ) は必ずゼロ点を持つゼロ点の総巻数はN + − ( − 1 ) b N − N_+ - (-1)^b N_- N + − ( − 1 ) b N − に等しい H ( k ) = ( h + ( k ) − μ Δ ⊥ ( k ) + Δ ∥ ( k ) [ Δ ⊥ ( k ) + Δ ∥ ( k ) ] † − h + T ( k ) + μ ) H(k) = \begin{pmatrix}
h_+(k) - \mu & \Delta_\perp(k) + \Delta_\|(k) \\
[\Delta_\perp(k) + \Delta_\|(k)]^\dagger & -h_+^T(k) + \mu
\end{pmatrix} H ( k ) = ( h + ( k ) − μ [ Δ ⊥ ( k ) + Δ ∥ ( k ) ] † Δ ⊥ ( k ) + Δ ∥ ( k ) − h + T ( k ) + μ )
系はCIノード類に属し、カイラル対称性を有する:
C ( k ) H ( k ) C † ( k ) = − H ( k ) C(k)H(k)C^\dagger(k) = -H(k) C ( k ) H ( k ) C † ( k ) = − H ( k )
節点超伝導は非ゼロのカイラル数により保護され、det [ h + ( k ) − μ − i Δ ( k ) U + ∗ ( k ) ] \det[h_+(k) - \mu - i\Delta(k)U_+^*(k)] det [ h + ( k ) − μ − i Δ ( k ) U + ∗ ( k )] のゼロ点を分析することで決定される。
厳密な平坦バンド近似下では、ゼロ温度超流動重みは:
[ D S F ] i j = − 2 e 2 V ∑ k , a ∂ A i ∂ A j μ 2 + λ + , k , a ( A ) ∣ A → 0 [D_{SF}]_{ij} = -\frac{2e^2}{V}\sum_{k,a} \partial_{A_i}\partial_{A_j}\sqrt{\mu^2 + \lambda_{+,k,a}(A)}\Big|_{A \to 0} [ D SF ] ij = − V 2 e 2 ∑ k , a ∂ A i ∂ A j μ 2 + λ + , k , a ( A ) A → 0
ここでλ + , k , a ( A ) = ( ∣ Δ ⊥ ( k , A ) ∣ + ( − 1 ) a − 1 ∣ Δ ∥ ( k , A ) ∣ ) 2 \lambda_{+,k,a}(A) = (|\Delta_\perp(k,A)| + (-1)^{a-1}|\Delta_\|(k,A)|)^2 λ + , k , a ( A ) = ( ∣ Δ ⊥ ( k , A ) ∣ + ( − 1 ) a − 1 ∣ Δ ∥ ( k , A ) ∣ ) 2 である。
捻れ角 : θ ∈ [ 1.05 ° , 1.15 ° ] \theta \in [1.05°, 1.15°] θ ∈ [ 1.05° , 1.15° ] BMモデルパラメータ : w 0 / w 1 = 0.8 w_0/w_1 = 0.8 w 0 / w 1 = 0.8 、v 0 = 5817 v_0 = 5817 v 0 = 5817 meV·Å、w 1 = 110 w_1 = 110 w 1 = 110 meV充填 : 1モアル単位胞あたり2.5個のホール相互作用強度 :
副格子内対形成: g 0 = 40 g_0 = 40 g 0 = 40 meV·(nm)² 副格子間対形成: g 1 = 70 g_1 = 70 g 1 = 70 meV·(nm)² 自己無撞着求解 : 平均場方程式を数値的に求解して対形成秩序パラメータを取得位相的解析 : 対形成ゼロ点の巻数を計算して理論予測を検証エネルギースペクトル計算 : ボゴリューボフ・ド・ジェンヌハミルトニアンのエネルギーギャップ構造を分析副格子内対形成 : 主な寄与は自明なチャネル(Δ⊥が支配的)から来ている副格子間対形成 : 主な寄与はオイラー阻害チャネル(Δ∥が支配的)から来ている副格子間対形成に対して:
Δ ∥ ( k ) \Delta_\|(k) Δ ∥ ( k ) はΓ M \Gamma_M Γ M 点で巻数1の2つのゼロ点を有するボゴリューボフ・ド・ジェンヌエネルギーギャップはゼロであり、節点超伝導の存在を確認した 比率r ∥ ⊥ = ⟨ ∣ Δ ∥ ∣ ⟩ / max ( ∣ Δ ⊥ ∣ ) > 1 r_{\|⊥} = \langle|\Delta_\||\rangle/\max(|\Delta_\perp|) > 1 r ∥⊥ = ⟨ ∣ Δ ∥ ∣ ⟩ / max ( ∣ Δ ⊥ ∣ ) > 1 はオイラー阻害チャネルが支配的であることを示す 数値計算によるTr [ D S F ] ≈ 10 8 \text{Tr}[D_{SF}] \approx 10^8 Tr [ D SF ] ≈ 1 0 8 H⁻¹ 理論的下界は同じ桁数の推定値を与える BKT転移温度T B K T ∈ [ 1 , 2 ] T_{BKT} \in [1,2] T B K T ∈ [ 1 , 2 ] Kは実験観測と一致している C 3 z C_{3z} C 3 z 不変対形成 : ゼロ点はK M K_M K M とK M ′ K'_M K M ′ に固定され、必ずしも節点超伝導を生成しない自発的ネマティック対形成 : 少なくとも1つのゼロ点がΓ M \Gamma_M Γ M にあり、微調整を除いて常に節点超伝導をもたらす本論文は3D半金属におけるオイラー阻害対形成の理論を2D モアレ系に拡張し、ゲージ不変な形式的記述を確立した。
既存のフォノン媒介、電子相関などのメカニズムと比較して、本論文は正常状態の位相的性質が超伝導特性に及ぼす根本的な影響を強調している。
均一なs波対形成に対する以前の下界を、一般的なC 2 z T C_{2z}T C 2 z T 不変対形成に拡張した。
普遍性 : オイラー阻害クーパー対形成はMATBGなどの位相的平坦バンド系における普遍的な現象であるメカニズム : ネマティックオイラー阻害対形成はMATBG節点超伝導を理解するための新しいメカニズムを提供する安定性 : ノードが存在しても、超流動重みは確定的な下界を有する平坦バンド近似 : 超流動重みの下界は厳密な平坦バンド近似に依存し、実際のMATBGでは完全には成立しない可能性がある平均場理論 : 量子揺らぎと強相関効果を無視しているパラメータ依存性 : 節点超伝導の出現は特定の相互作用パラメータ選択に依存するエネルギーバンド色分散が超流動重みに及ぼす補正を考慮する 他の対称性破れがオイラー阻害対形成に及ぼす影響を研究する オイラー阻害対形成を実験的に検証する方法を探索する 理論的革新性 : ゲージ不変なオイラー阻害対形成理論を確立し、重要な理論的価値を有する物理的洞察 : 位相と超伝導の間の深い関連性を解明し、MATBG超伝導の理解に新しい視点を提供する数学的厳密性 : 理論導出は厳密であり、数値計算は理論予測を十分に検証している実用的意義 : 予測されたBKT温度は実験観測と良好に一致している実験的検証の不足 : オイラー阻害対形成の存在を直接支持する実験的証拠が不足しているモデルの制限 : BMモデルは実際のMATBGの複雑性を完全に記述できない可能性があるパラメータ感度 : 結果は相互作用パラメータの選択に対して比較的敏感である本論文は位相材料における超伝導を理解するための重要な理論的枠組みを提供し、特にMATBGおよび他のモアレ超伝導体の研究に指導的意義を有する。理論予測の普遍性により、他の位相超伝導系にも適用可能である可能性がある。
非自明な位相を有する平坦バンド系 モアレ超伝導体(捻れた多層グラフェンなど) オイラー阻害を有する他の位相材料 論文は72篇の関連文献を引用しており、MATBGの実験的発見、理論モデル、位相超伝導理論など重要な研究を網羅し、研究に堅実な理論的基礎を提供している。