Effect of structure and composition on the electronic excitation induced amorphization of La 2 Ti 2−x Zr x O 7 ceramics
Michel Sassi, Tiffany Kaspar, Kevin M. Rosso, Steven R. Spurgeon
EEffect of structure and composition on the electronic excitation induced amorphization of La Ti Zr x O ceramics Michel Sassi , Tiffany Kaspar , Kevin M. Rosso , and Steven R. Spurgeon Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, USA *[email protected]
Abstract:
Understanding the response of ceramics operating in extreme environments is of interest for a variety of applications.
Ab initio molecular dynamic simulations have been used to investigate the effect of structure and B -site (=Ti, Zr) cation composition of lanthanum-based oxides (La B O ) on electronic-excitation-induced amorphization. We find that the amorphous transition in monoclinic layered perovskite La Ti O occurs for a lower degree of electronic excitation than for cubic pyrochlore La Zr O . While in each case the formation of O -like molecules drives the structure to an amorphous state, an analysis of the polyhedral connection network reveals that the rotation of TiO octahedra in the monoclinic phase can promote such molecule formation, while such octahedral rotation is not possible in the cubic phase. However, once the symmetry of the cubic structure is broken by substituting Ti for Zr, it becomes less resistant to amorphization. A compound made of 50% Ti and 50% Zr (La TiZrO ) is found to be more resistant in the monoclinic than in the cubic phase, which may be related to the lower bandgap of the cubic phase. These results illustrate the complex interplay of structure and composition that give rise to the radiation resistance of these important functional materials. ntroduction: Pyrochlores, with the general ideal formula A B O , represent an important class of complex oxide ceramics. With more than 500 different combinations of A - and B -site cations , pyrochlores exhibit a range of physical, chemical and electronic properties, making them useful for a variety of applications, including nuclear fuel and waste forms, solar energy conversion, electronics, and catalysis . While the A -site cation is generally occupied by tri- and tetravalent actinides and lanthanides, the B -site is usually occupied by transition or post-transition metals. The stability and structural flexibility of pyrochlores is governed by the radius ratio between the A - and B -site cations . A B O oxides with a r A /r B ratio ranging from 1.46 to 1.78 are isometric ( 𝐹𝑑3$𝑚 , Z=8, a=0.9-1.2 nm) and adopt the pyrochlore structure, which is related to ideal fluorite (AX ) but with ordered cation sites and anion vacancies. If the cation radius ratio is above 1.78, then a monoclinic layered perovskite structure ( 𝑃2 ( ) is favored, while for a ratio below 1.46, a distorted fluorite structure ( 𝐹𝑑3$𝑚 ) is preferred. Pyrochlores used in the nuclear fuel cycle will be exposed to extreme conditions potentially affecting their long-term performance and reliability. Pyrochlores are potential waste form candidates for the immobilization of plutonium and other minor actinides from spent nuclear fuel because significant amounts of radioisotopes can be incorporated into the A -site of the pyrochlore structure. To better understand the structural changes that can occur under irradiation, pyrochlores have been the subject of many experiments that focused on characterizing their response to electron , ion , and photon (pulsed laser) irradiation in terms of radiation-induced defect generation, damage processes, microstructural evolution, and phase transformations. Many studies were specifically designed to investigate the effect of α-decay in various pyrochlore compositions using keV-GeV ion-beam irradiation; in this range, ion stopping is mediated predominantly by interactions with the target electrons, inducing excitations and ionizations along the ion path. These experiments have highlighted the importance of structure, bond-type, and electronic configuration on their radiation stability, particularly for the titanate ( A Ti O ) and zirconate ( A Zr O ) pyrochlores . While pyrochlore structures display a wide range of behaviors in response to ion-beam irradiation as a function of composition, it is generally found that titanate pyrochlores, such as d Ti O , are less resistant to electronic excitations compared to their zirconate counterparts, such as Gd Zr O . Analysis of the thermal spike following 119 MeV U irradiation shows that Gd Ti O is readily amorphized, while Gd Zr O is transformed into a radiation-resistant anion-deficient fluorite structure. For other zirconate pyrochlores ( A Zr O ), a more complex response was observed, in which both a pyrochlore-fluorite phase transformation and amorphization occur. B -site mixed Gd Ti x Zr x O compounds showed a systematic increase in the resistance to ion-beam induced amorphization with increasing Zr content . This result highlights the important effect of the B -site cation on the response of pyrochlore materials to electronic excitations. Although zirconate pyrochlores are generally considered to be radiation-resistant, a systematic ion-beam study of A Zr O pyrochlores shows that the susceptibility to amorphization by electronic excitation increases with increasing cation radius ratio r A /r B15 . While theoretical studies investigating the effect of atomic collision processes, defect generation and phase transformation have been performed to understand the mechanisms leading to amorphization, the role of electronic excitation in amorphization has been investigated less. The simulation of electronic excitation effects is a challenging task, in which we must account for the dynamics of the electron-hole recombination . In order to mimic the effects of electronic excitation and ionization, a simplified methodology, consisting of removing several electrons from the higher valence band states, has been proposed and used in previous studies to investigate electronic excitations induced amorphization of Ge-Sb-Te alloys and titanate pyrochlores . It is worth noting that the effect of electronic excitation described by this method, and the resulting mechanisms leading to amorphization, are different from the intense electronic excitation induced by swift heavy ion irradiation, where local melting from the thermal spike leads to a quenched melt structure . Nevertheless, these simulations of electronic excitation provide guidance for intense electron and pulsed laser irradiation studies, and help interpret trends in the susceptibility to electronic excitation-induced amorphization. In this study, we report ab initio molecular dynamics simulations that focus on determining how the transition temperature from crystalline to amorphous phase is affected by the structure, composition and electronic excitation concentration in lanthanum-based compounds. In that regard, La Ti O and La Zr O are ideal materials to investigate since they espectively occur in monoclinic layered perovskite and cubic pyrochlore structures. Therefore, mixing the B -site cation modifies not only the composition of the material, but also its structure. In this study we also consider intermediate monoclinic La Ti Zr x O and cubic La Zr Ti x O compounds to follow the evolution of the amorphization temperature as function of the B -site composition. We find that differences in the octahedral connectivity of the cubic and monoclinic phases significantly influence the potential for defect generation and relative phase stability, with important implications for the design of nuclear energy systems. omputational Details: In order to investigate the effect of structure and composition on electronic excitation induced amorphization, several compounds with variable B -site cation composition have been created by using the special quasirandom structure (SQS) generation code available in the ATAT toolkit . Two series of SQS intermediate structures were created. One series uses the monoclinic layered perovskite La Ti O (LTO) crystal as its starting point and substitutes Zr for Ti, giving rise to La Ti Zr x O compounds (labeled Z-LTO) with a monoclinic lattice. The other series uses the cubic pyrochlore La Zr O (LZO) crystal as its starting point and substitutes Ti for Zr, leading to La Zr Ti x O compounds (labeled T-LZO) with a cubic lattice. For consistency with the end-member LTO and LZO crystals, only SQS structures containing a total of 88 atoms were generated, such that all of the structures possess 16 La sites, 16 B -sites occupied by either pure Ti, pure Zr, or a mix of both species, and 56 O atoms. Density functional theory (DFT) calculations were performed on all the structures using the VASP package . The exchange-correlation functional used the generalized gradient approximation (GGA) as parametrized by Perdew and Wang . All of the calculations accounted for spin-polarization, used the Vosko-Wilk-Nusair local density approximation scheme , and a cutoff energy for the projector augmented wave pseudo-potential of 400 eV. Prior to performing ab initio molecular dynamics simulations, the atomic coordinates and lattice parameters of each investigated structure were fully relaxed using a convergence criterion of 10 -5 eV/cell for the total energy and 10 -4 eV/Å for the force components. A 4×4×2 and 4×4×4 Monkhorst-Pack k -point sampling of the Brillouin zone was used for the monoclinic and cubic structures respectively. To mimic the effect of electronic excitation, several electrons have been removed from the higher valence band states and a jellium background has been used to compensate for the loss of charge. A similar procedure has been used previously to investigate the role of electronic excitation in the amorphization of Ge-Sb-Te alloys and titanate pyrochlores . The electronic excitation concentration has been defined as the percentage ratio of removed electrons with respect to the total number of electrons in the system. For the monoclinic LTO and cubic LZO crystals, the total number of electrons for a cell containing 88 atoms was 1728 and 2000, espectively, and the number of removed electrons varied from 7 to 35 for LTO and 8 to 40 for LZO, such that the excitation concentration investigated in the structures varied from 0.4% to 2% by steps of 0.4%. The electronic excitation concentration in the intermediate SQS structures followed a similar calculation. To determine the temperature at which the crystalline to amorphous transition occurs for each structure and electronic excitation conditions, ab initio molecular dynamic simulations have been performed at the Γ-point of the Brillouin zone in the NVT ensemble with a Nosé-Hoover thermostat for temperatures ranging from 200 K to 1000 K by step of 100 K. For each case, a three-step procedure was followed for a total simulation time of 27 ps with a time step of 2 fs. First, 6 ps were used to thermalize the crystal at the given temperature. Second, 15 ps were given to allow structural modifications associated with electronic excitation. This time lapse was tested in each structure and proven long enough to allow the system to reach an equilibrium state. Third, 6 ps were used to allow electron-hole recombination. esults and Discussion: The monoclinic LTO structure is a layered perovskite crystal containing TiO octahedra linearly connected, as shown in Fig. 1(a). In contrast, the cubic pyrochlore LZO crystal contains ZrO octahedra inter-connected by corners in a zigzag fashion. A comparison of the two crystals in their respective ground state structure indicates that the monoclinic LTO is 0.325 eV/atom less energetically favorable than the cubic LZO, as shown in Fig. 1(b). In the case where a monoclinic perovskite structure is made of 100% Zr atoms ( i.e. , monoclinic LZO) it would be 0.039 eV/atom less favorable than its ground state cubic LZO structure. Similarly, if a cubic pyrochlore crystal was made of 100% Ti atoms ( i.e. , cubic LTO), its energy would be 0.055 eV/atom less favorable than that of the monoclinic LTO crystal. The general trend shown in Fig. 1(b) is that the incorporation of Zr in both crystal structures leads to more energetically favorable phases. For both the monoclinic and cubic crystals, SQS structures have been generated to find compositionally intermediate structures with the most appropriate [Ti, Zr] distributions that mimic, for a small number of atoms, the first few, physically most relevant radial correlation functions of a perfectly random structure. Using this method enables us to investigate the effect of progressive [Ti, Zr] mixing on electronic excitation induced amorphization. Among the SQS structures generated for 25%, 50% and 75% of Zr (calculated from the B -site only), only those with the lowest energy have been reported in Fig. 1(b). The energy diagram shows that the phase transition from monoclinic perovskite to cubic pyrochlore should occur for 53% of Zr in the monoclinic LTO structure. It is worth noting that the structures at 50% [Ti, Zr] mixing have very similar energies but different crystal structures. In the following, we have focused our investigations on the most energetically stable structures at any given composition, mainly by considering only a monoclinic layered perovskite phase for compounds with a Zr fraction ≤50% and a cubic pyrochlore phase for compounds with a Zr fraction ≥50%. At a Zr fraction of 50% both phases have been investigated. Figure 1: (a) Representation of the monoclinic layered perovskite LTO and cubic pyrochlore LZO structures. (b) Evolution of the crystal energy as function of Zr fraction in the monoclinic and cubic lattices. To determine how the crystal structure, composition and electronic excitation concentration affect the transition temperature from crystalline to amorphous phase, the radial distribution function (RDF) has been calculated at the end of the simulation time at each temperature and for each structure investigated. A material is considered amorphous when the short-range order in its RDF is maintained but the long-range order is lost. Fig. 2(a) shows the effect of various electronic excitation concentrations on the RDF of cubic LZO at 300 K. The eatures of the RDF remain similar until an excitation concentration of 1.2%, at which point bonds between oxygen atoms start to form. The appearance of O -like molecules is associated with a modification of the intermediate and long-range ordering, suggesting the beginning of amorphization. At higher electronic excitation concentrations, the degree of amorphization increases with the number of O—O bonds formed. It is worth noting that a similar behavior has been obtained by Xiao et al. during the simulation of the effect of electronic excitation in several titanate pyrochlores. Under the influence of electronic excitation, the formation of O -like molecules in these materials is due to the electrons removed from the system, which mainly originate from oxygen 2 p valence orbitals located near the Fermi level. The anion disorder resulting from the formation of O molecules further induces cation displacements, together driving the structure toward an amorphous phase. Fig. 2(b) shows the evolution of the RDF features at different temperatures for an excitation concentration fixed at 0.8%. The RDFs show that the formation of O—O bonds starts at 600 K, which coincides with the beginning of amorphization. Performing similar calculations for other excitation concentrations and [Zr, Ti] mixing allows us to investigate how the B -site composition affect the transition temperature from crystalline to amorphous phase, and its interplay with electronic excitation concentration in the cubic pyrochlore lattice. For each composition investigated, the general trend, shown in Fig. 2(c), is a lowering of the transition temperature as the excitation concentration increases such that each cubic compound amorphizes at room temperature for 1.2% of excitation concentration. While cubic LZO is found to be the most resistant composition to electronic excitation, a compound made of 50% Zr (i.e. La ZrTiO ) is the easiest to amorphize as suggested by the value of its transition temperature, which is always the lowest across the range of electronic excitation concentration investigated. Interestingly, a similar trend with respect to B -site composition has been observed experimentally for the Gd (Zr x Ti ) O pyrochlore family, with the pure Gd Zr O phase being more resistant to electronic excitation than a phase made of 50% Zr fraction (i.e. Gd ZrTiO ) . The evolution of the transition temperature as function of Zr fraction for 0.4%, 0.8% and 1.2% of excitation concentration is shown in Fig. 2(d). While for 0.4% of excitation concentration, the transition temperature is highly dependent on the B -site cation composition, or 1.2% of excitation concentration the amorphization temperature starts to be composition-independent. Figure 2: Evolution of the radial distribution function (g(r)) of cubic LZO as function of (a) the electronic excitation concentration at 300 K and (b) the temperature for a fixed 0.8% of electronic excitation concentration. (c) Evolution of crystalline-to-amorphous transition temperature in cubic T-LZO compounds as function of the excitation concentration for 50%, 75% and 100% Zr fraction. (d) Evolution of crystalline-to-amorphous transition temperature as function of Zr fraction for 0.4%, 0.8% and 1.2% excitation concentration. The response of the monoclinic layered perovskite compounds to electronic excitation and the variations of the amorphization temperature with respect to B -site composition are very different than those of cubic pyrochlore compounds. The evolution of the RDF features for monoclinic LTO at 300 K for three excitation concentrations is shown in Fig. 3(a). As for the cubic compounds, the formation of O—O bonds is maintained, but at an excitation concentration of 0.8%, which is lower than the 1.2% determined for the cubic lattice. Although 0.4% excitation concentration is not enough to have created O—O bonds in the crystal, the RDF features have been slightly affected compared to monoclinic LTO without electronic excitation (0%). At 0.4% lectronic excitation, no Ti—O or La—O bonds have been broken, but their distances have varied by 0.1-0.2 Å on average, resulting in a slightly broadened and less intense RDF short-range peak for Ti—O. As shown in Fig. 3(b), the transition temperature tends to decrease in each monoclinic Z-LTO compound as the excitation concentration increases. Interestingly, the monoclinic layered perovskite lattice is most resistant to electronic excitation when the chemical mixing in the B -site is the highest, which is obtained for 50% of Zr fraction. This is in contrast to the cubic pyrochlore lattice, for which the most vulnerable compound is obtained for a 50% [Zr, Ti] mixing. For 1.2% of excitation concentration and above, the transition temperature is composition independent and all the monoclinic compounds are readily amorphized below room temperature. Figure 3: (a) Effect of electronic excitation concentration on the radial distribution function (g(r)) of the monoclinic LTO crystal at 300 K. (b) Evolution of transition temperature as function of excitation concentration and composition for monoclinic Z-LTO compounds. A direct comparison of the overall interplay between structure, composition and electronic excitation concentration on transition temperature is shown in Fig. 4. For both 0.4% and 0.8% excitation concentration, the monoclinic compound made of 50% [Zr, Ti] mixing is found to be more resistant than cubic compound of the same composition. Fig. 4 also shows that the variations of the transition temperature with Zr fraction are larger for the cubic phase than for he monoclinic phase. This behavior suggests that the introduction of Ti in the cubic LZO phase has a greater impact on its resistance than the introduction of Zr in the monoclinic LTO phase. Breaking the high symmetry of the cubic phase by mixing the B -site composition has a more destabilizing effect than having a mixed B -site composition in the already low symmetry monoclinic phase. Interestingly, since the transition from monoclinic to cubic phase occurs at 53% Zr fraction, Fig. 4 shows that the 50% [Zr, Ti] composition of the cubic lattice, which is the most vulnerable to electronic excitation, should not be involved. This has been symbolized by the black dashed lines in Fig. 4. Therefore, starting from a monoclinic layered perovskite LTO structure, the gradual substitution of Ti by Zr would yield structures that are increasingly more resistant to electronic excitations, with the cubic pyrochlore LZO being the most resistant to electronic excitation. This overall trend is in agreement with several experimental observations , which found that the susceptibility to amorphization by electronic excitation increases with cation radius ratio r A /r B and Zr content. In the case of cubic LZO and monoclinic LTO, the r A /r B ratio is 1.61 and 1.92 respectively, making the LZO pyrochlore phase more resistant to electronic excitation. It is worth noting that the increase in electronic excitation resistance shown in Fig. 4 is not monotonic; a monoclinic compound with 50% Zr is predicted to have a similar transition temperature to a cubic compound with 75% Zr. While only the monoclinic and cubic phases were considered to model the phase transition between LTO and LZO, the potential for other phases to be involved in the transition region could lead to some uncertainty in the amorphization behavior in this regime. We also note that the trend in the resistance to electronic excitation follows the variations of the band gap ( 𝐸 * ) of the materials, which increases with Zr fraction. In particular, the band gap of the monoclinic phase with 50% Zr is slightly larger than the band gap of the cubic phase with 50% Zr, such that in order, 𝐸 *MonoclinicLTO < 𝐸 *CubicT-LZO <𝐸 *MonoclinicZ-LTO < 𝐸 *CubicLZO . This correlation supports the idea that the radiation response cannot be solely described in terms of cation radius ratio and that other factors, such as defect formation energy, bond type (i.e. covalency and ionicity), electronic structure of the compounds, and energy of order-disorder transition, also play an important role in the amorphization resistance of ceramics to electronic excitation . Figure 4: Effect of structure, composition and excitation concentration on the crystalline-to-amorphous transition temperature from the monoclinic perovskite to cubic pyrochlore structure. A comparison of our results with similar theoretical calculations performed for titanate pyrochlores by Xiao et al. complements and deepens our understanding of the effect of structure and composition on electronic excitation induced amorphization in ceramics. Especially, Fig. 5(a) shows that the evolution of transition temperature with excitation concentration for the monoclinic layered perovskite LTO has a similar slope compared to the other titanate pyrochlores, while cubic LZO pyrochlore exhibits a steeper slope. This suggests that a chemical change of the B -site cation more strongly affects the slope of the transition temperature as a function of excitation concentration, compared to a chemical change of the A -site cation. Fig. 5(a) also shows that cubic pyrochlores are more resistant to electronic excitation than the monoclinic perovskite LTO structure, as evidenced by their higher transition temperatures. In the specific case of a comparison between La Ti O and Gd Ti O , experiments have shown that the perovskite structure amorphizes more rapidly than the pyrochlore structure . This is in agreement with LTO having a lower crystalline-to-amorphous transition temperature compared to Gd Ti O . An analysis of the structural modifications obtained at the end of the simulations for cubic LZO and monoclinic LTO allows us to identify the main structural features conferring resistance to electronic excitation. Fig. 5(b) shows an example of structural modifications obtained at 600 K and 0.8% excitation concentration for both cubic LZO and monoclinic LTO crystals. For the cubic ZO crystal, the zigzag corner sharing ZrO octahedral arrangement is mostly preserved, even if some O -like molecules have been formed in the structure. However, in the case of the monoclinic LTO crystal, the formation of O -like molecules more strongly disturbs the octahedral connectivity. Due to the breaking of Ti—O bonds, new La-La layers can form in addition to those initially present in the structure, and the formation of TiO square based pyramids is possible. The structural origin of these modifications is the possibility for TiO octahedra to rotate along the axis by which they are linearly inter-connected. Octahedral rotations allow the monoclinic LTO structure to accommodate electronic excitation with more flexibility, but they also facilitate the formation of O -like molecules, making the amorphization easier; ultimately, this results in the structure being less resistant to excitation. In the cubic pyrochlore structures, such octahedral rotation is prevented by the zigzag connectivity network between ZrO octahedra. Figure 5: (a) Crystalline-to-amorphous transition temperature as function of electronic excitation for the lanthanum-based pyrochlores and comparison with the results from Ref[31]. (b) Example of structural deformations obtained for cubic LZO and monoclinic LTO crystals at a temperature of 600 K and an excitation concentration of 0.8%. onclusions: Ab initio molecular dynamics simulations investigating the interplay between electronic excitation, structure and composition on the amorphization temperature of lanthanum-based ceramics have been performed. The effect of B -site composition is that the monoclinic layered perovskite phase is more resistant to excitation for 50% [Zr, Ti] mixing, while the cubic pyrochlore phase is more vulnerable for a high [Zr, Ti] mixing. A direct structural comparison between the monoclinic and cubic phases indicates that at 50% [Zr, Ti] mixing, the monoclinic perovskite structure is most resistant; however, the undoped cubic LZO crystal is found to be more resistant to excitation than the undoped monoclinic LTO crystal. A structural analysis suggests that the structural feature conferring excitation resistance in cubic pyrochlore crystals is the octahedral connectivity network with a zigzag pattern. In contrast, the rotation of linearly connected octahedra in the monoclinic LTO phase facilitates the formation of O -like molecules in the structure. Our results illustrate how crystal structure and octahedral distortions mediate the amorphization process, suggesting ways to improve the radiation response of important functional ceramics and waste forms. eferences: Subramanian, M. A., Aravamudan, G., and Rao, G. V. S. Oxide pyrochlores – A review.
Prog. Solid. St. Chem. , 55-143 (1983). 2 Chakoumakos, B. C. Systematics of the pyrochlores structure type, Ideal A B X Y. J. Solid State Chem. , 120-129 (1984). 3 Lian, J., Helean, K. B., Kennedy, B. J., Wang, L. M., Navrotsky, A., and Ewing, R. C. Effect of structure and thermodynamic stability on the response of lanthanide stannate pyrochlores to ion beam irradiation.
J. Phys. Chem. B , 2343-2350 (2006). 4
Ewing, R. C., Weber, W. J., and Lian, J. Nuclear waste disposal – pyrochlores (A B O ). Nuclear waste form for the immobilization of plutonium and “minor” actinides. J. Appl. Phys. , 5949 (2004). 5 Weber, W. J., Wald, J. W., and Matzke, H. Self-radiation damage in gadolinium titanate (Gd Ti O ). Mater. Lett. , 173 (1985). 6 Weber, W. J., Wald, J. W., and Matzke, H. Effects of self-radiation damage in Cm-doped Gd Ti O and CaZrTi O . J. Nucl. Mater. , 196-209 (1986). 7
Weber, W. J. et al.
Radiation effects in crystalline ceramics for the immobilization of high-level nuclear waste and plutonium.
J. Mater. Res. , 1434-1484 (1998). 8 Weber, W. J., and Roberts, F. P. A review of radiation effects in solid nuclear waste forms.
Nucl. Technol. , 178 (1983). 9 Weber, W. J. et al.
Radiation effects in glasses used for immobilization of high-level waste and plutonium disposition.
J. Mater. Res. , 1946 (1997). 10 Sickafus, K. E. et al.
Radiation tolerance of complex oxides.
Science , 748 (2000). 11
Sickafus, K. E. et al.
Radiation-induced amorphization resistance and radiation tolerance in structurally related oxides.
Nat. Mater. , 217 (2007). 12 Cheng, J. G., Wang, J., Dechakupt, T., and Trolier-McKinstry, S. Low-temperature crystallized pyrochlore bismuth zinc niobate thin films by excimer laser annealing.
Appl. Phys. Lett. , 232905 (2005). 13 Shamoto, S., , H., Ono, Y., Nakano, T., Nozue, Y., and Kajitani, T. Light-induced metal-insulator transition in Lu V O . J. Phys. Chem. Solids , 325-329 (2001). 4 Wang, S. X., Begg, B. D., Wang, L. M., Ewing, R. C., Weber, W. J., and Govidan Kutty, K. V. Radiation stability of gadolinium zirconate: A waste form for plutonium disposition.
J. Mater. Res. , 4470-4473 (1999). 15 Sattonnay, G. et al.
Effect of composition on the behavior of pyrochlores irradiated with swift heavy ions.
Nucl. Instrum. Methods B , 261-265 (2012). 16
Sattonnay, G. et al.
Phase transformation induced by high electronic excitations in ion-irradiated Gd (Zr x Ti ) O pyrochlores. J. Appl. Phys.
Wang, S. X., Wang, L. M., Ewing, R. C., Was, G. S., Lumpkin, G. R. Ion irradiation-induced phase transformation of pyrochlore and zirconolite.
Nucl. Instrum. Methods Phys. Res. B , 704-709 (1999). 18
Lian, J., Ewing, R. C., Wang, L. M. Ion-beam irradiation of Gd Sn O and Gd Hf O pyrochlore: bond-type effect. J. Mater. Res. , 1575-1580 (2004). 19 Ewing, R. C., Wang, L. M. Amorphization of zirconolite: alpha-decay event damage versus krypton ion irradiation.
Nucl. Instrum. Methods B. , 319-323 (1992). 20 Begg, B. D. et al.
Heavy-ion irradiation effects on structures and acid dissolution of pyrochlores.
J. Nucl. Mater. , 208-216 (2001). 21
Lian, J., Zu, X. T., Kutty, K. V. G., Chen, J., Wang, L. M., and Ewing, R. C. Ion-irradiation-induced amorphization of La Zr O pyrochlore. Phys. Rev. B , 054108 (2002). 22 Lian, J. et al.
Radiation-induced amorphization of rare-earth titanate pyrochlores.
Phys. Rev. B , 134107 (2003). 23 Lang, M. et al.
Single-ion tracks in Gd Zr Ti x O pyrochlores irradiated with swift heavy ions. Phys. Rev. B. , 224105 (2009). 24 Park, S. et al.
Response of Gd Ti O and La Ti O to swift-heavy ion irradiation and annealing. Acta Mater. , 1-11 (2015). 25 Chartier, A., Meis, C., Crocombette, J. P., Weber, W. J., and Corrales, L. R. Molecular dynamic simulation of disorder induced amorphization in pyrochlores.
Phys. Rev. Lett. , 025505 (2005). 26 Williford, R. E., Weber, W. J., Devanathan, R., and Gale, J. D. Effects of cation disorder on oxygen vacancy migration in Gd Ti O . J. Electroceramics . , 409-424 (1999). 7 Jiang, C., Stanek, C. R., Sickafus, K. E., and Uberuaga B. P. First-principles prediction of disordering tendencies in pyrochlore oxides.
Phys. Rev. B , 104203 (2009). 28 Chartier, A., Meis, C., Weber, W. J., and René Corrales, L. Theoretical study of disorder in Ti-substituted La Zr O . Phys. Rev. B , 134116 (2002). 29 Chartier, A., Catillon, G., Crocombette, J. P. Key role of the cation interstitial structure in the radiation resistance of pyrochlores.
Phys. Rev. Lett. , 155503 (2009). 30
Li, X. et al.
Role of electronic excitation in the amorphization of Ge-Sb-Te alloys.
Phys. Rev. Lett.
Zhang, Y., Bae, I. T., and Weber, W. J. Atomic collision and ionization effects in oxides.
Nucl. Instrum. Methods B , 2828-2833 (2008). 32
Xiao, H.Y., Weber, W.J., Zhang, Y., Zu, X. T., and Li, S. Electronic excitation induced amorphization in titanate pyrochlores: an ab initio molecular dynamic study.
Sci. Rep. Zhang-JMR-2010: Zhang, J. et al.
Nanoscale phase transition under extreme conditions within an ion track.
J. Mater. Res. , 1344 (2010). 34 Zunger, A., Wei, S.-H., Ferreira, L.G. and Bernard, J.E. Special quasirandom structures.
Phys. Rev. Lett.
353 (1990). 35 van de Walle, A. et al.
Efficient stochastic generation of special quasi random structures.
Calphad
Kresse, G. and Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.
Phys. Rev. B
Wang, Y. and Perdew, J.P. Correlation hole of the spin-polarized electron gas, with exact small-wave-vector and high-density scaling.
Phys. Rev. B
Vosko, S. H., Wilk, L. and Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis.
Can. J. Phys.
Blöchl, P.E. Projector augmented-wave method.
Phys. Rev. B
Monkhorst, H.J. and Pack, J.D. Special points for Brillouin-zone integrations.
Phys. Rev. B
Sattonnay, G. et al. Phase Transformations in Pyrochlores Irradiated with Swift Heavy Ions: Influence of Composition and Chemical Bonding.
Acta Phys. Pol. A , 862-866 (2013).
Acknowledgement:
This work was supported by the Nuclear Process Science Initiative (NPSI) at the Pacific Northwest National Laboratory (PNNL). Computational resources were provided by PNNL Institutional Computing (PIC). This work was performed at Pacific Northwest National Laboratory which is operated by the Battelle Memorial Institute for the U. S. Department of Energy under Contract No. DEAC05-76RL0-1830.
Author contributions:
M.S. designed and conduced the calculations. All the authors discussed the results, wrote and reviewed the manuscript.