[论文解读] Numerical Simulations of the Invar Effect in Fe-Ni, Fe-Pt, and Fe-Pd Ferromagnets

F. Liot, C. Hooley|arXiv (Cornell University)|Aug 14, 2012

Magnetic Properties and Applications被引用 5

一句话总结

本研究提出一种两阶段计算框架，用于研究铁磁性 Fe-Ni、Fe-Pt 和 Fe-Pd 合金中的 Invar 效应。该框架结合了用于描述温度依赖磁性的 Ising 模型与基于密度泛函理论（DFT）的自由能计算（作为晶格常数的函数），成功再现了实验中观测到的异常低热膨胀系数，包括在 x≈0.35 附近的 Fe-Ni 合金中著名的 Guillaume 效应。

ABSTRACT

The Invar effect in ferromagnetic Fe-Ni, Fe-Pt, and Fe-Pd alloys is investigated theoretically by means of a computationally efficient scheme. The procedure can be divided into two stages: study of magnetism and calculations of structural properties. In the first stage, an Ising model is considered and fractions of Fe moments which point up as a function of temperature are determined. In the second stage, density-functional theory calculations are performed to evaluate free energies of alloys in partially disordered local moment states as a function of lattice constant for various temperatures. Extensive tests of the scheme are carried out by comparing simulation results for thermal expansion coefficients of Fe1-xNix with x = 0.35, 0.4, ..., 0.8, Fe0.72Pt0.28, and Fe0.68Pd0.32 with measurements. The scheme is found to perform well, at least qualitatively, throughout the whole spectrum of test compounds. For example, the significant reduction of the thermal expansion coefficient of Fe1-xNix as x decreases from 0.55 to 0.35 near room temperature, which was discovered by Guillaume, is reliably reproduced. As a result of the overall qualitative agreement between theory and experiment, it appears that the Invar effect in Fe-Ni alloys can be investigated within the same computational framework as Fe-Pt and Fe-Pd.

研究动机与目标

• 开发一种计算高效的框架，用于研究铁基铁磁性合金中的 Invar 效应。
• 探究 Fe-Ni、Fe-Pt 和 Fe-Pd 中的 Invar 效应是否可由单一理论模型描述。
• 检验假设：在 x≈0.35 附近，Fe-Ni 合金中热膨胀系数的异常降低源于磁致结构稳定化。
• 通过将模拟的热膨胀系数与多种成分的实验数据对比，评估该方法的有效性。
• 确定尽管化学成分不同，Fe-Ni、Fe-Pt 和 Fe-Pd 中 Invar 效应是否由相同的物理机制主导。

提出的方法

• 采用两阶段计算方法：首先，利用 Ising 模型计算随温度变化的 Fe 自旋向上占据分数，以模拟磁性。
• 其次，通过 DFT 计算评估在无序局部磁矩（DLM）态下，合金自由能随晶格常数的变化，温度范围覆盖多个温度点。
• 该方法引入 Debye-Grüneisen 模型，以考虑晶格振动对自由能的贡献。
• 通过在每个温度下对晶格常数最小化自由能，确定平衡晶格常数。
• 该方案避免了对非共线磁性或部分化学有序性的显式处理，转而依赖有效磁矩和基于 DFT 的能量计算。
• 热膨胀系数由平衡晶格常数随温度的变化关系推导得出。

实验结果

研究问题

• RQ1 能否使用统一的计算框架描述 Fe-Ni、Fe-Pt 和 Fe-Pd 合金中的 Invar 效应？
• RQ2 当 Fe1-xNi x 中的 x 从 0.55 降低到 0.35 时，其在室温附近热膨胀系数的降低是否如模型预测的那样源于磁涨落？
• RQ3 DFT-Ising 方法在多大程度上再现了无序面心立方（fcc）铁基铁磁体中的实验热膨胀趋势？
• RQ4 随着温度升高，磁矩的减小是否是这些体系中 Invar 效应的关键驱动力？
• RQ5 尽管电子结构不同，该相同的物理机制是否能解释 Fe-Ni、Fe-Pt 和 Fe-Pd 中的 Invar 行为？

主要发现

• 该模型成功再现了 Fe0.65Ni0.35、Fe0.72Pt0.28 和 Fe0.68Pd0.32 在广阔温度范围内的异常低热膨胀系数。
• 当 x 从 0.55 降低到 0.35 时，Fe1-xNi x 在室温附近热膨胀系数的显著降低——这一现象最早由 Guillaume 发现——在模拟中被可靠再现。
• 在 x=0.35 时，Fe1-xNi x 在零温下的计算晶格常数与 Vegard 定律的偏差为 -0.14%，与实验值 -0.08% 非常接近。
• 该模型在 x=0.35 至 0.80 范围内以及 Fe0.72Pt0.28 和 Fe0.68Pd0.32 在不同温度下，与实验热膨胀系数表现出定性一致。
• 该框架成功捕捉了磁矩随温度的变化及其对结构稳定性的影响力，表明磁性在 Invar 效应中起核心作用。
• 尽管采用了简化假设（如忽略静态离子位移和非共线磁性），该方法的结果仍与实验观测高度一致，支持其在研究 Invar 效应中的实用性。