基于超轻玻色暗物质与自旋相互作用的轴子诱导自旋动力学模拟开源Python软件包axionbloch

《Universe》：axionbloch: An Open-Source Python Package for Simulating Axion-Induced Spin Dynamics

　　

超轻玻色暗物质与自旋的相互作用可被解释为作用于普通物质的赝磁场（pseudomagnetic field）。此类相互作用可利用自旋演化（Bloch）方程建模为常规磁相互作用。本文提出axionbloch，一个用于模拟由常规相互作用和奇异相互作用诱导的自旋动力学的开源Python软件包。数值模拟可作为推导轴子信号特征的工具，这对设计实验搜索和数据分析至关重要。模拟结果已对照理论预期进行校准，确保模拟信号的准确性。axionbloch已在GitHub上公开，允许研究人员在特定的轴子模型和实验配置下模拟赝磁信号。该软件包在Read the Docs上提供文档，并包含应用示例脚本。

研究背景方面，轴子最初被提出用于解决量子色动力学中的强CP问题，此后在各种与自发破缺整体对称性相关的标准模型扩展中提出了更一般的轴子样粒子（ALPs）。从宇宙学意义上看，轴子和ALPs也是有充分动机的暗物质候选者。轴子可能与普通标准模型粒子存在耦合，包括光子、胶子和费米子。关键在于，轴子场梯度可通过赝标量相互作用与费米子自旋耦合。轴子场梯度可被视为模仿真实磁场的赝磁场，由此引入轴子诱导赝磁场的概念。赝磁场诱导自旋动力学，为核磁共振（NMR）和共磁力仪等精密测量实验中的探测提供了途径。已有专门研究给出了轴子信号特征和实验灵敏度，但据研究人员所知，尚无利用数值方法验证所推导轴子信号特征的研究。另一方面，随着更多轴子模型和实验被提出，如宇宙学轴子背景、细粒度轴子和二次耦合轴子，需要一个工具来系统模拟各种物理条件和轴子场配置下的轴子诱导自旋动力学。

研究人员开展了axionbloch软件包的开发工作。该软件包接口用Python 3（3.10或更高版本）实现，底层数值积分用C++实现以提高效率。软件包为轴子场和实验配置建模提供了灵活框架，通过数值求解自旋演化方程，能够预测和分析实验可观测信号。研究得出以下结论：该软件包采用SI单位制，使输入输出可直接与实验参数和观测量比较；在旋转坐标系中采用四阶Runge–Kutta方法对Bloch方程进行数值积分，实现稳定高效的模拟；已进行全面校准以确保模拟准确性；并展示了来自局部维里化轴子暗物质晕的轴子诱导信号模拟示例。该软件包作为开源、有文档和测试的工具，旨在促进精密测量实验中轴子诱导自旋动力学的研究，并为方法开发和新轴子模型的纳入提供平台。此外，该软件包对NMR社区也具有价值，例如可用于教学目的，学生可指定脉冲序列进行模拟并查看可视化模拟结果。这项研究的重要意义在于填补了利用数值方法系统验证轴子信号特征的空白，为各种轴子模型和实验配置下的自旋动力学研究提供了标准化工具。

研究人员为开展研究用到的主要关键技术方法包括：采用Python 3实现上层接口，C++实现底层数值积分，通过pybind11实现Python与C++的接口；利用astropy库处理单位和物理量；采用四阶Runge–Kutta（RK4）方法在旋转坐标系（RCF）中对Bloch方程进行数值积分；利用OpenMP对独立赝磁场实现进行并行化加速；采用pytest框架进行自动化测试；使用Sphinx构建文档。样本方面，校准示例使用热极化质子NMR样品，轴子模拟示例选用129Xe样品，初始极化至50%。

研究结果方面：（1）单位、物理量与常数：软件包利用astropy处理单位和物理量，在axionbloch.constants中定义了轴子和NMR研究中常用的常数如旋磁比，所有常数均为astropy.units.Quantity实例，支持物理量的所有运算。（2）算法：软件包聚焦于自旋-磁共振，以磁化强度作为表征自旋系综的可观测量；在磁场或轴子诱导赝磁场存在下，磁化强度按Bloch方程演化，采用旋转坐标系去除偏置场引起的Larmor进动；横向弛豫时间包括T2弛豫和偏置场不均匀性引起的退相位，通过采样偏置场分布并求解个体Bloch方程模拟不均匀性；采用RK4方法进行数值积分，局部截断误差为O(dt^5)；采用面向对象编程模块化实现，将实验配置、轴子场建模和数值积分分离，类Sample、Magnet和MagField封装实验参数，轴子模型由专门类如MilkyWayAxionHalo实现，Simulation和Simulations类管理单次或多次模拟运行，RK4算法在专用C++库中实现，API为blochSimulation。（3）校准：通过相对L2误差量化模拟精度，激励校准中模拟横向磁化强度与理论预测比较，Rabi振荡示例显示误差随模拟速率增加而减小并在其他因素主导时趋于平台，在100 Hz中等模拟速率下误差一般低于1%；T1和T2弛豫校准同样显示误差随模拟速率增加而减小并趋于浮点精度决定的平台，误差远低于1%。（4）NMR模拟：脉冲NMR模拟中施加一个π/2脉冲和两个π脉冲，磁化强度表现出以Larmor频率振荡、被π/2脉冲翻转到x-y平面、无脉冲时以T2自由衰减、π脉冲后重聚焦（回波）、回波幅度以T2衰减等预期行为；连续波（CW）NMR模拟中，弱场共振时归一化横向磁化强度按正比于B1*t增长，模拟符合预期行为。（5）轴子模拟：MilkyWayAxionHalo类实现银河系暗物质晕轴子模型，轴子场被视为振荡场，频率由轴子质量决定，线宽由维里化暗物质晕的速度弥散决定；模拟的赝磁场和诱导横向磁化强度幅度显示信号在模拟初期因赝磁场而随时间增长，随后因T2弛豫而衰减；通过OpenMP对独立（赝）磁场实现进行并行化加速；示例中总积分步数为10^8，在配备Intel Core Ultra 7 155H CPU的个人计算机上运行时间约为1小时。

讨论部分总结：axionbloch作为模拟磁场或轴子场在各种物理条件下诱导自旋动力学的框架，采用SI单位制使输入输出可直接与实验参数和观测量比较，在旋转坐标系中采用RK4方法对Bloch方程进行数值积分实现稳定高效模拟，已进行全面校准确保准确性，并展示了来自局部维里化轴子暗物质晕的轴子诱导信号模拟示例。作为开源、有文档和测试的工具，axionbloch旨在促进精密测量实验中轴子诱导自旋动力学的研究，并为方法开发和新轴子模型纳入提供平台。该软件包对NMR社区也有价值，例如可用于教学目的。

研究结论翻译：本文提出了axionbloch，一个用于模拟磁场或轴子场在各种物理条件下诱导的自旋动力学的框架。该软件包采用SI单位制，使输入和输出可直接与实验参数和观测量比较。在旋转坐标系中采用四阶Runge–Kutta方法对Bloch方程进行数值积分，实现稳定高效的模拟。已进行全面校准以确保模拟的准确性。本文展示了来自局部维里化轴子暗物质晕的轴子诱导信号模拟示例，以证明该软件包的能力。作为一个开源、有文档和测试的工具，axionbloch旨在促进精密测量实验中轴子诱导自旋动力学的研究，并为方法开发和新轴子模型的纳入提供平台。需要注意的是，该软件包对NMR社区也具有价值，例如可用于教学目的——学生可以指定脉冲序列进行模拟并查看可视化模拟结果。