Symmetric multipartite Bell inequalities via Frank-Wolfe algorithms

技術要約：フランク・ウルフ法による対称マルチパーティ・ベル不等式

問題提起
本論文は、マルチパーティ量子状態、特にGHZ（Greenberger-Horne-Zeilinger）状態の非局所性を特徴付ける際の計算上の困難さ（計算不可能性）に対処している。マルチパーティのベル・シナリオでは、パラメータ空間がパーティ数（$N$）および測定設定数（$m$）に対して指数関数的にスケールする。その結果、量子状態が局所的な相関と区別できなくなる臨界可視性（$v$）を決定することや、ベル不等式の対応する局所境界を計算することは、標準的な線形計画法や網羅的な列挙を用いる場合、中程度の$N$および$m$の値であっても不可能となる。著者らは、特定のシナリオ、すなわち$N$個のパーティがブロッホ球の$XY$平面内で正多角形を形成する$m$個の平面測定を行うケースに焦点を当てている。

手法
著者らは、この指数関数的なスケーリングを克服するために、2つの主要なツールを組み合わせたフレームワークを提案している。

1. 対称性の活用： 著者らは、GHZ状態と正多角形測定設定に固有の対称性を特定し、これを活用している。これらの対称性は、シナリオに遍在するものではなく、その特定のインスタンスに埋め込まれたものである。著者らは、パーティの置換、測定インデックスの巡回シフト、および符号反転によって生成される群作用を定義している。レイノルズ作用素（群による平均化演算子）を適用することで、高次元の相関テンソルおよび局所多面体を、大幅に縮小された対称部分空間へと投影する。
   • これにより、探索空間の次元が$m^N$から$\lceil m/2 \rceil$へと減少する。
   • 決定論的戦略の列挙数が、$2^{(N-1)(m-1)}$から、バイナリ・ネックレスの数を含む組合せ項である$\binom{u_m + N - 2}{N-1}$（ここで$u_m$は巡回シフトの下での軌道の数）へと減少する。
2. 対称フランク・ウルフ・アルゴリズム： 違反を最大化するベル不等式（分離超平面）を見つけるために、著者らはフランク・ウルフ・アルゴリズム（条件付き勾配法）を採用している。局所多面体のすべての頂点を列挙する代わりに、アルゴリズムは現在の勾配に対する「ワーストケース」の決定論的戦略を見つけるために、線形最小化オラクル（LMO）を反復的に解く。
   • アルゴリズムは、縮小された対称部分空間内で動作する。
   • 「アクティブ・セット」を利用して、以前に見つかった頂点を再利用し、高コストなLMO呼び出しの回数を削減する。
   • LMOについては、大きな$N$に対してはヒューリスティック（交互最小化）を用い、小さな$N$（具体的には$N \ge 8$）に対しては、決定論的戦略の軌道への探索を制限するために正確な軌道列挙を用いる。

主な貢献

• アルゴリズムの加速： 論文は、対称性の低減がいかにベル不等式の計算と局所境界の計算を劇的に加速させるかを実証している。これにより、著者らは、$\sim 10^{21}$個の戦略を列挙し、$\sim 10^9$個の要素を持つテンソルを扱う必要があるインスタンス（$N=10, m=9$）を解くことが可能となった。
• 新しいベル不等式： この手法は、対称化された局所多面体のファセット（facet）を生成する。$N=3$から$10$において、これらの不等式はGHZ状態の非局所性堅牢性に関する既知の最良の上界を与える。
• $m=4$の場合の一般化： 著者らは、$m=4$の測定に関する不等式のパターンを特定した。これを任意のパーティ数$N$へと一般化し、局所境界の閉形式の表現を導出した。これらの境界はルーカス数列を形成し、ミン（Mermin）の不等式による漸近的なノイズ堅牢性を改善する。
• 検出効率の解析： 著者らは、これらの不等式に対する臨界検出効率の閾値を計算し、$N=4$の場合、$m$が増加するにつれて閾値が0.7706から0.7544へと改善されることを示している。

結果

• 臨界可視性： 著者らは、$N \in [3, 10]$および様々な$m$に対して、臨界可視性 $v_{GHZ_N}^m$ の数値的値を提供している。例えば、$N=10$かつ$m=9$の場合、臨界可視性は約$0.02301$である。
• 非単調性： 奇数$N$において、臨界可視性が$m$に対して非単調であることを示しており、これは非対称な構成において同一の正多角形測定を使用することの劣最適性に起因する。
• 堅牢性の下界： フランク・ウルフ・アルゴリズムを介して明示的な局所モデル（凸分解）を構築することにより、著者らは$XY$平面における射影測定の非局所性閾値（$v_{GHZ_N}^{XY}$）に対する厳密な下界を導出している。$N=10$の場合、約$0.01974$という下界を確立している。
• 非局所性の活性化： 導出された境界を用いて、著者らは$N=10$のスターネットワークにおける非局所性の活性化を実証している。初期の2量子ビット状態が局所的であるネットワークにおいて、有限の数の測定（$m=8$または$9$）によって非局所性が活性化できることを示しており、これは以前は無限の測定を必要とするとされていた結果を洗練させたものである。

意義と主張
本論文は、対称性が高次元の制約付き最適化問題を解くためのフランク・ウルフ・アルゴリズム内でいかに効果的に活用され得るかを示す概念実証として位置づけられている。著者らは、自身の手法を用いることで、線形計画法や従来のフランク・ウルフ法の適用が可能な標準的なシナリオを大幅に超えることができると主張している。

その意義は以下の点にある：

1. 最大10パーティまでのGHZ非局所性の堅牢性に関する、最もタイトな既知の境界を提供していること。
2. 有限の測定設定が、スターネットワークにおける非局所性活性化を証跡するのに十分であることを示し、無限の測定を必要としていた従来の知見を精緻化したこと。
3. 指数関数的な複雑さを、特定の対称的なインスタンスに対して管理可能なレベルへと減少させる計算フレームワークを提供したこと。

著者らは、$m=4$の不等式のパターンは見出したが、任意の$m$に対する一般的なケースは依然として未解決の課題であるとして、自身の不等式のパターンに関する一般化については控えめな姿勢をとっている。また、不等式の整数係数が非常に大きくなる可能性があることも指摘しており、これは既知の代数的性質であるが、文献によく見られる小さな整数とは対照的である。本研究は、グロテックスキー定数などの定数の境界をさらに改善できる可能性のある、他の対称的なケース（二部テンソルなど）を特定するためのステップとして提示されている。