--- title: "怎样读模拟给出的斯格明子拓扑荷和钉扎阈值" description: "逐图读 Kim 与 Mulkers 的论文（arXiv:2006.13336v2），再用同样的检查去读本课程实验关于拓扑荷、穿边消失和脱钉的记录。" canonical_url: "https://sandboxphysics.com/zh/stories/micromagnetics/reading-skyrmion-charge-and-pinning-results/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/micromagnetics/reading-skyrmion-charge-and-pinning-results/" last_modified: "2026-10-29" --- [微磁学 · 边学边做实验](https://sandboxphysics.com/zh/stories/micromagnetics/) L4 · 第 8 章 # 怎样读模拟给出的斯格明子拓扑荷和钉扎阈值 把 arXiv:2006.13336v2 打开放在旁边：J.-V. Kim 与 J. Mulkers，《On quantifying the topological charge in micromagnetics using a lattice-based approach》。论文用三个例子，区分拓扑荷的数值误差和真实变化。三个例子分别是：有温度而无边界的膜、产生斯格明子对的赛道、带边缘的圆盘。读完论文，再把同样的检查用到斯格明子实验的记录上；最后读钉扎阈值，这一步没有论文可作参照。 Sandbox Physics · 2026-10-29 这一章 [查看八章课程](https://sandboxphysics.com/zh/stories/micromagnetics/) 开始前 第 6、7 章；或者已经知道拓扑荷的两种定义和 Thiele 方程。 ## 式 1 与图 1：用导数算拓扑荷 式 1 就是第 6 章的连续拓扑荷，导数用格子上的有限差分计算。图 1(a) 是零温下的斯格明子，颜色表示 $m_z$；(b) 把每个格子映射到单位球面上，线段表示相邻格子之间的连接。这张连接网恰好把球面覆盖一次，式 1 给出 −0.99978290。 (c)、(d) 给每个磁矩加上随机无序。球面上的网变形了，但仍然覆盖球面一次，所以拓扑荷并没有变；式 1 却给出 −0.97115153。这张图说明：相邻磁矩差别大时，有限差分值会失准，即使纹理本身没变。 常见的误读是把 −0.971 当成斯格明子已经部分散开。(b) 和 (d) 的拓扑结构相同，变差的只是导数的估计。 ## 式 2、式 3 与图 2：用三角形来数 式 2、3 用相邻三个磁矩组成的三角形的带符号立体角来数拓扑荷，这是 Berg 与 Lüscher 的构造。图 2(a)、(b) 是把四个格子组成的方格分成两个三角形的两种分法；(c) 对两种分法取平均，每个三角形权重二分之一；(d) 处理有限形状：叉号表示空格，在弯曲边缘处只有一种分法放得下，那个孤立的三角形权重取 1。 这说明：晶格值不需要导数，相邻磁矩差别大时也不会变差。要避免的误读是认为它处处精确。只有三角形铺满一个封闭曲面时，比如周期性的膜，平滑纹理的总和才是整数；有开放边缘时不必是。 $$ \tan\frac{q_{ijk}}{2}=\frac{\mathbf m_i\cdot(\mathbf m_j\times\mathbf m_k)}{1+\mathbf m_i\cdot\mathbf m_j+\mathbf m_i\cdot\mathbf m_k+\mathbf m_j\cdot\mathbf m_k},\qquad Q=\frac{1}{4\pi}\sum_{\langle ijk\rangle}q_{ijk} $$ ## 图 3：100 到 400 K 下的无边界膜 第 3.1 节给出运行设置：边长 200 nm、厚 0.6 nm 的方膜，256 乘 256 个格子，面内周期性边界，$A=16\,\mathrm{pJ/m}$，$M_s=1.1\,\mathrm{MA/m}$，$D=2.7\,\mathrm{mJ/m^2}$，$\alpha=0.3$，不计偶极相互作用。拓扑荷每 5 ps 采样一次，共 100 ns。图 3 每个子图对应一个温度，插图是有限差分值的直方图。 有限差分曲线在起伏，温度越高散得越开，直方图的峰也不在 −1。晶格曲线几乎是平的；作者把残余的抖动（400 K 时从 −1.0000008 到 −0.9999996）归因于单精度运算。400 K 时晶格曲线还出现短暂的 −1/2 和 −1 状态，作者解读为热驱动的半子和斯格明子的产生与湮灭。 这说明：在无边界的膜里，晶格拓扑荷只在真有事情发生时才变化。误读是把有限差分曲线对时间平均，当作拓扑荷报告；这个平均值偏离 −1。 ## 图 4：四种网格下的赛道 第 3.2 节用一条 1000 乘 125 nm 的赛道，$D=0.1\,\mathrm{mJ/m^2}$，$\alpha=0.05$。一个 25 nm、面内各向异性的圆盘作为成核区，$25\,\mathrm{TA/m^2}$ 的电流（不含非绝热力矩）从那里不断甩出斯格明子与反斯格明子对。(a) 到 (d) 的格子数从 1024 乘 128 依次减到 384 乘 48。 最细的网格上，两种定义除了产生和湮灭时的少数采样点外都一致。格子变大后，有限差分曲线出现更多非整数值，进而是平滑的漂移，而晶格曲线在每个子图里都保持接近整数的平台。这说明：粗网格上平滑的非整数漂移来自导数。 要避免的误读：四个子图里的事件序列并不相同。作者解释说，圆形成核区在不同网格下被切成格子的方式不同。不同网格下事件不同，来自离散后的几何；两种拓扑荷在各自的网格上都算对了。 ## 第 3.3 节与图 5：有边缘的圆盘 现在边缘起作用了。在同样 256 乘 256 格子的 200 nm 圆盘里，零温下的斯格明子有限差分值为 −0.9304853，晶格值为 −0.9468352。作者把偏离 −1 的部分归因于手性边界条件造成的边缘磁矩倾斜，也就是 Rohart 与 Thiaville 推出的那个倾斜（第 6 章）。 图 5 中，有温度时两条曲线都在起伏，晶格曲线小一些，作者把这些起伏归因于边缘磁矩。300 K 和 400 K 时，滑动平均（1 ns 窗口，黑线）陡然降到零：斯格明子在边界处被消灭了。这个过程中平均曲线的小幅摆动，来自边缘处部分反转的区域。 这说明：在受限结构里，非整数拓扑荷是预期之中的，拓扑荷也可以从边缘流走。误读是把 −0.9468 当作数值误差，想靠加密网格消掉它。 完整例题 ## 把这些检查用到斯格明子实验的记录上 斯格明子实验可以把记录导出为 JSON，基准表给出了网格对比。读每个拓扑荷时只问一个问题：格子减半，它动不动？ 无边界膜里晶格拓扑荷为 −1，有限差分值在 2 nm 格子下为 −0.956，1 nm 格子下为 −0.989。格子减半，误差约缩小四倍，符合二阶差分，所以这点亏缺是数值误差。 60 nm 宽的赛道里，斯格明子居中时晶格拓扑荷为 −0.999935，离下边缘 23.6 nm 时为 −0.999631；1 nm 格子下分别为 −0.999937 和 −0.999636。格子减半，亏缺不动，所以它属于边缘。没有斯格明子的赛道，两种定义都给 0。 然后是最容易误读的情形：拓扑荷在边界处发生变化。$\beta=0$、$u=150\,\mathrm{m/s}$ 时，斯格明子被推出下边缘，2 nm 格子下 0.64 ns 后消失，1 nm 格子下 0.63 ns。其间每 10 ps 一次的采样中，分别有 16 个和 17 个落在 −0.9 与 −0.1 之间，相邻采样最大变化 0.23，最终拓扑荷与零相差不到 $10^{-7}$。拓扑荷是从边缘平滑地流出去的，两种网格下时间几乎相同。在周期性的膜里，平滑纹理的晶格拓扑荷是整数，只能跳变。拓扑荷经边缘平滑流出，又不随网格改变，记录下的就是一个物理事件：斯格明子离开了赛道。 ## 怎样读钉扎阈值 没有哪篇论文给出这些阈值的参考值，所以要从运行判据去读。第一个纳秒之后，如果中心在 0.5 ns 内移动慢于 0.2 m/s，记为被钉住；中心沿赛道越过 95 nm 或斯格明子消失，记为挣脱；15 ns 后仍未判定，记为未定。二分法在 2 nm 网格上从 30 到 50 m/s 走五步；1 nm 网格只在区间下方和上方各 8% 处各跑一次。 弱点在 41.875 m/s 时钉住斯格明子，42.5 m/s 时放走它；$\beta=0$ 时，边缘在 38.75 m/s 时挡住它，39.375 m/s 时挡不住。两者都是区间：0.625 m/s 的宽度是二分法的步长，不是误差棒。1 nm 的运行只说明网格让每个阈值移动不到 8%。阈值还取决于起点和运行历史，影响约 1 m/s。 Iwasaki 等说明了每个阈值由哪个力决定：自由斯格明子为 $|F_{\mathrm{pin}}|/|G|$，通道中为 $|F^{\mathrm{pin}}_x|/(\beta\mathcal D)$。这告诉你哪些参数重要，却给不出这条赛道的数值。误读是把 42 m/s 当作这层膜的性质；它是零温下、一个弱点、一个起点、按这套判据得到的结果。 ## 读下一个结果时的清单 读拓扑荷：用的是哪种定义；边界是周期性的还是开放的，DMI 是否让边缘倾斜；格子减半后怎样；温度、平均窗口和浮点精度。读阈值：被钉住和挣脱的判据、时间上限、区间宽度、网格检验、起点，以及 Thiele 方程说哪些力必须平衡。 带着问题动手 ## 先预测，再改变，最后解释 1. 运行之前，预测强电流把斯格明子推出边缘时，两条拓扑荷曲线会怎样离开 −1。 2. 在斯格明子实验中选“被推出边缘而消失”，看两条拓扑荷曲线，再导出记录，数一数落在 −0.9 与 −0.1 之间的采样点。 3. 解释：这里拓扑荷从 −1 平滑地变到 0，为什么不是数值误差？哪一次运行能证实这一点？ [沿赛道推动一个斯格明子](https://sandboxphysics.com/zh/labs/skyrmion-track-pinning/) 单畴实验把 muMAG 第 4 号标准问题的小条当作一个磁矩。畴壁和斯格明子实验使用 Thiaville 等的 Pt/Co 薄膜，杂散场按局域项处理。三个实验都在零温下运行。 ## 试着自己解释 ### 1. Kim 与 Mulkers 的图 4 中，四种网格上产生斯格明子对的时刻各不相同。是不是有一种拓扑荷定义算错了？ 展开答案与解释 不是。作者解释说，圆形成核区在每种网格下被切成格子的方式不同，所以事件不同。应当比较每条曲线的平台，而不是跨网格比较事件的时刻。 ### 2. 一个圆盘模拟在零温下给出晶格拓扑荷 −0.947。判定它是误差之前，你会检查什么？ 展开答案与解释 把格子减半，看它动不动；把边缘倾斜与 $\sin\theta_{\mathrm{edge}}=\Delta/\xi$ 比较；在周期性膜里算同一个斯格明子，那里晶格拓扑荷应为 −1。Kim 与 Mulkers 在圆盘里得到 −0.9468，并归因于边缘倾斜。 ### 3. 2 nm 网格上脱钉区间为 41.875 到 42.5 m/s，1 nm 网格在 38.8 m/s 被钉住、45.6 m/s 挣脱。你能说什么？ 展开答案与解释 在 2 nm 网格上，阈值落在这个区间内；在 1 nm 网格上，只知道它在 38.8 与 45.6 m/s 之间，所以网格让它移动不到 8%。没有独立的参考值，结论只适用于零温、这个弱点和这个起点。 ## 可以放在旁边读的来源 1. [Rohart & Thiaville, Skyrmion confinement in ultrathin film nanostructures in the presence of Dzyaloshinskii-Moriya interaction, Physical Review B 88, 184422 (2013), arXiv:1310.0666](https://arxiv.org/abs/1310.0666) 已读全文：Pt/Co/AlOx 参数和超薄膜中杂散场的局域处理，自由边缘的边界条件与边缘倾斜（式 3 至 8，以及图 1 中含完整偶极作用的对照），孤立斯格明子的方程（式 15、16）和近似半径（式 18）。 2. [Iwasaki, Mochizuki & Nagaosa, Current-induced skyrmion dynamics in constricted geometries, Nature Nanotechnology 8, 742 (2013), arXiv:1310.1655](https://arxiv.org/abs/1310.1655) 读了 arXiv 版全文：自旋转移力矩驱动下斯格明子的 Thiele 方程（式 4）、贴边滑行速度（式 5、6）、自由平面和窄通道中的钉扎阈值、边缘排斥以及在边缘湮灭（图 5）。文中模拟用的是晶格自旋模型，不是本课程的 Pt/Co 薄膜。期刊版没有打开。 3. [Kim & Mulkers, On quantifying the topological charge in micromagnetics using a lattice-based approach, IOP SciNotes 1, 025211 (2020), arXiv:2006.13336](https://arxiv.org/abs/2006.13336) 已读全文：有限差分拓扑荷（式 1）、由带号立体角求的晶格拓扑荷（式 2、3，图 2）、第 3 节的三个算例（图 3 至图 5），以及边缘倾斜使两种取值都不是整数的圆盘算例。 [上一章 电流把斯格明子斜着推走，边缘和弱点决定它停在哪里](https://sandboxphysics.com/zh/stories/micromagnetics/current-driven-skyrmions/) [回到课程 八章课程与全部实验](https://sandboxphysics.com/zh/stories/micromagnetics/)