Nd2Fe14B磁体因其高的剩磁（Br）、高矫顽力（Hcj）以及优异的最大能量积（BH）max而成为最出色的永磁材料之一[[1], [2], [3], [4]]。其独特的磁性能使其在高性能电机、磁共振成像（MRI）、传感器、电子设备和现代绿色能源技术中得到广泛应用[[5], [6], [7], [8], [9]]。自1983年Sagawa [10]和Croat [11]分别通过粉末冶金法和熔融旋压法成功制备出以Nd2Fe14B为主相的稳定Nd-Fe-B系列永磁材料以来，Nd2Fe14B磁体的研究和应用进入了快速发展阶段。这些材料的卓越性能不仅推动了永磁材料科学的发展，也为相关研究开辟了新的领域。
根据不同的制造工艺，Nd-Fe-B磁体主要可分为三类：烧结Nd-Fe-B磁体、粘结Nd-Fe-B磁体和热变形Nd-Fe-B磁体。烧结Nd-Fe-B磁体通过粉末冶金法制备，先熔化合金、制粉，然后在惰性气体或真空环境中压制成型；通常还需要进行老化热处理以提高结晶度[12,13]。粘结Nd-Fe-B磁体则是将Nd-Fe-B粉末与塑料或其他粘结材料混合后，通过注塑或压制成型[14]。热变形过程主要包括两个关键步骤：首先对纳米晶Nd-Fe-B薄片或粉末进行热压处理，形成各向同性的致密前驱体；随后在热变形过程中，材料发生模具变形，从各向同性转变为各向异性，因为外部压力方向与Nd2Fe14B晶粒的c轴晶体取向一致[15]。与其他制造技术相比，热变形工艺是诱导纳米晶硬磁材料织构的最可行方法。与传统烧结工艺相比，热变形工艺具有更高的矫顽力，并且能够生产出与最终产品形状和尺寸高度匹配的材料，从而减少了后续加工步骤，缩短了生产周期[16,17]。
近年来，人们对Nd-Fe-B永磁体热变形过程中的变形机制有了显著进展。众多研究[[18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29]]表明，热变形本质上是一个高温蠕变过程，晶界扩散和晶界滑移是控制塑性变形的两个主要机制。Grünberger等人[28]通过“固溶-沉淀扩散蠕变”模型系统地阐明了热变形的蠕变本质，指出高温下晶界相的熔化不仅增强了晶界扩散，还显著降低了晶界滑移的阻力。这一观点得到了多项实验研究的验证。例如，Yasuda等人[29]观察到Nd-Fe-B合金在拉伸变形过程中的优异伸长主要来源于晶界滑移；Liu等人[18]通过实验表明，在高温高压条件下，晶粒旋转行为主要由晶界滑移机制控制。Hioki等人的工作[20]通过系统实验观察进一步支持了这一观点，他们提出在热压过程中，晶界相首先熔化，随后在热变形阶段，施加的应力促进晶粒沿c轴（片状晶粒的宽表面）的各向异性生长和晶粒旋转。总体而言，这些研究确立了热变形的蠕变本质，并强调了晶界扩散和晶界滑移机制的关键作用；然而，在特定温度和压缩应力条件下哪种机制占主导地位仍不清楚，需要进一步研究。
分子动力学（MD）模拟是阐明金属和合金系统蠕变机制的强大工具[[30], [31], [32], [33], [34]]。Wang等人[[30], [31], [32]]通过MD模拟成功再现了单晶和多主元素合金的蠕变过程，利用经典的Mukherjee-Bird-Dorn本构方程进行多温度蠕变测试，确定了应力指数（n）和蠕变激活能等关键参数，明确了不同蠕变机制的主导作用，并建立了蠕变行为的应力-温度依赖的本构关系。类似地，Qian等人[34]通过MD模拟在多种应力水平下进行了系列蠕变实验，系统研究了不同蠕变机制之间的竞争关系。这些开创性研究为当前工作提供了重要的方法论参考。因此，本研究采用MD模拟来研究纳米晶Nd-Fe-B在热压过程中的蠕变行为。通过构建精确的原子模型，模拟高温高压条件下的热机械变形，我们旨在：阐明外部载荷和温度条件对蠕变机制转变的影响，并深入理解蠕变过程中的塑性变形机制转变。该研究将有助于丰富Nd-Fe-B在加工条件下的蠕变机制基础知识。
