学习自旋电子学的笔记07：根据微磁学基本能量密度公式推导有效场

阅读说明：由于CSDN自带的Markdown编辑器目前对大篇幅KATEX公式的支持性不太好，导致文章内容有了字数限制，一旦超过字数限制，就不能正常保存和发布了。所以，我将笔记07的部分内容转换为了图片的形式，方便大家阅读！原始的.md文档，里面的公式使用的是KATEX格式，和本文的PDF版本有需要的可以点击下载！

雁来音信无凭，路遥归梦难成。----李煜《清平乐别来春半》

---

########
本文链接：https://blog.csdn.net/qq_43572058/article/details/126234976
CSDN@搬砖工人_0803号
########

前言

在笔记06中为大家推荐了“ubermag”官网上的有关微磁学公式的推导笔记，按照老辈人的说法，微磁学相关公式的推导应该在我们入门后两个月内学会，但并不是每一个人都有着基本的数学素养，有的人面对这些数学式子就是头晕，比如我，，，好在近期闲来无事想着动纸笔比划一下，遇到的一些陌生的名词，陌生的导数，卡住的推导步骤等都会去网上翻一翻，于是总归学会了一点三脚猫的功夫了。

本文的内容就是根据几个微磁学基本能量密度公式来推导其有效场，这几个微磁学基本能量包括：交换能，塞曼能，界面DMI能，体DMI能，单轴各向异性能，立方各向异性能。这几个能量的能量密度公式是已知的，利用有效场的公式：

$$\overrightarrow{H_{eff}}=\frac{1}{-u_0{M}_s}\frac{\delta E_{any}}{\delta \overrightarrow{m}}$$

即可推导出每一项能量对应的有效场。式中：$u_0，M_s$ 分别表示真空磁导率和材料的饱和磁化强度；$\overrightarrow{m}=\overrightarrow{m(\vec{r})}=\frac{\overrightarrow{M(\vec{r})}}{M_s}$ 表示单位磁化强度（矢量场）；$E_{any}$ 表示以上能量的能量密度（标量场）；符号 $\delta$ 表示变分算子。

这里涉及到“变分”这个概念，我自己目前也不是很明白，但从网上查了一些言论描述如下：
设 $F=F(y)，y=y(x)$，有所谓的 函数的函数的导数 表示为：
$$F_y^{\prime }=\frac{\delta F}{\delta y}$$
即
$$\frac{\delta F}{\delta y}=\frac{dF}{dy}或\frac{\partial F}{\partial y}$$

于是，有效场的式子变为了能量密度（标量场）对单位磁化强度（矢量场）求偏导：

$$\overrightarrow{H_{eff}}=\frac{1}{-u_0{M}_s}\frac{\partial E_{any}}{\partial \overrightarrow{m}}$$

大家对这个式子都非常熟悉了，在大部分微磁学文章中，一般这个式子会出现在LLG方程的后面。不过，此处又出现了有难度的步骤，即标量场如何对矢量场求偏导呢？其实结合“梯度”的定义来看：
设有一个标量场 $\varphi =\varphi (x,y,z)$，一个矢量线元 $d\vec{l}$ ，则定义 $\varphi$ 的梯度为：

$$grad\varphi =\nabla \varphi =\frac{d\varphi }{d\vec{l}}=\frac{\partial \varphi }{\partial x}\hat{x}+\frac{\partial \varphi }{\partial y}\hat{y}+\frac{\partial \varphi }{\partial z}\hat{z}$$

梯度是一个标量场对矢量求导，从而不难得到“有效场等于能量密度对单位磁化强度的梯度”这个结论，这句话也经常出现在许多毕业论文第二章的开头，此时我才明白其中的含义。
于是借助梯度的概念将有效场的式子展开：

$$\overrightarrow{H_{eff}}=\frac{1}{-u_0{M}_s}{\nabla }_{\vec{m}}{E}_{any}=\frac{1}{-u_0{M}_s}(\frac{\partial E_{any}}{\partial m_x}\hat{x}+\frac{\partial E_{any}}{\partial m_y}\hat{y}+\frac{\partial E_{any}}{\partial m_z}\hat{z})$$

---

一、塞曼能的有效场

已知塞曼能的能量密度公式：

$$E_{zeeman}=-u_0{M}_s(\vec{m}\cdot \vec{H})$$

其中：$\vec{H}=\overrightarrow{H(\vec{r})}$ 为外加磁场（矢量场），单位 A/m。那么可以得到塞曼能的有效场： $\overrightarrow{H_{zeeman}}$

$=\frac{1}{-u_0{M}_s}\frac{\delta E_{zeeman}}{\delta \overrightarrow{m}}$

$$\begin{gathered} =\frac{1}{-u_0{M}_s}\frac{\partial E_{zeeman}}{\partial \overrightarrow{m}} \\ =\frac{1}{-u_0{M}_s}(\frac{\partial E_{zeeman}}{\partial m_x}\hat{x}+\frac{\partial E_{zeeman}}{\partial m_y}\hat{y}+\frac{\partial E_{zeeman}}{\partial m_z}\hat{z}) \end{gathered}$$

$=\frac{1}{-u_0{M}_s}[\frac{\partial }{\partial m_x}(-u_0{M}_s(\vec{m}\cdot \vec{H}))\hat{x}+\frac{\partial }{\partial m_y}(-u_0{M}_s(\vec{m}\cdot \vec{H}))\hat{y}+\frac{\partial }{\partial m_z}(-u_0{M}_s(\vec{m}\cdot \vec{H}))\hat{z})]$

$=\frac{\partial }{\partial m_x}(\vec{m}\cdot \vec{H})\hat{x}+\frac{\partial }{\partial m_y}(\vec{m}\cdot \vec{H})\hat{y}+\frac{\partial }{\partial m_z}(\vec{m}\cdot \vec{H})\hat{z}$ ******* (式子 1）

其中 $\vec{m}$ 和 $\vec{H}$ 的点乘： $\vec{m}\cdot \vec{H}=m_x{H}_x+m_y{H}_y+m_z{H}_z$

带回式子（1）的第一项：

$$\begin{gathered} \frac{\partial }{\partial m_x}(\vec{m}\cdot \vec{H})\hat{x} \\ =\frac{\partial }{\partial m_x}(m_x{H}_x+m_y{H}_y+m_z{H}_z)\hat{x} \\ =H_x\hat{x} \end{gathered}$$

同理，式子（1）的第二项：

$$\begin{gathered} \frac{\partial }{\partial m_y}(\vec{m}\cdot \vec{H})\hat{y} \\ =\frac{\partial }{\partial m_y}(m_x{H}_x+m_y{H}_y+m_z{H}_z)\hat{y} \\ =H_y\hat{y} \end{gathered}$$

同理，式子（1）的第三项：

$$\begin{gathered} \frac{\partial }{\partial m_z}(\vec{m}\cdot \vec{H})\hat{z} \\ =\frac{\partial }{\partial m_z}(m_x{H}_x+m_y{H}_y+m_z{H}_z)\hat{z} \\ =H_z\hat{z} \end{gathered}$$

最后合并这三项，式子（1）的结果为：

$\frac{\partial }{\partial m_x}(\vec{m}\cdot \vec{H})\hat{x}+\frac{\partial }{\partial m_y}(\vec{m}\cdot \vec{H})\hat{y}+\frac{\partial }{\partial m_z}(\vec{m}\cdot \vec{H})\hat{z}$

$$\begin{gathered} =H_x\hat{x}+H_y\hat{y}+H_z\hat{z} \\ =\vec{H} \end{gathered}$$

这个推导结果表示塞曼能的有效场： $\overrightarrow{H_{zeeman}}$ 就等于外加磁场 $\vec{H}$。

而且由塞曼能的能量密度表达式可知，若已知任意微磁能量的有效场 $\overrightarrow{H_{eff}^{any}}$ 的情况下，那么它对应的能量密度表达式即为：

$$E_{any}=-u_0{M}_s(\vec{m}\cdot \overrightarrow{H_{eff}^{any}})$$

二、交换能的有效场

$$\begin{gathered} =A\frac{\partial }{\partial m_z}[(\frac{\partial m_z}{\partial x}{)}^2+(\frac{\partial m_z}{\partial y}{)}^2+(\frac{\partial m_z}{\partial z}{)}^2] \\ =2A{\nabla }^2{m}_z \end{gathered}$$

将以上三项带回式子（1）合并得到交换能的有效场：$\overrightarrow{H_{ex}}$

$=\frac{1}{-u_0{M}_s}[\frac{\partial }{\partial m_x}(A[(\nabla m_x{)}^2+(\nabla m_y{)}^2+(\nabla m_z{)}^2])\hat{x}+\frac{\partial }{\partial m_y}(A[(\nabla m_x{)}^2+(\nabla m_y{)}^2+(\nabla m_z{)}^2])\hat{y}+\frac{\partial }{\partial m_x}(A[(\nabla m_z{)}^2+(\nabla m_y{)}^2+(\nabla m_z{)}^2])\hat{z}]$

$=\frac{2A}{-u_0{M}_s}({\nabla }^2{m}_x\hat{x}+{\nabla }^2{m}_y\hat{y}+{\nabla }^2{m}_z\hat{z})$

$=\frac{2A}{-u_0{M}_s}{\nabla }^2\vec{m}$

三、界面DMI能的有效场

$=D[0-0+\frac{\partial m_z}{\partial y}-m_z\frac{\partial }{\partial y}]$

$=D[\frac{\partial m_z}{\partial y}-m_z\frac{\partial }{\partial y}]$

其中式子（1）的第三项： $\frac{\partial }{\partial m_z}(D(m_x\frac{\partial m_z}{\partial x}-m_z\frac{\partial m_x}{\partial x}+m_y\frac{\partial m_z}{\partial y}-m_z\frac{\partial m_y}{\partial y}))$

$=D[\frac{\partial }{\partial m_z}(m_x\frac{\partial m_z}{\partial x})-\frac{\partial }{\partial m_z}(m_z\frac{\partial m_x}{\partial x})+\frac{\partial }{\partial m_z}(m_y\frac{\partial m_z}{\partial y})-\frac{\partial }{\partial m_z}(m_z\frac{\partial m_y}{\partial y})$

$=D[m_x\frac{\partial }{\partial x}-\frac{\partial m_x}{\partial x}+m_y\frac{\partial }{\partial y}-\frac{\partial m_y}{\partial y}]$

将以上三项带回式子（1）合并得到界面DMI能的有效场： $\overrightarrow{H_{DMI}^{inter}}$

$=\frac{1}{-u_0{M}_s}[D[\frac{\partial m_z}{\partial x}-m_z\frac{\partial }{\partial x}]\hat{x}+D[\frac{\partial m_z}{\partial y}-m_z\frac{\partial }{\partial y}]\hat{y}+D[m_x\frac{\partial }{\partial x}-\frac{\partial m_x}{\partial x}+m_y\frac{\partial }{\partial y}-\frac{\partial m_y}{\partial y}]\hat{z}]$

$=\frac{D}{-u_0{M}_s}[(\frac{\partial m_z}{\partial x}\hat{x}-m_z\frac{\partial \hat{x}}{\partial x})+(\frac{\partial m_z}{\partial y}\hat{y}-m_z\frac{\partial \hat{y}}{\partial y})+(m_x\frac{\partial \hat{z}}{\partial x}-\frac{\partial m_x}{\partial x}\hat{z}+m_y\frac{\partial \hat{z}}{\partial y}-\frac{\partial m_y}{\partial y}\hat{z})]$

$=\frac{D}{-u_0{M}_s}[(\frac{\partial m_z}{\partial x}\hat{x}-m_z\cdot 0)+(\frac{\partial m_z}{\partial y}\hat{y}-m_z\cdot 0)+(m_x\cdot 0-\frac{\partial m_x}{\partial x}\hat{z}+m_y\cdot 0-\frac{\partial m_y}{\partial y}\hat{z})]$ ********（$\hat{x}，\hat{y}，\hat{z}$ 都是（方向）常矢量，它们的偏导数为0）

$=\frac{D}{-u_0{M}_s}(\frac{\partial m_z}{\partial x}\hat{x}+\frac{\partial m_z}{\partial y}\hat{y}-\frac{\partial m_x}{\partial x}\hat{z}-\frac{\partial m_y}{\partial y}\hat{z})$

$=\frac{D}{u_0{M}_s}[(\nabla \cdot \vec{m})\hat{z}-\nabla m_z]$

四、体DMI能的有效场

已知体DMI能的能量密度公式：

$$E_{DMI}^{bulk}=D[\vec{m}\cdot (\nabla \times \vec{m})]$$

首先展开 $\vec{m}$ 的旋度 $\nabla \times \vec{m}$

$=\left|\begin{array}{ccc}\hat{x}& \hat{y}& \hat{z}\\ \frac{\partial }{\partial x}& \frac{\partial }{\partial y}& \frac{\partial }{\partial z}\\ m_x& m_y& m_z\end{array}\right|$

$=(\frac{\partial m_z}{\partial y}-\frac{\partial m_y}{\partial z},\frac{\partial m_x}{\partial z}-\frac{\partial m_z}{\partial x},\frac{\partial m_y}{\partial x}-\frac{\partial m_x}{\partial y})$

带回原式得到：
$$E_{DMI}^{bulk}=D[\vec{m}\cdot (\nabla \times \vec{m})]=D[(m_x,m_y,m_z)\cdot (\frac{\partial m_z}{\partial y}-\frac{\partial m_y}{\partial z},\frac{\partial m_x}{\partial z}-\frac{\partial m_z}{\partial x},\frac{\partial m_y}{\partial x}-\frac{\partial m_x}{\partial y})]$$

$$=D[m_x(\frac{\partial m_z}{\partial y}-\frac{\partial m_y}{\partial z})+m_y(\frac{\partial m_x}{\partial z}-\frac{\partial m_z}{\partial x})+m_z(\frac{\partial m_y}{\partial x}-\frac{\partial m_x}{\partial y})]$$

$=\frac{D}{-u_0{M}_s}((\frac{\partial m_z}{\partial y}\hat{x}-\frac{\partial m_y}{\partial z}\hat{x}+m_y\frac{\partial \hat{x}}{\partial z}-m_z\frac{\partial \hat{x}}{\partial y})+(-m_x\frac{\partial \hat{y}}{\partial z}+\frac{\partial m_x}{\partial z}\hat{y}-\frac{\partial m_z}{\partial x}\hat{y}+m_z\frac{\partial \hat{y}}{\partial x})+(m_x\frac{\partial \hat{z}}{\partial y}-m_y\frac{\partial \hat{z}}{\partial x}+\frac{\partial m_y}{\partial x}\hat{z}-\frac{\partial m_x}{\partial y}\hat{z}))$ **********（$\hat{x}，\hat{y}，\hat{z}$ 都是（方向）常矢量，它们的偏导数为0）

$=\frac{D}{-u_0{M}_s}((\frac{\partial m_z}{\partial y}\hat{x}-\frac{\partial m_y}{\partial z}\hat{x}+0-0)+(-0+\frac{\partial m_x}{\partial z}\hat{y}-\frac{\partial m_z}{\partial x}\hat{y}+0)+(0-0+\frac{\partial m_y}{\partial x}\hat{z}-\frac{\partial m_x}{\partial y}\hat{z}))$

$=\frac{D}{-u_0{M}_s}[(\frac{\partial m_z}{\partial y}-\frac{\partial m_y}{\partial z})\hat{x}+(\frac{\partial m_x}{\partial z}-\frac{\partial m_z}{\partial x})\hat{y}+(\frac{\partial m_y}{\partial x}\hat{z}-\frac{\partial m_x}{\partial y})\hat{z}]$

$=\frac{D}{-u_0{M}_s}(\nabla \times \vec{m})$

五、磁晶单轴各向异性能的有效场

已知磁晶单轴各向异性能的能量密度公式：

$$E_{uniaxial}=-K_u(\vec{m}\cdot \vec{u}{)}^2或E_{uniaxial}=-K_1(\vec{m}\cdot \vec{u}{)}^2-K_2(\vec{m}\cdot \vec{u}{)}^4（高阶形式）$$

其中 $\vec{u}$ 是磁晶易轴方向（单位矢量）， $K_u$ 是磁晶单轴各向异性常数，单位 $J/m^3$。那么可以得到磁晶单轴各向异性能的有效场： $\overrightarrow{H_{uniaxial}}$

$=\frac{1}{-u_0{M}_s}\frac{\delta E_{uniaxial}}{\delta \overrightarrow{m}}$

$=\frac{1}{-u_0{M}_s}\frac{\partial E_{uniaxial}}{\partial \overrightarrow{m}}$

$=\frac{1}{-u_0{M}_s}(\frac{\partial E_{uniaxial}}{\partial m_x}\hat{x}+\frac{\partial E_{uniaxial}}{\partial m_y}\hat{y}+\frac{\partial E_{uniaxial}}{\partial m_z}\hat{z})$

此处和推导塞曼能的有效场的步骤相同，将 $(\vec{m}\cdot \vec{u})$ 展开为 $m_x{u}_x+m_y{u}_y+m_z{u}_z$ 方便我们求导，带回上式：

$=\frac{1}{-u_0{M}_s}(\frac{\partial }{\partial m_x}[-K_u(m_x{u}_x+m_y{u}_y+m_z{u}_z{)}^2]\hat{x}+\frac{\partial }{\partial m_y}[-K_u(m_x{u}_x+m_y{u}_y+m_z{u}_z{)}^2]\hat{y}+\frac{\partial }{\partial m_z}[-K_u(m_x{u}_x+m_y{u}_y+m_z{u}_z{)}^2]\hat{z})$ ******* (式子 1)

其中式子（1）的第一项： $\frac{\partial }{\partial m_x}[-K_u(m_x{u}_x+m_y{u}_y+m_z{u}_z{)}^2]$

$=-2K_u[(m_x{u}_x+m_y{u}_y+m_z{u}_z)][u_x+0+0]$

$=-2K_u(\vec{m}\cdot \vec{u})u_x$

其中式子（1）的第二项： $\frac{\partial }{\partial m_y}[-K_u(m_x{u}_x+m_y{u}_y+m_z{u}_z{)}^2]$

$=-2K_u[(m_x{u}_x+m_y{u}_y+m_z{u}_z)][0+u_y+0]$

$=-2K_u(\vec{m}\cdot \vec{u})u_y$

其中式子（1）的第三项： $\frac{\partial }{\partial m_z}[-K_u(m_x{u}_x+m_y{u}_y+m_z{u}_z{)}^2]$

$=-2K_u[(m_x{u}_x+m_y{u}_y+m_z{u}_z)][0+0+u_z]$

$=-2K_u(\vec{m}\cdot \vec{u})u_z$

将以上三项带回式子（1）合并得到磁晶单轴各向异性能的有效场： $\overrightarrow{H_{uniaxial}}$

$=\frac{1}{-u_0{M}_s}(-2K_u(\vec{m}\cdot \vec{u})u_x\hat{x}-2K_u(\vec{m}\cdot \vec{u})u_y\hat{y}-2K_u(\vec{m}\cdot \vec{u})u_y\hat{z})$

$=\frac{-2K_u}{-u_0{M}_s}(\vec{m}\cdot \vec{u})(u_x\hat{x}+u_y\hat{y}-u_y\hat{z})$

$=\frac{2K_u}{u_0{M}_s}(\vec{m}\cdot \vec{u})\vec{u}$

六、磁晶立方各向异性能的有效场

已知磁晶立方各向异性能的能量密度公式：

$$E_{cubic}=-K_c[(\vec{m}\cdot \overrightarrow{u_1}{)}^2(\vec{m}\cdot \overrightarrow{u_2}{)}^2+(\vec{m}\cdot \overrightarrow{u_2}{)}^2(\vec{m}\cdot \overrightarrow{u_3}{)}^2+(\vec{m}\cdot \overrightarrow{u_3}{)}^2(\vec{m}\cdot \overrightarrow{u_1}{)}^2]$$

其中 $\overrightarrow{u_1}，\overrightarrow{u_2}，\overrightarrow{u_3}$ 是磁晶易轴方向（单位矢量），而且 $\overrightarrow{u_3}$ 是通过 $\overrightarrow{u_1}，\overrightarrow{u_2}$ 的叉乘确定的：$\overrightarrow{u_3}=\overrightarrow{u_1}\times \overrightarrow{u_2}$ 。$K_c$ 是磁晶立方各向异性常数，单位 $J/m^3$。那么可以得到磁晶立方各向异性能的有效场： $\overrightarrow{H_{cubic}}$

$$\begin{gathered} =\frac{1}{-u_0{M}_s}\frac{\delta E_{cubic}}{\delta \overrightarrow{m}} \\ =\frac{1}{-u_0{M}_s}\frac{\partial E_{cubic}}{\partial \overrightarrow{m}} \end{gathered}$$

$=\frac{1}{-u_0{M}_s}(\frac{\partial E_{cubic}}{\partial m_x}\hat{x}+\frac{\partial E_{cubic}}{\partial m_y}\hat{y}+\frac{\partial E_{cubic}}{\partial m_z}\hat{z})$

$=\frac{1}{-u_0{M}_s}(\frac{\partial }{\partial m_x}[-K_c[(\vec{m}\cdot \overrightarrow{u_1}{)}^2(\vec{m}\cdot \overrightarrow{u_2}{)}^2+(\vec{m}\cdot \overrightarrow{u_2}{)}^2(\vec{m}\cdot \overrightarrow{u_3}{)}^2+(\vec{m}\cdot \overrightarrow{u_3}{)}^2(\vec{m}\cdot \overrightarrow{u_1}{)}^2]]\hat{x}$

$+\frac{\partial }{\partial m_y}[-K_c[(\vec{m}\cdot \overrightarrow{u_1}{)}^2(\vec{m}\cdot \overrightarrow{u_2}{)}^2+(\vec{m}\cdot \overrightarrow{u_2}{)}^2(\vec{m}\cdot \overrightarrow{u_3}{)}^2+(\vec{m}\cdot \overrightarrow{u_3}{)}^2(\vec{m}\cdot \overrightarrow{u_1}{)}^2]]\hat{y}$

$+\frac{\partial }{\partial m_z}[-K_c[(\vec{m}\cdot \overrightarrow{u_1}{)}^2(\vec{m}\cdot \overrightarrow{u_2}{)}^2+(\vec{m}\cdot \overrightarrow{u_2}{)}^2(\vec{m}\cdot \overrightarrow{u_3}{)}^2+(\vec{m}\cdot \overrightarrow{u_3}{)}^2(\vec{m}\cdot \overrightarrow{u_1}{)}^2]]\hat{z})$

$=\frac{1}{-u_0{M}_s}(\frac{\partial }{\partial m_x}[-K_c[[(\vec{m}\cdot \overrightarrow{u_1})(\vec{m}\cdot \overrightarrow{u_2}){]}^2+[(\vec{m}\cdot \overrightarrow{u_2})(\vec{m}\cdot \overrightarrow{u_3}){]}^2+[(\vec{m}\cdot \overrightarrow{u_3})(\vec{m}\cdot \overrightarrow{u_1}){]}^2]]\hat{x}$

$+\frac{\partial }{\partial m_y}[-K_c[[(\vec{m}\cdot \overrightarrow{u_1})(\vec{m}\cdot \overrightarrow{u_2}){]}^2+[(\vec{m}\cdot \overrightarrow{u_2})(\vec{m}\cdot \overrightarrow{u_3}){]}^2+[(\vec{m}\cdot \overrightarrow{u_3})(\vec{m}\cdot \overrightarrow{u_1}){]}^2]]\hat{y}$

---

总结

本文总结了微磁学中通过几个基本的能量密度公式推导出其有效场，内容核心就是求偏导，各种各样的偏导。至于微磁学中的其它公式推导，暂时还没学到皮毛，只有等以后再说了。