All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

A Theoretical Calculation of Zero-Field Splitting Parameter for Mn2+ Ion in Cadmium Ammonium Phosphate Hexahydrate

Ram Kripal* and Manju Singh

EPR Laboratory, Department of Physics, University of Allahabad, Allahabad 211002, India

*Corresponding Author:
Ram Kripal
EPR Laboratory, Department of Physics, niversity of Allahabad, Allahabad 211002, India
+91 532 2470532
[email protected]

Received date: 22/07/2015 Accepte ddate: 01/09/2015 Published date: 03/09/2015

Visit for more related articles at Research & Reviews: Journal of Pure and Applied Physics


The crystal field parameters (CFP) and zero-field splitting (ZFS) parameter D of Mn2+ ion at axial symmetry site in cadmium ammonium phosphate hexahydrate (CAPH) single crystal are modeled. The calculation is carried out with the help of superposition model and the perturbation formulae. The theoretical ZFS parameter D is in good agreement with the experimental value. The theoretical modeling supports the conclusion drawn from experimental study that Mn2+ ion occupies substitutional site in CAPH. Ignoring local distortion, the value of D comes out as 306.1×10-4 cm-1; while considering small local distortion, the value of D is obtained as 308.4 ×10-4 cm-1. The experimental value of D is 305.5×10-4 cm-1.


Organic crystal, Crystal structure, Symmetry, Crystal fields, Spin-orbit effects, Electron paramagnetic resonance.


Electron paramagnetic resonance (EPR) studies yielding zero-field splitting (ZFS) parameters are important for the analysis of structural distortions and local site symmetry around the transition metal ions in crystals [1-7]. The microscopic spin-Hamiltonian (MSH) theory is widely used for the analysis of EPR spectra of transition metal ions. The spin Hamiltonian (SH) data from EPR spectra are correlated with optical spectroscopy and structural data with the help of MSH theory. Experimental and theoretical study indicated that the spin Hamiltonian parameters for transition metal (d5 (6S)) ions in crystals are affected very sensitively by the local structure. Consequently the SH theory has been extensively used in various types of crystal lattices.

The crystal-field parameters are evaluated by means of superposition model [8,9]. The ZFS parameters are determined through the electrostatic, the spin-orbit coupling, and the crystal-field parameters of the d5 ion in a crystal [10]. Out of iron group transition metal ions, Mn2+ is important for the reason that its ground state is 6S5/2 [11-15]. The electron spin is affected by high order interaction of crystalline electric field, and in the presence of external magnetic field orientation of spins become free [16]. The spin–lattice relaxation time is large for S state ions, due to which Mn2+ shows fine EPR spectra [2,15,17,18] at RT.

In this paper, the superposition model is used to obtain crystal field (CF) parameters and MSH theory together with CF parameters is used to predict the ZFS parameters for Mn2+ ions at the substitutional sites in cadmium ammonium phosphate hexahydrate (CAPH) crystal at room temperature (RT).

Crystal Structure

The crystal structure of CAPH is analogous to the naturally occurring biomineral struvite [19]. This crystal belongs to orthorhombic system having space group Pmn21, The unit cell dimensions are a = 6.941(2) Ǻ, b = 6.137(2) Ǻ and c = 11.199(4) Ǻ [20], and there exist two molecules per unit cell (Z = 2). Six water molecules form a distorted octahedral environment around the cadmium atom. The bond length of Cd-O lies within the range from 2.046 to 2.1076 Ǻ.

The Mn2+ ion in CAPH crystal substitutes at Cd2+ site [21]. The phenomenological SH of 3d5 ion configuration for tetragonal symmetry crystal field can be written as [2,22,23]




+image (1)

where the first term represents electronic Zeeman energy in which B is the external field, g is the spectroscopic splitting factor and μB is Bohr magneton. The second, third, and fourth terms are the second order axial, fourth-rank cubic and fourthrank axial ZFS terms [2]. The fifth term is the hyperfine interaction term. S, D, a, F represent the effective spin vector, second order axial, fourth-rank cubic and fourth-rank axial ZFS parameters, respectively. For the electronic Zeeman interaction isotropic approximation is used and usually it is valid for Mn2+ ions [2,24]. From the EPR spectra the direction of the maximum overall splitting is taken as the z axis and the direction of minimum splitting as the x axis [25]. The laboratory axis system (x, y, z) corresponds to the crystallographic axis system (CAS). The Z- axis of the local site symmetry axis i.e. the symmetry adopted axis (SAA) is along the metal-oxygen bond and the rest two axes X, Y are perpendicular to this.

The Hamiltonian for a d5 ion can be written as

H = H0+ Hcf + Hso (2)

in which Hcf=image

shows the crystal field Hamiltonian. H0 and Hso are free ion Hamiltonian and spin-orbit (SO) coupling, respectively. Because of spin-spin coupling being very small [26-28], its contribution has not been considered in Eq. (2). The total crystal field of SO interaction is taken as perturbation term [29,30]. The calculation in the strong-field scheme is formulated by Macfarlane for F-state ions [31]. In CAPH crystal, oxygen atoms form a distorted octahedron around the Cd atom [20] and yield axial local symmetry. For axial symmetry, the SO contributions to the ZFS parameter D for 3d5 ions can be written as [30]

image (3)

where P = 7(B+C), G = 10B+5C and D = 17B+5C. P, G, and D are the energy separations between the excited quartets and the ground sextet. Racah parameters B and C describe the electron-electron repulsion. In Eq. (3) only fourth order term is considered as other perturbation terms are very small [32]. By considering the average covalency parameter N, the parameter B, C and ξ are specified in terms of N as, B = N4B0, C = N4C0 and ξ = N2 ξ0, where B0 and C0 are the Racah parameters and ξ0 is the spin-orbit coupling parameter, for free ion [33,34]. B0= 960 cm-1, C0 = 3325 cm-1, ξ0 = 336 cm-1 [2] are used here.

From optical study of Mn2+ doped zinc ammonium phosphate (ZAPH) [35], B = 917 cm-1 and C = 2254 cm-1. As CAPH is isomorphous to ZAPH [35]; B, C and Dq are taken from this paper. By using these data in Eq. (4) [36] below,

image (4)

N can be obtained.

By applying the superposition model the crystal-field parameters for Mn2+ doped CAPH single crystal are calculated and using these parameters in Eq. (3) ZFS parameter D is determined. Similar approach has been used by several earlier workers [37] for estimating ZFS parameters.

Superposition Model

The superposition model is effectively used for explaining the crystal-field splitting. In recent times, this model is frequently applied for 3dn ions [8,32]. With the help of this model the crystal field parameters can be determined by the equations [38]

image 5(i)

image 5(ii)

image 5(iii)

where R0 is the reference distance. In general, it is chosen as the mean value of all six bond lengths (for axial symmetry). The mean value of two out of six Cd-O bond lengths are indicated by R10 and the mean value of remaining four bond lengths are indicated by R20, further ΔR1 and ΔR2 indicate the distortion parameters. Āk and tk are the intrinsic parameter and power law exponent, respectively [39,40].

For 6-fold cubic coordination image [10]. For 3d5 ions, image lies in the range 8-12 [32,40]. Using Dq= 756 cm-1 [35], imagecan be determined taking the value of ratio image. The power law exponent in the case of Mn2+ ion is adopted as t2 = 3, t4 = 7 [41]. For the intrinsic parameter values for the superposition model semi-ab initio calculations are done for other transition metal ions [42].

Results and Discussion

In the present study the bond length of six ligands Rj (Figure 1) is calculated as given in Table 1. The value of B and C are taken from the optical study [35] as 917 and 2254 cm-1, respectively. First of all ignoring the local distortion the value of D is calculated. For this, taking image = 8 and reference distance R0 = 0.28 nm, which is the mean value of six ligands distance, the value of D comes out to be D = 87342×10-4 cm-1. From EPR study, the experimental value is given as D = 305.5×10- 4 cm-1 [20]. Here it is seen that the theoretical value is very high. Now, if the ratio image is changed to 12, then D = 79282×10-4 cm-1, which is again high. Further, using the same ratio and R0 = 0.212 nm, which is the sum of ionic radii of Mn2+ = 0.08 nm and O2- = 0.132 nm, the result comes out as D = 1103 ×10-4 cm-1, which is also very large. If image = 8 is taken, then D becomes 467.4×10-4 cm-1, which is again large. Now by taking R0 = 0.210 nm, the value comes out as D = 306.1×10-4 cm-1, which is in reasonable agreement with the experimental value D = 305.5×10-4 cm-1.


Figure 1: Coordination about Mn2+ ion in CAPH.

Metal-Oxygen Metal-Oxygen bond distance Rj(nm)
Mn- O(2) 0.4414
Mn-O(3) 0.4117
Mn-O(H1) 0.1887
Mn-O(H2) 0.1992
Mn-O(H3) 0.2185
Mn-O(H4) 0.2234

Table 1: Metal-oxygen bond distances Rj for Mn2+ ion in cadmium ammonium phosphate hexahydrate single crystal.

Further, considering local distortions as ΔR1 = -0.0008 nm and ΔR2 = -0.0009 nm, reference distance R0 = 0.28 nm (the mean value of six ligands distance) and ratio image = 8, the value of D comes out to be 62759.9×10-4 cm-1. It is noticed that the theoretical value is very high. Again taking the same value of R0 (0.28 nm) but changed value of the ratio, image = 12, the value of D comes out as 55681.7 ×10-4 cm-1, which is also larger than the experimental value. Further, using the same ratio but R0 = 0.212 nm, which is the sum of ionic radii of Mn2+ = 0.08 nm and O2- = 0.132 nm, the result comes out as D = 1289.8 ×10-4, which is again very high and if image = 8 is taken with the same R0 = 0.212 nm, the result comes out as D = 37.0×10-4 cm-1, which is very low as compared to experimental value. At R0 = 0.166 nm, the value of D comes out as 308.4 ×10-4 cm-1, which is in reasonable agreement with the experimental result D = 305.5×10-4 cm-1 .

By using Eq. 5(i) - 5(iii), crystal-field parameters Bkq are obtained as: B20 = -19896 cm-1, B40 = 7468.1 cm-1, B44 = 10325.1 cm-1 when distortion parameters ΔR1 and ΔR2 are not considered in the calculation. Considering ΔR1 = -0.008 Ǻ and ΔR2 =-0.009 Ǻ, the Bkq parameters are: B20 = -9484 cm-1, B40 = 1329.8 cm-1, B44 = 1839.6 cm-1. The above theoretical investigation suggests that Mn2+ ion enters the host lattice substitutionally in CAPH. This supports the conclusion drawn by EPR study.


Zero-field splitting parameter D for Mn2+ doped CAPH in axial symmetry has been calculated with the help of superposition model and perturbation formulae. The theoretical value of D is in reasonable agreement with the experimental value. This analysis indicates that Mn2+ ion enters the host lattice substitutionally, which supports the experimental result obtained by EPR.


The authors are thankful to the Head, Department of Physics, University of Allahabad, Allahabad for providing departmental facilities.