Sound Attenuation at High Temperatures in Pt

Similar documents
Higher Order Elastic Constants of Thorium Monochalcogenides

Non linear properties of mixed divalent crystals

Temperature dependent elastic and ultrasonic properties of ytterbium monopnictides

Elastic and Acoustic Properties of XB 2 (X= V, Nb & Ta) Transition Metal Diborides

ULTRASONIC CHARACTERIZATION OF CERAMIC MATERIAL TITANIUM DIBORIDE

Quasi-Harmonic Theory of Thermal Expansion

Numerical calculation of the electron mobility in ZnS and ZnSe semiconductors using the iterative method

Anisotropic assessment of ultrasonic wave velocity and thermal conductivity in ErX (X:N, As)

Lecture 11 - Phonons II - Thermal Prop. Continued

Research Article Computations of Ultrasonic Parameters in Zr 100 X Sn X Alloys

Ultrasonic attenuation in rare-earth monoarsenides

Chapter 3. The complex behaviour of metals under severe thermal and stress

Micro-canonical ensemble model of particles obeying Bose-Einstein and Fermi-Dirac statistics

Lattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W

Equations of State. Tiziana Boffa Ballaran

First-principles calculation of defect free energies: General aspects illustrated in the case of bcc-fe. D. Murali, M. Posselt *, and M.

COMPARISON OF THERMAL CHARACTERISTICS BETWEEN THE PLATE-FIN AND PIN-FIN HEAT SINKS IN NATURAL CONVECTION

Asymptote. 2 Problems 2 Methods

Propagation of Acoustic Wave inside the Carbon Nanotube: Comparative Study with Other Hexagonal Material

Non-destructive characterization of XB 2 (X= V, Nb and Ta) transition metal diborides

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

16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity

Structure and Dynamics : An Atomic View of Materials

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

Title. Author(s)Tamura, S.; Sangu, A.; Maris, H. J. CitationPHYSICAL REVIEW B, 68: Issue Date Doc URL. Rights. Type.

Chapter 5 Phonons II Thermal Properties

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

Signal Loss. A1 A L[Neper] = ln or L[dB] = 20log 1. Proportional loss of signal amplitude with increasing propagation distance: = α d

A PSEUDOPOTENTIAL STUDY ON THE THERMODYNAMIC AND ELASTIC PROPERTIES OF Pd 39 Ni 10 Cu 30 P 21 BULK METALLIC GLASS

Phonon dispersion relation of liquid metals

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

Material properties and vibration isolation Technical information

Analysis of Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch. Abstract

Two Constants of Motion in the Generalized Damped Oscillator

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

Study of full band gaps and propagation of acoustic waves in two-dimensional piezoelectric phononic plates

Thermal properties of Engineering Materials

Professor, Institute of Engineering Mechanics, Harbin. China 2. Ph.D Student, Institute of Engineering Mechanics, Harbin. China 3

Acoustic forcing of flexural waves and acoustic fields for a thin plate in a fluid

Reliability of Axially Loaded Fiber-Reinforced-Polymer Confined Reinforced Concrete Circular Columns

8. INTRODUCTION TO STATISTICAL THERMODYNAMICS

Interatomic potentials for monoatomic metals from experimental data and ab initio calculations

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

UNIVERSITY COLLEGE LONDON EXAMINATION FOR INTERNAL STUDENTS

6.730 Physics for Solid State Applications

SCIENCE CHINA Physics, Mechanics & Astronomy

Phase Transitions in Strontium Titanate

Objective decomposition of the stress tensor in granular flows

Introduction to physical acoustics

Natural convection in a vertical strip immersed in a porous medium

Study of Superconducting State Parameters of Alloy Superconductors

ELASTICITY AND CONSTITUTION OF THE EARTH'S INTERIOR*

ULTRASONIC ABSORPTION COEFFICIENT AND PROPAGATION VELOCITY IN CRITICAL ΒΕΝΖONITRILE-ISOOCTANE MIXTURE

FATIGUE DURABILITY OF CONCRETE EXTERNALLY STRENGTHENED WITH FRP SHEETS

ENERGY ANALYSIS: CLOSED SYSTEM

Finite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements

The Boguslawski Melting Model

Plane waves in a rotating generalized thermo-elastic solid with voids

Introduction. Methods of vibration control

2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson

Nanoscale interfacial heat transfer: insights from molecular dynamics

Dispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson

Motivation. Confined acoustics phonons. Modification of phonon lifetimes Antisymmetric Bulk. Symmetric. 10 nm

Pressure derivatives of bulk modulus for materials at extreme compression

Colloque National sur les Techniques de Modélisation et de Simulation en Science des Matériaux, Sidi Bel-Abbès Novembre 2009

Structural stability of TiO and TiN under high pressure

Carrier concentration effect and other structure-related parameters on lattice thermal conductivity of Si nanowires

Lecture 25: Heat and The 1st Law of Thermodynamics Prof. WAN, Xin

RESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS

Temperature and pressure dependence of the Raman frequency shifts in anthracene

Analysis of Mobility of Intrinsic Germanium and Silicon near Room Temperature

Intrinsic Small-Signal Equivalent Circuit of GaAs MESFET s

Potentials, periodicity

Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals

Transverse vibration and instability of fluid conveying triple-walled carbon nanotubes based on strain-inertia gradient theory

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

AH 2700A. Attenuator Pair Ratio for C vs Frequency. Option-E 50 Hz-20 khz Ultra-precision Capacitance/Loss Bridge

Introduction to solid state physics

Quantum Condensed Matter Physics Lecture 5

( x) f = where P and Q are polynomials.

Supplementary Figures

Life Prediction Under Multiaxial Fatigue

Chapter 4. Improved Relativity Theory (IRT)

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

General Case for Deriving Four Pole Coefficients for Gas Pulsation

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

Monte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions

Curve Sketching. The process of curve sketching can be performed in the following steps:

PARAMETRIC PHONON-HELICON INTERACTION IN CENTROSYMMETRIC CRYSTALS WITH STRAIN-DEPENDENT DIELECTRIC CONSTANT

Non-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings

Impact of size and temperature on thermal expansion of nanomaterials

Digest Journal of Nanomaterials and Biostructures Vol. 5, No 1, March 2010, p

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004

Speeds of sound and isothermal compressibility of ternary liquid systems: Application of Flory s statistical theory and hard sphere models

Pre-AP Physics Chapter 1 Notes Yockers JHS 2008

Sound Radiation Of Cast Iron

Stoneley Waves at the Boundary Surface of Modified Couple Stress Generalized Thermoelastic with Mass Diffusion

STUDY OF VIBRATIONAL DYNAMICS OF Ca-BASED METALLIC GLASSES BY PSEUDOPOTENTIAL THEORY

Earlier Lecture. Basics of Refrigeration/Liquefaction, coefficient of performance and importance of Carnot COP.

Reflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress

Transcription:

Vol. 109 006) ACTA PHYSICA POLONICA A No. Sound Attenuation at High Temperatures in Pt R.K. Singh and K.K. Pandey H.C.P.G. College, Varanasi-1001, U.P., India Received October 4, 005) Ultrasonic attenuation due to phonon phonon interaction and thermoelastic loss was evaluated in VIII group transition metal Pt in a wide temperature range 100 K to 1500 K) or longitudinal and shear waves along 100, 110 and 111 directions and or shear waves polarised along dierent directions. Electrostatic and Born Mayer potentials were used to obtain second and third order elastic constants, taking nearest-neighbour distance and hardness parameter as input data. Second and third order elastic constants data obtained at dierent temperatures were used to obtain Gruneisen parameters and non-linearity or anisotropy parameters which in turn were used to evaluate and in Akhiezer regime. The results were discussed and it was ound that at lower temperatures increases )l )s rapidly with temperature and then rate o increase becomes very small. Contribution to the total attenuation due to thermoelastic loss is negligible so that due to phonon phonon interaction establishing that major part o energy rom sound wave is removed due to interaction with thermal phonons lattice vibrations). PACS numbers: 6.0.Dc, 6.65.+k ) 1. Introduction Platinum Z = 78) is a silvery white transition metal o VIII group in VI period with an outer electronic coniguration o 5d 9 6s 1 and shows high catalytic activity. It is highly malleable and ductile metal with high melting point 041.4 K). Platinum is corrosion resistant and its wear and tarnish resistance characteristics are well suited or making ine jewelry. Other distinctive properties include resistance to chemical attack, excellent high temperature characteristics and sta- corresponding author; e-mail: rksingh 17@redimail.com. Present address: Physics Department, B.H.U., Varanasi-1005, India. 19)

0 R.K. Singh, K.K. Pandey ble electrical properties. All these properties have been exploited or industrial applications. In recent past, elastic [1], thermal and mechanical [] properties o Pt have been investigated extensively at room temperature as well as in wide temperature region up to 1500 K). So ar as sound characterisation o platinum is concerned, no result is available or this metal in wide temperature range. Ultrasonic attenuation studies have been made in dierent types o solids viz. dielectric [3], mixed dielectric [4], semiconducting [5, 6] and metallic [7, 8] substances at room temperature as well as at higher temperatures. But results are available up to 500 K. Due to high melting point o platinum and its important high temperature uses, ultrasonic attenuation study has been made or Pt rom 100 K to 1500 K. Dierent causes may be attributed to the attenuation o ultrasonic waves propagating through solids crystalline or amorphous). O these important actors or sound attenuation are electron phonon interaction, phonon phonon p p) interaction, thermoelastic loss and loss due to screw and edge dislocations. Phonon phonon interaction is the principal cause o ultrasonic attenuation at high temperatures viz. 100 K and above in all types o solids i.e. metallic, semiconducting and dielectric [4 8]. In metals a small contribution about to 5%) o total attenuation coeicient occurs due to thermoelastic loss i.e. due to heat transer rom compressed to rareied parts o metals due to strain caused by applied sound waves. Sound attenuation studies oer the possibility to detect and characterize microstructural properties as well as laws in materials, controlling material behaviours, based on physical mechanisms to predict uture perormance o the material. Structural inhomogeneities, elastic parameters, non-linearity parameters are well connected with requency or temperature dependence o ultrasonic attenuation mechanisms. Starting rom electrostatic and Born Mayer repulsive potentials and taking nearest-neighbour distance and Born repulsive parameters as input data and considering interactions up to nd nearest neighbours, second and third order elastic constants have been obtained at dierent temperatures which in turn have been used to obtain non-linearity coupling parameters and hence ultrasonic absorption coeicients. The present communication deals with temperature variation o ultrasonic attenuation due to p p interaction rom 100 K to 1500 K along 100, 110 and 111 crystallographic directions o propagation or longitudinal and shear waves and or shear waves polarized along dierent directions. Ultrasonic attenuation due to thermoelastic loss has been evaluated and it has been ound that its contribution relative to contribution made by p p interaction is about 5%. Temperature variation o ultrasonic attenuation has been ound to increase rapidly at lower temperatures and aterwards rate o increase becomes small which nearly saturates beyond 1000 K along all the directions o propagation or longitudinal and shear waves.

Sound Attenuation at High Temperatures in Pt 1. Theory.1. Theory o second and third order elastic constants Second and third order elastic constants SOEC and TOEC), C 0 ij and C0 ijk, at 0 K have been obtained using electrostatic and Born Mayer potentials and ollowing Brugger s [9] deinition o elastic constants. Repulsive parameter and nearest-neighbour distance have been used as input data and interaction up to next nearest-neighbours has been considered. According to Brugger s deinition, n-th order elastic constant is deined as C ijklmn... = u/ ε ij ε kl ε mn...), 1) where u is the crystal ree-energy density and ε ij is strain tensor. In Voigt notation C IJK replaces C ijklmn... in which ij I... etc. For cubic crystals Pt has cc structure) three independent SOEC C 11, C 1, and C 44 ) and six independent TOEC C 111, C 11, C 144, C 166, C 456, and C 13 ) occur. Using the theory discussed in [9], SOEC and TOEC at 0 K viz. Cij 0 and C0 ijk obtained are given as C 0 11 = 1.56933 z e r 4 0 C 0 1 = C 0 44 = 0.344778 z e + 1 1 + 1 ) Qr 1 ) + 1 + 1 ) Qr ), qr 0 r 0 q qr 0 r0 q r 4 0 + 1 1 + 1 ) Qr ), qr 0 r0 q C111 0 = 10.6390 z e 3 r0 4 qr0 + 3 q + 1 ) r 0 q 3 Qr 1 ) 3 r 0 q + 3 ) q + r 0 q 3 Qr ), C11 0 = C166 0 = 1.0965 z e 3 r0 4 + r0 q + 3 ) q + r 0 q 3 Qr ), C13 0 = C144 0 = C456 0 = 1.086 z e 3 r0 4 r0 q + 3 q + 1 ) Qr ), r 0 q 3 Qr 1 ) = A exp r ) 0, Qr ) = A exp r ) 0. q q According to lattice dynamics developed by Leibried and Ludwig [10] temperature variation o SOEC and TOEC have been obtained by adding vibrational contribution to elastic constants. SOEC and TOEC at any temperature are given by and C ij T ) = C 0 ij + ij )

R.K. Singh, K.K. Pandey C ijk T ) = Cijk 0 + Cijk vib, 3) where Cij vib and Cijk vib are vibrational contributions to elastic constants where where 11 = 1,1 G 1 + G, 1 = 1,1 G 1 + G 1,1, 44 = G 1,1, 111 = 1,1,1 G 3 1 + 3,1 G G 1 + 3 G 3, 11 = 1,1,1 G 3 1 +,1 G 1,1 + G )G 1 + 3 G,1, 13 = 1,1,1 G 3 1 + 3,1 G 1 G 1,1 + 3 G 1,1,1, 144 =,1 G 1 G 1,1 + 3 G 1,1,1, 456 = 3 G 1,1, 166 =,1 G 1 G 1,1 + 3 G,1, G 3 = Qr 1 ) 30 + 30q 0 + 9q0 q0 3 q0) 4 H + G.1, G,1 = Qr ) 15 + 15q 0 + 9 q q 0 3 ) q0 4 H, 0 G = Qr 0 ) 6 6q 0 q0 + q0) 3 H + G1.1, G 1,1 = Qr ) 3 6 q 0 ) q0 + q0 3 H, G 1 = Qr 0 ) + q 0 q 0) + Qr ) + q0 q 0 ) H, where Qr 1 ) = A exp = 3 = ηw 0 8r 3 0 r ) 0 q coth X, and Qr ) = A exp r ) 0, q 1,1 =,1 = 1 ηw 0 r0 3 48 [ ηw0 1,1,1 = ηw 0 384r 3 0 X = ηw 0 kt, H = kt η = h π X sinh + coth X, X ) coth X 6 sinh X + X ] sinh X + coth X, and k is Boltzmann constant, [ ) r0 r0 Qr 1 ) + q q ) 1 Qr )].

Sound Attenuation at High Temperatures in Pt 3.. Theory o sound attenuation In the Akhiezer regime a sound wave passing through a solid can be attenuated by two processes [11]. First, i the wave is longitudinal, periodic contractions and dilations in the solid induce a temperature wave via thermal expansion. Energy is dissipated by heat conduction between regions o dierent temperatures. This is called thermoelastic loss. Second, dissipation occurs as the gas o thermal phonons tries to reach an equilibrium characterized by a local sound wave induced) strain. This is an internal riction mechanism. The physical basis or obtaining attenuation coeicient is that the elastic constants contributed by the thermal phonons relax [1]. The phonon contribution to the unrelaxed elastic constants is evaluated by taking into consideration the change in energy o the thermal phonons due to applied instantaneous strain. The requency o each mode ν i is changed by ν i ν i = γ j i S j, where γ j i is generalised Gruneisen parameter and S j is instantaneous strain. It is assumed that all the phonons o a given direction o propagation and polarization have equal change in requency. Then phonons o i-th branch and j-th mode suer a change in temperature T i T 0 = γ j i S j T is the temperature). A relaxed elastic constant is obtained ater there is phonon phonon coupling among various branches and the T i relax to a common temperature change, T, given by T T = γj i S j, where γ j i is the average value o γj i. Thermal relaxation time is τ = τ s = τ l = 3K C V v, 4) where K is thermal conductivity, C V is speciic heat per unit volume and v is Debye average velocity. According to Mason and Batemann [1], SOEC and TOEC which are measure o anharmonicity o the crystal) are related by Gruneisen parameter γ j i and hence by non-linearity parameter, D. Ultrasonic attenuation due to phonon phonon interaction in Akhiezer regime ωτ 1) [1] is given by ) = π DE 0 τ 3v 3, 5) p p where non-linearity coupling constant D = 9 γ j i ) 3 γj i C V T. 6) E 0 γ j i ) and γ j i are square average and average square Gruneisen parameters, v is sound wave velocity or longitudinal waves v l ) and or shear waves v s ) and d is density. Debye average velocity is given by 3 v 3 = 1 vl 3 + vs 3. 7) Propagation o sound wave through crystal produces compression and rareactions as a result heat is transmitted rom compressed region at higher

4 R.K. Singh, K.K. Pandey temperature) to rareied region at lower temperature) and hence thermoelastic loss occurs, which is given by [13] ) = 4π γ j i KT dvl 5. 8) th 3. Results and discussions Using the theory discussed above values o SOEC and TOEC obtained at absolute zero temperature are given in Table I. Thermal energy density, E 0, and speciic heat per unit volume, C V are evaluated as unction o the Debye temperature, θ D, using physical constants table [14]. TABLE I C 0 ij and C 0 ijk in [10 11 dyne/cm ] or Pt. C11 0 C1 0 C111 0 C11 0 C144 0 19.14 6.57 308.04 6.6 10.63 Thermal relaxation time has been obtained using Eq. 4), taking thermal conductivity rom physical constants table [15], Gruneisen parameter γ j i, square average Gruneisen parameter γ j i ), and non-linearity constants D l and D s along dierent crystallographic directions are evaluated using Eq. 6) and D l and D s are given in Table II. Sound absorption coeicients or longitudinal wave, and or shear )l wave, are evaluated using Eq. 5). Temperature variation o )s and )l along dierent crystallographic directions o propagation and polarization )s are shown in Figs. 1 3. Temperature variation o along dierent directions )th or longitudinal wave has been presented in Table III. Values o SOEC and TOEC at 0 K are given in Table I. Cagin et al. [] have evaluated SOEC or Pt using computer simulation technique at dierent temperatures. Their values at 300 K are C 11 = 8.96 10 11 dyne/cm, C 1 = 3.93 10 11 dyne/cm, and C 44 = 6.50 10 11 dyne/cm, while our values at the same temperature are C 11 = 19.19 10 11 dyne/cm, C 1 = 19.14 10 11 dyne/cm C 44 = 6.51 10 11 dyne/cm. Thus there is good agreement between values o C 11 and C 44. Deviations in value o C 11 and slight deviation in value o C 44 can be attributed to the approach adopted by them. They have used NVE molecular dynamics simulation, in which they have not considered certain parameters. Also our approach is very simple, which involves only two basic parameters, nearest- -neighbour distance and hardness parameter.

Sound Attenuation at High Temperatures in Pt 5 TABLE II Temperature variation o D l and D s along dierent directions or Pt. 100 110 111 Temp. [K] D l D s D l Ds Ds D l D s + D s ++ 001 1 10 110 11 100 8.3 4.08 38.67 17.11 54.87 79.73 8.7.59 00 30.38 4.01 41.9 17.00 54.08 84.75 7.87.35 300 30.73 3.94 41.43 16.88 53.6 86.68 7.46.11 400 9.73 3.87 39.65 16.77 5.47 86.35 7.05 1.88 500 9.04 3.81 38.33 16.66 51.70 86.38 6.67 1.66 600 8.30 3.74 36.99 16.55 50.96 86.6 6.9 1.44 700 7.57 3.68 35.68 16.45 50.4 86.09 5.93 1.30 800 6.84 3.63 34.40 16.34 49.55 85.87 5.4 1.03 900 6.16 3.57 33.1 16.5 48.87 85.65 5.4 0.84 1000 5.48 3.5 3.03 16.16 48. 85.40 4.9 0.65 1100 4.05 3.47 30.96 16.07 47.59 85.18 4.60 0.48 100 4. 3.4 9.90 15.98 46.98 84.93 4.30 0.30 1300 3.65 3.38 8.93 15.89 46.39 84.71 4.00 0.14 1400 3.08 3.34 8.00 15.81 45.81 84.49 3.74 19.98 1500.54 3.9 7.11 15.73 4.5 84.8 3.44 19.8 Average and square average Gruneisen parameters have the values as expected [16]. Non-linearity constants ratio along all the three directions o propagation is as expected [16]. Fig. 1. Variation o ) 100 direction or platinum. with temperature or longitudinal and shear waves along

6 R.K. Singh, K.K. Pandey Fig.. Variation o ) 110 direction or platinum polarisations are shown). with temperature or longitudinal and shear waves along Fig. 3. Variation o ) 111 direction or platinum polarisations are shown). with temperature or longitudinal and shear waves along From Figs. 1 3 it can be seen that temperature variation o along all the directions or longitudinal and shear waves. and )l ) )s is similar increase with temperature. At lower temperatures rate o increase is large, as temperature increases beyond ) 500 K, rate o increase becomes small. Beyond 500 K rate o variation o with temperature is small.

Rate o increase in Sound Attenuation at High Temperatures in Pt 7 )l with temperature is little larger or longitudinal waves compared to shear waves. This may be attributed to larger value o thermal relaxation time, τ or longitudinal waves compared to that o shear waves Eqs. 4) and 5)). Parameter in [10 18 db s cm 1 ] at dierent temperature or Pt. Temp. [K] )th TABLE III Direction o propagation 100 110 111 100 0.98 0.55 1.65 00.5 1.6 3.63 300 3.36 1.85 5.13 400 4.1.36 6.5 500 5.00.78 7.08 600 5.67 3.15 7.69 700 6.5 3.47 8.13 800 6.75 3.74 8.41 900 7.18 3.97 8.57 1000 7.54 4.16 8.6 1100 7.83 4.3 8.60 100 8.07 4.45 8.51 1300 8.7 4.55 8.36 1400 8.44 4.6 8.19 1500 8.51 4.67 7.9 From Table III and rom the values o given in Figs. 1 3, it can be )th seen that along all the directions is negligible in comparison to )th or )l )s. Although experimental results are not available in the temperature range studied, yet on the basis o the above discussions it can be concluded that the present approach is correct. Acknowledgments One o us R.K. Singh) grateully acknowledges inancial support F. No. 6-70/SPT/00MRP/NRCB) rom University Grants Commission, Government o India.

8 R.K. Singh, K.K. Pandey Reerences [1] D. Verma, M.L. Verma, A. Verma, Ind. J. Phys. A 78, 337 004). [] T. Cagin, G. Dereli, M. Uludogan, M. Tomak, Phys. Rev. B 59, 3648 1999). [3] S.K. Kor, R.K. Singh, Acta Phys. Pol. A 80, 805 1991). [4] R.K. Singh, J.P. Acad. Sci. 5, 80 1996). [5] S.K. Kor, R.K. Singh, J. Acoust. Soc. India 1, 03 1993). [6] S.K. Kor, R.K. Singh, Acustica 79, 83 1993). [7] S.K. Kor, R.K. Singh, Acustica 78, 9 1993). [8] R.R. Yadav, O.K. Pandey, Acta Phys. Pol. A 107, 933 005). [9] K. Brugger, Phys. Rev. A 133, 1611 1964). [10] G. Leibried, W. Ludwig, in: Solid State Physics, Eds. F. Seitz, D. Turnbull, Academic Press Inc., New York 1961. [11] J. Fabian, P.B. Allen, Phys. Rev. Lett. 8, 1478 1999). [1] W.P. Mason, Physical Acoustics, Vol. III B, Academic Press, New York 1965. [13] R.R. Yadav, D. Singh, Acoust. Phys. 49, 595 003). [14] AIP Handbook, Academic Press, New York 1981. [15] CRC Data Series, CRC Press, Florida 1999 000. [16] S.K. Kor, R.K. Yadav, Kailash, J. Phys. Soc. Japan 55, 8 1986).