[論文レビュー] Spin relaxation dynamics with a continuous spin environment: the dissipaton equation of motion approach

Wenxiang Ying, Su Yu+3 more|arXiv (Cornell University)|Feb 1, 2023

Spectroscopy and Quantum Chemical Studies被引用数 4

ひとこと要約

本稿は、連続的スピン浴におけるスピン緩和ダイナミクスのための散逸関数運動方程式（DEOM）手法を開発し、任意のスピン量子数 $S$ に対してフラクチュエーション・ディスシペーション定理（FDT）を導出し、時間領域プロニ適合（time-domain Prony fitting）を用いて指数的減衰基底関数を構築する。この手法により、標準的なスピン・ボソンモデルを超える非ガウス的・非調和的環境の高精度かつ高効率なシミュレーションが可能となる。

ABSTRACT

We present the quantum dynamics of a spin coupling to a bath of independent spins via the dissipaton equation of motion (DEOM) approach. The bath, characterized by a continuous spectral density function, is composed of spins that are independent level systems described by the su(2) Lie algebra. This represents an extreme class of anharmonic environment. Based on the conclusion drawn by Suarez and Silbey [J. Chem. Phys. 95, 9115 (1991)] and Makri [J. Chem. Phys. 111, 6164 (1999)] that the spin bath can be mapped to a Gaussian environment under its linear response limit, we derive the fluctuation-dissipation theorem (FDT) of the spin bath from a microscopic perspective, and generalize the discussion to the case of arbitrary bath spin quantum number S. Next, the time-domain Prony fitting decomposition scheme is applied to the bare-bath time correlation function (TCF) given by FDT to generate the exponential decay basis (or pseudo modes) for DEOM construction. The accuracy and efficiency of this strategy has been justified by a variety of numerical results. We envision this work provides new insights to extend the hierarchical equations of motion (HEOM) and DEOM approach to certain types of anharmonic enviroments with arbitrary TCF or spectral density

研究の動機と目的

• 2準位系が任意のスピン量子数 $S$ を持つ連続的スピン浴に結合する量子ダイナミクスフレームワークの構築。
• 微視的視点から、このようなスピン浴に対するフラクチュエーション・ディスシペーション定理（FDT）を導出することにより、非ゼロの $S$ への先行研究の一般化を図る。
• 裸浴の時間相関関数（TCF）の時間領域プロニ適合（t-PFD）による近似を用いて、指数的減衰基底を構築することで、DEOM形式を非調和的環境へ効果的に拡張すること。
• さまざまな浴パラメータと温度領域において、手法の正確性と効率性を検証すること。

提案手法

• 任意のスピン量子数 $S$ を持つ独立した $\mathfrak{su}(2)$-代数的レベル系で構成されるスピン浴に対するFDTを導出し、有効スペクトル密度 $J_{\text{eff}}(\omega;\beta,S)$ を得る。
• 時間領域プロニ適合分解（t-PFD）を用いて、裸浴のTCFを指数的減衰の和として近似し、DEOM構築の基底を形成する。
• 得られた擬似モード表現を用いて、非マークフィアンのダイナミクスを記述するための散逸関数階層を構築し、非調和的スピン浴の下でのダイナミクスを記述する。
• 線形応答極限と有効スペクトル密度マッピングを活用することで、任意のスペクトル密度および非ガウス的環境へとアプローチを一般化する。
• 既知のベンチマークと比較することで、数値的検証を実施し、特にスピン浴がガウス的（ボソン的）浴に還元される高スピン極限 $S \gg 1$ での有効性を確認する。
• 浴のカットオフ周波数や結合強度が異なるモデル系にフレームワークを適用し、収束性と計算コストを評価する。

Figure 1: Population dynamics of the zero-temperature SSB models. The spin baths are parameterized as (a) $\alpha=0.5$ , $\omega_{c}/\Delta=1$ . (b) $\alpha=0.1$ , $\omega_{c}/\Delta=6$ . (c) $\alpha=0.2$ , $\omega_{c}/\Delta=10$ . (d) $\alpha=0.5$ , $\omega_{c}/\Delta=10$ . (e) $\alpha=0.75$ , $\om

実験結果

リサーチクエスチョン

• RQ1 任意のスピン量子数 $S$ を持つスピン浴に対して、どのように微視的視点からフラクチュエーション・ディスシペーション定理を導出できるか？
• RQ2 時間領域プロニ適合戦略が、異なる $S$ と温度におけるスピン浴の裸浴時間相関関数をどの程度正確に表現できるか？
• RQ3 DEOM形式は、標準的なボソン的浴の仮定を超えて、非ガウス的・非調和的環境（例：スピン浴）を効果的に記述できるか？
• RQ4 t-PFDに基づくDEOMの正確性は、指数的項の数や浴パラメータ $\alpha$ および $\omega_c/\Delta$ にどのように依存するか？
• RQ5 有限温度スピン浴に対して、異なるt-PFDスキームを用いる際の計算コストと正確性のトレードオフはいかなるものか？

主な発見

• スピン浴のFDTは $ C(t) = \frac{1}{\pi} \int_{-\infty}^{\infty} d\omega \, e^{-i\omega t} \frac{J_{\text{eff}}(\omega;\beta,S)}{1 - e^{-\beta\omega}} $ として導出され、ここで $ J_{\text{eff}}(\omega;\beta,S) = J(\omega) \zeta(\omega;\beta,S) $ であり、$ \zeta(\omega;\beta,S) $ はスピン統計を反映する。
• $ S = 1/2 $ の場合、$ \zeta(\omega;\beta,1/2) = \tanh(\beta\omega/2) $ となり、既存の文献で知られている結果と一致する。
• 高スピン極限 $ S \gg 1 $ において、ホルシュタイン＝プリマコフ変換によりスピン浴はボソン的浴と同型となり、標準的なスピン・ボソンモデルが回復される。
• $ 4+4 $ 項の t-PFD 戦略は、有限温度モデルのTCFの実部および虚部を高精度に適合可能であり、視覚的検証により高い正確性が確認された。
• 中程度の指数的項数で高い正確性を達成しており、非調和的環境への計算的効率性とスケーラビリティが示された。
• 数値的結果から、高カットオフ周波数や強い結合（$ \alpha = 10 $）を含むさまざまな浴パラメータに対しても、手法のロバスト性が確認され、有限温度でも正確性を維持していることが示された。

Figure 2: Convergence test of the localization model ( $\omega_{c}/\Delta=1,\ \alpha=10$ ) using $t$ -PFD strategies $2+2$ , $3+3$ and $4+4$ . Comparisons are made against ML-MCTDH (digitized from Ref. 33 ).