永磁铁是现代生活背后许多环节的隐形引擎，为从电动汽车电机到发电风力涡轮机的各种设备提供动力。在这些最强大的磁体核心，是一种被称为钕铁硼的材料。虽然磁体的本体是由这种材料的大晶体组成的，但其强度的秘密在于这些晶体相遇时形成的薄且肉眼不可见的边界。这些边界就像砖块之间的砂浆；如果它们脆弱或紊乱，整个结构就会失去抵抗退磁的能力，尤其是在受热时。近年来，科学家们发现了一种通过分解并重新组装旧磁体来进行回收的方法，但这个过程往往会打乱这些边界的精细排列，导致新磁体比旧的更弱。研究人员发现，在回收过程中添加微量的锆和硼可以帮助解决这个问题，使微小的针状颗粒恰好在晶体交界处形成。这些颗粒似乎能将结构固定在一起，但直到现在，人们仍不知道为什么自然界选择在这些特定位置构建它们，而不是将它们随机散布在整个材料中。
为了解开这个谜团，一个研究小组转向了一种强大的计算机模拟技术，这种技术允许他们观察原子的行为，就像它们是真实存在的一样，而无需物理实验室。他们专注于一个具体问题：为什么仅存在于极微量中的锆和硼会聚集在晶体边缘，而不是留在中间混合在一起。使用一种能够解释所涉及原子复杂磁行为的方法，他们构建了磁体内部及其边界的数字模型。随后，他们测试了如果将硼、锆和其他元素的单个原子放置在这些不同位置会发生什么。目标是测量将这些原子保持在晶体中间与将其保持在边界所需的能量。在原子的世界里，系统自然寻求最低能量状态，就像球向山谷滚动一样。通过计算这些能量水平，研究人员可以预测原子更倾向于在哪里沉降。
模拟结果揭示了一个清晰且令人惊讶的偏好等级。研究人员发现，硼原子有一种非常强烈的、自然的趋向，即奔向边界，通过这样做显著降低其能量。锆原子表现出类似的、虽然稍弱一点的向边缘移动的欲望。当研究人员将这两种元素放在一种模拟真实磁体中粒子结构的配置中时，组合体仍然偏好边界。这意味着，即使磁体中的锆总量极少——仅约0.1%——材料也会自然地将其引导至晶体边缘，使其能够与硼汇合。计算机模型显示，一旦这些元素在边界相遇，它们就会形成一种稳定的、低能量的排列，有助于锁定结构。这解释了为什么含有极少量锆的磁体仍能产生独特的、集中的颗粒于其边界：边界本身就像是一个针对这些特定原子的磁铁，将它们从主体材料中拉出并聚集在一起。
该研究还观察了这些原子变化如何影响磁体的整体强度。磁体科学中的一个主要担忧是，添加外来元素可能会破坏铁原子的磁对齐，从而削弱磁体。然而，模拟显示，当锆和硼聚集在边界时，它们并不会干扰周围晶体的磁有序性。晶体内部的铁原子保持强烈的对齐，从而保留了磁体的力量，而边界区域则变成了一个稳定且经过加固的区域。这一发现将原子的微观行为与磁体的宏观性能联系了起来。它表明，回收磁体之所以能重新变强，不仅是因为有了新的粒子，更是因为材料的热力学特性自然地驱动了必要的成分到达正确的位置。边界充当了一个优选的储库，收集稀缺元素并将其组织成一种结构，从而稳定晶界，防止晶体生长得过大并失去其磁力。
虽然计算机模型提供了一个关于这些原子为何聚集的清晰图景，但研究人员指出，实际粒子的形成涉及一个随时间变化的复杂热量与运动过程，而模拟并未完全捕捉到这一点。该研究专门侧重于设定这一过程基础的能量偏好。它证实了边界是这些元素在能量上更有利的场所，但并未计算它们移动的速度或最终粒子结晶的具体步骤。研究结果为实验室中观察到的现象提供了一个根本性的解释：为什么少量的锆可以带来磁体结构在组织上的巨大提升。通过理解边界是这些元素的“热力学汇”，科学家可以更好地设计依赖这种自然趋向来生产高质量磁体的回收工艺，确保其中的珍贵材料不仅被回收，而且其潜力得以完全恢复。
技术摘要：回收 Nd2Fe14B 中 Zr–B 在晶界处偏析的第一性原理热力学研究
问题陈述
Nd2Fe14B 永磁体的矫顽力和热稳定性关键取决于其晶界（GB）和晶间（IG）区域的化学成分与结构。在回收的 Nd2Fe14B 永磁体中，仅回收硬磁 Nd2Fe14B 相往往无法恢复实现高性能所需的最佳晶间微观结构。最近关于含有约 0.1 at.% Zr 的回收 Nd2Fe14B 永磁体的实验研究表明，在晶界和三叉晶界处出现了纳米级 ZrB2 沉淀物。这些沉淀物与抑制晶粒粗化和增强矫顽力有关。然而，实验数据无法解释这种优先定域化的原子尺度起源：为什么在块体 Zr 浓度极低且均匀的情况下，Zr 和 B 会在晶界区域聚集并形成类 ZrB2 结构？现有的关于 Nd–Fe–B 界面的第一性原理研究通常局限于单一化学成分、理想化界面，或使用开壳层近似处理稀土 4f 态，这导致在一致的电子框架下理解多组分（特别是 B、Zr 和 Dy）耦合偏析的热力学方面存在空白。
方法论
为了填补这一空白，作者采用了带有哈伯德修正的自旋极化密度泛函理论（DFT+U）来计算缺陷偏析的原子热力学。
- 模型： 本研究利用了代表块体 Nd2Fe14B 和理想化 Σ1 (001) 旋转晶界的 136 原子超胞。晶界模型通过绕 [001] 轴将块体单元的上半部分旋转 90° 来构建。
- 缺陷构型： 分析了五种成分匹配的块体/晶界对：原始态、间隙硼（Bint）、置换锆（ZrNd）、置换镝（DyNd）以及代表类 ZrB2 基元的局部 Zr–B 复合物（Zr–2B）。
- 电子结构： 计算使用 VASP 软件包，采用 PBE 广义梯度近似和旋转不变 DFT+U 公式（Dudarek 等人），其中对 Nd 和 Dy 的 4f 态应用了 Ueff=6 eV。该方法在所有构型中保持了关联电子描述，避免了开壳层近似。
- 验证： 执行了严格的磁态一致性检查。研究验证了所有含缺陷的结构都保持在与原始参考态相同的高矩 Fe 子格范围内，确保能量差反映的是真实的位点偏好，而非不同磁基态导致的伪影。
- 热力学指标： 计算了偏析能（Eseg）以确定位点偏好，并计算了晶界过剩能（γGB）以评估界面稳定性。
关键结果
计算确立了 B、Zr、Dy 和 Zr–B 复合物偏析的明确热力学层级：
- 强硼偏析： 间隙硼表现出最强的晶界偏好，其偏析能为 Eseg(B)=−1.59 eV/缺陷。与块体中的正值（+1.01 eV/B）相比，B 在边界处的形成能变为负值（$-0.58$ eV/B），表明边界是过量 B 的热力学汇。
- 独立的锆偏好： Zr 也表现出直接的界面热力学偏好，其 Eseg(Zr)=−0.77 eV/缺陷。
- Zr–B 配位结构的稳定性： 当 Zr 和 B 以类 ZrB2 的局部构型结合时，界面偏好性得以保留（Eseg=−0.32 eV/缺陷），尽管驱动力弱于孤立物种。这表明随着 Zr–B 配位的形成，界面仍然是有利的。
- 镝的对比： 与 Zr 和 B 不同，在这一理想化的零温模型中，Dy 表现出正的偏析能（Eseg=+1.90 eV/缺陷），这表明在这些特定条件下，Dy 在基体中而非边界处具有能量稳定性。
- 界面稳定性： B、Zr 以及 ZrB2 类基元的存在降低了相对于原始界面的晶界过剩能（分别降低了 $-0.33、-0.16和-0.20Jm^{-2}$），而 Dy 则增加了过剩能。
- 磁性兼容性： 所有定域于边界的缺陷对归一化磁化强度的扰动均小于其块体对应物。Fe 子格矩和自旋分辨 d 态保持与高矩基态的一致性，证实了 Zr–B 偏析不会破坏周围晶粒的磁性特征。
意义与主张
本文声称提供了对回收 Nd2Fe14B 永磁体中观察到的 Zr 和 B 再分布现象的第一性原理热力学解释。主要贡献在于确立了一条热力学路径，即晶界作为硼和锆的优先储库。这解释了为何全球 Zr 浓度极低（~0.1 at.%）却能产生局部浓缩的 ZrB2 沉淀物：边界在能量上倾向于同时积累这两种组分，从而促进了它们的共定域化及随后的类 ZrB2 配位。
作者强调，这些结果阐明了回收 Nd2Fe14B 磁体微观结构演化的“前置原子步骤”。虽然计算并未模拟沉淀物生长的动力学过程或直接计算矫顽力，但它们证明了 ZrB2 富集边界结构的形成是由组分偏析驱动的热力学过程。这种偏析稳定了晶间区域，并且在磁性上与高矩 Nd2Fe14B 相兼容，为含有 Zr 的回收磁体中观察到的增强微观结构稳定性和矫顽力提供了微观依据。该研究将这种热力学驱动力与最终决定如 Dy 等元素分布的复杂加工环境（有限温度、扩散梯度）区分开来。
每周获取最佳 materials science 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题，再试一次？
无垃圾邮件，随时退订。