共振的核心：自旋锁定磁共振中的微扰有限时间响应

《Journal of Magnetic Resonance》：A kernel of resonance: Perturbative finite-time responses in spin-lock magnetic resonance

　　

摘要

已经提出使用自旋锁定（Spin-lock）方法通过磁共振来检测微弱的低频磁场，包括神经磁场。目前大多数分析模型仅限于正弦波目标或长时间有效描述，因此它们不适用于重建任意时变波形。我们将自旋锁定检测视为一个有限时间的微扰观测问题。从旋转参考系的布洛赫方程（Bloch equation）出发，我们将弱目标场与参考自旋锁定动态分离，并应用一阶戴森展开（first-order Dyson expansion）。对于未调制自旋锁定，纵向磁化扰动是目标波形的线性函数，等效于一个有限时间的卷积核响应。我们将这种方法扩展到自旋锁定相位调制主动共振激励（Spin-lock Phase-modulated Active Resonance, SPARX）。调制引起的伪磁场被包含在参考动态中，而只有弱目标场被视为扰动。数值布洛赫模拟验证了一阶响应端点对于高斯包络载波、汉宁窗（Hann-window）sinc包络载波以及来自开放式快速不可见频率标记脑电图数据集的波形的有效性。我们还在0到期间，对完整布洛赫解和一阶解进行了比较。对于所有三个磁化分量，振幅扫描可以识别微扰区域和高阶效应的出现。对于SPARX，我们比较了一个基于有限阻尼三角函数的伴随数值核。我们展示了伪场调制如何相对于未调制自旋锁定改变响应带。这种前向模型将自旋锁定脉冲设计与弱磁场波形的时间序列重建联系起来。

引言

在功能神经成像中，直接通过磁共振检测微弱的神经磁场仍然很重要但具有挑战性。自旋锁定方法对此非常有用，因为射频场可以锁定磁化，使其对旋转参考系中的低频扰动敏感。然而，目标场通常不是纯正弦波。生理磁场和由刺激引起的磁场会随时间变化。实际问题在于，如何用有限的几次自旋锁定测量来约束这个时间信号。

这种方法解决了MRI的强空间编码与传统功能MRI的间接时间机制之间的差距。血氧水平依赖性MRI是一种主要的神经成像方法，但其对比度反映了神经活动的血管和代谢效应，而不是电生理电流本身[1]，[2]。自旋锁定方法作为一种非血流动力学替代方案被提出，因为振荡磁场可以在旋转参考系中与磁化发生共振。刺激诱导的旋转饱和（Stimulus-Induced Rotary Saturation, SIRS）概念最早由该方法确立[3]。后续研究在 phantom 和生物系统中检测到了电流或振荡场[4]，[5]，[6]。进一步的发展包括体内自旋锁定成像神经电磁振荡[7]、稳态和超低场方法[8]、为鲁棒性设计的自旋锁定方案[9]，[10]、基于旋转激励的验证[11]，以及应用于非正弦脉冲场[12]。因此，自旋锁定MRI已从基于共振的感测概念发展成为一个用于场检测、鲁棒性和实验验证的脉冲设计家族。

分析研究主要描述了受控振荡目标的磁化动态。封闭形式的布洛赫方程解已被用来解释SIRS类型的动态，包括松弛过程，并减少了重复的数值积分[13]，[14]。低场phantom研究也将理论响应与实验上可检测到的微弱振荡场联系起来[15]。同时，如SPARX这样的相位调制自旋锁定设计在旋转参考系中引入了一个可控的伪磁场，从而控制了参考信号水平[16]。然而，大多数当前的分析方法仍然侧重于正弦波目标。重建不仅需要一个频率的响应，还需要一个前向算子，将准备过程中的任何弱波形映射到准备后的测量磁化。

本研究从数学和数值上推导并验证了这样的有限时间响应模型。我们将布洛赫动态分为参考自旋锁定系统和弱目标场扰动。一阶戴森展开将测量的纵向磁化扰动表示为目标波形的积分和由参考轨迹及传播器定义的核函数。我们首先为未调制自旋锁定准备推导出这种关系。在这种情况下，参考系统的系数是常数，核函数形式简单。然后我们将其扩展到SPARX。其调制产生的伪磁场故意设计得较大，因此属于参考系统而非扰动部分。

主要贡献不在于为单个波形提供另一个布洛赫解或替代数值布洛赫积分。直接积分总是可以预测特定的实验结果。相反，我们提供了一个可解释的线性前向模型，将脉冲设计、核函数形状和端点对比度联系起来。该模型可以通过数值方法评估，比较未调制和相位调制准备方案，并组装成用于时间序列重建的观测矩阵。我们通过完整布洛赫模拟对其进行测试，确定一阶近似何时失效，并为后续的反问题研究确定响应度量方法。

章节片段

理论

本节将自旋锁定检测描述为一个有限时间的微扰观测问题。首先我们将布洛赫方程写成一个仿射线性系统。然后我们将动态分为参考系统和弱目标场扰动。参考轨迹及其均匀传播器定义了响应核函数。我们首先为未调制自旋锁定准备推导出这个核函数，然后将相同的方法扩展到相位调制的SPARX。

方法

本节首先给出所有分析共有的数值条件，然后描述五个验证实验。依次评估了一阶戴森响应对不同持续时间目标波形、所有磁化分量及其振幅依赖性、SPARX核函数、频率响应和一阶精度，以及伪场振幅对响应带的影响。结果和讨论按照相同的顺序进行。

结果

本节按照与方法相同的顺序报告了五个验证实验。每个小节在开头指出了相应的图表。主要数值结果总结在图表内的表格中。文本仅报告观察结果；其解释留待讨论部分。

讨论

在有限持续时间波形验证中，无目标参考轨迹与直接布洛赫积分结果一致，所有三种具有不同特性的波形的端点相对误差都很小。这些结果支持了参考矩阵、传播器和常微分方程（ODE）实现的一致性。它们还表明，当目标场仍处于微扰范围内时，相同的一阶核函数形式适用于不同的有限持续时间波形。核函数的积分是

结论

我们推导并数值验证了一个用于检测弱磁场波形的有限时间微扰响应模型。对于未调制自旋锁定，一阶戴森展开给出了一个核函数积分，将任何弱目标波形映射到端点纵向磁化。当参考动态中包含相位调制引起的伪磁场时，相同的框架也适用于SPARX。该模型将自旋锁定动态与弱场逆向重建联系起来。

手稿准备过程中生成式AI和AI辅助技术的声明

在准备本作品期间，作者使用了OpenAI ChatGPT和Codex来辅助语言编辑、LaTeX起草和代码辅助。使用这些工具后，作者根据需要审阅和编辑了内容，并对手稿的内容负全责。

利益冲突声明

作者声明以下可能被视为潜在利益冲突的财务利益/个人关系：上田弘之（Hiroyuki Ueda）报告获得了日本科学技术机构（Japan Science and Technology Agency）的财务支持；上田弘之还报告获得了日本学术振兴会（Japan Society for the Promotion of Science）的财务支持。如果有其他作者，他们声明没有已知的可能影响工作的财务利益或个人关系。

上田弘之

日本京都大学工程研究生院电气、电子与数字科学系，京都