SciELO - Scientific Electronic Library Online

 
vol.64 número1Numerical solutions of the Maung-Norbury-Kahana equation with the coulomb potential in momentum spaceSpace time geometry in the atomic hydrogenoid system. Approach to a dust relativistic model from causal quantum mechanics índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.64 no.1 México ene./feb. 2018

 

Research

Research in Condensed Matter

Theoretical studies on the electron paramagnetic resonance parameters for the tetragonal VO2+ center in CaO-Al2O3-SiO2 system

J.Z. Lina  * 

B.F. Zhanga 

1Department of physics, Suzhou Institute of Technology, Jiangsu University of Science and Technology, Zhangjiagang, 215600, P. R. China.


Abstract

In this work, the high order perturbation formulae of the electron paramagnetic resonance (EPR) parameters, g & A for a VO2+ 3d1 ion in tetragonal symmetry (compressed octahedron) are adopted to theoretically study these parameters of the tetragonal VO2+ center in CaO-Al2O3-SiO2 system. Because the spin-orbit (SO) coupling parameter ζp (151 cm-1) of ligand O2- is close to that (248 cm-1) of the central V4+ ion in the cluster [VO6]8- with strong covalence, the effect of the SO coupling parameter ζp on EPR parameters are considered on the basis of cluster approach. The obtained EPR parameters show good agreement with the experimental values and need fewer adjustable parameters than previous studies. And the local structure of the tetragonal VO2+ center in CaO-Al2O3-SiO2 system due to Jahn-Teller effect is obtained from the analysis of the EPR parameters. The results are discussed.

Keywords: VO2+; CaO-Al2O3-SiO2; electron paramagnetic resonance (EPR); crystal-fields; spin Hamiltonians

PACS: 71.70.Ch; 71.70.Ej; 75.10.Dg

1. Introduction

It is well-known that electron paramagnetic resonance (EPR) technique is a powerful tool to study the environmental symmetry of the transition-metal ions doped into the crystals (or glasses) which is produced by the ligands around the metal ion and can also give valuable information about the electric fields of these impurity ions 1. The vanadyl ion (VO2+) with 3d1 configuration is probably the most stable among the molecular paramagnetic transition metal ions and is often used as an impurity probe for detecting the site symmetry of the central ion and the bonding nature with the ligands in the EPR studies 2,3 Due to the strong V4+-O2- covalence bonding in VO2+, many VO2+ complexes in crystals (or glasses) possess C 4v symmetry with the axial symmetry for g and A values 4,5,6. For example, EPR investigations for VO2+ center in CaO-Al2O3-SiO2 glass-ceramic system were performed by Farah and the anisotropic g factors (g, g) and the hyperfine structure constants (A, A) were also obtained for this center 6. According to the EPR results, the oxygen octahedron around the central ion V4+ was found to suffer tetragonally compressed distortion 6. According to the EPR results, the oxygen octahedron around the central ion V4+ was found to suffer tetragonally compressed distortion 6. Based on the traditional crystal-field theory, theoretical calculations for EPR spectra of V4+ center in CaO-Al2O3-SiO2 were carried out through the complete diagonalization energy matrix method (CDM) and perturbation theory method (PTM) by Wei et al, 7 and the theoretical results show good agreement with the observed values (see Table I) 7. However, there may be some imperfections in their calculations. First, many adjusted parameters (i.e., the crystal field parameters B 20, B 40 and B 44, the orbital-reduction factor k and the core polarization constant κ) were introduced in their studies. Second, the contributions to EPR parameters from the ligand orbitals and SO coupling interactions were ignored for the strong covalent system CaO-Al2O3-SiO2: VO2+ characterized by the small orbital-reduction factor k(0.71) 7. Since the SO coupling coefficient ζp(151 cm-14) of the ligand O2- is close to that (248 cm-11) of the central V4+ ion, these contributions from ζp and the p- (s-) orbitals of the ligands may be somewhat important for this center due to strong covalency. Finally, No information about the local structures has been determined for the octahedral cluster [VO6]8- in CaO-Al2O3-SiO2 system, since previous computations failed to correlate with the defect structures around the impurity. Due to the reason that analysis of EPR spectra can provide useful information about local defect structures of the cluster [VO6]8- in CaO- Al2O3-SiO2 system, and this would be helpful in understanding the properties of this glass-ceramic system. Thus, further studies on the above EPR parameters are of scientific and practical meanings. So, in this paper, the EPR parameters are quantitatively analyzed based on high-order perturbation formulas of these parameters for a d1 ion in tetragonally compressed octahedral, the contributions from the tetragonally compressed distortion as well as the SO coupling and the orbitals of the ligands are considered in these formulas.

Table I. The EPR parameters of the tetragonal VO2+ center in CaO-Al2O3-SiO2 system 

aRef. 7.

bCalculations by using Eqs. (9)-(12) and including the tetragonal distortion in equation (20), but without considering the ligand orbital contributions (i.e., taking λγ=0, ζp=0 and then ζ'=ζ=Nζd, k = k' = N and P = P' = NP 0) in this work.

cCalculations by using Eqs. (9)-(12) and including both the tetragonal distortion in Eq. (20) and the ligand orbital contributions in this work.

dRef. 6.

2. Calculation

As for nd1 ions in tetragonally compressed octahedra, the higher orbital doublet 2eg in original cubic symmetry would split into two orbital singlets 2a1θ and 2b1ε. Meanwhile, the ground orbital triplet 2t2g in cubic symmetry can be divided into an orbital singlet 2 b 2 (ζ) and a doublet 2e(η,ξ), with the former lying lowest 8. Thus, the energy intervals can be written as follows:

E1=E(2e)-E(2b2)=-3Ds+5Dt (1)

E2=E(2b1)-E(2b2)=10Dq (2)

E3=E(2a1)-E(2b2)=10Dq-4Ds-5Dt (3)

Here, D q , D s and D t are the cubic field parameter and the tetragonal field parameters, respectively. They can be calculated from the local geometry of the studied impurity centers and the point-charge model 9,10,11, and be expressed as follows:

Dq=-eqr4/(6R¯5) (4)

Ds=2eqr2(1/R3-1/R3)/7 (5)

Dt=-2eqr4(1/R5-1/R5)/21 (6)

Here q(= -2e) is the effective charge of the oxygen ligand. r2 and r4 denote the expectation values of the radial wavefunction of d1 orbital in cluster, which can be expressed in terms of the corresponding free-ion values and average covalency factor N with respect to the admixture between the central metal ion and the ligand orbitals 9,10,11:

r23.3222N(a.u.),r421.0698N(a.u.) (7)

R¯ denotes the average metal-ligand (V-O) bonding length for the studied system (here, we take R¯1.985 Å for the VO2+ in cubic case 13). R and R denote the metal-ligand distances parallel and perpendicular to the tetragonal axis, respectively. For the VO2+ center in CaO-Al2O3-SiO2 system forming tetragonal [VO6]8- cluster, due to Jahn-Teller effect (note: the Jahn-Teller ion V4+ can suffer the Jahn-Teller effect under octahedral environments, which may slightly modify the immediate environment by stretching or contraction of the impurity-ligand bonding lengths), the parallel and perpendicular metal-ligand lengths can be expressed in terms of the average distance and the tetragonal compression ΔR (see Fig. 1) as:

RR¯-2ΔR,RR¯+ΔR (8)

Figure 1. The local structure of [VO6]8- cluster in CaO-Al2O3-SiO2 system. 

As mentioned before, previous studies 7 did not consider the contributions to EPR parameters from the s-& p-orbitals and the SO coupling coefficients of the ligands. Thus, using the routines provided in Refs. 14, the high-order (third-order) perturbation formulae of the g tensor and A tensor for 3d1 ion in tetragonally compressed octahedra including above contributions can be obtained based on cluster approach:

g=ge-8k'ζ'E2-2k'ζ2E12-4k'ζζ'E1E2-2ζ2E12 (9)

g=ge-2kζE1+kζ2E12-2kζ'2E1E2+2k'ζ'ζE1E2-4ζ'2E22-ζ2E12 (10)

A=P-κ-47+P'(g-ge)+37(g-ge) (11)

A=P-κ-27+P'1114(g-ge) (12)

Where ge2.0023 is the spin-only value. The energy denominators E1 and E2 stand for the energy divisions between the excited 2b1, 2e and the ground 2b2 states, as mentioned in Eqs. (1-3). P and P' are the dipolar hyperfine constant related to the interaction within t2g states and the interaction between t2g and eg states.

The SO coupling parameters ζ, ζ' the orbital reduction factors k, k' and the dipolar hyperfine constant P, P' in Eqs. (9-12) may be written as 15,16:

ζ=Nt(ζd+λt2ζp/2)       ζ'=(NtNe)1/2(ζd-λtλeζp/2) (13)

k=Nt1+λt2            k' =(NtNe)12[1-λtλe+λsA2] (14)

P=NtP0                  P'=(NtNe)1/2P0 (15)

ζd and ζp are the SO coupling coefficients of the free 3d1 and ligand ions, respectively. P0(128×10-4 cm-117) is the dipolar hyperfine structure parameter of the free VO2+ ion. A denotes the integral Rns/ynpy with the impurity-ligand distance R of the studied system. Nγ and λγ (or λs) are, respectively, the normalization factors and the orbital mixing coefficients for the cubic (O h ) irreducible representations γ(=eg and t2g) and can be determined from the normalization conditions 15,16

Nt(1-2λtSdpt+λt2)=1 (16)

Ne(1-2λeSdpe-2λsSds+λe2+λs2)=1 (17)

And the approximate relationships 15,16

N2=Nt2[1+λt2Sdpt2-2λtSdpt] (18)

N2=Ne2[1+λe2Sdpe2+λs2Sds2-2λeSdpe-2λsSds] (19)

Here Sdpγ (and S ds ) are the group overlap integrals. N is the average covalency factor, characteristic of the covalency effect (or reduction of the SO coupling coefficient and the dipolar hyperfine structure parameter) for the central ion in the studied system. generally, the mixing coefficients increase with increasing the group overlap integrals, and one can approximately adopt the proportional relationship between the mixing coefficients and the related group overlap integrals, i.e., λe/Sdpeλs/Ss within the same irreducible representation eg. Thus, the integrals Sdpt0.04553, Sdpe0.1198, Sds0.0896 and A1.0283 can be calculated from the Slater-type SCF functions 18,19 and the metal-ligand bonding length R (1.985 Å 13). According to the optical spectra parameter for VO2+ doped in glasses and crystals, 20,21 the average covalency factor N(0.81) can be adopted here. Then, the molecular orbital coefficients Nt0.8289, Ne0.8819, λt0.5021, λe0.4359 and λs0.3271 can be obtained from Eqs. (16-19). And the parameters ζ221 cm-1, ζ'198 cm-1k0.9334 and k'0.6893 can be determined from Eqs. (13-15) and the free-ion values ζd248 cm-1 for V4+5 and ζp151cm-1 for O2-4.

In Eqs. (11) and (12), the core polarization constant k is taken as 0.84, which is the same as κ(0.84) adopted in the previous work, 7. Thus, there are only one unknown parameters ΔR in the formulae of the EPR parameters. Substituting the related values (i.e., ζ,  ζ', k and k') into Eqs. (9)-(12) and fitting the calculated results to the experimental data, one can obtain:

ΔR0.1Å (20)

The theoretical results are shown in Table I. For comparison, the calculated results based on the formulae of the EPR parameters within the scheme of the conventional crystal-field theory and various fitted parameters in Refs. 7 are also given in the table.

3. Discussion

From Table I, we find that the theoretical results (Cal c ) based on the high-order perturbation formulae including the contributions from the ligand orbitals and SO coupling interactions show better agreement with the experimental values than those (Cal b ) neglecting the above contributions. This means that the contributions to the EPR parameters from the ligand orbitals and SO coupling interactions are important and should be considered in the explanations of these parameters for present system. In addition, these results exhibit a slight improvement compared with those (Cal a ) based on the formulae of the traditional crystal-field theory and various fitted parameters in previous work 7, although fewer adjustable parameters (the tetragonal distortion parameter ΔR, the core polarization constant κ ) than those (the crystal field parameters B 20, B 40 and B 44, the orbital reduction factor k, the core polarization constant κ) in Ref. 7 are adopted in this work.

  • (1) The studied system CaO-Al2O3-SiO2:VO2+ exhibits strong covalency and significant impurity-ligand orbital admixtures, characterized by the small covalency factor N(0.81<1) and the considerable mixing coefficients (Δt0.5021, λe0.4359 and λs0.3271). In view of the high valence state of V4+, the covalency and hence the significant admixture between the metal and ligand orbitals can be expected. In addition, the validity of the average covalency factor N can be further demonstrated by the relationship: N21-h(L)k(M) 22, where the parameter h(L) and k(M) are the characteristic of the ligand and central metal ions, respectively. From the values h(O 2-) and k(V2+)0.1 and k(V3+) 0.2 22 one can reasonably obtain k(V4+) 0.3 here by extrapolation. Thus, the value N0.84 is obtained which is very close to that adopted (N0.81) in this work.

  • (2) The positive sign of the tetragonal compression ΔR(0.1 Å) obtained in this work suggest that the oxygen octahedron around impurity V4+ ion exhibits tetragonally compressed distortion along C 4 axis, which is consistent with the ground state 2 b 2 based on the experimental EPR results g<g<ge and A>A). The local structure distortion can be explained by the Jahn-Teller effect. Because in compression, the ground state is nondegenerate singlet 2 b 2, while in elongation, the ground state is degenerate doublet 2 e, it can be expected that the latter is unstable. Similar distortions caused by Jahn-Teller effect are also found for V4+ (or other d1 ions, such as Nb4+23 and W5+24) in many other glasses and crystals 20,25,26. In addition, the crystal field parameters Dq(1685 cm-1), Ds(-3257 cm-1) and Dt(876 cm-1) obtained in present calculations are similar as V4+ in other crystals and glasses 2,20,21. So, the tetragonal compression ΔR obtained in this work is reliable.

  • (3) The calculated results show that the hyperfine structure constants A and A are negative as given in Table I, and the absolute values of them are in good agreement with the experimental values. It should be pointed out that the signs of the hyperfine structure constants are very difficult to determine. Thus, many experiments give them as absolute ones 1,2,6,21. Our calculations indicate that both signs of A and A are negative, this point is in accordance with the Ref. 7 and also be supported by the theoretical results for V4+ ions in many crystals 1,13,26 and glasses 20,25. In fact, from the formulae of A constants (Eqs. (11) and (12)), one can find that the negative signs of A constants are determined by the negative anisotropic g-shifts (g-ge, g-ge) and the negative terms related to core polarization constant κ which indicates the contribution to A constants by the unpaired s-electron. And the large value of κ (0.84) leads to large magnitudes of the hyperfine structure constants. It seems that the displacement toward central ion for the two ligands along C4 axis is larger than that (see Fig. 1) for the planar ligands may lead to more V4+ 3d-3s admixture and hence the larger κ. In addition, many previous works suggest that κ is within the range 0.6-1.0 for V4+ (or VO2+) ions in various crystals and glasses 1,2,13,20,25,26. The value of κ (0.84) adopted in this work lies within the range and can be regarded as rational.

  • (4) The calculation errors of this work can be discussed as follows. First, the approximation of the theoretical model (the cluster approach) and third-order perturbation formulae adopted here can lead to some errors. Second, the covalency factor N obtained from those of the similar [VO6]8- cluster may also introduce some errors into the calculated results. In fact, if one assumes an increase of the factor N by about 0.3%, the calculated results can be improved and show better agreement with the observed values (especially, the constants A and A) than Cal c in Table I. Third, the bond length of VO2+ in cubic field is applied in the calculations of the group overlap integrals, which is not the actual metal-ligand distance. This would introduce small errors into the integrals and hence bring forward an influence on the final EPR parameters. The errors of the final g factors are estimated to be no more than 1% as the integrals change by 10%. Finally, the contributions to the EPR parameters from charge-transfer (CT) mechanism are not involved. Fortunately, these contributions are very small for the cluster [VO6]8-27 and can be reasonably ignored.

4. Conclusions

The defect structure and EPR parameters for the tetragonal VO2+ center in CaO-Al2O3-SiO2 system are quantitatively studied from the high order (third order) perturbation formulae of these parameters including both the tetragonal distortion and the ligand orbitals contributions. The [VO6]8- cluster is found to suffer the tetragonal compression by about 0.1 Å along the C4 axis due to the Jahn-Teller effect. The negative signs of the constants A and A are suggested in this work.

Acknowledgments

This work was supported by the Foundation of the College Natural Science Research Project of Jiangsu Province (Grant No.15KJD140003), the PhD Startup Fund of Jiangsu University of Science and Technology, and the Youth Foundation of Zhangjiagang campus of Jiangsu University of Science and Technology.

References

1. E. Kalfaoglu, B. karabulut, J. Magn. Magn. Mater 324 (2012) 1593. [ Links ]

2. R. Kirpal, M.G. Misra, Appl. Magn. Reson. 44 (2013) 759. [ Links ]

3. R. Biyik, Physica B 404 (2009) 3483. [ Links ]

4. N.S. Reddy, R. Viswanatha, R.P.S. Chakradhar, J. Alloy. Compd 695 (2017) 1386. [ Links ]

5. R.R. Kumar, A.S. Rao, B.C.V. Reddy, Opt. Mater 4 (1995) 723. [ Links ]

6. H. Farah, J. Mater. Sci 38 (2003) 727. [ Links ]

7. Q. Wei, P.X. Zhang, D.Y. Zhang, J.H. Zhou, Pramana-J. Phys 73 (2009) 1087. [ Links ]

8. A. Abragam, B. Bleanely, Electron Paramagnetic Resonance of Transition Ions., (Oxford University Press, London, 1970). [ Links ]

9. M.G. Zhao, Chin. Phys. Lett 20 (2003) 2229. [ Links ]

10. M.G. Zhao, J. Chem. Phys 109 (1998) 8003. [ Links ]

11. M.G. Zhao, Y. Lei, J. Phys: Condens. Matter 9 (1997) 529. [ Links ]

12. F. Kubec and Z. Šroubek, J. Chem. Phys . 57 (1972) 1660. [ Links ]

13. J.T. Gong, L.J. Wang, W.L. Feng, Spectros. Lett. 43 (2010) 306. [ Links ]

14. D. Gourier, L. Colle, A.M. Lejus, D. Vivien, and R. Moncorge, J. Appl. Phys. 63 (1988) 1144. [ Links ]

15. H.M. Zhang, X. Wan, Z.H. Zhang, Physica B 407 (2012) 2201. [ Links ]

16. S.Y. Wu, Y.X. Hu, X.F. Wang, C.J. Fu, Radiat. Eff. Defect S. 165 (2010) 298. [ Links ]

17. R. Kripal, S. Shukla, Chin. Phys. Lett 26 (2009) 087601-1-4. [ Links ]

18. E. Clementi, D.L. Raimondi, J. Chem. Phys . 38 (1963) 2686. [ Links ]

19. E. Clementi, D.L. Raimondi, W.P. Reinhardt, J. Chem. Phys 47 (1967) 1300. [ Links ]

20. C.Y. Li, X.K. Huang, Q. Tu, W.Y. Wang, X.M. Zheng, Spectrosc. Spect. Anal 35 (2015) 1808. [ Links ]

21. B. Karabulut, I. Ilkin, R. Tapramaz, Z. Naturforsch 60a (2005) 95. [ Links ]

22. K.H. Karlsson, T. Perander, Chemica Scripa 3 (1973) 201. [ Links ]

23. X.X. Wu, X.P. Yu, W.C. Zheng, Acta. Phys. Pol. A 128 (2015) 72. [ Links ]

24. R.M. Peng, Y.M, W.C. Zheng, C.F. Wei, Physica B 461 (2015) 106. [ Links ]

25. W.L. Feng, Philos. Mag. 89 (2009) 1391. [ Links ]

26. Y.P. Huang, Z. Zhang, B.L. Yu, Optik 124 (2013) 4533. [ Links ]

27. L.X. Zhu, M.J .Wang, Dalton. Trans 42 (2013) 16289. [ Links ]

Received: June 19, 2017; Accepted: September 26, 2017

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License