DOI:
10.1039/D2RA02419E
(Paper)
RSC Adv., 2022,
12, 16717-16722
Mechanisms of temperature-dependent oxygen absorption/release and appearance of intermediate phase in κ-Ce2Zr2O8: study based on oxygen vacancy formation energy computations†
Received
14th April 2022
, Accepted 30th May 2022
First published on 6th June 2022
Abstract
This study clarified the mechanisms of the temperature-dependent oxygen absorption/release properties and appearance of the intermediate phase for κ-Ce2Zr2O8, which is known to have a high oxygen storage/release capacity (OSC). First-principle computations revealed that the vacancy formation energies depend on the number of vacancies and can be categorized into two groups: low-energy and high-energy. The intermediate phase observed experimentally was assigned to the state after all the oxygen vacancies in the low-energy group were formed. We also found that the mechanism of the improved OSC performance by Ti substitution could be explained in terms of the vacancy formation energies.
1 Introduction
Oxygen storage/release materials have been actively investigated in a wide variety of fields, including environmental protection1 and energy conversion.2 In these fields, ceria (CeO2)-based materials are known to have high oxygen storage/release capacities (OSCs) and are used in industrial applications such as automobile three-way catalysts,3–5 water–gas shift reactions,6 and solid oxide fuel cells.7 In automobile three-way catalysts, ceria-based materials are used as oxygen storage/release materials to adequately control the oxygen concentration on the catalyst surfaces by storing oxygen in an oxidized atmosphere and releasing it under a reduced atmosphere.3,5,8 Among the ceria-based materials, cation-ordered Ce2Zr2O8 (κ-Ce2Zr2O8) has been reported to have an outstanding OSC, where almost all the Ce atoms contribute to the redox conversion.9–11 The mechanism of the high OSC performance of κ-Ce2Zr2O8 has been studied both experimentally12,13 and theoretically.14,15 These studies have clarified that oxygen vacancies are formed at oxygen sites coordinated only with Zr atoms, even though the binding strength of Zr4+–O2− is usually stronger than that of Ce4+–O2−. The oxygen vacancy formation energy is known to be dominated by the structural relaxations that occur after the formation of oxygen vacancies. In κ-Ce2Zr2O8, when the Zr4+–O2− bond is broken and an oxygen vacancy is formed, the remaining oxygen atoms move toward the vacancy sites as a result of the driving force caused by the differences in the electrostatic forces and ionic radii of Zr and Ce. The oxygen vacancies at the Zr sites are stabilized by this relaxation. Simultaneously, the valency of Zr is retained at 4+ by electron transfer via the redox reaction of Ce from 4+ to 3+. This electron transfer prevents the significant movement of Zr atoms, and therefore, the energetically unfavorable lattice distortion is minimized. This latter feature enables the topotactic transformation of κ-Ce2Zr2O8 into Ce2Zr2O7 pyrochlore by maximizing the locality of the structural relaxation.14 In theoretical studies, the oxygen vacancy formation energies were examined using density functional theory (DFT) calculations and used to understand the atomic-scale structural changes and oxygen vacancy formation properties.14,15 Experimental analyses also clarified the presence of the intermediate phase, β-Ce2Zr2O7.5, in the space group F
3m between the κ-Ce2Zr2O8 (Fm
m) and Ce2Zr2O7 pyrochlore (Fd
m) states.13,16 The existence of the intermediate state results in the two-step oxygen absorption/release behavior observed during heating under an oxidized/reduced atmosphere.16 Recently, Ti substitution in Ce0.5TixZr0.5−xO2 was found to greatly enhance the OSC performance.17 Ce0.5Ti0.1Zr0.4O2 was found to have approximately twice the OSC value of κ-Ce2Zr2O8 at a lower temperature (473 K). A bond valence sum (BVS) analysis18,19 showed that this OSC enhancement could be explained by the formation of weakly bound oxygen atoms.17 However, despite the efforts made in this field, the mechanisms of the appearance of the intermediate phase and the two-step oxygen absorption/release behavior have not yet been fully explained. The mechanism of the OSC enhancement by Ti substitution is also unclear.
In this study, we systematically investigated the formation of oxygen vacancies from a fully oxidized state (κ-Ce2Zr2O8) to a fully reduced state (Ce2Zr2O7 pyrochlore) using DFT, including the Hubbard-U term (DFT+U).20,21 The same method was applied to Ti-substituted systems.
The systematic computations of the oxygen vacancy formation energies revealed that they depended on the number of vacancies and could be categorized into two groups: low- and high-energy. The intermediate phase obtained experimentally was assigned to the state after all the oxygen vacancies in the low-energy group were formed. The theoretical conclusion was verified by comparing the simulated X-ray diffraction (XRD) spectra with the experimental results. This new finding reasonably explained the two-step oxygen absorption/release behavior observed in the experiments. The OSC enhancement at low temperatures by Ti substitution was also reasonably explained by the oxygen vacancy formation energy computations.
2 Method
2.1 First-principle computations
The periodic DFT method was used for the electronic structure computations. The computations were performed using the vienna ab initio simulation package (VASP).22,23 The generalized gradient approximation (GGA) of the Perdew–Burke–Ernzerhof (PBE) functional24 was used to describe the exchange–correlation interaction of electrons. The projector-augmented wave (PAW) method25 was used to describe the core–valence interaction. GGA cannot accurately describe the electronic structure of reduced ceria-based materials because of the presence of the self-interaction error.20,26–28 Therefore, an on-site Coulomb correction was included using the DFT+U method29 to properly incorporate the behavior of Ce 4f electrons. The previous study showed that significant delocalisation because of the self-interaction error that remained for U < 5 eV, while the results were converged for U ≥ 5 eV where the localisation was found to be independent of the choice of U.20,26 The similar oxygen vacancy formation energies of CeO2 have been reported for U greater than 5 eV (3.15, 3.19, 3.20 eV for U = 5, 6.1, 6.3 eV, respectively).15,20,26 Therefore, the value of U was set to 5 eV in this study. Spin-polarized calculations were performed to correctly treat the sub-stoichiometric oxygen. The cut-off energy of the plane wave was set to 400 eV, which has been validated for the CeZrO2 system.15 The k-point sampling was performed using a 2 × 2 × 2 Monkhorst-Pack mesh, except for the calculation of gas-phase oxygen molecules, which was performed using the Γ point and a cubic unit cell with 10 Å per side. Both the lattice and atomic positions were fully relaxed to obtain a stable structure. No symmetries were applied to the structural optimization. To investigate the OSC mechanism in detail, Bader charge analysis was performed using the code of Henkelman et al.30
2.2 Model systems
The 96 atom super cell (Ce16Zr16O64) shown in Fig. 1 was used for the κ-Ce2Zr2O8 computations. The 64 oxygen atoms in the 96 atom super cell can be classified into three types: 8 Oa (coordinated only with Ce), 8 Ob (coordinated only with Zr), and 48 Oc (coordinated with both Ce and Zr atoms),10–12 as shown in Fig. 1. Oxygen vacancies are known to preferentially form at Ob sites.12,14 When all 8 Ob atoms are removed (complete reduction), the remaining 88 atom unit cell represents Ce2Zr2O7 pyrochlore.
![image file: d2ra02419e-f1.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-f1.gif) |
| Fig. 1 Super cell of Ce16Zr16O64 (96 atoms) for κ-Ce2Zr2O8. The cations (Ce4+ and Zr4+) were ordered. Examples of three types of oxygen atoms, Oa, Ob, and Oc, are indicated. | |
To model the Ti-substituted system, we selected a Ce16Ti2Zr14O64 system. The Ti cation ratio of this system (6.25 mol%) was close to the composition that exhibited a significant OSC improvement in an experiment (9.7 mol%).17 A Ti-substituted system was experimentally confirmed to have a cation-ordered structure.17 The preferential vacancy formation sites were considered to be the same as those in the κ-Ce2Zr2O8 system. The formation of oxygen vacancies was expected to be activated by Ti substitution because the Ti–O bond is weaker than the Zr–O bond as a result of the greater electronegativity of Ti4+ (the electronegativities estimated by the empirical model31 were 1.73 and 1.61 for Ti4+ and Zr4+, respectively). Therefore, 8 Ob sites were chosen as vacancy formation sites. The remaining task was to determine the Ti substitution sites.
Several combinations of Zr sites can be replaced by Ti atoms. We examined the calculated total energies of all the symmetrically independent combinations of Ti substitution sites to determine the most stable combination (see ESI† for the calculation result). This examination indicated that in the most stable state, with two Ti substitutions, the two substituted Ti atoms were adjacent to each other (see Fig. 2). Because the ionic radius of Ti4+ is smaller than that of Zr4+ (Zr4+: 0.84 Å, Ti4+: 0.74 Å),32 the substitution caused a structure distortion around the substituted sites. This distortion could be localized using the configuration shown in Fig. 2 because the induced strain could be cancelled by adjacent substituted Ti atoms. This was considered to be the reason that two substituted Ti atoms energetically preferred the adjacent substitution sites. The oxygen vacancy formation energies of Ce16Ti2Zr14O64 were computed using this atomic configuration.
![image file: d2ra02419e-f2.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-f2.gif) |
| Fig. 2 Optimized Ti-substituted structure for Ce16Ti2Zr14O64 system. | |
2.3 Oxygen vacancy formation energy
The oxygen vacancy formation energy, Efv, relative to the gas-phase O2 molecule was computed using the following equation: |
![image file: d2ra02419e-t1.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-t1.gif) | (1) |
Here, N is the number of oxygen vacancies, E[N] is the energy of the bulk with N oxygen vacancies, and E[Ogas2] is the energy of the oxygen molecule in vacuum. As previously described, in the bulk energy calculations, the cell volume, shape, and atomic positions were fully relaxed. Unlike previous studies, where Efv(N) was computed only for a discrete number of vacancies (N = 1, 4, and 8 in ref. 14), we systematically computed Efv(N) for all cases (N = 1–8), from the fully oxidized state (κ-Ce16Zr16O64) to the fully reduced state (Ce16Zr16O56 pyrochlore).
The oxygen vacancy sites in Ce16Zr16O64 were assumed to be oxygen atoms coordinated only to Zr atoms (Ob, see Fig. 1). Although all the Ob atoms were symmetrically equivalent for N = 1, some of the remaining Ob atoms became inequivalent after the vacancy formation. In this study, the oxygen vacancy formation energies were examined for all symmetrically inequivalent combinations of Ob sites, and the most stable combinations were selected for Efv(N). In the case of Ce16Ti2Zr14O64, the Ob atoms, some of which were coordinated to both Zr and Ti atoms, were not symmetrically equivalent even for N = 1. The oxygen vacancy formation energies were computed for all the combination patterns of the Ob sites, and the most stable combinations were selected for Efv(N).
2.4 Symmetry analysis and XRD pattern simulation
The XRD patterns were simulated using the VESTA program33 to assign the theoretically obtained structure to the experimentally observed intermediate phase. The FINDSYM program34 was used to examine the symmetry of the theoretically obtained crystals and compare them with the space group of the intermediate phase. The tolerance parameters of the program were set as follows: lattice: 0.01, position: 0.1, and rotational moment: 0.1.
3 Results and discussion
Fig. 3 shows the computed oxygen vacancy formation energies, Efv(N), for κ-Ce16Zr16O64 and Ti-substituted Ce16Ti2Zr14O64. For both κ-Ce16Zr16O64 and Ti-substituted Ce16Ti2Zr14O64, the dependence of the formation energy on N appeared to be a sigmoid function, indicating that the Efv(N) values could be categorized into two groups: low-energy (first four) and high-energy (last three) groups. The vacancy formation energies of the high-energy group were less than 0.5 eV. These were significantly smaller than those of other materials (e.g., 2.09 eV for tetragonal Ce2Zr2O8 (ref. 14)). This indicated that κ-Ce2Zr2O8 can release oxygen at relatively low temperatures, making it an excellent OSC material. The oxygen vacancy formation energies at N = 5 were between those of the two groups. The locations of the low-energy and fifth vacancy sites are shown in Fig. 4. A close examination of the structure indicated that the four vacancy sites in the low-energy group were at locations that maximized their distances. As reported in the literature,14 the sacrificial reduction of the surrounding Ce atoms can produce oxygen vacancies with little charge fluctuation of the Zr atoms in the κ-Ce2Zr2O8 system. The maximized distances allow Ce atoms to efficiently influence the vacancy sites via the sacrificial reduction mechanism. It should also be noted that the lattice distortions caused by vacancy formation can be minimized by maximizing the distance. Fig. 4 shows that the fifth oxygen vacancy site was located at the center of the four low-energy group sites. This special location was considered to suppress the lattice distortion and make the vacancy formation energy relatively small. However, when additional oxygen vacancies were generated at the high-energy group sites, the symmetry was broken, and the additional vacancy sites had to be adjacent to the existing vacancy sites. In these cases, the distance between vacancy sites became significantly shorter. Therefore, the additional vacancy formations could not benefit from the sacrificial reduction of Ce atoms, and the lattice distortions became significantly large. Accordingly, the two distinct energy groups were attributed to the symmetry of the vacancy formation sites.
![image file: d2ra02419e-f3.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-f3.gif) |
| Fig. 3 Vacancy formation energies for κ-Ce16Zr16O64 and Ce16Ti2Zr14O64 systems. | |
![image file: d2ra02419e-f4.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-f4.gif) |
| Fig. 4 Locations of vacancy sites for low-energy group and 5th vacancy site. | |
In the experiment,13 a stable intermediate phase with space group F
3m at a composition of Ce2Zr2O7.5 was observed. In this study, the composition of the intermediate phase (Ce16Zr16O60) was assigned to N = 4. Therefore, the experimentally observed intermediate phase was considered to be the state after all the oxygen vacancies in the low-energy group were formed. This assignment was supported by the simulated XRD patterns of the structures of N = 4 (Ce16Zr16O60), N = 0 (κ-Ce16Zr16O64), N = 5 (Ce16Zr16O59), and N = 8 (Ce16Zr16O56 pyrochlore). The simulated XRD patterns are shown in Fig. 5. The peaks with a Miller index of (200), which are characteristic of the experimentally observed intermediate phase,13 are visible in the XRD patterns of Ce16Zr16O60 and Ce16Zr16O59, but do not appear in the XRD patterns of Ce16Zr16O56 and Ce16Zr16O64. For more precise assignments, the space groups were examined for the computed structures of Ce16Zr16O60 and Ce16Zr16O59. The space group of Ce16Zr16O60 was identified as F
3m, which corresponded to the experimentally detected space group of the intermediate phase, whereas that of Ce16Zr16O59 was identified as P
3m. Hence, the experimentally observed intermediate phase was assigned to the state after all the oxygen vacancies were formed at the low-energy group sites. The theoretically obtained formation energy gap between the low- and high-energy groups, as shown in Fig. 3, explains the stability of this intermediate phase.
![image file: d2ra02419e-f5.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-f5.gif) |
| Fig. 5 Simulated XRD patterns for (a) Ce16Zr16O64, (b) Ce16Zr16O60, (c) Ce16Zr16O59, and (d) Ce16Zr16O56. The peaks corresponding to a Miller index of (200), which is characteristic of the intermediate phase, are highlighted by arrows. | |
As shown in Fig. 3, the oxygen vacancy formation energies in the low-energy group were reduced by Ti substitution, whereas no significant change was observed in the high-energy group. This indicated that the release temperatures of the oxygen atoms at the low-energy group sites decreased. In other words, the amount of oxygen released at low temperatures increased. This feature was consistent with experimental results.17 Unlike κ-Ce16Zr16O64, the eight oxygen sites, Ob, in the Ce16Ti2Zr14O64 system are not equivalent because of the Ti substitution. Based on the electronegativities of Ti (1.54) and Zr (1.33), we expected that oxygen atoms directly coordinated to Ti would be more easily released than oxygen atoms coordinated only to Zr because oxygen cannot take as many electrons from Ti as Zr. This expectation could be true if vacancies were formed without any structural relaxation, as verified by our preliminary DFT+U calculations on frozen structures. However, when the lattice relaxation was switched on, the vacancy formation energy for the oxygen atoms coordinated only to Zr was smaller than that of Ti. This stemmed from the large strain caused by the removal of oxygen coordinated directly by Ti, which has a smaller ionic radius than Zr.
Although oxygen vacancies were formed at the Zr sites regardless of the Ti substitution, the oxygen vacancy formation energies were significantly reduced by the Ti substitution. To investigate the mechanism of this reduction, Bader charge analyses were conducted for Ce16Ti2Zr14O64 and κ-Ce16Zr16O64. The electron numbers obtained by these Bader charge analyses for all the atoms before and after the first oxygen vacancy formation are compared in Fig. 6. The sacrificial reduction of Ce atoms was observed in the κ-Ce16Zr16O64 system, where the electron numbers of three specific Ce atoms increased when an oxygen vacancy was formed. The reduced Ce sites formed a triangle around the vacancy sites (inset in Fig. 6). However, the electron numbers of Zr atoms, which involved Zr atoms directly coordinating the oxygen vacancy sites, did not change significantly. This result was consistent with that of a previous study,14 although two Ce atoms instead of three were reduced during the oxygen vacancy formation in past calculations, where the supercell was fixed during the calculations. This number difference was attributed to the relaxation of the supercell: two Ce atoms were reduced when the supercell was fixed in our preliminary calculations. In contrast, in Ce16Ti2Zr14O64 systems, the electron amounts for the Ce atoms were already increased with the Ti substitution before the vacancy formation, as shown in Fig. 6. These electrons were confirmed to originate from the oxygen atoms (the Bader charge for the oxygen atoms was decreased by the Ti substitution). However, the number of electrons in the Zr atoms was not affected. The decrease in the number of oxygen electrons resulted in a decrease in the strength of the ionic bonds. This was considered to be the cause of the decrease in the oxygen vacancy formation energies.
![image file: d2ra02419e-f6.tif](/image/article/2022/RA/d2ra02419e/d2ra02419e-f6.gif) |
| Fig. 6 Bader charge amounts for each cation before and after oxygen vacancy formation for κ-Ce16Zr16O64 and Ce16Ti2Zr14O64 systems. Three specific reduced Ce atoms are indicated by blue arrows, and their positions are shown in the inset. | |
4 Conclusion
In this study, we systematically investigated the formation of oxygen vacancies from the fully oxidized state (κ-Ce2Zr2O8) to the fully reduced state (Ce2Zr2O7 pyrochlore) using first-principles computations based on the DFT+U method. Systematic computations of the oxygen vacancy formation energies of cation-ordered Ce2Zr2O8 revealed that the arrangement of vacancy sites was determined so that the distance between vacancy sites was as large as possible, and the symmetry was as high as possible. This ordered arrangement minimized the lattice distortion and maximized the benefit of Ce reduction during the vacancy formation. This feature resulted in the appearance of two oxygen vacancy formation energy groups: low- and high-energy groups. The intermediate phase obtained experimentally was assigned to the state after all the oxygen vacancies in the low-energy group were formed on the basis of the simulated XRD spectra. The theoretically obtained oxygen vacancy energies explained the two-step oxygen absorption/release behavior observed experimentally during heating under oxidized and reduced atmospheres. In addition, this ordered vacancy formation was considered to be one of the reasons for the possibility of topotactic transformation between κ-Ce2Zr2O8 and Ce2Zr2O7 pyrochlore. The OSC enhancement by Ti substitution at lower temperatures was also reproduced by the theoretical oxygen vacancy formation energies. The charge transfer from oxygen atoms to Ce atoms caused by the Ti substitution was considered to be the main factor for the improved OSC.
Conflicts of interest
The authors declare that they have no conflicts of interest.
Acknowledgements
The authors wish to acknowledge their useful discussions with Dr Akira Morikawa (Toyota Central R&D Labs, Inc.).
References
- S. Carrettin, P. Concepción, A. Corma, J. M. Lopez Nieto and V. F. Puntes, Nanocrystalline CeO2 increases the activity of Au for CO oxidation by two orders of magnitude, Angew. Chem., Int. Ed., 2004, 43(19), 2538–2540 CrossRef CAS PubMed.
- G. A. Deluga, J. R. Salge, L. D. Schmidt and X. E. Verykios, Renewable hydrogen from ethanol by autothermal reforming, Science, 2004, 303(5660), 993–997 CrossRef CAS PubMed.
- A. F. Diwell, R. R. Rajaram, H. A. Shaw and T. J. Truex, The role of ceria in three-way catalysts, in Studies in Surface Science and Catalysis, Elsevier, 1991, vol. 71, pp. 139–152 Search PubMed.
- J. Kašpar, P. Fornasiero and M. Graziani, Use of CeO2-based oxides in the three-way catalysis, Catal. Today, 1999, 50(2), 285–298 CrossRef.
- J. Kašpar and P. Fornasiero, Nanostructured materials for advanced automotive de-pollution catalysts, J. Solid State Chem., 2003, 171(1–2), 19–29 CrossRef.
- Q. Fu, H. Saltsburg and M. Flytzani-Stephanopoulos, Active nonmetallic Au and Pt species on ceria-based water-gas shift catalysts, Science, 2003, 301(5635), 935–938 CrossRef CAS PubMed.
- Z. Xu, Z. Qi and A. Kaufman, Effect of oxygen storage materials on the performance of proton-exchange membrane fuel cells, J. Power Sources, 2003, 115(1), 40–43 CrossRef CAS.
- P. Li, X. Chen, Y. Li and J. W. Schwank, A review on oxygen storage capacity of CeO2-based materials: Influence factors, measurement techniques, and applications in reactions related to catalytic automotive emissions control, Catal. Today, 2019, 327, 90–115 CrossRef CAS.
- Y. Nagai, T. Yamamoto, T. Tanaka, S. Yoshida, T. Nonaka, T. Okamoto, A. Suda and M. Sugiura, X-ray absorption fine structure analysis of local structure of CeO2–ZrO2 mixed oxides with the same composition ratio (Ce/Zr=1), Catal. Today, 2002, 74(3–4), 225–234 CrossRef CAS.
- H. Kishimoto, T. Omata, S. Otsuka-Yao-Matsuo, K. Ueda, H. Hosono and H. Kawazoe, Crystal structure of metastable κ-CeZrO4 phase possessing an ordered arrangement of Ce and Zr ions, J. Alloys Compd., 2000, 312(1–2), 94–103 CrossRef CAS.
- A. Suda, Y. Ukyo, H. Sobukawa and M. Sugiura, Improvement of oxygen storage capacity of CeO2–ZrO2 solid solution by heat treatment in reducing atmosphere, J. Ceram. Soc. Jpn., 2002, 110(1278), 126–130 CrossRef CAS.
- T. Yamamoto, A. Suzuki, Y. Nagai, T. Tanabe, F. Dong, Y. Inada, M. Nomura, M. Tada and Y. Iwasawa, Origin and dynamics of oxygen storage/release in a Pt/ordered CeO2–ZrO2 catalyst studied by time-resolved XAFS analysis, Angew. Chem., Int. Ed., 2007, 46(48), 9253–9256 CrossRef CAS PubMed.
- T. Sasaki, Y. Ukyo, K. Kuroda, S. Arai, S. Muto and H. Saka, Crystal Structure of Ce2Zr2O7 and β-Ce2Zr2O7.5, J. Ceram. Soc. Jpn., 2004, 112(8), 440–444 CrossRef CAS.
- H.-F. Wang, Y.-L. Guo, G.-Z. Lu and P. Hu, Maximizing the Localized Relaxation: The Origin of the Outstanding Oxygen Storage Capacity of κ-Ce2Zr2O8, Angew. Chem., Int. Ed., 2009, 48(44), 8289–8292 CrossRef CAS PubMed.
- H.-T. Chen and J.-G. Chang, Oxygen vacancy formation and migration in Ce1−xZrx O2 catalyst: A DFT+U calculation, J. Chem. Phys., 2010, 132(21), 214702 CrossRef PubMed.
- T. Sasaki, Y. Ukyo, A. Suda, M. Sugiura, K. Kuroda, S. Arai and H. Saka, Oxygen Absorption Behavior of Ce2Zr2O7+x and Formation of Ce2Zr2O7.5, J. Ceram. Soc. Jpn., 2003, 111(1294), 382–385 CrossRef CAS.
- Y. Goto, A. Morikawa, M. Iwasaki, M. Miura and T. Tanabe, Enhanced oxygen storage capacity of cation-ordered cerium–zirconium oxide induced by titanium substitution, Chem. Commun., 2018, 54(28), 3528–3531 RSC.
- I. D. Brown and D. Altermatt, Bond-valence parameters obtained from a systematic analysis of the inorganic crystal structure database, Acta Crystallogr., Sect. B: Struct. Sci., 1985, 41(4), 244–247 CrossRef.
- N. E. Brese and M. O'keeffe, Bond-valence parameters for solids, Acta Crystallogr., Sect. B: Struct. Sci., 1991, 47(2), 192–197 CrossRef.
- M. Nolan, S. C. Parker and G. W. Watson, The electronic structure of oxygen vacancy defects at the low index surfaces of ceria, Surf. Sci., 2005, 595(1–3), 223–232 CrossRef CAS.
- S. Fabris, G. Vicario, G. Balducci, S. de Gironcoli and S. Baroni, Electronic and atomistic structures of clean and reduced ceria surfaces, J. Phys. Chem. B, 2005, 109(48), 22860–22867 CrossRef CAS PubMed.
- G. Kresse and J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54(16), 11169 CrossRef CAS.
- J. Hafner, Materials simulations using VASP—a quantum perspective to materials science, Comput. Phys. Commun., 2007, 177(1–2), 6–13 CrossRef CAS.
- J. P. Perdew, K. Burke and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett., 1996, 77(18), 3865 CrossRef CAS PubMed.
- P. E. Bliichl, Projector augmented-wave method, Phys. Rev. B: Condens. Matter Mater. Phys., 1994, 50(24), 17953–17979 CrossRef.
- M. Nolan, S. Grigoleit, D. C. Sayle, S. C. Parker and G. W. Watson, Density functional theory studies of the structure and electronic structure of pure and defective low index surfaces of ceria, Surf. Sci., 2005, 576(1–3), 217–229 CrossRef CAS.
- H.-T. Chen, J.-G. Chang, H.-L. Chen and S.-P. Ju, Identifying the O2 diffusion and reduction mechanisms on CeO2 electrolyte in solid oxide fuel cells: a DFT+ U study, J. Comput. Chem., 2009, 30(15), 2433–2442 CrossRef CAS PubMed.
- G. Vicario, G. Balducci, S. Fabris, S. de Gironcoli and S. Baroni, Interaction of hydrogen with cerium oxide surfaces: a quantum mechanical computational study, J. Phys. Chem. B, 2006, 110(39), 19380–19385 CrossRef CAS PubMed.
- M. Cococcioni and S. De Gironcoli, Linear response approach to the calculation of the effective interaction parameters in the LDA+ U method, Phys. Rev. B: Condens. Matter Mater. Phys., 2005, 71(3), 035105 CrossRef.
- W. Tang, E. Sanville and G. Henkelman, A grid-based Bader analysis algorithm without lattice bias, J. Phys.: Condens. Matter, 2009, 21(8), 084204 CrossRef CAS PubMed.
- J. Portier, G. Campet, J. Etourneau and B. Tanguy, A simple model for the estimation of electronegativities of cations in different electronic states and coordinations, J. Alloys Compd., 1994, 209(1–2), 285–289 CrossRef CAS.
- R. D. Shannon, Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides, Acta Crystallogr., Sect. A: Cryst. Phys., Diffr., Theor. Gen. Crystallogr., 1976, 32(5), 751–767 CrossRef.
- K. Momma and F. Izumi, VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data, J. Appl. Crystallogr., 2011, 44(6), 1272–1276 CrossRef CAS.
- H. T. Stokes and D. M. Hatch, FINDSYM: program for identifying the space-group symmetry of a crystal, J. Appl. Crystallogr., 2005, 38(1), 237–238 CrossRef CAS.
Footnote |
† Electronic supplementary information (ESI) available: The calculated total energies of the Ce16Ti2Zr14O64 systems for the symmetrically independent combinations of Ti substitutions. See https://doi.org/10.1039/d2ra02419e |
|
This journal is © The Royal Society of Chemistry 2022 |