ENTCALC: Toolkit for calculating geometric entanglement in multipartite quantum systems

「ENTCALC: 多粒子量子系における幾何学的エンタングルメント計算のためのツールキット」テクニカル・サマリー

問題提起
多粒子系における量子エンタングルメントの定量化は、計算負荷の高い課題である。エンタングルメント尺度（エンタングルメント・エントロピーやネガティビティなど）はいくつか存在するが、それらにはしばしば限界がある。ある尺度は純粋状態にのみ適用可能であり、別の尺度は特定のエンタングル状態（PPTエンタングル状態など）を検出できない。また、可分性の判定はNP困難な問題である。具体的には、一般的な混合多粒子状態に対する幾何学的エンタングルメントの計算は、可分状態の集合に対する最適化を伴うため困難である。既存のアプローチは、純粋状態に限定したり、ペアワイズ（二体）の部分系のみに焦点を当てたり（全域的な多粒子構造を無視）、あるいは定量的な尺度を提供しないエンタングルメント・ウィットネスに依存したりすることが多い。したがって、解析解が得られない場合でも、純粋状態および混合多粒子状態の両方に対して幾何学的エンタングルメントの高精度な推定値を提供できるツールが必要とされている。

手法
著者らは、Python (entcalcpy) および MATLAB (entcalc) の両方で実装されたソフトウェアパッケージ ENTCALC を提示する。このツールキットは、下界と上界の両方を計算することで幾材的エンタングルメント ($E_G$) を推定し、真の値が存在する区間を特定する。

1. 幾何学的エンタングルメントの定義:

   • 純粋状態 $|\psi\rangle$ に対して、$E_G$ は $1 - \max_{|\phi\rangle \in S} |\langle\phi|\psi\rangle|^2$ と定義される。ここで $S$ は完全可分状態の集合である。
   • 混合状態 $\rho$ に対して、それは凸屋根拡張（convex-roof extension）を介して $E_G(\rho) = \min_{\{p_i, |\psi_i\rangle\}} \sum p_i E_G(|\psi_i\rangle)$ と定義される。
   • 量子フィデリティを用いた等価な定式化は、$E_G(\rho) = 1 - \max_{\sigma \in S} F(\rho, \sigma)$ である。

2. 純粋状態のためのアルゴリズム:

   • SDPアプローチ: 本パッケージは、PPT（部分転置正値）緩和に基づく半定値計画法（SDP）の定式化を使用する。三体状態の場合、下界はPPT状態とのオーバーラップを最大化することから導かれる。同時に、誤差推定（$4\sqrt{\epsilon}$、ここで $\epsilon$ は最適なPPT状態の最大固有値に関連する）を伴う上界も計算される。
   • 勾配降下法: 最も近い可分状態を見つけるために反復アルゴリズムが使用され、これが上界を提供する。この手法は混合状態のサブ関数としても用いられる。

3. 混合状態のためのアルゴリズム:

   • 上界: 対象の状態の可分分解を反復的に洗練させ、フィデリティを最大化する勾配降下アルゴリズムを用いて計算される。
   • 下界: ユーザーが精度と計算コストのバランスを選択できるよう、4つの異なる下界手法を実装している。
     • 下界 1 (SDP): PPT状態の集合上でフィデリティを最大化することに基づく。効率的であるが、PPTエンタングル状態に対してはゼロとなる。
     • 下界 2 (k-対称拡張): $k$-対称拡張を取り入れることで精度を向上させる。これは、$k \to \infty$ において可分状態の集合へと収束する、よりタイトな緩和の階層を形成する。
     • 下界 3 (純粋度-エンタングルメント相補性): 熱力学とエンタングルメントの関係から導出され、混合状態の純化（purification）を利用し、補助変数に対して最適化を行う。
     • 下界 4 (グローバル純度-エンタングルメント): 下界3と同様であるが、純化されたシステム全体に対して作用する。計算次元は高くなるが、よりタイトな下界を提供する。

主な貢献

• デュアル・バウンド推定: 単一の境界やウィットネスのみを提供する従来のメソッドとは異なり、ENTCALCは下界と上界を同時に計算する。これにより、真の幾何学的エンタングルメントの狭い区間を特定することができ、推定の信頼性を大幅に高める。
• 混合状態およびPPT状態の処理: 本ツールキットは、ネガティビティに基づく尺度では検出できないPPTエンタングル状態（バウンド・エンタングルメント）を含む、混合状態におけるエンタングルメントを定量化することができる。
• 新しい理論的拡張: 著者らは、純度-エンタングルメント相補性関係に基づく、多粒子シナリオのための新しい下界を導出するために既存の理論的枠組みを拡張した。
• ソフトウェア実装: 標準的な科学ライブラリ（NumPy, QuTiP, CVXPY, SciPy, CVX, QETLAB）を活用した、PythonおよびMATLABの両方でのオープンソースパッケージのリリースにより、高精度なエンタングルメント計算へのアクセスを容易にした。

結果と応用
著者らはENTCALCをいくつかの代表的な例に適用し、下界と上界の差が常に非常に小さい（多くの場合 $< 10^{-5}$）ことを観察しており、これは高い精度を示している。

1. PPTエンタングル状態: 本パッケージは、$3 \otimes 3$ ホロデツキ（Horodecki）PPTエンタングル状態の幾何学的エンタングルメントを、境界の差が $\approx 5.15 \times 10^{-6}$ という精度で定量化した。
2. 低ランク状態: 3, 4, 5量子ビット系におけるGHZ状態とW状態の混合物に対するエンタングルメントを定量化した。
3. スピン鎖における熱的エンタングルメント:
   • XXモデル: ペアワイズのエンタングルメントがゼロである場合でも、三量子ビット・ハイゼンベルクXX鎖の熱状態における多粒子エンタングルメントを検出した。強磁性的（$J>0$）および反強磁性的（$J<0$）な結合の両方において、バウンド・エンタングルメント（PPTエンタングル状態）の領域を特定した。
   • XXXモデル（磁場あり）: 磁場強度および逆温度の関数としての幾何学的エンタングルメントの不連続性を通じて、量子相転移のシグネチャーを観察した。
   • 六角形スピン鎖: 外部磁場を調整することで、直接相互作用しない非隣接サイト間（例：6量子ビット環におけるスピンB、D、F）のエンタングルメントを活性化できることを示した。これはペアワイズ解析では捉えられない現象である。
4. ノイズチャネル: GHZ状態およびW状態に対する振幅減衰およびデポラリゼーション・ノイズに対する堅牢性を定量化した。結果として、高ノイズレベルではW状態の方が一般に振幅減衰に対して頑健であるが、低ノイズレベルではGHZ状態の方が初期のエンタングルメントが高いことが確認された。

意義と主張
本論文は、ENTCALCが、多粒子混合状態の幾何学的エンタングルメントを高精度に計算できる最初のパッケージであると主張している。下界と上界の間のタイトな区間を提供することで、本ツールは解析解が困難な理論的研究において信頼できる数値的手法を提供する。著者らは、テストケースにおける境界間の差が極めて小さいことが、採用された境界のタイトさを裏付けていると強調している。彼らは、本パッケージを、バウンド・エンタングルメント、量子相転移、またはノイズ環境下でのエンタングルメント・ダイナミクスといった複雑な現象を調査するための価値あるリソースとして位置づけており、ペアワイズの尺度が見逃してしまう特徴（非局所的な多粒子エンタングルメントなど）を明らかにできる点に注目している。本研究は、NP困難な可分性問題を一般的に解決することを目的とするものではなく、特定の物理的に関連のある状態に対して、実用的かつ高精度な近似を提供するものである。