非局所性、可積分性、およびベル演算子のスペクトルにおける量子カオス

背景と学術的系譜

起源と学術的系譜

本稿で取り上げる問題は、ベルの不等式[1]によって初めて有名に記述された非局所性という現象を理解しようとする長年の探求に端を発する。非局所性は、古典物理学を超えた高度な情報処理タスクに不可欠である[2]が、単純なスピン1/2粒子よりも複雑な系において、これらの非局所相関がどのように現れるかについての我々の理解は、依然として非常に限られている。具体的には、この分野は「高スピン」系における非局所性の特徴付けに苦慮してきた。高スピン系では、各量子サブシステムは、電子のようなスピン1/2粒子の2つの結果とは対照的に、3つ以上の可能な結果を持つことができる。

過去の研究における重要な「痛点」は、ベルシナリオの指数関数的なスケーリングであった。相互作用する量子粒子、測定設定、および可能な結果の数が増加するにつれて、ベル不等式の解析の複雑さは法外に増大する[6]。スピン1/2系では、研究者たちは対称性を利用し、より単純な相関に焦点を当てることでこの問題を回避する方法を見出した[7-11]。これにより、膨大な数の粒子を持つ系を研究することが可能になった[12, 13]。しかし、これらの方法は高スピン粒子に対しては「はるかに開発が進んでいない」。これは、そのような系におけるベル相関の理論的理解と実験的実証の両方に大きなギャップがあることを意味する。本稿は、これらのより豊かな3準位系における非局所性が量子カオスのような他の複雑な量子現象とどのように関連するかについての深い探求を妨げてきた、これらの高次元系（qutrit）のための局所射影測定の滑らかなパラメータ化という基本的なタスクでさえ、より複雑であることを強調している。このような適切な方法論の欠如は、これらのより豊かな3準位系における非局所性が量子カオスのような他の複雑な量子現象とどのように関連するかについての深い探求を妨げてきた。著者らは、特にSU(3)モデルを用いて、このギャップを埋めるために本稿を執筆する必要性を感じた。SU(3)モデルは、動的な複雑さと非局所相関の間の相互作用を探求するのに自然に適している。

直感的な領域用語

• ベル非局所性: 遠く離れたアリスとボブがそれぞれ特別なコインを持っていると想像してほしい。もし二人が同時にコインを投げ、コミュニケーションや事前の取り決めが不可能であったにもかかわらず、常に反対の結果（一方が表、他方が裏）を得たとしたら、それはベル非局所性に似ている。これは、局所的な影響を超えた「不気味な」つながりを示唆する、古典的な説明では説明できないほど強い相関である。
• ベル演算子: これはベル非局所性ゲームの数学的な「スコアカード」と考えてほしい。アリスとボブが行う任意の測定セットに対して、この演算子は特定の値を計算する。この値がある閾値を下回った場合、その量子系がベル非局所性を示していることを意味する。本稿では、この演算子を系のエネルギーとその進化を記述する数学的表現である「有効ハミルトニアン」として扱っている。
• キュートリット: 標準的なライトスイッチにはオンとオフの2つの状態がある。キュートリットは、オフ、薄暗い、明るいの3つの異なる状態を持つ特別なライトスイッチのようなものである。これは量子情報の基本的な単位であり、キュービット（2つの状態）に似ているが、自由度が1つ増えているため、それが記述する系はより複雑で興味深いものになる。
• 量子カオス: ビリヤード台を想像してほしい。もしそれが完璧な長方形であれば、ボールの軌道は予測可能（可積分）かもしれない。しかし、台が非常に不規則で複雑な形状をしている場合、ボールは信じられないほど予測不可能で「カオス的」な方法で跳ね返るだろう。量子カオスとは、この種の複雑で予測不可能な振る舞いが量子系のエネルギー準位にどのように現れるかについてのことであり、しばしばそれらをランダムな数字の寄せ集めのように見せる。
• 可積分性: ビリヤード台の例えに倣って、可積分性とは、台が完璧な形状をしており、ボールの動きが非常に予測可能で、単純で規則的なパターンに従う場合を指す。量子系では、これはエネルギー準位が無相関で、単純な分布（ポアソン）からランダムに抽出された数値のように、予測可能な方法で広がることを意味する。

記法表

記法	説明

問題定義と制約

中核問題の定式化とジレンマ

本稿が取り組む中心的な問題は、高次元多体系における量子非局所性の理解における大きなギャップに起因する。ベルの不等式はスピン1/2（キュービット）系における非局所性の特徴付けに役立ってきたが、サブシステムあたり3つ以上の結果を持つ系（キュートリットまたは高スピン粒子）への適用は、依然としてほとんど開発されていない。

現状は、これらの高次元系において非局所相関がどのように現れるか、そしてより重要なことに、それらが量子カオスのような複雑な物理現象とどのように関連するかについての理解が限定的であることである。過去の研究は、主に対称性を利用することによってキュービットにおけるベルシナリオの指数関数的なスケーリングを軽減することに焦点を当ててきた。しかし、キュートリットに対する同様の方法論ははるかに成熟しておらず、そのような系におけるベル相関の実験的実証は事実上存在しない。さらに、ハミルトニアンのスペクトル特性は量子カオスの確立された診断法であるが、この「スペクトルレンズ」はベル演算子にめったに適用されていない。

望ましい終着点は、多体スピン1系におけるベル非局所性と量子カオスの間の明確なつながりを確立することである。これには、適切な置換対称なベル不等式の導入、対応するベル演算子の構築、そして様々な測定構成下でのこの演算子のスペクトル特性の分析が含まれる。最終的な目標は、ベル不等式の違反を最大化する特定の測定設定を特定し、これらの最適な構成に関連する量子カオス（または可積分性）の性質を理解することである。

本稿が埋めようとしている欠けているリンクは、高次元量子系における最大ベル非局所性と関連するベル演算子のスペクトル統計との間の正確な数学的および物理的関係である。具体的には、本稿は最大非局所性がカオス的（ウィグナー・ダイソン）または可積分的（ポアソン）なスペクトル挙動に関連しているかどうか、そしてどのような根本的なメカニズムがこのつながりを駆動しているかを理解しようとしている。

歴史的に研究者を閉じ込めてきた痛みを伴うトレードオフまたはジレンマであり、本稿が強調しているのは、これらの系で最大ベル不等式違反をもたらす条件が、期待される複雑なカオス的ダイナミクスをもたらさないという直感に反する発見である。代わりに、本稿は、非局所性を最大化する最適な測定設定が、驚くべきことにベル演算子がポアソン的レベル統計（可積分性の兆候）を示すことを明らかにする。逆に、一般的なまたはわずかに摂動された測定は、量子カオスの特徴であるウィグナー・ダイソン統計をもたらす。これはジレンマを提示する：最も深遠な量子非局所性を明らかにする条件そのものが、系のスペクトルダイナミクスを単純化するように見え、この可積分性の微調整された脆い性質を示唆している。

制約と失敗モード

多体高スピン系における非局所性と量子カオスの問題を理解することは、いくつかの厳しく現実的な制約に満ちている。

• 計算スケーリング: 最も重要な壁は、「パーティー、測定設定、および結果の数に対するベルシナリオの指数関数的なスケーリング」である。$n$個のパーティー、それぞれが2つの測定設定と3つの結果（キュートリット）を持つ場合、複雑さはすぐに扱えなくなる。本稿は、各サイト$i \in [n]$に対して異なる測定パラメータを持つベル演算子の構築が「計算上要求が厳しすぎる」と明示的に述べている。この計算上の限界により、考慮された置換対称なベル不等式（PIBI）の分析は$n=32$パーティーに制限された。
• 測定パラメータ化の困難さ: キュートリットの局所射影測定を滑らかにパラメータ化することは、キュービットの場合よりもはるかに複雑である。主な困難は、ベルシナリオに必要なスペクトル特性を維持する連続的なユニタリ演算子の族を構築することにある。パウリ行列の直接的な一般化は、線形結合に対してユニタリ性を維持しないため、測定設定の最適化は自明なタスクではない。
• データの希少性と実験的検証: 高スピン粒子については、「そのような系におけるベル相関の実験的実証はまだ報告されていない」。このような系に対する実際のデータや実験プラットフォームの欠如は、理論モデルと計算シミュレーションが探求の主な手段であることを意味し、直接的な経験的検証を制限する。
• 可積分性の脆さ: 最大ベル違反で観察される可積分挙動は、極めて脆い。これは、最適な測定設定の周りの消滅的に小さい近傍内でのみ現れる。最適な測定設定からの「わずかな逸脱でさえ、ポアソンからウィグナー・ダイソン統計への遷移を誘発する」。さらに、「ポアソン様RCS分布をもたらす測定設定の体積は、$n$の増加とともに減少する」ことは、系サイズが増加するにつれて、この可積分領域を見つけ維持することがますます困難で微調整されることを意味する。これは、観察された可積分性がベル演算子の一般的な特性ではなく、最適な構成の特異な特徴であることを示唆しており、その検出と研究を困難にしている。
• 有限サイズ効果: 最適測定が非ゼロのフィットされたRCSパラメータ（$\lambda > 0$）をもたらす一般則（図2のオレンジ色の点）からの明白な例外は、「対応する既約表現のヒルベルト空間次元の制限に起因する有限サイズ効果」として解釈される。これは、系サイズが小さい場合、スペクトル統計が漸近的な挙動を完全に反映しない可能性があり、可積分性またはカオスの誤解につながる可能性があることを意味する。

なぜこのアプローチなのか

選択の必然性

多体3準位系（キュートリット）のための新しい置換対称なベル不等式（PIBI）を導入し、その関連ベル演算子のスペクトル特性を分析することを中心とした選択されたアプローチは、問題の範囲を考えると、単なる好みではなく必然であった。従来の「SOTA」手法、例えば標準的な畳み込みニューラルネットワーク（CNN）、拡散モデル、またはトランスフォーマーは、この領域とは完全に直交している。それらは、画像生成、自然言語処理、またはデータ分類のようなタスクのために設計された機械学習パラダイムであり、演算子スペクトルにおける量子非局所性とカオスに関する基本的な理論的研究のためではない。したがって、これらの手法は考慮されず、ここでは適用できない。

著者らは、ギャップの大きさ、「しかし、スピン1/2粒子を超える系で非局所相関がどのように現れるか、そしてそれらが複雑な物理的挙動とどのように関連するかについての我々の理解は限られている」と強調されているように、既存の方法の不十分さを当初から認識していた（2ページ）。具体的には、高スピン粒子（キュートリットなど）については、「同様の方法論は…はるかに開発が進んでおらず、そのような系におけるベル相関の実験的実証はまだ報告されていない」（2ページ）。これは、キュービット（スピン1/2）のための確立されたフレームワークが、高次元ヒルベルト空間とSU(3)対称性の固有の複雑さのために、キュートリット（スピン1）に容易に拡張できなかったことを示している。ベルシナリオのパーティー、測定設定、および結果に対する指数関数的なスケーリング[6]は、この複雑さを管理しながら物理的な洞察を維持する、調整されたアプローチの必要性をさらに強調した。キュートリットのための特定のベル演算子の構築は、スペクトル分析を通じてベル非局所性と量子カオスとの間のギャップを埋めるための唯一の方法として解釈された。

比較優位性

この方法の定性的な優位性は、ベンチマークで既存のアルゴリズムを上回るのではなく、以前は解決不可能であった問題空間に取り組む能力と、新しい基本的なつながりを明らかにする能力に由来する。

1. 高次元系への拡張: このアプローチは、多体スピン1系（キュートリット）におけるベル非局所性と量子カオスを調査することを独自に可能にする。指摘されているように、そのような高スピン粒子に対する方法論は「はるかに開発が進んでいなかった」（2ページ）。SU(3)系は、「通常、外部駆動を必要とするSU(2)モデルとは対照的に、固有の量子カオスダイナミクスを既にホストしている」（3ページ）ため、この探求の自然なプラットフォームとして特に魅力的である。
2. 対称性の利用による計算上の実現可能性: ベルシナリオにおける中心的な課題は、その「指数関数的なスケーリング」（2ページ）である。この方法論は、PIBIに固有の置換対称性とシューア・ワイル双対性を利用して、ベル演算子を「多項式サイズの各ブロックにブロック対角化する」（4ページ、セクションB）ことで、これを克服する。これは、$O(N^2)$から$O(N)$への削減として明示的に述べられていないが、このブロック対角化は、各ブロック内での計算複雑性を指数関数的なものから多項式的なものに劇的に減らし、最大$n=32$パーティーの系での分析を可能にする（5ページ）。
3. 量子カオスに対する演算子ベースの視点: 以前の研究は通常、固定ハミルトニアン下での量子状態の進化を通じて量子カオスを調査していたのに対し、この方法論は「根本的に異なる」（7ページ、考察）演算子ベースの視点を提供する。これは、ベル演算子自体のスペクトル統計が異なる測定構成によって変化し、可積分性（ポアソン統計）が「最大ベル違反と一致する微調整された、対称性強化された点」（7ページ、考察）から現れることを明らかにする。この構造的な利点は、非局所性とカオスとの間に、より深く、より本質的なつながりを提供する。

制約との整合性

選択された方法論は、問題の暗黙的な制約と完全に整合し、「問題の厳しい要件と解決策のユニークな特性との結婚」を形成する。

1. 多体スピン1系への焦点: 方法論の核心は、「多体スピン1系に合わせた置換対称な多体ベル不等式」（3ページ）の導入である。これは、この文脈では以前はほとんど探求されていなかったキュートリット系の研究という制約に直接対処する。
2. 指数関数的スケーリングへの対処: 置換対称性とシューア・ワイル双対性を利用してベル演算子をブロック対角化する（4ページ、セクションB）方法論への依存は、多体ベルシナリオにおける指数関数的スケーリングの計算上の課題に直接対処する。これにより、問題はかなりの数のパーティーに対して扱えるようになる。
3. 量子カオスと非局所性の調査: 比率の連続レベル間隔（RCS）を使用したベル演算子のスペクトル特性の分析（4ページ、セクションC）は、可積分性（ポアソン統計）とカオス（ウィグナー・ダイソン統計）を診断するための直接的なメカニズムである。これは、非局所性と量子カオスの相互作用を探求するという目標を完全に満たす。
4. 効率的な測定パラメータ化: 本稿は、「ベルシナリオの構造を尊重しながら、局所射影測定を効率的にパラメータ化する方法」の必要性を明示的に述べている（4ページ、セクションB）。採用されたユニタリベースのパラメータ化（式7）は、ユニタリ性を保証し、ベルシナリオ構造を維持し、効率的な最適化を可能にする。これは、最適な測定設定を見つけるために不可欠である。
5. 置換対称性の維持: PIBI自体は置換対称であり、最適化戦略は「すべてのパーティーが同じ測定ペアを共有する」（$n$個のパーティーすべてに対して$\theta_i^x = \theta^x$）と仮定している（4ページ、セクションB）。これは、最適な測定値の探索を単純化するために、この対称性を完全に活用している。

代替案の却下

著者らは、量子カオスと相関に関する先行研究と比較して、その演算子ベースのフレームワークのユニークな利点を強調することによって、代替アプローチを暗黙のうちに却下している。著者らは次のように述べている。「我々のアプローチは、カオスと量子相関を結びつける以前の努力とは根本的に異なる。以前の研究は通常、エンタングルメント成長[43]、量子ディスコード[44]、またはカオスダイナミクスにおけるレゲット・ガーグ不等式の違反[45]のような量子状態の特性を通じてカオスを調査するが、これらはすべて固定ハミルトニアン下での状態進化に依存している。」（7ページ、考察）。

これらの状態進化ベースの代替案を却下する理由は、測定設定の関数としてベル演算子自体の固有のスペクトル特性を直接調査することを可能にしないからである。これらの先行研究はハミルトニアンを固定し状態の進化を観察するのに対し、本稿の方法論はベル演算子を、測定構成によってスペクトル統計が変化する動的な実体として扱う。これにより、「可積分性が最大ベル違反と一致する微調整された、対称性強化された点から現れる」（7ページ、考察）という発見が可能になる。これは、状態進化中心の分析ではアクセスできない発見である。さらに、前述のように、スピン1/2系のための既存の方法論は、高スピン粒子に対して「はるかに開発が進んでいなかった」（2ページ）ため、この研究の中心であるキュートリット系には適していなかった。

FIG. 3. Histogram of RCS parameters λ resulting from fitting the RCS distribution of the Bell operator constructed from the PIBI (1) with random projectors. Here n = 25 and (p, q) = (25, 0), i.e. the lowest point in Fig. 2, but other irreps show a similar behaviour de- spite having a lower fraction of points exhibiting nonlocality

数学的および論理的メカニズム

マスター方程式

本稿の数学的エンジンの絶対的な核心は、ベル演算子$B$である。これはエルミート演算子であり、その期待値は特定の置換対称ベル不等式（PIBI）に対応する。本稿では、この演算子を式(3)で次のように定義している。

$$ B = \sum_{i \in [n]} \sum_{a \in \{0,1\}} E_{a|x=0}^{(i)} + \sum_{i \neq j \in [n]} \sum_{a \in \{0,1\}} (E_{a|x=0}^{(i)} E_{a|x=1}^{(j)}) - 2 \sum_{i \neq j \in [n]} (E_{0|x=0}^{(i)} E_{1|x=1}^{(j)} + E_{0|x=1}^{(i)} E_{1|x=0}^{(j)}) $$

この演算子は、その最小固有値がベル不等式の最大量子違反を直接定量化するため、中心的なものである。この演算子のスペクトル特性、特にそのエネルギー準位の統計が、量子カオスと可積分性を特徴付けるために分析される。

項ごとの解剖

マスター方程式（式3）を項目ごとに分解してみよう。

• $B$: これはベル演算子そのものである。

  • 数学的定義: $n$個のキュートリットの$3^n$次元ヒルベルト空間上で作用するエルミート演算子。
  • 物理的/論理的役割: これは有効な多体ハミルトニアンとして機能し、その期待値$\text{Tr}[\rho B]$は、与えられた量子状態$\rho$と測定設定に対するPIBI（式1）の違反を定量化する。負の期待値は非局所性を示す。
  • この構造である理由: 演算子は、PIBI（式1）の集団確率の構造を反映して、1体および2体測定項の和として構築されている。

• $\sum_{i \in [n]}$: これは$n$個のパーティー（サブシステム）すべてにわたる和である。

  • 数学的定義: $i=1$から$n$までの和。
  • 物理的/論理的役割: これはベル不等式と演算子の置換対称性を反映している。すべてのパーティーは対称的に扱われる。
  • 和である理由: 個々のサブシステムまたはサブシステムのペアからの寄与をシステム全体にわたって集計するため。

• $\sum_{a \in \{0,1\}}$: これは測定結果$a=0$と$a=1$にわたる和である。

  • 数学的定義: 指定された2つの結果にわたる和。
  • 物理的/論理的役割: これらはPIBI（式1）で考慮される特定の可能な結果である。本稿では、結果$a=2$は、非信号制約$E_{2|x} = I - E_{0|x} - E_{1|x}$を通じて暗黙的に処理されると述べている。
  • 和である理由: 異なる測定結果に関連する確率または演算子を組み合わせるため。

• $E_{a|x}^{(i)}$: これはパーティー$i$に対する局所的な正演算子値測定（POVM）演算子を表す。

  • 数学的定義: $E_{a|x}^{(i)} = \mathbb{I}^{(i-1)} \otimes E_{a|x} \otimes \mathbb{I}^{(n-i)}$。ここで、$E_{a|x}$は単一キュートリットに対する局所POVMであり、$\mathbb{I}^{(k)}$は$k$個のサブシステム上の恒等演算子である。局所$E_{a|x}$は、ユニタリ演算子$U(\theta_x)$から逆フーリエ変換のような表現を用いて導出される。本稿の例（式9および$P_{00}$に関する記述）に基づくと、結果$a=0$の場合、$E_{0|x} = \frac{1}{3} (U(\theta_x)^3 + U(\theta_x)^2 + U(\theta_x))$、結果$a=1$の場合、$E_{1|x} = \frac{1}{3} (U(\theta_x)^3 + \zeta U(\theta_x)^2 + \zeta^2 U(\theta_x))$となる。
  • 物理的/論理的役割: これらは、サブシステム$i$で測定設定$x$を実行し結果$a$を得ることに対応する基本的な量子演算子である。これらは、ベル不等式における古典的な確率を量子力学的な観測量に変換する。
  • テンソル積である理由: 他のサブシステムをそのままにしておく一方で、特定のサブシステム$i$に作用する局所測定を記述するため。

• $E_{a|x=0}^{(i)}$: これは測定設定$x=0$と結果$a$を持つパーティー$i$の局所POVMである。

  • 物理的/論理的役割: 局所測定の特定の選択を表す。

• $E_{a|x=1}^{(j)}$: これは測定設定$x=1$と結果$a$を持つパーティー$j$の局所POVMである。

  • 物理的/論理的役割: 局所測定のもう1つの特定の選択を表す。

• $E_{a|x=0}^{(i)} E_{a|x=1}^{(j)}$: これは異なるパーティー$i$と$j$の2つの局所POVMの積である。

  • 数学的定義: これは局所POVMのテンソル積であり、実質的に$( \mathbb{I}^{(i-1)} \otimes E_{a|x=0} \otimes \mathbb{I}^{(j-i-1)} \otimes E_{a|x=1} \otimes \mathbb{I}^{(n-j)} )$である。
  • 物理的/論理的役割: パーティー$i$が設定$x=0$を、パーティー$j$が設定$x=1$を使用する同時測定を表し、両方とも結果$a$を得る。これらの項は2体相関に対応する。
  • 乗算（テンソル積）である理由: 別々の、通信しないサブシステム上の同時イベントを記述するため。

• $\sum_{i \neq j \in [n]}$: これはパーティーのすべての異なるペア$(i, j)$にわたる和である。

  • 数学的定義: $i, j \in \{1, ..., n\}$で$i \neq j$となる和。
  • 物理的/論理的役割: システム内のすべての可能な2体相関を考慮し、再び置換対称性を反映する。
  • 和である理由: すべての可能なサブシステムのペアからの寄与を集計するため。

• 係数$-2$: これは重み係数である。

  • 数学的定義: スカラー乗数。
  • 物理的/論理的役割: これは、不等式を定義するために特定の相関項がこの特定の係数で減算されるPIBI（式1）の構造から直接由来する。
  • この値である理由: これは、古典的な境界を設定するために設計された、この特定のベル不等式の特徴的な係数である。

• $\zeta = e^{-2\pi i/3}$: これは3の3乗根である。

  • 数学的定義: 局所POVMの定義で使用される複素数。
  • 物理的/論理的役割: 3準位系（キュートリット）の射影測定を定義するために使用される逆フーリエ変換から生じる。異なる測定結果を定義するために不可欠である。

• $U(\theta_x)$: これは測定設定をパラメータ化するユニタリ演算子である。

  • 数学的定義: $U(\theta) := e^{i g(\theta)} D e^{-i g(\theta)}$。ここで、$g(\theta) = g_0 + \sum_{l=1}^M \theta_l (g_l - g_0)$はエルミート行列$\{g_k\}$の基底から構築されたエルミート演算子であり、$D = \text{diag}(1, \zeta, \zeta^2)$は対角行列である。$\theta = (\theta_1, ..., \theta_M)$は実パラメータのベクトルである。
  • 物理的/論理的役割: このパラメータ化により、局所測定が射影的かつユニタリであることが保証され、量子力学的な構造が維持されつつ、測定選択の連続的な変化と最適化が可能になる。

正直なところ、PIBI（式1）からベル演算子（式3）への正確なマッピングについては完全には確信が持てない。本稿では式(3)が関連するベル演算子であると述べているが、式(1)の確率を演算子（式2を使用）に直接項ごとに変換すると、式(1)に存在する一部の1体および2体項が欠落した、わずかに異なる演算子構造が得られるだろう。式(3)のこの形式につながる暗黙の単純化または参照された方法論[27, 28]からの特定の文脈がある可能性がある。さらに、テキスト中の$P_{a|x}(\theta_x)$の一般的な式は、$P_{00}$および$P_{01|01}$の特定の例と比較して、$U(\theta_x)$のべき乗と$\zeta$係数にタイプミスがあるように見える。私は、提供された例に基づいた局所POVMの定義を解釈した。

ステップバイステップの流れ

ここでは、単一の抽象的なデータポイント、つまり量子測定の概念が、ベル演算子を構築および分析するための数学的機構を通過する様子を想像してみよう。

1. 測定設定入力: プロセスは、測定パラメータのセット、全体として$\theta = (\theta_0, \theta_1)$で表されるものから始まる。これらのパラメータは、2つの設定、$x=0$と$x=1$に対する局所測定の選択を定義する。

2. エルミート演算子構築（式6）: 各測定設定$\theta_x$に対して、エルミート演算子$g(\theta_x)$が構築される。これは、基底エルミート行列（$g_0, g_1, \dots, g_M$）をパラメータ$\theta_l$によって与えられた重みで混合するようなものである。このステップにより、後続のユニタリ演算子が望ましい特性を持つことが保証される。

3. ユニタリ演算子生成（式7）: エルミート演算子$g(\theta_x)$は、ユニタリ演算子$U(\theta_x)$を生成するために使用される。これには、$i g(\theta_x)$の指数関数を取り、それを対角行列$D$（単位根を含む）で共役し、その後$e^{-i g(\theta_x)}$で逆共役することが含まれる。これにより、$U(\theta_x)$が固定スペクトルを維持することが保証され、射影測定の定義に不可欠である。

4. 局所POVM定義: 各$U(\theta_x)$から、結果$a \in \{0,1,2\}$の局所POVM $E_{a|x}$が導出される。これは、逆フーリエ変換のような操作を使用して行われる。例えば、結果$a=0$の場合、$E_{0|x}$は$U(\theta_x)$のべき乗（具体的には、$U(\theta_x)^3 + U(\theta_x)^2 + U(\theta_x)$）を$1/3$でスケーリングして合計することで形成される。他の結果については、特定の単位根$\zeta$が係数として導入される。

5. グローバルPOVM展開: 各局所POVM $E_{a|x}$は、$n$個のキュートリットのヒルベルト空間全体に作用するように展開される。パーティー$i$の場合、$E_{a|x}^{(i)}$は、$i$番目のキュートリットに$E_{a|x}$が作用し、他の$n-1$個のキュートリットには恒等演算子$\mathbb{I}$が作用することを意味する。これにより、集合ベル演算子の構成要素が作成される。

6. ベル演算子組み立て（式3）: これらのグローバルPOVM $E_{a|x}^{(i)}$は、ベル不等式の構造に従って組み合わせられる。

   • 最初の項は、設定$x=0$と結果$a=0,1$のすべての単一パーティー測定を合計する。
   • 2番目の項は、パーティー$i$が設定$x=0$を、パーティー$j$が設定$x=1$を使用し、両方とも同じ結果$a=0$または$a=1$を得るすべての2パーティー同時測定を合計する。
   • 係数2で減算される3番目の項は、パーティー$i$が設定$x=0$を、パーティー$j$が設定$x=1$を使用するが、特定の結果が異なる（$0,1$または$1,0$）2パーティー同時測定を合計する。
     この全体の組み立てにより、ベル演算子$B(\theta)$が得られる。

7. 固有値計算: $B(\theta)$が構築されると、その固有値が計算される。これらの固有値は、この有効ハミルトニアンの可能な「エネルギー準位」を表す。最小固有値、$\lambda_{min}$は、最大量子違反を示すため、特に興味深い。

8. RCS計算とフィッティング（式5）: 整列された固有値は、連続レベル間隔の比（RCS）を計算するために使用される。この分布$P(r)$は、パラメータ$\lambda$を抽出するために、補間式（式5）にフィットされる。この$\lambda$値は、このメカニズムの部分の最終的な出力であり、システムがポアソン統計（$\lambda=0$、可積分性）またはウィグナー・ダイソン統計（$\lambda=1$、量子カオス）を示すかどうかを示す。

最適化ダイナミクス

メカニズムの学習と収束は、ベル不等式違反を最大化する最適な測定設定$\theta$を見つけることに中心がある。これは、ベル演算子の最小固有値、$\lambda_{min}(B(\theta))$を最小化することに対応する。

本稿では、これを通常、数値シーソーや確率的勾配降下法のような技術を用いて取り組まれる非凸最適化問題として説明している。本質的に、システムは$\theta$の高次元パラメータ空間を反復的に探索する。

1. 損失地形: ここでの「損失地形」は、$\lambda_{min}(B(\theta))$の関数である。目標は、この地形の最も深い「谷」を見つけることである。本稿の主な発見は、この地形の最適な点（最大違反）が非常に特定のスペクトル特性、すなわちポアソン的レベル統計（$\lambda=0$）に対応することである。これは可積分性を示す。

2. 反復更新:

   • 勾配ベースの方法: 確率的勾配降下法を使用する場合、アルゴリズムはパラメータ$\theta$に対する$\lambda_{min}$の勾配を計算する。この勾配は、最も急な上昇（または下降、最小化している場合）の方向を示す。その後、パラメータ$\theta$は、負の勾配方向に移動することによって反復的に更新され、徐々に損失地形の谷に下降していく。
   • 数値シーソー: これは、他のすべてのパラメータを固定したまま、1つのパラメータ$\theta_l$を最適化し、次に次のパラメータに進み、収束するまで繰り返すことを含む。これは、完全な勾配の計算が困難または計算コストが高い非凸問題によく使用される。

3. 収束と地形の形状: 本稿は、可積分挙動が「微調整された」特徴であることを明らかにする。これは、最適な領域が非常に小さく脆いことを意味する。最適な点周りの「ポアソン様挙動の体積」は、$n$が増加するにつれて急速に減少する。これは、可積分な最小値が孤立しており、カオス的（ウィグナー・ダイソン）統計をもたらす広大な領域に囲まれていることを示唆している。これにより、最適設定からのわずかな摂動がすぐに系をカオス領域に押し込むため、最適化は繊細なタスクとなる。

4. 対称性の出現: 最適化ダイナミクスに関する重要な洞察は、最大量子違反の点の近くでベル演算子にパリティ対称性が現れることである。これは、最適化プロセスが成功した場合、実質的にベル演算子$B(\theta_{opt})$が特定のパリティ演算子と可換になるような測定パラメータのセットを見つけることを意味する。この可換性は、パリティ固有値基底における$B(\theta_{opt})$のブロック対角化構造につながる（図5に示す）。このブロック対角化は固有値レベルを分離し、スペクトル全体にわたるレベル反発を抑制する。これは、最適な違反におけるポアソン的統計（可積分性）の観察された出現の構造的な説明である。したがって、最適化メカニズムは、成功した場合、単に数値的な最小値を見つけるだけでなく、パラメータ空間内の点を見つける。そこでは、演算子のスペクトル特性を根本的に変更する、基本的な対称性が現れる。

結果、限界、および結論

実験設計とベースライン

彼らの主張を厳密に検証するために、著者らは多体3準位系、すなわちキュートリットのための置換対称な多体ベル不等式（PIBI）を中心とした実験フレームワークを設計した。中心的なアイデアは、有効ハミルトニアンとして解釈できる関連ベル演算子$B(\theta)$を定義し、その後、様々な測定構成下でのそのスペクトル統計を分析することであった。パラメータ$\theta$は、各パーティーが選択した局所測定設定をカプセル化する。

キュートリットのような高次元系における重要な課題は、ユニタリ性とベルシナリオ構造を維持しながら、局所射影測定を滑らかにパラメータ化することである。著者らは、ユニタリベースのパラメータ化を採用し、パラメータベクトル$\theta_x$から量子射影子$P_{a|x}^{(i)}(\theta_x)$を構築することで、これを解決した。最適化問題を単純化するために、彼らはPIBIの置換対称性を活用し、最適な違反はすべてのパーティーが同じ測定ペア、すなわちすべての$i \in [n]$に対して$\theta_x^{(i)} = \theta_x$を共有する場合に達成されると仮定した。これにより、最適化はグローバルパラメータセット$(\theta_0, \theta_1)$に削減された。

量子カオスを診断するための主要なツールは、スペクトル展開の複雑さを回避する、連続レベル間隔の比（RCS）$P(r, \lambda)$であった。RCS分布のパラメータ$\lambda$は、重要な補間因子として機能する：$\lambda=0$はポアソン統計を示し、可積分挙動を示唆する一方、$\lambda=1$はウィグナー・ダイソン統計（具体的には、ガウス直交アンサンブルまたはGOE）に対応し、量子カオスの特徴である。

彼らの数学的主張がテストされたベースラインモデルには、以下が含まれる。
1. ランダム測定設定: 一般的な測定選択がカオス的挙動をもたらすことを示すために、著者らはSU(3)から行列をサンプリングすることによって$10^3$以上のランダム射影子を生成し、結果のベル演算子のRCS分布を計算した。これは、可積分性が一般的な特性ではないことを示すための制御として機能した。
2. 最適点周りの摂動: 観察された可積分挙動の堅牢性を評価するために、彼らは最適な測定設定を体系的に摂動させ、RCS分布がポアソンからウィグナー・ダイソン統計にどれだけ早く遷移するかを観察した。
3. 異なるSU(3)既約表現（irrep）: 分析は、特定の単一部分空間に固有の発見ではないことを保証するために、$(p,q)$ペアによって特徴付けられる様々な既約表現全体で実行された。

実験は、$n=8$（彼らのPIBIが非局所性を検出する最小の系サイズ）から最大$n=32$パーティーまで実行された。これは彼らの計算上の限界であった。結果は、この範囲全体で一貫して観察され、彼らの発見の妥当性を強化した。

証拠が証明すること

中核メカニズム—最適な量子測定、非局所相関、および可積分性との間の相互作用—が実際に現実で機能したという決定的な、否定できない証拠は、説得力があり多面的である。

第一に、最も顕著な証拠は、RCS分布の直接比較にある。図1に示すように、$n=25$パーティーと$(21,2)$既約表現の場合、最適な測定設定（青い曲線）で得られたベル演算子は、可積分性の明確な兆候であるポアソン分布（$\lambda=0$）に完全に適合する。対照的に、ランダム測定設定（赤い曲線）が使用されると、RCS分布はGOEウィグナー・ダイソン統計（$\lambda=1$）に密接に一致し、紛れもなく量子カオスを示唆する。この視覚的な証拠は、測定設定の選択のみに基づいて可積分挙動とカオス的挙動の遷移を示す強力なデモンストレーションである。

第二に、図2は$n=25$の様々な既約表現全体にわたる包括的な要約を提供する。これは、最大量子違反を置換対称性の尺度である$\eta = p/(p+q)$に対してプロットする。決定的に、最大ベル違反が達成される既約表現を表す青い点は、一貫してポアソンRCS統計（$\lambda=0$）を示す。これは、最大非局所性と可積分性を直接結びつける。一部の「オレンジ色の点」は中間的な$\lambda$値（$0 < \lambda < 1$）を示すが、これらは有限サイズ効果として解釈され、真の反例ではなくクロスオーバー領域を示唆しており、大きな$n$の漸近的な極限でポアソン統計に収束すると予想される。

第三に、堅牢性分析はこれらの主張をさらに確固たるものにする。図3は、$(p,q)=(25,0)$のランダム射影子に対する$\lambda$パラメータのヒストグラムであり、$\lambda$値が1の近くに強く集中していることを示しており、一般的な測定設定が実際にウィグナー・ダイソン統計をもたらすことを確認している。興味深いことに、ランダム設定であっても、高い割合（この場合75.78%）が依然として非局所性を検出したが、カオス的なスペクトル署名を伴っていた。これは、可積分領域のユニークな性質を強調している。さらに、図4は、ポアソン様RCS分布をもたらす最適な点周りの「体積」が$n$の増加とともに急速に減少することを示している。これは、可積分挙動の脆く微調整された性質を実証しており、これは一般的な特性ではなく、最適な設定の消滅的に小さい近傍内でのみ現れる。

最後に、著者らは最大量子違反の点の近くでパリティ対称性の出現を発見した。図5は、$n=15$キュートリットの最適なベル演算子をパリティ固有値基底で示しており、明確なブロック対角化構造を示している。パリティ演算子と可換なベル演算子に起因するこのブロック対角化は、各ブロック内でのポアソン統計の出現を自然に説明する。ブロックの分離は、スペクトル全体にわたるレベル反発を抑制する。最大違反固有状態が一貫して特定の対称セクター（奇数$n$の場合はeeo、偶数$n$の場合はeee）に存在する事実は、これが真の出現対称性であり、数値的アーティファクトではないことを強く示唆しており、これらの構成の特別な性質を強化している。

限界と将来の方向性

これらの発見は、非局所性、可積分性、および量子カオスとの間に深遠なつながりを示しているが、固有の限界を認識し、将来の探求の方向性を考慮することが重要である。

1つの明確な限界は、パーティーの数、$n$に対する計算上の制約である。現在の分析は、ベルシナリオの指数関数的なスケーリングにより、$n=32$に限定されている。これは、特に「クロスオーバー」領域（図2のオレンジ色の点）のような一部の観察結果が有限サイズ効果として解釈されることを意味する。この解釈はもっともらしいが、決定的な確認には、現在計算上不可能であるより大きな$n$への推進が必要となるだろう。将来の研究では、より大きな系のためのより効率的な数値的手法または解析的近似を探索して、漸近的な挙動を確認できる可能性がある。

もう1つの側面は、研究されているベル不等式の特異性である。本稿は、特定の2入力、3出力の置換対称ベル不等式に焦点を当てている。最大違反における可積分性のスペクトル署名が、他の種類のベル不等式またはより複雑な多体ベルシナリオに一般化されるかどうかは未解決の問題である。将来の研究では、他の3出力置換対称ベル不等式[27]を体系的に調査し、分析を異なる入力数または出力数に拡張できる可能性がある。

すべてのパーティーにわたる同一の測定設定（$\theta_x^{(i)} = \theta_x$）の仮定は、最適化問題を大幅に単純化する。しかし、各パーティーが異なるランダム測定を設定できるより一般的なシナリオは、ポアソンRCS分布からのさらなる顕著な逸脱をもたらす可能性が高いと本稿は述べている。このより一般的なケースを探求することは、可積分性が現れる、または破壊される条件についてのより深い洞察を提供する可能性がある。

将来に向けて、これらの発見はいくつかのエキサイティングな議論のトピックと研究の方向性を開く。

1. より深い解析的関係と対称性: 図2で観察された興味深いパターンは、同じハイパーチャージを持つ既約表現が明確な曲線に沿って整列していることを示しており、最大PIBI違反とSU(3)部分セクターの置換対称性の度合いとの間のより深い解析的関係を示唆している。これらの数学的なつながりを解き明かすことは、量子相関とカオスのより根本的な理解につながる可能性がある。
2. 量子多体系の自己テスト: 出現するパリティ対称性と最大量子違反近くの可積分挙動は、量子多体系の自己テスト[46]のための強力なツールとなる可能性がある。これらの構成が、基底となる量子状態と測定構造の単純化された特徴付けを可能にする場合、それは量子デバイスを検証するための新しい実験プロトコルへの道を開く可能性がある。
3. 量子カオスに対する演算子ベースの視点: 本稿は、ベル演算子自体のスペクトル統計が測定構成によって変化する、量子カオスをプローブするための新しい演算子ベースの視点を導入している。これは、固定ハミルトニアンに対する状態進化に依存する従来の、状態進化中心のアプローチとは対照的である。この視点をさらに探求することは、量子カオスの性質とその基本的な量子情報特徴との関連性についての新しい洞察をもたらす可能性がある。このフレームワークは、カオスの他の演算子ベースの診断とどのように関連しており、量子情報における他の種類の演算子に一般化できるか？
4. 可積分性とカオスの工学: 可積分性が脆く微調整された特徴であるため、特定のタスクのために可積分挙動を維持または誘発するように量子系または測定プロトコルを設計できるか？逆に、堅牢で複雑なダイナミクスが望まれるアプリケーションのためにカオスへの遷移を利用できるか？これは、量子制御と量子コンピューティングに影響を与える可能性がある。
5. 高スピン系の役割: スピン1（キュートリット）系への焦点は、超冷却原子プラットフォームで自然に発生し、固有の量子カオスダイナミクスをホストするため、特に適切である。これらの発見がさらに高スピン系（$d > 3$のクディット）にどのようにスケールするかについてのさらなる調査は、高次元に固有の普遍的な原理または新しい現象を明らかにする可能性がある。

FIG. 1. Ratio of Consecutive level Spacings (RCS) for the Bell oper- ators associated with the PIBI (1) for n = 25, obtained with optimal (blue) and random (red) measurement settings. Spectra are shown for the irreducible representation (p, q) = (21, 2), chosen for illus- tration, which has 825 eigenvalues in the symmetric subspace. Solid curves are fits to the interpolating RCS function in Eq. (5), which spans Poisson statistics (λ = 0, indicating integrability) to Wigner- Dyson statistics (λ = 1, indicating chaos). In the Wigner-Dyson case, the Gaussian Orthogonal Ensemble (GOE) seems to provide a better fit than the Gaussian Unitary Ensemble (GUE), suggesting that time- reversal symmetry is preserved in the chaotic regime FIG. 2. Maximal quantum violation of the PIBI (1) for n = 25 parties, restricted to the (p, q) irrep subsector of SU(3). The classical bound is βc = 0, so ⟨B⟩< 0 certifies nonlocality. The parameter η = p/(p+q) quantifies the degree of permutation invariance of each irrep, with η = 1 being the fully symmetric case. Blue points cor- respond to irreps (p, q) whose Bell operator exhibits Poisson RCS statistics, signalled by λ = 0, indicative of integrability. Orange points correspond to irreps for which the fitted RCS parameter is non- zero (0 < λ < 1). The histograms of the orange cases are fitted with values λ = 0.11 for the irrep (15, 2) and λ = 0.456 for (9, 8), sug- gesting a crossover regime between the Poisson and GOE limits, an interpretation further supported by the significant bin weight in their RCS histograms at small spacing values (see Supplementary Figures 1, 2 and Supplementary Tables I and II [28] for explicit values of all λ’s obtained and some illustrative histograms). Irreps shown in gray do not detect nonlocality and are included only for completeness. Dashed lines connect irreps with same p + 2q and p −q. The largest violations occur in the fully symmetric sector, as expected from the permutationally invariant structure of the Bell operator [3, 31] FIG. 5. Optimal Bell operator B for n = 15 qutrits in the parity eigen- basis, revealing a block-diagonal structure with four non-empty par- ity sectors (for odd n these are ooo, eeo, eoe, and oee). Gray boxes enclose the sectors as a visual guide. The color scale indicates ma- trix element magnitudes; zero entries are shown in white to highlight their sparsity