[论文解读] Advanced Optimal Control Methods for Spin Systems

David L. Goodwin|arXiv (Cornell University)|Jan 1, 2017

Advanced NMR Techniques and Applications参考文献 273被引用 11

一句话总结

本论文提出了一种基于牛顿型算法与解析计算的二阶方向导数的自旋系统最优控制方法，实现二次收敛，显著加快磁共振应用中的收敛速度。该方法通过正则化确保海森矩阵条件良好，并在固态NMR、电子自旋共振及磁化强度到单重态转移实验中实现高保真度控制，同时大幅减少功能评估次数。

ABSTRACT

Work within this thesis advances optimal control algorithms for application to magnetic resonance systems. Specifically, presenting a quadratically convergent version of the gradient ascent pulse engineering method. The work is formulated in a superoperator representation of the Liouville-von Neumann equation. A Newton-grape method is developed using efficient calculation of analytical second directional derivatives. The method is developed to scale with the same complexity as methods that use only first directional derivatives. Algorithms to ensure a well-conditioned and positive definite matrix of second directional derivatives are used so the sufficient conditions of gradient-based numerical optimisation are met. A number of applications of optimal control in magnetic resonance are investigated: solid-state nuclear magnetic resonance, magnetisation-to-singlet pulses, and electron spin resonance experiments.

研究动机与目标

• 通过利用二阶导数信息，开发一种收敛速度更快的磁共振系统最优控制算法。
• 通过在优化过程中保持海森矩阵条件良好且正定，确保数值稳定性和收敛性。
• 将该方法应用于具有挑战性的自旋系统，包括固态NMR、磁化强度到单重态转移以及ESR实验。
• 通过高效的数值方法，减少控制脉冲优化所需的实验评估次数。
• 为实际实验环境中的应用，实现鲁棒且平滑的控制脉冲。

提出的方法

• 基于保真度泛函对控制场的解析二阶方向导数，开发牛顿-GRAPPA（Newton-GRAPE）方法。
• 采用超算符表示的李维-冯诺依曼方程，以建模受控的开放量子系统。
• 应用乔列斯基分解与信赖域方法对海森矩阵进行正则化，保持其正定性。
• 利用Krylov子空间技术与传播子重用方法，加速矩阵指数与方向导数的计算。
• 通过分段常数近似参数化控制脉冲，并应用范数平方与溢出惩罚项，以满足功率与鲁棒性约束。
• 集成反馈控制回路，采用Nelder-Mead与多方向搜索算法，在实验环境中优化脉冲形状。

实验结果

研究问题

• RQ1 与一阶方法相比，二阶优化方法是否能显著加速自旋系统最优控制中的收敛？
• RQ2 如何高效计算并应用解析二阶方向导数，以在不增加计算复杂度的前提下实现量子控制？
• RQ3 哪些正则化技术可确保在自旋系统中使用基于海森矩阵的优化时的数值稳定性和收敛性？
• RQ4 该方法在磁化强度到单重态转移与固态NMR实验中，能在多大程度上提升控制保真度？
• RQ5 在真实世界反馈控制场景中，面对实验噪声与硬件限制，该方法的性能如何？

主要发现

• Newton-GRAPE方法实现了二次收敛优化，其收敛速度显著优于一阶梯度上升方法。
• 在OOP-ESEEM实验中，仅通过119次实验评估即完成WURST脉冲的优化，达到归一化目标值1.28。
• 在双脉冲Hahn自旋回波实验中，Nelder-Mead算法在341次实验后收敛，最终目标值为0.98，表明信号增强显著。
• 通过乔列斯基分解与信赖域方法进行海森矩阵正则化，确保了正定性与稳定收敛，即使在病态问题中亦然。
• 该方法实现了鲁棒且平滑的磁化强度到单重态转移脉冲，功率与溢出极低，经实验反馈优化验证。
• 传播子重用与Krylov子空间技术降低了计算成本，使该方法的复杂度与基于一阶导数的方法相当。