research papers
Structure model of γ-Al2O3 based on planar defects
aInstitute of Materials Science, TU Bergakademie Freiberg, Gustav-Zeuner-Straße 5, D–09599 Freiberg, Germany
*Correspondence e-mail: m.rudolph@iww.tu-freiberg.de
The defect structure of γ-Al2O3 derived from boehmite was investigated using a combination of (SAED) and powder X-ray diffraction (XRD). Both methods confirmed a strong dependence of the diffraction line broadening on the diffraction indices known from literature. The analysis of the SAED patterns revealed that the dominant structure defects in the spinel-type γ-Al2O3 are antiphase boundaries located on the lattice planes , which produce the shifts . Quantitative information about the defect structure of γ-Al2O3 was obtained from the powder XRD patterns. This includes mainly the size of γ-Al2O3 crystallites and the density of planar defects. The correlation between the density of the planar defects and the presence of structural vacancies, which maintain the stoichiometry of the spinel-type γ-Al2O3, is discussed. A computer routine running on a fast graphical processing unit was written for simulation of the XRD patterns. This routine calculates the atomic positions for a given kind and density of planar defect, and simulates the diffracted intensities with the aid of the Debye scattering equation.
Keywords: γ-alumina; microstructure defects; antiphase boundaries; rotational boundaries; selected-area electron diffraction; powder X-ray diffraction; Debye equation; anisotropic broadening.
1. Introduction
γ-Alumina (γ-Al2O3) is one of the intermediate aluminium oxides that accompany the temperature-induced transition of boehmite [γ-AlO(OH)] towards the thermodynamically stable corundum (α-Al2O3), see Euzen et al. (2002):
Because of its large surface area and high γ-Al2O3 is often used as a functional constituent of catalytic converters or surface coatings. However, as γ-Al2O3 is only stable below 700oC (Paglia et al., 2004; Rudolph et al., 2017), it must be stabilized for high-temperature applications.
metastableAll intermediate alumina phases possess a slightly distorted cubic close-packed (c.c.p.) et al., 1982; Wilson, 1979). It survives the topotactic and is rearranged only by the formation of α-Al2O3 at approximately 1200oC (Euzen et al., 2002). Aluminium cations occupy octahedral and tetrahedral sites in the crystal structures of most intermediate alumina phases. However, in contrast to the fully occupied oxygen the aluminium contains structural vacancies, which balance the [Al]/[O] ratio. In consecutive intermediate alumina phases, these vacancies show different degrees of ordering that increases in general from γ-Al2O3 to θ-Al2O3.
of oxygen anions. This is already present in the octahedral double layers of boehmite (ChristensenThe different distribution of vacancies on the cation sites is one of the reasons why different structural descriptions of γ-Al2O3 can be found in the literature for differently prepared samples, for samples with different crystallite size etc. Thus, there is a need for a generalized structure model of γ-Al2O3, which would take into account the presence of structural vacancies with partially correlated positions and which would allow for quantitative description of these structure defects. Such a structure model is a first step towards understanding the successive phase transformation in metastable alumina phases and a prerequisite for targeted manipulation of their thermal stability.
Zhou & Snyder (1991) described the of γ-Al2O3 as a defective spinel structure. The term `defective' stands for a small tetragonal distortion of the cubic lattice, for small static displacements of atoms from their ideal positions, for the presence of vacancies in the cation and for the displacement of some cations into the non-spinel sites. In an ideal spinel structure with the and with the , the O2− anions occupy the Wyckoff positions 32e within the c.c.p. the A2+ cations occupy the tetrahedral sites 8a and the B3+ cations occupy the octahedral sites 16d (Sickafus et al., 2004). From the cation valency and from the cation-to-anion ratio (2/3) in Al2O3, it becomes evident that the trivalent Al3+ cations have to be located on both, tetrahedral and octahedral, sites. Moreover, a non-integer number of cation sites have to remain vacant to obey the stoichiometry of Al2O3.
In the structure models of γ-Al2O3 based on the analysis of integral intensities obtained from the X-ray or neutron powder diffraction patterns, the structural vacancies are typically randomly distributed over the standard tetrahedral and octahedral spinel sites (8a and 16d). In order to be able to reproduce the measured integral intensities more accurately, the sites 8b, 16c and 48f, which are empty in ideal spinel structures, were assumed to be partially occupied by Al3+ cations (Verwey, 1935; Ushakov & Moroz, 1984; Zhou & Snyder, 1991; Paglia et al., 2003; Smrčok et al., 2006). As a random distribution of vacancies at the cation sites can lead to an energetically unfavorable clustering of Al3+ cations (Cowley, 1953), specific ordering of vacancies is usually assumed when the local structure is investigated, e.g. in ab initio calculations (Digne et al., 2004; Menéndez-Proupin & Gutiérrez, 2005).
However, none of the above ). Alternative microstructure models of γ-Al2O3 are based on the assumption that the anisotropy of the line broadening observed is caused by planar defects, which are known to broaden the diffraction lines differently even for equivalent hkl (Warren, 1990; Guinier, 1994). Typical examples of the planar defects reported in conjunction with the microstructure models of γ-Al2O3 are antiphase boundaries (Dauger & Fargeot, 1983) and stacking faults (Cowley, 1953; Fadeeva et al., 1977; Paglia et al., 2006; Tsybulya & Kryukova, 2008). In these models, a coalescence of vacancies at the planar defects was assumed in order to avoid unwanted occupation of neighboring cation lattice sites (Kryukova et al., 2000).
models can explain the observed dependence of the X-ray diffraction (XRD) line broadening on the diffraction indices, which was already reported by Zhou & Snyder (1991Based on these considerations, Tsybulya & Kryukova (2008) and Pakharukova et al. (2017) developed a three-dimensional model of the real structure of γ-Al2O3, which consists of small unperturbed nano-sized domains of γ-alumina with spinel-like The nano-domains are terminated by the lattice planes {001}, {011} and {111} and possess polyhedral shapes. Neighboring nano-domains are mutually shifted along the gliding planes {001}, {011} and {111}. The shift vectors are summarized in Table 1.
|
The defects from Table 1 act as stacking faults for the cation but they keep the anion intact. Consequently, they do not cause any broadening of the diffraction lines, which stem predominantly from the scattering of X-rays on the oxygen whereas the diffraction lines produced by the scattering on the aluminium are strongly broadened. The capability of this model to explain the observed anisotropy of the XRD line broadening was recently illustrated by Pakharukova et al. (2017) for γ-Al2O3 derived from a boehmite-based aerogel precursor.
In the present study, we extend the microstructure models invented by Tsybulya & Kryukova (2008) and Pakharukova et al. (2017) by considering not only fractional lattice translations but also selected lattice rotations, and discuss the influence of resulting microstructure defects on (SAED) and XRD patterns. In all cases, only the planar defects are considered, which do not affect the atomic ordering within the oxygen sublattice.
The SAED patterns were simulated using JEMS (Stadelmann, 2012) and DIFFaX (Treacy et al., 1991). For simulation of the XRD patterns, a fast algorithm based on the Debye equation (Debye, 1915) was written for a graphical processing unit. The simulated diffraction patterns were compared with the diffraction patterns measured on γ-Al2O3 that was obtained by annealing of highly crystalline boehmite. As the γ-alumina derived from boehmite exhibits a small tetragonal distortion, the defect is described in I41/amd instead of (Paglia et al., 2003). For the orientation relationship between cubic and tetragonal γ-alumina,
the equivalent lattice planes and Wyckoff positions are summarized in Table 2.
|
2. Experimental
The γ-Al2O3 under study was prepared by heating highly crystalline boehmite powder (Actilox® B20, Nabaltec) in air for 20 h at 600oC. The heating rate was 2°C min−1. In the original state, the boehmite particles had the shape of thin, almost rhombic platelets (Fig. 1). Their size varies between a few hundred nanometres and several micrometres (Rudolph et al., 2017). After annealing, this particle shape was preserved, but the corresponding crystallographic directions changed according to the orientation relationship
The SAED patterns were recorded using a JEM-2200 FS field-emission transmission electron microscope (TEM) that was operated at the acceleration voltage of 200 kV. The TEM was equipped with a high-resolution objective lens ( = 0.5 mm) and a corrector of the spherical aberration that was located in the primary beam. The XRD experiments were performed on a Bragg–Brentano diffractometer (URD6 from Seifert/FPM) using Cu Kα radiation ( = 1.5406 Å, = 1.5444 Å) and a The XRD patterns were collected in a 2θ range between 15 and 70° with a 0.02° step size and a counting time of 12 s per step. For calculation of the XRD patterns the same wavelengths and 2θ range were utilized.
3. Results
3.1. Models of planar defects in γ-Al2O3 obtained from analysis of the SAED pattern
SAED pattern of γ-Al2O3 derived from boehmite (Fig. 2) consists of narrow diffraction spots and streaks, which are elongated in the c* direction. This kind of line broadening indicates the presence of planar defects on the lattice planes , which was already described by Cowley (1953). The observed dependence of the line broadening on the diffraction indices confirms different degrees of disorder in individual sublattices. As it can be seen from equation (7), the strongly elongated diffraction spots (111, 113, 220, 224, 331, 333 etc.) stem from the disordered cation while the narrow diffraction spots (004, 222, 440, 444 etc.) are dominated by the scattering on a fully occupied and well ordered c.c.p. anion (Lippens & de Boer, 1964; Paglia et al., 2004; Tsybulya & Kryukova, 2008).
The difference in shape of the elongated diffraction spots with even (220, 224 etc.) and odd (111, 113, 331, 333 etc.) indices stems from different arrangements of Al3+ cations on regular (8a and 16d) and irregular (8b, 16c and 48f) Wyckoff sites. It follows from Table 4 and equation (7) that the structure factors of the reflections with odd diffraction indices depend on the occupancy of the regular Wyckoff positions 8a, 16d and the irregular Wyckoff positions 8b and 16c, while the structure factors of the reflections with even diffraction indices depend on the occupancy of the regular 8a and on the occupancy of the irregular positions 8b and 48f.
The observed dependence of the line broadening on diffraction indices (Fig. 2) confirms that in γ-Al2O3 derived from highly crystalline boehmite, the anion remains intact, while the cation is disordered in such a way that the tetrahedral sites 8a, 8b and 48f with x = 0.375 are more disordered than the octahedral positions 16d and 16c, as already proposed by Zhou & Snyder (1991) and Paglia et al. (2004).
In our microstructure models derived from non-distorted γ-Al2O3 (Fig. 3a), the defect structures were produced by antiphase boundaries (APBs). The examples of the microstructure models include simple glide (Figs. 3b and 3f), which is connected with a shift of the atoms within the glide plane, APBs producing out-of-plane shifts of the atoms (Figs. 3c and 3e) and rotational boundaries (Fig. 3d). The simple glide, which corresponds to the model suggested by Tsybulya & Kryukova (2008), produces conservative APBs. The out-of-plane shift, which corresponds to the model proposed by Dauger & Fargeot (1983), produces non-conservative APBs. Rotational boundaries (RBs) can produce either conservative or non-conservative APBs. The model presented in Fig. 3(d) corresponds to non-conservative APBs. Conservative APBs do not change the local chemical composition, while some of the non-conservative APBs change the local stoichiometry of Al2O3. The SAED patterns that correspond to the models from Figs. 3(a)–3(f) are shown in Figs. 4(a)–4(f).
The SAED patterns were simulated using the computer program JEMS (Stadelmann, 2012) for orthorhombic supercells, which contained six APBs of the respective kind. The size of the respective was about 5.6 Å along the electron beam and approximately 50.5 ×46.9 Å2 in the in-plane directions [110] and [001], respectively. Thus, the mean distances between the planar defects were approximately 6.7 Å for the APBs from Figs. 3(b)–3(d), 7.2 Å for APBs from Fig. 3(e) and 9.9 Å for APBs from Figs. 3(e) and 3(f). Such a high defect density suppressed automatically the formation of reflections.
For simulations of the SAED patterns, the specific APBs from Figs. 3(a)–3(f) were complemented by the antiphase boundaries with crystallographically equivalent shifts and gliding or rotation planes. A full set of APBs that belong to the APB system is presented as an example in Table 5. The positions of individual atoms within the were generated as described in Appendix A.
The simulations of the SAED patterns (Figs. 4c and 4d) confirmed that only non-conservative APBs with the gliding planes can produce the streaks in the c* direction and the anisotropic broadening of the diffraction spots, which were observed in the experimental SAED pattern (Fig. 4g). Conservative APBs (Fig. 4b) can be excluded as a possible source of the observed line broadening, as these planar defects do not affect the shape of the reflections with even diffraction indices (220, 224 etc.).
The APBs located on the lattice planes (hh0), (hhh) and broaden the diffraction spots along the normal direction to the respective lattice plane (Figs. 4e and 4f), thus the presence of these defects can be excluded as well. Analogously, planar defects on other lattice planes can also be eliminated, as they would not broaden the diffraction spots in the c* direction.
3.2. Possible models of planar defects in γ-Al2O3 from the point of view of powder XRD
Although the structural information contained in the SAED pattern of a single crystal is much more comprehensive than the information contained in the powder XRD pattern, the powder XRD patterns were additionally consulted in order to obtain statistically relevant information about the density of microstructure defects. Furthermore, powder XRD provides direct access to all reflections in hhl (in cubic notation).
in contrast to the SAED experiment. Without laborious sample preparation, the plate-like particles were transparent in the electron beam only in the vicinity of the direction, because the facets of the particles have a specific orientation with respect to the crystallographic axes. Consequently, our SAED experiments were limited to the diffraction spotsThe XRD pattern measured in the powder sample of γ-Al2O3 is displayed in Fig. 5 together with the XRD patterns simulated for the APBs discussed above. The diffraction patterns in Figs. 5(a), 5(b), 5(c), 5(e) and 5(f) were simulated using the same structure models as the SAED patterns in Figs. 4(a), 4(b), 4(c), 4(e) and 4(f). Additional XRD patterns were simulated for conservative APBs , and suggested by Tsybulya & Kryukova (2008), and for non-conservative APBs . In analogy with the SAED simulations, the powder XRD patterns were simulated by taking into account the crystallographically equivalent APBs. For the APB system , the distribution of the APBs is shown as an example in Fig. 8 and the individual shift vectors are listed in Table 5.
It follows from the comparison of simulated and measured XRD patterns (Fig. 5) that conservative APBs , and can be excluded from further considerations, because these defects do not reproduce the experimentally observed broadening of diffraction lines 220, 422 and/or 511/333 satisfactorily. On the contrary, the conservative APBs and as well as the non-conservative APBs , and reproduce the observed anisotropic line broadening quite well. The simulated XRD patterns show differences mainly in the intensities of the diffraction lines 111, 311 and 222. The best agreement was achieved for the conservative APBs . The simulations of the powder diffraction patterns for rotational boundaries (RBs) (Fig. 6) revealed that RB { h00} is the most probable defect of this kind.
However, it must be kept in mind that the APBs , and were excluded by SAED simulations, as they broaden the diffraction spots in a wrong direction (Figs. 4e and 4f) or as they do not broaden the diffraction spots 220, 224 etc. (Fig. 4b). Thus, the only promising candidates for description of the planar defects in γ-Al2O3 derived from boehmite are the non-conservative APBs and the RBs { h00}.
3.3. Quantitative description of planar defects in γ-Al2O3 based on the analysis of the XRD patterns
The powder XRD patterns in Fig. 5 were simulated for cuboidal crystallites with the edge length of 10.8 nm, which contain 24 equally distributed planar defects of the respective kind. The mean distances between the defects were approximately 1.2, 3.1 and 2.7 nm, which corresponds to the planar defect densities of 16.7, 4.5 and 8.3% if the planar defects are located on the planes { 400}, { 440} and { 222}, respectively. The size of the cuboidal crystallites (10.8 nm) was determined from the broadening of the `narrow' diffraction lines 222, 400, 440 and 444.
For calculation of the XRD patterns, a fast computer routine (cuDebye) was written. It is based on the Debye equation (Debye, 1915) and runs on a graphic processing unit. Individual atomic positions are generated as described in Appendix A. The overall isotropic temperature factor [B in equation (9)] was ∼2 Å2. As an idealized structure was used for the simulations, this factor accounts mainly for static displacements in the structure of γ-Al2O3. The performance of the simulation routine can be illustrated by the short computing time, which was less than 5 min when calculating the XRD powder pattern with 6500 data points for a structure model consisting of 1.2×106 atoms. A drawback of this routine is that it produces intensity oscillations in the low-angle region ( 30°) and near the tails of the Bragg peaks (Dopita et al., 2013).
4. Consequences of the proposed structure model of γ-Al2O3
4.1. Effect of the parameters of the proposed structure model on the XRD patterns
From the point of view of the kinematical diffraction theory (Warren, 1990; Guinier, 1994), one reason for the anisotropic line broadening is a limited coherence of the in certain crystallographic directions, which leads to the broadening of the points along the corresponding directions in the In γ-Al2O3 derived from boehmite, the coherence of the is interrupted by planar defects like non-conservative APBs and RBs { h00}, as they modify the atomic ordering and fragment the aluminium into very small coherent domains. From the crystallographic point of view, the APBs and RBs in γ-Al2O3 produce disordered with broad diffraction maxima instead of distinct satellites. However, as the change in the atomic ordering concerns only the cation (Al3+) sites, the diffraction lines stemming from the scattering on oxygen anions are always unaffected.
The change in the atomic ordering is caused by the shift vectors belonging to the APBs. For instance, the non-conservative APBs shift some of the Al3+ cations from the regular positions 8a and 16d to 8b and 16c or even to 48f, which corresponds basically to the structure models published by Verwey (1935), Ushakov & Moroz (1984), Zhou & Snyder (1991), Paglia et al. (2003) and Smrčok et al. (2006). However, the crystallographic description of these structure defects implies equal occupancies of crystallographically equivalent atomic positions, while the APBs shift only some of the Al3+ cations from their regular positions to the irregular positions. The positions of oxygen anions that occupy the Wyckoff sites 32e remain unaffected in all cases.
Frequently, a sequence of APBs is required to achieve the specific shift vector. The shift of the tetrahedral cations from the Wyckoff sites 8a to 8b and the shift of the octahedral cations from the Wyckoff sites 16d to 16c are associated with one of the translation vectors , which result, for example, from the sum of the shift vectors and . The migration of tetrahedrally coordinated Al3+ cations from the Wyckoff sites 8a to some of the sites 48f is simply produced by a single non-conservative shift of the type , cf. Table 6. The shift of the Al3+ cations from the 8a sites to the remaining 48f sites is associated with one of the translation vectors that results, for instance, from the sum of the shift vectors and .
In comparison with the non-conservative shifts, the conservative shifts of the Al3+ cations produce less disorder in the The conservative counterpart of the APBs discussed above causes a shift of the Al3+ cations from the tetrahedrally coordinated Wyckoff sites 8a to 16 of the 48f sites only. The 32 atomic positions , , , , , , and , where z is an odd multiple of , are inaccessible (Fig. 3a). This example shows that the non-conservative shifts within the APB system { h00} move the Al3+ cations to the `non-spinel' positions with a higher probability than the conservative ones.
A consequence of the above atomic shifts is a change in the phase of the scattered wave that primarily modifies the
of a crystallite with APBs and consequently the intensities of diffraction lines. For APBs, the phase shift can be concluded easily from the scalar product of the shift vector of the respective APB with the diffraction vector. The change in the at the APB is given by the multiplicative factor:where h = (hkl) is the vector of diffraction indices and RT is the transposed shift vector. When , the shift vector does not change the phase of the at the APB. The domains separated by such APBs are fully coherent and the corresponding diffraction lines are not affected by these defects. The factor changes the sign of the `local' the magnitude of the calculated over the whole crystallite and finally the diffracted intensity.
An overview of the values calculated for heavily broadened diffraction lines 111, 220 and 311 and for possible shift vectors is given in Table 3. For the narrow diffraction lines 222, 004, 440 and 444, the multiplicative factors calculated for these shift vectors are always equal to unity. It also follows from equation (7) that the structure factors of the diffraction lines 222, 004, 440 and 444 are not affected by the displacement of the aluminium cations from the Wyckoff positions 16d to 16c and from the Wyckoff positions 8a to 8b or 48f.
|
As the phase shift and the fragmentation of the cation g and 4h). As can be seen from Table 3, all crystallographically equivalent shifts producing non-conservative APBs on the lattice planes yield negative for the reflection 220, while only half of the values are negative for the reflections 111 and 113.
are two concurrent consequences of the presence of specific APBs that limit the coherence of the from the point of view of the diffraction methods, the probability of the change in the multiplicative factor can be related to the broadening of individual diffraction lines. This feature can be illustrated in the different broadening of the reflections 220 and 113 (Figs. 4At equal probabilities of the crystallographically equivalent shifts, only half of the non-conservative APBs contribute to the broadening of the reflections 111 and 311, while more than half of the crystallographically equivalent shifts `broaden' the reflections of the type 220. The particular reflections 220 and are broadened by all non-conservative shifts of the APBs , see Table 3. In other words, the diffraction on the lattice planes { 111} and { 311} does not recognize half of the non-conservative APBs as domain boundaries, while the fraction of these APBs that are active as domain boundaries for the diffraction lines 220 is higher. This phenomenon is apparent in the measured and simulated SAED patterns (Figs. 4c, 4g and 4h), where the reflection 220 is approximately two times broader than the reflections 111 and 113.
In powder XRD patterns, the line broadening induced by APBs is additionally modified by the integration of the diffracted intensities over a certain reciprocal-space volume that is usually performed using the powder pattern power theorem (Warren, 1990). In the case of the APBs { h00}, which broaden the points along the a*, b* and c* directions, this integration makes strongly broadened diffraction lines asymmetric. Another phenomenon, which additionally affects the broadening and the shape of individual diffraction lines, is the partial coherence of neighboring nano-domains (Rafaja et al., 2004) that are separated from each other by APBs. Nevertheless, the powder pattern power theorem and the partial coherence of the nano-domains are considered a priori, if the Debye equation is utilized for calculation of the diffracted intensities.
The comparison of measured XRD patterns with the XRD patterns simulated for non-conservative APBs (Fig. 5c) and RBs { h00} (Fig. 6d), both having the mean distance of 1.2 nm, shows that the corresponding models can reproduce the defect structure of γ-Al2O3 derived from boehmite quite well. The only noteworthy discrepancy can be seen for the diffraction line 111, which is extremely weak in the simulated pattern. Partial extinction of this diffraction line is caused mainly by the almost equidistant distribution of APBs and RBs in the microstructure model used for the diffraction pattern simulation. Fragments of the separated by a planar defect with a phase change factor scatter with opposite phases, thus their structure factors (scattered amplitudes) are fully subtracted if these regions have the same size.
An irregular distance between planar defects with negative phase-change factors counteracts this extinction and consequently increases the diffracted intensity. This effect is further enhanced by the partial coherence of neighboring regions, which is the strongest for short diffraction vectors and becomes weaker with increasing magnitude of the diffraction vector (Rafaja et al., 2004). In this context, it is worth noting that the sequence of the phase-change factors is different for 111 and 311 (Table 3) and that the magnitude of the diffraction vector is almost twice as large for the diffraction line 311 than for 111, as it follows from . From this point of view, the intensity of the diffraction line 111 is much more sensitive to the regularity of the APB or RB distribution than the intensity of the diffraction line 311.
4.2. Influence of the proposed antiphase boundaries on the local stoichiometry of γ-Al2O3
From the crystallographic point of view, γ-Al2O3 crystallizes in a defect spinel structure, in which the Wyckoff positions 32e are fully occupied by oxygen and the Wyckoff positions 8a and 16d partially occupied by aluminium. The partial occupation of the lattice sites is needed to retain the Al2O3 stoichiometry. On the local scale, where microstructure defects like APBs and RBs are considered, the lattice sites cannot be occupied partially. They can either be fully occupied or stay vacant.
In the APB model discussed above, the formation of `virtual' vacancies is a necessary consequence of the local lattice shift. As illustrated in Fig. 7, the non-conservative APB moves the tetrahedral cations located at the Wyckoff positions 8a near the antiphase boundary very close to the octahedral cations 16d of the original Thus, the distance between both kinds of the cations would be solely , where a is the lattice parameter of the cubic structure.
In order to avoid such atomic proximity, one of the two close cations located at the APB must be removed (Cowley, 1953; Dauger & Fargeot, 1983; Tsybulya & Kryukova, 2008). The removal of the cations located within the APBs implies a clustering of vacancies in the vicinity of these defects. Based on the ab inito calculations performed by Wolverton & Hass (2000), the cations at the octahedral sites are preferentially replaced by vacancies in our model. This replacement produces two vacancies per unit area of the APB ( a2). In Fig. 7, these vacancies are located on a (001) plane. In a general case, two octahedral sites on an { h00} plane remain vacant.
The number of vacancies related to a single APB can be used to estimate the number of APBs () in a single crystallite that are necessary to maintain the stoichiometry of γ-Al2O3. In the defect spinel structure of γ-alumina, the stoichiometry ratio [Al]/[O] is given by
where nV (= 2) is the number of vacancies per unit area of APB. N3 is the number of elementary cells in a cuboidal γ-Al2O3 crystallite. is the number of the Wyckoff positions 8a and 16d that accommodate Al3+ cations and is the number of the Wyckoff positions 32e that are occupied by the O2− anions. It follows from equation (5) that the number of APBs scales with the linear size of the crystallite (N) as
The simulations of the SAED and XRD patterns were performed with cuboidal γ-Al2O3 crystallites that had the edge length of 10.8 nm and that consisted of 13 to 14 cubic elementary cells with the lattice parameter . Consequently, the volume of a single crystallite contained about 18 non-conservative APBs .
If the APBs would be equally distributed on all crystallographically equivalent lattice planes { h00}, i.e. on (h00), (0k0) and , their mean distances would be equal to 1.53 nm. The mean distances of the APBs located only on the lattice planes would be 0.56 nm. The theoretical mean distance between the APBs distributed over the crystallographically equivalent lattice planes { h00}, i.e. 1.53 nm, is very close to the experimentally determined mean distance of the APBs from Section 3.3, which was about 1.2 nm. This comparison indicates that the APBs are located with almost the same probability on all crystallographically equivalent lattice planes { h00}.
Finally, it should be noted that the mean distance of APBs is slightly affected by the size of the cuboidal crystallites. As the APBs cannot be located on the crystallite surface in the model their distances become smaller with decreasing crystallite size although their density (number per crystallite volume) is kept constant. For an infinitely large crystallite, the mean distance of APBs distributed over the lattice planes { h00} and approaches and , respectively.
4.3. γ-Al2O3 δ-Al2O3
With increasing temperature and consequently higher cation mobility, the cation vacancies in γ-Al2O3 start to rearrange, which leads to the formation of a periodic sequence of APBs that is interpreted as formation of the δ-Al2O3 (Dauger & Fargeot, 1983). The periodic ordering of cation vacancies in δ-Al2O3 leads to the formation of reflections, which gradually replace the strongly broadened reflections in γ-Al2O3. The first indications of the formation are already visible in Fig. 2. The progress of the γ-Al2O3 δ-Al2O3, the formation of several intermediate states between the metastable phases γ-Al2O3 and δ-Al2O3, and the completed formation of δ-Al2O3 at 975oC were reported by Rudolph et al. (2017).
An indirect consequence of our model of γ-Al2O3 is that the atomic ordering in δ-Al2O3 (as well as the and the size of the δ-Al2O3 unit cell) may depend on the density and distribution of the planar defects in original γ-Al2O3. As already pointed out by Tsybulya & Kryukova (2008) and Pakharukova et al. (2017), the distribution of planar defects in γ-Al2O3 can depend on the synthesis procedure. Therefore, it is not surprising that different unit cells were found for plasma sprayed δ-Al2O3 and for δ-Al2O3 derived from boehmite (Levin & Brandon, 1998).
5. Conclusions
Based on the results of γ-Al2O3 was proposed. It was shown that in γ-Al2O3 prepared from boehmite, the dominant defects are non-conservative antiphase boundaries and rotational boundaries { h00}, which affect the occupancy of the cation sites in the spinel-like structure of γ-Al2O3. Both kinds of defect disrupt the coherence of the cation for diffraction and consequently broaden certain diffraction lines. Furthermore, it was shown that the antiphase boundaries and rotational boundaries { h00} induce a shift in the Al3+ cations located in the nano-domains from the Wyckoff sites 8a and 16d to some of the Wyckoff sites 8b, 16c and 48f, while the positions of the O2− anions remain unaffected.
and powder X-ray diffraction, a defect-structure model ofIn contrast to the previously published b, 16c and 48f are occupied statistically, the presented model occupies the Wyckoff positions 8b, 16c and 48f non-uniformly.
models, in which the Wyckoff sites 8The shift of atoms caused by the planar defects leads to a convergence of the cations in the vicinity of the planar defects, which is avoided by replacing half of the unfavorably coordinated cations with structural vacancies. The presence of these vacancies retains the stoichiometry of Al2O3. For quantification of these planar defects, a computer routine based on the Debye scattering equation was written. The comparison of measured and simulated X-ray diffraction patterns revealed that the relevant planar defects in γ-Al2O3 obtained from boehmite annealed for 20 h at 600oC have distances of approximately 1.2 nm and that they are located on all crystallographically equivalent lattice planes (h00), (0k0) and .
APPENDIX A
Modeling of planar defects and simulation of X-ray diffraction patterns
The presented structure model of γ-Al2O3 is based on an idealized cubic spinel structure with the , in which the Wyckoff positions 8a and 16d are occupied by aluminium cations and the Wyckoff positions 32e by oxygen anions. The fractional coordinates of the idealized structure are given in Table 4. The partial occupancy of the 16d balances the stoichiometry of Al2O3.
|
The Wyckoff positions 8b, 16c and 48f can be found with a marginal occupancy in several structure reports (Verwey, 1935; Ushakov & Moroz, 1984; Zhou & Snyder, 1991; Paglia et al., 2003; Smrčok et al., 2006). Partial occupancy of these Wyckoff positions modifies the structure factors (and intensities) of the diffraction lines with odd indices, and the structure factors and intensities of the diffraction lines 220, 422, 442 etc. see the following equation:
In this equation, and are the atomic scattering factors of aluminium and oxygen, and O is the occupancy of the corresponding Wyckoff site. For the modeling of the planar defects, a Matlab routine was written. In this routine, the planar defects are generated on the respective lattice planes i.e. { h00}, { hh0} or { hhh}. The planar defects are distributed randomly with the desired probability over the respective lattice planes, as it is illustrated in Fig. 7 for lattice planes { h00}. In the next step, the respective shift vectors were applied, as shown in the example in Fig. 8 for conservative antiphase boundaries . In the case of non-conservative shifts, the overlapping parts of the structure are considered only once.
The APBs from Fig. 8 move the aluminium cations from the Wyckoff positions 16d (small filled circles) to the Wyckoff positions 16c (empty positions between the oxygen anions). The Wyckoff positions 32e occupied by oxygen anions (large open circles) are unaffected by all planar defects under consideration.
The APBs were chosen, because they can be depicted quite demonstratively in 2D. The entire set of the shift vectors belonging to the respective lattice plane is shown in Table 5. The kind of APB (conservative or non-conservative) is distinguished by the respective upper index. The APBs located on other lattice planes are classified in Table 6.
|
|
Successive application of such γ-Al2O3 crystallites into small domains (thin lines in Fig. 9), which can scatter both in the same phase or in the antiphase, as discussed above and illustrated in Table 3. When the rotational boundaries are modeled, the lattice planes defining the positions of these planar defects are set up in a similar way to the antiphase boundaries. However, the respective rotation axis, i.e. two-, three- and fourfold for the RBs { hh0}, { hhh} and { h00}, is applied instead of the shift vector.
translations tessellates the original | Figure 9 of |
The atomic positions within the nanocrystallites were simulated for the crystallite size of 10.8 nm, which was determined from the width of the `narrow' diffraction lines 222, 400, 440 and 444. This size corresponds to the size of the coherent oxygen a = 7.942 Å, which was afterwards tetragonally distorted in order to arrive at the measured c/a ratio of 0.985.
For the typical defect densities, such a crystallite contains between 100 and 1000 nano-domains of different shape and size. The planar defects were modeled for a cubic structure with the elementary cell size ofThe XRD patterns were calculated with the aid of the general scattering equation (Debye, 1915)
which intrinsically accounts for the powder pattern power theorem and for the partial coherence of nano-domains in a nanocrystallite. In equation (8), q is the magnitude of the diffraction vector (), is the diffraction angle and λ is the wavelength of the X-rays. rij denotes the distances between the atoms i and j within the crystallite. The factor Cij contains the atomic scattering factors f, the occupancies O and the temperature factors B for each atomic pair
The temperature factors were assumed to be isotropic and equal for all atoms. For comparison with the measured XRD patterns, the intensities calculated using equation (8) were corrected for polarization. Other effects including instrumental broadening were neglected (Dopita et al., 2015).
In order to account for a random distribution of the planar defects in powder samples, in which several crystallites can diffract concurrently, the planar defects in a crystallite were generated ten times at random positions but with the same type of defect and with the same defect density, and summed. The resulting simulated intensities were averaged, before they were compared with the experimental data.
The intensity calculation was speeded up by dividing the atoms in equation (8) into subgroups with the same properties ( f,O,B) and distances ( rk):
nk means the frequency of the distance rk. Moreover, as the oxygen remains unchanged in the faulted structure of γ-Al2O3, only the contribution of the aluminium to the diffracted intensities has to be recalculated for each kind and density of microstructure defect, when this procedure is utilized.
Consequently, the most time-consuming step is the calculation of the distribution of the atomic distances. However, when running on a CUDA capable graphic processing unit (GTX 980 TI from NVIDIA), our C++/Cuda routine cuDebye needs less than 5 min to simulate the X-ray powder diffraction pattern from nanocrystallites containing approximately 1.2 ×106 atoms. The experimental code of cuDebye that was written in C++ using Visual Studio Community 2015 and the CUDA Toolkit 8.0 is deposited at GitHub (https://github.com/Martin-Rudolph/cuDebye).
Supporting information
Information about the supporting structure and video files for Figure 2. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup1.txt
Video file of Vesta model for Figure 3a. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup2.mp4
Video file of Vesta model for Figure 3b. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup3.mp4
Video file of Vesta model for Figure 3c. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup4.mp4
Video file of Vesta model for Figure 3d-RB. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup5.mp4
Video file of Vesta model for Figure 3d-APB. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup6.mp4
Video file of Vesta model for Figure 3e. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup7.mp4
Video file of Vesta model for Figure 3f. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup8.mp4
Structure input file for Vesta for Figure 3a. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup9.txt
Structure input file for Vesta for Figure 3b. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup10.txt
Structure input file for Vesta for Figure 3c. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup11.txt
Structure input file for Vesta for Figure 3d-RB. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup12.txt
Structure input file for Vesta for Figure 3d-APB. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup13.txt
Structure input file for Vesta for Figure 3e. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup14.txt
Structure input file for Vesta for Figure 3f. DOI: https://doi.org/10.1107/S2052252518015786/ct5007sup15.txt
Acknowledgements
The authors would like to thank Professor Andreas Leineweber and Marius Wetzel for fruitful discussions.
Funding information
This work was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project Number 169148856 – Sonderforschungsbereich (SFB, Collaborative Research Centre) 920.
References
Christensen, A. N., Lehmann, M. S. & Convert, P. (1982). Acta Chem. Scand. 36a, 303–308. CrossRef Google Scholar
Cowley, J. M. (1953). Acta Cryst. 6, 53–54. CrossRef IUCr Journals Google Scholar
Dauger, A. & Fargeot, D. (1983). Radiat. Eff. 74, 279–289. CrossRef CAS Google Scholar
Debye, P. (1915). Ann. Phys. 351, 809–823. CrossRef Google Scholar
Digne, M., Sautet, P., Raybaud, P., Euzen, P. & Toulhoat, H. (2004). J. Catal. 226, 54–68. CrossRef CAS Google Scholar
Dopita, M., Emmel, M., Salomon, A., Rudolph, M., Matěj, Z., Aneziris, C. G. & Rafaja, D. (2015). Carbon, 81, 272–283. CrossRef CAS Google Scholar
Dopita, M., Rudolph, M., Salomon, A., Emmel, M., Aneziris, C. G. & Rafaja, D. (2013). Adv. Eng. Mater. 15, 1280–1291. Web of Science CrossRef CAS Google Scholar
Euzen, P., Raybaud, P., Krokidis, X., Toulhoat, H., Le Loarer, J.-L., Jolivet, J.-P. & Froidefond, C. (2002). In Handbook of Porous Solids, Vol. 3, pp. 1591–1677. Weinheim: Wiley-VCH Verlag GmbH. Google Scholar
Fadeeva, V. I., Kagan, A. S. & Zevin, L. S. (1977). Acta Cryst. A33, 386–389. CrossRef CAS IUCr Journals Google Scholar
Guinier, A. (1994). X-ray Diffraction In Crystals, Imperfect Crystals, and Amorphous Bodies, pp. 209–237, 270–274. New York: Dover Publications. Google Scholar
Kryukova, G., Klenov, D., Ivanova, A. & Tsybulya, S. (2000). J. Eur. Ceram. Soc. 20, 1187–1189. CrossRef CAS Google Scholar
Levin, I. & Brandon, D. (1998). J. Am. Ceram. Soc. 81, 1995–2012. CrossRef CAS Google Scholar
Lippens, B. C. & de Boer, J. H. (1964). Acta Cryst. 17, 1312–1321. CrossRef IUCr Journals Web of Science Google Scholar
Menéndez-Proupin, E. & Gutiérrez, G. (2005). Phys. Rev. B, 72, 035116. Google Scholar
Paglia, G., Božin, E. S. & Billinge, S. J. L. (2006). Chem. Mater. 18, 3242–3248. Web of Science CrossRef ICSD CAS Google Scholar
Paglia, G., Buckley, C. E., Rohl, A. L., Hunter, B. A., Hart, R. D., Hanna, J. V. & Byrne, L. T. (2003). Phys. Rev. B, 68, 144110. Web of Science CrossRef ICSD Google Scholar
Paglia, G., Buckley, C. E., Rohl, A. L., Hart, R. D., Winter, K., Studer, A. J., Hunter, B. A. & Hanna, J. V. (2004). Chem. Mater. 16, 220–236. CrossRef CAS Google Scholar
Pakharukova, V., Yatsenko, D., Gerasimov, E. Y., Shalygin, A., Martyanov, O. & Tsybulya, S. (2017). J. Solid State Chem. 246, 284–292. CrossRef CAS Google Scholar
Rafaja, D., Klemm, V., Schreiber, G., Knapp, M. & Kužel, R. (2004). J. Appl. Cryst. 37, 613–620. Web of Science CrossRef CAS IUCr Journals Google Scholar
Rudolph, M., Salomon, A., Schmidt, A., Motylenko, M., Zienert, T., Stöcker, H., Himcinschi, C., Amirkhanyan, L., Kortus, J., Aneziris, C. G. & Rafaja, D. (2017). Adv. Eng. Mater. 19, 1700141. CrossRef Google Scholar
Sickafus, K. E., Wills, J. M. & Grimes, N. W. (2004). J. Am. Ceram. Soc. 82, 3279–3292. CrossRef Google Scholar
Smrčok, L., Langer, V. & Křesťan, J. (2006). Acta Cryst. C62, i83–i84. CrossRef IUCr Journals Google Scholar
Stadelmann, (2012). Jems v3.76. Saas-Fee, Switzerland. Google Scholar
Treacy, M. M. J., Newsam, J. M. & Deem, M. W. (1991). Proc. R. Soc. London Ser. A, 433, 499–520. CrossRef Web of Science Google Scholar
Tsybulya, S. V. & Kryukova, G. N. (2008). Phys. Rev. B, 77, 024112. CrossRef Google Scholar
Ushakov, V. A. & Moroz, E. M. (1984). React. Kinet. Catal. Lett. 24, 113–118. CrossRef CAS Google Scholar
Verwey, E. (1935). Z. Kristallogr. Cryst. Mater. 91, 65–69. CAS Google Scholar
Warren, B. E. (1990). X-ray Diffraction, pp. 31–35, 216–227, 262–298. New York: Dover Publications. Google Scholar
Wilson, S. J. (1979). J. Solid State Chem. 30, 247–255. CrossRef CAS Google Scholar
Wolverton, C. & Hass, K. C. (2000). Phys. Rev. B, 63, 024102. CrossRef Google Scholar
Zhou, R. S. & Snyder, R. L. (1991). Acta Cryst. B47, 617–630. CrossRef CAS IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.