Interacting Chern insulator transition on the sphere: revealing the Gross-Neveu-Yukawa criticality

技術要約：球面上における相互作用するチャーン絶縁体転移

問題の所在
二次元における自明なバンド絶縁体からチャーン絶縁体への転移は、対称性の破れを伴わずにトポロジーが変化する、典型的な「ランドー理論を超えた（beyond-Landau）」臨界現象を象徴している。この低エネルギー物理は、単一のディラックコーンによって支配されている。非相互作用極限（$u_0=0$）は自由フェルミオン共形場理論（CFT）によって記述されるが、電子相互作用（$u_0 > 0$）は系をグロス・ネーヴ・ユーカワ（GNY）固定点へと駆動する。この多重臨界点は、フェルミオンと質量演算子が同時にギャップレスとなる点であり、高エネルギー物理学および物性物理学の両面において広く関心の高いベンチマーク問題である。

しかし、相互作用する単一ディラックコーンを忠実に数値的に扱うことは、主に以下の4つの障害により極めて困難である：

1. パリティ・アノマリー： 局所的で時間反転対称性を持つ格子レギュラリゼーションは、パリティ・アノマリーに苦しむことが多く、意図した普遍性クラスを不明瞭にするために、フレーバー数の増大や対称性の破れを必要とする場合がある。
2. 強相関： 臨界相互作用ディラックフェルミオンは、自由ディラックCFTからは摂動的に到達できない。なぜなら、四フェルミオン相互作用は自由な固定点において無次元（irrelevant）だからである。
3. 対称性の制限： $SO(3)$ 対称性が欠如している場合、ほとんどの微視的な計算は、全演算子の内容を決定するための状態・演算子対応を利用できず、秩序パラメータからのみ臨界指数を抽出することしかできない。
4. 分散の問題： $SO(3)$ を保存する標準的なファジー球構成は、通常、球の中心にあるモノポールによるランダウ準位を生成してしまい、対象とする線形ディラック分散を再現できない。

手法
著者らは、これらの障害を回避するために、ディラックコーンを球（$S^2$）の表面上に直接定式化するフレームワークを提案している。

• 幾何学的セットアップ： 半径 $R=1$ の球面上にディラック演算子を配置することで、必要なスピン接続を曲率によって供給する。ディラック演算子は $D = \sigma \cdot L + 1$ （ここで $L$ は軌道角運動量、$\sigma$ はスピン）の形をとる。このセットアップは、完全な回転 $SO(3)$ 対称性と離散対称性（時間反転 $T$、電荷共役 $C$、パリティ $P$）を保持する。
• ハミルトニアンの構成： 系はハミルトニアン $H = H_0 + H_{int}$ で記述される。相互作用項は、半径方向の局所スピン演算子 $\sigma_\Omega = \hat{n}_\Omega \cdot \sigma$ を用いることで、$SO(3)$ および離散対称性を保存する。相互作用ポテンシャルは、ハルダン・シュードポテンシャル $U_l$ を用いて展開され、これにより接触相互作用（$U_0$）や高次項（$U_1, U_2$）の制御が可能となる。
• 厳密対角化（ED）： 著者らは、有限の角運動量カットオフ $j_{max}$ を持つ系に対してEDを実行する。ヒルベルト空間の急速な増大を考慮し、本研究では $j_{max} = 3/2$ および $j_{max} = 5/2$（電子数が10個未満の系に対応）に焦点を当てる。
• スペクトル解析： $S^2 \times \mathbb{R}$ 上のエネルギー・スペクトルは、径方向量子化を用いて、CFT演算子のスケーリング次元 $\delta E = v \Delta$ へと写像される。速度 $v$ は、最低次スカラー演算子 $\sigma$ のスケーリング次元をコンフォーマル・ブートストラップの値（$\Delta_\sigma = 0.66$）に一致させることで較正される。シュードポテンシャル $(U_0, U_1, U_2)$ は、ディラックスピノル $\psi$ とそのデセンダントのスケーリング次元をブートストラップの予測値に一致させるように微調整される。

主要な貢献と結果

1. GNY臨界性の実現： 本研究は、球面上における控えめなヒルベルト空間において、堅牢なGNY臨界点が実現可能であることを示している。相互作用パラメータを調整することで、著者らは系が期待通りの共形対称性を示す臨界領域を特定している。
2. 共形スペクトルの抽出： EDスペクトルは、電荷セクター $q=0, 1, 2$ にわたる8つの異なる共形マルチプレットを明らかにする。質量演算子 $\sigma$、ディラックスピノル $\psi$、およびシングレット・ペアリング演算子 $\sigma_T$ の主要なスケーリング次元は、最先端のコンフォーマル・ブートストラップの境界値と数パーセント以内の誤差で一致している（例：$j_{max}=5/2$ における $\sigma_T$ の誤差は約3%）。
3. 新たなプライマリー演算子の発見： 著者らは、GNYスペクトルにおいてこれまで未発見であった2つのプライマリー演算子を特定した：
   • $\tilde{J}_\mu$：保存電流 $J_\mu$ と同じ量子数を持つが、微分デセンダントを許容する非保存電流。
   • $D_{\mu\nu}$：$q=2$ セクターにおけるスピン2の演算子であり、スケーリング次元 $\Delta \approx 3.03$（$j_{max}=5/2$ 時）を持つ。この値は、$(2+1)$次元CFTにおけるスピン-$\ell$ の演算子のユニタリティ境界 $\Delta = \ell + 1$ をほぼ飽和させている。
4. 弱く破れた高次スピン対称性の証拠： $D_{\mu\nu}$ によるユニタリティ境界の近飽和、およびより低いスピンのデセンダントの不在は、$D_{\mu\nu}$ がほぼ保存された電流として機能していることの直接的な数値的証拠である。著者らは、これが $p$ 波ペアリング演算子 $\psi^T R (\gamma_\mu \partial_\nu + \gamma_\nu \partial_\mu) \psi$ に対応すると推測しており、これはGNY CFTで予想される弱く破れた高次スピン対称性を表している。
5. 球面EDの検証： 低エネルギー励起が特殊共形群の表現を形成していることが確認され、3次元共形対称性の出現が検証された。この手法は、保存された演算子の制限されたデセンダント構造（$\partial_\mu J^\mu = 0$ が低スピンの微分を禁止する）およびスピノル・マルチプレットの特徴的な構造を正確に再現することに成功している。

意義
本論文は、曲がった幾何学上にディラックコーンを配置することが、強く相互作用するトポロジカル転移を研究するための「驚くほど単純かつ正確なセットアップ」を提供すると主張している。相互作用するディラックフェルミオンにファジー球の枠組みを一般化することで、このアプローチはGNY固定点のためのコンパクトで透明な数値的顕微鏡を提供する。これは、パリティ・アノマリーを自然に回避し、完全な演算子内容を抽出するために必要な完全な回転対称性を保持することにより、単一フレーバー・ディラックフェルミオンのシミュレーションにおける長年の困難を解決している。新たなプライマリーの特定と、弱く破れた高次スピン対称性の確認は、GNY普遍性クラスのCFTデータセットを豊かにし、$\epsilon$ 展開や大規模 $N$ 手法に依存していた理論的解析に対する独立した数値的ベンチマークを提供するものである。