⚛️ quantum physics

Preparation of the single-spinon wave function on a quantum computer

本論文は、ハイゼンベルク模型およびハルデン・シャストリー模型の両方における単一スピンン波関数のエネルギーを準備し解析するための、ユニタリ演算子の線形結合に基づく量子アルゴリズムを提示し、計算リソースと性能の観点から3つの異なる戦略を比較するものである。

原著者： D. Faílde, A. Gómez, J. Fernández-Rossier

公開日 2026-09-02

📖 1 分で読めます🧠 じっくり読む

CC BY 4.0

原著者： D. Faílde, A. Gómez, J. Fernández-Rossier

原論文は CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

量子材料の微視的な世界において、原子はしばしば長い一次元の鎖状に配列し、隣接する粒子と相互作用する小さな磁石のように振る舞います。数十年にわたり、物理学者たちはこれらの中性子鎖を研究することで、最も根本的なレベルで磁性がどのように機能するかを理解しようとしてきました。この分野における重要な謎の一つは、相互作用が完璧に均衡し、量子的な無秩序状態を作り出している特定のタイプの磁性鎖に関するものです。このような系において、集団励起（鎖が乱された際の反応）は、音波や光波のような馴染みのある波としては振る舞いません。その代わりに、理論は、これらの励起がより小さく独立した断片へと分裂することを示唆しています。これらの断片は「スピノン」と呼ばれます。完全な磁気波とは異なり、スピノンは単一の原子が持つ磁気強度の半分しか持ちません。これは古典的な世界では不可能に思える性質ですが、量子の領域では極めて正常なことです。これらのスピノンがどのように動き、エネルギーを運ぶのかを理解することは、エキゾチックな材料の挙動を把握する上で極めて重要ですが、それらは現実の世界において孤立した物体として存在しないため、物理的な実験で直接観察することは非常に困難です。

この問題を解決するために、研究チームは量子コンピューティングの台頭する力に目を向けました。物理的な原子の鎖を構築する代わりに、彼らはデジタル量子プロセッサを使用して、これらの磁性鎖の振る舞いをシミュレートし、単一のスピノンの数学的な設計図を構築しました。研究者たちは、原子がすぐ隣の隣人とのみ相互作用するモデルと、距離に関係なくすべての原子が鎖内の他のすべての原子と相互作用するモデルという、2つの特定の磁性鎖モデルに焦点を当てました。彼らの目標は、単一のスピノンを表す量子状態を準備する方法を見出すことでした。これは、これまで理論的な計算や古典的なコンピュータによるシミュレーションによって制限されてきた課題でした。デジタル量子プロセッサ上でこの状態を生成することに成功したことで、彼らは、物理的なラボでそれらを孤立させる必要なく、これらの捉えどころのない粒子の特性（エネルギーや運動量など）を探求するための新しい手法を実証しました。

彼らが開発したプロセスは、磁性鎖の基底状態を準備することから始まる多段階の手順です。この基底状態とは、システムが保持できる最低エネルギーの構成であり、スピンが均衡している安定した配置です。この安定した状態が量子コンピュータ上に確立されると、研究者は追加の磁気ユニットを導入し、事実上、鎖に原子を一つ追加します。この追加が完璧な均衡を乱し、スピノンが存在するための条件を作り出します。しかし、単に原子を追加するだけでは不十分です。研究者は、この追加ユニットが定義された運動量を持つように、鎖全体に広がった特定のパターンを作り出す必要がありました。これを達成するために、彼らは「ユニタリ演算子の線形結合」と呼ばれる手法を用いました。簡単に言えば、これは量子コンピュータに対して、最終的な状態を一度のステップで強制的に作り出そうとするものではありません。代わりに、追加のスピンが鎖の異なる位置に存在する多くの可能性の重ね合わせを作成し、それらの可能性を精密な数学的手法で組み合わせることで、最終的なスピノンの状態を形成したのです。

チームはこの手法を2種類の異なる磁性鎖でテストしました。原子が隣接する相手としか対話しない鎖については、基底状態を見つけるために「変分量子固有値ソルバー」として知られる柔軟で反復的なアプローチを用いました。長距離相互作用を持つ鎖については、特定の数学的射影に基づく、より直接的で厳密な手法を用いました。どちらの場合においても、彼らは単一スピノンの状態を生成し、その特性を測定することに成功しました。彼らは、スピノンのエネルギーが運動量に応じて変化する様子が、数十年前の理論式と完全に一致することを発見しました。これにより、彼らのデジタル構築が正確であることが確認されました。また、彼らは「ノルム」と呼ばれる量を測定しました。これは、状態がいかに「実在的」で明確であるかを示すものです。彼らの結果は、スピノンの状態が特定の運動量の範囲内にのみ存在し、その範囲外では消失することを示しており、これはスピノンが量子の風景における特定の領域に限定されているという理論的理解と一致しています。

直接的な状態生成の手法は有効ですが、機能するために「アンシラ量子ビット」として知られるかなりの数の追加量子ビットを必要とし、操作の成功は確率的である（つまり、毎回必ず成功するわけではない）という側面があります。これに対処するため、研究者たちはより少ないリソースを必要とする代替戦略を模索しました。一つのアプローチは、問題を管理可能な小さな断片に分解し、それらを個別に測定することで、スピノンのエネルギーと特性を算出することでした。これは、状態全体を一度に作成しようとするのではなく、単一の追加量子ビットを用いてシステムを探索する「アダマール・テスト」と呼ばれる手法です。このアプローチは、同じ量のデータを収集するために非常に多くの個別の回路を実行する必要がありますが、必要なハードウェアの観点からははるかに効率的です。研究者たちはこれらの異なる戦略を比較し、直接的な手法は強力であるものの、量子ビット間の接続が限られている現在の世代の量子コンピュータには、ステップバイステップの測定アプローチの方が適していることを見出しました。

この研究はすべてコンピュータ・シミュレーションを通じて行われました。つまり、研究者たちは量子プロセッサの振る舞いを模倣する古典的なコンピュータ上でアルゴリズムを実行しました。彼らは実際の物理的な量子ハードウェア上でこれらの実験を行ったわけではありませんが、シミュレーションは手法の実現可能性をテストするのに十分な厳密さを備えていました。結果は、提案された手法が実行可能であり、単一スピノンの波動関数を正確に準備し分析できることを示しています。研究者たちはさらに、より複雑な長距離相互作用を持つ「ハルデン・シャストリー・モデル」へと研究を広げ、同じ手法が成功裏に適用できることを見出しました。単純な短距離鎖にうまく機能する数学的ツールが、より複雑な長距離系に適用されると失敗することが多い中で、これは重要な進展です。量子アルゴリズムが両方のタイプのモデルを扱えることを示すことで、チームはより複雑な磁性材料の研究への道を開きました。

究極的には、本研究は量子コンピュータ上で複雑な量子粒子を研究するためのロードマップを提供しています。単一のスピノンの状態を準備し測定する能力は、これらの粒子が二次元材料の中でどのように振る舞うかといった、より複雑な現象を理解するための基礎的なステップです。現実の世界では、ベテ・アンザッツや密度行列繰り込み群のような一次元鎖に有効な数学的手法は、二次元系に適用されると非効率になるか、あるいは完全に機能しなくなります。量子シミュレーションは、このボトルネックに対する潜在的な解決策を提供します。研究者たちは、これらの手法がシミュレーション環境において有効であることを示すことで、量子コンピュータが将来的に、高次元におけるスピノンの波動関数やダイナミクスを探求するための主要なツールとなり、古典的なコンピュータでは到達不可能な量子物質の振る舞いへの新しい窓を開くことを示しました。

技術要約：量子コンピュータにおける単一スピンオン波動関数の準備

問題提起
本論文は、量子コンピュータ上で1次元 $S=1/2$ スピンモデルの単一スピンオン波動関数を準備するという課題に取り組んでいる。反強磁性的相互作用を持つハイゼンベルク模型は、ギャップレスなスピン励起スペクトルを持つが、これらは孤立したマグノンではなく、スピンオン（ $S=1/2$ ）と呼ばれる複合粒子として記述されることが多い。Kulkaら[1]によって提案された単一スピンオン波動関数のヒューリスティックなアンザッツは、 $L$ サイトの鎖の基底状態構成に $S_z=1/2$ のサイトを1つ追加することで、 $L+1$ サイトの鎖における $S=1/2$ の状態を構築するものである。しかし、このアンザッツを量子ハードウェア上で実装することは容易ではない。なぜなら、得られる状態は正規化されておらず、単一のユニタリ演算子によって生成することができず、ユニタリの線形結合（LCU）を必要とするためである。さらに、このアンザッツはXXZ模型に対しては検証されているものの、長距離相互作用を持つハルダン・シャストリー（Haldane-Shastry）模型への適用可能性は確立されていなかった。

手法
著者らは、図1に示される3つの主要な段階からなるプロトコルを提案している。

基底状態の準備： ターゲットとなるスピン模型（ハイゼンベルクまたはハルダン・シャストリー）の基底状態（ $|\Psi_{GS}\rangle$ $∣ Ψ_{GS} ⟩$ ）を $L$ $L$ 量子ビットで準備する。
- ハイゼンベルク模型の場合： 変分量子固有値ソルバー（VQE）を利用する。彼らは、バレンスボンド結晶ではなく、ギブス回転（Givens rotations）によって得られた自由フェルミオン状態を初期状態とする、浅いパラメータ化量子回路（PQC）を用い、その後に $R_{ZZ}$ および $R_{XX+YY}$ ゲートの層を適用する。
- ハルダン・シャストリー模型の場合： 基底状態がグッツウィラー射影された自由フェルミオン海であることを利用する。彼らはギブス回転を用いて自由フェルミオン状態を準備し、一連のCNOTゲートを通じてグッツウィラー射影を適用する。このステップは非決定的であり、成功確率はシステムサイズに対して指数関数的に減少する。
単一スピンオン状態の構築： 格子を1サイト拡張し（ $|\uparrow\rangle$ 状態のアドアンサ量子ビットを追加）、スピンオン・アンザッツを適用する。アンザッツは、拡張された構成のフーリエ変換として定義される： $|\Psi(q)\rangle = \frac{1}{\sqrt{L+1}} \sum_{m=0}^{L} e^{iqm} |\Psi(m)\rangle$ 。ここで、 $|\Psi(m)\rangle$ は、余剰のスピンがサイト $m$ に挿入された基底状態を表す。
エネルギーおよびノルムの推定： 著者らは、運動量依存のノルム $N(q)$ $N (q)$ とエネルギー分散 $\epsilon(q)$ $ϵ (q)$ を計算するために、3つの異なる戦略を評価している：
- ユニタリの線形結合（LCU）： $L$ 個のアドアンサ量子ビットを用いて、直接重ね合わせ状態 $|\Psi(q)\rangle$ を準備する非決定的アルゴリズム。ノルムは、アドアンサレジスタを $|0\rangle^{\otimes L}$ で測定する確率によって推定される。
- SWAPストリングの期待値： アドアンサを用いずに、基底状態上で直接オーバーラップ $\langle \Psi(n) | \Psi(m) \rangle$ を計算する決定的アプローチ。しかし、これには指数関数的な数のパウリ文字列（ $O(4^L)$ ）の測定を必要とする。
- アダマールテスト： オーバーラップと遷移振幅を項ごとに推定するために、単一のアドアンサ量子ビットを使用する決定的手法。これは $O(L^2)$ 回の回路を必要とするが、パウリ文字列の指数関数的なスケーリングや、LCU法のような高いアドアンサ数を回避できる。

主な貢献

ハルダン・シャストリー模型への拡張： 著者らは、もともとXXZ模型のために開発されたKulkaらの単一スピンオン・アンザッツが、ハルダン・シャストリー模型におけるスピンオンを正確に記述することを実証した。彼らは、このアンザッツが解析的なエネルギー分散 $\epsilon(q) = \frac{J}{2}(q_0^2 - q^2)$ と、ノルムが非ゼロとなる運動量領域（ $-\pi/2 < q < \pi/2$ ）を正しく再現することを確認した。
量子実装戦略： 本論文は、量子ハードウェア上で単一スピンオン・アンザッツを実装するための詳細なロードマップを提供している。彼らは、LCU、SWAPストリング、およびアダマールテストの手法の計算リソース（量子ビット、ゲート、回路の深さ、回路数）を比較している。
イン・シリコ（in-silico）検証： 提案された手法は、システムサイズ最大 $L=16$ （ハイゼンベルク）および $L=8$ （ハルダン・シャストリー）に対して数値シミュレーション（in-silico）を通じてテストされている。シミュレーションにより、アダマールテストおよびLCUアプローチが厳密対角化法と一致する結果を与えることが確認された。

結果

基底状態の準備： ハイゼンベルク模型に対するVQEアプローチは、中程度の鎖長（ $L \leq 8$ ）において低い不完全度（infidelity）を達成しており、性能は $L$ の増加に伴いわずかに低下するものの、サンプリングノイズの閾値以下に留まっている。ハルダン・シャストリー模型については、グッツウィラー射影法は高い忠実度を得られるが、成功確率が指数関数的に減少する（ $p \sim 2^{-\sqrt{2}L/2}$ ）ため、中程度の鎖長への実用的な適用が制限される。
スピンオン分散： 計算されたエネルギー分散 $\epsilon(q)$ およびノルム $N(q)$ は、両方の模型において解析的な予測と一致している。ノルム $N(q)$ は、 $L$ が大きくなるにつれて $|q| > \pi/2$ で正しく消失し、スピンオンの運動量領域を確認している。
手法の比較：
- LCU法は直接的な状態準備を可能にするが、 $L$ 個のアドアンサ量子ビットと非決定的なポストセレクションステップを必要とする。その成功確率は $N(q)/(L+1)$ に対して線形にスケールするため、基底状態のグッツウィラー射影よりもスケーラブルである。
- アダマールテストは、単一のアドアンサ量子ビットのみを必要とし、決定的である。これは $O(L^2)$ 回の回路実行を要求するが、直接的なSWAPストリング法に見られるパウリ項の指数関数的なスケーリングや、LCUによる高い量子ビットオーバーヘッドを回避できる。著者らは、ハルダン・シャストリー模型においては、基底状態準備の成功確率の指数関数的減少が誤差予算を支配するため、精度を維持するためにショット数の指数関数的な増加が必要になると指摘している。

意義および主張
本論文は、物理的な量子系において単一スピンオン波動関数を実現するための実用的なフレームワークを提供することを主張している。これは、古典的な手法（ベテ解法やDMRGなど）が非効率になる領域、特に高次元における準粒子の研究において不可欠なステップである。著者らは、これらの手法は現在「in-silico」でテストされていることを明示しており、物理的なハードウェア上での実験的実現を主張しているのではない。主な貢献は、特定の非自明なスピン波動関数を準備するためのアルゴリズムの提案とリソース分析であり、これは、厳密解が存在しない二次元モデルへの将来的な実装に向けた前駆的な研究として機能する。本研究は、近年の量子ハードウェアの制約という文脈において、決定的な項別推定（アダマールテスト）と確率的な状態準備（LCU）の間のトレードオフを浮き彫りにしている。

関連論文