磁矩、磁场和运动方程
铁磁体里电子自旋叠加成磁化强度 M=Msm。它的大小即饱和磁化强度 Ms,由材料决定;微磁学只追踪它的方向,也就是单位矢量 m。磁场 H 的单位是安培每米,文中用毫特斯拉表示 μ0H。
电子自旋带有角动量,所以磁场转动它的方式,和在蜂蜜里拨动指南针不一样:力矩改变的是角动量,而角动量就是磁矩本身。不计阻尼,得到 Landau-Lifshitz-Gilbert 方程的第一项(Abert 2019,式 117)。旋磁比取 muMAG 标准问题采用的 γ0=2.211×105mA−1s−1。
磁矩为什么绕圈
叉乘同时垂直于两个因子,因此 m 的变化垂直于 m,长度保持为 1;又垂直于 H,两者夹角 θ 保持不变。m 的端点只能沿半径为 sinθ 的圆绕磁场跑,这就是进动。
端点速率是 ∣dm/dt∣=γ0Hsinθ,除以半径得角速度 ω=γ0H,与倾角无关。换算成频率,每特斯拉的 μ0H 对应 28.0 GHz。
阻尼把圆圈变成螺旋
真实磁体会把能量交给晶格和传导电子。Gilbert 加了一项阻碍运动的项 αm×m˙,其中无量纲阻尼 α 很小。把 m˙ 解出来,得到下面的显式形式(Abert 2019,式 118),第二项从 m 指向磁场。
磁矩在磁场中的能量密度是 E=−μ0Msm⋅H。代入显式形式可知它只降不升,速率为 E˙=−αγ0μ0Ms∣m×H∣2/(1+α2)(Abert 2019,式 121 至 124)。α=0 时能量守恒,磁矩永远转下去。
完整例题
停下之前转几圈
让磁场沿 z,用倾角 θ 和方位角 φ 描述 m。显式形式的第一项只改变 φ,速率 ω=γ0H/(1+α2);第二项只改变 θ,θ˙=−αωsinθ。
分离变量:dθ/sinθ=−αωdt。1/sinθ 的积分是 lntan(θ/2),于是 tan(θ/2) 按时间常数 1/(αω) 指数衰减。这段时间里磁矩转了 ω/(2π)×1/(αω)=1/(2πα) 圈,与磁场大小无关。
标准问题的阻尼 α=0.02 对应 7.96 圈。在 25 mT 的磁场中,进动频率 0.700 GHz,时间常数 11.4 ns。单畴实验的第一项检验就是积分这条螺旋,它的四阶积分器复现这个公式,实测收敛阶为 4.00。
实验里的磁体为什么不一样
单畴实验研究的是 muMAG 第 4 号标准问题中的那块纳米磁体:长 500 nm、宽 125 nm、厚 3 nm,Ms=800kA/m,α=0.02,并把它当成一个整体转动的磁矩。这个磁矩感受到的不只是外加的均匀磁场:磁体自己会产生很强的磁场,大小随磁矩指向而变。下一章算出这个场,并解释为什么不加外场,磁矩也会以 4.66 GHz 振荡。
带着问题动手
先预测,再改变,最后解释
- 阻尼从 0.02 提高到 0.1,球面上的螺旋会怎样变?停下之前转的圈数变多还是变少?每一圈会不会变慢?
- 在单畴实验中选“轻轻一推,按 Kittel 频率振荡”,看一遍,再把“Gilbert 阻尼”设为 0.1,再看一遍。
- 用每衰减到 1/e 所转的圈数 1/(2πα),以及进动速率里的因子 1/(1+α2),解释你看到的现象。
单畴实验把 muMAG 第 4 号标准问题的小条当作一个磁矩。畴壁和斯格明子实验使用 Thiaville 等的 Pt/Co 薄膜,杂散场按局域项处理。三个实验都在零温下运行。
试着自己解释
1. 磁场 μ0H=100mT 时进动频率是多少?
展开答案与解释
约 2.80 GHz:每特斯拉 28.0 GHz 乘以 0.1 T。α=0.02 时阻尼因子 1/(1+α2) 的影响不到千分之一。
2. α=0.01 时,tan(θ/2) 缩小到 1/e 期间磁矩转几圈?
展开答案与解释
约 15.9 圈,来自 1/(2πα)。阻尼减半,圈数加倍。
3. 证明不计阻尼时,恒定磁场中的能量 E=−μ0Msm⋅H 不会变化。
展开答案与解释
变化率为 E˙=−μ0MsH⋅m˙=μ0Msγ0H⋅(m×H)。叉乘垂直于每个因子,所以它与 H 的点积为零。