## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733

### Braz. J. Phys. vol.40 no.3 São Paulo Sept. 2010

#### http://dx.doi.org/10.1590/S0103-97332010000300016

**Theoretical studies of the local structure and the EPR parameters for substitutional Mo ^{5+} in TiO_{2} **

**Ji-Zi Lin ^{*}**

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

**ABSTRACT**

The electron paramagnetic resonance(EPR) parameters g-factors g_{i}(i = x,y,z) and the hyperfine structure constants A_{i} for the substitutional Mo^{5+} in rutile are theoretically studied from the perturbation formulas of these parameters for a d^{1} ion in a rhombically compressed octahedron. On the basis of the studies, The oxygen octahedron around Mo^{5+} is found to transform from the original elongation on the host Ti^{4+} site to compression in the impurity center due to the Jahn-Teller effect. The calculated EPR parameters based on the above local structure in this work are in good agreement with the experimental data.

**Keywords: **Electron paramagnetic resonance (EPR); Crystal-fields and spin hamiltonians; Mo^{5+}; Rutile(TiO_{2}).

**1. INTRODUCTION**

Rutile (TiO_{2}) doped with transition-metal ions has attracted interest extensively because of its unique phase transition [1] and magnetic [2-4], optical [5,6] and catalytic [7-9] properties. In general, these properties are sensitively dependent upon the local structures around the impurity ions in the host, It is well known that electron paramagnetic resonance (EPR) is a useful technique to study the electronic and local structure properties of paramagnetic impurities in crystals. For example, EPR investigations were performed decades ago, and one rhombic center was observed for the substitutional Mo^{5+} in rutile by means of EPR technique, and the EPR parameters anisotropic g factors g* _{i}*(i = x,y,z) and the hyperfine structure constants A

_{i}were also measured at 77ºK [10].

Till date, however, the above experimental results have not been satisfactorily interpreted. In the treatments of the g factors based on angular overlap model and the simple perturbation formulas in Refs [11]. , the contributions from the ligand orbitals and spin-orbit coupling interactions were ignored. In fact, the covalency effect and hence the admixtures between the orbitals of the metal and the ligand ions may become significant, due to the high valence state of the impurity Mo^{5+}. Meanwhile, the energy denominators in the formulas of g factors were not correlated with the local structure around the impurity center. Actually, the local structure properties around Jahn-Teller ions(such as Ni^{3+}, Ti^{3+},V^{4+}) in crystals may be dissimilar to those of the hosts due to the local modification of bond-lengths and angles arising from the Jahn-Teller effect [12-14], In view that the studies on the Jahn-Teller behaviours of Mo^{5+} in rutile would be useful to understand the properties of this material with dopants. In addition, the hyperfine structure constants were not interpreted.

In order to make further investigations on the EPR parameters for the substitutional Mo^{5+} center in rutile, the EPR parameters are quantitatively analyzed from the second-order perturbation formulas of these parameters for a d^{1} ion in rhombic symmetry, by considering reasonable local lattice (angular) distortion around the impurity ion. In these formulas, the contributions from the rhombic distortion as well as the spin-orbit coupling and the orbitals of the ligands are taken into account based on the cluster approach.

**2. CALCULATIONS**

In rutile, there are one substitutional and one interstitial octahedral site which can be occupied by the doping transition-metal ions. The oxygen octahedra in both sites exhit significant low symmetrical distortions. For the substitutional site, the ligand octahedron is slightly elongated with two longer bond-lengths *R*_{||}( ≈ 1.99 *Å*) parallel with the [110] axis and four coplanar shorter bond-lengths* R*⊥( ≈ 1.94 *Å*) perpendicular to the axis, with the planar bonding angle q of about 81º [15] .The above low symmetrical distortions can be characterized by the axial distortion angle Δα = α – α_{0}, with α ≈ tan^{–1}(*R*_{⊥}/*R*_{//}) and α_{0} ≈ 45º of the cubic symmetry. Thus, we have Δα ≈ –0.73º for the host substitutional site [15]. Additionally, the ligand octahedron for this site also suffers large perpendicular (or rhombic) distortion characterized by the deviation Δθ of the plannar bonding angle θ related to the cubic value θ_{0} ( = 90º), i.e., Δθ = θ – θ_{0} = –9º [15]. Thus, the ligand octahedron has large rhombic distortion for the host substitutional site . When Mo^{5+} is doped into the lattice of TiO_{2}, it may occupy the substitutional site [10]. and the local structure around the ion Mo^{5+} may exhibit deformation because of the Jahn-Teller effect, the impurity-ligand bonding lengths and angles may be unlike those in the host and lead to local lattice (angular) distortion.

For a d^{1}(Mo^{5+}) ion in rhombically-distorted octahedra, its higher orbital doublet ^{2}E* _{g}* of the original cubic case would split into two orbital singlets

^{2}A

_{1g}( |z

^{2}> ) and

^{2}A'

_{1g}( |xy > ), while the original lower orbital triplet

^{2}T

_{2g}would be separated into three orbital singlets

^{2}B

_{2g}( |xz > ),

^{2}B

_{3g}|yz > ) and

^{2}B

_{1g}( |x

^{2}-y

^{2}> ), with the latter lying lowest [11]. Note that in the above notations xy > and |x

^{2}-y

^{2}> are interchanged due to a rotation of the frame of axes in the XY plane.

According to Refs. [11], the perturbation formulas of the g factors for a d^{1} ion in rhombic symmetry were established based on the conventional crystal-field theory. The experimental g factors were investigated from the above formulas by adjusting various parameters (i.e.,the energy separations between the ground and the excited states). Nevertheless, the contributions from the ligand orbitals and spin-orbit coupling interactions were ignored in these treatments. By including the above contributions and the similar perturbation procedure in Refs. [16], the perturbation formulas of the g factors and the hyperfine structure constants for the d^{1} ion in rhombically-distorted octahedra can be obtained on the basis of the cluster approach:

where g_{s} ( = 2.0023) is the spin-only value. *P*_{0} is the dipolar hyperfine structure parameter of the free d^{1} ion. κ is the isotropic core polarization constant. *N* is the average covalency factor, characteristic of the covalency effect of the studied system. The energy denominators E_{1}, E_{2} and E_{3} stand for the energy separations between the excited ^{2}B_{2g}, ^{2}B_{3g} and ^{2}A'_{1g} and the ground ^{2}B_{1g} states and can be obtained from the strong cubic field approach:

Here *D _{q}* is the cubic field parameter, and

*D*

_{s}, D_{t}, D_{ξ}and

*D*

_{η}the rhombic ones.

Based on the cluster approach, the spin-orbit coupling coefficients ζ, ζ' and the orbital reduction factors *k*, *k*' in equation (1) may be written as

where ζ_{d} and ζ_{p} are the spin-orbit coupling coefficients of the d^{1} and the ligand ions in free states, respectively. *A* denotes the integral , where *R* is the impurity-ligand distance in the studied system. *N*_{γ} and λ_{γ} (or λ_{s}) are, respectively, the normalization factors and the orbital mixing coefficients for the cubic (*O _{h}*) irreducible representations γ( =

*e*and

_{g}*t*

_{2g}). They are usually determined from the normalization conditions

and the approximate relationships

Here *N* is the average covalency factor. *S _{dp}*

_{γ}(and

*S*) are the group overlap integrals. In general, 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., λ

_{ds}*≈ λ*

_{e}/S_{dpe}_{s}

*/S*within the same irreducible representation

_{s}*e*.

_{g}From the superposition model [17] and the geometrical relationship of the studied impurity center, and rhombic field parameters can be determined as follows:

Here θ'(=θ_{0} + Δθ') denotes the planar bonding angle for the substitutional site in rutile, . _{2} (*R*_{0}) and _{4}(*R*_{0}) are the intrinsic parameters (with the reference bonding length *R*_{0}), while *t*_{2} ( ≈ 3) and *t*_{4}( ≈ 5) are the power-law exponents [18]. The relationships _{4}(*R*_{0}) = *D _{q}* and

_{2}(

*R*

_{0}) = 9 ~ 12

_{4}(

*R*

_{0}) can be obtained from Refs. [18] and we take

_{2}(

*R*

_{0}) = 12

_{4}(

*R*

_{0}) here. The reference bonding length

*R*

_{0}can be reasonably estimated from the empirical relationship

*R*

_{0}≈ + (

*r*)/2, = (

_{i }– r_{h}*R*

_{//}+ 2

*R*

_{⊥})/3 ≈ 1.957

*Å*(the average metal-ligand bonding length),

*r*( ≈ 0.66

_{i}*Å*)and

*r*( ≈ 0.68

_{h}*Å*) denote the ionic radius for the impurity Mo

^{5+}and the host Ti

^{4+}[19]. From the distance

*R*

_{0}( ≈ 1.947

*Å*) and the Slater-type SCF functions [20,21] the integrals

*S*≈ 0.047,

_{dpt}*S*≈ 0.1168,

_{dpe}*S*≈ 0.0935 and

_{ds}*A*≈ 1.2673 are obtained. As mentioned before, because of the Jahn -Teller effect, the impurity-ligand bonding lengths

*R'*

_{//}and

*R'*

_{⊥}, and hence the angle α' ≈ tan

^{-1}(

*R'*

_{⊥}/

*R'*

_{//}) may be different from those (

*R*

_{//},

*R*

_{⊥}and α) in the host. The impurity-ligand bonding lengths can be written in terms of the local distortion angle >Δα'( = α' – α

_{0}), the average distance as:

For TiO_{2}:Mo^{5+}, to our knowledge, no optical spectra data were reported. However, the value *D _{q}* ≈ 2500cm

^{–1}was obtained for the [MoO

_{4}]

^{3–}cluster in CaWO

_{4}:Mo

^{5+}[22] .thus, the value

*D*≈ 3750cm

_{q}^{–1}can be approximately estimated for the studied octahedral [MoO

_{6}]

^{7–}cluster here.

Thus, by using the formulas of the g factors (Eq. (1)), and fitting the calculated g factors to the experimental data, we have

The values *N _{t}* ≈ 0.7550,

*N*≈ 0.8124, λ

_{e}_{t}≈ 0.6187, λ

_{e}≈ 0.5098 and λ

_{s}≈ 0.4081 are calculated from equations (4) and (5). Then the parameters ζ ≈ 799cm

^{–1}, ζ' ≈ 788cm

^{–1},

*k*≈ 0.899 and

*k'*≈ 0.534 can be determined from equation (3) and the free-ion values ζ

_{d}≈ 1030cm

^{-1}for Mo

^{5+}[22] and ζ ≈ 151cm

^{–1}for O

^{2–}[23].

In the formulas of the hyperfine structure constants, the free-ion value *P*_{0} is about –68.2 × 10^{–4} cm^{–1} for Mo^{5+} [24]. By fitting the calculated hyperfine structure constants to the observed values, the core polarization constant is determined: κ ≈ 0.42.

The corresponding EPR parameters are shown in table 1. For comparisons, the calculated results based on the formulas of the g factors within the scheme of the conventional crystal-field theory and various fitted parameters in Refs. [11] are also collected in Table 1.

**3. DISCUSSION**

From table 1, one can find that the calculated EPR parameters based on the perturbation formulas in Eq. (1) by considering the local lattice deformation as well as the ligand orbital and spin-orbit coupling contributions show better agreement with the experimental data than those in the absence of the above contributions. This means that the formulas and the related parameters adopted in this work can be regarded as suitable.

The local axial and rhombic distortion angles Δα' and Δθ' for the impurity center are different from those in the host sites. The ligand octahedron in the substitutional center is found to be slightly (Δα' ≈ 0.995º) compressed while elongated (Δα ≈ –0.73º) in the host Ti^{4+} site, and the perpendicular distortion (Δθ' ≈ –0.305º) are much smaller than that (Δθ ≈ –9.0º) in the replaced cation site. The above distortions can be described as the contraction (by about 0.08*Å*) and relaxation (by about 0.04*Å*) of the parallel and perpendicular bond lengths, and the angle θ' is much closer to that ( ≈ 90º) of cubic case. As a result, the ligand octahedron is slightly compressed. The reason may be due to the Jahn-Teller effect which may lend to compress the ligand octahedron, as Ni^{3+} on the octahedral Al^{3+} site in LaSrAlO_{4} [12]. Interestingly, similar local lattice modifications around impurity ions in crystals are also found for substitutional V^{4+} centres in rutile-type crystals [13]. Thus, the local structure properties around Mo^{5+} ion in rutile obtained in this work can be regarded as reasonable and the anisotropy of the experimental g factors can be understood, by considering the local axial and rhombic distortion angles Δα' and Δθ'.

Seen from Table 1 one can find that the agreement between the theoretical g factors and the experimental data in present work is better than that in Refs. [11]. Based on the formulas of the g factors within the scheme of the conventional crystal-field theory in Refs. [11], the observed g factors were explained by various fitted parameters and the energy separations were not determined from the local structure of the impurity center. In addition, the contributions from the SO coupling coefficient and the orbitals of the ligands were not taken into account. In view of the high valence state of Mo^{5+}, the distinct covalency and hence, significant admixture of the metal and ligand orbitals can be expected. This point may be illustrated by the small covalency factor *N* ( ≈ 0.733 < 1) in equation (5) and the obvious mixing coefficients (λ_{t} ≈ 0.6187, λ_{e} ≈ 0.5098 and λ_{s} ≈ 0.4081) obtained in this work. Based on the studies, inclusion of the contributions from the SO coupling coefficient and the orbitals of the ligands leads to the variations of about 0.04, 0.05 and 0.016 for the calculated *g _{x}, g_{y}* and

*g*, respectively. Obviously, for the impurity ion Mo

_{z}^{5+}with high valence state in the octahedral cluster, the above contributions could not to be ignored.

There may be some errors in the above calculations. First, the approximation of the perturbation formulas and the cluster approach adopted here can induce some errors. Second, the cubic field parameter Dq extrapolated from the [MoO_{4}] ^{3–} cluster in CaWO_{4} would lead to some errors. Finally, the impurity-ligand distance is obtained from the empeirical relationship related to the difference between the ionic radius for imputity and host ions, which may bring forward some errors in the average distance for present system. It means that the actual average distance can be slightly different from . This would bring small errors to the group overlap integrals in the calculations. In fact, the errors in the final g factors and distortion angles in Eqs. (1) and (7) due to the above variation in the average distance turn out to be no more than 0.4%, since is mainly related to the magnitude of the g factors.

**4. CONCLUSIONS**

The electron paramagnetic resonance (EPR) parameters g-factors g_{i}(i=x,y,z) and the hyperfine structure constants A* _{i}* for the substitutional Mo

^{5+}in rutile are theoretically investigated in this work. Based on the studies, The oxygen octahedron around Mo

^{5+}is found to transform from the original elongation on the host Ti

^{4+}site to compression in the impurity center due to the Jahn-Teller effect. The theoretical EPR parameters based on the above local structure in this work are in good agreement with the experimental data.

[1] H.L. Ma, J.Y. Yang, Y. Dai, Y.B. Zhang, B. Lu, G. H. Ma. Appl. Surf. Sci., 253, 7497 (2007). [ Links ]

[2] G.Y. Gao, K.L. Yao, Z.L. Liu. Phys. Lett. A, 359, 523 (2006). [ Links ]

[3] N. Akdogan, B.Z. Rameev, R.I.K. haibullin, A. Westphalen, L.R. Tagirov, B. Aktas, H. Zabel. J. Magn. Mater., 300, e4 (2006). [ Links ]

[4] G.Y. Gao, K.L. Yao, Z.L. Liu, Y.L. Li, J.L. Jiang, Y.C. Li. Physica B, 382,14 (2006). [ Links ]

[5] T.N. Ishikawa, T. Nakajima ,Y. Shinohara. J. Mol. Struct-Theochem., l545, 67 (2001). [ Links ]

[6] F. Gravier , B. Boulanger. Opt. Mater. 30, 33 (2007). [ Links ]

[7] J. Wang, Y.H. Lv, L.Q. Zhang, B. Liu, R.Z. Jiang, G.X. Han, R. Xu, X.D. Zhang. Ultrason Sonochem, 17, 642 (2010). [ Links ]

[8] X.W. Zhang, L.C. Lei. J. Haza. Mater.,153, 827 (2008). [ Links ]

[9] Y. Dongsun, K. Iigon, K. Sangsoo, H. H. Chang, L. Changyu, C. Seongjin. Appl. Surf. Sci., 253, 3888 (2007). [ Links ]

[10] R. T. Kyi. Phys. Rev., 128, 151 (1962). [ Links ]

[11] P. DE Montgolfier, Y. Boudeville. Chem. Phys. Lett., 37 (1), 97 (1976). [ Links ]

[12] Y.V. Yablokov, T.A. Ivanova. Coord .Chem. Rev., 190-192, 1255 (1999). [ Links ]

[13] S.Y. Wu, L.H. Wei, Z.H. Zhang, X.F. Wang, J.Z. Lin. Z. Naturforsch. 63a, 523 (2008). [ Links ]

[14] S. Lenjer, O. F. Schirmer, H. Hesse, Th. W. Kool. Phys. Rev B, 66, 165106 (2002). [ Links ]

[15] H. Hikita, K. Takeda, Y. Kimura. Phys. Rev. B, l46 (22), 14381 (1992). [ Links ]

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

[17] D.J. Newman, B. Ng. Rep. Prog. Phys., 52, 699 (1989). [ Links ]

[18] W. L. Yu, X. M. Zhang, L. X. Yang, B. Q. Zen. Phys. Rev. B, 50, 6756 (1994). [ Links ]

[19] K. Chandrasekharan, V.S. Murty, F.J. Kumar, P. Ramaswamy, P.S. Rao. J. Mater. Sci., 35, 1561 (2000). [ Links ]

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

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

[22] A.S. Chakravarty. Introuduction to the Magnetic Properties of Solids. New York: Wiley -Interscience Publication.1980. [ Links ]

[23] E.K. Hodgson, I. Fridovich. Biochem. Biophys. Res. Commun., 54, 270 (1973). [ Links ]

[24] B.R. McGarvey. J. Phys. Chem., 71 (1), 51 (1967). [ Links ]

(Received on 17 April, 2010)

* Electronic address: linjizi13@126.com