Acoustic study of nano-crystal embedded PbO P 2 O 5 glass

Similar documents
STRONG CONFIGURATIONAL DEPENDENCE OF ELASTIC PROPERTIES OF A CU-ZR BINARY MODEL METALLIC GLASS

Chapter 3 Chapter 4 Chapter 5

T. Egami. Model System of Dense Random Packing (DRP)

12. Spectral diffusion

Chapter 12. Nanometrology. Oxford University Press All rights reserved.

Kinetically arrested long-range magnetic ordered phase. Alok Banerjee

The effect of rigidity on torsional vibrations in a two layered poroelastic cylinder

International Journal of Pure and Applied Sciences and Technology

Sound Attenuation at High Temperatures in Pt

STUDIES ON ZnS - CuS NANOPARTICLE SYSTEM.

Physics of disordered materials. Gunnar A. Niklasson Solid State Physics Department of Engineering Sciences Uppsala University

Optical and Photonic Glasses. Lecture 39. Non-Linear Optical Glasses III Metal Doped Nano-Glasses. Professor Rui Almeida

Thermal properties of Engineering Materials

Kinetics. Rate of change in response to thermodynamic forces

Preparation of Silver Nanoparticles and Their Characterization

Raman spectroscopy of self-assembled InAs quantum dots in wide-bandgap matrices of AlAs and aluminium oxide

Chapter 4. The Effect of Elastic Softening and Cooperativity on the Fragility of

Shape transition in ZnO nanostructures and its effect on blue-green photoluminescence

International Journal of Pure and Applied Sciences and Technology

Theoretical approach to the Poisson s ratio behaviour during structural changes in metallic glasses. Abstract

An Introduction to Polymer Physics

Introductory Nanotechnology ~ Basic Condensed Matter Physics ~

Relaxation in Glass: Review of Thermodynamics. Lecture 11: Thermodynamics in the Glass Transition Region

Dielectric behaviour and a.c. conductivity in Cu x Fe 3 x O 4 ferrite

Optical and Photonic Glasses. Lecture 30. Femtosecond Laser Irradiation and Acoustooptic. Professor Rui Almeida

Pre-yield non-affine fluctuations and a hidden critical point in strained crystals

4. Thermal properties of solids. Time to study: 4 hours. Lecture Oscillations of the crystal lattice

Ultrasonic investigations of some bismuth borate glasses doped with Al 2 O 3

arxiv:cond-mat/ v2 [cond-mat.mtrl-sci] 26 Nov 1999

Nanophysics: Main trends

Impact of Magnetic Impurities on Transient Propagation of Coherent Acoustic Phonons in II-VI Ternary Semiconductors

Structure and Dynamics : An Atomic View of Materials

Bandgap engineering through nanocrystalline magnetic alloy grafting on. graphene

Concepts in Surface Physics

Potential-dependent dynamic fracture of nanoporous gold

Electric field dependent sound velocity change in Ba 1 x Ca x TiO 3 ferroelectric perovskites

Optimisation using measured Green s function for improving spatial coherence in acoustic measurements

The electronic structure of materials 1

Local Electronic Structures and Chemical Bonds in Zr-Based Metallic Glasses

Thermally Functional Liquid Crystal Networks by Magnetic Field Driven Molecular Orientation

GeSi Quantum Dot Superlattices

Dynamics of Supercooled Liquids The Generic Phase Diagram for Glasses

Lecture 11 - Phonons II - Thermal Prop. Continued

Quasi-Harmonic Theory of Thermal Expansion

Impedance Analysis and Low-Frequency Dispersion Behavior of Bi 4 Ti 3 O 12 Glass

STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY

File name: Supplementary Information Description: Supplementary Figures, Supplementary Tables and Supplementary References

Widely Tunable and Intense Mid-Infrared PL Emission from Epitaxial Pb(Sr)Te Quantum Dots in a CdTe Matrix

Melting and solidi cation of Pb nanoparticles embedded in an Al matrix as studied by temperature-modulated di erential scanning calorimetry

Chapter 5: Molecular Scale Models for Macroscopic Dynamic Response. Fluctuation-Dissipation Theorem:

SAMPLE ANSWERS TO HW SET 3B

Relevance of jamming to the mechanical properties of solids Sidney Nagel University of Chicago Capri; September 12, 2014

Effective Temperatures in Driven Systems near Jamming

Conductivity studies of lithium zinc silicate glasses with varying lithium contents

Structural and Mechanical Properties of Nanostructures

Crystalline Surfaces for Laser Metrology

Sub -T g Relaxation in Thin Glass

dc Conductivity and dielectric properties of V 2 O 5 -Bi 2 O 3 -ZnO glass

Scanning Force Microscopy II

Flow of Glasses. Peter Schall University of Amsterdam

Chapter 6 ELECTRICAL CONDUCTIVITY ANALYSIS

Acoustic Velocity, Impedance, Reflection, Transmission, Attenuation, and Acoustic Etalons

Tb 2 Hf 2 O 7 R 2 B 2 7 R B R 3+ T N

Introduction to thermodynamics

Intensity / a.u. 2 theta / deg. MAPbI 3. 1:1 MaPbI 3-x. Cl x 3:1. Supplementary figures

AM11: Diagnostics for Measuring and Modelling Dispersion in Nanoparticulate Reinforced Polymers. Polymers: Multiscale Properties.

Chapter 10. Nanometrology. Oxford University Press All rights reserved.

Solid-State Diffusion and NMR

An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation

Saveetha Engineering College, Thandalam, Chennai. Department of Physics. First Semester. Ph6151 Engineering Physics I (NOV/DEC 2014)

arxiv: v1 [cond-mat.other] 11 Sep 2008

Chapter 2. Dielectric Theories

Chapter 2. Thermodynamics of pure metals and alloys

Temperature-Modulated Differential Scanning Calorimetry Analysis of High- Temperature Silicate Glasses

2008,, Jan 7 All-Paid US-Japan Winter School on New Functionalities in Glass. Controlling Light with Nonlinear Optical Glasses and Plasmonic Glasses

Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi

Elastic and Inelastic Scattering in Electron Diffraction and Imaging

Physical and Dielectric Properties of Silver Lithium Niobate Mixed Ceramic System

G. R. Strobl, Chapter 5 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). J. Ferry, "Viscoelastic Behavior of Polymers"

7.2.1 Seismic waves. Waves in a mass- spring system

Supplementary Figure 1: Approach to the steady state. Energy vs. cycle for different

Supplementary Figures:

Indentation cracking of glass II: Relationship between crack resistance and equibiaxial flexural strength

Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).

Doppler echocardiography & Magnetic Resonance Imaging. Doppler echocardiography. History: - Langevin developed sonar.

SUPPLEMENTARY INFORMATION

Higher Order Elastic Constants of Thorium Monochalcogenides

Chapter 4 Section 4.3

Metal nanoparticle-doped coloured films on glass and polycarbonate substrates

Quiz 1 Introduction to Polymers

3.320 Lecture 23 (5/3/05)

ELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS

Synthesis and Characterization of Exfoliated Graphite (EG) and to Use it as a Reinforcement in Zn-based Metal Matrix Composites

Microscopical and Microanalytical Methods (NANO3)

Surface Studies by Scanning Tunneling Microscopy

Conduction. Metals, Semiconductors and Interconnects. Fig 2.1

Differential Acoustic Resonance Spectroscopy Analysis of Fluids in Porous Media

VIII. Rubber Elasticity [B.Erman, J.E.Mark, Structure and properties of rubberlike networks]

Welcome to Physics 460 Introduction to Solid State Physics

THz nanocrystal acoustic vibrations from ZrO 2 3D supercrystals

Transcription:

Bull. Mater. Sci., Vol. 9, No. 4, August 6, pp. 357 363. Indian Academy of Sciences. Acoustic study of nano-crystal embedded PbO P O 5 glass SUDIP K BATABYAL, A PAUL, P ROYCHOUDHURY and C BASU* Department of Solid State Physics, Indian Association for the Cultivation of Science, Kolkata 7 3, India Department of Physics, New Alipore College, Kolkata 7 53, India MS received 6 December 5 Abstract. Nano-crystal embedded PbO P O 5 glass has been prepared and characterized by XRD and TEM measurements. The ultrasonic velocity and attenuation measured within the temperature range 8 3 K show significant structure and interesting feature with the presence of nano-crystalline region. The glass samples were prepared by melt-quench method and nano-crystals of different sizes were produced by heat treatment of the glasses for different durations of heating. All the processes were carried out at or above glass transition temperature. A theoretical model that takes account of the effects of thermally activated relaxation, anharmonicity as well as microscopic elastic inhomogeneities arising out of fluctuations has been successfully applied to interpret the variation of ultrasonic velocity and attenuation data. An interesting outcome of this application has been to propose a method for the determination of the size of nano-crystals from the ultrasonic attenuation data. Keywords. Glasses; acoustical properties; nanostructured materials; glass ceramic. 1. Introduction During the last two decades, studies of different physical properties of nano-crystal embedded glass matrix have attracted attention of technologists as well as scientists for fabrication of glass ceramic through controlled crystallization. The studies usually involve X-ray diffraction measurement (XRD (Komatsu et al 1991; Balaya and Sunandana 1993, electrical conductivity measurement (Chatterjee and Chakravorty 199; Maity et al 1996, optical (Kraevskii and Solinov 3 and thermal (Balaya and Sunandana 1993 properties measurements etc. In addition, the differential scanning calorimetry measurements give us information about the mechanism of crystallization and activation energy for crystal growth in glass. Further, we may also mention that after different thermal treatments, crystallization behaviour of TeO LiNbO glasses by determining the phases using X-ray diffraction was reported (Komatsu et al 1991. The magnetization measurement of vitroceramics containing lithium ferrite crystal showed closed link to the crystal sizes by annealing of LiFe 5 O 8 glass (Chaumont and Benier 1983. The effects of crystallization on electrical conductivity of 3Li O 7TeO glass were also studied (Balaya and Sunandana 1993. The main aspects of such studies of nano particles are two folds, viz. (i to observe the change of different physical properties due to nano size of the materials and (ii to develop various processes of tuning these properties with *Author for correspondence (sspcb@iacs.res.in various sizes of particles. Heat treatment of glassy materials at or above glass transition temperature is relatively an easy process to produce nano-crystals within the glass matrix. In this process one can control the size of the produced nano particles by controlling the heat treatment temperature and duration of thermal treatment. A large number of workers have generated nano-crystals by this process and studied their different physical properties. But ultrasonic characterization of such materials is yet to be done. In this communication, we present the different phases with varying sizes of nano-crystals developed during heat treatment of (PbO 5 (P O 5 5 glass and investigate the variation of ultrasonic velocity and attenuation with the size of nano-crystals and finally propose an empirical relation between attenuation parameter and size of nano-crystals.. Experimental The (PbO 5 (P O 5 5 glass samples were prepared by melt quench method as described earlier (Mukherjee et al 1993; Maiti et al 1996; Paul et al 1997, 1998. The glass transition temperatures (T g of the samples were estimated from differential thermal analysis (DTA. The glass ceramic phases were developed within the samples by heat treatment. For this purpose, two temperatures were considered, one the glass transition temperature, 443 K and the other, 3 K above the glass transition temperature i.e. 473 K. Keeping the temperatures fixed at 443 K and 473 K, the samples were heated for different durations of time: 357

358 Sudip K Batabyal et al 3 min, 1 h and 9 min. The development of crystal phases was confirmed from XRD. To ensure the complete formation of crystal phases one of the glass samples was heattreated for a period of h at 673 K which is above the crystallization temperature. It was then taken as the reference polycrystalline sample for a comparison of particle size and degree of crystallization by transmission electron microscopy (TEM and XRD, respectively of the samples undergoing different heat treatments. Figure 1 shows XRD pattern of all the samples including the reference polycrystalline one. TEM micrographs were also taken for all glass ceramic samples and are presented in figure. For ultrasonic measurements bulk samples with exactly parallel and polished opposite faces are prepared. Ultrasonic velocity and attenuation measurements were carried out at a frequency of 1 MHz in the temperature range 8 3 K. The polycrystalline sample was very much fragile and so the ultrasonic measurements could not be performed. Details of the experimental setup and procedure for measurements are given elsewhere (Mukherjee et al 1993; Maiti et al 1996; Paul et al 1997, 1998. 3. Results and discussion From TEM micrographs we observed nano sized crystals developed by heat treatment within the glass matrix for all the samples excepting the one heat treated at 443 K for 3 min. The average sizes of the nano-crystals have been calculated from TEM micrographs. The particle sizes of the samples heat treated at different temperatures are given in table 1. From the data it is clear that the size of the crystals increases with the duration of heat treatment. Again the XRD spectra show a hump around θ = 5 for the samples undergoing heat treatment up to 1 h at different fixed temperatures. But when the heat treatment time is 9 min, a large number of crystalline peaks indicating Intensity (a.u. 11 (P 4 O 7 11 (P O 5 1 (P O 5 1 (P O 5 41 (P 4 O 7 1 3 4 5 6 7 θ Figure 1. XRD spectra of the glass ceramic samples. Reference S 6 S 5 S 4 S 3 S S 1 aggregate formation have been observed. XRD data, however, do not lead to any conclusive information about the structure of the nano aggregates formed. The XRD profiles S3 and S6 indicate the presence of P O 5 in the studied samples. Thus, in addition to entering into the aggregate formation, P O 5 exists also in the unreacted form. Obviously XRD data do not suggest any definite composition. Characteristic features of the composition and related structural details by suitable methods may reveal the true picture. Variations of longitudinal ultrasonic velocity and attenuation coefficient with temperature have been presented in figures 3 and 4, respectively. The variation of longitudinal velocity with temperature of all the samples shows a negative slope. The value of ultrasonic velocity is almost same at different temperatures for all the samples except the one heat treated at 473 K for 9 min in which case it is slightly lower. A broad peak has been observed for all the samples in the temperature variation of ultrasonic attenuation. The peaks become more flat and the values become lower than that of the sample heat-treated at 443 K for 3 min where no crystals were observed. To explain the variation of ultrasonic velocity and attenuation with temperature for nano-crystal embedded glass, we consider a model which includes the effects of thermally activated relaxation in a double well created by the glassy disorder, anharmonicity as well as microscopic elastic inhomogeneities generated by freezing of different natural fluctuations. The scenario envisaged by the model is that due to inherent disorder in glass, atoms or a group of atoms constituting localized structural defects can have different configurational states and may move in a double well potential having two minima corresponding to two local equilibrium configurations. At low temperature these particles can only move from one well to the other by quantum mechanical tunneling. At high temperature (range of our study the two defect configurations should actually be treated as two thermodynamical states involving a thermodynamical average over some vibrational degrees of freedom and the double well description is interpreted as free energy curve rather than potential energy curve. When a sound wave passes through such material, the strain field causes a relative change in free energy of the two states corresponding to the two possible configurations and thus produces thermal disequilibrium. The change in free energy due to the strain field is characterized by a deformation potential, D, which may be defined as the relative energy shift of the two states for a strain field of unit strength. A relaxation process follows to restore the equilibrium configuration by surmounting the barrier height via a thermally activated process. We can calculate the change, dv, of the sound velocity due to relaxation effect from its value in absence of relaxation (i.e. at K and the dissipation factor, often called internal friction, Q 1 (related to the attenuation coefficient, α, by Q 1 = α v/ω arising from the relaxation effect and obtained the following expressions (Jäckle et al 1976

Acoustic study of nano-crystal embedded PbO P O 5 glass 359 Figure. TEM micrographs of glass ceramic samples. Table 1. Sample specifications and particle size of the nano-crystals. Sample Temperature of heat Duration of heat Particle size as name treatment (K treatment obtained from TEM (nm S1 443 3 min S 443 1 h 13 6 S3 443 9 min 17 S4 473 3 min 14 S5 473 1 h S6 473 9 min 49

36 Sudip K Batabyal et al v l (1 3 m/s 3.3 s1 s s3 3.7 s4 s5 3. s6 3.17 3.1 5 1 15 5 3 Figure 3. Variation of longitudinal ultrasonic velocity with temperature (The solid lines are fit to aid the eye. δv = nd v 8ρvkT 1+ ωτ Q 1 = nd 1 vkt + ωτ 4ρ 1 ωτ, (1, ( where τ (relaxation time = τ e V/kT, n the number of relaxing particles per unit volume, ρ the density of material, v the sound velocity in absence of relaxation (i.e. at K and D the deformation potential. The above relations refer to one kind of defect particle only with a single value of τ and D. In glass there will be several kinds of defects having several values of τ (hence different values of activation energy, V, i.e. barrier height between the two wells in a double well and different values of deformation potential, D. The asymmetry parameter,, is defined by the difference in energy of the minima of the two wells. The distribution of represented by f ( is considered to be constant (denoted by f (Gilroy and Phillips 1981, an assumption valid below the glass transition temperature. Further the distribution of V is taken in the form, g(v = 1/V exp( V/V, where V is the energy value of the order that corresponds to T g (Gilroy and Phillips 1981. Introducing the above distributions in (1 and (, an approximate solution in closed form is obtained as follows (Gilroy and Phillips 1981 δv D v 4ρv Q = f 1 = D 4ρ v kt V [1 ( ωτ ] kt kt V π ( ωτ. V f, (3 The theory described above, which accounts for the attenuation behaviour, is found to be inadequate to explain the velocity variation with temperature. It should be noted (4 that ωτ << 1 and from (3 δv decreases as T increases. Evidently, the relaxation effect alone cannot explain the high temperature variation of sound velocity in some glasses where the velocity rises after attaining a minimum. For such cases an improvement can be made by introducing the effect of anharmonicity as well as the microscopic elastic inhomogeneities that exist in glass. Using an extension of the quasiharmonic continuum model of Garber and Granato (1975 to isotropic material (Claytor and Sladek 1978, the anharmonic effect on elastic constant, C, at any temperature, T, can be expressed as C = C [1 ΓF(T/Θ], (5 where C is the elastic constant at K, the factor Γ determines the anharmonic effect, Θ the Debye temperature at K and the function, F(T/Θ, is given by 4 Θ/ T 3 T T x dx F = 3. Θ Θ x e 1 (6 Then the anharmonic effect of sound velocity is immediately obtained using (5 as follows 1 T vt = v 1 ΓF. Θ (7 The microscopic inhomogeneities referred to earlier are associated with the freezing of natural fluctuation of the order parameter of liquid glass at temperature close to T g during glass formation. When the rate of cooling is high for glass formation the order fluctuations are not dispersed but are frozen at T g and are retained in the solid state. Details of the calculation are given in our earlier works (Paul et al 1998, 1999. Considering the effect of relaxation, anharmonicity and frozen-in fluctuation, the expression for velocity variation becomes R = v T v8 v 8 fd 1 P 8 P T H 8 H T ρv 4 [ ] + Γ[ ] + A B A B S 8 + S 8 S T + S T 1 v v v v ρ fd 1 fd 1 1 P 8 ΓH8 1 P T ΓH 4ρv 4ρv = A B S 8 + S 8 1 1 v v 8 8 ρv ρ 1 P 8 ΓH8 4ρv fd 1 P Γ H 4 fd 1 T, (8

Acoustic study of nano-crystal embedded PbO P O 5 glass 361 1. 1. S 1 S.. 1. 1 15 5 3 5 1 15 5 3 1. S 3 S 4.. 5 1 15 5 3 1. 5 1 15 5 3 S 5 S 6.. 5 1 15 5 3 5 1 15 5 3 Figure 4. Variation of longitudinal ultrasonic attenuation coefficient with temperature (The solid lines are fit to aid the eye. Table. Values of the theoretical parameters for different glass ceramic samples. V f D /v A/v B/v Sample τ (s (1 3 J (1 3 Jm 5 s Γ (1 5 (1 7 Q s S 3 1 1 83 1 893 S3 3 1 1 99 9 9799 S4 1 13 3 1 5 1 37 19 9938 S5 3 79 1 6 99 37 9859 S6 3 7 1 64 1 9893

36 Sudip K Batabyal et al B / where P = 1 ( ωτ K T V B8/, P = 1 ( ωτ K V, H T = T 8 F(T/Θ; H 8 = F(8/Θ, S T = (T g T and S 8 = (T g 8. We would like to mention here that in case of attenuation the major contribution comes from the thermally activated relaxation process. In an earlier communication we have estimated the contribution of anharmonicity and the frozen-in fluctuation towards attenuation and have shown that these are negligible (Paul et al 1999. Transforming our attenuation data into internal friction, Q 1 Q -1 (1 3 8 4 1 15 5 3 Figure 5. Theoretical curve (solid line of Q 1 (longitudinal against temperature and experimental points ( are also shown for the sample, S4. and taking ( f D /v, τ and V of (3 and (4 as parameters, we fit the experimental data by similar procedure as elaborated in our earlier publications (Paul et al 1997, 1999. Besides the parameters ( f D /v, τ and V already evaluated from attenuation data we treat Γ, A/v and B/v in (8 as additional parameters to fit the velocity data. The evaluated parameters are given in table. Figures 5 and 6 show a representative fit of experimental data with the theory for attenuation and velocity, respectively. It is expected that the mechanical stress at the crystal glass interface is likely to influence sound propagation. Also, the wave propagation will depend on crystalline fraction and the number of nano-crystals formed. However, the changes in the ultrasonic velocities for the samples, S1 to S5, throughout the range of temperature studied are rather small and do not appear to follow a regular order. While the marked change is observed for sample 6 there is overlap between the curves of other samples, so with this limited observation the analysis of the velocity data is not attempted. From the analysis of the Q 1 data it was found that V is the same for all the samples whereas the values of f D / v go on decreasing with the size of the nano-crystals as shown in figure 7. The constant value of V indicates that the distribution of barrier height has not been affected by the presence of small quantities of nano-crystals. The change in the values of the ultrasonic velocities of different samples is very small throughout the entire range of temperature studied, so it is expected that the change in v is also small for all the samples. This implies that the change in deformation potential, if any, is too small to be perceptible. Hence the change in values of the parameter, f D /v, is mostly due to the change in the values of f alone, from sample to sample. As the duration of heat treatment increases the size of the nano-crystals increases R - -4-6 -8 8 13 18 3 8 Figure 6. Fit to the relative change in longitudinal sound velocity (dashed line, theoretical curve considering the relaxation process alone; (dotted line, theoretical curve considering the relaxation and anharmonicity; (solid line, theoretical curve considering the relaxation and anharmonicity and frozen-in fluctuation and (, experimental points of the sample, S4. f o D / v o 1.9 1.7 1.5 1.3 1.1.9.7.5 8 13 18 3 8 33 38 43 48 53 particle size (nm Figure 7. Variation of the parameter, f D /v, with particle size of the nano-crystals.

Acoustic study of nano-crystal embedded PbO P O 5 glass 363 in the system, and the parameter, f D /v, shows a decreasing trend. It is found with the increasing size of the nanocrystals the value of the parameter, f D /v, decreases exponentially and follows the relation y = y + A exp [ x/b], where y, A and b are constants (see figure 7. For this glass system, y = 67578, A = 13 715 and b = 4 5447. We could expect that once we evaluate the values of these constants for a particular glass system in which nanocrystals are embedded, the size of the crystals within the glass matrix can be evaluated from the attenuation data. Hence this method may be used to determine size of nanocrystals using ultrasonic as a tool. References Balaya P and Sunandana C S 1993 J. Non-Cryst. Solids 16 53 Chatterjee A and Chakravorty D 199 J. Phys. D: Appl. Phys. 3 97 Chaumont C and Benier J C 1983 Solid State Commun. 48 357 Claytor T N and Sladek R J 1978 Phys. Rev. B18 584 Garber J A and Granato A V 1975 Phys. Rev. B11 399 Gilroy K S and Phillips W A 1981 Philos. Mag. B43 735 Jäckle J, Piché L, Arnold W and Hunklinger S 1976 J. Non- Cryst. Solids 365 Komatsu T, Tawarayama H, Mohri H and Matusita K 1991 J. Non-Cryst. Solids 135 15 Kraevskii S L and Solinov V F 3 J. Non-Cryst. Solids 316 37 Maiti A, Mukherjee S, Ghosh U S and Basu C 1996 J. Non- Cryst. Solids 194 13 Maity A K, Nath D and Chakravorty D 1996 J. Phys.: Condens. Matter 8 5717 Mukherjee S, Maiti A, Ghosh U S and Basu C 1993 Philos. Mag. B67 83 Paul A, Ghosh U S and Basu C 1997 J. Non-Cryst. Solids 1 65 Paul A, Chattopadhyay A K and Basu C 1998 J. Appl. Phys. 84 513 (Erratum: 1999 J. Appl. Phys. 85 465 Paul A, Maiti A and Basu C 1999 J. Appl. Phys. 86 3598