MuMax3

定义

MuMax3是一款开源的、以图形处理器加速的微磁学模拟程序，用于计算纳米至微米尺度磁性材料中磁化强度随时间和空间的变化^[3]。它以有限差分方法离散空间，在二维或三维正交网格上求解微磁学方程，属于微磁学领域常用的数值计算工具之一^[4]。

该软件由根特大学的DyNaMat研究组开发并维护，源代码在GPLv3许可证下发布，可自由修改与再分发^[1][2]。程序以Go语言和CUDA编写，运行需要英伟达显卡及其专有驱动程序，支持Linux、Windows和Mac平台^[3][1]。

原理

MuMax3的核心是求解Landau-Lifshitz-Gilbert方程，即计算约化磁化强度单位矢量m随时间的变化率^[3]。方程中磁矩受到有效场的作用，该有效场由外加磁场、静磁场、海森堡交换场、Dzyaloshinskii-Moriya交换场、磁晶各向异性场和热涨落场等分量叠加而成^[3]。除进动项外，程序还实现了Zhang-Li与Slonczewski两种自旋转移矩项，用于描述电流对磁矩的作用^[1]。

在空间离散上，程序采用正交网格的有限差分格式：磁化强度与有效场等体量定义在网格单元中心，交换耦合等界面量定义在相邻单元的界面上^[3]。为节省显存，材料的空间依赖参数并不逐格存储，而是为每个单元分配0至256之间的区域编号，具体参数存放在以区域编号索引的查找表中；单元之间的交换耦合则以三角矩阵存储，从而支持任意两种材料之间的交换作用^[3]。

静磁场通过对磁化强度与退磁核做离散卷积求得，并借助快速傅里叶变换加速，卷积定理是其计算基础^[3]。程序还提供可选的一维、二维或三维周期性边界条件，通过在模拟盒两侧添加若干镜像来计入重复磁化图像产生的静磁场^[3]。几何建模使用构造实体几何方法，形状由判断点是否在体内的函数表示，可经平移、旋转、缩放并通过布尔运算组合^[3]。

发展历程

MuMax3的设计、验证与性能测试结果由Arne Vansteenkiste、Jonathan Leliaert、Mykola Dvornik、Felipe Garcia-Sanchez和Bartel Van Waeyenberge等人于2014年发表，论文刊登在《AIP Advances》上，文中报告程序的各个组成部分均通过与解析解及微磁学标准问题的比对加以验证^[5][3]。

此后该软件由根特大学DyNaMat研究组持续维护，并被收纳进Debian等发行版的软件仓库，版本号随更新不断推进^[2]。围绕主程序还出现了若干扩展版本，例如面向合成反铁磁体、带有神经形态计算功能的分支，以及用于模拟腔量子电动力学中磁振子与光子相互作用的Mumax3-cQED^[6][7]。

应用

MuMax3广泛用于自旋电子学器件的研究，特别是依赖自旋波、磁畴壁和斯格明子动力学的低功耗自旋器件的设计与优化^[4]。使用者在给定材料参数、各向异性项、几何形状、外加磁场与电流以及自旋转移矩后，即可对微磁系统进行数值求解^[4]。

该程序也用于较专门的材料与结构研究，如软硬磁多层膜的交换耦合与最大磁能积计算、FeGe纳米柱中的磁化反转与自旋织构、自旋阀中的磁液滴孤子，以及Co/Pt铁磁-重金属双层中的斯格明子等^[8]。在教学中，它也常被作为GPU微磁学模拟的入门工具，配合批处理脚本和Jupyter笔记本演示模拟流程^[4]。

局限

MuMax3采用规则网格的有限差分方法，在几何形状的灵活度上不如有限元类微磁学软件，后者可以更好地贴合复杂边界，但通常以牺牲计算性能为代价^[3]。其计算精度受单精度浮点与快速傅里叶变换噪声底限制，例如通过卷积求得的场在约0.01微特斯拉量级以下就不再可靠^[3]。

在物理建模范围上，该软件的验证基于标准微磁学问题，其结果的物理有效性以模型不涉及纠缠效应、也不涉及自旋与声子、光子等准粒子的显著相互作用为前提^[4]。此外，程序运行依赖英伟达显卡，需要计算能力5.0或以上的GPU，因而在硬件选择上受到限制^[1]。

参见

• 微磁学 —— MuMax3所模拟的物理对象与理论基础。

• Landau-Lifshitz-Gilbert方程 —— MuMax3数值求解的核心运动方程。

• 斯格明子 —— 常借助该软件研究的拓扑磁结构。

• 自旋电子学 —— 该软件的主要应用方向之一。

• 图形处理器 —— MuMax3实现加速所依赖的计算硬件。