Advances in Condensed Matter Physics

Volume 2013 (2013), Article ID 675410, 8 pages

http://dx.doi.org/10.1155/2013/675410

Research Article

Electronic Structure Calculations of A₂Ti₂O₇ (A = Dy, Ho, and Y)

H. Y. Xiao

School of Physical Electronics, University of Electronic Science and Technology of China, Chengdu 610054, China

Received 8 July 2013; Accepted 30 July 2013

Academic Editor: Liang Qiao

Copyright © 2013 H. Y. Xiao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Ab initio calculations have been performed on titanate pyrochlores A₂Ti₂O₇ (A = Dy, Ho, and Y) to investigate their electronic structures. The generalized gradient approximation (GGA) + formalism has been used to correct the strong onsite Coulomb repulsion between the localized 4f electrons. The effects of effective values on the structural and electronic properties of A₂Ti₂O₇ (A = Dy, Ho, and Y) have been discussed. It is shown that Dy₂Ti₂O₇ and Ho₂Ti₂O₇ exhibit different electronic structures from Y₂Ti₂O₇. The strong interaction between Dy and Ho 4f electrons and O 2p orbitals may increase the covalency of and bonds and decrease their irradiation resistance.

1. Introduction

Materials with the A₂B₂O₇ pyrochlore structure have wide ranges of composition that lead to remarkable properties and wide variations in ionic and electronic conductivity, catalytic activity, and electrooptic and piezoelectric behavior [1]. Because they can be used to immobilize actinides [2–5], the pyrochlores have attracted significant attention both theoretically and experimentally [6–21]. In A₂B₂O₇ pyrochlore structure, the A and B cations occupy the 16 (0.5, 0.5, 0.5) and 16 (0, 0, 0) sites, respectively, and the oxygens are in the 48 (, 0125, 0.125) and 8 (0.375, 0.375, 0.375) positions (using the Wyckoff notation) [22]. Single crystals of the A₂Ti₂O₇ pyrochlores (A = Sm to Lu and Y) have been irradiated by 1 MeV Kr⁺ ions, and their microstructural evolutions as a function of increasing radiation dose have been characterized [7]. A slight deviation from the monotonic trend of critical amorphization temperature versus the ionic radius was observed for Y₂Ti₂O₇. is frequently used as a measurement of the resistance of a material to amorphization, and lower values of are the result of substantial dynamic annealing occurring during irradiation, which allows the materials to remain in the crystalline state. Generally, pyrochlores that are closer to the ideal fluorite structures are more susceptible to the radiation-induced pyrochlore-to-defect fluorite structural transition. The defect fluorite structure results from disordering of the A- and B-site cations, as well as the anion vacancies. Thus, pyrochlore compositions that are more easily disordered to the defect fluorite structure are more “resistant” to ion-beam-induced amorphization [7]. On the other hand, theoretical investigations reported by Sickafus et al. [18] have demonstrated that compounds with very dissimilar cationic radii (e.g., closer to the ideal pyrochlore structure, ) should exhibit the greatest susceptibility to structural destabilization (e.g., amorphization), whereas compounds with more similar radii (e.g., closer to the ideal fluorite structure) should behave more robustly in a radiation environment. For A₂Ti₂O₇ (A = Dy, Ho, and Y), the ionic radius ratio () decreases monotonically. According to Sickafus’ point of view, the compounds should become more radiation resistant as the cation A varies from Dy to Ho. However, the critical amorphization temperatures were measured to be 910, 850, and 780 K for Dy₂Ti₂O₇, Ho₂Ti₂O₇, and Y₂Ti₂O₇ compounds, respectively [7], meaning that of these three compounds, Y₂Ti₂O₇ is the most radiation resistant. It is clear that cation radius ratio alone cannot be used to explain the different responses of these compositions to radiation.

To resolve the discrepancy between experimental observations and theoretical predictions, it is important to fundamentally understand the electronic structures of A-site elements and their effects on the stability of the pyrochlore structure. In the present study, ab initio total energy calculations [17, 19, 23] based on density functional theory have been performed on A₂Ti₂O₇ (A = Dy, Ho, and Y) pyrochlores. GGA + formalism has been used to account for the strong on-site Coulomb repulsion among the 4f electrons in Dy₂Ti₂O₇ and Ho₂Ti₂O₇. How the structural and electronic properties of A₂Ti₂O₇ (A = Dy, Ho, and Y) pyrochlores are affected by the choice of is discussed. These calculations will provide significant insight into the effects of electronic configuration on the thermochemical stability of pyrochlore of different compositions.

2. Calculational Method

All the calculations have been completed using the VASP code [24] with spin-polarized effects taken into account. A primitive unit cell containing 22 atoms was used for the present investigations, with a Monkhorst-Pack -point mesh. The ion-electron interaction was described by PAW pseudopotentials with the following atomic valence configurations: Ti (3s², 3p⁶, 3d², and 4s²), Y (4s², 4p⁶, 4d¹, and 5s²), Dy (5s², 5p⁶, 4f¹⁰, and 6s²), and Ho (5s², 5p⁶, 4f¹¹, and 6s²). The PBE functional within the generalized gradient approximation was used to describe the exchange-correlation potential energy [25, 26], with the basis set for valence electrons consisting of plane waves with a cut-off energy of 400 eV. The calculations were performed based on ferromagnetic ordering of the magnetic moments and antiferromagnetic ordering, and spin ice model [27, 28] is not considered in the present work. The Hubbard correction was introduced using the method proposed by Dudarev et al. [29], in which the Hubbard parameter reflecting the strength of onsite Coulomb interaction and parameter adjusting the strength of exchange interaction are combined into a single parameter .

3. Results and Discussion

3.1. Atomic and Electronic Structure of Dy₂Ti₂O₇

The pyrochlore structure can be completely described by the a-cell edge, , and the 48 oxygen positional parameter, . The dependence of the lattice parameter , , and band gap of the Dy₂Ti₂O₇ on is shown in Figures 1(a), 1(b), and 1(c), respectively. As shown in Figure 1, the structural parameter decreases slightly for  eV. Above this value, the parameter changes more significantly. A of 2.6 eV yields a lattice constant of 10.12 Å, in excellent agreement with experiments [7]. For O48f positional parameter , it deviates from experimental value of 0.3275 with increasing values; that is, introducing makes the deviation from the experiment larger. Concerning the band structure of Dy₂Ti₂O₇, a sharp increase of band gap value with growing is observed. For a certain value of around 3.5 eV, the band gap value of 2.85 eV matches the calculated optical band gaps performed by Nemoshkalenko et al. [30]. It is shown that the band gap of Dy₂Ti₂O₇ presents stronger dependence on the value than lattice parameters.

(a)

(b)

(c)

Figure 1: Dependence of the lattice parameter (a), (b), and band gap (c) of the Dy₂Ti₂O₇ on .

Figure 2 displays the total density of state (DOS) distribution of Dy₂Ti₂O₇ at , 2.0, and 3.0 eV. The Fermi level is set to 0 eV. At the pure GGA level, the f band does not split but shows a large peak around the Fermi level, leading to a metallic ground state, which disagrees with experiments [31]. Obviously, the pure GGA calculation without modifying the intra-atomic Coulomb interaction gives wrong results. It is suggested that introduction of a penalty function which corrects the intraband Coulomb interaction by the Hubbard parameter is necessary for strongly correlated systems [32] such as Dy₂Ti₂O₇. If the is increased to 2 eV, the f bands undergo splitting, and a semiconducting solution with a finite separation of the occupied and unoccupied f band is found. The obtained band gap between the valence band edge (contributed by Dy 4f) and the bottom of the conduction band (contributed by Dy 4f) is 1.53 eV. In the case of  eV, the unoccupied f bands shift toward higher energy level, resulting in a larger band gap value of 2.36 eV. A notable difference between the cases of  eV and  eV is that there is no mixture between the occupied f bands and the O 2p orbitals in the valence region at  eV, whereas the occupied f bands become hybridized with the O 2p orbitals in the valence region at  eV.

(a)  eV

(b)  eV

(c)  eV

Figure 2: Total DOS distribution for Dy₂Ti₂O₇. (a)  eV; (b)  eV; and (c)  eV. The upper and lower panels show the spin-up and spin-down channels, respectively.

3.2. Atomic and Electronic Structure of Ho₂Ti₂O₇

Figure 3 shows the lattice parameters and band gaps of the Ho₂Ti₂O₇ as a function of . Generally, Ho₂Ti₂O₇ shows similar dependence on the effective values to the case of Dy₂Ti₂O₇. The lattice constants decrease with increasing , and the reason for this behavior is a slight hybridization of Ho 4f and O 2p orbitals. The experimental lattice constant [7] of 10.104 Å is obtained for  eV. For the O48f positional parameter, it deviates from experimental value of 0.3285 with increasing values. As compared with Dy₂Ti₂O₇, the band gap value increases less significantly with . At  eV, the band gap of Ho₂Ti₂O₇ is still 0.82 eV smaller than experimental measurement of 3.2 eV [33].

(a)

(b)

(c)

Figure 3: Dependence of the lattice parameter (a), (b), and band gap (c) of the Ho₂Ti₂O₇ on .

Figure 4 presents the total DOS of Ho₂Ti₂O₇ at , 2.0, and 3.0 eV. At  eV, the pure GGA calculation without modifying the intra-atomic Coulomb interaction yields a metallic ground state, in contrast to experiments [33]. As the is increased to 2 eV, the f bands splits, and the occupied and unoccupied f band contributes significantly to the valence bands and conduction bands, respectively. The corresponding band gap between the valence band edge and the bottom of the conduction band is 1.65 eV. Different from the total DOS distribution of Dy₂Ti₂O₇ at  eV, the occupied f bands hybridize with O 2p orbitals in the valence region. At  eV, it is noted that the main effect of increasing value is to push the unoccupied f bands toward higher energy level, resulting in a larger band gap value of 2.38 eV.

(a)  eV

(b)  eV

(c)  eV

Figure 4: Total DOS distribution for Ho₂Ti₂O₇. (a)  eV; (b)  eV; and (c)  eV. The upper and lower panels show the spin-up and spin-down channels, respectively.

3.3. Atomic and Electronic Structure of Y₂Ti₂O₇

Yttrium titanate pyrochlore is an important member of pyrochlore family, and it is often served as a model system for pyrochlores because of its simple electronic structure [12]; that is, no f electrons exist in Y₂Ti₂O₇. For this composition, GGA + method has also been employed to study if the intra-atomic electron correlations are important for Y-4d states. In Figure 5, we plot the equilibrium lattice parameters and band gap value of Y₂Ti₂O₇ as functions of . As shown in the figure, the lattice constant and O48f positional parameter deviate from experimental values of 10.1 Å and 0.33 [7] with increasing values. Specially, the choice of values has almost no effects on the band gap values. It is indicated that the intra-atomic electron correlations are negligible for Y-4d states in Y₂Ti₂O₇.

(a)

(b)

(c)

Figure 5: Dependence of the lattice parameter (a), (b), and band gap (c) of the Y₂Ti₂O₇ on .

The total DOS distribution of Y₂Ti₂O₇ at  eV is presented in Figure 6. The valence bands mainly consist of O 2p orbitals with small contribution from Ti 3d states, and the conduction bands are mainly composed of Ti 3d states hybridized with O 2p orbitals. The corresponding band gap is 2.82 eV, in excellent agreement with the calculated value of 2.84 eV reported by Jiang et al. [34]. This value is 0.32 eV larger than our previous work on Y₂Ti₂O₇ [10], as a result of different pseudopotential of Y employed.

Figure 6: Total DOS distribution for Y₂Ti₂O₇ at  eV.

3.4. Comparison of the Electronic Properties of Dy₂Ti₂O₇, Ho₂Ti₂O₇, and Y₂Ti₂O₇

The following discussions are based on the results obtained by  eV for Dy,  eV for Ho, and  eV for Y. Comparison of the partial DOS distribution of these three compositions is presented in Figure 7. It is noted that the Dy₂Ti₂O₇ and Ho₂Ti₂O₇ have similar DOS results; in contrast, Y₂Ti₂O₇ shows a very different character. This is probably due to the fact that 4f electrons play an important role in Dy₂Ti₂O₇ and Ho₂Ti₂O₇. The electronic structures of a series of titanate oxides A₂Ti₂O₇ (A = Sm-Er, Yb, and Lu) have been studied by Nemoshkalenko et al. [30] using X-ray photoelectron, emission spectroscopy, as well as the first-principles band structure calculations, where the lanthanide 4f states are assumed not to be hybridized with the other states. The calculations in the present work show that A-site 4f electrons do take part in the chemical bonding. For Dy₂Ti₂O₇ and Ho₂Ti₂O₇, the most striking features are the hybridization of Dy 4f and Ho 4f orbitals with O 2p orbitals, as shown in Figure 7. Especially, a strong hybridization occurs for Dy₂Ti₂O₇ in its upper valence region. The 4f electrons also contribute greatly to the lower conduction bands of Dy₂Ti₂O₇ and Ho₂Ti₂O₇. Unlike Dy₂Ti₂O₇ and Ho₂Ti₂O₇, interaction is more significant in Y₂Ti₂O₇, since its valence and conduction bands are mainly contributed by O 2p states hybridized with Ti 3d states and Ti 3d orbitals hybridized with O 2p orbitals, respectively.

(a)

(b)

(c)

Figure 7: Partial DOS distributions of (a) Dy₂Ti₂O₇; (b) Ho₂Ti₂O₇; and (c) Y₂Ti₂O₇.

Under irradiation, Y₂Ti₂O₇ is the most radiation resistant and Dy₂Ti₂O₇ is the least [7]. Since the pyrochlore compositions that are more easily disordered to the defect fluorite structure are more “resistant” to ion-beam-induced amorphization [18], Dy₂Ti₂O₇ is the least probable to transform into defect-fluorite structure and Y₂Ti₂O₇ is the most. This means that of the three compositions investigated, Dy₂Ti₂O₇ is the most stable thermodynamically and Y₂Ti₂O₇ is the least. The cation radius ratio criteria, as proposed by Sickafus et al. [18], clearly cannot be used to explain the different responses of these compositions to radiation. The radiation tolerance of nonmetallic solids has been correlated with the nature of the chemical bond in earlier work [35–38]. They demonstrated that the more covalently bonded materials are more readily amorphized at lower temperatures under heavy ion irradiation. For pyrochlores, the less covalently bonded compositions are more easily disordered to defect-fluorite structures [7], which are highly radiation resistant and remain crystalline at extreme radiation dose. In the present work, the strong interaction between Dy and Ho 4f electrons and O 2p orbitals may increase the covalency of and bonds and decrease the irradiation resistance of Dy₂Ti₂O₇ and Ho₂Ti₂O₇.

4. Conclusions

The electronic structures of A₂Ti₂O₇ (A = Dy, Ho, and Y) have been investigated using GGA + method. The effects of effective values on the structural and electronic properties of pyrochlores have been studied. It is shown that for strongly correlated systems such as Dy₂Ti₂O₇ and Ho₂Ti₂O₇, it is necessary to correct the intraband Coulomb interaction by the Hubbard parameter. We suggest that the electronic structure can be reasonably described with of  eV for Dy₂Ti₂O₇ and Ho₂Ti₂O₇.

Dy₂Ti₂O₇ and Ho₂Ti₂O₇ have similar DOS distribution; in contrast, Y₂Ti₂O₇ shows a very different character. The DOS distributions of Dy₂Ti₂O₇ and Ho₂Ti₂O₇ show that A-site 4f electrons hybridize significantly with O 2p orbitals in the valence region. Since Dy₂Ti₂O₇ and Ho₂Ti₂O₇ are less radiation resistant than Y₂Ti₂O₇, it is suggested that the strong interaction between Dy and Ho 4f electrons and O 2p electrons may increase the covalency of and bonds and decrease the irradiation resistance of Dy₂Ti₂O₇ and Ho₂Ti₂O₇.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 11004023) and by the Project Sponsored by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.

References

B. J. Wuensch, K. W. Eberman, C. Heremans et al., “Connection between oxygen-ion conductivity of pyrochlore fuel-cell materials and structural change with composition and temperature,” Solid State Ionics, vol. 129, no. 1, pp. 111–133, 2000.

R. C. Ewing, “Materials science: displaced by radiation,” Nature, vol. 445, no. 7124, pp. 161–162, 2007.

W. J. Weber and R. C. Ewing, “Plutonium immobilization and radiation effects,” Science, vol. 289, no. 5487, pp. 2051–2052, 2000.

R. C. Ewing, W. J. Weber, and J. Lian, “Nuclear waste disposal—pyrochlore (A₂B₂O₇): nuclear waste form for the immobilization of plutonium and “minor” actinides,” Journal of Applied Physics, vol. 95, no. 11 I, pp. 5949–5971, 2004.

G. R. Lumpkin, K. P. Hart, P. J. McGlinn, T. E. Payne, R. Giere, and C. T. Williams, “Retention of actinides in natural pyrochlores and zirconolites,” Radiochimica Acta, vol. 66-67, pp. 469–474, 1994.

A. Meldrum, C. W. White, V. Keppens, and L. A. Boatner, “Irradiation-induced amorphization of Cd₂Nb₂O₇ pyrochlore,” Physical Review B, vol. 63, no. 10, Article ID 104109, pp. 1041091–10410911, 2001.

J. Lian, J. Chen, L. M. Wang et al., “Radiation-induced amorphization of rare-earth titanate pyrochlores,” Physical Review B, vol. 68, no. 13, Article ID 134107, pp. 1341071–1341079, 2003.

J. Chen, J. Lian, L. M. Wang, R. C. Ewing, R. G. Wang, and W. Pan, “X-ray photoelectron spectroscopy study of disordering in Gd₂(Ti1-xZrx)₂O₇ pyrochlores,” Physical Review Letters, vol. 88, Article ID 105901, 4 pages, 2002.

H. Y. Xiao, X. T. Zu, F. Gao, and W. J. Weber, “First-principles study of energetic and electronic properties of A₂Ti₂O₇ (A=Sm, Gd, Er) pyrochlore,” Journal of Applied Physics, vol. 104, Article ID 073503, 6 pages, 2008.

Z. L. Zhang, H. Y. Xiao, X. T. Zu, F. Gao, and W. J. Weber, “First-principles calculation of structural and energetic properties for A₂Ti₂O₇ (A=Lu, Er, Y, Gd, Sm, Nd, La),” Journal of Materials Research, vol. 24, no. 4, pp. 1335–1341, 2009.

H. Y. Xiao, L. M. Wang, X. T. Zu, J. Lian, and R. C. Ewing, “Theoretical investigation of structural, energetic and electronic properties of titanate pyrochlores,” Journal of Physics, vol. 19, Article ID 346203, 2007.

H. Y. Xiao, F. Gao, and W. J. Weber, “Ab initio investigation of phase stability of Y₂Ti₂O₇ and Y₂Zr₂O₇ under high pressurePhysical Review B,” vol. 80, Article ID 212102, 4 pages, 2009.

Z. J. Chen, H. Y. Xiao, X. T. Zu et al., “Structural and bonding properties of stannate pyrochlores: a density functional theory investigation,” Computational Materials Science, vol. 42, no. 4, pp. 653–658, 2008.

N. Li, H. Y. Xiao, X. T. Zu et al., “First-principles study of electronic properties of La₂Hf₂O₇ and Gd₂Hf₂O₇,” Journal of Applied Physics, vol. 102, Article ID 063704, 2007.

H. Y. Xiao, F. X. Zhang, F. Gao, M. Lang, R. C. Ewing, and W. J. Weber, “Zirconate pyrochlores under high pressure,” Physical Chemistry Chemical Physics, vol. 12, no. 39, pp. 12472–12477, 2010.

H. Y. Xiao and W. J. Weber, “Pressure induced structural transformation in Gd₂Ti₂O₇ and Gd₂Zr₂O₇,” Journal of Physics, vol. 23, Article ID 35501, 2011.

J. M. Pruneda and E. Artacho, “First-principles study of structural, elastic, and bonding properties of pyrochlores,” Physical Review B, vol. 72, Article ID 085107, 8 pages, 2005.

K. E. Sickafus, L. Minervini, R. W. Grimes et al., “Radiation tolerance of complex oxides,” Science, vol. 289, no. 5480, pp. 748–751, 2000.

R. Terki, H. Feraoun, G. Bertrand, and H. Aourag, “Full potential linearized augmented plane wave investigations of structural and electronic properties of pyrochlore systems,” Journal of Applied Physics, vol. 96, no. 11, pp. 6482–6487, 2004.

A. Chartier, C. Meis, W. J. Weber, and L. R. Corrales, “Theoretical study of disorder in Ti-substituted La₂Zr₂O₇,” Physical Review B, vol. 65, Article ID 134116, 11 pages, 2002.

R. E. Williford, W. J. Weber, R. Devanathan, and J. D. Gale, “Effects of cation disorder on oxygen vacancy migration in Gd₂Ti₂O₇,” Journal of Electroceramics, vol. 3, no. 4, pp. 409–424, 1999.

X. J. Wang, H. Y. Xiao, X. T. Zu, Y. Zhang, and W. J. Weber, “Ab initio molecular dynamics simulations of ion–solid interactions in Gd₂Zr₂O₇ and Gd₂Ti₂O₇,” Journal of Materials Chemistry C, vol. 1, p. 1665, 2013.

W. R. Panero, L. Stixrude, and R. C. Ewing, “First-principles calculation of defect-formation energies in the Y₂(Ti,Sn,Zr)₂O₇ pyrochlore,” Physical Review B, vol. 70, Article ID 054110, 11 pages, 2004.

G. Kresse and D. Joubert, “From ultrasoft pseudopotentials to the projector augmented-wave method,” Physical Review B, vol. 59, no. 3, pp. 1758–1775, 1999.

J. P. Perdew, J. A. Chevary, S. H. Vosko et al., “Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation,” Physical Review B, vol. 46, no. 11, pp. 6671–6687, 1992.

J. A. White and D. M. Bird, “Implementation of gradient-corrected exchange-correlation potentials in Car-Parrinello total-energy calculations,” Physical Review B, vol. 50, no. 7, pp. 4954–4957, 1994.

M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske, and K. W. Godfrey, “Geometrical frustration in the ferromagnetic pyrochlore Ho₂Ti₂O₇,” Physical Review Letters, vol. 79, no. 13, pp. 2554–2557, 1997.

A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, and B. S. Shastry, “Zero-point entropy in ‘spin ice’,” Nature, vol. 399, no. 6734, pp. 333–335, 1999.

S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, “Electron-energy-loss spectra and the structural stability of nickel oxide: an LSDA+U study,” Physical Review B, vol. 57, no. 3, pp. 1505–1509, 1998.

V. V. Nemoshkalenko, S. V. Borisenko, V. N. Uvarov et al., “Electronic structure of the R₂Ti₂O₇ (R=Sm-Er, Yb, Lu) oxides,” Physical Review B, vol. 63, no. 7, Article ID 075106, pp. 0751061–0751068, 2001.

A. K. Pandit, T. H. Ansari, R. A. Singh, and B. M. Wanklyn, “Electrical conduction in Dy₂Ti₂O₇ single crystal,” Materials Letters, vol. 11, no. 1-2, pp. 52–58, 1991.

V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, “First-principles calculations of the electronic structure and spectra of strongly correlated systems: the LDA + U method,” Journal of Physics Condensed Matter, vol. 9, no. 4, pp. 767–808, 1997.

S. T. Bramwell and M. J. P. Gingras, “Spin ice state in frustrated magnetic pyrochlore materials,” Science, vol. 294, no. 5546, pp. 1495–1501, 2001.

Y. Jiang, J. R. Smith, and G. Robert Odette, “Prediction of structural, electronic and elastic properties of Y₂Ti₂O₇ and Y₂TiO₅,” Acta Materialia, vol. 58, no. 5, pp. 1536–1543, 2010.

H. M. Naguib and R. Kelly, “Criteria for bombardment-induced structural changes in non-metallic solids,” Radiation Effects, vol. 25, no. 1, pp. 1–12, 1975.

K. Trachenko, “Understanding resistance to amorphization by radiation damage,” Journal of Physics Condensed Matter, vol. 16, no. 49, pp. R1491–R1515, 2004.

K. Trachenko, J. M. Pruneda, E. Artacho, and M. T. Dove, “How the nature of the chemical bond governs resistance to amorphization by radiation damage,” Physical Review B, vol. 71, no. 18, Article ID 184104, 2005.

R. C. Ewing, W. J. Weber, and J. Lian, “Nuclear waste disposal-pyrochlore (A₂B₂O₇): nuclear waste form for the immobilization of plutonium and “minor” actinides,” Journal of Applied Physics, vol. 95, no. 11 I, pp. 5949–5971, 2004.