# Computations of Ultrasonic Parameters inZr100-XSnX Alloys

Authors: Pramod Kumar Yadawa
Journal: Physics Research International
 (2011)

---

## Abstract

The ultrasonic properties like ultrasonic attenuation, sound velocity in the hexagonalZr100-XSnX alloys have been studied along unique axis at room temperature. The second- and third-order elastic constants (SOEC & TOEC) have been calculated for these alloys using Lennard-Jones potential. The velocities VL and VS1 have minima and maxima, respectively, at 45° with unique axis of the crystal, while VS2 increases with the angle from unique axis. The inconsistent behaviour of angle-dependent velocities is associated to the action of second-order 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.


---

## Body

## 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 high-pressure experimental works indicate that, at room temperature and ambient pressure, Zr is a hexagonal close-packed (hcp) structure (α phase). At high temperature and zero pressure, it transforms martensitically into the body-centered 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 muffin-tin orbital (FP-LMTO) method [8], the tight-binding (TB) approach [9], the embedded-atom (EAM) method [10], and the modified embedded-atom (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 relationshipV=C/ρ, where C 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 higher-order 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 hexagonal-structuredZr100-XSnX alloys (Zr82Sn18: alloy 1; Zr80Sn20: alloy 2; Zr78Sn22: alloy 3; Zr76Sn24: alloy 4) at room temperature. The ultrasonic attenuation coefficient, acoustic coupling constants, higher-order 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 hexagonal-structured materials and alloys.


## 2. Theory

In the present investigation, the theory is divided into three parts.


### 2.1. Second- and Third-Order Elastic Constants

The elastic energy density (U) is function of the strain components.(1)U=F(exx;eyy;ezz;eyz;ezx;exy)=F(e1;e2;e3;e4;e5;e6),
where eij(iorj=x,y,z) is component of strain tensor.The second- (CIJ) and third- (CIJK) order elastic constants of material are defined by following expressions. (2)CIJ=∂2U∂eI∂eJ,IorJ=1,…,6,CIJK=∂3U∂eI∂eJ∂eK,IorJorK=1,…,6.The elastic energy density is well related to interaction potentialΦ(r) between atoms. Let the interaction potential be the Lennard-Jones Potential or many-body interaction potential, which is formulated as:(3)Φ(r)=-a0rm+b0rn,
where a0, b0 are constants, and m, n are integers. The definition of higher order elastic constants (2) with this  potential (3) under equilibrium and symmetric condition six second- and ten third-order elastic constants (SOEC and TOEC) for the hexagonal closed packed structured materials [13, 17] (4)C11=24.1p4C′,C12=5.918p4C′,C13=1.925p6C′,C33=3.464p8C′,C44=2.309p4C′,C66=9.851p4C′,C111=126.9p2B+8.853p4C′,C112=19.168p2B-1.61p4C′,C113=1.924p4B+1.155p6C′,C123=1.617p4B-1.155p6C′,C133=3.695p6B,C155=1.539p4B,C144=2.309p4B,C344=3.464p6B,C222=101.039p2B+9.007p4C′,C333=5.196p8B,
where p=c/a; axial ratio; C′=χa/p5; B=ψa3/p3. The rest second-order 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:(5)Φ=Φ0+χ∑i=12[Δri2]2+ψ∑i=12[Δri2]3.According to (5), χ and ψ can be written as:(6)χ=12![d2Φ(r)d(r2)2],ψ=13![d3Φ(r)d(r2)3].
In solving (6) for hexagonal closed packed structured materials:(7)χ=(18)[{nb0(n-m)}{an+4}],ψ=-χ{6a2(m+n+6)}.The parametersχ and ψ can be calculated using (7) with appropriate values of m, n, and b0 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 hexagonal-structured crystal [19, 20]. The ultrasonic velocities as a function of angle between direction of propagation and unique axis for hexagonal-structured materials are [21] (8)VL2={13(1V13+1V23+1V33)C33cos⁡2θ+C11sin2⁡θ+C44+{[C11sin2θ-C33cos⁡2θ+C44(cos⁡2θ-sin2θ)]2+4cos⁡2⁡θsin2θ(C13+C44)2}1/213(1V13+1V23+1V33)}/2ρ,VS12={13(1V13+1V23+1V33)C33cos⁡2θ+C11sin2⁡θ+C44-{[C11sin2⁡θ-C33cos⁡2θ+C44(cos⁡2θ-sin2⁡θ)]2+4cos⁡2θsin2⁡θ(C13+C44)2[C11sin2⁡θ-C33cos⁡2θ+C44(cos⁡2θ-sin2⁡θ)]2}1/2}/2ρVS22={C44cos⁡2⁡θ+C66sin2⁡θ}ρ,
where VL, VS1 and VS2 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 (TD) is an important physical parameter for the characterization of materials, which is well related to the Debye average velocity (VD): (9)TD=ℏVD(6π2na)1/3kB,
here
(10)VD={13(1V13+1V23+1V33)}-1/3,
where ℏ is quantum of action and is equal to Planck’s constant divided by 2π;kB is Boltzmann Constant; na 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 phonon-phonon interaction (Akhieser loss) and thermoelastic relaxation mechanisms. The ultrasonic attenuation coefficient(𝒜Akh) due to phonon-phonon interaction and thermoelastic relaxation mechanisms is given by the following expression [17, 22]:(11)(Af2)Akh=4π2(3E0〈(γij)2〉-〈γij〉2CVT)τ2ρV3,(12)(Af2)Th=4π2〈γij〉2kT2ρVL5,
where, f: frequency of the ultrasonic wave; V: ultrasonic velocity for longitudinal and shear wave; VL: longitudinal ultrasonic velocity; E0: thermal energy density; γij: Grüneisen number (i,j are the mode and direction of propagation).The Grüneisen number for hexagonal-structured crystal along〈001〉 orientation or θ=0° is direct consequence of second- and third-order elastic constants. D=3(3E0〈γij〉2-〈γij〉2CVT)/E0 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:(13)τ=τS=τL2=3kCVVD2.HereτL and τS are the thermal relaxation time for longitudinal and shear wave. k and CV are the thermal conductivity and specific heat per unit volume of the material, respectively.


## 2.1. Second- and Third-Order Elastic Constants

The elastic energy density (U) is function of the strain components.(1)U=F(exx;eyy;ezz;eyz;ezx;exy)=F(e1;e2;e3;e4;e5;e6),
where eij(iorj=x,y,z) is component of strain tensor.The second- (CIJ) and third- (CIJK) order elastic constants of material are defined by following expressions. (2)CIJ=∂2U∂eI∂eJ,IorJ=1,…,6,CIJK=∂3U∂eI∂eJ∂eK,IorJorK=1,…,6.The elastic energy density is well related to interaction potentialΦ(r) between atoms. Let the interaction potential be the Lennard-Jones Potential or many-body interaction potential, which is formulated as:(3)Φ(r)=-a0rm+b0rn,
where a0, b0 are constants, and m, n are integers. The definition of higher order elastic constants (2) with this  potential (3) under equilibrium and symmetric condition six second- and ten third-order elastic constants (SOEC and TOEC) for the hexagonal closed packed structured materials [13, 17] (4)C11=24.1p4C′,C12=5.918p4C′,C13=1.925p6C′,C33=3.464p8C′,C44=2.309p4C′,C66=9.851p4C′,C111=126.9p2B+8.853p4C′,C112=19.168p2B-1.61p4C′,C113=1.924p4B+1.155p6C′,C123=1.617p4B-1.155p6C′,C133=3.695p6B,C155=1.539p4B,C144=2.309p4B,C344=3.464p6B,C222=101.039p2B+9.007p4C′,C333=5.196p8B,
where p=c/a; axial ratio; C′=χa/p5; B=ψa3/p3. The rest second-order 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:(5)Φ=Φ0+χ∑i=12[Δri2]2+ψ∑i=12[Δri2]3.According to (5), χ and ψ can be written as:(6)χ=12![d2Φ(r)d(r2)2],ψ=13![d3Φ(r)d(r2)3].
In solving (6) for hexagonal closed packed structured materials:(7)χ=(18)[{nb0(n-m)}{an+4}],ψ=-χ{6a2(m+n+6)}.The parametersχ and ψ can be calculated using (7) with appropriate values of m, n, and b0 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 hexagonal-structured crystal [19, 20]. The ultrasonic velocities as a function of angle between direction of propagation and unique axis for hexagonal-structured materials are [21] (8)VL2={13(1V13+1V23+1V33)C33cos⁡2θ+C11sin2⁡θ+C44+{[C11sin2θ-C33cos⁡2θ+C44(cos⁡2θ-sin2θ)]2+4cos⁡2⁡θsin2θ(C13+C44)2}1/213(1V13+1V23+1V33)}/2ρ,VS12={13(1V13+1V23+1V33)C33cos⁡2θ+C11sin2⁡θ+C44-{[C11sin2⁡θ-C33cos⁡2θ+C44(cos⁡2θ-sin2⁡θ)]2+4cos⁡2θsin2⁡θ(C13+C44)2[C11sin2⁡θ-C33cos⁡2θ+C44(cos⁡2θ-sin2⁡θ)]2}1/2}/2ρVS22={C44cos⁡2⁡θ+C66sin2⁡θ}ρ,
where VL, VS1 and VS2 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 (TD) is an important physical parameter for the characterization of materials, which is well related to the Debye average velocity (VD): (9)TD=ℏVD(6π2na)1/3kB,
here
(10)VD={13(1V13+1V23+1V33)}-1/3,
where ℏ is quantum of action and is equal to Planck’s constant divided by 2π;kB is Boltzmann Constant; na 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 phonon-phonon interaction (Akhieser loss) and thermoelastic relaxation mechanisms. The ultrasonic attenuation coefficient(𝒜Akh) due to phonon-phonon interaction and thermoelastic relaxation mechanisms is given by the following expression [17, 22]:(11)(Af2)Akh=4π2(3E0〈(γij)2〉-〈γij〉2CVT)τ2ρV3,(12)(Af2)Th=4π2〈γij〉2kT2ρVL5,
where, f: frequency of the ultrasonic wave; V: ultrasonic velocity for longitudinal and shear wave; VL: longitudinal ultrasonic velocity; E0: thermal energy density; γij: Grüneisen number (i,j are the mode and direction of propagation).The Grüneisen number for hexagonal-structured crystal along〈001〉 orientation or θ=0° is direct consequence of second- and third-order elastic constants. D=3(3E0〈γij〉2-〈γij〉2CVT)/E0 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:(13)τ=τS=τL2=3kCVVD2.HereτL and τS are the thermal relaxation time for longitudinal and shear wave. k and CV are the thermal conductivity and specific heat per unit volume of the material, respectively.


## 3.  Results and Discussion

### 3.1. Higher-Order Elastic Constants

The unit cell parameters “a” (basal plane parameter) and “p” (axial ratio) for Zr100-XSnX alloys are 3.17 Å, 3.23 Å, 3.21 Å, 3.21 Å, and 1.622, 1.599, 1.603, 1.600, respectively, [23]. The values of m and n for chosen materials are 6 and 7. The value of b0 is 1.79×10-64 erg cm7 for all alloys. The SOEC and TOEC have been calculated for these alloys using (4) and are presented in Table 1.Table 1
SOEC, TOEC, and Bulk Modulus (B) in the unit of 1010 Nm−2 of ZrSn alloys at room temperature.
(a)AlloysC11C12C13C33C44C66B122.4165.5054.71222.3125.6528.79110.775219.1164.6943.90417.9624.6837.5009.003320.2154.9644.14719.1624.9747.9289.551420.4245.0154.17619.2405.0098.0099.628Zr [24]14.106.766.4316.692.583.68Sn [25]7.013.813.68(b)AlloysC111C112C113C123C133C344C144C155C222C3331−365.54−57.96−12.09−15.36−75.37−70.66−17.90−11.93−289.23−278.902−311.73−49.42−10.01−12.73−60.67−56.88−14.83−9.88−246.65−218.153−329.65−52.26−10.64−13.52−64.72−60.68−15.75−10.50−260.83−233.734−333.06−52.81−10.71−13.61−64.99−60.92−15.86−10.57−263.52−233.95Published 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 Table1, 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 (B) for these alloys can be calculated with the formula B=2(C11+C12+2C13+C33/2)/9. The evaluated B for these alloys is presented in Table 1.


### 3.2. Ultrasonic Velocity and Allied Parameters

The density (ρ) and thermal conductivity (k) at room temperature have been taken from the literature [23]. The value of CV and E0 are evaluated using tables of physical constants and Debye temperature. The quantities ρ,CV,E0,k, and calculated acoustic coupling constants (DL & DS) are presented in Table 2.Table 2
Density (ρ: in 103 kg m−3), specific heat per unit volume (CV: in 106 Jm−3K−1), thermal energy density (E0: in 108Jm−3), thermal conductivity (k: in Wm−1K−1), and acoustic coupling constant (DL, DS) of ZrSn alloys.
AlloysρCVE0kDLDS129.977.3315.125.4757.361.11228.617.0514.925.0856.411.19329.287.1415.044.5656.301.16429.537.1515.064.2156.371.18The calculated orientation-dependent ultrasonic wave velocities and Debye average velocities at 300 K are shown in Figures1, 2, 3 and 4. Figures 1–3 show that the velocities VL and VS1 have minima and maxima, respectively, at 45° with the unique axis of the crystal while VS2 increases with the angle from the unique axis. The combined effect of SOEC and density is reason for abnormal behaviour of angle-dependent velocities.Figure 1
VL versus angle with unique axis of crystal.Figure 2
VS1 versus angle with unique axis of crystal.Figure 3
VS2 versus angle with unique axis of crystal.Figure 4
VD versus angle with unique axis of crystal.The nature of the angle-dependent velocity curves in the present work is found similar as that for heavy rare earth metals, Laves-phase 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 (VD) of these alloys are increasing with the angle and have maxima at 55° at 300 K (Figure 4). Since VD is calculated using VL, VS1, and VS2 [13, 22]; therefore, the angle variation of VD is influenced by the constituent ultrasonic velocities. The maximum VD at 55° is due to a significant increase in longitudinal and pure shear (VS2) wave velocities and a decrease in quasishear (VS1) 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 Figure5. The angle-dependent thermal relaxation time curves follow the reciprocal nature of VD as τ∝3K/CVVD2. This implies that τ for these are mainly affected by the thermal conductivity. The τ for hexagonal-structured 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 θ=55° 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 quasi-identical, but in Figure 5, their relaxation times differ.Figure 5
Relaxation time versus angle with unique axis of crystal.


### 3.3. Ultrasonic Attenuation

In the evaluation of ultrasonic attenuation, it is supposed that wave is propagating along the unique axis (〈001〉 direction) of these alloys. The attenuation coefficient over frequency square (𝒜/f2)Akh for longitudinal (𝒜/f2)L and shear wave  (𝒜/f2)S are calculated using (11) under the condition ωτ≪1 at room temperature. The thermoelastic loss over frequency square (𝒜/f2)th is calculated with (12). The values of (𝒜/f2)L, (A/f2)S, (𝒜/f2)th, and total attenuation (𝒜/f2)Total, are presented in Table 3.Table 3
Ultrasonic attenuation coefficient (in 10−17 Nps2m−1) of ZrSn alloys.
Alloys1234(𝒜/f2)Th0.1150.5100.4040.373(𝒜/f2)L362.491529.937417.917385.668(𝒜/f2)s27.49742.02432.51330.404(𝒜/f2)Total390.103572.471450.834416.445Table3 indicates that the thermoelastic loss is very small in comparison to Akhieser loss, and ultrasonic attenuation for longitudinal wave (𝒜/f2)L is greater than that of shear wave (𝒜/f2)S. This reveals that ultrasonic attenuation due to phonon-phonon interaction along longitudinal wave is governing factor for total attenuation ((𝒜/f2)Total=(𝒜/f2)Th+(𝒜/f2)L+(𝒜/f2)S). 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𝒜∝V-3 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 third-group nitrides (AlN: 4.441×10-17 Nps2m−1, GaN: 14.930×10-17 Nps2m−1, and InN: 20.539×10-17 Nps2m−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 third-group nitrides.


## 3.1. Higher-Order Elastic Constants

The unit cell parameters “a” (basal plane parameter) and “p” (axial ratio) for Zr100-XSnX alloys are 3.17 Å, 3.23 Å, 3.21 Å, 3.21 Å, and 1.622, 1.599, 1.603, 1.600, respectively, [23]. The values of m and n for chosen materials are 6 and 7. The value of b0 is 1.79×10-64 erg cm7 for all alloys. The SOEC and TOEC have been calculated for these alloys using (4) and are presented in Table 1.Table 1
SOEC, TOEC, and Bulk Modulus (B) in the unit of 1010 Nm−2 of ZrSn alloys at room temperature.
(a)AlloysC11C12C13C33C44C66B122.4165.5054.71222.3125.6528.79110.775219.1164.6943.90417.9624.6837.5009.003320.2154.9644.14719.1624.9747.9289.551420.4245.0154.17619.2405.0098.0099.628Zr [24]14.106.766.4316.692.583.68Sn [25]7.013.813.68(b)AlloysC111C112C113C123C133C344C144C155C222C3331−365.54−57.96−12.09−15.36−75.37−70.66−17.90−11.93−289.23−278.902−311.73−49.42−10.01−12.73−60.67−56.88−14.83−9.88−246.65−218.153−329.65−52.26−10.64−13.52−64.72−60.68−15.75−10.50−260.83−233.734−333.06−52.81−10.71−13.61−64.99−60.92−15.86−10.57−263.52−233.95Published 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 Table1, 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 (B) for these alloys can be calculated with the formula B=2(C11+C12+2C13+C33/2)/9. The evaluated B for these alloys is presented in Table 1.


## 3.2. Ultrasonic Velocity and Allied Parameters

The density (ρ) and thermal conductivity (k) at room temperature have been taken from the literature [23]. The value of CV and E0 are evaluated using tables of physical constants and Debye temperature. The quantities ρ,CV,E0,k, and calculated acoustic coupling constants (DL & DS) are presented in Table 2.Table 2
Density (ρ: in 103 kg m−3), specific heat per unit volume (CV: in 106 Jm−3K−1), thermal energy density (E0: in 108Jm−3), thermal conductivity (k: in Wm−1K−1), and acoustic coupling constant (DL, DS) of ZrSn alloys.
AlloysρCVE0kDLDS129.977.3315.125.4757.361.11228.617.0514.925.0856.411.19329.287.1415.044.5656.301.16429.537.1515.064.2156.371.18The calculated orientation-dependent ultrasonic wave velocities and Debye average velocities at 300 K are shown in Figures1, 2, 3 and 4. Figures 1–3 show that the velocities VL and VS1 have minima and maxima, respectively, at 45° with the unique axis of the crystal while VS2 increases with the angle from the unique axis. The combined effect of SOEC and density is reason for abnormal behaviour of angle-dependent velocities.Figure 1
VL versus angle with unique axis of crystal.Figure 2
VS1 versus angle with unique axis of crystal.Figure 3
VS2 versus angle with unique axis of crystal.Figure 4
VD versus angle with unique axis of crystal.The nature of the angle-dependent velocity curves in the present work is found similar as that for heavy rare earth metals, Laves-phase 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 (VD) of these alloys are increasing with the angle and have maxima at 55° at 300 K (Figure 4). Since VD is calculated using VL, VS1, and VS2 [13, 22]; therefore, the angle variation of VD is influenced by the constituent ultrasonic velocities. The maximum VD at 55° is due to a significant increase in longitudinal and pure shear (VS2) wave velocities and a decrease in quasishear (VS1) 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 Figure5. The angle-dependent thermal relaxation time curves follow the reciprocal nature of VD as τ∝3K/CVVD2. This implies that τ for these are mainly affected by the thermal conductivity. The τ for hexagonal-structured 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 θ=55° 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 quasi-identical, but in Figure 5, their relaxation times differ.Figure 5
Relaxation time versus angle with unique axis of crystal.


## 3.3. Ultrasonic Attenuation

In the evaluation of ultrasonic attenuation, it is supposed that wave is propagating along the unique axis (〈001〉 direction) of these alloys. The attenuation coefficient over frequency square (𝒜/f2)Akh for longitudinal (𝒜/f2)L and shear wave  (𝒜/f2)S are calculated using (11) under the condition ωτ≪1 at room temperature. The thermoelastic loss over frequency square (𝒜/f2)th is calculated with (12). The values of (𝒜/f2)L, (A/f2)S, (𝒜/f2)th, and total attenuation (𝒜/f2)Total, are presented in Table 3.Table 3
Ultrasonic attenuation coefficient (in 10−17 Nps2m−1) of ZrSn alloys.
Alloys1234(𝒜/f2)Th0.1150.5100.4040.373(𝒜/f2)L362.491529.937417.917385.668(𝒜/f2)s27.49742.02432.51330.404(𝒜/f2)Total390.103572.471450.834416.445Table3 indicates that the thermoelastic loss is very small in comparison to Akhieser loss, and ultrasonic attenuation for longitudinal wave (𝒜/f2)L is greater than that of shear wave (𝒜/f2)S. This reveals that ultrasonic attenuation due to phonon-phonon interaction along longitudinal wave is governing factor for total attenuation ((𝒜/f2)Total=(𝒜/f2)Th+(𝒜/f2)L+(𝒜/f2)S). 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𝒜∝V-3 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 third-group nitrides (AlN: 4.441×10-17 Nps2m−1, GaN: 14.930×10-17 Nps2m−1, and InN: 20.539×10-17 Nps2m−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 third-group nitrides.


## 4. Conclusions

Based on the above discussion it is worthwhile to state the following:(i)
The present method to evaluate second- and third-order 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 velocitiesVL, VS1, and VS2, hence, a good resemblance in VD implies that our calculated velocities are correct.(iv)
The〈001〉 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θ=55° 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 third-group 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 phonon-phonon interaction mechanism is predominant over total attenuation as a governing factor of thermal conductivity.(viii)
Zr82Sn18 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 well-known physical properties, may expand future prospects for the application and study of these materials.



---

*Source: 350540-2011-07-03.xml*