Maximal Orbital Entanglement from Non-Interacting Fermions

技術要約：非相互作用フェルミオンからの最大軌道もつれ

問題提起
量子多体系物理学および量子化学において、軌道回転は通常、量子状態の複雑性を少数の関連する軌道へと圧縮するために（テンソルネットワーク法やハートリー＝フォック理論などで）、軌道相関を最小化する目的で用いられる。逆に、量子情報処理（QIP）は、リソースとしてもつれ（エンタングルメント）を最大化しようとする。自由状態（非相互作用フェルミオン）における軌道相関の最小化は十分に理解されている（それらの自然軌道基底において消失する）が、それらの軌道もつれの最大化については、未だほとんど探索されておらず、特定の状態、平均的な性質、あるいは区別可能な粒子やボゾン系に対してのみ議論されてきた。

本研究は、根本的な問いに取り組む。非相互作用フェルミオン（具体的にはスレーター行列式および一般的な自由状態）は、最適な軌道回転の下で、どの程度の軌道もつれおよび全相関を保持できるのか？

手法
著者らは、$N$フェルミオン状態における軌道相関の最大化問題を、1粒子簡約密度行列（1RDM）の空間におけるスペクトル最適化問題へと再定式化した。

1. 1RDMへの還元： 自由状態（スレーター行列式やガウス状態を含む）は、その1RDMによって完全に決定されるため、多体ヒルベルト空間における最適化は、$d \times d$のエルミテート行列（ここで$d$は格子サイト数）に対する最適化へと還元される。
2. ブロックスペクトルと主要化（Majorization）： 問題はさらに、格子の分割によって誘導されるブロックスペクトル（1RDMの対角ブロックの固有値）の集合を特徴付ける問題へと還元される。著者らは、最適のブロックスペクトルを特定するために、主要化（$\prec$）の概念を利用する。彼らは、目的関数が厳密にシューア凹（Schur-concave）であることを確立している。したがって、最大化因子は、許容されるブロックスペクトルの集合の中で、主要化において最小となる（すなわち、集合内の他のすべての要素によって主要化される）要素に対応する。
3. 新しい数学的ツール： 一般的な問題を解くために、著者らは、行列のブロックのスペクトルに関する新しい主要化関係（定理1および定理4）を導出し、これによって行列のスペクトルと、その対角ブロックのスペクトルとの関係を規定する。

主な貢献と結果

• 結果1：スレーター行列式の最大もつれ
  $d = d_A + d_B$ 個のサイトを持つ格子上の $N$ フェルミオン系を、2つの領域に分割した場合、スレーター行列式（SD）によって達成可能な最大軌道エンタングルメント・エントロピーは以下の通りである：
  $$E_{\text{max}} = \min\{d_A, d_B, N, d-N\} \times \log 2$$
  これは、軌道ペアの「ビームスプリッター型」の回転によって達成される（具体的には、軌道 $j$ と $d+1-j$ の混合）。粒子数と正孔数が十分な場合（$\min\{N, d-N\} \ge \min\{d_A, d_B\}$）、この値はあらゆるフェルミオン状態において可能な最大値に等しく、非相互作用状態がQIPにおいて相互作用状態と同等のリソースとなり得ることを示唆している。

• 結果2：任意の自由状態における最大相関（二部分割）
  著者らは、結果1を、与えられた自然占有数 $n = (n_1, \dots, n_d)$ を持つ任意の自由状態（混合状態を含む）へと一般化する。任意の二部分割 $\lambda = (\lambda_1, \lambda_2)$ に対して、ブロックスペクトルの最大化因子 $b_{\text{max}}$ の閉形式を導出する。$b_{\text{max}}$ の成分は、自然占有数のペアの平均値（例：$(n_1+n_d)/2, (n_2+n_{d-1})/2, \dots$）である。最大全相関は $C_{\text{max}} = H(b_{\text{max}}) - H(n)$ となる。

• 結果3：任意の純粋状態に対するエンタングルメント境界
  導出された主要化関係を用いることで、著者らは、任意のパリティ超選択規則（SSR）を満たすフェルミオン純粋状態（相互作用のある基底状態を含む）の軌道エンタングルメントに対する上限を、その自然占有数を用いて確立する：
  $$E(|\Psi\rangle) \le \frac{1}{2} H(b_{\text{max}})$$
  この境界は、標準的な次元に基づく境界よりも厳密に強いものである。

• 結果4：最も細かい分割（自明な多部分割）
  各格子サイトが独立した部分系となる最も細かい分割（$\lambda = (1, \dots, 1)$）では、最大化因子は、すべての成分が平均充填率 $\eta = (\sum n_j)/d$ に等しい、最大混合ベクトルとなる。この最適ユニタリは、軌道の離散フーリエ変換（DFT）に対応する。

• 結果5：一般的な多部分割（$d \le 5$）
  著者らは、最大5サイトまでの格子のすべての分割に対して、明示的な閉形式の解を提供する。彼らは、最大化因子の関数形式が $n$ の具体的な値に依存し、異なる「レジーム」を定義することに注目している。解かれたすべてのケースにおいて、最大化因子は主要化最小であり、その成分は自然占有数の部分集合の平均である。

• 予想1および2： $d=12$ までの数値的証拠に基づき、著者らは、任意の次元および分割において、最適なブロックスペクトルは常に許容集合内の主要化最小の要素であり、その成分は、インデックスの順序付けられた分割によって定義される自然占有数の部分集合の平均である、という旨の予想を述べている。

• 拡張：

  • 超選択規則（SSR）： 著者らは、局所的な粒子数に関するSSRの影響を分析している。SSRはアクセス可能な相関を減少させるものの、SSRなしで導出された最適ユニタリは効果的であり（厳密な最適解ではない場合でも）、相関の損失は粒子数が大きい極限において消失することを見出した。
  • 純粋ガウス状態： 固定された粒子数を持たない可能性のある一般的な純粋フェルミオンガウス状態について、バランスの取れた $2k$ 個の軌道による $k$ ebits の最大エンタングルメントは、BCS型状態によって達成されることが示されている。

意義と主張
本論文は、多体系物理学において複雑性を低減するための出発点と見なされがちな平均場状態（スレーター行列式および自由状態）を、量子情報処理のための実行可能なリソースへと変貌させることを主張している。

• リソースの特定： 非相互作用フェルミオンは、粒子相関がゼロであるにもかかわらず、最適な回転の下で最大級の軌道もつれを保持できることを示している。
• 操作の実現可能性： 最適な回転は、要素的な2軌道ユニタリ（ビームスプリッター型）で構成されており、これらは光格子、量子ドット、または超伝導プロセッサにおいて実験的に実現可能である。
• 理論的完結： 本研究は、軌道相関の最小化として確立されているものに対し、その最大化の側面を提供することで、同一粒子における軌道の相対性の全体像を完成させている。
• フェルミオン特異性： 自由状態における軌道もつれが、真にフェルミオン的なリソースであることを特定している（自由なボゾン系は常に分離可能であるため）。

著者らは、任意の多部分割に対する一般的な解については、控えめな立場を維持しており、$d \le 5$ および $d \le 12$（数値的）の事例を、すべての $d$ に対する証明された定理としてではなく、自身の予想を裏付ける強力な証拠として提示している。また、導出された主要化関係は、行列解析における独立した関心を持つ新しい結果であることも明記している。