SciELO - Scientific Electronic Library Online

 
vol.15 issue6Determination of the numerical parameters of a continuous damage model for the structural analysis of clay brick masonryCement mixtures containing copper tailings as an additive: durability properties author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Materials Research

Print version ISSN 1516-1439

Mat. Res. vol.15 no.6 São Carlos Nov./Dec. 2012 Epub Oct 11, 2012

http://dx.doi.org/10.1590/S1516-14392012005000135 

Microstructural changes and effect of variation of lattice strain on positron annihilation lifetime parameters of zinc ferrite nanocomposites prepared by high enegy ball-milling

 

 

Abhijit BanerjeeI; Srinjoy BidII; Hema DuttaIII; Sandeep ChaudhuriIV; Dipankar DasV; Swapan Kumar PradhanI, *

IDepartment of Physics, The University of Burdwan, Golapbag, Burdwan-713104, West Bengal, India
IIDepartment of Physics, Academy of Science, Aidco Nagar, Hooghly, West Bengal, India
IIIDepartment of Physics, Vivekananda College, Burdwan-713103, West Bengal, India
IVDepartment of Physics, University of Surrey, Guildford, Surrey, UK
VUGC-DAE Consortium for Scientific Research, Kolkata 700 098

 

 


ABSTRACT

Zn-ferrite nanoparticles were synthesized at room temperature by mechanical alloying the stoichiometric (1:1 mol%) mixture of ZnO and α-Fe2O3 powder under open air. Formation of both normal and inverse spinel ferrite phases was noticed after 30 minutes and 2.5 hours ball milling respectively and the content of inverse spinel phase increased with increasing milling time. The phase transformation kinetics towards formation of ferrite phases and microstructure characterization of ball milled ZnFe2O4 phases was primarily investigated by X-ray powder diffraction pattern analysis. The relative phase abundances of different phases, crystallite size, r.m.s. strain, lattice parameter change etc. were estimated from the Rietveld powder structure refinement analysis of XRD data. Positron annihilation lifetime spectra of all ball milled samples were deconvoluted with three lifetime parameters and their variation with milling time duration was explained with microstructural changes and formation of different phases with increase of milling time duration.

Keywords: nano-Zn-ferrites, ball-milling, Rietveld method, positron annihilation


 

 

Introduction

Ferrites are a group of technologically important materials used in magnetic, electronic and microwave fields. Magnetic nanocrystalline materials hold great promise for atomic engineering of materials with functional magnetic properties1-3. Many magnetic nanocrystals show superparamagnetism in single domain particles below a certain critical size. Magnetic nanocrystals have been extensively applied in magnetic recording medium, information storage, bio-processing and magneto-optical devices4,5. The sulfur absorption capacity of milled zinc ferrite increases with decreasing crystallite size due to structure-reactivity relationship at high temperature.

Crystalline ZnFe2O4 (cubic, a = 0.8441 nm, space group: Fdm, Z = 8; ICDD PDF #22-1012) is a normal spinel at room temperature. Spinel structure is made up by a cubic close-packed array of oxygen atoms with tetrahedral (T) and octahedral (O) cavities. In the normal 2-3 spinels, one-eighth of the T sites and one-half of the O sites are filled by the divalent (A) cations (Mg2+, Zn2+, Mn2+, Cd2+, etc.) and the trivalent (B) cations (Al3+, Fe3+, Cr3+ etc.), respectively, in the ratio AB2O4. The structural formula of zinc ferrite is usually written as (Zn1-δ2+ Feδ3+)[Znδ2+ Fe2-δ3+]O42-, where round ( ) and square [ ] brackets denote T and O sites of co-ordination, respectively, and δ represents the degree of inversion (defined by the fraction of the T sites occupied by B cations). There are two ordered configuration stable at low temperature, the one with δ = 0 (normal spinel) and the other with δ = 1 (inverse spinel). When the temperature increases, disorder takes place, since A and B cations undergo increasing intersite exchange over the three cation sites per formula unit (one T and two O sites). The completely random distribution of A and B cations over the three cation sites corresponds to δ = 2/3, which is asymptotically approached at very high temperatures. Same type of cation distribution is also observed in ball-milled samples6-7. The change in temperature or the change in milling parameter may result in change in the degree of inversion. It has been found that a metastable nanoscale structural state of mechanosynthesized ZnFe2O4 is characterized by a substantial displacement of Fe3+ cations to tetrahedral sites and of Zn2+ cations into octahedral sites. The inverse-normal transition in mechanosynthesized zinc ferrite proceeds rapidly in the temperature range 885-1073 K and the activation energy of the transition is E ~72 kJ.mol-1[6].

Formation of nanocrystalline ZnFe2O4 as normal and inverse spinel structures was noticed after ball-milling the stoichiometric mixture of ZnO and α-Fe2O3 powders in open air for different lengths of time. Formation of nanocrystalline materials in the process of ball-milling leads to significant amount of structural and microstructural defects which can be characterized by X-ray diffraction and positron annihilation spectroscopy studies.

The powder patterns of almost all the ball-milled materials, milled at different milling time are composed of a large number of overlapping reflections of α-Fe2O3, ZnO and ZnFe2O4 (normal and inverse spinel) phases. The Rietveld analysis based on structure and microstructure refinement8-9 was adopted for precise determination of several microstructural parameters as well as relative phase abundance of such multiphase material.

The purpose of the present work is to characterize defect states of the zinc ferrite nanocomposites by positron annihilation lifetime spectroscopy (PALS) and establish their relationship with microstructure parameters obtained from X-ray analysis. PALS is a powerful technique to characterize defects in solid materials10. A nanostructured material contains unlimited grain boundaries and these boundaries are rich in lattice defects. Thus, PALS method can be effectively used to characterize nanostructures in terms of lattice imperfections. To the best of our knowledge, this type of analysis will help to understand the nature of deformation generated in the process of mechanical alloying in nanocrystalline ferrite powders.

 

Experimental Method

High-energy ball-milling of ZnO (M/s Merck, 99% purity) and α-Fe2O3 (M/s Glaxo, 99% purity) mixture in 1:1 mol% was conducted in a planetary ball mill (Model P5, M/S Fritsch, GmbH, Germany). The time of milling was varied from 30 minutes to 10 hours depending upon the rate of formation of zinc ferrite phase. The step-scan data (of step size 0.02° 2θ and counting time 5 seconds) for the entire angular range (15-80° 2θ) of the unmilled mixture and all ball-milled samples were recorded using Ni-filtered CuKα radiation from a highly stabilized and automated Philips X-ray generator (PW 1830) operated at 35 kV and 25 mA.

For PALS measurements, about 12 µCi 22Na activity was deposited and dried on a thin aluminium foil and was covered with an identical foil. This assembly was used as the positron source. The source correction was determined using a properly annealed defect- free aluminium sample. The PALS system used was a standard fast-fast coincidence set-up with two identical 1-inch tapered off BaF2 scintillator detectors fitted with XP 2020Q photomultiplier tubes. The time resolution obtained using 60Co source with 22Na gates was 285 ps. All lifetime spectra were analysed using PATFIT 88[11] programme.

 

Method of Analysis

The Rietveld's powder structure refinement analysis12-15 of X-ray powder diffraction data using the Rietveld software MAUD 2.26[9] revealed the refined structural parameters, such as atomic coordinates, occupancies, lattice parameters, thermal parameters etc. and microstructural parameters, such as crystallite size and r.m.s. lattice strain etc. The experimental profiles were fitted with the most suitable pseudo-Voigt (pV) analytical function15 because it takes individual care for both the crystallite size and strain broadening of the experimental profiles. Positron annihilation lifetime data were deconvoluted with three lifetime components using the PATFIT programme. A total source correction of 10% had been deducted while analysing the spectra.

 

Results and Discussion

X-ray diffraction analysis

Figure 1 shows the X-ray powder diffraction patterns of unmilled and ball milled mixture of ZnO and α-Fe2O3 powders milled for different durations. The powder pattern of unmilled mixture contained only the individual reflections of ZnO and α-Fe2O3 phases and the precursors were free from impurities. It was noticed that in the course of ball milling the mixture, ZnFe2O4 phase was formed and its amount increased gradually with increasing milling time. After 30 minutes milling, the formation of ZnFe2O4 was noticed clearly with the appearance of (220) (2θ = 29.95°) and strongest (311) (2θ = 35.3°) reflections in the XRD pattern. It may also be noticed that the content of ZnO phase was reduced to a large extent in comparison to α-Fe2O3 phase. It indicated that the ZnO phase was much prone to deformation fault as all the reflections were sufficiently broadened due to reduction in particle size and accumulation of lattice strain in the course of milling. As a result, solid-state diffusion between ZnO and α-Fe2O3 nanoparticles enhanced extensively with increasing milling time. In the course of further milling, broadening as well as degree of overlapping of neighbouring reflections were increased with increasing milling time. After 2.5 hours of milling, except the strongest (104) (2θ = 33.18°) reflection, all other reflections of α-Fe2O3 were disappeared completely in the XRD pattern. It was the indication of either (i) significant reduction in content of the phase or (ii) significant increase in peak-broadening due to reduction in particle size and accumulation of lattice strain aroused from the high energy impact of milling or due to both these effects. Concurrently, an incredible change in intensities of ZnFe2O4 reflections was observed in the XRD pattern. The rate of mechanosynthesis of ZnFe2O4 was then increased rapidly in course of milling. After 6.5 hours milling, all reflections of starting precursors were completely disappeared and it appeared that the ZnFe2O4 phase was completely grown up as the intensity distribution of the XRD pattern resembled perfectly in accordance with the ICDD PDF # 22-1012. It was reported earlier that the present process of mechanosynthesis of zinc ferrite at room temperature by ball milling may also lead to the formation of a metastable inverse spinel structure. The inverse spinel zinc ferrite was derived by a substantial displacement of Fe3+ cations to tetrahedral (T) sites and equal amount of Zn2+ cations into octahedral (O) sites of the cubic close-packed anionic sublattice. After 2.5 hours milling, when the formation of ZnFe2O4 was almost completed, verification for formation of inverse spinel structure along with normal spinel was made because without considering the inverse spinel, fitting quality of XRD powder data was poor.

Figure 2 shows the comparison of the quality of profile fitting in the 2.5 hours and 10 hours-milled samples with and without consideration of inverse spinel structure of ZnFe2O4. It is evident from the figure that the inclusion of the inverse spinel with normal one improved the profile fitting quality significantly. It indicated towards the coexistence of both the normal and inverse spinel structures of ZnFe2O4 in the samples prepared with higher milling time.

Figure 3 shows the dependence of relative phase abundances of different phases with increasing milling time. The content (wt. (%)) of ZnO decreased rapidly but that of α-Fe2O3 decreased slowly with increasing milling time. After 2.5 and 6.5 hours milling, ZnO and α-Fe2O3 phases disappeared completely from the respective XRD patterns. The formation of normal spinel phase was noticed after 30 minutes of milling and its content increased continuously (9.4-41.3 wt. (%)) up to 2.5 hours milling. Simultaneously, almost equal amount (40.1 wt. (%)) of inverse spinel phase was noticed to form. Content of normal phase remained almost unchanged up to 10 hours milling but that of inverse phase increased to a large extent after 6.5 hours milling when α-Fe2O3 was completely utilized for making exact stoichiometric (1:1 mol%) Zn-ferrite phase. This sudden increase in inverse phase content was essentially due to occupancy of octahedral vacancies by Fe3+ cations of inverse spinel structures16.

 

 

The nature of variation of crystallite sizes (D) of ZnO, α-Fe2O3 and ferrite phases are shown in Figure 4a. Crystallite sizes of all the phases were considered to be isotropic. The crystallite size of ZnO decreased rapidly to ~11 nm within 30 minutes milling and then decreased slowly with increasing milling time. Crystallite size α-Fe2O3 phase was reduced from ~75 nm to ~13 nm within 30 minutes ball milling and then decreased slowly to a saturation value ~7 nm. Normal ZnFe2O4 phase was formed after 30 minutes of milling with ~19 nm crystallite size and then reduced to ~10 nm after 2.5 hours of milling. This decrease in normal spinel crystallite size was manifested in the growth of very small crystallite size of ~6 nm of inverse spinel phase after 2.5 hours of milling. After complete formation of both spinel phases, heat energy produced by high energy impacts was utilised to release the accumulated strain inside the nanoparticles and as a result, the crystallite size increased to a small extent. This effect may also be explained as the agglomeration of nanometric particles by re-welding mechanism.

The nature of variation of r.m.s lattice strain with increasing milling time is shown in Figure 4b. It is obvious that the strain value of normal spinel lattice increased very rapidly within 30 minutes and then increased very slowly with increasing milling time, but after 2.5 hours milling, lattice strain suddenly started to release and after 5 hours milling it reached a saturation value and remained almost unchanged with increasing milling time. This nature of variation corroborates the variation of particle size with increasing milling time.

Positron annihilation spectroscopy

Positron annihilation lifetime spectra (PALS) of all the ball-milled samples were deconvoluted using three lifetime-components τ1, τ2 and τ3 with corresponding intensities I1, I2 and I3 respectively. In general, for bulk polycrystalline samples, the shortest positron lifetime (τ1) is assigned to the free annihilation of positrons at defect-free sites, the intermediate lifetime (τ2) is assigned to the lifetime of positrons trapped at the defect sites (mono or di-vacancies) and the longest one (τ3) to the pick-off annihilation of o-Ps atoms formed inside large voids. But, in the case of polycrystalline samples with crystallite size of the order of a few tens of nanometers, the assignments of lifetime parameters are different. When crystallite sizes become smaller compared to the mean positron diffusion length (~100 nm), positrons pass through the grains and annihilate mainly in the grain boundary regions. In this case, the shortest lifetime (τ1) represents the weighted average of positron lifetimes at the grain boundary defects (mono or di-vacancies) and the lifetime corresponding to annihilation with free electrons residing at the grain boundaries. The intermediate lifetime τ2 corresponds to annihilations at the triple junctions. Triple junctions are open volumes present at the intersection of three or more grain boundaries and their size is of the order of 8-10 missing atoms. The fraction of positrons that gets annihilated inside a grain depends on the positron 'trapping capability' of the defects present at the grain boundaries. The more is the positron trapping at grain boundary defects, the less is the probability of annihilation with free electrons inside a grain.

Figure 5 shows the variation of mean positron lifetime (τm) for ZnO + α-Fe2O3 nanocomposites as a function of ball-milling duration. τm is related to average defect density and is defined as

 

 

For a given type (size) of defect, higher the defect density larger will be the value of τm. It is clear that τm remains more or less constant throughout the ball-milling process. This implies that the overall defect density, as seen by the positrons remains more or less constant although ball-milling is expected to introduce several changes in structure and phase content of the nanocomposite. To obtain an insight into the positron capture mechanism at various trapping centers, the variation of the individual lifetime parameters with ball-milling duration is to be seen.

Figure 6 shows the variation of positron lifetimes τ1, τ2 and intensity I2 with milling duration. Before milling, i.e. at milling duration t = 0 hours, lifetime τ1 may be ascribed to free annihilation of positrons in the bulk α-Fe2O3 and ZnO crystals. Positron trapping can be also expected in the grain boundary defects of α-Fe2O3 as its average crystallite size is ~75 nm (Figure 4a), which is less than positron diffusion length in solids. τ2 may be ascribed to the lifetime of positrons trapped at the defect sites viz., grain boundary defects in ZnO and triple-junctions of α-Fe2O3 nanocrystals.

After 30 minutes of milling, the XRD analysis shows that the sample contains nanocrystals of ZnO, α-Fe2O3 and ZnFe2O4 (Figures 1 and 3). Hence the shortest lifetime (τ1) after milling the mixture for 30 minutes, may be assigned to the mixed lifetime of positrons trapped at defects (mono or di-vacancies) present at the grain boundaries of all the nanocrystals (ZnO, α-Fe2O3 and ZnFe2O4) present in the sample and also to the free annihilation inside these nanocrystalline grains. The intermediate lifetime (τ2) has been assigned to the lifetime of positrons trapped at the triple junction formed at the intersection of three or more nanocrystalline grains. Similar assignment of lifetimes was seen to successfully explain the results obtained in case of ZnO nanocrystals17. τ1 and τ2 in this case can be expressed as follows

and,

where i stands for species of nanocrystalline grains viz., ZnO, α-Fe2O3 or ZnFe2O4, GB stands for grain-boundary, B stands for bulk i.e. the defect free sites inside a grain, n stands for the fraction of positrons annihilating in a particular site (GB or B) of a particular type of grain i, and TP stands for triple junction.

Upto a milling duration of 3.5 hours, τ1 shows a slight increasing trend. This is because of the fact that the crystallite sizes become smaller as the ball-milling duration is increased which in turn increases the surface area to volume ratio. As a result of this, the contribution from the annihilation at grain boundary increases resulting in the increase in the value of τ1. However τ2 shows a decreasing trend as milling hour changes from 30 minutes to 2.5 hours and remains more or less same up to milling time 3 hours. The fall in τ2 may be ascribed to the decrease in weight percentage of α-Fe2O3 (from 58.6 to 14.9%) with increase in milling time. Disappearance of α-Fe2O3 from the sample in the course of milling indicates decrease in the contribution of triple junctions to τ2. It can also be seen that the corresponding intensity I2, which gives an idea of the defect concentration related to τ2, maintains a steady value up to 3.5 hours of milling. Another point to be noted is that though the wt. (%) of α-Fe2O3 phase changes, τ1 still remains more or less constant in the region of 2.5 to 3.5 hours of milling. A decrease in wt. (%) could have decreased the contribution of annihilation at triple junction of α-Fe2O3 nanocrystals. This may be explained as following manner. Since the crystallite size of α-Fe2O3 were also seen to decrease in this region, the fraction of positrons annihilating at the grain boundaries becomes more, thereby compensating the decrease in the value of τ1. However τ1 was seen to decrease after milling the sample for 5 hours. The phase content of the sample at this stage includes α-Fe2O3, normal spinel ZnFe2O4 and inverse spinel ZnFe2O4 phases. Their content and grain sizes do not change much in this milling stage compared to the previous milling stage (3.5 hours). The only change is a substantial decrease in r.m.s. strain value of the normal spinel ZnFe2O4 structure18. Lattice strain depends on the lattice imperfection and lattice imperfections may arise from the non-stoichiometric composition of the crystal structure. In this case, decrease in r.m.s. strain value of the normal spinel ferrite structure strongly indicates the removal of vacancy type defects, which arises out of non-stoichiometric composition in the ferrite structure. The removal of defects, from the normal spinel structure, leads to a decrease in contribution of defect related lifetime to the lifetime component τ1, leading to a reduction in its value. After 6.5 hours of milling, τ1 again increases. At this milling stage, it was seen that the α-Fe2O3 phase completely vanishes from the composite sample. Simultaneously, the wt. (%) of inverse spinel component also increases after milling the sample for 6.5 hours. This indicates a reduction in the contribution of the free annihilation of positrons in the α-Fe2O3 crystals as well as of the annihilation at the grain boundary defect sites in the inverse spinel structure resulting in an overall increase in τ1 value. However, as the complete solid solution of α-Fe2O3 results in increase the wt. (%) of inverse spinel ZnFe2O4, the value of τ2 did not change at this stage.

Figure 7 shows the variation of τ3 and I3 as a function of ball-mill duration. Except an initial increase, value of τ3 remains more or less steady. However, I3 shows some fluctuations in its behavior. These fluctuations do not correlate themselves anyhow to the structural changes and hence are not discussed in this article.

 

Conclusion

The quantitative analysis of the XRD data evaluated on the basis of Rietveld's powder structure refinement method yielded detailed information about the structure and microstructure of mechanosysthesized nanoscale zinc ferrite as well as the distribution of cations in the spinel ferrite. The main feature of the structural disorder of mechanosynthesized zinc ferrite was the defect induced inverse spinel phase transition. The degree of inversion increased rapidly with increasing milling time and then wt. (%) of inverse phase approached towards a saturation value. Positron annihilation lifetime data shows that the mean lifetime τm does not change much with ball milling durations. This implies that the overall defect density, as seen by the positrons remains more or less constant during milling. The variation of individual lifetimes τ1 and τ2 and corresponding intensities I1 and I2 shows evaluation of different phases with milling duration and confirms the formation of inverse spinel ferrite structure.

 

Acknowledgement

One of the authors, SKP, wishes to thank the University Grant Commission (UGC) India, for granting CAS-I programme under the thrust area "Condensed Matter Physics including Laser applications" to the Department of Physics, The University of Burdwan under the financial assistance of which the work has been carried out. We are also grateful to UGC-DAE-CSR, Kolkata centre for providing Positron annihilation spectroscopy facility.

 

References

1. Goldman A. Modern Ferrite Technology. Van Nostrand Reinhold: New York; 1990.         [ Links ]

2. Berskovsky BM, McDvcdcv VF and Krakov MS. Magnetic Fluids: Engineering Applications. Oxford: Oxford University Press; 1993.         [ Links ]

3. Ayala RE and Marsh DW. Characterization and long range reactivity of zinc ferrite in high temperature desulfurization processes. Industrial & Engineering Chemistry Research. 1991; 30:55-60. http://dx.doi.org/10.1021/ie00049a009        [ Links ]

4. Bissett LA and Strickland LD. Analysis of a fixed-bed gasifier IGCC configuration. Industrial & Engineering Chemistry Research. 1991; 30:170-176. http://dx.doi.org/10.1021/ie00049a025        [ Links ]

5. Sack RO and Ghiorso MS. An internally consistent model for the thermodynamic properties of Fe_Mg titanomagnetite aluminate spinels. Contributions to Mineralogy and Petrology. 1991; 106:474-505. http://dx.doi.org/10.1007/BF00321989        [ Links ]

6. Druska P, Steinike U and Sepelak V. Surface structure of mechanically activated and of mechanosynthesized zinc ferrite. Journal of Solid State Chemistry. 1999; 146:13-21. http://dx.doi.org/10.1006/jssc.1998.8284        [ Links ]

7. Chinnasamy CN, Narayanasamy A, Ponpandian N, Chattopadhyay K, Guerault H, and Greneche J-M. Magnetic properties of nanostructured ferrimagnetic zinc ferrite. Journal of Physics Condensed Matter. 2000; 12:7795-7805. http://dx.doi.org/10.1088/0953-8984/12/35/314        [ Links ]

8. Lutterotti L, Scardi P and Maistrelli P. LSI- a computer program for simultaneous refinement of material structure and microstructure. Journal of Applied Crystallography. 1992; 25:459-462. http://dx.doi.org/10.1107/S0021889892001122        [ Links ]

9. Lutterotti L. MAUD. version 2.26. 2011. Available from: <http://www.ing.unitn.it/~maud>         [ Links ].

10. Krause-Rehberg R and Leipner HS. Positron Annihilation in Semiconductors: Solid-State Sciences. Berlin: Springer; 1999. p. 127.         [ Links ]

11. Kirkegaard P and Eldrup M. POSITRONFIT: a versatile program for analysing positron lifetime spectra. Computer Physics Communications. 1972; 3:240-255. http://dx.doi.org/10.1016/0010-4655(72)90070-7        [ Links ]

12. Rietveld HM. Line profiles of neutron powder-diffraction peaks for structure refinement. Acta Crystallographica. 1967; 22:151-152. http://dx.doi.org/10.1107/S0365110X67000234        [ Links ]

13. Rietveld HM. A profile refinement method for nuclear and magnetic structures. Journal of Applied Crystallography. 1969; 2:65-71. http://dx.doi.org/10.1107/S0021889869006558        [ Links ]

14. Sain S, Patra S and Pradhan SK. Microstructure and optical band-gap of mechano-synthesized CdxZn1-xS quantum dots. Journal of Physics D. 2011; 44: Article ID 075101.         [ Links ]

15. Young RA and Willes DB. Profile shape functions in Rietveld refinements. Journal of Applied Crystallography. 1982; 15:430-438. http://dx.doi.org/10.1107/S002188988201231X        [ Links ]

16. Wei Q-M, Li J-B and Chen Y-J. Cation distribution and infrared properties of NixMn1 -xFe2O4 ferrites. Journal of Materials Science. 2001; 36:5115-5118. http://dx.doi.org/10.1023/A:1012473207424        [ Links ]

17. Mishra AK, Chaudhuri SK, Mukherjee S, Priyam A, Saha A and Das D. Characterization of defects in ZnO nanocrystals: Photoluminescence and positron annihilation spectroscopic studies. Journal of Applied Physics. 2007; 102:103514-103514-6. http://dx.doi.org/10.1063/1.2817598        [ Links ]

18. Bid S and Pradhan SK. Preparation of zinc ferrite by high-energy ball-milling and microstructure characterization by Rietveld's analysis. Materials Chemistry and Physics. 2003; 82:27-37. http://dx.doi.org/10.1016/S0254-0584(03)00169-X        [ Links ]

 

 

Received: February 20, 2012
Revised: July 20, 2012

 

 

* e-mail: skp_bu@yahoo.com