Ground state energy and phase transitions of Long-range XXZ using VQE

本論文は、異なる相における基底状態エネルギー誤差の感度を分析することにより、長距離XXZ鎖における無限次相転移を検出するための、制約付きアンザッツ回路を用いた変分量子固有値ソルバー（VQE）の手法を提案し、同時に、厳密対角化の結果に対して当該手法の妥当性を検証するものである。

要約

あなたは、広大で霧に包まれた風景の中で、最も低い地点を探そうとしているところだと想像してください。物理学の世界では、この「最も低い地点」は基底状態エネルギーと呼ばれ、それが粒子（原子磁石など）のシステムが完全に静止したときに、どのように落ち着こうとするのかを教えてくれます。

通常、このような複雑なシステムの最低地点を見つけることは、人間の脳や、あるいは世界で最も強力なスーパーコンピュータにとっても、あまりに巨大なパズルを解くようなものです。ここで、この論文の著者たちは新しいツールを導入しました。それが**変分量子固有値ソルバー（VQE）**です。VQEを、量子コンピュータを使って最低地点の「推測」を行い、古典的なコンピュータを使ってその推測を真実にできるだけ近づくまで洗練させていく、スマートなハイブリッド・ロボットだと考えてください。

課題：2種類の「境界」

研究者たちは、長距離XXZ鎖と呼ばれる特定のモデルを研究していました。小さな磁石（スピン）が互いに影響を与え合えるラインを想像してみてください。通常、磁石は隣接する相手としか対話しませんが、このモデルでは、ライン全体を飛び越えて叫び合う（長距離相互作用）ことができます。

チームは、これらの磁石の振る舞いが劇的に変化する「境界」を見つけたいと考えていました。これらは相転移と呼ばれます。彼らは2種類の境界を見つけました。

1. 「崖」（一次転移）： これは、急な崖から転落するようなものです。エネルギーが突然、鋭く変化します。地面がガクンと落ちるため、非常に見つけやすいものです。
2. 「緩やかな斜面」（無限次転移）： これは、ずっと難しいものです。非常に緩やかで滑らかな丘を登っているようなものです。突然の落下や崖はありません。変化があまりに緩やかに行われるため、標準的なツールではその境界を見落としてしまいます。通常、科学者はこれを見つけるために特別な「グローバルな地図」（複雑な秩序変数）を必要としますが、これは計算が困難です。

秘密兵器：「下手な推測」戦略

ここがこの論文の巧妙な点です。通常、科学者はVQEを単に最低エネルギーの「正確な数値」を得るために使用します。しかし、著者たちは興味深いことに気づきました。**「推測の質は、あなたがどこにいるかに依存する」**ということです。

彼らは、量子ロボット（アンザッツ回路）を、特定のルールに基づいて設計しました。それは、**「全スピン（磁化）を一定に保つ」**というルールです。

• 「正しい」領域にいる場合： ロボットが、磁石が自然にその一定スピン状態を好むような相（フェーズ）にいる場合、ロボットは素晴らしい推測を行います。誤差（ロボットの推測と真の答えとの差）は極めて小さくなります。
• 「間違った」領域にいる場合： もしロボットが、磁石がその状態を望まない相にいる場合、ロボットは苦戦します。ロボットは磁石に間違った形を強制しようとし、その結果、誤差は巨大になります。

「コンパス」の比喩

目に見えない境界を見つけるために、著者たちは単に誤差の大きさを見たのではありません。彼らは、誤差の**「方向」**を見ました。

あなたが森の中を歩いていて、一歩ごとにコンパスを落としているところを想像してください。

• 森のある部分（相A）では、コンパスの針はバラバラな方向を向き、激しく回転しています。
• 別の部分（相B）では、コンパスの針がすべて同じ方向に整然と向いています。

著者たちは、この現象を測定するために**「方向性コヒーレンス（Directional Coherence）」**と呼ばれる手法を用いました。彼らは何千もの地点で「誤差」を計算し、その変化の方向を観察しました。

• コンパスの針が混沌としていたとき、彼らはそこが一つの相であることを知りました。
• 針が突然整列したとき、彼らは境界を越えたことを知りました。

これにより、彼らは簡単な「崖」の境界と、隠れた「緩やかな斜面」の境界の両方を、ロボットの失敗の仕方を観察するだけで特定することができました。彼らは新しい複雑な地図を必要としたのではなく、ただロボットがいかに躓（つまず）くかを観察する必要があったのです。

結果

• 簡単な境界（一次転移）に対して： 誤差が崖のように突然跳ね上がる様子が見られました。
• 難しい境界（無限次転移）に対して： 「コンパスの針（誤差の勾配）」が混沌とした状態から組織化された状態へと変化する様子が見られました。これにより、標準的な手法が見逃していた境界が明らかになりました。
• 精度： ロボットの「脳」（回路のレイヤー）を少し深くすることで、困難なケースであっても、システムの正確なエネルギーを非常に高い精度（誤差3〜7%以内）で計算することができました。

結論

この論文は、システムの変わり目を見つけるために、必ずしも完璧な計算が必要なわけではないと主張しています。時には、「計算がどのように失敗するか」を研究することが、計算そのものよりも、システムの構造について多くのことを教えてくれることがあります。彼らはこの「誤差分析」の手法を用いて、長距離磁性鎖における単純かつ複雑な相転移の両方を描き出すことに成功しました。これは、量子アルゴリズムが単に問題を解くだけでなく、物質の振る舞いの隠れたルールを「発見」するために使用できることを証明しています。

技術概要

技術要約：変分量子固有値ソルバーを用いた長距離XXZ鎖の基底状態エネルギーと相転移

問題提起
本論文は、量子多体系における相転移、特に長距離XXZ（LRXXZ）鎖における相転移の特定という課題に取り組んでいる。変分量子固有値ソルバー（VQE）は、解析解のないハミルトニアンの基底状態エネルギーを見つけるために広く確立されているが、相境界（特に無限次相転移）を検出する際の有用性については、まだ十分に探索されていない。相転移を検出するための従来の手法は、基底状態エネルギーや二点相関関数の導関数における特異点に依存している。これらの手法は有限次の転移には有効であるが、無限次相転移には対応できないことが多い。無限次相転移は、通常、グローバルな範囲を持つ秩序変数や、量子モンテカルロ法のような複雑な手法を必要とする。著者らは、エネルギー最小化のみに伝統的に使用されるVQEを、システムの位相に対するエネルギー誤差の挙動を分析することによって、一次および無限次の相転移境界の両方を特定するように適応させることができることを示すことを目的としている。

手法
本研究では、以下のハミルトニアンで定義される1次元LRXXZ鎖に対してVQEアルゴリズムを採用している：
$$H = -J \sum_{i \neq j} \frac{S^x_i S^x_j + S^y_i S^y_j + \Delta S^z_i S^z_j}{|i - j|^\alpha}$$
ここで、$J$は結合定数、$\Delta$は異方性、$\alpha$は相互作用の範囲を定義する。

1. アンザッツ回路設計： 著者らは、ハミルトニアン・バリエーショナル・アンザッツ（HVA）に着想を得つつ、効率化のために修正を加えた特定の変分アンザッツを設計した。

   • 初期化： システムは、$|000\dots\rangle$状態に交互にパウリXゲートを適用することで準備されるネール状態（$|0101\dots\rangle$）から始まる。
   • 制約： 回路は、最適化の全過程を通じて一定の正味スピン（z方向のゼロ磁化）を維持するように構築されている。
   • 構造： ハミルトニアンは偶数および奇数の結合項に分解される。共有パラメータを持つ標準的なHVAとは異なり、著者らはすべてのゲートを個別にパラメータ化している。これにより、総パラメータ数$4(N-1)$を維持しながら、低い回路深さ（$p=2$）を実現している。
   • ゲートの修正： 標準的な$R_Z$ゲートを用いた初期の実装は効果的ではないことが判明したため、これらは基底状態への収束を改善するために制御$R_Z$（CRZ）ゲートに置き換えられた。

2. 相転移検出戦略：

   • コアとなる仮説は、基底状態エネルギーの精度は、初期アンザッツ状態に対するシステムの位相に依存するということである。
   • 著者らは、厳密対角化（ED）の結果とVQEの結果との間のエネルギー差（誤差）、$E_d = E_{exact} - E_{VQE}$を計算する。
   • 方向的コヒーレンス： 相境界をマッピングするために、著者らはパラメータ空間（$\Delta, \alpha$）におけるこのエネルギー差（$E_d$）の勾配ベクトルを分析する。勾配角度の平均の正弦（sine）と余弦（cosine）の大きさを、移動窓を用いて計算し、「方向的コヒーレンス」を決定する。
   • 論理： アンザッツが適切である位相（例：ネール初期化による反強磁性相）では、誤差は低く、勾配ベクトルは系統的なパターンを示す。アンザッツが適合しない位相（例：強磁性または無秩序な常磁性相）では、誤差は高くなり、勾配ベクトルはランダムになる。これらの挙動の間の境界が相転移をマークする。

主な結果
研究は、$N=12$サイト、開境界条件のシステムに対して実施され、$J=1$（深さ1および2）と$J=-1$（深さ2）の結果を比較した。

1. 一次転移（PMからFMへ）：

   • $\Delta = 1$における常磁性（PM）から強磁性（FM）への転移について、本手法は境界を正常に特定した。
   • 初期ネール状態は正味の磁化がゼロであるため、VQEアンザッツはFM基底状態に到達できない。その結果、システムがFM相に入ると、$E_d$が急激に増加する。
   • 方向的コヒーレンスのプロットは、$\Delta = 1$において勾配ベクトルの向きの急激な変化を示し、平均場理論の予測と一致する形で転移を明確に描き出した。

2. 無限次転移（PMからAFMへ）：

   • 常磁性（PM）相から反強磁性（AFM）相への転移は無限次相転移であり、基底状態エネルギーの導関数のみでは検出が極めて困難である。
   • 著者らは、この転移の間で誤差$E_d$は緩やかに増加するものの、勾配ベクトルの「方向的コヒーレンス」が明確な区別を提供することを観察した。
   • AFM領域では勾配ベクトルはランダムに整列しており、PM領域ではコヒーレントであった。0.7という方向的コヒーレンスの閾値が、これらの領域を正常に分離した。
   • 特定された臨界点は$\Delta = -3$において$\alpha_c \approx 1.75$であり、これは幾何学的エンタングルメント（$\alpha \approx 1.78$）から得られた先行研究の結果とよく一致している。
   • 決定的なことに、著者らはこの転移が標準的な基底状態エネルギーの勾配やエネルギーギャップのプロットでは識別できなかったことを指摘しており、彼らの誤差ベースのアプローチの必要性を検証している。

3. 基底状態エネルギーの精度：

   • 深さ2の回路を用いることで、著者らは基底状態エネルギー計算において高い精度を達成した。
   • $J < 0$の場合、最大相対誤差は約3%であった。
   • $J > 0$の場合、最大相対誤差は7%であった。著者らは、$J > 0$での誤差が高い理由として、追加の修正（ヒルベルト空間を拡張するための$R_y$ゲートの適用など）なしではオプティマイザが基底状態に到達するのが困難であるためであるとしているが、手法自体は依然として有効であった。

意義と主張
本論文は、無限次相転移の境界を特定するためにVQEを具体的に利用した最初の研究であると主張している。主要な貢献は、単なるエネルギー最小化の成功に頼るのではなく、特定の相におけるアンザッツの「失敗」を利用して相境界をマッピングする手法である。

• 新規性： 本手法は、エネルギーの勾配と方向的コヒーレンスの観点からエネルギーデータを分析することで、標準的なエネルギーベースの診断法では通常アクセス不可能な無限次相転移を探索できることを示している。
• 効率性： 提案されたアンザッツは、低い回路深さ（$p=2$）と管理可能なパラメータ数で正確な結果を実現しており、近未来の量子デバイスに適している。
• 汎用性： 著者らは、このアプローチがLRXXZモデルに限定されるものではなく、基底状態の対称性が境界を越えて変化する他のシステムにおける相転移を特定するための新しいパラダイムを提供するものであると示唆している。この手法は、評価される特定の位相に対して敏感なアンザッツを設計することに依存している。

著者らは控えめなトーンを維持しており、ノイズの影響は考慮されていないこと、および方向的コヒーレンスの指標は誤差挙動を分析するためのツールであり、根本的な物理法則ではないことを認めている。彼らは、手法の成功がアンザッツの具体的な設計と初期状態の選択に依存していることを強調している。