Mumax3 微磁学模拟这事我估计劝退了不少人。软件本身免费开源、GPU 算得快、支持脚本批量跑但问题是——你打开一个官方示例脚本第一眼看到的是一堆SetGridSize、Msat 800e3、relax()这样的代码旁边连个像样的注释都没有新手上手基本靠猜。我当年也是从改别人的代码模板开始的踩了无数坑之后才慢慢总结出一套从基础设置到材料参数优化都顺手的写法。这篇就把我常用的 Mumax3 代码模板拆开讲讲从第一行代码怎么搭到参数怎么批量扫描再到跑挂了怎么排查一次性说清。适合刚接触微磁模拟的科研小白也适合给课题组搭公共模板的老手参考。1. Mumax3代码模板的整体设计与核心思路1.1 为什么Mumax3适合“模板化”工作流Mumax3 这类脚本式模拟软件和传统的 GUI 界面软件有个本质区别你的输入不是一串鼠标点击而是一个可以复用的文本文件。这意味着什么意味着你第一次调好的那套几何、材料、外场和输出配置可以原封不动作为下一组实验的起点。这种复用性在科研里尤其值钱因为你做的绝大多数模拟都是从“上次那个能跑的脚本”改出来的。我见过不少实验室新人来了先抄师兄的脚本但抄过来的脚本往往带着一堆冗余设置甚至还有隐藏 bug。真正高效的做法是把一份经过验证的脚本沉淀成代码模板参数区集中在开头输出设置统一注释规范。这样做的好处很直接——换材料、改尺寸、调参数只需要动模板里那十几个数值不用每次从零开始猜测语法。代码模板的另一个价值是“可追溯”。科研讲究可重复性如果你三个月后要复现去年的一组数据靠脑子记参数是不可能的但如果你当时基于模板跑的模板本身加上运行记录基本能完整还原所有条件。这也是我强烈建议用模板而不是随手写脚本的原因。1.2 一份mx3脚本的标准骨架Mumax3 脚本用的是它自己的 mx3 语言语法上继承了 Go 的简洁风格。一份完整脚本按功能可以拆成五个区块网格与几何、材料参数、初始状态、演化与求解、输出监控。我第一次写出这份骨架的时候感觉像突然拿到了地图。// 区块1网格与几何 Nx : 256 Ny : 128 Nz : 1 SetGridSize(Nx, Ny, Nz) SetCellSize(2e-9, 2e-9, 2e-9) // 区块2材料参数 Msat 800e3 Aex 1.3e-11 alpha 0.02 // 区块3初始状态 m Uniform(1, 0, 0) // 区块4外场与求解 B_ext vector(0, 0, 0) relax() // 区块5输出 tableadd(M) tableautosave(1e-12) save(m)区块顺序不能乱。网格没定义后面所有参数都没有物理载体材料参数没设初始化磁化结构就无从谈起输出不配置跑完了拿不到数据等于白算。理解了这五个区块你看任何官方示例脚本都不会再眼花缭乱。有一点特别提醒mx3 里:和是有区别的。:用来定义局部变量才是给模拟器内置参数赋值。初学经常把Msat 800e3写成Msat : 800e3结果脚本直接报错。这种细节模板里最好用注释标出来。1.3 做模板前先问自己三个问题很多新手拿到模板就改参数改完就跑跑完发现结果不对再回头检查浪费大量算力。我建议动手写模板前先花十分钟过三个问题。第一研究对象是什么几何薄膜、纳米线、纳米圆盘还是三维颗粒几何决定了网格划分方式和是否需要设置周期边界。第二关心的物理量是什么你关心的是静态平衡态比如磁畴结构还是动态演化比如磁化反转过程这决定了脚本里是直接用relax()收敛还是用run(t)跑一段时间。第三特征尺寸是多少特征尺寸直接决定网格大小上限网格太粗会丢物理太细会浪费 GPU 显存。这三个问题想清楚了模板才能写得有的放矢。我见过太多人把二维薄膜的模板硬套到三维结构上结果不是算得慢就是结果完全背离物理。2. 基础设置从空脚本到第一个能跑的仿真2.1 网格与几何第一行代码决定一切SetGridSize和SetCellSize是我认为整份脚本里最关键的两行。前者决定把仿真区域切成多少个格子后者决定每个格子代表的实际尺寸。二者相乘就是你的仿真区域物理大小。这里有个概念必须时刻记住交换长度。它是衡量磁化方向能发生多大空间尺度变化的物理量公式是lex sqrt(2 * Aex / (mu0 * Msat^2))其中mu0 4 * pi * 1e-7。以坡莫合金为例Msat 800e3 A/mAex 1.3e-11 J/m算出来lex大约 5.7 nm。网格尺寸必须小于这个值一般取 2 到 3 nm 才稳妥。如果网格比交换长度还大模拟会丢掉短波长的磁结构涡旋、斯格明子这类现象直接算不出来。实际写代码时我通常先设定物理尺寸再反推网格数。比如要做 500 nm × 250 nm 的矩形薄膜网格尺寸取 2 nm那么SetGridSize(250, 125, 1)SetCellSize(2e-9, 2e-9, 2e-9)。注意网格数尽量取偶数某些傅里叶变换操作在偶数网格上更稳。2.2 材料参数的“标准填法”材料参数区是模板的“灵魂”也是隐性 bug 最多的区域。我整理了一张常用参数表照着填一般不会错参数名物理含义单位常见值示例Msat饱和磁化强度A/m800e3坡莫合金Aex交换刚度常数J/m1.3e-11坡莫合金alpha吉尔伯特阻尼系数无量纲0.01~0.1Ku1单轴各向异性常数J/m³应材料而定anisU各向异性轴方向无量纲矢量vector(0, 0, 1)Dind界面DMI强度J/m²1e-3重金属/铁磁界面Dbulk体DMI强度J/m³1e-3B20结构特别强调单位。Mumax3 内部统一用 SI 单位制长度是米磁化强度是安培每米交换刚度是焦耳每米。比如纳米尺寸要写成2e-9而不是2。饱和磁化强度要写成800e3而不是800。我见过一个项目Aex 多写了一个数量级结果模拟出来的磁畴尺寸明显偏大查了一整天才发现是单位换算错了。各向异性能。如果你的材料是垂直各向异性的比如 CoFeB 薄膜通常要同时设置Ku1和anisU后者指向易轴方向。这个方向要和各向异性定义一致否则算出来的矫顽力完全不靠谱。磁性参数不太确定时可以先去查文献里同类材料的参数不要自己拍脑袋。2.3 初始磁化与外场让模拟贴近物理初始磁化状态的选择直接影响模拟能否收敛到正确的平衡结构。最基础的是Uniform(dx, dy, dz)把所有格点的磁化设成同一个方向更贴近真实的是RandomMag()为每个格点生成随机方向适合模拟从高温冷却下来的无序磁化过程。对初学者我的建议是如果只是算静态平衡态用Uniform作为初态就够如果一个初态算出来的结果只有局部能量最小可以尝试多种初态对比均匀 x、均匀 -x、随机看看最终是否收敛到同一个结构。这个办法我用得最多能有效避开“选错初态导致结果失真”的问题。外场通过B_ext设置单位是特斯拉。比如沿 x 轴加 50 mT 场写作B_ext vector(50e-3, 0, 0)。如果要跑磁滞回线通常会写一个循环逐步改变外场强度Mumax3 里既可以在脚本里用for循环也可以由外部脚本多次调用。relax()是快速阻尼收敛到局部能量最小的操作适合求静态结构。run(t)则是真实的时间演化适合看磁矩动态过程。二者的区别在于时间复杂度——relax()相当于把阻尼设得很大迅速收敛run()则按 LLG 方程逐步推进。搞错这两者你会得到看起来“正常”但其实完全错误的动态结果。2.4 输出设置后处理从文件类型开始很多人轻视输出设置跑完发现数据文件不是自己想要的格式再补跑一遍时间全浪费了。我通常在模板里固定加这几条tableadd(M) tableadd(E) tableautosave(1e-12) save(m)tableadd(M)是把总磁化矢量记录到表格tableadd(E)记录总能量tableautosave设置自动保存间隔。运行结束后目录下会生成一个以脚本名命名的.txt表格文件里面每一列对应一个记录量方便后续画曲线。磁化结构的空间分布要保存成 OVF 格式。默认情况下save(m)输出的是文本 OVF体积比较大。如果你存储空间吃紧可以转换成二进制 OVF体积小很多读取也快。最终可视化可以用开源的 ParaView 或 Python 的 ubermag 工具链OVF 文件在这两个工具里都能直接打开。有一点要提前规划OVF 文件是“快照”默认每一帧覆盖旧的还是追加新的取决于文件名是否带时间戳。做长时间演化时如果不加控制一个仿真可能产生几十 GB 的文件非常容易把磁盘塞爆。模板里我就养成了习惯保存时给文件名加时间戳或索引比如save(m, m_ str(i))方便后期回溯。3. 材料参数优化的实操方法3.1 先搞懂哪些参数值得扫材料参数优化不是把所有参数都扔进去跑一遍网格搜索那既不经济也不科学。我习惯先按“对模拟结果的影响方式”把参数分成静态主导和动态主导两类。静态主导的参数包括Msat、Aex、Ku1、DMI。这些参数直接决定平衡态的结构和能量景观。比如你要调饱和磁化强度看对矫顽力的影响Msat就是核心变量你要研究斯格明子稳定性DMI 和Aex的比值是核心。动态主导的参数主要是alpha阻尼大小影响磁矩进动衰减速度和驰豫时间但对静态稳定结构基本没影响。做优化前先明确目标量。比如目标是提高矫顽力Hc那你追踪的是M-H曲线中的过零场点目标是稳定斯格明子你追踪的是斯格明子半径随外场的变化。目标不明确扫描再多的参数也提炼不出有用结论。我在这一块吃过亏早期什么都想扫结果数据一堆能写进论文的没几个。3.2 用Mumax3内置功能做参数扫描Mumax3 的脚本语言本身支持循环和条件判断所以最简单的参数扫描方式就是直接在脚本里写一个for循环。for i : 0; i 10; i { B_ext vector(0, 0, -(float32(i) * 0.01 0.05)) run(1e-10) save(m, loop_ str(i)) }这个脚本会从 50 mT 开始步进 10 mT逐步变到 150 mT每次都运行一段真实时间并保存快照。看起来很直接但有个隐藏坑如果要求每一步都达到平衡态直接循环可能不够准因为步长太大会导致磁化来不及响应步长太小又会浪费大量时间。我的经验是先用大步长跑一遍找出磁化发生突变的场区再在那个区间加密步长效率和准确率都能兼顾。新版 Mumax3 里还提供了一些更高级的参数扫描机制但不同版本语法兼容性不太一样。如果你用的版本较新建议直接看官方文档里 “Sweep” 相关章节。不过说实在的我日常批量调参更多用外部脚本控制原因在下一节讲。3.3 更通用的办法批量脚本加参数模板标签如果你的参数扫描量级在几十到上百组Mumax3 脚本里的循环会变得笨重单次跑太久、中途崩了不好续、之后也不容易单独分析某一组的结果。我更推荐的做法是写一个通用模板脚本用固定格式的标签标记参数位置然后让 Shell 或 Python 脚本批量替换标签并逐个调用 mumax3。在 mx3 模板里可以这样标记// PARAM:MSAT Msat 800e3 // PARAM:AEX Aex 1.3e-11然后写一个简单的 Shell 脚本用sed替换后再运行for msat in 600e3 700e3 800e3 900e3; do sed s/PARAM:MSAT/Msat $msat/ template.mx3 run_$msat.mx3 mumax3 -o out_$msat run_$msat.mx3 done这样每跑完一组参数结果都在独立目录里不会互相干扰。等全部跑完再用一个后处理脚本汇总所有目录的角度曲线或能量曲线。我实验室的模拟基本都是这种流程维护起来比一个超长脚本清晰得多。更进阶一点可以用 Python 生成参数组合并直接调用命令行。Python 在科学计算群体里普及度高用subprocess调mumax3也不复杂。还有个附带好处你可以把参数组合、运行时间、输出文件路径都记录在一个 CSV 里以后写论文方法部分时数据集怎么来的一查就有。3.4 结果判读从Table数值到OVF可视化跑完不是终点会读数据才算真正做完。Table 文件里的M是总磁化矢量单位通常是 A/m 或者相对值具体看你的配置E是总能量单位是 J。画曲线前建议先把表里数据归一化把磁化定标到Msat上这样不同材料参数的曲线才能横向比较。看 OVF 快照时要重点留意磁矩空间分布是否连续、有没有出现尖锐的边界或莫名的纹理。如果出现了教科书上没有的怪异结构先别急着下结论检查一下网格是否够细、初态是否合理、是否真的收敛到平衡。我在调试斯格明子模拟时经常发现某个参数下斯格明子半径反复跳动后来确认是收敛条件不够严增加驰豫时间后结构就稳定了。另外提一句不要只盯着磁化分布能量曲线是判断收敛的好帮手。如果能量曲线在运行末尾还在明显下降那说明还没有到达平衡态。碰到这种情况要么把relax()时长加长要么检查阻尼是否过小。4. 常见问题与排查技巧实录4.1 报错速查表Mumax3 报错信息本身不算友好但我把这几年常遇见的错误整理成了一个小表格照着排查能省下很多时间现象可能原因处理方法启动时报GPU相关错误显存不足缩小网格数、降低Nz、使用-o指定输出目录避免主目录积压脚本第一行就报语法错:和混用检查变量定义与参数赋值运行中输出NaN时间步长过大或参数发散细化网格、降低外场强度、增大阻尼alpha模拟结果出现大量高频条纹网格太粗或交换长度没满足按lex sqrt(2*Aex/(mu0*Msat^2))算网格上限生成文件巨大导致磁盘满保存频率过高降低tableautosave频率二进制 OVF 替代文本 OVF数据结果和文献对不上单位换算或材料参数数量级错误逐一核对 SI 单位用交换长度做自检这张表我直接放在实验室共享目录里新同学报错了先查表解决不了再找我。效果还挺好至少同类问题不会反复占用我时间。4.2 单位与数量级最常见的隐形杀手如果说我只让新同学记一条经验那就是单位。Mumax3 没有图形界面帮你检查单位写错了也不会立刻报错而是给你一份看起来合理、实际上完全偏离物理的结果。数量级出错的典型场景就是Aex坡莫合金一般在1e-11量级有人写成1e-12交换长度直接缩水三倍多磁畴结构当然不对。怎么自查我推荐一个简单方法每次设完Msat和Aex手动算一下交换长度lex再检查你的网格尺寸是否小于等于lex的 1/2 到 1/3。如果网格比lex大不管后面怎么调都白搭。再比如时间尺度磁化动力学的特征时间在皮秒到纳秒量级如果你run(1e-6)相当于跑了一微秒对微磁模拟来说通常长得离谱很可能绕了很多圈振荡。各向异性强度也会有坑。Ku1的单位是 J/m³数值通常从几千到几百万不等取决于材料体系。如果同时有静磁能和交换能三者量级关系直接决定磁畴结构。遇到结果不合常理先核对这三项能量在输出里的数值哪一个异常问题就出在哪。4.3 收敛、发散与伪物理结果模拟发散不一定全是代码错很多时候是物理设定本身矛盾。比如初始状态和外场方向冲突太大磁矩在最初的几步里会发生剧烈进动如果阻尼又小就有可能在数值上冲过头出现NaN。我的习惯是复杂模型第一次跑之前先把外场设成零跑一个简单的relax()看看基础状态是否自洽没问题再加外场。伪物理结果比发散更隐蔽。比如随机初态得出来的多畴结构可能只反映数值噪声而不是真实的能量最低态。体面的做法是换几个随机种子或者用均匀初态算一遍两者结果一致才敢信。再比如边界效应尺寸不大的结构边界会显著改变磁矩分布看结果时要把边界区域和非边界区域区别开。还有一个常见问题是“收敛过度设成绝对化”。不同目标量收敛的判定标准不一样求能量最小值时能量变化趋于零就够但要算矫顽力得保证磁化曲线整体稳定后再扫场。多打印能量和磁化随迭代步数的变化能帮你判断什么时候该停。4.4 批量任务的工程化管理参数优化一旦进入批量阶段问题就不只是物理还有项目管理。我早期用脚本一次性提交 30 组任务跑到第 15 组时磁盘满了前面 14 组结果也被连带搞乱。所以现在我的流程是每组任务一个独立输出目录文件名带参数标签定期清理大体积 OVF只保留关键快照和 Table 文件。建议在模板自带的说明注释里写明这个模板适用的模型、验证过的材料参数、已知的注意事项。这样一来模板不只是代码更像一本“运行手册”。新手拿到后照着注释改参数出错概率大幅下降。课题组里统一这么维护个人的经验也就自然转成了团队的资产。一个小技巧分享每次仿真完我会把脚本、参数、输出摘要压缩成同一个编号的归档目录并在 Table 文件头部记录运行时间和参数组合。几个月后你再回来看不用猜当初用的什么参数文件自己会说话。走完这一步Mumax3 的代码模板才算真正形成了从设置到优化的闭环。