Accelerating off-lattice kinetic Monte Carlo simulations to predict hydrogen vacancy-cluster interactions in $\alpha$-Fe

Abstract

We present an enhanced off-lattice kinetic Monte Carlo (OLKMC) model, based on a new method for tolerant classification of atomistic local-environments that is invariant under Euclidean-transformations and permutations of atoms. Our method ensures that environments within a norm-based tolerance are classified as equivalent. During OLKMC simulations, our method guarantees to elide the maximum number of redundant saddle-point searches in symmetrically equivalent local-environments. Hence, we are able to study the trapping/detrapping of hydrogen from up to five-vacancy clusters and simultaneously the effect hydrogen has on the diffusivity of these clusters. These processes occur at vastly different timescales at room temperature in body-centred cubic iron. We are able to predict the diffusion pathways of clusters/complexes without a priori assumptions of their mechanisms, not only reproducing previously reported mechanisms but also discovering new ones for larger complexes. We detail the hydrogen-induced changes in the clusters' diffusion mechanisms and find evidence that, in contrast to mono-vacancies, the introduction of hydrogen to larger clusters can increase their diffusivity. We compare the effective hydrogen diffusivity to Oriani's classical theory of trapping, finding general agreement and some evidence that hydrogen may not always be in equilibrium with traps, when the traps are mobile. Finally, we are able to compute the trapping atmosphere of meta-stable states surrounding non-point traps, opening new avenues to better understand and predict hydrogen embrittlement in complex alloys.

References in corpus (8)

• Hydrogen-enhanced local plasticity in aluminum: an ab initio study
• The Kinetic Activation-Relaxation Technique: A Powerful Off-lattice On-the-fly Kinetic Monte Carlo Algorithm
• Kinetic Activation Relaxation Technique
• Ab-initio molecular dynamics simulation of hydrogen diffusion in $\alpha$-iron
• Understanding long-time vacancy aggregation in iron: a kinetic activation-relaxation technique study
• A fully quantum mechanical calculation of the diffusivity of hydrogen in iron using the tight binding approximation and path integral theory
• New off-lattice Pattern Recognition Scheme for off-lattice kinetic Monte Carlo Simulations
• Extended Pattern Recognition Scheme for Self-learning Kinetic Monte Carlo (SLKMC-II) Simulations

Cited by in corpus (2)

• Multiscale modeling of hydrogen diffusion in iron considering the effect of dislocations
• Atomistic Simulations of H-Cu Vacancy Cosegregation and H Diffusion in Cu Grain Boundary