在向上的磁畴附近做小倾斜
取向上的磁畴,m≈(mx,my,1),mx、my 都很小。把有效场(Abert 2019,第 2 节)中的每一项展开到倾斜的一阶。外场和有效各向异性都沿 z,合起来是 H+HK0,其中 μ0HK0=2K0/Ms=927mT。实验里的链在磁畴中也保留了第 4 章那个小的难轴项,所以 mx 还多感受到 −HKmx。交换作用贡献 (2A/μ0Ms)∇2m。
对平面波 m∝ei(kx−ωt),拉普拉斯算符变成 −k2,交换作用就像一个把倾斜拉回去的附加场 Hex=2Ak2/(μ0Ms)。波越短,磁化弯得越急,回复场也越强。
界面 DMI 在这一阶不出现。它的能量密度是 D(mx∂xmz−mz∂xmx)。因为 mz≈1−(mx2+my2)/2,第一项是倾斜的二阶量;第二项在一阶下是 −D∂xmx,一个全导数,积分为零。这个推导由我们给出,只适用于垂直磁畴中沿 x 传播的波,不涉及其他几何。
完整例题
色散关系和能隙
不计阻尼,LLG 方程化为 m˙x=−γ0H2my 和 m˙y=γ0H1mx,两个刚度场见下式。对第一式再求一次导,代入第二式,得 m¨x=−γ02H1H2mx:倾斜以 ω=γ0H1H2 进动。
k=0、零场时,μ0H1=958mT,μ0H2=927mT,最低频率 f0=ω/2π=26.4GHz。沿磁畴方向加 40 mT,能隙升到 27.5 GHz。低于能隙时,天线只能在身边激起一点很快消失的倾斜。
反过来,从频率求波长。零场、60 GHz 时,由 (γ0H1)(γ0H2)=ω2 解出 Hex:需要 μ0Hex=1.20T,于是 k=0.203nm−1,波长 30.9 nm。10 nm 的波需要 11.5 T 的交换场,频率是 348 GHz。
网格对波做了什么
在格子大小为 a 的链上,二阶差分把 k2 换成 (2−2coska)/a2,短波时偏小。用 1 nm 的格子,60 GHz 的波长是 30.86 nm,连续介质给 30.92 nm,差 0.2%;到 150 GHz 差别增大到 0.6%。这就是第 3 章那个二阶网格误差,只是这里表现为波长的偏移。
实验用链自己的离散色散来检验链:一条 512 格的周期链,加一点小扰动,不计阻尼运行,再对空间和时间做傅里叶变换。每个谱峰都落在公式值半个频率格(1.95 GHz)以内。
完整例题
阻尼让波数变成复数
天线以一个实频率驱动薄膜。有 Gilbert 阻尼时,进动方程里的每个 γ0H 都变成 γ0H−iαω,条件变为 (γ0H1−iαω)(γ0H2−iαω)=ω2。令 X=γ0Hex−iαω,这是关于 X 的实系数二次方程,X 为实数,于是 k2 多出一个正比于 αω 的虚部。
k 的实部给出波长 λ=2π/Rek;虚部给出振幅衰减为 1/e 所走的距离 ℓ=1/Imk。
在 1 nm 的链上,零场、60 GHz:α=0.1 时 λ=30.7nm,ℓ=55.1nm;换成薄膜本身的 α=0.5,λ=28.6nm,ℓ=11.9nm,不到半个波长就衰减掉了。自旋波预设因此把阻尼降到 0.1,这是为了让波看得见而假设的值。
实验中天线在 60 GHz 激发的波,链上测得 30.75 nm,带阻尼的公式给 30.74 nm。基准测试不用闭式解,而用牛顿迭代解同一个条件,在 40、60、100、150 GHz 下与链的波长相差 0.07% 以内,衰减长度相差 0.6% 以内。
这个图像省略了什么
链只把杂散场保留为局域项,这和 Thiaville 等处理超薄膜时用的近似相同。波自身杂散场的长程部分被略去了,本课程没有估算它对色散的影响。天线是 5 nm 宽的均匀场,不是真实微波传输线的模型;也没有温度。