Research Article  Open Access
Pramod Kumar Yadawa, "Computations of Ultrasonic Parameters in Alloys", Physics Research International, vol. 2011, Article ID 350540, 7 pages, 2011. https://doi.org/10.1155/2011/350540
Computations of Ultrasonic Parameters in Alloys
Abstract
The ultrasonic properties like ultrasonic attenuation, sound velocity in the hexagonal alloys have been studied along unique axis at room temperature. The second and thirdorder elastic constants (SOEC & TOEC) have been calculated for these alloys using LennardJones potential. The velocities and have minima and maxima, respectively, at 45Β° with unique axis of the crystal, while increases with the angle from unique axis. The inconsistent behaviour of angledependent velocities is associated to the action of secondorder elastic constants. Debye average sound velocities of these alloys are increasing with the angle and has maximum at 55Β° with unique axis at room temperature. Hence, when a sound wave travels at 55Β° with unique axis of these alloys, then the average sound velocity is found to be maximum. The mechanical and ultrasonic properties of these alloys will be better than pure Zr and Sn due to their high SOEC and ultrasonic velocity and low ultrasonic attenuation. The comparison of calculated ultrasonic parameters with available theoretical/experimental physical parameters gives information about classification of these alloys.
1. Introduction
Group IV transition metal zirconium (Zr) and its alloys are very important materials both from scientific and technological points of view. Scientifically, the electronic transfer between the broad sp band and the narrow d band is the driving force behind many structural and electronic transitions in these materials [1β3]. Technologically, these alloys have applications in the aerospace industry due to their light weight, static strength, and stiffness; they do not degrade rapidly as the temperature increases, and they also show oxidation resistance [4]. The mechanical properties of these alloys can be greatly improved by controlling the crystallographic phases present.
Static highpressure experimental works indicate that, at room temperature and ambient pressure, Zr is a hexagonal closepacked (hcp) structure ( phase). At high temperature and zero pressure, it transforms martensitically into the bodycentered cubic (bcc) structure ( phase) before reaching the melting temperature [5] while at room temperature and under pressure, the hcp phase transforms into another hexagonal structure.
The first form of tin is called grey tin (Sn) and is stable at temperature below 286.2βK. This form crystallizes in the diamond structure [6], in common with the lighter elements in column IVβA of the periodic table. At temperature above 286.2βK, grey tin slowly turns into tinβs second form, white tin. This second form (Sn) is a tetragonal distortion of diamond with two atoms per unit cell [7]. Tin is used to form many useful alloys.
Up to now, only a few theoretical methods have been applied successfully to calculate the elastic constants of Zr, such as the full potential linear muffintin orbital (FPLMTO) method [8], the tightbinding (TB) approach [9], the embeddedatom (EAM) method [10], and the modified embeddedatom (MEAM) method [11, 12]. However, elastic constants of ZrSn alloys have not yet been reported.
Ultrasonic properties offer the possibility to detect and characterize microstructural properties, as well as flaws in materials, controlling materials behaviour based on physical mechanism to predict future performance of the materials. Various investigators have shown considerable interest on ultrasonic properties of metals and alloys. Wave propagation velocity is key parameter in ultrasonic characterization and can provide information about crystallographic texture. The ultrasonic velocity is directly related to the elastic constants by the relationship , where is the relevant elastic constants, and is the density of that particular material. Also ultrasonic attenuation is very important physical parameter to characterize the material, which is well related to several physical quantities like thermal conductivity, specific heat, thermal energy density and higherorder elastic constants [13, 14]. The elastic constants provide valuable information about the bonding characteristic between adjacent atomic planes and the anisotropic character of the bonding and structural stability [15, 16].
Therefore, in this work I predict the ultrasonic properties of hexagonalstructured alloys (: alloy 1; : alloy 2; : alloy 3; : alloy 4) at room temperature. The ultrasonic attenuation coefficient, acoustic coupling constants, higherorder elastic constants, thermal relaxation time, and ultrasonic wave velocities for these alloys for each direction of propagation of wave are calculated at 300βK. The calculated ultrasonic parameters are discussed with related thermophysical properties for the characterization of the chosen alloys. The obtained results are analyzed in comparison to other hexagonalstructured materials and alloys.
2. Theory
In the present investigation, the theory is divided into three parts.
2.1. Second and ThirdOrder Elastic Constants
The elastic energy density () is function of the strain components. where is component of strain tensor.
The second () and third () order elastic constants of material are defined by following expressions.
The elastic energy density is well related to interaction potential between atoms. Let the interaction potential be the LennardJones Potential or manybody interaction potential, which is formulated as: where , are constants, and , are integers. The definition of higher order elastic constants (2) with this potential (3) under equilibrium and symmetric condition six second and ten thirdorder elastic constants (SOEC and TOEC) for the hexagonal closed packed structured materials [13, 17] where ; axial ratio; ; . The rest secondorder elastic constants have zero value because under 180Β° rotation, they have equal and opposite value for the same stress. The harmonic and anharmonic parameters ( and ) can be calculated using one experimental SOEC [13, 17]. In the present study, I have expanded the theory for theoretical evaluation of parameters and . The potential energy can be expanded in the powers of changes in the squares of distances. The expansion up to cubic term can be written as:
According to (5), and can be written as: In solving (6) for hexagonal closed packed structured materials:
The parameters and can be calculated using (7) with appropriate values of , , and so that the calculated values of elastic constants justify the experimental data.
2.2. Ultrasonic Velocity
The anisotropic behaviour of the material can be understood with the knowledge of ultrasonic velocity because the velocity is related to the second order elastic constants [17]. On the basis of mode of atomic vibration, there are three types of velocities (longitudinal, quasishear, and shear) in acoustical region [18]. These velocities vary with the direction of propagation of wave from the unique axis of hexagonalstructured crystal [19, 20]. The ultrasonic velocities as a function of angle between direction of propagation and unique axis for hexagonalstructured materials are [21] where , and are longitudinal, quasi shear and pure shear wave ultrasonic velocities. Variables and represent the density of the material and angle with the unique axis of the crystal, respectively. The Debye temperature () is an important physical parameter for the characterization of materials, which is well related to the Debye average velocity (): here where is quantum of action and is equal to Planckβs constant divided by ; is Boltzmann Constant; is atom concentration.
The above formulae have been used for the evaluation of ultrasonic velocity and related parameters for the selected materials.
2.3. Ultrasonic Attenuation and Allied Parameters
The predominant causes for the ultrasonic attenuation in a solid at room temperature are phononphonon interaction (Akhieser loss) and thermoelastic relaxation mechanisms. The ultrasonic attenuation coefficient due to phononphonon interaction and thermoelastic relaxation mechanisms is given by the following expression [17, 22]: where, : frequency of the ultrasonic wave; : ultrasonic velocity for longitudinal and shear wave; : longitudinal ultrasonic velocity; : thermal energy density; : GrΓΌneisen number ( are the mode and direction of propagation).
The GrΓΌneisen number for hexagonalstructured crystal along orientation or is direct consequence of second and thirdorder elastic constants. is known as acoustic coupling constant, which is the measure of acoustic energy converted to thermal energy. When the ultrasonic wave propagates through crystalline material, the equilibrium of phonon distribution is disturbed. The time for reestablishment of equilibrium of the thermal phonon distribution is called thermal relaxation time () and is given by following expression:
Here and are the thermal relaxation time for longitudinal and shear wave. and are the thermal conductivity and specific heat per unit volume of the material, respectively.
3. Results and Discussion
3.1. HigherOrder Elastic Constants
The unit cell parameters βaβ (basal plane parameter) and βpβ (axial ratio) for alloys are 3.17ββ«, 3.23ββ«, 3.21ββ«, 3.21ββ«, and 1.622, 1.599, 1.603, 1.600, respectively, [23]. The values of and for chosen materials are 6 and 7. The value of is βergβcm^{7} for all alloys. The SOEC and TOEC have been calculated for these alloys using (4) and are presented in Table 1.
(a)  
 
(b)  

Published theoretical/experimental data on the elastic constants for the chosen alloys are not available for comparison. Comparison for the present values to theoretical/experimental studies of pure metals, as shown in Table 1, demonstrates that results obtained from the present investigation are higher than those of the pure metals as reported by other investigators [24, 25]. The obtained values of the SOEC and TOEC are of the same order as previous experimental and theoretical studies of other metallic alloys and metals [19, 21, 26β28]. Therefore, it is concluded that the theoretical approach to evaluate the SOEC and TOEC is valid for the chosen alloys. The bulk modulus () for these alloys can be calculated with the formula . The evaluated for these alloys is presented in Table 1.
3.2. Ultrasonic Velocity and Allied Parameters
The density () and thermal conductivity () at room temperature have been taken from the literature [23]. The value of and are evaluated using tables of physical constants and Debye temperature. The quantities , and calculated acoustic coupling constants ( & ) are presented in Table 2.

The calculated orientationdependent ultrasonic wave velocities and Debye average velocities at 300βK are shown in Figures 1, 2, 3 and 4. Figures 1β3 show that the velocities and have minima and maxima, respectively, at 45Β° with the unique axis of the crystal while increases with the angle from the unique axis. The combined effect of SOEC and density is reason for abnormal behaviour of angledependent velocities.
The nature of the angledependent velocity curves in the present work is found similar as that for heavy rare earth metals, Lavesphase compounds, and other hexagonal wurtzite structured materials (GaN, AlN, InN) [13, 17, 19, 21, 22, 28]. Thus, the computed velocities for these alloys are justified.
Debye average velocities () of these alloys are increasing with the angle and have maxima at 55Β° at 300βK (Figure 4). Since is calculated using , , and [13, 22]; therefore, the angle variation of is influenced by the constituent ultrasonic velocities. The maximum at 55Β° is due to a significant increase in longitudinal and pure shear () wave velocities and a decrease in quasishear () wave velocity. Thus, it can be concluded that when a sound wave travels at 55Β° with the unique axis of these alloys then the average sound wave velocity is maximum.
Thus, the present average sound velocity directly correlates with the Debye temperature, specific heat, and thermal energy density of theses materials.
The calculated thermal relaxation time is visualised in Figure 5. The angledependent thermal relaxation time curves follow the reciprocal nature of as . This implies that for these are mainly affected by the thermal conductivity. The for hexagonalstructured materials is the order of picoseconds [19]. Hence, the calculated justifies the hcp structure of chosen alloys at room temperature. The minimum for wave propagation along implies that the reestablishment time for the equilibrium distribution of thermal phonons will be minimum for propagation of wave along this direction. The value of thermal conductivity of alloy: 3 is greater than alloy: 4. So, in Figures 1β4, the ultrasonic velocities of alloy 3 and alloy 4 seem quasiidentical, but in Figure 5, their relaxation times differ.
3.3. Ultrasonic Attenuation
In the evaluation of ultrasonic attenuation, it is supposed that wave is propagating along the unique axis ( direction) of these alloys. The attenuation coefficient over frequency square for longitudinal and shear wave are calculated using (11) under the condition at room temperature. The thermoelastic loss over frequency square is calculated with (12). The values of , , , and total attenuation , are presented in Table 3.

Table 3 indicates that the thermoelastic loss is very small in comparison to Akhieser loss, and ultrasonic attenuation for longitudinal wave is greater than that of shear wave . This reveals that ultrasonic attenuation due to phononphonon interaction along longitudinal wave is governing factor for total attenuation . The total attenuation is mainly affected by thermal energy density and thermal conductivity. Thus, it may predict that at 300βK the alloy 1 behaves as its purest form and is more ductile as evinced by minimum attenuation while alloy 2 is least ductile. Therefore impurity will be least in the alloy 1 at room temperature.
Since and velocity is the largest for alloy 1 among alloy 2; thus, the attenuation should be smallest and material should be most ductile. The minimum ultrasonic attenuation for alloy 1 justifies its quite stable hcp structure state. The total attenuation of these alloys are much larger than thirdgroup nitrides (AlN: βNps^{2}m^{β1}, GaN: βNps^{2}m^{β1}, and InN: βNps^{2}m^{β1}) due to their large thermal conductivity and acoustic coupling constants [19, 29]. This implies that the interaction between acoustical phonon and quanta of lattice vibration for these alloys is large in comparison to thirdgroup nitrides.
4. Conclusions
Based on the above discussion it is worthwhile to state the following:(i)The present method to evaluate second and thirdorder elastic constants involving many body interaction potential for hexagonal wurtzite crystal structured materials (alloys) is correct.(ii)All elastic constants and density are mainly the affecting factor for anomalous behaviour of ultrasonic velocity in these alloys.(iii)When a sound wave travels at 55Β° with the unique axis of these crystals, then the average sound wave velocity is maximum. Since the Debye average velocity is calculated using the constituent velocities , , and , hence, a good resemblance in implies that our calculated velocities are correct. (iv)The direction is the direction of symmetry for these alloys as they have the same quasishear and pure shear wave velocities.(v)The order of thermal relaxation time for these alloys is found in picoseconds, which justifies their hcp structure at 300βK. The reestablishment time for the equilibrium distribution of thermal phonons will be minimum for the wave propagation along due to being the smallest value of along this direction.(vi)The acoustic coupling constant of these alloys for longitudinal wave are found five times larger than thirdgroup nitrides. Hence, the conversion of acoustic energy into thermal energy will be large for these alloys. This shows general suitability of chosen alloys.(vii)The ultrasonic attenuation due to phononphonon interaction mechanism is predominant over total attenuation as a governing factor of thermal conductivity.(viii)Zr_{82}Sn_{18} alloy is more suitable for industrial and other uses, as it has the highest elastic constants as well as wave velocity in comparison to other chosen alloys.(ix)The mechanical and ultrasonic properties of these alloys will be better than pure Zr and Sn due to their high SOEC and ultrasonic velocity and low ultrasonic attenuation.
Thus, obtained results in the present work can be used for further investigations, general and industrial applications. Our theoretical approach is valid for ultrasonic characterization of these materials at room temperature. The ultrasonic behavior in these alloys as discussed above shows important microstructural characteristic feature, which is well connected to thermoelectric properties of the materials. These results, together with other wellknown physical properties, may expand future prospects for the application and study of these materials.
Acknowledgments
The author wishes to express sincere thanks to Dr. Ashok K. Chauhan, President, AMITY School of Engineering and Technology (ASET), New Delhi, Prof. B. P. Singh, Senior Director, ASET New Delhi, and Prof. S. K. Kor and Prof. Raja Ram Yadav, Department of Physics, Allahabad University, for encouraging and providing necessary facilities in carrying out the work.
References
 J. S. Gyanchandani, S. C. Gupta, S. K. Sikka, and R. Chidambaram, βStructural stability of hafnium under pressure,β Journal of Physics, vol. 2, no. 30, Article ID 011, pp. 6457β6459, 1990. View at: Publisher Site  Google Scholar
 J. S. Gyanchandani, S. C. Gupta, S. K. Sikka, and R. Chidambaram, βThe equation of state and structural stability of titanium obtained using the linear muffintin orbital bandstructure method,β Journal of Physics, vol. 2, no. 2, pp. 301β305, 1990. View at: Publisher Site  Google Scholar
 A. L. Kutepov and S. G. Kutepova, βCrystal structures of Ti under high pressure: theory,β Physical Review B, vol. 67, no. 13, article 132102, 2003. View at: Google Scholar
 R. Ahuja, L. Dubrovinsky, N. Dubrovinskaia et al., βTitanium metal at high pressure: synchrotron experiments and ab initio calculations,β Physical Review B, vol. 69, no. 18, article 184102, 2004. View at: Publisher Site  Google Scholar
 D. Young, Phase Diagrams of the Elements, University of California Press, Berkeley, Calif, USA, 1991.
 P. Villars and L. D. Calvert, Eds., Pearson's Handbook of Crystallographic Data for Intermetallic Phases, American Society for Metals, Materials Park, Ohio, USA, 2nd edition, 1991.
 G. Simmons and H. Wang, Single Crystal Elastic Constants and Calculated Aggregate Properties: A Handbook, MIT Press, Cambridge, Mass, USA, 2nd edition, 1971.
 E. Clouet, J. M. Sanchez, and C. Sigli, βFirstprinciples study of the solubility of Zr in Al,β Physical Review B, vol. 65, no. 9, article 094105, 2002. View at: Google Scholar
 I. Schnell, M. D. Jones, S. P. Rudin, and R. C. Albers, βTightbinding calculations of the elastic constants and phonons of hcp Zr: complications due to anisotropic stress and longrange forces,β Physical Review B, vol. 74, no. 5, article 054104, 2006. View at: Publisher Site  Google Scholar
 D. J. Oh and R. A. Johnson, βRelationship between c/a ratio and point defect properties in HCP metals,β Journal of Nuclear Materials, vol. 169, pp. 5β8, 1989. View at: Google Scholar
 Y. M. Kim, B. J. Lee, and M. I. Baskes, βModified embeddedatom method interatomic potentials for Ti and Zr,β Physical Review B, vol. 74, no. 1, article 014101, 2006. View at: Google Scholar
 M. I. Baskes and R. A. Johnson, βModified embedded atom potentials for HCP metals,β Modelling and Simulation in Materials Science and Engineering, vol. 2, no. 1, pp. 147β163, 1994. View at: Publisher Site  Google Scholar
 P. K. Yadawa, D. Singh, D. K. Pandey, and R. R. Yadav, βElastic and acoustic properties of heavy rareearth metals,β The Open Acoustics Journal, vol. 2, no. 1, pp. 61β67, 2009. View at: Google Scholar
 D. Singh and P. K. Yadawa, βEffect of platinum addition to coinage metals on their ultrasonic properties,β Platinum Metals Review, vol. 54, no. 3, pp. 172β179, 2010. View at: Publisher Site  Google Scholar
 P. Ravindran, L. Fast, P. A. Korzhavyi, B. Johansson, J. Wills, and O. Eriksson, βDensity functional theory for calculation of elastic properties of orthorhombic crystals: application to TiSi_{2},β Journal of Applied Physics, vol. 84, no. 9, pp. 4891β4904, 1998. View at: Google Scholar
 L. Louail, D. Maouche, A. Roumili, and F. A. Sahraoui, βCalculation of elastic constants of 4d transition metals,β Materials Letters, vol. 58, no. 24, pp. 2975β2978, 2004. View at: Publisher Site  Google Scholar
 A. K. Yadav, R. R. Yadav, D. K. Pandey, and D. Singh, βUltrasonic study of fission products precipitated in the nuclear fuel,β Materials Letters, vol. 62, no. 1718, pp. 3258β3261, 2008. View at: Publisher Site  Google Scholar
 C. Oligschleger, R. O. Jones, S. M. Reimann, and H. R. Schober, βModel interatomic potential for simulations in selenium,β Physical Review B, vol. 53, no. 10, pp. 6165β6173, 1996. View at: Google Scholar
 D. K. Pandey, D. Singh, and R. R. Yadav, βUltrasonic wave propagation in IIIrd group nitrides,β Applied Acoustics, vol. 68, no. 7, pp. 766β777, 2007. View at: Publisher Site  Google Scholar
 S. K. Verma, R. R. Yadav, A. K. Yadav, and B. Joshi, βNondestructive evaluations of gallium nitride nanowires,β Materials Letters, vol. 64, no. 15, article 11226, 2010. View at: Publisher Site  Google Scholar
 D. K. Pandey, P. K. Yadawa, and R. R. Yadav, βAcoustic wave propagation in Lavesphase compounds,β Materials Letters, vol. 61, no. 25, pp. 4747β4751, 2007. View at: Publisher Site  Google Scholar
 D. K. Pandey, P. K. Yadawa, and R. R. Yadav, βUltrasonic properties of hexagonal ZnS at nanoscale,β Materials Letters, vol. 61, no. 30, pp. 5194β5198, 2007. View at: Publisher Site  Google Scholar
 D. Biswas, A. K. Meikap, S. K. Chattopadhyay, and S. K. Chatterjee, βLowtemperature electrical transport properties of disordered Zr_{100x}Sn_{x} alloys,β Physica B, vol. 382, no. 12, pp. 51β57, 2006. View at: Publisher Site  Google Scholar
 Y.J. Hao, L. Zhang, X.R. Chen, Y.H. Li, and H.L. He, βPhase transition and elastic constants of zirconium from firstprinciples calculations,β Journal of Physics, vol. 20, no. 23, article 235230, 2008. View at: Publisher Site  Google Scholar
 A. Berroukche, B. Soudini, and K. Amara, βMolecular dynamics simulation study of structural, elastic and thermodynamic properties of tin below 286K,β International Journal of Nanoelectronics and Materials, vol. 1, no. 1, pp. 41β51, 2008. View at: Google Scholar
 B. MagyariKope, L. Vitos, and G. Grimvall, βPerturbative quantumstate estimation anomalous behaviour of lattice parameters and elastic constants in hcp AgZn alloys,β Physical Review B, vol. 77, article 052102, 2004. View at: Google Scholar
 Y. Matsuo, βElastic behavior and phase stability in the H.C.P. $\epsilon $phase of AgZn alloy,β Journal of the Physical Society of Japan, vol. 53, no. 4, pp. 1360β1365, 1984. View at: Google Scholar
 D. K. Pandey, D. Singh, and P. K. Yadawa, βUltrasonic study of osmium and ruthenium,β Platinum Metals Review, vol. 53, no. 2, pp. 91β97, 2009. View at: Publisher Site  Google Scholar
 H. Landolt and R. Bornstein, LandoltBornstein: Numerical Data and Function Relationship in Science and Technology. Group III, vol. 11, Springer, Berlin, Germany, 1979.
Copyright
Copyright © 2011 Pramod Kumar Yadawa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.