Effect of Vacancies on Hydrogen Mobility and Trapping in Elemental Fe and Cr: A DFT and kMC Study

This study employs a combined DFT and kMC approach to demonstrate that vacancies significantly reduce hydrogen mobility and increase activation energy in BCC Fe and Cr, with a more pronounced trapping effect observed in Cr due to stronger hydrogen-vacancy interactions.

The Gist

The Big Picture: Hydrogen as a "Tiny, Unwanted Guest"

Imagine a metal structure (like a steel beam in a bridge or a part of a nuclear reactor) as a giant, crowded dance floor made of atoms. Usually, everyone is dancing in perfect, organized rows. But sometimes, a tiny, hyperactive guest named Hydrogen sneaks in.

Hydrogen is very small and moves incredibly fast. While it might seem harmless, if it gets stuck in the wrong spots, it can make the metal brittle and prone to cracking (a problem called "hydrogen embrittlement").

This study asks a specific question: What happens when there are empty spots (vacancies) on this dance floor? Do these empty spots act like traps that catch the fast-moving Hydrogen, or do they let it slip away? The researchers looked at two specific types of metal floors: Iron (Fe) and Chromium (Cr).

The Tools: Two Different Ways to Look at the Problem

To solve this, the scientists used a "multiscale" approach, which is like using two different cameras to film the same event:

1. The Microscope (DFT): They used a super-powerful computer simulation (Density Functional Theory) to zoom in on the atomic level. This let them see exactly how much energy it takes for a Hydrogen atom to jump from one spot to another, or how tightly it gets stuck in an empty spot.
2. The Time-Lapse Camera (kMC): Because atoms move too fast to watch in real-time, they used a Kinetic Monte Carlo (kMC) simulation. Think of this as a time-lapse video that speeds up time by billions of times. This allowed them to watch how Hydrogen moves across a large area over a long period, seeing where it gets stuck and how fast it travels.

Key Findings: The "Trap" Analogy

1. The Empty Seat (Vacancy)

In a perfect metal crystal, every seat is taken. But sometimes, a seat is missing. This is a vacancy.

• The Discovery: Hydrogen loves these empty seats. It gets attracted to them like a magnet.
• The Capacity: Just like a small car can only fit a certain number of people, a single vacancy can only hold a limited number of Hydrogen atoms. The study found that up to six Hydrogen atoms can crowd into the space around one vacancy.

2. Iron vs. Chromium: The "Velcro" Difference

The researchers compared how well Iron and Chromium hold onto these Hydrogen guests.

• Iron (Fe): Think of Iron's vacancy as a piece of light tape. It holds the Hydrogen, but it's not super sticky. The Hydrogen can still wiggle free relatively easily.
• Chromium (Cr): Think of Chromium's vacancy as super-strong Velcro. It grabs the Hydrogen much tighter. The study showed that Hydrogen is trapped more strongly in Chromium than in Iron. In fact, the "stickiness" (binding energy) is higher in Chromium, meaning it's harder for the Hydrogen to escape.

3. The "Crowded Room" Effect

As more Hydrogen atoms pile into the vacancy (up to six), the rules change.

• The Trend: Usually, as the room gets more crowded, it gets easier for the last person to leave because they are being pushed out by the others. The study confirmed that generally, the energy needed to escape (detrapping) goes down as more Hydrogen arrives.
• The Surprise: Previous studies suggested that the sixth Hydrogen atom in Iron would just fall out effortlessly (barrierless). However, this study found that even the sixth atom in Iron still has to fight a little bit to get out. It's not a free exit; there is still a small "door" it has to push through.

4. The Traffic Jam (Diffusion)

Finally, the researchers looked at the big picture: How fast does Hydrogen travel through the metal?

• The Result: When there are many vacancies (empty seats), the Hydrogen gets stuck more often. It's like a highway where cars keep getting pulled into side parking spots. The more parking spots (vacancies) there are, the slower the traffic moves.
• The Difference: This traffic jam is much worse in Chromium than in Iron. Because Chromium's "Velcro" is so strong, the Hydrogen gets stuck for longer, making the metal much less permeable to Hydrogen. In Iron, the Hydrogen moves faster, but it still slows down significantly if there are many vacancies.

Summary

This paper is essentially a detailed investigation into how "empty seats" in metal affect the movement of tiny Hydrogen atoms.

• Vacancies act as traps.
• Chromium is a much stronger trap than Iron.
• More vacancies mean slower movement for the Hydrogen.
• Even the last Hydrogen atom in an Iron vacancy has to work to escape, correcting some previous ideas that it would just fall out easily.

By understanding these tiny interactions, scientists can better predict how metals will behave in harsh environments, helping to prevent materials from becoming brittle and breaking.

Technical Summary

Technical Summary: Effect of Vacancies on Hydrogen Mobility and Trapping in Elemental Fe and Cr

Problem Statement
Hydrogen embrittlement (HE) significantly degrades the mechanical properties of metallic alloys, particularly those based on iron (Fe) and chromium (Cr), which are critical for nuclear, aerospace, and automotive applications. While hydrogen diffusion in pristine Fe has been extensively studied, the behavior of hydrogen in Cr, and specifically the role of vacancies in influencing diffusion and trapping in both elements, remains comparatively limited. Understanding how vacancy defects alter hydrogen mobility, trapping energetics, and detrapping mechanisms is essential for modeling hydrogen transport and predicting embrittlement in structural materials.

Methodology
The study employs a multiscale approach combining first-principles Density Functional Theory (DFT) and Kinetic Monte Carlo (kMC) simulations:

• DFT Calculations: Performed using the Vienna Ab initio Simulation Package (VASP) with the PBE-GGA functional and PAW potentials. 54-atom supercells were used for BCC Fe and Cr. The climbing-image nudged elastic band (CI-NEB) method determined migration and detrapping energy barriers. Zero-point energy (ZPE) corrections were calculated but found to have a negligible effect on binding and detrapping trends, and thus were excluded from the final barrier calculations.
• kMC Simulations: Used to model hydrogen diffusion over extended time and length scales. Hydrogen atoms were distributed on tetrahedral sites in 50×50×50 supercells. The model incorporated vacancy trapping with a defined 2 Å radius, allowing up to six hydrogen atoms per vacancy. H–H interactions were neglected due to high repulsion, and vacancy migration was ignored due to its significantly higher energy barrier compared to hydrogen. Simulations were averaged over 18 runs at various temperatures (300 K to 1100 K) and vacancy concentrations (40, 60, and 80 ppm).

Key Contributions and Results

1. Vacancy Formation and Binding Energies:

   • Hydrogen preferentially occupies tetrahedral sites in pristine Fe and Cr. Upon the introduction of a vacancy, hydrogen migrates toward the vacancy center, stabilizing near octahedral sites.
   • Up to six hydrogen atoms can occupy the vicinity of a single vacancy in both metals.
   • Binding Strength: Cr exhibits stronger hydrogen-vacancy binding than Fe. Bader charge analysis shows hydrogen gains more charge in Cr (0.35 e) than in Fe (0.30 e). Consequently, the vacancy formation energy decreases more significantly with hydrogen accumulation in Cr.
   • Charge Density: QTAIM analysis indicates higher electron density at the bond critical point for Fe-H compared to Cr-H, yet the overall trapping tendency is stronger in Cr due to the higher solute formation energy in Cr, driving hydrogen more strongly toward the vacancy.

2. Migration and Detrapping Barriers:

   • Migration: Hydrogen diffusion is slightly faster in Fe (0.05 eV) than in Cr (0.059 eV).
   • Detrapping Trends: Generally, detrapping energy decreases as hydrogen occupancy at the vacancy increases. However, a non-monotonic trend is observed where the fourth hydrogen atom has a lower detrapping energy than the fifth, attributed to the symmetric planar arrangement of the first four atoms versus the out-of-plane position of the fifth.
   • Sixth Hydrogen Discrepancy: Contrary to previous reports suggesting barrierless detrapping for the sixth hydrogen atom in Fe, this study finds a finite detrapping barrier (0.126 eV). In Cr, the sixth atom also exhibits a finite barrier (0.290 eV).

3. Diffusion Coefficients and Mobility:

   • Vacancy defects significantly reduce hydrogen mobility and increase the effective activation energy for diffusion.
   • Concentration Effect: In Fe, increasing vacancy concentration from 0 to 20 ppm reduces the diffusion coefficient at 300 K by approximately 67.8%. In Cr, the reduction is more severe, reaching approximately 97.1% for the same vacancy concentration increase.
   • Trapping Efficiency: Cr acts as a stronger trapping medium than Fe. At 300 K, the number of trapped hydrogen atoms is higher in Cr systems with equivalent vacancy concentrations.

Significance and Claims
The paper claims that the combined DFT–kMC framework provides a detailed, systematic understanding of hydrogen trapping and detrapping mechanisms in BCC metals. The study highlights distinct differences between Fe and Cr, specifically noting that Cr exhibits stronger trapping capabilities and a more pronounced suppression of hydrogen mobility in the presence of vacancies.

The authors emphasize that their findings correct previous assumptions regarding the sixth hydrogen atom in Fe, demonstrating that it does not undergo barrierless detrapping as some earlier studies using embedded atom method (EAM) potentials suggested. By quantifying the variation in activation energy and diffusion coefficients with vacancy concentration, the work offers critical insights into the mechanisms governing hydrogen transport and potential embrittlement in structural materials containing these elements.