[论文解读] 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

一句话总结

本文提出了一种精确的解析方法，用于在任意时变波形下求解量子自旋系统的时间演化，方法结合代数图论与路径求和形式化。通过将自旋动力学表示为图上的行走路径，并利用收敛的连分数计算时序指数，该方法可获得闭式解，其在单组分、双组分及三组分自旋系统中，均优于传统的常微分方程求解器和分段常数近似方法，在精度与效率方面表现更优。

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.

研究动机与目标

• 开发一种通用的解析框架，用于求解在任意时变波形下的量子自旋系统的时间演化问题。
• 克服数值常微分方程求解器与分段常数近似方法在模拟复杂波形下自旋动力学时的局限性。
• 利用图论与组合方法，为自旋系统中的时序指数提供精确的闭式解。
• 展示该方法在核磁共振（NMR）与自旋-自旋耦合系统中计算性能与精度方面的优越性。
• 将路径求和形式化方法拓展至具有多种哈密顿结构（包括标量耦合）的多自旋系统。

提出的方法

• 该方法采用路径求和形式化，将哈密顿量的时间有序指数表示为以哈密顿量为邻接矩阵的图上所有行走路径之和。
• 使用伏尔泰拉复合运算替代求解式中的标准矩阵乘法，从而通过连分数实现精确解。
• 通过时间离散化的伏尔泰拉复合运算实现数值计算，仅需对条件良好且为三角形的矩阵进行乘法与求逆运算。
• 对于大规模系统，应用Lanczos路径求和方法通过三对角化降低矩阵规模，同时保持时间有序指数的精度。
• 通过在NMR中使用六参数啁啾脉冲对方法进行验证，并与ode45及分段常数传播器近似（PCPA）进行比较。
• 直接计算完整的3×3布洛赫方程系统，从而无需重新运行模拟即可评估任意初始态下的结果。

实验结果

研究问题

• RQ1 能否利用代数图论，为任意时变波形下的量子自旋动力学推导出精确的解析解？
• RQ2 与标准常微分方程求解器及分段常数近似方法相比，路径求和方法在精度与计算效率方面表现如何？
• RQ3 路径求和形式化在多自旋系统中，特别是具有复杂耦合结构的系统中，其通用性可扩展到何种程度？
• RQ4 通过矩阵预处理（如Lanczos路径求和）能否在降低计算成本的同时保持高精度？
• RQ5 在涉及啁啾脉冲与标量耦合的真实NMR场景中，路径求和方法的性能如何？

主要发现

• 路径求和方法在精度与计算时间方面始终优于ode45与PCPA，尤其在高精度要求下表现更优。
• 在单组分系统中，路径求和的梯形与辛普森规则在7秒内实现相对误差低于10⁻⁸，而PCPA在10⁻⁸精度下需超过12秒。
• 在双组分系统中，路径求和方法在10⁻⁸相对误差下将计算时间减少高达90%，其中PS-Simpson仅需1.69秒，而PCPA需12.47秒。
• 路径求和方法提供了以‹-解析式表示的精确解析表达式，可通过收敛的诺伊曼级数展开实现任意精度计算。
• 该方法具备尺度不变性，并可与Lanczos预处理兼容，从而在较大自旋系统中实现高效计算。
• 与纯数值方法（如常微分方程积分或PCPA）不同，路径求和形式化可提供对时间演化的解析洞察。