-
Transferability of Zr-Zr interatomic potentials
Authors:
Oliver G. Nicholls,
Dillion Frost,
Vidur Tuli,
Jana Smutna,
Mark R. Wenman,
Patrick A. Burr
Abstract:
Tens of Zr inter-atomic potentials (force fields) have been developed to enable atomic-scale simulations of Zr alloys. These can provide critical insight in the in-reactor behaviour of nuclear fuel cladding and structural components exposed, but the results are strongly sensitive to the choice of potential. We provide a comprehensive comparison of 13 popular Zr potentials, and assess their ability…
▽ More
Tens of Zr inter-atomic potentials (force fields) have been developed to enable atomic-scale simulations of Zr alloys. These can provide critical insight in the in-reactor behaviour of nuclear fuel cladding and structural components exposed, but the results are strongly sensitive to the choice of potential. We provide a comprehensive comparison of 13 popular Zr potentials, and assess their ability to reproduce key physical, mechanical, structural and thermodynamic properties of Zr. We assess the lattice parameters, thermal expansion, melting point, volume-energy response, allotropic phase stability, elastic properties, and point defect energies, and compare them to experimental and ab-initio values. No potential was found to outperform all others on all aspects, but for every metric considered here, at least one potential was found to provide reliable results. Older embedded-atom method (EAM) potentials tend to excel in 2-3 metrics each, but at the cost of poorer transferability. The two highest-performing potentials overall, with complementary strengths and weaknesses, were the 2021 angular-dependent potential of Smirnova and Starikov (Comp. Mater. Sci. 197, 110581) and the 2019 embedded-atom method potential of Wimmer et al (J. Nucl. Mater. 532, 152055). All potentials trained through machine learning algorithms proved to have lower overall accuracy, and less transferability, than simpler and computationally faster potentials available. Point defect structures and energies is where the greatest divergence and least accuracy is observed. We created maps that will help modellers select the most suitable potential for a specific application, and which may help identify areas of improvement in future potentials.
△ Less
Submitted 17 May, 2023; v1 submitted 7 July, 2022;
originally announced July 2022.
-
Structural and phase evolution in U$_3$Si$_2$ during steam corrosion
Authors:
Jiatu Liu,
Patrick A. Burr,
Joshua T. White,
Vanessa K. Peterson,
Pranesh Dayal,
Christopher Baldwin,
Deborah Wakeham,
Daniel J. Gregg,
Elizabeth S. Sooby,
Edward G. Obbard
Abstract:
U$_3$Si$_2$ nuclear fuel is corroded in deuterated steam with in situ neutron diffraction. Density functional theory is coupled with rigorous thermodynamic description of the hydride including gas/solid entropy contributions. H absorbs in the 2$b$ interstitial site of U$_3$Si$_2$H$_x$ and moves to 8$j$ for $x\ge 0.5$. Hydriding forces lattice expansion and change in $a/c$ ratio linked to site pref…
▽ More
U$_3$Si$_2$ nuclear fuel is corroded in deuterated steam with in situ neutron diffraction. Density functional theory is coupled with rigorous thermodynamic description of the hydride including gas/solid entropy contributions. H absorbs in the 2$b$ interstitial site of U$_3$Si$_2$H$_x$ and moves to 8$j$ for $x\ge 0.5$. Hydriding forces lattice expansion and change in $a/c$ ratio linked to site preference. Rietveld refinement tracks the corrosion reactions at 350-500 °C and preference for the 8j site. Above 375 °C, formation of UO$_2$, U$_3$Si$_5$ and USi$_3$ take place in the grain boundaries and bulk. Hydriding occurs in bulk and precedes other reactions.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
Importance of elastic finite-size effects: neutral defects in ionic compounds
Authors:
P. A. Burr,
M. W. D. Cooper
Abstract:
Small system sizes are a well known source of error in DFT calculations, yet computational constraints frequently dictate the use of small supercells, often as small as 96 atoms in oxides and compound semiconductors. In ionic compounds, electrostatic finite size effects have been well characterised, but self-interaction of charge neutral defects is often discounted or assumed to follow an asymptot…
▽ More
Small system sizes are a well known source of error in DFT calculations, yet computational constraints frequently dictate the use of small supercells, often as small as 96 atoms in oxides and compound semiconductors. In ionic compounds, electrostatic finite size effects have been well characterised, but self-interaction of charge neutral defects is often discounted or assumed to follow an asymptotic behaviour and thus easily corrected with linear elastic theory. Here we show that elastic effects are also important in the description of defects in ionic compounds and can lead to qualitatively incorrect conclusions if inadequately small supercells are used; moreover, the spurious self-interaction does not follow the behaviour predicted by linear elastic theory. Considering the exemplar cases of metal oxides with fluorite structure, we show that numerous previous studies, employing 96-atom supercells, misidentify the ground state structure of (charge neutral) Schottky defects. We show that the error is eliminated by employing larger cells (324, 768 and 1500 atoms), and careful analysis determines that elastic effects, not electrostatic, are responsible. The spurious self-interaction was also observed in non-oxide ionic compounds and irrespective of the computational method used, thereby resolving long standing discrepancies between DFT and force-field methods, previously attributed to the level of theory. The surprising magnitude of the elastic effects are a cautionary tale for defect calculations in ionic materials, particularly when employing computationally expensive methods or when modelling large defect clusters. We propose two computationally practicable methods to test the magnitude of the elastic self-interaction in any ionic system.
△ Less
Submitted 6 September, 2017;
originally announced September 2017.
-
Defect processes in Be$_{12}$X Beryllides
Authors:
Matthew L Jackson,
Patrick A Burr,
Robin W Grimes
Abstract:
The stability of intrinsic point defects in Be$_{12}$X intermetallics (where X = Ti, V, Mo or W) are predicted using density functional theory simulations and discussed with respect to fusion energy applications. Schottky disorder is found to be the lowest energy complete disorder process, closely matched by Be Frenkel disorder in the cases of Be$_{12}$V and Be$_{12}$Ti. Antitisite and X Frenkel d…
▽ More
The stability of intrinsic point defects in Be$_{12}$X intermetallics (where X = Ti, V, Mo or W) are predicted using density functional theory simulations and discussed with respect to fusion energy applications. Schottky disorder is found to be the lowest energy complete disorder process, closely matched by Be Frenkel disorder in the cases of Be$_{12}$V and Be$_{12}$Ti. Antitisite and X Frenkel disorder are of significantly higher energy. Small clusters of point defects including Be divacancies, Be di-interstitials and accommodation of the X species on two Be sites were considered. Some di-interstitial, divacancy and X$_{2Be}$ combinations exhibit negative binding enthalpy (i.e. clustering is favourable), although this is orientationally dependent. None of the Be$_{12}$X intermetallics are predicted to exhibit significant non-stoichiometry, ruling out non-stoichiometry as a mechanism for accommodating Be depletion due to neutron transmutation.
△ Less
Submitted 6 April, 2017;
originally announced April 2017.
-
Experimental and DFT investigation of (Cr,Ti)3AlC2 MAX phases stability
Authors:
Patrick A. Burr,
Denis Horlait,
William E. Lee
Abstract:
Using a synergistic combination of experimental and computational methods, we shed light on the unusual solubility of (Cr,Ti)3AlC2 MAX phase, showing that it may accommodate Cr only at very low concentrations (<2at%) or at the exact Cr/(Cr+Ti) ratio of 2/3, even when the ratio of reactants is far from this stoichiometry (1/2 < Cr/(Cr+Ti) < 5/6). In both phases Cr exclusively occupies the 4f sites,…
▽ More
Using a synergistic combination of experimental and computational methods, we shed light on the unusual solubility of (Cr,Ti)3AlC2 MAX phase, showing that it may accommodate Cr only at very low concentrations (<2at%) or at the exact Cr/(Cr+Ti) ratio of 2/3, even when the ratio of reactants is far from this stoichiometry (1/2 < Cr/(Cr+Ti) < 5/6). In both phases Cr exclusively occupies the 4f sites, bridging carbide layers with the Al layer. Despite this, the peculiar stability of (Cr2/3Ti1/3)3AlC2 is attributed to the formation of strong, spin-polarised Cr-C bonds, which result in volume reduction and marked increase in c/a ratio.
△ Less
Submitted 8 August, 2016; v1 submitted 20 July, 2016;
originally announced July 2016.
-
Solubility and partitioning of impurities in Be alloys
Authors:
Patrick A. Burr,
Simon C. Middleburgh,
Robin W. Grimes
Abstract:
The most energetically favourable accommodation processes for common impurities and alloying elements in Be metal and Be-Fe-Al intermetallics were investigated using atomic scale simulations. Fe additions, combined with suitable heat treatments, may scavange Al and Si through their incorporation into the FeBe5 intermetallic. In the absence of Fe, Al and Si will not be associated with Be metal. Li…
▽ More
The most energetically favourable accommodation processes for common impurities and alloying elements in Be metal and Be-Fe-Al intermetallics were investigated using atomic scale simulations. Fe additions, combined with suitable heat treatments, may scavange Al and Si through their incorporation into the FeBe5 intermetallic. In the absence of Fe, Al and Si will not be associated with Be metal. Li and Mg are also not soluble, but may react with other impurities if present (such as Al or H). Mg may also form the MgBe13 intermetallic phase under certain conditions. He and H exhibit negligible solubility in all phases investigated and whilst He will tend to form bubbles, H can precipitate as BeH2. Similarly, C additions will form the stable compound Be2C. Finally, oxygen exhibits a strong affinity to Be, exhibiting both some degree of solubility in all phases considered here (though especially metallic Be) and a highly favourable energy of formation for BeO.
△ Less
Submitted 6 July, 2016;
originally announced July 2016.
-
Resolving the structure of TiBe$_{12}$
Authors:
M. L. Jackson,
P. A. Burr,
R. W. Grimes
Abstract:
There has been considerable controversy regarding the structure of TiBe$_{12}$, which is variously reported as hexagonal and tetragonal. Lattice dynamics simulations based on density functional theory show the tetragonal phase space group $I4/mmm$ to be more stable over all temperatures, while the hexagonal phase exhibits an imaginary phonon mode, which, if followed, would lead to the cell adoptin…
▽ More
There has been considerable controversy regarding the structure of TiBe$_{12}$, which is variously reported as hexagonal and tetragonal. Lattice dynamics simulations based on density functional theory show the tetragonal phase space group $I4/mmm$ to be more stable over all temperatures, while the hexagonal phase exhibits an imaginary phonon mode, which, if followed, would lead to the cell adopting the tetragonal structure. We then report the predicted ground state elastic constants and temperature dependence of the bulk modulus and thermal expansion for the tetragonal phase.
△ Less
Submitted 14 March, 2016;
originally announced March 2016.
-
Accommodation of Tin in Tetragonal ZrO2
Authors:
B. D. C. Bell,
S. T. Murphy,
P. A. Burr,
R. W. Grimes,
M. R. Wenman
Abstract:
Atomic scale computer simulations using density functional theory were used to investigate the behaviour of tin in the tetragonal phase oxide layer on Zr-based alloys. The $Sn_{Zr}^{\times}$ site defect was shown to be dominant across most oxygen partial pressures, with $Sn_{Zr}^{"}$ charge compensated by fully charged oxygen vacancies occurring at partial pressures below $10^{-31}$ atm. Insertion…
▽ More
Atomic scale computer simulations using density functional theory were used to investigate the behaviour of tin in the tetragonal phase oxide layer on Zr-based alloys. The $Sn_{Zr}^{\times}$ site defect was shown to be dominant across most oxygen partial pressures, with $Sn_{Zr}^{"}$ charge compensated by fully charged oxygen vacancies occurring at partial pressures below $10^{-31}$ atm. Insertion of additional positive charge into the system was shown to significantly increase the critical partial pressure at which $Sn_{Zr}^{"}$ is stable. Recently developed low-Sn nuclear fuel cladding alloys have demonstrated an improved corrosion resistance and a delayed transition compared to Sn-containing alloys, such as Zircaloy-4. The interaction between the positive charge and the tin defect is discussed in the context of alloying additions, such as niobium and their influence on corrosion of cladding alloys.
△ Less
Submitted 24 February, 2015;
originally announced February 2015.
-
From solid solution to cluster formation of Fe and Cr in $α$-Zr
Authors:
P. A. Burr,
M. R. Wenman,
B. Gault,
M. P. Moody,
M. Ivermark,
M. J. D. Rushton,
M. Preuss,
L. Edwards,
R. W. Grimes
Abstract:
To understand the mechanisms by which Fe and Cr additions increase the corrosion rate of irradiated Zr alloys, a combination of experimental (atom probe tomography, x-ray diffraction and thermoelectric power measurements) and modelling (density functional theory) techniques are employed to investigate the non-equilibrium solubility and clustering of Fe and Cr in binary Zr alloys. Cr occupies both…
▽ More
To understand the mechanisms by which Fe and Cr additions increase the corrosion rate of irradiated Zr alloys, a combination of experimental (atom probe tomography, x-ray diffraction and thermoelectric power measurements) and modelling (density functional theory) techniques are employed to investigate the non-equilibrium solubility and clustering of Fe and Cr in binary Zr alloys. Cr occupies both interstitial and substitutional sites in the α-Zr lattice, Fe favours interstitial sites, and a low-symmetry site that was not previously modelled is found to be the most favourable for Fe. Lattice expansion as a function of alloying concentration (in the dilute regime) is strongly anisotropic for Fe additions, expanding the $c$-axis while contracting the $a$-axis. Defect clusters are observed at higher solution concentrations, which induce a smaller amount of lattice strain compared to the dilute defects. In the presence of a Zr vacancy, all two-atom clusters are more soluble than individual point defects and as many as four Fe or three Cr atoms could be accommodated in a single Zr vacancy. The Zr vacancy is critical for the increased solubility of defect clusters, the implications for irradiation induced microstructure changes in Zr alloys are discussed.
△ Less
Submitted 10 October, 2015; v1 submitted 27 January, 2015;
originally announced January 2015.
-
Crystal structure, thermodynamics, magnetics and disorder properties of Be-Fe-Al intermetallics
Authors:
Patrick A. Burr,
Simon C. Middleburgh,
Robin W. Grimes
Abstract:
The elastic and magnetic properties, thermodynamical stability, deviation from stoichiometry and order/disorder transformations of phases that are relevant to Be alloys were investigated using density functional theory simulations coupled with phonon density of states calculations to capture temperature effects. A novel structure and composition were identified for the Be-Fe binary ε phase. In abs…
▽ More
The elastic and magnetic properties, thermodynamical stability, deviation from stoichiometry and order/disorder transformations of phases that are relevant to Be alloys were investigated using density functional theory simulations coupled with phonon density of states calculations to capture temperature effects. A novel structure and composition were identified for the Be-Fe binary ε phase. In absence of Al, FeBe_5 is predicted to form at equilibrium above ~ 1250 K, while the ε phase is stable only below ~ 1650 K, and FeBe_2 is stable at all temperatures below melting. Small additions of Al are found to stabilise FeBe_5 over FeBe_2 and ε, while at high Al content, AlFeBe_4 is predicted to form. Deviations from stoichiometric compositions are also considered and found to be important in the case of FeBe_5 and ε. The propensity for disordered vs ordered structures is also important for AlFeBe_4 (which exhibits complete Al-Fe disordered at all temperatures) and FeBe_5 (which exhibits an order-disorder transition at ~ 950 K).
△ Less
Submitted 18 March, 2015; v1 submitted 18 July, 2014;
originally announced July 2014.
-
The thermodynamics of hydride precipitation: the importance of entropy, enthalpy and disorder
Authors:
S. C. Lumley,
R. W. Grimes,
S. T. Murphy,
P. A. Burr,
A. Chroneos,
P. R. Chard-Tucke,
M. R. Wenman
Abstract:
The thermodynamics of H/α-Zr solid solution and zirconium hydride phases were studied using density functional theory. Disorder in ζ, γ and δ hydrides and solid solutions were modelled using a statistically significant number of randomly generated structures in combination with special quasi-random structures and solid solutions with a range of concentrations. This is used in conjunction with a ca…
▽ More
The thermodynamics of H/α-Zr solid solution and zirconium hydride phases were studied using density functional theory. Disorder in ζ, γ and δ hydrides and solid solutions were modelled using a statistically significant number of randomly generated structures in combination with special quasi-random structures and solid solutions with a range of concentrations. This is used in conjunction with a calculation of thermodynamic parameters of the system, including the temperature dependent sensible enthalpy, configurational entropy and vibrational entropy of the different crystals in the system, developed from phonon density of states. It was found that precipitation of hydrides is not thermodynamically favourable for Zr-H solid solutions containing less than 300 ppm H, suggesting that a mechanism must cause local concentration of H atoms to a greater amount than found globally in experimental samples containing hydrides. Temperature drives the reaction in the direction of solution, primarily due to entropic effects. Generally, γ hydride is the most stable phase, although it is very close in energy to the δ-phase. The sensible enthalpy of precipitation assists in stabilising HCP hydrides, and the configurational entropy change during precipitation favours FCC hydrides. None of the thermodynamic contributions are found to be negligible in driving precipitation.
△ Less
Submitted 22 January, 2014;
originally announced January 2014.
-
Hydrogen solubility in zirconium intermetallic second phase particles
Authors:
P. A. Burr,
S. T. Murphy,
S. C. Lumley,
M. R. Wenman,
R. W. Grimes
Abstract:
The enthalpies of solution of H in Zr binary intermetallic compounds formed with Cu, Cr, Fe, Mo, Ni, Nb, Sn and V were calculated by means of density functional theory simulations and compared to that of H in α-Zr. It is predicted that all Zr-rich phases (formed with Cu, Fe, Ni and Sn), and those phases formed with Nb and V, offer lower energy, more stable sites for H than α-Zr. Conversely, Mo and…
▽ More
The enthalpies of solution of H in Zr binary intermetallic compounds formed with Cu, Cr, Fe, Mo, Ni, Nb, Sn and V were calculated by means of density functional theory simulations and compared to that of H in α-Zr. It is predicted that all Zr-rich phases (formed with Cu, Fe, Ni and Sn), and those phases formed with Nb and V, offer lower energy, more stable sites for H than α-Zr. Conversely, Mo and Cr containing phases do not provide preferential solution sites for H. In all cases the most stable site for H are those that offer the highest coordination fraction of Zr atoms. Often these are four Zr tetrahedra but not always. Implications with respect to H-trapping properties of commonly observed ternary phases such as Zr(Cr,Fe)2, Zr2(Fe,Ni) and Zr(Nb,Fe)2 are also discussed.
△ Less
Submitted 29 July, 2013;
originally announced July 2013.
-
Hydrogen accommodation in Zr second phase particles: Implications for H pick-up and hydriding of Zircaloy-2 and Zircaloy-4
Authors:
P. A. Burr,
S. T. Murphy,
S. C. Lumley,
M. R. Wenman,
R. W. Grimes
Abstract:
Ab-initio computer simulations have been used to predict the energies associated with the accommodation of H atoms at interstitial sites in α, β-Zr and Zr.M intermetallics formed with common alloying additions (M = Cr, Fe, Ni). Intermetallics that relate to the Zr2(Ni,Fe) second phase particles (SPPs) found in Zircaloy-2 exhibit favourable solution enthalpies for H. The intermetallic phases that r…
▽ More
Ab-initio computer simulations have been used to predict the energies associated with the accommodation of H atoms at interstitial sites in α, β-Zr and Zr.M intermetallics formed with common alloying additions (M = Cr, Fe, Ni). Intermetallics that relate to the Zr2(Ni,Fe) second phase particles (SPPs) found in Zircaloy-2 exhibit favourable solution enthalpies for H. The intermetallic phases that relate to the Zr(Cr,Fe)2 SPPs, found predominantly in Zircaloy-4, do not offer favourable sites for interstitial H. It is proposed that Zr(Cr,Fe)2 particles may act as bridges for the migration of H through the oxide layer, whilst the Zr2(Ni,Fe)-type particles will trap the migrating H until these are dissolved or fully oxidised.
△ Less
Submitted 29 July, 2013;
originally announced July 2013.