An exact semidefinite characterization of the stoquasticity cone of permutationally invariant Bell operators

技術要約：置換不変なベル演算子のストークスティシティ・コーンに関する厳密な半正定値記述

問題設定
大規模な原子アンサンブル（数百から数十万個の粒子に及ぶ）におけるベル相関は、置換不変（PI）なベル演算子を用いて検出されることが多い。この領域における重要な観察結果は、多くの実験的に実現された演算子において、最大違反状態が非負の振幅を持つ基底で表現できることである。この基底において、ベル演算子は「ストークスティック（stoquastic）」、すなわち、その非対角要素が実数かつ非正であるという性質を持つ。ベル係数がPI演算子をストークスティックにする集合は、**ストークスティシティ・コーン（stoquasticity cone）**と呼ばれる凸領域を形成する。

$K \geq 3$ の演算子におけるこのコーンの特性を明らかにすることは、未解決の問題であった。2体演算子（$K=2$）の場合、コーンはシステムサイズ $n$ に依存しない固定された数の超平面によって定義される多面体錐であることが既知である。しかし、高次の演算子は制約の増殖を伴う。具体的には、ストークスティシティの条件は、すべての非対角要素が非正であることを要求する。このとき、非対角要素の数はシステムサイズ $n$ とともにスケールするため、コーンを定義する線形不等式の数は $n$ とともに増加し、コーンのコンパクトな閉形式による記述を妨げる。この増殖は、量子・古典ギャップを最大化するためのベル演算子の最適化（そのような最適化にはコーンの双対（頂点）記述が必要となる）を阻害している。

手法
著者らは、対称部分空間（ハミルトニアンの最大のブロック）におけるPIベル演算子の構造的特性を利用している。彼らのアプローチは、以下の3つの段階で進行する：

1. 非対角要素の多項式構造： 著者らは、任意の特定の非対角成分 $d$（主対角からの距離 $d$）に沿って、行列要素 $(B)_{k, k+d}$ は任意ではなく、スケーリングされた行インデックス $k$ の単変数多項式の値になることを確立している。具体的には、$(B)_{k, k+d} = \sqrt{f_d(n,k)} h_d(n, k; \vec{\alpha})$ であり、ここで $h_d$ は $k$ に関する次数高々 $K-d$ の多項式であり、係数はベル係数 $\vec{\alpha}$ に対して線形である。
2. 熱力学的極限への還元： 各非対角成分に対して $O(n)$ 個の離散的な不等式をチェックする代わりに、著者らは熱力学的極限（$n \to \infty$）へと移行する。$k = (n-d)x$（ただし $x \in [0, 1]$）と置換することで、離散的な制約の集合は、多項式 $h_d(x, \vec{\alpha})$ が連続区間 $[0, 1]$ 上のすべての $x$ に対して非正であるという連続的な条件へと置き換えられる。これにより、$O(n)$ 個の制約が、単一の連続的な条件 $\sup_{x \in [0,1]} h_d(x, \vec{\alpha}) \leq 0$ へと削減される。
3. 半正定値記述： 区間上の単変数多項式が非正であることを保証する問題は、マルコフ・ルカシュ（Markov–Lukács）の定理を用いて解決される。この定理は、ある単変数多項式が $[0, 1]$ 上で非負であるための必要十分条件は、正定値グラム行列を伴う特定の平方和（SOS）分解を持つことであると述べている。したがって、$\sup h_d \leq 0$ という条件は、線形行列不等式（LMI）の実行可能性と同値である。これにより、ストークスティシティ・コーンは射影スペクトラヘドロン（projected spectrahedron）（またはスペクトラヘドラル・シャドウ）として特徴付けられる。

主要な貢献および結果

• 厳密な半正定値定式化： 本論文は、任意の次数 $K$ のPIベル演算子に対する、熱力学的極限におけるストークスティシティ・コーンの厳密な半正定値記述を提供している。多変数多項式の非負性問題はしばしば階層的な近似（Positivstellensatz）を必要とするが、非対角要素の単変数的な性質により、本手法は単一の半正定値プログラム（SDP）と近似階層なしでの厳密な表現を可能にしている。
• 2体ケース（$K=2$）： 本手法は、2体ストークスティシティ・コーンが3つの超平面によって定義される多面体錐であるという既知の結果を再現する。これは、本アプローチの妥当性を検証するものである。
• 3体ケース（$K=3$）： 3体演算子について、著者らは新たな幾何学的構造を導出している。第2および第3の非対角成分は線形制約（超平面）を与える一方で、第1の非対角成分は二次制約を与える。内部領域におけるコーンの境界は、単一の補助変数を持つ標準的なLMIによって記述される二次錐（quadric cone）（曲面の境界）である。
  • 著者らは、厳密な熱力学的極限において、この曲面の境界が対の超平面（多面体形状）へと退化することを示す。
  • しかし、有限の $n$ においては、曲率は $O(1/n)$ の効果として残る。論文ではこの誤差を定量化しており、熱力学的極限による近似が非常に正確であり（例：$n \approx 5 \times 10^5$）、相対的な偏差は $1/n$ でスケールすることを示している。
• パリティと対称性： 著者らは、極限多項式におけるパリティ特性 $h^\infty_d(1-x) = (-1)^{K-d} h^\infty_d(x)$ を特定した。この対称性は、$K-d$ が奇数の場合、非対角成分の制約が厳密な極限において等式制約へと退化することを意味し、$K-d$ が偶数の場合は、区間上の非負多項式を定義することを意味する。

意義および主張
本論文は、多くの体数のPIベル演算子における不等式の「増殖」という問題に対し、その複雑さは表層的なものに過ぎないことを示すことで、この問題を解決したと主張している。鍵となる洞察は、各非対角成分に沿ったストークスティシティ条件が単一の単変数多項式によって支配されているため、全実行可能領域を単一のSDPで捉えることができるという点にある。

• 統一： 本アプローチは、既知の2体演算子の多面体記述と、これまで未知であった高次演算子の構造を統一する。
• 計算の実行可能性： ベル演算子がストークスティックであるかどうかの判定は、単一の半正定値実行可能性問題へと簡略化される。所望の特性を持つストークスティック・ベル演算子を設計することは、このコーン上の半正定値最適化問題となる。
• 物理的関連性： 本研究の結果は、大規模な原子アンサンブル（ボース＝アインシュタイン凝縮体など）を用いた現在の実験に直接適用可能である。論文では、有限系に対する熱力学的極限記述の忠実度を定量化し、実験データの解析におけるその有用性を確認している。
• 限界と今後の課題： 著者らは、SDPによる特徴付けは厳密であるものの、$K \geq 3$ の場合のコーンの明示的な頂点記述（極限レイ）は依然として未解決の問題であることを指摘している。さらに、qutrit（スピン1）系への拡張は、対称セクターが2次元になるためより困難な問題であると特定している。これは、多変数多項式となり、階層なしでは厳密なSOS表現が一般には不可能になるためである。

要約すると、本論文は、制約の数がシステムサイズに対して線形にスケールする問題を、固定サイズの半正定値プログラムへと変換することにより、置換不変なベル演算子のストークスティシティを分析するための、厳密かつ効率的なフレームワークを提供している。