Quantum annealing through a first-order phase transition: field theory approach

技術要約：一次相転移を通じた量子アニーリング：場の手法によるアプローチ

問題提起
量子アニーリングは、ハミルトニアンを緩やかに変化させる際に、システムがその瞬間の基底状態に従うことに依存している。二次相転移における欠陥生成は、キブブル・ズレック・メカニズム（励起数のべき乗則スケーリング $n_{ex} \propto \beta^\nu$ を予測する）を通じてよく理解されているが、実験的な観測では、スイープ速度 $\beta \to 0$ としても消えない、エラー率における上昇した一定のオフセット $n_0$ がしばしば見られる。本論文は、この持続的な欠陥生成の背後にあるメカニズムを解明することを目的としており、それが一次量子相転移に起因すると仮定している。すなわち、ポテンシャルの障壁が消失する前に、システムがメタステーブルな「偽の真空」にトラップされ、結果として大きな集団振動と計算エラーを引き起こすという仮説である。

手法
著者らは、欠陥の動力学をモデル化するために、量子場理論の手法を用いている。

1. 場理論の定式化: 長距離相互作用を持つ相互作用イジングスピンの動力学を、外部の時間依存力を持つ $0+1$ 次元の実スカラー $\phi^4$ 理論へと還元している。ラグランジアンは $L = \frac{m}{2}\dot{\phi}^2 - V(\phi, t)$ で定義され、ここでポテンシャル $V(\phi, t) = -\beta t \phi - \frac{k}{2}\phi^2 + \frac{g}{4}\phi^4$ である。$k$ は転移の性質を決定し、$k>0$ は一次転移、$k=0$ は臨界境界、$k<0$ は劣臨界通過に対応する。
2. モデルのマッピング: $N$ 個のスピンを代表する最小のモデルとして、駆動型リプキン・メシュコフ・グリック・モデル（LMGm）を用いる。著者らは、大 $N$ 極限におけるLMGmのハミルトニアンパラメータ（$N, D, B_x$）と、場理論のパラメータ（$m, k, g$）の間の対応関係を導出している。
3. 半古典的解析: 大 $S$ （スピンサイズ）極限において、量子動力学は時間依存ポテンシャル中の粒子の古典運動へとマップされる。励起数 $n_{ex}$ は、$t \to -\infty$ から $t \to +\infty$ までの断熱不変量 $I = \frac{1}{2\pi} \oint p \, dx$ の変化を通じて計算される。
4. 解析および数値的手法:
   • 劣臨界の場合（$k<0$）、移動する最小値付近での調和振動子近似から生じる積分に対して、最急降下法を適用する。
   • 超臨界の場合（$k>0$）、ポテンシャルの最小値が消失する臨界点付近の動力学を分析する。著者らは、臨界時刻付近において、運動方程式が可積分なパレレーヴェ・I（Painlevé-I）方程式に帰着することを見出した。
   • LMGm ($N=60$) および古典場方程式の数値シミュレーションを行い、解析的なべき乗則と前係数の妥当性を検証している。

主な貢献と結果
本論文は、スイープ速度 $\beta$ および臨界点からの距離 $k$ の関数としての励起数 $n_{ex}$ に関する具体的なべき乗則を導出している：

1. 劣臨界領域 ($k < 0$):
   励起数は $\beta \to 0$ において指数関数的に抑制される：
   $$n_{ex} \propto \beta^{\nu_-} e^{-\Delta/\beta}, \quad \text{ここで } \Delta \propto |k|^{\mu_-}$$
   導出された指数は $\nu_- = 1$ および $\mu_- = 2$ である。これは $\mu_- = 3/2$ を示唆していた平均場理論の予測とは対照的である。

2. 臨界境界 ($k = 0$):
   システムは純粋なべき乗則を示す：
   $$n_{ex} \propto \beta^\nu, \quad \text{ただし } \nu = 3/4.$$

3. 一次転移 ($k > 0$):
   励起数は一定のオフセットを伴うスケーリング則に従う：
   $$n_{ex} \approx n_0 + c\beta^{\nu_+}$$

   • オフセット $n_0$ は $n_0 \propto k^{\mu_+}$ であり、$\mu_+ = 3/2$ である。
   • 指数 $\nu_+$ は、一定のオフセットが支配的な遅い転移の極限において、約 $0.8$（具体的には$4/5$）であることが判明した。これは、臨界境界や高速スイープで見られる$3/4$ という値とは異なる。
   • 著者らは、$\nu_+ \approx 0.8$ という指数は、システムが臨界点付近で異常に長い時間を過ごすこと（臨界スローイングダウン）により、システムが新しいグローバル最小値に落下する前にポテンシャルの井戸が修正されるために生じると指摘している。

4. 可積分性: 著者らは、完全な量子ハミルトニアンおよび一般的なラグランジアンは非可積分であるが、一次転移における臨界モーメント付近の有効な場動力学は、可積分なパレレーヴェ-I方程式によって記述されることを観察している。

意義と主張
著者らは、自らの場理論が、量子アニーリング実験（D-Waveによるものなど）で観察される「一定のオフセット」に対する普遍的な説明を提供すると主張している。これは、実験的な不完全さのみによるものではなく、一次臨界点に起因するものであるとしている。

• シグネチャの特定: べき乗指数の急激な変化と、$\beta \to 0$ における非零の $n_0$ の存在は、一次臨界点のシグネチャとして機能する。
• アルゴリズムへの影響: 量子アニーリング・アルゴリズムの性能は、それが通過する相転移の普遍性クラスによって特徴付けられる可能性がある。一次臨界点（例えば、不要なエネルギー障壁を生み出すパラメータの不確定性を抑制することによって）を特定し回避することで、エラーのスケーリングを一定のオフセットから指数関数的に抑制される領域へと変換し、計算の収束性を向上させることができる。
• 普遍性: 結果は、イジングスピン・ハミルトニアンにマップされる計算問題に共通する特徴である、「よく混合された」（長距離相互作用を持つ）系に適用可能であるとして提示されている。

本論文は、既存の実験的な異常を説明できるものの、これらの知見に基づいた新しい、より高速なアルゴリズムの開発は、依然として投機的な将来の方向性であると控えめに結論付けている。