First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion
- Univ. of Illinois, Urbana-Champaign, IL (United States). Dept. of Materials Science and Engineering
- Univ. of New South Wales, Sydney, NSW (Australia). School of Mechanical and Manufacturing Engineering; Imperial College London, London (United Kingdom). Materials Dept.
Zirconium alloys are the most widely used nuclear fuel cladding materials for light water power reactors where irradiation damage causes solute redistribution, leading to degradation of alloy properties such as corrosion resistance. Designing radiation-tolerant zirconium alloys requires a thorough understanding of the atomic-scale transport behavior of the alloying elements in Zr. In this work, we perform density function theory calculations to investigate the diffusion of Sn, Cr, Fe, Be, Al, and Ni in the hexagonal close-packed (HCP) Zr matrix. We develop a methodology to accurately model the metastable vacancy states along the basal migration path, known to occur in group IV metals. We compute the vacancy-mediated solute diffusion coefficients and drag ratios using the kinetic Monte Carlo method and an analytic Green's function method—the agreement between the two validates our methodology. The computed diffusion coefficients of Sn and Al show good agreement with the experimental data and we expect these solutes to diffuse via the vacancy-mediated mechanism. We use a Green's function approach, parameterized with data from density functional theory calculations, to compute the interstitial diffusion coefficients of Cr, Fe, Be, and Ni in the HCP Zr lattice. The computed diffusion coefficients of Cr, Ni, and Be agree with the experimental measurements within one order of magnitude, while those of Fe are within two orders of magnitude of the experimental measurements. The drag ratios for Cr, Fe, Be, and Ni are positive up to 1235 K, which suggests that nonequilibrium vacancy fluxes could drag these solutes toward sinks such as dislocation loops and grain boundaries. We also compute the transport coefficients without including the metastable states, and using the eight- and thirteen-frequency model. Our results show significant differences in drag ratio for the eight- and thirteen-frequency model predictions compared with the Green's function methodology, but smaller errors in the solute diffusivity. Combining interstitial and vacancy-mediated diffusivities, we predict the unusual result that increased vacancy concentration slows down solute diffusivity, while a sufficiently high vacancy concentration can change the dominant mechanism to an accelerated vacancy-mediated diffusion.
- Research Organization:
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
- Sponsoring Organization:
- USDOE Office of Nuclear Energy (NE), Nuclear Energy University Program (NEUP); USDOE
- Grant/Contract Number:
- IRP-12-4728
- OSTI ID:
- 1530403
- Alternate ID(s):
- OSTI ID: 1546350
- Journal Information:
- Physical Review Materials, Vol. 3, Issue 3; ISSN 2475-9953
- Publisher:
- American Physical Society (APS)Copyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Atomic-scale modeling of the thermodynamic and kinetic properties of dilute alloys driven by forced atomic relocations
|
journal | December 2019 |
Atomic scale modeling of the thermodynamic and kinetic properties of dilute alloys driven by forced atomic relocations | text | January 2019 |
Similar Records
Predicting the temperature dependence of self-diffusion behavior in Ni-Cr alloys via molecular dynamics
A Simple Kinetic Model of Zircaloy Zr(Fe,Cr){sub 2} Precipitate Amorphization During Neutron Irradiation