skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion

Abstract

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. 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 withinmore » 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.« less

Authors:
 [1];  [2];  [1]
  1. Univ. of Illinois, Urbana-Champaign, IL (United States). Dept. of Materials Science and Engineering
  2. Univ. of New South Wales, Sydney, NSW (Australia). School of Mechanical and Manufacturing Engineering; Imperial College London, London (United Kingdom). Materials Dept.
Publication Date:
Research Org.:
Lawrence Berkeley National Laboratory-National Energy Research Scientific Computing Center (NERSC)
Sponsoring Org.:
USDOE
OSTI Identifier:
1530403
Resource Type:
Journal Article
Journal Name:
Physical Review Materials
Additional Journal Information:
Journal Volume: 3; Journal Issue: 3; Journal ID: ISSN 2475-9953
Country of Publication:
United States
Language:
English

Citation Formats

Jain, Abhinav C. P., Burr, Patrick A., and Trinkle, Dallas R. First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion. United States: N. p., 2019. Web. doi:10.1103/PhysRevMaterials.3.033402.
Jain, Abhinav C. P., Burr, Patrick A., & Trinkle, Dallas R. First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion. United States. doi:10.1103/PhysRevMaterials.3.033402.
Jain, Abhinav C. P., Burr, Patrick A., and Trinkle, Dallas R. Fri . "First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion". United States. doi:10.1103/PhysRevMaterials.3.033402.
@article{osti_1530403,
title = {First-principles calculations of solute transport in zirconium: Vacancy-mediated diffusion with metastable states and interstitial diffusion},
author = {Jain, Abhinav C. P. and Burr, Patrick A. and Trinkle, Dallas R.},
abstractNote = {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. 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.},
doi = {10.1103/PhysRevMaterials.3.033402},
journal = {Physical Review Materials},
issn = {2475-9953},
number = 3,
volume = 3,
place = {United States},
year = {2019},
month = {3}
}

Works referenced in this record:

A review of microstructure evolution in zirconium alloys during irradiation
journal, October 1988


Redistribution of the alloying elements during Zircaloy-2 oxidation
journal, September 1997


Tin clustering and precipitation in the oxide during autoclave corrosion of Zircaloy-2
journal, January 2015


Corrosion of Zirconium Alloys Used for Nuclear Fuel Cladding
journal, July 2015


The effect of Sn concentration on oxide texture and microstructure formation in zirconium alloys
journal, October 2015


Barrier oxide chemistry and hydrogen pick-up mechanisms in zirconium alloys
journal, January 2016


Atom probe tomography study of alloying element distributions in Zr alloys and their oxides
journal, November 2013


EELS and atom probe tomography study of the evolution of the metal/oxide interface during zirconium alloy oxidation
journal, July 2015


Light-element diffusion in Mg using first-principles calculations: Anisotropy and elastodiffusion
journal, August 2016


Exact Model of Vacancy-Mediated Solute Transport in Magnesium
journal, March 2017


Ab initio magnesium-solute transport database using exact diffusion theory
journal, May 2018


Performance of iron–chromium–aluminum alloy surface coatings on Zircaloy 2 under high-temperature steam and normal BWR operating conditions
journal, March 2016


First principles calculations of beryllium stability in zirconium surfaces
journal, January 2017


Ab initio atomic-scale determination of point-defect structure in hcp zirconium
journal, February 2005


Self-diffusion in the hexagonal structure of Zirconium and Hafnium: computer simulation studies
journal, December 2005


Atomistic calculation of point-defect diffusion anisotropy and irradiation growth in α -zirconium
journal, June 2007


Anisotropy of the Vacancy Migration in Ti, Zr and Hf Hexagonal Close-Packed Metals from First Principles
journal, November 2007


Issues in the ab-initio assessment of hcp transition metals self-diffusion
journal, August 2012


Self-interstitial configurations in hcp Zr: a first principles analysis
journal, February 2013


Point defect modeling in materials: Coupling ab initio and elasticity approaches
journal, October 2013

  • Varvenne, Céline; Bruneval, Fabien; Marinica, Mihai-Cosmin
  • Physical Review B, Vol. 88, Issue 13
  • DOI: 10.1103/PhysRevB.88.134102

Pressure effect on stabilities of self-Interstitials in HCP-Zirconium
journal, July 2014

  • Peng, Qing; Ji, Wei; Lian, Jie
  • Scientific Reports, Vol. 4, Issue 1
  • DOI: 10.1038/srep05735

Hydrogen and vacancy clustering in zirconium
journal, January 2016


Effect of alloying elements on the properties of Zr and the Zr–H system
journal, February 2014


First-principles investigation on diffusion mechanism of alloying elements in dilute Zr alloys
journal, August 2018


Anomalous energy pathway of vacancy migration and self-diffusion in hcp Ti
journal, June 2011


A first-principles study of the diffusion coefficients of alloying elements in dilute α-Ti alloys
journal, January 2016

  • Xu, W. W.; Shang, S. L.; Zhou, B. C.
  • Physical Chemistry Chemical Physics, Vol. 18, Issue 25
  • DOI: 10.1039/C6CP01899H

Diffusion calculations with two atomic models in h.c.p Zr-Nb diluted alloys
journal, August 2017


New Technique for Calculating Correlations for Vacancy Diffusion
journal, June 1962


Diffusion kinetics in dilute binary alloys with the h.c.p. crystal structure
journal, June 2014


Initial observation of grain boundary solute segregation in a zirconium alloy (ZIRLO) by three-dimensional atom probe
journal, August 2009


Enrichment of Fe and Ni at metal and oxide grain boundaries in corroded Zircaloy-2
journal, December 2012


The effect of matrix chemistry on dislocation evolution in an irradiated Zr alloy
journal, May 2017


Nano-scale chemical evolution in a proton-and neutron-irradiated Zr alloy
journal, April 2017


Redistribution of alloying elements in Zircaloy-2 after in-reactor exposure
journal, November 2014


51Cr diffusion in α-Zr single crystals
journal, November 1983


Chromium diffusion in α-zirconium, Zircaloy-2 and Zr-2.5Nb
journal, March 1993


Diffusion of 59 Fe in single-crystal α-Zr
journal, September 1988


Ultra-fast solute diffusion in α-Ti and α-Zr
journal, August 1972


Ni Diffusion in Zr and Zr Alloys
journal, January 1987


Anomalous self-diffusion in a single crystal of α-zirconium
journal, November 1984


Diffusion and vacancy properties of α-Zr
journal, July 1986


Self-diffusion in α-Zr single crystals
journal, May 1991


α-Zr self-diffusion anisotropy
journal, May 1995


From solid solution to cluster formation of Fe and Cr in α -Zr
journal, December 2015


First-principles appraisal of solute ultra-fast diffusion in hcp Zr and Ti
journal, March 2013


Ab initio approach to the effect of Fe on the diffusion in hcp Zr II: The energy barriers
journal, July 2009


Microstructure Evolution in Zr Alloys during Irradiation: Dose, Dose Rate, and Impurity Dependence
journal, January 2008

  • Griffiths, M.; Limback, Magnus; Kammenzind, Bruce
  • Journal of ASTM International, Vol. 5, Issue 1
  • DOI: 10.1520/JAI101146

Radiation-Induced Dislocation and Growth Behavior of Zirconium and Zirconium Alloys – a Review
journal, June 2013


The Kinetic Monte Carlo method: Foundation, implementation, and application
journal, July 2008

  • Battaile, Corbett C.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issue 41-42
  • DOI: 10.1016/j.cma.2008.03.010

Frequency factors and isotope effects in solid state rate processes
journal, January 1957


Generalized Gradient Approximation Made Simple
journal, October 1996

  • Perdew, John P.; Burke, Kieron; Ernzerhof, Matthias
  • Physical Review Letters, Vol. 77, Issue 18, p. 3865-3868
  • DOI: 10.1103/PhysRevLett.77.3865

Projector augmented-wave method
journal, December 1994


Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set
journal, October 1996


High-precision sampling for Brillouin-zone integration in metals
journal, August 1989


Special points for Brillouin-zone integrations
journal, June 1976

  • Monkhorst, Hendrik J.; Pack, James D.
  • Physical Review B, Vol. 13, Issue 12, p. 5188-5192
  • DOI: 10.1103/PhysRevB.13.5188

A ‘missing’ family of classical orthogonal polynomials
journal, January 2011


A climbing image nudged elastic band method for finding saddle points and minimum energy paths
journal, December 2000

  • Henkelman, Graeme; Uberuaga, Blas P.; Jónsson, Hannes
  • The Journal of Chemical Physics, Vol. 113, Issue 22, p. 9901-9904
  • DOI: 10.1063/1.1329672

Thermal equations of state of the α , β , and ω phases of zirconium
journal, May 2005


First-principles study of atomic hydrogen adsorption and initial hydrogenation of Zr(0001) surface
journal, January 2013

  • Zhang, Peng; Wang, Shuangxi; Zhao, Jian
  • Journal of Applied Physics, Vol. 113, Issue 1
  • DOI: 10.1063/1.4772675

Vacancy clustering in zirconium: An atomic-scale study
journal, October 2014


Elemental vacancy diffusion database from high-throughput first-principles calculations for fcc and hcp structures
journal, January 2014


The vacancy properties of α-Zr
journal, February 1981


Tracer diffusion in α-Zr
journal, April 1974


Beryllium diffusion in zirconium
journal, March 1976


Mechanism of Diffusion of Copper in Germanium
journal, November 1956


Vacancy Mediated and Interstitial Solute Transport in Zr from Density Functional Theory Calculations
dataset, January 2018

  • Jain, Abhinav C. P.; Burr, Patrick A.; Trinkle, Dallas R.
  • Materials Data Facility
  • DOI: 10.18126/M2R04K

First principles methods using CASTEP
journal, January 2005

  • Clark, Stewart J.; Segall, Matthew D.; Pickard, Chris J.
  • Zeitschrift für Kristallographie - Crystalline Materials, Vol. 220, Issue 5/6
  • DOI: 10.1524/zkri.220.5.567.65075

Soft self-consistent pseudopotentials in a generalized eigenvalue formalism
journal, April 1990


Representations of quasi-Newton matrices and their use in limited memory methods
journal, January 1994

  • Byrd, Richard H.; Nocedal, Jorge; Schnabel, Robert B.
  • Mathematical Programming, Vol. 63, Issue 1-3
  • DOI: 10.1007/BF01582063

Relaxation of Crystals with the Quasi-Newton Method
journal, February 1997

  • Pfrommer, Bernd G.; Côté, Michel; Louie, Steven G.
  • Journal of Computational Physics, Vol. 131, Issue 1
  • DOI: 10.1006/jcph.1996.5612

A generalized synchronous transit method for transition state location
journal, October 2003