# Elastic-stiffness coefficients of titanium diboride.

Using resonance ultrasound spectroscopy, we measured the monocrystal elastic-stiffness coefficients, the Voigt [C.sub.ij], of [TiB.sub.2]. With hexagonal symmetry, [TiB.sub.2] exhibits five independent [C.sub.ij]: [C.sub.11], [C.sub.33], [C.sub.44], [C.sub.12], [C.sub.13]. Using Voigt-Reuss-Hill averaging, we converted these monocrystal values to quasiisotropic (polycrystal) elastic stiffnesses. Briefly, we comment on effects of voids. From the [C.sub.ij], we calculated the Debye characteristic temperature, the Gruneisen parameter, and various sound velocities. Our study resolves the enormous differences between two previous reports of [TiB.sub.2]'s [C.sub.ij].Key words: Debye temperature; elastic constants; Gruneisen parameter; hexagonal crystal; monocrystal; sound velocities; titanium diboride; voids

Accepted: September 10, 2008

Available online: http://www.nist.gov/jres

1. Introduction

Currently, we see a renascence of research on transition-metal diborides. For example, well-known are many current studies on [ReB.sub.2] and [OsB.sub.2] as superhard materials, rivaling even diamond.

Titanium diboride's physical and mechanical properties received review elsewhere (1), (2). This compound is known well for low mass density, high hardness, high melting point, low electrical resistivity, good thermal conductivity, and good chemical inertness. However, the problems in preparing full-dense, high-quality monocrystals preclude their extensive study. Akimitsu and colleagues (3) reported superconductivity at 40 K in a same-crystal-structure companion compound: [MgB.sub.2].

The importance of elastic-stiffness coefficients for both science and technology also received review (4).

Titanium diboride's monocrystal elastic stiffnesses appeared in two reports. The values of Gilman and Roberts (5) depart strongly from the more recent report of Spoor and colleagues (6), mainly in the off-diagonal [C.sub.12] and [C.sub.13]. The first measurements were made by a pulse-echo method, the second by resonance ultrasonic spectroscopy, but by a very different mechanical setup than used in the present study.

Here, we report a third set of measurements, which agree closely with the Spoor et al. results. From our [C.sub.ij] measurements, we estimate the Debye characteristic temperature, the quintessential harmonic-lattice property. Also, we estimate the Gruneisen parameter, the quintessential anharmonic property.

Figure 1 shows the [TiB.sub.2] crystal structure, obviously a layered structure suggesting strong elastic anisotropy, not observed as shown below, creating a conundrum.

[FIGURE 1 OMITTED]

2. Measurements

2.1 Crystal

An oriented parallelepiped specimen 1.8 mm x 2.3 mm x 3.9 mm was prepared from a larger crystal grown by a floating zone method (7), (8). Fluorescent x-ray analysis revealed no significant impurities. The exact stoichiometry was [TiB.sub.1.97]. The crystal faces with respect to the above dimensions were [2 - 1 - 10], [- 1100], [0001]. These directions represent an a-axis, the c-axis, and the direction orthogonal to both. Laue x-ray diffraction confirmed these orientations within 1[degrees].

2.2 Method

To measure the [C.sub.ij], we used resonance ultrasound spectroscopy, summarized in Fig. 2 (9-11). Briefly, one clamps lightly a regular-shape (cube, cylinder, cube, parallelepiped) specimen between two piezoelectric transducers. One transducer is swept through frequency and the second transducer detects the specimen's macroscopic vibration frequencies (Fig. 3). Frequencies of a specimen are determined by five factors: (1) shape, (2) size, (3) mass or mass density, (4) elastic-stiffness coefficients [C.sub.ij], and (5) crystal-axis orientation relative to macroscopic shape. Thus, by measuring the resonance frequencies [f.sub.n], one can determine by an inverse calculation the stiffnesses [C.sub.ij]. The problem is strongly overdetermined: about one hundred [f.sub.n] to determine five [C.sub.ij]. Not reported here, one can also determine the complete internal-friction tensor [Q.sub.ij.sup.1] = [[DELTA][f.sub.n]/[f.sub.n]] = [[C*.sub.ij]/[C.sub.ij]], where [DELTA][f.sub.n] denotes resonance-peak width, [C*.sub.ij] the imaginary part, and [C.sub.ij] the real part of the total [C.sub.ij] tensor. Well-described elsewhere (12), the inverse problem involves Lagrangean minimization, the Rayleigh-Ritz method, and a least-squares procedure for measured and deduced [f.sub.n] values.

[FIGURE 2 OMITTED]

We determined mass density by careful mass and size measurements: [rho] = 4.502 [+ or -] 0.016 g/[cm.sup.3]. This compares with reported x-ray mass densities of 4.504 to 4.53. We ascribe the differences to heavier impurities on the Ti sites. From this mass density, we concluded our specimen contains no significant void content. From the sharp resonance peaks (Fig. 3), we concluded our specimen contains few cracks.

[FIGURE 3 OMITTED]

2.3 Errors

Errors arise from many sources: crystal orientation, crystal dimensions, nonparallelism, the inverse problem (measured-frequency to [C.sub.ijkl], conversion), specimen-transducer interactions. Several authors described these errors elsewhere, especially the first reference (13-15). The effect of the slight (1 %) departure from stoichiometry is hard to estimate. Because of the strong covalent bonding within the boron planes, we conjecture a very small error arising from an occasional missing boron atom.

3. Results and Discussion

Table 1 shows our principal results: the [C.sub.ij] and their uncertainties. (Because the two previous reports omitted error estimates, very detailed comparisons are precluded.) Only five [C.sub.ij] are independent because [C.sub.66] = [([C.sub.11] - [C.sub.12])/2]. As usual (they fail to correspond directly to a phonon), the off-diagonal [C.sub.ij] ([C.sub.12] and [C.sub.13]) show the largest uncertainties. Table 1 also gives the Spoor et al. results (6), which differ from ours by an average of 2.2 %. The average uncertainty in our six [C.sub.ij] is 0.7 %. The close agreement between our results and the Spoor et al. results suggests strongly that the earlier Gilman-Roberts results (5) are wrong, the largest discrepancies occurring in the off-diagonal [C.sub.ij]: [C.sub.12] and [C.sub.13]. Further support for the correctness of the Ledbetter-Tanaka/Spoor et al. results arises from ab initio calculations that yielded B = 251 GPa and v = 0.12 (16), versus 417 and 0.32 for Gilman-Roberts.

Table 1. Monocrystal and polycrystal elastic constants of titanium diboride. Unless specified, all units are GPa, except for the Poisson ratios, V. which arc dimensionless Present Gilman - Roberts Spoor et al. (b) [rho](g/[cm.sup.3]) 4.502 [+ or -] 0.016 Monocrystal elastic constants [C.sub.11] 654.4 [+ or -] 1.9 690 660 [C.sub.33] 458.1 [+ or -] 1.4 440 432 [C.sub.44] 262.6 [+ or -] 0.3 250 260 [C.sub.66] 302.7 [+ or -] 0.7 140 306 [C.sub.12] 48.98 [+ or -] 1.4 410 48 [C.sub.13] 95.25 [+ or -] 0.55 320 93 [E.sub.11] 633.3 [+ or -] 2.2 389 639 [E.sub.33] 432.3 [+ or -] 1.7 254 408 [V.sub.12] 0.0460 [+ or -] 0.0022 0.3877 0.0437 [V.sub.13] 0.1984 [+ or -] 0.0015 0.4453 0.2059 [V.sub.31] 0.1354 [+ or -] 0.0010 0.2909 0.1314 Debye characteristic temperature [[THETA].sub.D] (K) 1217 [+ or -] 6 989 1211 Voigt-Reuss-Hill-average quasiisotropic (polycrystal) elastic constants [C.sub.L] 599.9 642 593 B 247.5 417 244 G 264.3 169 262 E 584.7 446 579 V 0.1063 0.3219 0.1037 (a) Ref. 5. (b) Ref. 6.

Table 1 also shows the principal Young moduli [E.sub.ii] computed by

[E.sub.ij] = [S.sub.ij.sup.-1]. (1)

Here [S.sub.ij] denotes the elastic-compliance tensor, the tensor inverse of the [C.sub.ij]. As we expect from [C.sub.11] > [C.sub.33], [E.sub.11] exceeds [E.sub.33]; that is, [TiB.sub.2] is much stiffer within the basal plane than along the c-axis. Obviously, this relates to the crystal structure where the covalent-bonded boron atoms lie in the plane perpendicular to [x.sub.3].

Table 1 shows also the three principal Poisson ratios [V.sub.ij] computed by

[v.sub.ij] = -[[S.sub.ij] / [S.sub.ii]]. (2)

Within the boron plane, the Poisson ratio [v.sub.12] is extremely low, reflecting strong covalent bonding and the strong resistance of boron atoms to change their bond angles. The [v.sub.13] Poisson ratio is only slightly below normal, indicating weaker bonds out of the boron-atom planes than those within the planes.

The shear elastic anisotropy of hexagonal crystals can be expressed in various ways. The simplest is [[C.sub.66]/[C.sub.44]], 1.15 for [TiB.sub.2], thus weak elastic anisotropy (the isotropic case corresponding to 1.00). Because the Young modulus depends so strongly on the shear modulus, a Young-modulus variation with direction also serves as an effective shear-anisotropy indicator. Spoor and colleagues showed a Young-modulus representation surface; it is nearly spherical (6), From the alternating-layer boron-titanium crystal structure (Fig. 1), one expects higher elastic anisotropy than one finds.

The lower part of Table 1 gives the averaged-over-direction quasiisotropic elastic constants obtained from the [C.sub.ij] by a Voigt-Reuss-Hill average (17). These constants include longitudinal modulus [C.sub.L], bulk modulus B, shear modulus G, Young modulus E, and Poisson ratio v. These are the elastic constants appropriate to a full-density nontextured polycrystalline aggregate in which the grain boundaries cause no softening. Grain size produces no effect on elastic constants if it is small relative to specimen size. (A reviewer pointed out that nanosize grains may soften the elastic constants.) The most unusual feature of our results is the high GIB ratio, thus low Poisson ratio. Our nonpublished results on a [YB.sub.66] monocrystal gave an averaged-over-direction (Kroner method) Poisson ratio of 0.13. For boron, a handbook value is v= 0.089 (18). Both the bulk and shear moduli of [TiB.sub.2] exceed considerably (by 30 % to 40 %) the handbook values of boron: 248/179 = 1.39 and 264/203 = 1.30, showing the strong interatomic bonds in [TiB.sub.2], both for extension-compression and for shear-bending. Beside implying a low Poisson ratio, the very high G/B ratio holds implications for several other crystal-bonding properties: covalency, Cauchy-relation departure, many-body forces, and others. We plan to discuss all these elsewhere.

The elastic stiffnesses yield three useful sound velocities: the longitudinal velocity

[v.sub.l] = [([[C.sub.L] / [rho]]).sup.[1/2]] = [[[[B + (4/3)G]/[rho]]].sup.[1/2]], (3)

the shear or transverse velocity

[v.sub.s] = [(G/[rho]).sup.[1/2]], (4)

and the mean velocity (as defined by Debye)

[3/[v.sub.m.sup.3]] = [2/[v.sub.t.sup.3]] + [1/[v.sub.l.sup.3]]. (5)

We found [v.sub.1], = 1.54, [v.sub.s] = 0.766, and [v.sub.m] = 0.835 cm/[mu]s.

From the [C.sub.ij] and the mass density, we can compute the acoustic Debye characterisitic temperature [[THETA].sub.D]. At zero temperature, the acoustic [[THETA].sub.D] becomes identical with the calorimetric [[THETA].sub.D] (19). [[THETA].sub.D] is proportional to the mean sound velocity:

[[THETA].sub.D] = [Kv.sub.m], (6)

where [v.sub.m] denotes mean sound velocity and K is given by

K = [h/k][(3/[4[pi][V.sub.a]]).sup.[1/3]]. (7)

Here h denotes Planck's constant, k Boltzmann's constant, and [V.sub.a] atomic volume. The velocity [v.sub.m] comes from the integration over all directions:

[3v.sub.m.sup.-3] = [integral] [summation over ([alpha] = 1,3)] [v.sub.[alpha].sup.-3] [[d[OMEGA]]/[4[pi]]]. (8)

Here [v.sub.1] denotes the quasilongitudinal wave velocity, [v.sub.2] and [v.sub.3] the quasitransverse wave velocities, and d[OMEGA] an increment of solid angle. Equation (8) can not be integrated analytically, and numerous numerical and approximate methods have been used for its solution. Phase velocities [v.sub.a] are roots of the Christoffel equations:

det([C.sub.ijkl][x.sub.j][x.sub.k] - [rho][v.sup.2][[delta].sub.il]) = 0. (sum over repeated indices) (9)

This expression follows from equations of motion for plane, monochromatic waves, where [rho] denotes mass density, [C.sub.ijkl] fourth-rank elastic-stiffness tensor, [x.sub.i] components of the unit wave vector relative to the cubic axes, and [[delta].sub.il] the Kronecker operator. Equation (9) usually yields three distinct real roots [rho][V.sup.2].

We computed an exact [[THETA].sub.D] by using a distribution of 489 vectors proposed by Overton and Schuch (20), which, we distributed over the usual 48 [100]-[l 10]-[111] stereographic unit triangles, thus a total of 23 472 directions. We obtained [[THETA].sub.D] = 1217K. This result differs enormously from the handbook calorimetric value, 1576 K. Often, calorimetric values contain large errors because of large extrapolations to zero temperature and/or large uncertainties in the lowest-temperature lattice specific heat.

From B and from handbook values of heat capacity C, volume thermal expansivity [beta], and volume V, we computed the effective Gruneisen parameter:

[gamma] = [[[B.sub.S][beta]V]/[C.sub.P]]. (10)

We obtained [gamma] = 1.71, the handbook values for boron and titanium being 1.85 and 1.33. Alternative gammas can be computed when one knows the third-order elastic-stiffness coefficients, the [C.sub.ijklmn]] (21). The few gammas known for compounds similar to [TiB.sub.2] preclude any comparisons. Following Pearson's reasoning (22), the Ti-B bonds would lead to [gamma] = 2.

Finally, because some authors suspect voids/cracks in [TiB.sub.2], we want to describe briefly how the above full-dense quasiisotropic elastic constants would change with voids. Focusing on [A1.sub.2][O.sub.3], Ledbetter, Lei, and Datta (23) gave a theory for void effects on elastic constants. Principal results include the following: Voids soften the bulk modulus more than the shear modulus. In the dilute limit, for spherical voids, our results agree with the classical results of Mackenzie (24):

[[[DELTA]B]/B] = -3c, ([v.sub.0] = [1/3]) (11)

[[[DELTA]G]/G] = -[15/8]c. ([v.sub.0] = [1/3]) (12)

Here, c denotes void volume fraction and [v.sub.0] the void-free-state Poisson ratio. However, we emphasize that the often-used rule-of-thumb that elastic stiffness varies as mass density is true only in the dilute limit. Void shape plays a key role, especially if the voids possess an oblate-spheroid (disc) shape. Dunn and Ledbetter (25) focused on the interesting, unexpected effects of voids and cracks on the Poisson ratio. Other authors addressed this problem using various approximations (26).

4. Conclusions

1. Our elastic-stiffness-coefficient measurements on titanium diboride support the report of Spoor and colleagues rather than the older (now handbook) values of Gilman and Roberts. The principal differences in the two previous [C.sub.ij] sets lie in [C.sub.12] and [C.sub.13].

2. From the [C.sub.ij] we computed several additional useful physical properties:

(i) sound velocities; (ii) Debye temperature;

(iii) Gruneisen parameter.

3. Our computed acoustic Debye temperature, 1217 K, is about 20 % lower than the handbook calorimetric value.

4. Our computed Gruneisen parameter, 1.71, suggests the importance of Ti-B bonds along with B-B bonds.

5. Using the [[C.sub.66]/[C.sub.44]] ratio as a shear-mode elasticanisotropy criterion, [TiB.sub.2] shows low shearmode anisotropy, 1.15. A surprise because of the B-Ti-B ... layer crystal structure.

6. Our elastic constants estimated for a full-dense polycrystal depart strongly from those given in Munro's review (2), especially in the uncertainties. For example, Munro proposed a 70 % uncertainty in the Poisson ratio and a 24 % uncertainty in the bulk modulus. Our polycrystal results differ from the Spoor et al. result by only 2 %. This finding agrees with the well-known fact that modern measurement methods give the elastic-stiffness coefficients easily within one percent.

7. Finally, we are surprised by the low elastic shear anisotropy shown by such a strongly layered crystal structure. Reference (8) gives a possible explanation. Part of the boron-atom p-obitals lie in a p-d hybridized band, weakening the [rho]-electron contribution to B-B bonding. If the bonding weakens, the anisotropy decreases.

Acknowledgments

Thanks especially to Dr. Edwin Fuller (NIST-Gaithersburg) for an incisive critical reading. S. Kim (NIST-Boulder) helped with measurements. Most of this study proceeded while H. Ledbetter was at NIST-Boulder. F. Drymiotis (Clemson U., Physics) contributed Fig. 1.

5. References

(1) J. Hoard and R. Huges, Chemistry of Boron and its Compounds (Wiley, New York, 1967).

(2) R. Munro, Material properties of titanium diboride. J. Res. Natl. Inst. Stand. Technol. 105, 709-720 (2000).

(3) J. Nagamatsu, N. Kakagawa, T. Muranaka, Y. Zenitani, and J. Akimitsu, Superdonductivity at 39 K in magnesium diboride, Nature 410, 63-64(2001).

(4) H. Ledbetter, Elastic properties, in Materials at Low Temperatures (Amer. Soc. Metals. Metals Park, Ohio, 1983) pp. 1-45.

(5) J. Gilman and B. Roberts, Elastic constants of TiC and Ti[B.sub.2], J. Appl. Phys. 32, 1405(1961).

(6) P. Spoor, J. Maynard, M. Pan, D. Green, J. Hellmann, and T. Tanaka, Elastic constants and crystal anisotropy of titanium diboride, Appl. Phys. Lett. 70, 1959-1961 (1997).

(7) Y. Ishizawa and T. Tanaka, in Science of Hard Materials, E. A. Almond, C. A. Brooks, and R. Warren, eds. (Adam Hilger, Boston 1986) pp. 29-43.

(8) T. Tanaka and Y. Ishizawa, The deHaas-van Alphen effect in [TiB.sub.2], J. Phys. C 13, 6671-6676 (1980).

(9) O. Anderson, Rectangular parallelepiped resonance--A technique of resonance ultrasound and its applications to the determination of elasticity at high temperatures, J. Acoust. Soc. Amer. 91, 2245-2253 (1982).

(10) A. Migliori and J. Sarrao, Resonant Ultrasound Spectroscopy (Wiley-Interscience, New York, 1997).

(11) J. Maynard, Resonant ultrasound spectroscopy, Phys. Today 49 (1), 26-31 (1996).

(12) P. Heyliger, A. Jilani, H. Ledbetter, R. Leisure, and C.-L. Wang, Elastic constants of isotropic cylinders using resonant ultrasound, J. Acoust. Soc. Amer. 94, 1482-1487 (1993).

(13) P. Spoor, P. White, and J. Maynard, An investigation of computational problems associated with resonant ultrasound spectroscopy, J. Acoust. Soc. Am. 99, 2492-2603 (1996).

(14) R. Leisure and F. Willis, Resonant ultrasound spectroscopy, J. Phys.: Condens. Matter 9, 6001-6029 (1997).

(15) R. Schwarz and J. Vuorinen, Resonant ultrasound spectroscopy: applications, current status, and limitations, J. Alloys and Compounds 310, 243-250 (2000).

(16) K. Panda and K. Chandran, Determination of elastic constants of titanium diboride ([TiB.sub.2]) from first principles using FLAPW implementation of the density-functional theory, Computational Mater. Sci. 35, 134-150 (2006).

(17) H. Ledbetter, Monocrystal-polycrystal elastic-constant models, in Handbook of Elastic Properties of Solids, Liquids, and Gases, Volume III (Academic, New York 2000) pp. 313-324.

(18) K. Gschneidner, Physical properties and interrelationship of metallic and semimetallic elements, Solid State Phys. 16, 275-426(1964).

(19) G. Leibfried and W. Ludwig, Theory of anharmonic effects in crystals, Solid State Phys. 12, 275-444 (1961).

(20) W. Overton and A. Schuch, Method for highly accurate numerical spatial integration: Applications to lattice dynamics of cubic-type crystals, Los Alamos Lab. Rept. 3162 (1965).

(21) H. Ledbetter and S. Kim, Elastic Gruneisen parameters of cubic elements and compounds, in Handbook of Elastic Properties of Solids, Liquids, and Gases, Volume II (Academic, New York, 2001) pp. 107-122.

(22) R. Pearson, Chemical Hardness (Wiley, New York, 1992) p. 188.

(23) H. Ledbetter, M. Lei, and S. Datta, Elastic constants of porous ceramics, in Handbook of Elastic Properties of Solids, Liquids, and Gases, Volume II (Academic, New York, 2001) pp. 463-471.

(24) J. Mackenzie. Elastic constants of a solid containing spherical holes, Proc. Phys. Soc. Lond. B63, 2-11 (1950).

(25) M. Dunn and H. Ledbetter, Poisson ratio of porous and micro-cracked solids, J. Mater. Res. 10, 2715-2722 (1995).

(26) M. Kachanov, On the effective elastic constants of cracked solids, Int. J. Fract. 146, 295-399 (2007).

Hassel Ledbetter

Mechanical Engineering

University of Colorado

Boulder, Colorado 80309

Takaho Tanaka

National Institute for Materials Science,

Sengen, Namiki, Tsukuba,

Ibaraki 305-0047

Japan

hassel.ledbetter@colorado.edu

tanaka.takaho@nims.go.jp

About the authors: Hassel Ledbetter is a materials scientist in Mechanical Engineering at the University of Colorado in Boulder, Colorado. He retired from the Materials Reliability Division at NIST-Boulder, where he established the world's most complete elastic-coefficient-measurement laboratory, especially for measurements to liquid-helium temperature. Then he spent four years at Los Alamos National Laboratory, at the National High Magnetic Field Laboratory, studying low-temperature physical properties of various materials, including actinides such as plutonium. Takaho Tanaka is a materials scientist at the National Institute for Materials Science (NIMS). He retired from the Boride Research Group at NIMS, where he established the world's highest-temperature crystal-growth laboratory, especially for growing crystals of refractory materials such as rare-earth borides, transition metal borides, and carbides. He also discovered and characterized many icosahedral boron-rich rare-earth borides.

The National Institute of Standards and Technology is an agency of the U.S. Department of Commerce.

Printer friendly Cite/link Email Feedback | |

Author: | Ledbetter, Hassel; Tanaka, Takaho |
---|---|

Publication: | Journal of Research of the National Institute of Standards and Technology |

Article Type: | Report |

Geographic Code: | 9JAPA |

Date: | Nov 1, 2009 |

Words: | 3431 |

Previous Article: | Comparison of the NIST and PTB air-kerma standards for low-energy x-rays. |

Next Article: | DAVE: a comprehensive software suite for the reduction, visualization, and analysis of low energy neutron spectroscopic data. |

Topics: |