Lattice vacancy migration barriers in Fe-Ni alloys, and why Ni atoms diffuse slowly: An ab initio study

This ab initio study reveals that nickel atoms diffuse significantly slower than iron atoms in Fe-Ni alloys due to a coupling between local lattice distortions and spin-polarized electronic structure, which causes iron to relax into vacancies while nickel remains rigidly fixed.

The Gist

The Big Picture: The "Slow-Motion" Mystery of Nickel

Imagine you have a crowded dance floor made of two types of dancers: Iron (Fe) and Nickel (Ni). They are holding hands in a specific pattern. Sometimes, a dancer leaves the floor, creating an empty spot (a vacancy). To keep the dance going, a neighbor has to step into that empty spot.

Scientists have known for a long time that in this specific alloy (Fe-Ni), the Iron dancers are fast, while the Nickel dancers are incredibly slow. It's like Iron is a sprinter and Nickel is a turtle.

This matters because there is a special, super-strong magnetic material called Tetrataenite (found in meteorites) that scientists want to make in factories to replace rare-earth magnets in things like electric cars and wind turbines. To make this material, the Iron and Nickel atoms need to swap places and organize themselves into a neat, ordered pattern. But because the Nickel atoms are so slow, this process takes forever (or requires extreme heat that ruins the material).

The Question: Why is Nickel so much slower than Iron?

The Answer (from this paper): It's not just about size; it's about their "magnetic personalities" and how they react to an empty space.

---

The Investigation: A Digital "Push and Pull" Experiment

The researchers used a supercomputer to simulate this dance floor. They didn't just guess; they used a method called NEB (Nudged Elastic Band).

Think of the NEB method like a movie reel.

1. Frame 1: An atom is in its original spot.
2. Frame 10: The atom is halfway into the empty spot.
3. Frame 20: The atom has fully swapped places.

The computer calculates exactly how much "energy" (effort) is needed to push the atom from Frame 1 to Frame 20. This effort is called the migration barrier.

The Result:

• Iron atoms: The barrier was low. It was easy to push them.
• Nickel atoms: The barrier was 40% higher. It was much harder to push them.

This confirmed the experimental observation: Nickel is sluggish. But why?

---

The "Aha!" Moment: The Rubber Band vs. The Rock

The researchers discovered the secret lies in how the atoms react when a neighbor leaves the room.

1. The Iron Atom: The "Relaxing" Dancer

When an Iron atom sees an empty spot next to it, it doesn't just sit there. It stretches and relaxes into the empty space.

• The Analogy: Imagine Iron is wearing a rubber band suit. When a neighbor leaves, the rubber band snaps, and the Iron atom stretches out to fill the gap.
• The Physics: As Iron stretches into the empty space, its magnetic field gets stronger. This is a "happy" state for the Iron atom; it lowers its energy. Because it moves easily and gains energy by moving, it's very mobile.

2. The Nickel Atom: The "Rigid" Dancer

When a Nickel atom sees an empty spot, it does nothing. It stays stiff.

• The Analogy: Imagine Nickel is a rock or a statue. Even when a neighbor leaves, the Nickel atom refuses to budge. It stays rigidly fixed in its original spot.
• The Physics: Nickel's electronic structure (its internal "magnetic wiring") doesn't allow it to stretch into the empty space without paying a huge energy penalty. It's energetically "expensive" for Nickel to move. So, it stays put.

The Consequence: The "Long Walk" Problem

Because Iron atoms stretch into the empty spot before they even start their journey, the actual distance they need to travel to swap places is shorter.

• Iron: "I'm already halfway there because I stretched out. Easy swap!" (Low energy cost).
• Nickel: "I am stuck in my original spot. I have to drag my whole heavy body all the way across the gap." (High energy cost).

The paper shows a direct link: The further an atom has to physically travel to swap places, the more energy it takes. Since Nickel refuses to stretch, it has to travel the full distance, making the "ticket price" (energy barrier) very high.

Why Does This Matter?

1. Meteorites vs. Factories: This material (Tetrataenite) forms naturally in meteorites over millions of years because the atoms move so slowly. In a factory, we need to make it in hours or days.
2. The Bottleneck: The "slow Nickel" is the bottleneck. If we can't figure out how to make Nickel move faster (or find a way around it), we can't mass-produce these powerful, rare-earth-free magnets.
3. Future Tech: Understanding this "magnetic stiffness" of Nickel helps scientists design better ways to manufacture these magnets, perhaps by using pressure or specific magnetic fields to force the Nickel to relax.

Summary in One Sentence

Iron atoms are like flexible rubber bands that stretch into empty spaces to make swapping easy, while Nickel atoms are like rigid rocks that refuse to move, making them incredibly slow to diffuse and hard to organize into the super-strong magnets we need for the future.

Technical Summary

1. Problem Statement

The diffusion of solute elements in alloys is critical for metallurgical processes such as chemical ordering, segregation, and precipitation. In the Fe-Ni binary system, experimental observations have long established that Nickel (Ni) atoms diffuse significantly slower than Iron (Fe) atoms, particularly in the ferromagnetic regime.

This sluggish diffusion of Ni presents a major challenge for the bulk synthesis of the $L1_0$ ordered phase (tetrataenite), a tetragonal, atomically ordered structure with high magnetocrystalline anisotropy. This phase is a promising candidate for rare-earth-free permanent magnets. However, the low atomic ordering temperature and the slow diffusion kinetics of Ni prevent the efficient laboratory synthesis of bulk $L1_0$ FeNi. While previous studies (both experimental and computational) have confirmed the mobility difference, a clear atomic-scale physical mechanism explaining why Ni atoms are less mobile than Fe atoms in ferromagnetic Fe-Ni alloys has been lacking.

2. Methodology

The authors employed ab initio electronic structure calculations based on Density Functional Theory (DFT) to investigate lattice vacancy migration barriers.

• System Configuration: The study covered ferromagnetic $Fe_xNi_{1-x}$ alloys ($0.4 \le x \le 0.6$) in three structural states:
  • Atomically disordered face-centered cubic (A1) phase.
  • Partially ordered structures.
  • Fully ordered tetragonal $L1_0$ phase.
• Supercell Generation: Representative alloy configurations were generated using Metropolis Monte Carlo (MC) simulations based on a Bragg-Williams Hamiltonian with effective pair interactions (EPIs) derived from KKR-CPA (Korringa-Kohn-Rostoker Coherent Potential Approximation) theory. This ensured accurate sampling of chemical disorder and short-range order.
• Barrier Calculation: The Nudged Elastic Band (NEB) method was implemented within the CASTEP code to calculate the energy barriers for lattice vacancy migration.
  • Process: A vacancy was created by removing an atom. A neighboring atom (either Fe or Ni) was selected to swap with the vacancy.
  • Optimization: Both initial and final states were geometry-optimized to find local energy minima.
  • Dataset: A total of 192 migration barriers were calculated (103 Fe-vacancy interchanges and 89 Ni-vacancy interchanges) across various compositions and ordering degrees.
• Electronic Structure Analysis: The study analyzed local structural relaxations, spin-polarized electronic structures, and magnetic moments to correlate atomic movement with electronic changes.

3. Key Contributions

The paper provides the first comprehensive atomic-scale explanation for the disparity in diffusion rates between Fe and Ni in ferromagnetic Fe-Ni alloys. The primary contribution is the identification of a coupling between local lattice distortions and spin-polarized electronic structure.

Specifically, the authors demonstrate that:

1. Fe atoms "relax" into vacancies: When a vacancy forms, Fe atoms significantly displace toward the empty site to optimize bonding and lower energy.
2. Ni atoms remain rigid: Ni atoms show minimal displacement, remaining fixed near their original lattice positions.
3. Magnetic Origin: This structural difference is driven by the electronic structure. The relaxation of Fe atoms increases their local magnetic moment (spin polarization), which lowers the system's total energy. Ni atoms, due to their filled minority-spin $t_{2g}$ states, lack the electronic flexibility to undergo similar relaxation without an energy penalty.

4. Key Results

A. Migration Barrier Heights

• Ni is significantly less mobile: The average lattice vacancy migration barrier for Ni-vacancy interchanges is $1.13 \pm 0.12$ eV.
• Fe is more mobile: The average barrier for Fe-vacancy interchanges is $0.78 \pm 0.11$ eV.
• Magnitude: The Ni barrier is approximately 40% higher than the Fe barrier. According to the Arrhenius equation, this substantial difference in activation energy explains the exponentially slower diffusion rates of Ni observed experimentally.
• Consistency: This trend holds across all compositions ($x=0.41$ to $0.59$) and ordering states (disordered A1 to ordered $L1_0$).

B. Structural Relaxation

• Displacement Analysis: Upon vacancy creation, Fe atoms relax by an average of $0.0428$ Å (in the ordered $L1_0$ phase), whereas Ni atoms relax by only $0.0127$ Å.
• Correlation: There is a positive correlation between the distance an atom travels during the interchange and the barrier height. Fe atoms, which move shorter distances (due to pre-relaxation into the vacancy), encounter lower barriers. Ni atoms, which do not relax, must traverse the full nearest-neighbor distance, resulting in higher barriers.

C. Electronic and Magnetic Mechanism

• Spin Polarization: The introduction of a vacancy alters the local electron density. Fe atoms adjacent to vacancies experience an increase in magnetic moment (e.g., from $2.805 \mu_B$ to $2.937 \mu_B$ in the $L1_0$ phase).
• Orbital Physics: In the minority-spin channel, Fe has partially filled $t_{2g}$ states near the Fermi level, allowing them to form new orbital combinations and relax into the vacancy space. Ni has largely filled $t_{2g}$ states well below the Fermi level, restricting its ability to rearrange electronically and structurally.
• Total Magnetization: The geometry optimization of vacancy-containing cells results in an average increase in total supercell magnetization of $0.2 \pm 0.1 \mu_B$, primarily accruing on Fe atoms near the vacancy.

D. Influence of Ordering

• In the ordered $L1_0$ phase, barriers differ slightly between in-layer and inter-layer swaps, suggesting potential anisotropy in diffusion coefficients. However, the fundamental difference between Fe and Ni mobility persists regardless of the degree of atomic order.

5. Significance and Implications

• Fundamental Physics: The study resolves a long-standing question in materials science regarding the origin of diffusion asymmetry in Fe-Ni alloys, linking it directly to spin-polarized electronic structure and local magnetic moment changes.
• Materials Design: Understanding that Ni diffusion is kinetically hindered by high energy barriers due to its electronic rigidity provides a target for future processing strategies. To synthesize bulk $L1_0$ FeNi, methods must be developed to either lower these barriers (e.g., via strain, external fields, or defect engineering) or bypass the kinetic limitations.
• Computational Resource: The authors generated a high-quality dataset of ~200 NEB calculations. This dataset serves as a critical training and validation set for developing Machine-Learned Interatomic Potentials (MLIPs). Since semi-empirical potentials often fail to capture these specific magnetic and electronic barrier heights accurately, MLIPs trained on this ab initio data will enable large-scale kinetic Monte Carlo (kMC) simulations to model the disorder-order transition in Fe-Ni alloys with unprecedented accuracy.

In conclusion, the paper establishes that the sluggish diffusion of Ni in Fe-Ni alloys is not merely a geometric or mass effect but is fundamentally driven by the electronic rigidity of Ni atoms in the minority-spin channel, which prevents the energy-lowering structural relaxations that Fe atoms readily undergo.