Temperature Changes in the Elastic Properties of Sodium Chloride Single Crystals
Abstract Generated abstract
This paper investigates the low-temperature elastic behavior of sodium chloride single crystals as part of a systematic study of alkali-metal halides. Using a pulsed ultrasonic relative-measurement method, the authors measured longitudinal and transverse sound velocities in natural and artificially grown NaCl crystals from 300 K to 120 K, then calculated elastic constants, elastic compliances, anisotropy factors, compressibility, bulk modulus, Young’s modulus, shear modulus, and Poisson’s ratio. The results show increasing longitudinal velocities and elastic stiffness with decreasing temperature, comparatively stable transverse velocities, and better measurement reproducibility in artificially grown crystals than in natural rock salt. Extrapolated elastic constants at 0 K agree well with previously reported values by Overton and Swim.
Full Text
Temperature-Induced Changes in the Elastic Properties of Sodium Chloride Single Crystals
A. V. Sharko, A. A. Botaki
(Recommended for publication by the Department of General Physics of TPI)
The study of the elastic properties of crystals is of great importance for the development of the microscopic theory of the properties of solids. Despite the extensive work carried out by both theorists [1, 2] and experimentalists, the regularities in the change of mechanical properties with the chemical composition of a series of alkali-metal halides have as yet been little studied. The data available in the literature from experimental studies of temperature-induced changes in the elastic properties of single crystals of alkali-metal halides are of little use for revealing such dependences, since they were obtained by different authors and under different conditions (the only exception being the works of A. V. Stepanov and S. P. Nikanorov with co-workers [3—6], who measured temperature-induced changes in the elastic properties of most single crystals of alkali-metal halides at temperatures above room temperature).
The authors of the present work set themselves the aim of systematically studying temperature-induced changes in the elastic properties of alkali-metal-halide single crystals at temperatures below room temperature. The present work is devoted to the study of the elastic properties of sodium chloride single crystals by the pulsed ultrasonic method. The velocities of propagation of ultrasonic vibrations were measured with a DUK-6V pulsed ultrasonic flaw detector. Using the method of relative measurements, the following were determined: \(C_{l\infty}\), the velocity of propagation of longitudinal waves in an unbounded medium; \(C_l\), the velocity of propagation of longitudinal waves in a thin rod; and \(C_\perp\), the velocity of propagation of transverse waves. A medium is considered unbounded if its transverse dimensions are much greater than the wavelength of the ultrasonic vibrations, and a specimen is considered thin if its transverse dimensions are smaller than the wavelength of the elastic waves in the specimen.
The essence of the relative-measurement method is as follows. From the scale of the depth-measuring device of the flaw detector, calibrated using a reference specimen with a known ultrasonic propagation velocity \(C_0\), the apparent length \(l'\) of the specimen under test was determined; its true length \(l_0\) was determined by direct measurements. Denoting by \(C\) and \(C_0\), respectively, the velocities of propagation of elastic waves in the investigated and reference specimens, and by \(n\) the number of the multiple reflection of the echo pulse, the desired value of the velocity \(C\) can be determined from the relation
\[ l' = n \frac{C_0}{C} l_0 . \tag{1} \]
Specimens for measuring the elastic constants were made from natural single crystals of sodium chloride and from crystals grown artificially by the Kyropoulos method in the crystallization laboratory of TPI from chemically pure starting material. The specimens had the form of rectangular bars with faces coinciding with the cleavage planes of the single crystal. When recording the temperature dependence of the ultrasound propagation velocity, the specimen was glued to the probe head of the flaw detector with a mixture consisting of 88% ED-5 epoxy resin and 12% polyethylene polyamine and was placed in a cryostat. In preparing the specimens for measurements, special attention was paid to the quality and parallelism of the working faces of the specimen, as well as to their coincidence with the \(\{100\}\) planes. The working faces had a finely ground surface followed by additional polishing. In order to eliminate internal stresses and the effects of mechanical treatment, before the final finishing of the specimen surfaces the latter were annealed at a temperature of \(600^\circ\) for 24 hours, followed by slow cooling. The specimens for measurements of \(C_{l\infty}\) and \(C_{\perp}\) had dimensions \(80\times 40\times 40\) mm, and for measurements of \(C_l\), \(6\times 6\times 25\) mm. The density of sodium chloride single crystals for different temperatures [7] was determined from the formula
\[ \rho(t)=2.1680\left[1-1.12\cdot 10^{-4}t-5\cdot 10^{-8}t^2\right], \]
where \(t\) is the temperature in \({}^\circ\mathrm{C}\).
The change in the true length of the specimen \(l(t)\) as a function of temperature was calculated by the formula [8]
\[ l(t)=l_0\left[1+\int_t^{t_0}\alpha(t)\,dt\right], \]
where \(\alpha(t)\) is the coefficient of linear expansion, which, according to Henglein [9], can be determined from the formula
\[ \alpha(t)=\frac{\rho(t)-\rho_0}{\rho(t)\cdot t}, \tag{3} \]
where \(\rho(t)\) and \(\rho_0\) are the densities at temperatures \(t\) and \(0^\circ\mathrm{C}\), respectively. The Young’s modulus in the \(\langle 100\rangle\) direction and the elastic constant \(C_{11}\) were determined from the measured values of \(C_{l\infty}\) and \(C_l\) using the known formulas
\[ C_{l\infty}=\sqrt{\frac{C_{11}}{\rho}};\qquad C_l=\sqrt{\frac{E_{100}}{\rho}}. \tag{4} \]
To find the third independent constant—the shear modulus \(G\)—pulses of transverse waves were sent into the specimen; for this purpose prismatic probe heads at a frequency of \(2.5\) MHz with an ultrasonic-vibration entry angle of \(50^\circ\) were used, which made it possible to transform longitudinal waves into shear waves. In the general case, when an ultrasonic wave is incident on a plane-parallel plate at a certain angle to the normal, both longitudinal and transverse waves arise in the specimen. At the chosen angle of incidence of the ultrasound, total internal reflection of the longitudinal waves occurs, and only transverse waves propagate in the medium under investigation, directed at a certain angle with respect to the crystallographic axes. Taking this angle into account made it possible to determine the orientation factor \(F(\gamma)\) [10]
\[ F(\gamma)=\gamma_1^2\gamma_2^2+\gamma_1^2\gamma_3^2+\gamma_2^2\gamma_3^2, \tag{5} \]
where \(\gamma\) are the direction cosines of the angles formed by the direction of propagation of the transverse waves with the crystallographic axes.
Using the formula \(C_{\perp}=\sqrt{G/\rho}\), from the measured values of \(C_{\perp}\) we find \(G\) in the direction of propagation of the shear waves.
The calculation of the elastic constants \(C_{ik}\) and the elastic compliances \(S_{ik}\) was carried out using the well-known formulas [11]:
\[ S_{11}=\frac{1}{E_{100}}; \]
\[ C_{11}=\frac{S_{11}+S_{12}}{(S_{11}-S_{12})(S_{11}+2S_{12})}; \]
\[ \frac{1}{G}=S_{44}+4\left[(S_{11}-S_{12})-\frac{1}{2}S_{44}\right]F(\gamma); \]
\[ C_{44}=\frac{1}{S_{44}}; \tag{6} \]
\[ C_{12}=\frac{-S_{12}}{(S_{11}-S_{12})(S_{11}+2S_{12})}. \]
The Poisson coefficient was found from the expression [12]
Table 1
| T, °K | \(\rho\), g/cm³ | \(C_{l\infty}\) | \(C_l\) | \(C_{\perp}\) |
|---|---|---|---|---|
| 300 | 2,1613 | 4697 | 4410 | 2426 |
| 290 | 2,1639 | 4710 | 4427 | 2426 |
| 280 | 2,1663 | 4724 | 4447 | 2427 |
| 270 | 2,1687 | 4737 | 4471 | 2427 |
| 260 | 2,1711 | 4759 | 4492 | 2427 |
| 250 | 2,1732 | 4771 | 4512 | 2427 |
| 240 | 2,1759 | 4790 | 4551 | 2427 |
| 230 | 2,1782 | 4803 | 4556 | 2427 |
| 220 | 2,1806 | 4813 | 4574 | 2427 |
| 210 | 2,1829 | 4829 | 4594 | 2427 |
| 200 | 2,1851 | 4848 | 4615 | 2427 |
| 190 | 2,1874 | 4866 | 4635 | 2428 |
| 180 | 2,1896 | 4873 | 4662 | 2428 |
| 170 | 2,1919 | 4882 | 4673 | 2428 |
| 160 | 2,1943 | 4895 | 4693 | 2428 |
| 150 | 2,1962 | 4909 | 4714 | 2428 |
| 140 | 2,1984 | 4921 | 4731 | 2429 |
| 130 | 2,2005 | 4935 | 4751 | 2429 |
| 120 | 2,2025 | 4951 | 4772 | 2429 |
\[ \sigma=\frac{1-2\left(\dfrac{C_{\perp}}{C_{l\infty}}\right)^2}{2-2\left(\dfrac{C_{\perp}}{C_{l\infty}}\right)^2}. \tag{7} \]
The measurement data and the elastic constants and anisotropy factors calculated from them, \((S_{11}-S_{12})-\frac{1}{2}S_{44}\) and \(C_{44}; \frac{C_{11}-C_{12}}{2}\), as functions of temperature are given in Tables 1, 2, and 3. In addition, the tables give the temperature changes in the compressibility \(\chi=3(S_{11}+2S_{12})\), the bulk modulus \(\frac{C_{11}+2C_{12}}{3}\), and the measure of resistance to deformation caused by a shearing stress applied in the \(\{110\}\) plane in the \(\langle 100\rangle\) direction, \(\frac{C_{11}-C_{12}}{2}\); plots of these dependences are shown in Fig. 1. Extrapolation of the experimental values of the temperature dependences of the elastic constants to \(0^\circ\)K gives the following values:
Fig. 1.
\[ E_{100}=5.309\cdot10^{11}\ \frac{\mathrm{dyn}}{\mathrm{cm}^{2}}; \qquad G_{100}=C_{44}=1.310\cdot10^{11}\ \frac{\mathrm{dyn}}{\mathrm{cm}^{2}}; \]
\[ C_{11}=5.625\cdot10^{11}\ \frac{\mathrm{dyn}}{\mathrm{cm}^{2}}; \qquad C_{12}=0.973\cdot10^{11}\ \frac{\mathrm{dyn}}{\mathrm{cm}^{2}}; \]
\[ S_{11}=0.1874\cdot10^{-11}\ \frac{\mathrm{cm}^{2}}{\mathrm{dyn}}; \qquad S_{12}=0.0290\cdot10^{-11}\ \frac{\mathrm{cm}^{2}}{\mathrm{dyn}}, \]
\[ S_{44}=0.7368\cdot10^{-11}\ \frac{\mathrm{cm}^{2}}{\mathrm{dyn}}, \]
which is in good agreement with the values for \(0^\circ\)K obtained by Overton and Swim [8].
The determination of the elastic constants, as well as the recording of their low-temperature dependences, was carried out in parallel on natural and artificial—
Table 2
| \(T,\ ^\circ\mathrm{K}\) | Elastic constants \(\times 10^{-11}\ \mathrm{dyn}/\mathrm{cm}^{2}\): \(C_{11}\) | Elastic constants \(\times 10^{-11}\ \mathrm{dyn}/\mathrm{cm}^{2}\): \(C_{12}\) | Elastic constants \(\times 10^{-11}\ \mathrm{dyn}/\mathrm{cm}^{2}\): \(C_{44}\) | Elastic coefficients \(\times 10^{11}\ \mathrm{cm}^{2}/\mathrm{dyn}\): \(S_{11}\) | Elastic coefficients \(\times 10^{11}\ \mathrm{cm}^{2}/\mathrm{dyn}\): \(-S_{12}\) | Elastic coefficients \(\times 10^{11}\ \mathrm{cm}^{2}/\mathrm{dyn}\): \(S_{44}\) |
|---|---|---|---|---|---|---|
| 300 | 4.769 | 1.314 | 1.273 | 0.2380 | 0.0514 | 0.7855 |
| 290 | 4.801 | 1.306 | 1.274 | 0.2357 | 0.0504 | 0.7849 |
| 280 | 4.834 | 1.298 | 1.276 | 0.2334 | 0.0491 | 0.7837 |
| 270 | 4.868 | 1.291 | 1.277 | 0.2306 | 0.0484 | 0.7831 |
| 260 | 4.908 | 1.282 | 1.279 | 0.2284 | 0.0473 | 0.7819 |
| 250 | 4.948 | 1.271 | 1.280 | 0.2263 | 0.0464 | 0.7812 |
| 240 | 4.994 | 1.260 | 1.282 | 0.2238 | 0.0452 | 0.7800 |
| 230 | 5.025 | 1.252 | 1.283 | 0.2212 | 0.0443 | 0.7794 |
| 220 | 5.052 | 1.244 | 1.285 | 0.2192 | 0.0433 | 0.7782 |
| $T$, °K | $C_{11}$ | $C_{12}$ | $C_{44}$ | $S_{11}$ | $-S_{12}$ | $S_{44}$ |
|---|---|---|---|---|---|---|
| \multicolumn{3}{c}{Elastic constants $\times 10^{-11}$ dyn/cm²} | \multicolumn{3}{c}{Elastic coefficients $\times 10^{11}$ cm²/dyn} | |||||
| 210 | ,091 | 1,233 | 1,286 | 0,2169 | 0,0423 | 0,7776 |
| 200 | 5,134 | 1,218 | 1,287 | 0,2148 | 0,0414 | 0,7770 |
| 190 | 5,178 | 1,202 | 1,289 | 0,2127 | 0,0402 | 0,7756 |
| 180 | 5,201 | 1,191 | 1,290 | 0,2105 | 0,0393 | 0,7752 |
| 170 | 5,224 | 1,185 | 1,292 | 0,2088 | 0,0386 | 0,7740 |
| 160 | 5,258 | 1,169 | 1,294 | 0,2068 | 0,0376 | 0,7726 |
| 150 | 5,294 | 1,155 | 1,295 | 0,2049 | 0,0367 | 0,7722 |
| 140 | 5,323 | 1,145 | 1,297 | 0,2032 | 0,0360 | 0,7710 |
| 130 | 5,359 | 1,128 | 1,298 | 0,2012 | 0,0350 | 0,7702 |
| 120 | 5,399 | 1,111 | 1,299 | 0,1992 | 0,0340 | 0,7689 |
in artificially grown NaCl single crystals. The measurement results showed comparatively good agreement on both specimens, with a certain lowering of the elastic moduli in the rock-salt single crystal, which is apparently explained by the higher dislocation density in the latter. Repeated measurements gave better stability of the measurement results on the artificially grown crystal than on rock salt; therefore the table gives the results of measurements on the artificially grown single crystal.
Table 3
| $T$, °K | $\dfrac{C_{11}-C_{12}}{2}$ $\times 10^{-11}$ dyn/cm² |
$\dfrac{C_{44}}{C_{12}}$ | $\dfrac{C_{11}+2C_{12}}{3}$ $\times 10^{-11}$ dyn/cm² |
$\dfrac{2C_{44}}{C_{11}-C_{12}}$ | $\left[(S_{11}-S_{12})-\dfrac{1}{2}S_{44}\right]\times 10^{11}$ cm²/dyn |
$E_{100}\times 10^{-11}$ dyn/cm² |
$\times 10^{11}$ cm²/dyn |
$\sigma$ |
|---|---|---|---|---|---|---|---|---|
| 300 | 1,728 | 0,969 | 2,466 | 0,737 | 0,1033 | 4,201 | 0,4056 | 0,3181 |
| 290 | 1,747 | 0,976 | 2,471 | 0,729 | 0,1063 | 4,242 | 0,4047 | 0,3194 |
| 280 | 1,768 | 0,983 | 2,476 | 0,722 | 0,1090 | 4,285 | 0,4038 | 0,3207 |
| 270 | 1,783 | 0,989 | 2,484 | 0,714 | 0,1126 | 4,335 | 0,4026 | 0,3221 |
| 260 | 1,813 | 0,997 | 2,491 | 0,706 | 0,1152 | 4,380 | 0,4014 | 0,3242 |
| 250 | 1,838 | 1,007 | 2,497 | 0,696 | 0,1179 | 4,422 | 0,4005 | 0,3254 |
| 240 | 1,876 | 1,018 | 2,501 | 0,687 | 0,1210 | 4,467 | 0,3993 | 0,3273 |
| 230 | 1,886 | 1,024 | 2,510 | 0,680 | 0,1242 | 4,520 | 0,3984 | 0,3286 |
| 220 | 1,904 | 1,033 | 2,513 | 0,675 | 0,1266 | 4,561 | 0,3975 | 0,3295 |
| 210 | 1,929 | 1,043 | 2,519 | 0,667 | 0,1296 | 4,609 | 0,3969 | 0,3311 |
| 200 | 1,958 | 1,057 | 2,523 | 0,657 | 0,1323 | 4,655 | 0,3953 | 0,3327 |
| 190 | 1,998 | 1,073 | 2,527 | 0,648 | 0,1350 | 4,700 | 0,3957 | 0,3336 |
| 180 | 2,003 | 1,081 | 2,529 | 0,644 | 0,1378 | 4,750 | 0,3954 | 0,3348 |
| 170 | 2,020 | 1,091 | 2,531 | 0,640 | 0,1396 | 4,788 | 0,3951 | 0,3357 |
| 160 | 2,044 | 1,107 | 2,533 | 0,633 | 0,1420 | 4,835 | 0,3948 | 0,3369 |
| 150 | 2,069 | 1,121 | 2,535 | 0,626 | 0,1445 | 4,880 | 0,3945 | 0,3381 |
| 140 | 2,088 | 1,131 | 2,539 | 0,621 | 0,1463 | 4,922 | 0,3939 | 0,3389 |
| 130 | 2,115 | 1,150 | 2,529 | 0,614 | 0,1499 | 4,970 | 0,3939 | 0,3401 |
| 120 | 2,114 | 1,111 | 2,541 | 0,606 | 0,1517 | 5,020 | 0,2936 | 0,3409 |
References
- Huntington. Solid State Phys., 7, 213, 1958.
- K. A. Aleksandrov, T. V. Ryzheva. Kristallografiya, 6, 289, 1961.
- A. V. Stepanov and I. M. Eidus. ZhETF, 29, 669, 1955.
- S. P. Nikanorov, A. V. Stepanov. FTT, 6, 1987, 1964.
- S. P. Nikanorov, A. V. Stepanov. ZhETF, 37, 1814, 1959.
- Yu. M. Chernov, A. V. Stepanov. FTT, 3, 2872, 1961.
- International Critical Tables of Numerical Data (Mc—Graw—Hill Company N. V.), vol. II, p. 43, 1928.
- Overton Jr. and Swim R. T. Phys. Rev., 84, 758, 1951.
- Fr. A. Henglein. Zeit Physik Chemie, 115, 91, 1925.
- B. Boren. Defects in Crystals. Moscow, IL, p. 27, 1969.
- K. Zener. Elasticity and Anelasticity of Metals. Collection edited by S. V. Vonsovsky. Moscow, IL, p. 22, 1954.
- Nondestructive Testing (handbook edited by R. McMaster). Moscow, IL, vol. 2, p. 269, 1965.