[論文レビュー] Exact solutions for the time-evolution of quantum spin systems under arbitrary waveforms using algebraic graph theory

Pierre-Louis Giscard, Mohammadali Foroozandeh|arXiv (Cornell University)|May 10, 2022

Spectroscopy and Quantum Chemical Studies参考文献 38被引用数 2

ひとこと要約

本稿では、代数的グラフ理論とパス和形式を用いて、任意の時間依存波形下での量子スピン系の時間発展を正確に解析的に解く手法を提示する。スピンダイナミクスをグラフ上のウォークとして表現し、収束する連分数を用いて時間順序指数を計算することで、単一部、二部、三部スピン系の全範囲で、従来のODEソルバーや区分定数近似法よりも精度と効率に優れた閉形式解を提供する。

ABSTRACT

A general approach is presented that offers exact analytical solutions for the time-evolution of quantum spin systems during parametric waveforms of arbitrary functions of time. The proposed method utilises the \emph{path-sum} method that relies on the algebraic and combinatorial properties of walks on graphs. A full mathematical treatment of the proposed formalism is presented, accompanied by an implementation in extsc{Matlab}. Using computation of the spin dynamics of monopartite, bipartite, and tripartite quantum spin systems under chirped pulses as exemplar parametric waveforms, it is demonstrated that the proposed method consistently outperforms conventional numerical methods, including ODE integrators and piecewise-constant propagator approximations.

研究の動機と目的

• 任意の時間依存波形下での量子スピン系の時間発展を解く一般化された解析的フレームワークの構築を目的とする。
• 複雑な波形を伴うスピンダイナミクスのシミュレーションにおいて、数値的ODEソルバーや区分定数近似法の限界を克服することを目的とする。
• グラフ理論的および組合せ的手法を用いて、スピン系における時間順序指数の正確な閉形式解を提供することを目的とする。
• NMRおよびスピン-スピン結合系において、計算性能と精度の面で本手法の優位性を示すこと。
• スカラー結合を含む多スピン系や多様なハミルトニアン構造を有する系へのパス和形式の適用範囲を拡張すること。

提案手法

• 本手法はパス和形式を採用しており、ハミルトニアンを隣接行列とするグラフ上のすべてのウォークの和として、ハミルトニアンの時間順序指数を表現する。
• リゾルベント定式化における標準的行列乗算の代わりにボルテラ合成を用い、連分数による正確な解を可能にする。
• 数値的実装では、時間量子化されたボルテラ合成を用い、良好な条件数を持つ三角行列の乗算と逆行列計算のみを必要とする。
• 大規模系では、三重対角化を介して行列サイズを縮小しながらも、時間順序指数の精度を保持するラントツパス和法を適用する。
• NMRにおける6パラメータのチープドパルスを用いた妥当性の検証を行い、ode45および区分定数プロパゲータ近似（PCPA）と比較した。
• 全3×3ブロッホ方程式系に対して直接解を計算することで、初期状態を再計算せずに任意の初期状態の評価が可能となる。

実験結果

リサーチクエスチョン

• RQ1 代数的グラフ理論を用いて、任意の時間依存波形下での量子スピンダイナミクスに対して正確な解析的解を得ることは可能か？
• RQ2 標準的ODEソルバーや区分定数近似法と比較して、パス和法の精度と計算効率はいかがなっているか？
• RQ3 複雑な結合構造を有する多スピン系へ、パス和形式をどの程度一般化できるか？
• RQ4 ラントツパス和法のような行列前処理を用いることで、計算コストを低減しつつも高い精度を維持できるか？
• RQ5 チープドパルスおよびスカラー結合を含む現実的なNMR状況において、パス和アプローチの性能はいかがなものか？

主な発見

• パス和法は、特に高精度要件下において、ode45およびPCPAを常に上回る精度と計算時間の両面で優位性を示した。
• 単一部系では、パス和の台形則およびシンプソン則が7秒未塔で相対誤差10⁻⁸未満を達成したのに対し、PCPAは10⁻⁸精度で12秒以上を要した。
• 二部系では、相対誤差10⁻⁸において、パス和法がPCPA比で計算時間を最大90％まで短縮した。PS-Simpsonは1.69秒、PCPAは12.47秒であった。
• パス和アプローチは、収束するネウマン級数展開により任意の精度で計算可能な‹-リゾルベントの関数として正確な解析的表現を提供する。
• スケール不変性とラントツ前処理との相性の良さにより、より大きなスピン系に対しても効率的な計算が可能である。
• パス和形式は、ODE統合やPCPAといった純粋な数値的手法とは異なり、時間発展に関する解析的洞察を提供する。