🔬 materials science

Development of a magnetic interatomic potential for cubic anti-ferromagnets: the case of NiO

本文提出了一种通过整合海森堡交换模型与内尔模型来开发立方反铁磁体磁性原子间势能的新颖方法，并成功将其应用于 NiO 以创建经过验证的势能，从而能够实现对耦合磁弹性现象的大规模分子动力学模拟。

原作者： Ievgeniia Korniienko, Pablo Nieves, Jakub Sebesta, Roberto Iglesias, Dominik Legut

发布于 2026-07-20

📖 1 分钟阅读☕ 轻松阅读

原作者： Ievgeniia Korniienko, Pablo Nieves, Jakub Sebesta, Roberto Iglesias, Dominik Legut

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下固体内部那个隐形的微观世界，就像一座由原子组成的微型繁华城市。在这座城市里，原子并非静止不动；它们在不断地抖动、碰撞，并以特定的方式“手拉手”。但在某些特殊的材料中，这些原子拥有一种秘密超能力：它们表现得像微型磁铁。通常，当我们想到磁铁时，脑海中浮现的是吸附在金属上的冰箱贴。但在一种被称为“反铁磁体”的不同类型的材料中，邻居们就像是在玩一场拔河游戏，每个人都向相反的方向拉扯。一个原子指向“上”，它的邻居指向“下”，它们彼此抵消，使得整个物体对外界表现不出磁性。

科学家们对这些材料非常感兴趣，因为当施加磁场或挤压它们时，它们可以改变形状。这被称为“磁弹性效应”。这就像是你捏一个压力球，它突然决定指向北方，或者你拿着一块磁铁在它附近挥动，它就缩减了体积。为了理解这一现象，研究人员使用强大的计算机模拟技术。他们构建了一个被称为“原子间势函数”的“配方”，这本质上是一套指令，告诉计算机原子应该如何推、拉以及对磁力做出反应。困难之处在于，对于这种“拔河”性质的材料，其配方一直缺少一个关键成分：磁性规则。没有这些规则，计算机模拟就像是在不知道风向的情况下预测天气。

这篇论文旨在为一种名为氧化镍（NiO）的特定材料编写那本缺失的“配方书”，它是这种“拔河”磁体的典型代表。来自捷克和西班牙的研究团队开发了一种新方法，将磁性规则与原子移动和拉伸的规则结合起来。他们创建了两个不同版本的配方。一个版本是基于离子固体（如食盐）行为的经典理论，另一个则是利用高水平计算机计算从零开始构建的全新定制配方。

团队通过模拟一个拥有超过8.5万个原子的庞大“城市”来测试他们的新配方。他们检查了当磁序发生变化时，模拟出的材料是否会以恰当的幅度收缩或膨胀，以及它在面对压力时的反应是否正确。结果令人振奋：两种配方都以极高的准确度预测了材料的行为，与超级计算机进行的复杂计算相吻合。具体而言，他们发现材料在磁自旋有序化时会收缩约0.14%，并且成功重现了材料在不同方向上承受挤压时的抗性。

然而，作者谨慎地指出，目前这些结果仅适用于零温（一个原子停止抖动的理论状态）下的模拟。虽然这些模型在计算机中表现完美，但科学家们建议，若要了解这些材料在真实的、温暖的世界中的表现，仍需开展更多工作。他们还指出，虽然他们的新型定制配方（称为RF-MEAM）在预测材料抗扭转能力方面比旧配方更好，但关于该材料的磁性拉伸强度与其他研究相比究竟有多强，目前仍存在一些谜团。

简而言之，这篇论文并不声称已经解开了所有关于磁性材料的谜团，但它为未来的探索者搭建了一座非常坚实且可靠的桥梁。通过成功地将磁性规则与氧化镍的运动规则相结合，他们为更大规模、更详细的模拟打开了大门。这意味着在未来，科学家们或许能够设计出更好的传感器、执行器，甚至是利用这些聪明的、能改变形状的磁性材料开发的量子技术，而这一切，都源于他们终于掌握了一份关于原子如何共同起舞的良好地图。

问题陈述
原子间势函数是磁性材料分子动力学（MD）模拟的基础；然而，将磁性特征整合到复杂反铁磁（AFM）材料的势函数中仍然是一个重大挑战。虽然结合了经典自旋与分子动力学的自旋-晶格模型（SD-MD）已在立方铁磁（FM）材料中得到成功构建和验证，但将其扩展到反铁磁系统却非常困难。其核心难点在于，该模型需要同时在可变原子间距离下重现磁动力学，并在存在磁相互作用的情况下，由原子间势函数控制弹性性质。目前，针对原型立方反铁磁材料氧化镍（NiO）的可得原子间势函数可以描述电荷和电场响应，但缺乏磁相互作用。这一局限性阻碍了对金属氧化物中耦合电场与磁场现象（如磁弹性效应和磁电效应）的研究。

方法论
作者提出了一种通过在标准原子间势函数中添加磁哈密顿量，从而将磁性属性整合进立方反铁磁材料原子间势函数的方案。总哈密顿量（ $H_{sl}$ ）包括动能、非磁性原子间势（ $V$ ）和磁性势（ $V_{mag}$ ）。

磁哈密顿量： 磁项包括海森堡交换相互作用和内尔（Néel）相互作用（ $H_{N\acute{e}el}$ ）。内尔相互作用通过偶极（ $l$ ）项和四极（ $q$ ）项来解释磁各向异性，这两者分别描述了各向异性磁致伸缩和磁晶各向异性（MCA）。
贝特-斯莱特（Bethe-Slater）参数化： 交换作用（ $J$ ）、偶极（ $l$ ）和四极（ $q$ ）参数的空间依赖性使用贝特-斯莱特曲线进行建模。针对反铁磁情况，作者基于具有反平行自旋的简单立方（SC）结构的内尔模型，推导出了这些参数的具体表达式。
NiO 特性： 鉴于 NiO 晶胞的复杂性（面心立方 Ni 子晶格，包含混合的平行和反平行最近邻原子），该模型通过仅关注 Ni 原子第二近邻（NNN）来简化相互作用。选择这一做法的理由是，NiO 中的 NNN 交换作用比最近邻相互作用更强，且 NNN Ni 子晶格呈现出一种所有自旋均为反平行的简单立方结构，从而简化了参数推导。
势函数构建： 构建了两种不同的原子间势函数以验证该方法论：
- SD-MD 1： 基于离子固体玻恩（Born）模型（Fisher 和 Matsubara），利用 Buckingham 势处理短程相互作用。
- SD-MD 2： 基于无参考修正嵌入原子方法（RF-MEAM），拟合自通过含 Hubbard 校正的 45 个扭曲 NiO 晶体结构的密度泛函理论（DFT）数据集。
输入参数： 用于贝特-斯莱特曲线的关键参数（内尔温度、弹性常数、自发体积磁致伸缩、各向异性磁弹性常数和 MCA 常数）是通过结合实验数据和从头算计算（使用 VASP 和 MAELAS 代码）推导而来的。

主要贡献

方法论框架： 本文建立了一套通用的程序，用于开发包含海森堡交换作用和内尔相互作用的立方反铁磁材料原子间势函数，特别解决了建模反平行自旋配置的挑战。
双势函数验证： 作者成功地将该框架应用于 NiO，创建了两个包含磁弹性相互作用的不同模型（SD-MD 1 和 SD-MD 2）。
参数推导： 该工作提供了反铁磁系统贝特-斯莱特参数的详细推导过程，将这些参数与宏观性质（如内尔温度、自发体积磁致伸缩和各向异性磁弹性常数）联系起来。
定制势函数开发： 研究展示了如何为 NiO 开发定制的 RF-MEAM 势函数，该势函数在重现 DFT 计算的弹性常数（特别是 $C_{44}$ ）方面，比传统的玻恩模型势函数表现更好。

结果
开发的模型在零温下针对 DFT 计算和实验数据进行了验证：

力学与磁学性质： 两种 SD-MD 模型在力学性质（弹性常数 $C_{11}, C_{12}, C_{44}$ 、体积模量）和磁学性质（内尔温度、自发体积磁致伸缩 $\omega_s$ 和磁晶各向异性常数 $K_1$ ）方面均与 DFT 及实验值高度一致。
磁致伸缩： 模型准确重现了输入的体积磁致伸缩（ $\omega_s \approx -0.0014$ ）和各向异性磁致伸缩系数（ $\lambda_{100}, \lambda_{111}$ ）。
压力依赖性： 模型正确捕捉到了磁晶各向异性常数 $K_1$ 对静水压力的线性依赖关系，与 DFT 计算所得的斜率相吻合。
势函数对比： 与玻恩模型（SD-MD 1）相比，基于 RF-MEAM 的模型（SD-MD 2）在重现剪切弹性常数 $C_{44}$ 及其产生的磁致伸缩系数方面表现出更优越的性能，解决了后者在处理 NiO 时的已知局限性。

意义
本文提出了开发专门针对立方对称反铁磁材料并包含磁学性质的原子间势函数的“首次尝试”。作者声称这项工作是该方法论在零温度下的概念验证。其主要意义在于能够实现大规模的磁弹性现象分子动力学模拟，而此前由于缺乏合适的势函数，这一能力受到了限制。作者指出，这些模型为研究金属氧化物中的耦合电场与磁场效应、磁弹性与磁电耦合、磁性纳米结构的形状效应以及磁振子-声子动力学开辟了道路。该工作在声明上保持了适度，指出对于有限温度研究的结果应谨慎解读，因为这需要进一步的详细研究。

类似论文