Print version ISSN 0103-9733
Braz. J. Phys. vol.40 no.3 São Paulo Sept. 2010
Zhi-Hong ZhangI; Shao-Yi WuI, II, *; Pei XuI; Li-Li LiI
IDepartment of Applied Physics, University of Electronic Science and Technology of China, Chengdu 610054, P.R. China
IIInternational Centre for Materials Physics, Chinese Academy of Sciences, Shenyang 110016, P.R. China
The electron paramagnetic resonance (EPR) parameters (g factors and the hyperfine structure constants) for Ni2+ and Co+ in MgO are theoretically studied from the perturbation formulas of these parameters for a 3d8 ion in octahedral crystal-fields. In the computations, the ligand orbital and spin-orbit coupling contributions are taken into account using the cluster approach. The calculated EPR parameters are in good agreement with the experimental data. The larger g factor and the smaller magnitude of the hyperfine structure constant for Ni2+ as compared with those for Co+ can be attributed to the higher spin-orbit coupling coefficient and the lower dipolar hyperfine structure parameter of the former, respectively.
Keywords: Impurity and defects; Electron paramagnetic resonance; Ni2+; Co+; MgO.
Magnesium oxide (MgO) is usually regarded as a model system to investigate magnetic [1,2], adsorption [3,4], structure [5,6] and optical properties  of doped transition-metal impurities. Particularly, MgO containing Ni2+ and Co+ can exhibit unique catalytic [8-10] and tunable laser properties [11,12]. Normally, these properties are closely correlated with the electronic states of the transition-metal ions in the hosts, which can be investigated by means of electron paramagnetic resonance (EPR) technique. For example, EPR experiments were carried out for Ni2+ and Co2+ doped MgO, and the EPR parameters (i.e., the isotropic g factors and the hyperfine structure constants) were also measured for the cubic Ni2+ and Co+ centers [13,14].
Up to now, however, the above EPR experimental results have not been satisfactorily explained. On the other hand, the EPR spectra and magnetic properties have been extensively investigated for Ni2+ in various chlorides by considering only the central ion orbital and spin-orbit coupling contributions . Nevertheless, the contributions to the EPR parameters from the ligand orbital and spin-orbit coupling interactions were not taken into account in the previous studies. In fact, for the 3d8 ions in oxides, the systems may still show some covalency and impurity-ligand orbital admixtures. In addition, the EPR spectra for MgO:Co+ have not been interpreted until now. Considering that (i) investigations on the EPR parameters for Ni2+ and Co+ in MgO can reveal useful information about electronic structures which would be helpful to understand the properties of these systems and that (ii) 3d8 ions can be treated as model systems containing only two unpaired holes, further quantitative studies on the EPR spectra for the Ni2+ and Co+ centers are of scientific and practical significance. In this work, the improved perturbation formulas of the EPR parameters based on the cluster approach are applied to the theoretical analysis of the Ni2+ and Co+ centers in MgO. In the calculations, the ligand orbital and spin-orbit coupling contributions are considered in a uniform way. The results are discussed.
Judging from the observed isotropic g factors and the hyperfine structure constants, the experimental EPR signals [13,14] can be assigned to the substitutional cubic Ni2+ and Co+ centers in MgO. When a 3d8 ion locates on an octahedral (Oh) site, the free-ion configuration 3F would be separated into two orbital doublets 3T1g and 3T2g and one singlet 3A2g, with the latter lying lowest corresponding to the isotropic g and A signals . As for the g factor of a 3d8 ion in octahedra, the perturbation formula was established using the conventional crystal-field model [17,18], by including only the contributions from the central ion orbital and spin-orbit coupling interactions. In order to study the EPR spectra of the 3d8 centers more exactly, the ligand orbital and spin-orbit contributions may be taken into account. Thus, the improved g formula based on the cluster approach is applied here. Meanwhile, the perturbation formula of the hyperfine structure constant for a 3d8 ion in regular octahedra can be similarly derived. Thus, we have  :
Here gs ( ≈ 2.0023) is the spin-only value. ζ and ζ´ are the spin-orbit coupling coefficients, k and k´ are the orbital reduction factors, and P and P´ are the dipolar hyperfine structure parameters for a 3d8 ion in crystals. κ is the core polarization constant. The energy denominators Ei(i = 1 ~ 3) denote the energy separations between the excited 3T2, 1T2a and 1T2b and the ground 3A2 states [17-19]. They can be described in terms of the cubic field parameter Dq and the Racah parameters B and C for the 3d8 ion in crystals: E1 ≈ 10 Dq, E2 ≈ 10 Dq + 12 B and E3 ≈ 8 B + 2C + 10 Dq [17-19]. From the cluster approach containing the ligand p- and s-orbital contributions , the spin-orbit coupling coefficients ζ and ζ´, the orbital reduction factors k and k' and the dipolar hyperfine structure parameters P and P' can be expressed as
Here and are the spin-orbit coupling coefficients of the free 3d8 and the ligand ions, respectively. A denotes the integral , where R is the impurity-ligand distance of the studied systems. Nγ and λγ (here the subscripts γ = t and e denote the irreducible representations T2g and Eg of Oh group, respectively) are the normalization factors and the orbital admixture coefficients. They are determined from the approximation relationships 
and the normalization conditions 
Here N is the average covalency factor, characteristic of the covalency of the studied systems. Sdpγ (and Sds) are the group overlap integrals. In general, the orbital admixture coefficients increase with increasing the group overlap integrals, and one can approximately adopt the proportional relationship λs/λe ≈ Sds/Sdpe between the orbital admixture coefficients and the related group overlap integrals within the same irreducible representation Eg. Obviously, omission of the ligand contributions (i.e., ζ´ = ζ = N, k´ = k = N, P´ = P = P0N), the above g formula returns to that of the previous work based on the conventional crystal-field model [17,18].
Usually, the impurity-ligand distance R is different from the host cation-anion distance RH in a pure crystal due to the difference between the ionic radius ri of impurity and the radius rh of host ions. Fortunately, studies based on experimental superhyperfine parameter and extended X-ray absorption fine structure (EXAFS) measurements have verified that the empirical formula R ≈ RH + (ri rh)/2 is approximately valid for impurity ions in crystals . From RH ( ≈ 2.105 Å  ) for MgO, ri ( ≈ 0.69 and 0.82 Å  ) for Ni2+ and Co+ as well as rh ( ≈ 0.66 Å  ) for Mg2+, the distances R are obtained and listed in Table 1. From the distances R and the Slater-type self-consistent field (SCF) wave functions [23,24], the group overlap integrals are calculated and shown in Table 1. According to the optical spectra for Ni2+ in MgO [25,26], the cubic field parameter Dq ≈ 860 cm-1 and the covalency factor N ≈ 0.90 can be obtained. Since the isoelectronic monovalent Co+ suffers weaker crystal-fields (i.e., lower Dq) and covalency effect (i.e., higher N)  when coordinated to the same oxygen ligands, the spectral parameters Dq ≈ 780 cm-1 and N ≈ 0.91 may be estimated for Co+ in MgO. Then the Racah parameters are determined from the relationships B ≈ N2B0 and C ≈ N2 C0 , with the corresponding free-ion values B0 ≈ 1208 and 878 cm-1 and C0 ≈ 4459 and 3828 cm-1  for Ni2+ and Co+, respectively. Using Eqs. (3) and (4), the molecular orbital coefficients Nγ and λγ (and λs) can be calculated. From the free-ion values ≈ 649 and 456 cm-1  for Ni2+ and Co+ and ≈ 151 cm-1  for O2-, the parameters in Eq. (2) are obtained and shown in Table 1. In the formula of the hyperfine structure constant, the dipolar hyperfine structure parameters are P0 ≈ 112 × 10-4 cm-1 and 228 × 10-4 cm-1  for Ni2+ and Co+, respectively. The core polarization constant can be determined from the empirical relationship κ ≈ -2χ/(3 < r-3 > ), where χ is characteristic of the density of unpaired spins at the nucleus of the central ion and < r-3 > the expectation value of the inverse cube of the radial wave function of the 3d8 orbital . Applying < r-3 > ≈ 7.094 and 5.388 a.u.  and χ ≈ -3.15 and -3.31 a.u.  for Ni2+ and Co+ in oxides, one can obtain κ ≈ 0.3 and 0.41 for MgO:Ni2+ and MgO:Co+, respectively. Substituting the above values into the formulas of the EPR parameters, the corresponding theoretical results (Cal. b) are calculated and shown in Table 2. To clarify the importance and the tendency of the covalency and the ligand contributions for Ni2+ and Co+ in MgO, the results (Cal. a) based on omission of the ligand contributions are also collected in Table 2 for comparison.
Table 2 indicates that the theoretical EPR parameters (Cal. b) for Ni2+ and Co+ in MgO based on the improved formulas containing the ligand contributions are in better agreement with the observed values than those (Cal. ª) in the absence of the above contributions.
1) The experimental isotropic g and A factors are attributed to the cubic Ni2+ and Co+ centers on the substitutional Mg2+ site in MgO. It is noted that there are some low symmetrical 3d8 centers in other oxides, e.g., the trigonal Ni2+ and Cu3+ centers in α Al2O3 . The anisotropic g factors g// and g⊥  for the trigonal centers may be similarly analyzed from the perturbation formulas  for a trigonally distorted octahedral 3d8 cluster, and the contributions from the trigonal distortion can be quantitatively considered in the calculations of the trigonal crystal-field parameters from the superposition model  and the local geometrical relationship of the impurity centers. Interestingly, the larger g factors  for Cu3+ than those for Ni2+ are attributable to the higher spin-orbit coupling coefficient ( ≈ 876 cm-1  ) of the former than that ( ≈ 649 cm-1  ) of the latter. In addition, Co+ replacing the host Mg2+ in MgO may induce some means of charge compensation due to the fewer charge of the impurity. On the other hand, larger size of Co+ than Mg2+ can lead to enhancement of the local metal-ligand interactions around the impurity and thus make Co+ stable on Mg2+ site. Further, local charge compensation (e.g., oxygen vacancy nearby) would break the original cubic symmetry of the ideal Mg2+ site and yield anisotropic EPR parameters. In view of the observed isotropic EPR signals [13,14], the charge compensation may occur in the outer ligand spheres far away from the impurity Co+, and the possible disturbance of the local structure of this center can be regarded as very small and negligible for simplicity. Of course, further experimental investigations of possible charge compensation for Co+ in MgO seem necessary and meaningful.
2) From Eqs. (1) and (2), the g factor largely depends upon the spin-orbit coupling coefficient ζ' (related to the central ion spin-orbit coupling coefficient). Thus, the larger g factor for MgO:Ni2+ than that for MgO:Co+ can be illustrated by the higher spin-orbit coupling coefficient ( ≈ 649 cm-1) of the former than that ( ≈ 456 cm-1) of the latter. On the other hand, the hyperfine structure constant is sensitively related to the dipolar hyperfine structure parameter P0 and the dominant contribution proportional to the core polarization constant κ. So, the larger magnitude of the hyperfine structure constant for MgO:Co+ than that for MgO:Ni2+ can be understood in view of the higher values of P0 and κ in the former.
3) The average covalency factors N ( ≈ 0.90 and 0.91 for Ni2+ and Co+) in this work still show some influences of the covalency on the EPR parameters, although the spin-orbit coupling coefficient ( ≈ 151 cm-1  ) of the oxygen ligand is much smaller than that ( ≈ 649 or 456 cm-1  ) of the impurity Ni2+ or Co+. Thus, omission of the ligand contributions yields larger g factors and slightly lower hyperfine structure constants in magnitude (Cal. ª). It seems that the improved formulas of the EPR parameters adopted in this work are superior to the previous ones [17,18] based on the conventional crystal-field model. Further, the covalency and the strength of the crystal-fields exhibit the decreasing tendency from Ni2+ to Co+ in the same MgO host, i.e., N (Ni2+) < N (Co+) and Dq (Ni2+) < Dq (Co+). This point is consistent with the lower valence state and hence weaker covalency and impurity-ligand interactions of the latter.
The EPR parameters for Ni2+ and Co+ in MgO are satisfactorily explained from the perturbation formulas based on the cluster approach. Inclusion of the ligand contributions yields better theoretical results as compared with those in the absence of these contributions. The larger g factor and the smaller magnitude of the hyperfine structure constant for Ni2+ than those for Co+ can be attributed to the higher spin-orbit coupling coefficient and the lower dipolar hyperfine structure parameter of the former.
 L.J. Shi, Phys. Lett. A 374 (2010) 1292. [ Links ]
 L.X. Ye, C.M. Lee, J.H. Lai, A. Canizo-Cabrera, W.J. Chen, T.H. Wu, J. Magn. Magn. Mater. 322 (2010) L9. [ Links ]
 S. Fernandez, A. Markovits, C. Minot, Chem. Phys. Lett. 463 (2008) 106. [ Links ]
 Y.F. Zhukovskii, E. A. Kotomin, G. Borstel, Vacuum 74 (2004) 235. [ Links ]
 Y. Wang, placeE. Florez, F. Mondragon, T. N. Truong, Surface Science 600 (2006) 1703. [ Links ]
 T. Mineva, V. Alexiev, C. Lacaze-Dufaure, E. Sicilia, C. Mijoule, N. Russo, J. Mol. Structure: THEOCHEM 903 (2009) 59. [ Links ]
 A. S. Shalabi, N. Abdullah, W. S. A. Halim, K. M. Eid, M. A. Kamel, H. O. Taha, Physica B 403 (2008) 316. [ Links ]
 D.L. Li, Y.Y. Zhan, K. Nishida, Y. Oumi, T. Sano, T. Shishido, K. Takehira, Appl. Catalysis A 363 (2009) 169. [ Links ]
 T. Furusawa, T. Sato, M. Saito, Y. Ishiyama, M. Sato, N. Itoh, N. Suzuki, Appl. Catalysis A 327 (2007) 300. [ Links ]
 X.H. Guo, Y. Li, R.J. Shi, Q.Y. Liu, E.S. Zhan, W.J. Shen, Appl. Catalysis A 371 (2009) 108. [ Links ]
 P.F. Moulton, Proc. IEEE 80 (1992) 348. [ Links ]
 D. Ma, X. Ma, J. Chen, Y. Liu, Phys. Rev. B 56 (1997) 1780. [ Links ]
 J.S. Thorp, M.D. Hossain, J. Magn. Magn. Mater. 22 (1981) 311. [ Links ]
 J.W. Orton, P. Auzins, J.E. Wertz, Phys. Rev. 119 (1960) 1691. [ Links ]
 M.G. Zhao, M.L. Du, G.Y. Sen, J. Phys. C 20 (1987) 5557. [ Links ]
 A. Abragam, B. Bleaney, Electron Paramagnetic Resonance of Transition Ions (Oxford University Press, London, 1970). [ Links ]
 A K. Petrosyan, A.A. Mirzakhanyan, Phys. Status Solidi B 133 (1986) 315. [ Links ]
 R.M. Macfarlane, Phys. Rev. B 1 (1970) 989. [ Links ]
 S.Y. Wu, W.Z. Yan, X.Y. Gao, Spectrochim. Acta A 60 (2004) 701. [ Links ]
 S.Y. Wu, X.Y. Gao, H.N. Dong, J. Magn. Magn. Mater. 301 (2006) 67. [ Links ]
 M. Moreno, M.T. Barriuso, J.A. Aramburu, Appl. Magn. Reson. 3 (1992) 283. [ Links ]
 R.C. Weast, CRC Handbook of Chemistry and Physics (CRC Press, Boca Raton, 1989) p. B196; F187. [ Links ]
 E. Clementi, D.L. Raimondi, J. Chem. Phys. 38 (1963) 2686. [ Links ]
 E. Clementi, D.L. Raimondi, W.P. Reinhardt, J. Chem. Phys. 47 (1967) 1300. [ Links ]
 R. Pappalardo, D.L. Wood, R.C. Jr. Linares, J. Chem. Phys. 35 (1961) 1460. [ Links ]
 S. Minomura, H.G. Drickamer, J. Chem. Phys. 35 (1961) 903. [ Links ]
 J.S. Griffith, The Theory of Transition-metal Ions (Cambridge University press, London,1964). [ Links ]
 M.G. Zhao, J.A. Xu, G.R. Bai, H.S. Xie, Phys. Rev. B 27 (1983) 1516. [ Links ]
 E.K. hodgson, placeI. fridovich, Biochem. Biophys. Res. Commun. 54 (1973) 270. [ Links ]
 B.R. McGarvey, J. Phys. Chem. 71 (1967) 51. [ Links ]
 P.R. Locher, S. Geschwind, Phys. Rev. Lett. 11 (1963) 333. [ Links ]
 Z.H. Zhang, S.Y. Wu, L.H. Wei, X.F. Wang, Y.X. Hu, Radiation Effects & Defects in Solids 164 (2009) 493. [ Links ]
 D.J. Newman, B. Ng, Rep. Prog. Phys. 52 (1989) 699. [ Links ]
(Received on 21 June, 2010)