Extracting conserved operators from a projected entangled pair state

この論文は、無限テンソル積状態（iPEPS）の量子幾何学と忠実度感受性を用いて、その状態を固有状態とする局所的な保存演算子（ハミルトニアン）を高精度で抽出する手法を提案し、短距離 RVB 状態やトポロジカルな励起状態を含む新しいハミルトニアンの発見に成功したことを報告しています。

要約

この論文は、**「量子コンピュータが描いた『絵』から、その絵を描いた『ルール（法則）』を逆算して見つける方法」**を提案した研究です。

少し専門的な用語を、わかりやすい比喩に置き換えて解説します。

1. 背景：逆探偵小説のような問題

通常、物理学では「ルール（ハミルトニアン）」を決めて、そのルールに従って「状態（波関数）」がどうなるかを計算します。
しかし、この研究は逆です。
「すでに完成された『状態（絵）』しか手元にない。じゃあ、この絵を描いた『ルール』は何だったのか？」という逆探偵小説のような問題に挑んでいます。

• 例え話：
  • 料理の味（状態）だけを知っている。
  • その味を作った「レシピ（ルール）」を、味を分析して逆算したい。
  • でも、料理は複雑すぎて、レシピの全貌はわからない。

2. 従来の方法の限界

以前から、この「逆算」を試みる方法（親ハミルトニアンの構築）はありましたが、いくつかの欠点がありました。

• 完璧すぎるルール： 従来の方法は、「この状態が必ず『一番低いエネルギー（地面）』になるように」という厳格なルールを作ろうとしました。
• 複雑すぎる： その結果、ルールが非常に複雑になり、現実の装置で再現するのが難しかったり、必要な相互作用の範囲（局所性）が広すぎて扱いづらかったりしました。

3. この論文の新しいアプローチ：「静かな振動」を聞く

著者たちは、**「静的構造因子（Static Structure Factors）」**という新しい指標を使いました。

• 比喩：
  • 状態（絵）を「静かな湖」だと想像してください。
  • 何かのルール（ハミルトニアン）がその湖に作用すると、湖は「波（振動）」を起こします。
  • もし、そのルールが「保存量（守恒量）」、つまり湖の性質を壊さないものなら、湖は**全く波立たない（振動がゼロ）**はずです。
  • この論文は、**「波立たない（振動がゼロになる）ルール」**を、湖の表面の微妙な揺らぎ（相関関係）を分析することで見つけ出します。

4. 具体的な手法：「変形」させてチェックする

彼らは、以下のような手順でルールを見つけました。

1. パラメータで少し変形する： 手元の「絵（状態）」を、小さなパラメータを使って少しだけ歪ませます（生成関数という手法）。
2. 微分する： 「どのくらい歪んだか？」を計算します。
3. ゼロを探す： もし「歪んでも、元の状態と全く同じ振る舞いをする（エネルギーが変わらない）」ルールがあれば、その計算結果はゼロになります。
4. 行列を解く： この「ゼロになる条件」を満たすルールを、数学的な行列計算（核空間の求解）として見つけ出します。

5. 驚くべき成果：どんな「絵」でもルールが見つかる

この方法は、以下のような難しいケースでも成功しました。

• AKLT 状態（完璧な例）： 理論的に完璧に解ける状態から、正しいルールを再現しました。
• 近似された状態（現実的な例）： 実際の計算では「完璧な絵」ではなく、「少しノイズの混じった近似の絵」しか作れません。通常、ノイズがあるとルールは壊れるはずですが、この方法は**「近似された絵」からも、非常に精度の高いルールを抽出できました。**
• RVB 状態（超伝導の鍵）： 高温超伝導に関わる「RVB 状態」という複雑な絵から、これまでにない**「4 つの点だけで相互作用するシンプルなルール」**を見つけ出しました。従来の方法だと「8 つの点」が必要だったものが、これなら「4 つ」で済むので、実験的に作りやすいルールです。
• 量子多体傷（Quantum Many-Body Scars）： 通常、エネルギーの低い状態（基底状態）がルールに従いますが、この方法で見つけたルールでは、「エネルギーの真ん中にある状態（励起状態）」も、同じルールに従って動いているという奇妙な現象（量子多体傷）を発見しました。これは、量子情報が壊れにくい状態を作るヒントになります。

6. まとめ：なぜこれがすごいのか？

• 「不完全な情報」から「完璧な法則」を引ける： 計算機が近似して作った「少しぼやけた絵」から、物理法則を高精度で逆算できます。
• シンプルで実用的： 従来の複雑なルールではなく、「局所的（近くの粒子同士だけ）」でシンプルなルールを見つけられます。
• 新しい物理の発見： 単に既存のルールを復元するだけでなく、**「新しいタイプの物理現象（量子多体傷など）」**を発見するツールとしても機能します。

一言で言うと：
「量子状態という『完成された作品』を、その微細な揺らぎを分析することで、その作品を支配する『シンプルな法則』を逆探知する、新しい強力なレンズを開発した」という研究です。これにより、量子シミュレータの検証や、新しい量子物質の設計が格段に進むことが期待されます。

技術概要

論文「Extracting conserved operators from a projected entangled pair state」の技術的サマリー

本論文は、与えられたテンソルネットワーク状態（特に無限投影エンタングルメントペア状態：iPEPS）から、その状態が固有状態となる保存演算子（ハミルトニアンを含む）を抽出する新しい手法を提案するものです。従来の「親ハミルトニアンの構築」が抱える課題を克服し、近似計算の誤差が存在する状況下でも高精度な保存演算子を特定できることを示しています。

以下に、問題設定、手法、主要な貢献、結果、および意義について詳述します。

1. 問題設定 (Problem)

量子多体系において、基底状態や特定のテンソルネットワーク状態が与えられたとき、その状態を固有状態とするハミルトニアン（またはより一般的な保存演算子）をどのように復元するかは重要な逆問題です。

• 既存の課題: 1 次元の MPS からは効率的にハミルトニアンを学習する手法が存在しますが、2 次元の PEPS への拡張は困難です。PEPS の正確な縮約は #P 完全であり、物理的に意味のある計算には CTMRG（コーナー転送行列再正規化群）などの近似縮約法が用いられます。これにより相関関数や静的構造因子に制御された誤差が生じ、従来の「親ハミルトニアン構築法（frustration-free 構築）」では、近似誤差の影響で正確な保存演算子を得られない、あるいは局所性が悪化するという問題がありました。
• 目的: 近似縮約法で計算可能な静的構造因子を入力とし、誤差が存在する状況下でも、与えられた iPEPS を（近似）固有状態とする幾何学的に k-局所的な保存演算子を抽出するアルゴリズムの開発。

2. 手法 (Methodology)

提案手法の核心は、静的構造因子行列の核（Kernel）を解くことと、生成関数（Generating Function）の微分による多サイト演算子の評価にあります。

A. 静的構造因子と保存演算子の関係

保存演算子 $\hat{H}$ は、状態 $|\Psi\rangle$ における量子揺らぎがゼロ（または最小）となる演算子として定義されます。

• 幾何学的に k-局所的なエルミート演算子の基底 $\{\hat{o}^\alpha_x\}$ を定義し、その静的構造因子行列 $S^{\alpha,\beta}(q)$ を計算します。
  $$S^{\alpha,\beta}(q) = \sum_x e^{-iq\cdot x} \left( \frac{\langle \Psi | \hat{o}^\alpha_x \hat{o}^\beta_0 | \Psi \rangle}{\langle \Psi | \Psi \rangle} - \frac{\langle \Psi | \hat{o}^\alpha_x | \Psi \rangle \langle \Psi | \hat{o}^\beta_0 | \Psi \rangle}{\langle \Psi | \Psi \rangle^2} \right)$$
• この行列 $S$ の固有値問題 $S\mathbf{h} = s\mathbf{h}$ を解きます。固有値 $s$ がゼロ（または近似ゼロ）に対応する固有ベクトル $\mathbf{h}$ が、保存演算子 $\hat{H} = \sum_x e^{iq\cdot x} \hat{h}_x$ の係数を決定します。
• 量子幾何との関連: 生成関数で定義されるテンソルネットワークの多様体において、静的構造因子行列は量子計量テンソルに比例し、保存演算子は忠実度感受性（Fidelity Susceptibility）がゼロとなる方向に対応します。

B. 多サイト演算子の評価：生成関数法

iPEPS における多サイト演算子の静的構造因子を高精度に評価するために、生成関数法を採用・拡張しました。

• 演算子 $\hat{O} = \sum_x e^{-iq\cdot x} \hat{o}_x$ を、パラメータ $\mu$ を用いた生成関数 $\langle G(\mu, q)| = \langle \Psi | \prod_x (1 + \mu e^{-iq\cdot x} \hat{o}_x)$ の微分として表現します。
• 多サイト演算子の場合、生成関数は iPEPO（無限投影エンタングルメントペア演算子）となり、その結合次元が増加します。
• 数値的微分（中央差分法、5 点公式）を用いて $\mu=0$ での微分値を計算し、CTMRG による iPEPS の縮約と組み合わせることで、静的構造因子を算出します。
• 安定性の工夫: 自動微分（AD）は多サイト演算子において数値的不安定性を示すため、高精度な有限差分法を採用しました。また、トポロジカル秩序状態や猫状態（Cat state）のような非単射（non-injective）iPEPS においては、対称性の破れを適切に扱うための環境固定点の対称化処理も考慮しています。

3. 主要な貢献と結果 (Key Contributions & Results)

A. 2D AKLT 状態（厳密解）へのベンチマーク

• 2 次元 AKLT 状態（スピン 2）に対して、1 サイトおよび 2 サイトの静的構造因子を計算しました。
• 厳密な縮約と CTMRG 近似の両方で、ゼロ固有値を持つベクトルから、既知の対称性生成子（スピン演算子）および 2 サイトのフラストレーションフリーな親ハミルトニアンの完全な集合を正確に抽出することに成功しました。

B. 変分 iPEPS からの非フラストレーションフリー演算子の抽出

• 2 次元 XX モデル（階段状 Z 場）の変分最適化で得られた基底状態（iPEPS）を解析しました。
• 自発的対称性の破れがある相（強磁性相）においても、U(1) 対称性の生成子や、秩序パラメータに対応する保存演算子を高精度で検出しました。
• 臨界点近傍では、元のハミルトニアンの係数を 4 桁の精度で再現する非フラストレーションフリーな保存演算子を抽出しました。これは、従来のフラストレーションフリー構築法では不可能な成果です。

C. RVB 状態からの近似局所親ハミルトニアンの発見

• 短距離共鳴価結合（RVB）状態（臨界状態、結合長発散）に対して、4 サイト・プランケット（4 格子点）の局所演算子を探索しました。
• 得られたハミルトニアンは、RVB 状態を近似固有状態とし、そのエネルギー期待値が変分最適化で得られる基底状態エネルギーと極めて一致していました。
• 重要点: 従来のフラストレーションフリー構築では少なくとも 8 サイトの演算子が必要とされるのに対し、本手法では4 サイトの局所性で良好な親ハミルトニアンを構築できました。

D. 量子多体傷（Quantum Many-Body Scars）の候補

• 変形されたイジング波動関数（および双対的な変形トリックコード状態）に対して、4 サイトプランケットの静的構造因子を解析しました。
• 得られたハミルトニアンは、任意のひずみパラメータ $\beta$ に対して、変形イジング状態（または変形トリックコード状態）を励起固有状態（基底状態ではなく、スペクトルの中央に位置する状態）として持ちます。
• この状態はエンタングルメントエントロピーが面積則を満たすため、**量子多体傷（Quantum Many-Body Scars）**の候補として特定されました。

4. 意義と展望 (Significance & Outlook)

• 近似誤差への頑健性: 本手法は、iPEPS の近似縮約に伴う誤差や、変分最適化による不完全な基底状態であっても、保存演算子を高精度に抽出できることを実証しました。
• 局所性の向上: 従来の親ハミルトニアン構築法に比べて、より高い局所性（短い相互作用範囲）を持つ演算子を発見できます。これは、実験的な量子シミュレータや NISQ デバイスでのモデル実装や検証に極めて有用です。
• 新しい物理の発見: 自発的対称性の破れがある系や臨界系、さらには量子多体傷のような非平衡的な現象を持つモデルの自動発見を可能にします。
• 計算コスト: 静的構造因子の評価コストは iPEPS の結合次元に対して多項式時間であり、実用的な規模で適用可能です。

結論:
本論文は、テンソルネットワーク状態から「データ駆動型」で保存演算子を抽出する強力な枠組みを提供しました。これは、量子物質の低エネルギー有効理論の導出、新しい量子モデルの設計、および量子シミュレータのモデル検証において、重要なツールとなり得ます。特に、近似計算の限界を克服し、より局所的で物理的に意味のあるハミルトニアンを特定できる点は、2 次元量子多体系の研究において画期的な進展です。