DMI 让畴壁变便宜,但能量仍为正
第 4 章给出了单位面积的畴壁能。不计那个小的难轴项,手性正确的 Néel 壁能量为 (Rohart 与 Thiaville 2013,式 10)。这层膜的 ,所以 时畴壁只要每平方米 2.0 mJ。
畴壁能在 处降到零。超过它,畴壁反而降低能量,薄膜会碎成摆线状结构。低于它,均匀磁化才是基态。这时斯格明子是一种能量更高、但稳定的纹理:小的扰动去不掉它,可它也不是能量最低的状态。Rohart 与 Thiaville 提到,这类斯格明子被称为孤立斯格明子。实验把 限制在每平方米 2.8 到 3.4 mJ,低于 。
完整例题
把畴壁首尾接成一圈
把半径为 的斯格明子想成一段弯成圆圈的直畴壁,膜厚为 。畴壁本身的能量是 ;弯曲要额外付出与 成正比的交换能,Rohart 与 Thiaville 估计为 (式 17)。前一项希望圈小,后一项希望圈大。
对 求导并令其为零,得 。代入 ,就是他们的式 18。 时预测 9.47 nm。
这个估计假设圈比畴壁宽得多,。用完整的径向剖面(数值求解他们的式 16a)来检验:每平方米 3.0 mJ 时为 9.40 nm,与估计接近;2.6 mJ 时为 5.77 nm,估计却是 7.42 nm,此时圈只比畴壁宽一点;3.4 mJ 时为 17.3 nm,估计为 15.5 nm。实验用的是完整剖面。
拓扑荷数的是什么
把膜上每一点映射到它的磁矩方向,也就是单位球面上的一点。均匀磁化的膜全部映到北极。斯格明子的磁矩从核心处的南极出发,经过畴壁里赤道上的所有方向,在远处回到北极:它把球面恰好覆盖一次。拓扑荷数的就是覆盖的次数,并按取向带上正负号(Kim 与 Mulkers 2020,式 1)。
对绕方位角转一圈的径向剖面 ,积分化为 。向上膜中的向下核心给出 ,向下膜中的向上核心给出 。
Abert 指出了它的后果:斯格明子无法通过磁化的连续变化变成均匀态。在网格上,这种保护只是近似成立的,本章余下部分讲的就是这一点。
完整例题
在网格上数的两种办法
直接的办法是把式 1 中的导数换成有限差分。只要相邻磁矩差别不大,它就可用,误差随格子尺寸的平方减小。
晶格办法不用导数,Kim 与 Mulkers 取自 Berg 与 Lüscher。把每四个相邻格子中心组成的方格分成两个三角形,再把每个三角形的三个磁矩在球面上张成的带符号立体角加起来。方格有两种分法,Kim 与 Mulkers 对两者取平均,每个三角形权重为二分之一(式 2、3)。
拿沿 、、 的三个磁矩试一下:分子 ,分母为 1,所以 ,,正好是整个球面 的八分之一。把 换成 ,符号翻转为 。
在无边界的膜(周期性边界)中,实验给出的斯格明子晶格拓扑荷为 −1,有限差分拓扑荷在 2 nm 格子下为 −0.956,1 nm 格子下为 −0.989。格子减半,有限差分的误差缩小约四倍,符合二阶方法。Kim 与 Mulkers 的图 1 给出同样的对比:平滑斯格明子为 −0.99978,加入随机无序后为 −0.971,尽管覆盖方式并没有变。
开放边缘让拓扑荷偏离整数,这并不是出错
有 DMI 时,自由边缘上的磁矩无法保持垂直。Rohart 与 Thiaville 推出了边界条件(式 4b、6),并给出直边处的倾斜 ,(式 8)。 时 ,边缘磁矩倾斜到 ,向膜内约经过 的距离消失。
倾斜的边缘磁矩覆盖了球面的一部分,所以有开放边缘时,两种定义都不必给出整数。在斯格明子实验 60 nm 宽的赛道里,斯格明子居中时晶格拓扑荷为 −0.999935,离下边缘 23.6 nm 时为 −0.999631;换成 1 nm 格子,分别是 −0.999937 和 −0.999636。格子减半,这点亏缺不变,说明它来自边缘而不是网格。没有斯格明子的赛道两种定义都给 0,因为直边沿整个长度倾斜方式都一样。
Kim 与 Mulkers 在一个边缘弯曲的 200 nm 圆盘里看到更大的偏离:零温下有限差分为 −0.9305,晶格方法为 −0.9468,他们把它归因于边缘倾斜。他们的结论是:在带手性边缘的受限结构里,非整数的拓扑荷要谨慎解读。本课程把这一点定为规则:靠近边缘时拓扑荷偏离整数,首先说明的是边缘;其中是否有一部分属于数值误差,只有网格测试才能判断。