Numerical solution of Bloch‘s equation for neutron spin precession

1 相关物理

1.1极化向量 P

p是一个三维向量，$P_x,P_y,P_z$表示中子在三个方向上的极化率

1.2布洛赫方程

程序中用数值的方法模拟了布洛赫方程

布洛赫方程比较难计算出数值解，本文（程序）的主要工作就是找到另一种方法，模拟布洛赫方程
将B与P的关系模拟出来

1.3磁通量的积分 S

1.4旋磁比

$$\gamma =\frac{\omega }{B_0}$$
本项目内容：
在研究极化中子成像技术和装置的基础上，设计与实现一个对应的中 子成像系统，采用高性能优化算法（计算过程必须快速、准确和精确），通过跟踪中子传输过程中的中子自旋，对大量的中子自旋历史进行平均，解决布洛赫的 进动方程的任意磁感应配置问题，包括时间依赖性等，并通过与已知解析表达式 的简单情况的比较来评估精度，给出非均匀磁场情况下的精度和执行时间

2 方法

研究极化向量P与磁场$B(t)$的关系，但又要避开布洛赫方程。

1. 首先根据距离切片，
2. 切片内按$B$的方向旋转坐标系，使$P$和$z$轴夹角为拉莫尔进动角$\theta$，便于建立微分方程
3. 引入磁通量积分$S$，问题变为：切片内$P$与$S(t)$的关系
4. P可用$\theta ,\varphi$表示，给出了$\theta ,\varphi$关于$S(t)$的微分方程
5. 用龙格库塔方法来求解微分方程的数值解，用龙贝格积分计算$S(t)$

3 过程

1. 将磁场切片，近似认为切片内磁场为稳恒磁场
   
   其中，以第一个切片的磁场方向为基准坐标系

2. 在切片内，根据切片内的磁场方向，旋转坐标系，为每个切片确定旋转矩阵R
   

   旋转的过程对应的旋转矩阵
   
   旋转坐标系后，将P分解
   根据布洛赫方程，P的方向会受B影响（P大小不变），所以P与B之间的夹角$\theta$也会随着B的变换而变换。P除了与切片内磁场相关，还与轨迹相关

3. 所以$S(t)$可以作为P的自变量，跟踪中子自旋的过程变为，跟踪P如何随着S变换

   因为：$P_x^2+P_y^2+P_z^2=1$ ，P的自由度为2
   即：用$\varphi ,\theta$就可以将P向量确定下来
   $P_z=cos\theta$
   $P_x=\sqrt{1-P{b}_z^2}cos\varphi$
   $P_y=\sqrt{1-P{b}_z^2}sin\varphi$

4. 根据两个公式求解微分方程$\frac{dcos\theta }{dS},\frac{d\varphi }{dS}$的数值解，进而得出P

5. 数学
   由上面的过程知，求p的微分方程转为了求$\varphi ,\theta$的微分方程
   论文中给出了公式【18】【20】
   
   
   

通过求解微分方程，最终得到$\varphi ,\theta$的数值解

可以解出$P_x、P_y、P_z$，最后每个切片内再乘以$R^{-1}$来还原坐标系后就得到每个切片内的极化向量P了

5 程序