Oxygen vacancy (VO), one of the representative native defects in oxides, has received much attention because of the important role played thereby in determining the physical properties of materials [1-3]. Various tools have been used for the qualitative and quantitative analysis of VO. These include photoluminescence , ultraviolet-visible (UV-vis) , Raman , and X-ray photoelectron spectroscopic (XPS) techniques . Quantitative analyses of VO have been successfully performed via Rutherford backscattering spectroscopy  or X-ray absorption spectroscopy using synchrotron radiation . Additionally, X-ray diffraction (XRD) is a simple and useful tool for analysis of VO because it reveals the crystal structure and the electron density distribution of periodic arrays of atoms .
Fitting of X-ray diffraction data using the Rietveld refinement has been attempted for the quantitative analyses of VO in terms of oxygen site occupancy [10,11]. However, additional corrections and the use of neutron or synchrotron X-rays were required [11-13]. Electron density profiling using the maximum entropy method (MEM) is also a suitable tool for analysis of VO because it uses the more precise Rietveld refinement that resolves summation-terminated errors and affords a better structural model [14,15]. Furthermore, MEM introduces negligible modeling errors via least-biased electronic reconstruction of X-ray diffraction patterns in real space .
We sought to confirm whether MEM analysis could be used for analysis of VO. In this study, we applied such analysis to VO that changed as a function of the Co-doping concentration in ZnO. Co-doped ZnO is a good candidate room-temperature magnetic semiconductor and has been intensively studied in terms of intrinsic ferromagnetism. However, the origin of such ferromagnetism remains controversial, and the presence thereof limits applications of the semiconductor. VO was regarded, until recently, to explain the ferromagnetism and was reported to be affected by Co-doping concentration . Herein, we analyzed the change in VO as a function of Co-doping concentration and compared the results with XPS data. A method of analysis of VO is proposed, using conventional XRD and MEM techniques.
ZnO and Co-doped ZnO (Zn1−xCoxO, x = 0.01, 0.05, 0.1) powder samples were fabricated by sol-gel methods [16,17]. Zinc acetate dihydrate (Sigma-Aldrich, St. Louis, MO, USA) and cobalt acetate tetrahydrate (Sigma-Aldrich, St. Louis, MO, USA), used as starting materials, were dissolved in 2-methoxyethanol (Sigma-Aldrich, St. Louis, MO, USA) and stabilized by monoethanolamine (Sigma-Aldrich, St. Louis, MO, USA). To exclude the possibility of external contamination, the dissolution and drying processes were performed under a pure argon gas (99.999% purity) atmosphere, and under vacuum, each for 10 h, respectively. The organic residuals in samples were completely removed via an intermediate heat treatment at 300°C and a subsequent final heat treatment at 800°C under vacuum for 10 h . The samples used in this study were characterized using XRD, and we found a high degree of crystallinity, which was comparable to that of commercially available high-quality powder samples (ZnO; CAS 1314-13-2, Sigma-Aldrich, St. Louis, MO, USA). We also characterized the samples using synchrotron radiation, and we found high sample quality. A characterization study using synchrotron radiation will be submitted to a specialized journal soon. XRD (Empyrean Series 2, PANalytical) experiments were performed to analyze the crystal structures and electron density distributions of the powder samples. The Rietveld and MEM analyses were performed using a published technique . The MEM calculation was performed using ENIGMA software  with 66 × 66 × 104 pixels. The electron density distribution was reconstructed using the VESTA visualization program . An X-ray photoelectron spectrometer (model: Theta Probe (Thermo Electron Co., Waltham, MA, USA), Korean Basic Science Institute, Busan Center) was used for atomic composition analysis. Magnetic-field-dependent magnetization was measured using a vibrating sample magnetometer (VSM) equipped with a physical property measurement system (PPMS; Model 6000, Quantum Design, San Diego, CA, USA).
Results and discussion
Figure 1 shows the XRD patterns of the ZnO and ZnCoO powder samples. The Bragg peaks correspond to the wurtzite ZnO structure and other possible secondary species [16,21]. The diffraction intensity of each sample is presented logarithmically to better observe the background levels. Figure 1 shows that the XRD patterns of all samples had pronounced diffraction peaks that corresponded to wurtzite ZnO, with no additional peaks corresponding to secondary phases (at least within the detection limits). This indicates that the Co2+ ions were well-substituted during doping into Zn sites, without creation of secondary species [16,22]. Doping Co into ZnO tended to decrease the full-width at half-maximum (FWHM) values of XRD peaks. To quantitatively analyze this trend, the grain sizes of each powder sample were calculated using the Scherrer equation . These sizes were estimated to be 93, 107, 115, and 143 nm for ZnO, Zn0.99Co0.01O, Zn0.95Co0.05O, and Zn0.9Co0.1O, respectively. Co doping thus improved the crystallinity of the samples, in proportional to the Co-doping concentration, in our experimental range . The crystallinity of ZnO can be greatly affected by lattice strain  or native defects such as VO [26-28] and interstitial zinc (Zni) [27,28]. The contribution of lattice strain (resulting from the Co-doping concentration) to differences in crystallinity was assumed to be negligible because all powder samples were post-annealed at both 300°C and 800°C [17,29], which would relieve any lattice strain. Co doping was thought to induce insignificant amounts of lattice strain that might result from a difference in ionic radii, because the ionic radius of Zn2+ (74 pm) is similar to that of Co2+ (72 pm) . Consequently, the observed enhancement in the crystallinity of the samples was attributed to a decrease in the density of crystal defects caused by Co doping. A quantitative analysis of VO and Zni was performed to explore this hypothesis; we refined the XRD patterns of the samples.
Figure 2 shows the Rietveld refinement results for the ZnO and ZnCoO samples performed using single-phase wurtzite ZnO. The reliability factors based on the Bragg intensity (R I) and structure factor (R F) were below 1.5% and 0.8%, respectively, for all samples (Table 1). The refinement results revealed that the lattice constants steadily increased with increasing Co-doping level. It is known that formation of Zni leads to an increase in lattice constants , while formation of VO induces a decrease in lattice constants . The observed reduction of internal defects, and the increased lattice constants with increasing Co-doping level, indicated that the main defects embedded in the ZnO sample were VO, and hence, their density decreased as the Co-doping concentration increased.
An additional Rietveld refinement of the oxygen occupancies in each finalized Rietveld refinement was performed to examine the validity of this concept. Table 1 also lists the change in the oxygen occupancy [Δg O = (g O(ZnO) − g O(ZnCoO))/g O(ZnO)] for each sample, where g O refers to oxygen site occupancy. The oxygen occupancy increased for Zn0.99Co0.01O (i.e., negative Δg O was obtained) but decreased for Zn0.95Co0.05O and Zn0.9Co0.1O. Also, the oxygen vacancy was greater for Zn0.95Co0.05O than for Zn0.9Co0.1O. MEM analysis was used to further examine this trend.
It was expected that the presence of VO in ZnO and ZnCoO would induce changes in the electron density distributions of the oxygen atoms. Hence, the electron densities of the oxygen atoms as a function of Co-doping concentration were investigated via MEM/Rietveld analysis . Figure 3a,b show electron density maps of the ZnO and Zn0.9Co0.1O samples in the (110) plane, respectively. The Zn sites of the ZnO sample and the Zn(Co) sites of the Zn0.9Co0.1O sample did not noticeably differ in electron density. This is because of the low-dopant concentration used, and because Zn and Co exhibited similar electron density distributions, as Zn2+ and Co2+ have similar total numbers of electrons [33,34]. However, the electron densities at the central O atoms of ZnO and Zn0.9Co0.1O were clearly different, thus, 20.65 and 21.91 e/A3, respectively. This indicated that the oxygen sites in the wurtzite ZnO structure became increasingly occupied by oxygen atoms after Co doping. Thus, Co doping decreased the VO content.
Figure 3c shows the electron density line profiles along the direction of the O-Zn bond. These profiles enable precise analysis of oxygen occupancy as a function of the Co content. The lines were normalized to the electron density at the Zn(Co) atomic position to allow comparison of VO with Zn occupancy. The electron density at the O atomic position increased in the order Zn0.99Co0.01O < ZnO < Zn0.95Co0.05O < Zn0.9Co0.1O, in agreement with the Rietveld refinement results. The sample with 1% Co doping exhibited significantly lower oxygen electron density, which did not agree with the Rietveld refinement data.
XPS was used to rationalize the conflicting VO results. Figure 4 shows the percentage change in oxygen occupancy as a function of Co concentration as determined by the Rietveld refinement, MEM, and XPS techniques. The MEM data were calculated from the integrals of the electron density profiles of the Zn(Co) and O atom positions shown in Figure 3c and the ratio of the oxygen electron density area to the Zn(Co) electron density area. The change in oxygen occupancy revealed by XPS was obtained from the integrals of the areas of the Zn 2p, Co 2p, and the O 1s peaks. The calculation of change in oxygen occupancy was performed on the basis of a pure ZnO sample without Co dopant, as per the data of Table 1. The MEM and XPS results both showed that the sample with 1 mol% Co doping exhibited an abrupt drop in oxygen content; the samples with more doping had increased oxygen contents, in agreement with both the Rietveld refinement and the MEM results. Consequently, the XPS profile results were in better agreement with the MEM results than with the data from the Rietveld refinement. This indicated that MEM analysis, via repeated error correction, provided more accurate structural information. These findings suggest that reliable defect analysis may be possible using easily accessed laboratory X-ray data.
The magnetic field dependences of magnetization (M-H curves) were measured for all samples (Figure 5a). Pure ZnO is diamagnetic, and Co-doped ZnO exhibits paramagnetic behavior because of the 3d electron of Co2+. We reconfirmed that the ZnCoO samples were not intrinsically ferromagnetic, regardless of Co concentration [6,16]. Figure 5b shows the magnetic susceptibilities, which are the slopes of the M-H curves; these increased nonlinearly with increased Co-doping level. With increasing Co concentration, not all Co spins behave paramagnetically; some spins assume configurations differing in alignment. Considering the absence of secondary phases in the above structural analysis, we conjecture that increasing numbers of Co atoms assuming positions neighboring oxygen atoms created an antiferromagnetic configuration via superexchange interaction.
The observed trend, that creation of VO was suppressed with increasing Co-doping level, is attributable to differences in the Zn-O and Co-O bond strengths; the O2− ions in the wurtzite ZnO structure are tetrahedrally coordinated and thereby form four Zn-O bonds . Doping of Co2+ ions into ZnO creates Co-O bonds, the diatomic bond dissociation energy of which is higher than that of the Zn-O bond by 84 kJ/mol (Zn-O: 284 kJ/mol, Co-O: 368 kJ/mol) . This indicates that the Co-O bonds created by Co doping enhanced the average bond strength between oxygen ions and neighboring cations, i.e., doping decreased the possibility of oxygen-cation bond dissociation during sample fabrication or post-treatment processing . The experimental results indicate that VO decreased at high-level Co doping (i.e., above 5 mol%). The supporting analyses suggest that Co doping can impede creation of VO. However, the abrupt increase of VO at 1 mol% of Co doping is not well-understood and warrants additional study.
The Rietveld refinement results of the X-ray diffraction patterns of the ZnCoO system indicated that increased Co doping of ZnO tended to decrease the VO, but the VO increased slightly upon 1 mol% of Co doping. The MEM results were in better agreement with the XPS data, which indicated that MEM analysis could be a reliable tool for the study of VO. Additional research is needed to explain the anomalous behavior at 1 mol% of Co doping. More advanced X-ray electron density studies using synchrotron radiation would provide more precise and reliable data, but nevertheless, our present work shows that MEM is a reliable technique for the analysis of defects in materials characterized by XRD, which is a readily accessible tool in the material scientist laboratory. This approach will be of particular value in early-stage studies of oxide systems.
Ji Hun Park and Yeong Ju Lee contributed equally to this work.
Competing interests The authors declare that they have no competing interests.
Authors’ contributions JHP designed and planned the experiment. JHP and YJL carried out Rietveld refinement and MEM analysis. JSB contributed to experiments and analysis for XPS. BSK performed powder sample fabrication. YCC assisted with sample characterization and manuscript preparation. CM and YK contributed to MEM analysis. SL performed data analysis and guided this work. JHP, SL, and SYJ wrote the main manuscript with help from the co-authors. SYJ is responsible for project planning. All authors discussed the results and reviewed the manuscript. All authors read and approved the final manuscript.
Authors’ information JHP, BSK and YJL are graduate students of the Department of Cogno-Mechatronics Engineering, Pusan National University, Republic of Korea. JSB is a senior researcher of Korea Basic Science Institute (Busan Center), Republic of Korea. YCC is a research professor at the Crystal Bank Institute, Pusan National University, Republic of Korea. CM and YK are professors at the Department of Physical Science, Hiroshima University, Japan. SL, the corresponding author, is a research scholar at the Department of Materials Science and Engineering, University of Maryland. SYJ, the corresponding author, is a professor at the Department of Cogno-Mechatronics Engineering, Pusan National University, Republic of Korea.
This work was supported for 2 years by a Pusan National University Research Grant.
Janotti, A; Van de Walle, CG: Oxygen vacancies in ZnO, Appl Phys Lett., 2005, 87, S. 122102Crossref
Zhao, Q; Xu, XY; Song, XF; Zhang, XZ; Yu, DP: Enhanced field emission from ZnO nanorods via thermal annealing in oxygen, Appl Phys Lett., 2006, 88, S. 033102Crossref
Park, JH; Lee, S; Kim, B-S; Kim, W-K; Cho, YC; Oh, MW: Effects of Al doping on the magnetic properties of ZnCoO and ZnCoO:H, Appl Phys Lett., 2014, 104, S. 052412Crossref
Yang, H; Wang, YQ; Wang, H; Jia, QX: Oxygen concentration and its effect on the leakage current in BiFeO3 thin films, Appl Phys Lett., 2010, 96, S. 012909Crossref
Ciatto, G; Trolio, AD; Fonda, E; Alippi, P; Testa, AM; Bonapasta, AA: Evidence of cobalt-vacancy complexes in Zn1 − xCoxO dilute magnetic semiconductors, Phys Rev Lett., 2011, 107, S. 127206CrossrefPubMed
Charles K. Crystal structure. In: Stuart J, Patricia M, Martin B, editors. Introduction to solid state physics 8th edition. John Wiley & Sons, Inc; 2004. P. 3–22.Google Scholar
Gržetaa, B; Tkalčecb, E; Goebbertb, C; Takedac, M; Takahashic, M; Nomurad, K: Structural studies of nanocrystalline SnO2 doped with antimony: XRD and Mössbauer spectroscopy, J Phys Chem Solids., 2002, 63, S. 765–72Crossref
Dann, SE; Weller, MT: Structure and oxygen stoichiometry in Sr3Co2O7-y (0.94 ≤ y ≤ 1.22), J Solid State Chem, 1995, 115, S. 499–507Crossref
Yamazaki, S; Toraya, H: Rietveld refinement of site-occupancy parameters of Mg2-xMnxSiO4 using a new weight function in least-squares fitting, J Appl Cryst., 1999, 32, S. 51–9Crossref
Itoh, T; Nishida, Y; Tomita, A; Fujie, Y; Kitamura, N; Idemoto, Y: Determination of the crystal structure and charge density of (Ba0.5Sr0.5)(Co0.8Fe0.2)O2.33 by Rietveld refinement and maximum entropy method analysis, Solid State Commun., 2009, 149, S. 41–4Crossref
Kitaura, R; Kitagawa, S; Kubota, Y; Kobayashi, TC; Kindo, K; Mita, Y: Formation of a one-dimensional array of oxygen in a microporous metal-organic solid, Science., 2002, 298, S. 2358–61CrossrefPubMed
Takata, M; Nishibori, E; Sakata, M: Charge density studies utilizing powder diffraction and MEM. Exploring of high Tc superconductors, C60 superconductors and manganites, Z Kristallogr, 2001, 216, S. 71–86Crossref
Lee, S; Kim, B-S; Seo, S-W; Cho, YC; Kim, SK; Kim, JP: A study of the correlation between hydrogen content and magnetism in ZnCoO, J Appl Phys, 2012, 111, S. 07C304
Kim, SJ; Cha, SY; Kim, JY; Shin, JM; Cho, YC; Lee, S: Ferromagnetism in ZnCoO due to hydrogen-mediated Co − H − Co complexes: how to avoid the formation of Co metal clusters?, J Phys Chem C., 2012, 116, S. 12196–202Crossref
Tanaka, H; Takata, M; Nishibori, E; Kato, K; Iishi, T; Sakata, M: ENIGMA: maximum-entropy method program package for huge systems, J Appl Cryst., 2002, 35, S. 282–6Crossref
Momma, K; Izumi, F: VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data, J Appl Cryst., 2011, 44, S. 1272–6Crossref
Koch, U; Fojtik, A; Weller, IH; Henglein, A: Photochemistry of semiconductor colloids. Preparation of extremely small ZnO particles, fluorescence phenomena and size quantization effects, Chem Phys Lett, 1985, 122, S. 507–10Crossref
Yuhas, BD; Zitoun, DO; Pauzauskie, PJ; He, R; Yang, P: Transition-metal doped zinc oxide nanowires, Angew Chem., 2006, 118, S. 434–7Crossref
Patterson, AL: The Scherrer formula for X-ray particle size determination, Phys Rev., 1939, 56, S. 978–82Crossref
Kim, K-K; Song, J-H; Jung, H-J; Choi, W-K; Park, S-J; Song, J-H: The grain size effects on the photoluminescence of ZnO/α-Al2O3 grown by radio-frequency magnetron sputtering, J Appl Phys., 2000, 87, S. 3573–6Crossref
Yang, P; Yan, H; Mao, S; Russo, R; Johnson, J; Saykally, R: Controlled growth of ZnO nanowires and their optical properties, Adv Funct Mater., 2002, 12, S. 323–31Crossref
Kang, HS; Kang, JS; Kim, JW; Lee, SY: Annealing effect on the property of ultraviolet and green emissions of ZnO thin films, J Appl Phys., 2004, 95, S. 1246–50Crossref
Kim, K-K; Kim, H-S; Hwang, D-K; Lim, J-H; Park, S-J: Realization of p-type ZnO thin films via phosphorus doping and thermal activation of the dopant, Appl Phys Lett., 2003, 83, S. 63–5Crossref
Reynolds, DC; Look, DC; Jogai, B; Jones, RL; Litton, CW; Harsch, W: Optical properties of ZnO crystals containing internal strains, J Lumin., 1999, 82, S. 173–6Crossref
Yan, L; Ong, CK; Rao, XS: Magnetic order in Co-doped and „Mn, Co… codoped ZnO thin films by pulsed laser deposition, J Appl Phys., 2004, 96, S. 508–11Crossref
Khalid, M; Ziese, M; Setzer, A; Esquinazi, P; Lorenz, M; Hochmuth, H: Defect-induced magnetic order in pure ZnO films, Phys Rev B., 2009, 80, S. 035331Crossref
Lee, C-R; Wang, C-C; Chen, K-C; Lee, G-H; Wang, Y: Bond characterization of metal squarate complexes [MII(C4O4)(H2O)4; M = Fe, Co, Ni, Zn], J Phys Chem A., 1999, 103, S. 156–65Crossref
SenKov, ON; Miracle, DB: Effect of the atomic size distribution on glass forming ability of amorphous metallic alloys, Mater Res Bull., 2001, 36, S. 2183–98Crossref
Wang, ZL: Nanostructures of zinc oxide, Mater Today., 2004, 7, S. 26–33Crossref
Tan, K; Nijem, N; Canepa, P; Gong, Q; Li, J; Thonhauser, T: Stability and hydrolyzation of metal organic frameworks with paddle-wheel SBUs upon hydration, Chem Mater., 2012, 24, S. 3153–67Crossref