Published in physic.philica.com
In this paper we discuss and compare theoretical and experimental data of thermal conductivity in crystalline Silicon (natural and enriched). We focus in particular on the role that the isotopic mass defects have in determining the value of this conductivity.
Thermal Conductivity of the Crystalline Silicon
Amelia Carolina Sparavigna
Department of Applied Science and Technology, Politecnico di Torino, Torino, Italy
Abstract: In this paper we discuss and compare theoretical and experimental data of thermal conductivity in crystalline Silicon (natural and enriched). We focus in particular on the role that the isotopic mass defects have in determining the value of this conductivity.
Keywords: Solid State physics, Materials Science, Thermal Conductivity, Crystalline Silicon, Isotopic Mass Defects, Phonon Scattering.
In the theoretical study of the thermal conductivity in non-metallic solids, a widely used approach is that of solving the Boltzmann equation of the phonon thermal transport to obtain the related conductivity [1-3]. Many researchers are using methods of solution based on the relaxation time approximation, of which we proposed a simplified discussion in [4,5]. Other researchers prefer the use of methods which are involving more specific atomic models [6-12]. In general, all the used methods are able of describing the behavior of the thermal conductivity, when they are involving a proper formulation of the normal and umklapp phonon-phonon scattering and of the other scattering mechanisms due to the presence of defects in the crystal.
If we consider a solid where only phonons are acting in the thermal transport, the sample turns out to be like a box containing them. At very low temperatures, we have a few phonons interacting only with the walls of the box. As the temperature increases, the population of phonons increases too and consequently their thermal conductivity. However, with a further increase of temperature, also the interactions among phonons become relevant. Therefore, after the thermal conductivity has reached a maximum, its value starts decreasing with a further increase of the temperature.
Among the processes which are influencing the values of the thermal conductivity we find the scattering of phonons due to the presence of isotopic mass defects. Actually, the natural crystals of several elements are composed by different isotopes; that is, in these crystals with natural composition we have the sites of the crystalline lattice occupied by atoms of the same element having different masses. This disorder of the lattice reduces the thermal conductivity.
For an experimental study of the thermal conductivity without the isotope effect, the sample needs to be made of an isotopically enriched material. Among the materials which have received specific enrichment for the analysis of thermal conductivity we find diamond-like Carbon, Silicon and Germanium [13-17]. As resulting from experiments, the isotopic purity of the solid is strongly influencing the thermal conductivity. This is important for technological applications of these materials, such as for a better comprehension of the thermal transport.
Here we discuss the case of crystalline Silicon. Specifically, we discuss some of the literature on it, which is providing experimental data. Theoretical works had been proposed in the past by the author, who solved the phononic Boltzmann equation in crystals in which different isotopes are present. Let us start for these theoretical works on the thermal conductivity of Silicon.
An iterative approach
In [18,19] we proposed an iterative approach to the solution of the Boltzmann equation. In the calculation of the thermal conductivity, we inserted competing scattering mechanism due to isotope impurities and the three-phonon collisions, without resorting to the systematic use of the relaxation time approximation. The model that we used was not so specific as in the following works [6,7], but it was able to give a good agreement with experimental data for silicon, here shown in the Figure 1.
Figure 1: Thermal conductivity of Silicon as a function of temperature at room pressure. The curves A and B are the results of calculations proposed in , with and without isotopic mass defects, in the case of absence of the boundary scattering. This scattering is considered for curve B’. The experimental points are from Holland and Neuringer .
Let us note that the results given in the Figure 1 had not been obtained by using adjustable parameters. Actually, the use of adjustable parameters is rather common in the methods involving the relaxation time approximation. A criticism to this use was made by Srivastava too .
In  we have also considered the isotope effect in Germanium. And for it too we observed the relevant role of the purity of the sample. For our further discussion on Silicon, let us extrapolate from the Figure 1 the thermal conductivity at 200 K. For the natural Silicon it is of 2.2 W/cm-K and for the pure sample, it is of 3.8 W/cm-K.
Following the methods given in [6,7], the effect of isotopes on the thermal conductivity of Silicon was investigated again in . Actually, in ,  and , we applied the iterative method on a revised model of the solid: we used a crystal described by its real Brillouin zone and Hamiltonian to determine properly the three-phonon normal and umklapp processes.
Figure 2: Thermal conductivity of Silicon as a function of temperature at room pressure. The curves A and B are the results of calculations proposed in , with and without isotope scattering. The experimental points are from Refs. 20 and 23.
Again, for the further discussion, let us extrapolate from the Figure 2 the thermal conductivity at 200 K. For the natural Silicon it is of 2. W/cm-K and for the pure sample, it is of 3. W/cm-K .
As we can see from the values obtained in  and , we have that the ratio of the thermal conductivity of enriched and natural Silicon, at 200 K, is 1.72 in the case of , and it is 1.5 in the case of . It means that, using a more realistic mode, the effect of the isotopic mass defects is reduced. From , we estimate that the thermal conductivity at 200 K is reduced of 30% at the most.
At the time we investigated the use of the iterative approach to solve the Boltzmann equation, we had not the opportunity of having a service, such as Google Scholar or other web supports, to find easily the literature on the subject. As the reader can see, we had at our disposal only the access to a few experimental data for Silicon. Now, let us examine the further data that we find on the web today.
First, we consider the paper by Capinski, Maris, Bauser, et al. . The experimental data are given in the Figure 3.
Fig. 3: Thermal conductivity of isotopically pure 28Si (squares) and Si of natural isotopic abundance (triangles) . The circles are the data from Ref.25.
The data from Ref.24 tell us that, at 200 K, the thermal conductivity of the natural Silicon is of about 2.8 W/cm-K and that of the isotopically pure Silicon of about 4 W/cm-K . The ratio of the thermal conductivity is 1.42. The isotope scattering is reducing the thermal conductivity of 30%. These results seem in good agreement with the previously given theoretical data.
In 2000, new data on the thermal conductivity of Silicon had been published . Here in the Figure 4, we give some of the experimental points.
Fig. 4: Thermal conductivity of isotopically pure 28Si and Si of natural isotopic abundance. Some data from Ref.26. The open dots represent the measurements on sample SIN1 of natural silicon. As told in  the “plus” symbols denote the “standard” curve for natural Si, as given by .
Let us consider the data from Ref.26 at 200 K. The thermal conductivity of the natural Silicon is of about 2.3 W/cm-K and that of the isotopically pure Silicon of about 4 W/cm-K (Figure 3). In , the authors discussed the result in the following manner. “In principle, the thermal conductivity enhancement, which is most prominent at the maximum of kappa, where umklapp processes are frozen out, should vanish rapidly at higher temperatures. However, we find an almost constant enhancement in kappa of about 60% between 100 and 300 K and estimate this effect to persist up to at least 400 K. This observation also contradicts the predictions of a recent calculation for the thermal conductivity enhancement in isotopically pure silicon”. The authors are referring to article .
The behavior at high temperature
Data published in  had been discussed in the work by Inyushkin in . He noted, about the phonon free path at low temperatures, that a “tremendous discrepancy” between theoretical and experimental values exists, and it “cannot be rationalized in the framework of the current views of heat conduction in crystals in the boundary scattering regime”. Actually, also the lack of results for sample SIN1 (natural Silicon) at low temperatures was stressed.
The authors of  considered Inyushkin’s analysis and reconsider their results. “In view of the interest, in this problem” , “we have recently remeasured the same and other similar samples with a different apparatus around T = 300 K. We consistently found an increase in the k of 28Si with respect of that of natSi of about 10 %, whereas the values of this enhancement measured below 100 K are the same as those in the original data” .
Fig. 5: The same as in the Figure 4 with the red squares representing the correction reported in  for the data of pure 28Si.
The authors, in the new measurements, found the thermal conductivity of the natural Silicon the same as before, but the “values reported earlier for 28Si are systematically 90 W/m-K higher than the present ones”. “This difference … below 100 K … reduces the enhancement … reported earlier for 28Si from 60% to about 10%. The reason for the discrepancy is not known, but we believe the new data to be more reliable”.
Therefore we can revise accordingly the data of the Figure 4, obtaining Figure 5.
Let us consider the value of the thermal conductivity of the enriched Silicon at 200 K. From the Figure 4, we have a value of 400 W/m-K: if we reduce it of 90 W/m-K, we obtain a value of 3 W/cm-K (Figure 5). Let us note that this value is in good agreement with the result of the theoretical model in . From the plot, the conductivity of natural Silicon seems being 2.2 W/cm-K. The ratio of the conductivity of enriched and natural Silicon seems being 1.36. Then, the enrichment increases, at 200 K, the conductivity of 36% .
Let us consider again the results given in . In the Figure 6, we can see side by side the theoretical results given in the Ref.22 and the experimental data given in the Figure 5, with the corrections for the enriched sample (red squares).
Figure 6: On the left, the theoretical curves from , on the right the experimental data [26,30]. The red lines had been added to help the reader in the comparison.
We can discuss the comparison. First of all, I have removed the experimental data of the enriched Silicon below 3 K. These data depend on the boundary scattering and therefore on the size of the sample. For the enriched sample, the specific value of the thermal conductivity at the peak is different, however the thermal conductivity obtained from theory between 3 K and 10 K and between 30 K and 200 K is in good agreement with experimental data. For the natural sample, above 10 K, the overall behaviors are in agreement. Moreover, the theoretical curve is in agreement with the open circles from .
In , the researchers measured the thermal conductivity of isotopically enriched 28Si “independently in three laboratories by high precision experiments on a total of four samples of different shape and degree of isotope enrichment in the range from 5 to 300 K with particular emphasis on the range near room temperature. The results obtained in the different laboratories are in good agreement with each other. They indicate that at room temperature the thermal conductivity of isotopically enriched 28Si exceeds the thermal conductivity of Si with a natural, unmodified isotope mixture by 10±2%. This finding is in disagreement with an earlier report by Ruf et al.”. The results given in  are in agreement with the data from .
Let us note that the results for natural Silicon from , at 200 K, are rather different from data given in , and also from  and .
Fig.7: Data from .
From the Figure 7, at 200 K, the conductivity is comprised between 2.5 and 2.9 W/cm-K. The ratio of the conductivity is 1.16. Actually, the new measurements had increased the value of the thermal conductivity in the natural sample of 16%.
Let us stress this fact, that it is a good set of measurements on the natural sample which is relevant for any model of the phonon scattering by the isotopic mass defects. For what concerns the calculation in , the given model based on a real representation of the crystalline lattice and of the scattering processes produces values of the thermal conductivity in good agreement for the pure crystal. For the crystal with natural isotopic mass defects, the formula used for the scattering produced by these defects is giving, at 200 K, a conductivity 25% lower than the new experimental value. In the case of Diamond and Germanium, that were studied in  too, the agreement between theory and experiments is good on all the range of temperature.
In this paper we have presented some theoretical curves and experimental data on the thermal conductivity of crystalline Silicon. Let us stress that a good overall agreement of the theoretical curves with experimental values exists, in particular for enriched Silicon . The theoretical curves had been obtained from a realistic model, without the use of adjustable parameters. For what concerns the natural Silicon, the most recent data here considered give a thermal conductivity at 200 K which is about 10% larger than the previously determined value. It would be interesting to understand the origin of this difference in the experimental values.
 Ziman, J.M. (1962). Electrons and phonons: The theory of transport phenomena in solids. Clarendon Press. Reprint Edition in 2001.
 Gurevich, V.L. (1986). Transport in phonon systems. North-Holland.
 Srivastava, G.P. (1990). The physics of phonons. Taylor and Francis.
 Sparavigna, A. C. (2016). The Boltzmann equation of phonon thermal transport solved in the relaxation time approximation – I – Theory. Mechanics, Materials Science & Engineering Journal, 2016 (3), 34-45. Available ?hal-01280946?
 Sparavigna, A. C. (2016). The Boltzmann Equation Of Phonon Thermal Transport Solved In the Relaxation Time Approximation – II – Data Analysis. Mechanics, Materials Science & Engineering Journal, 2016 (3), 57-66. Available ?hal-01280952?
 Omini, M., & Sparavigna, A. (1996). Beyond the isotropic-model approximation in the theory of thermal conductivity. Physical Review B, 53 (14), 9064-9073.
 Omini, M., & Sparavigna, A. (1997). Heat transport in dielectric solids with diamond structure. Nuovo Cimento D, 19 (10), 1537-1564.
 Broido, D.A., Ward, A., & Mingo, N. (2005). Lattice thermal conductivity of silicon from empirical interatomic potentials. Physical Review B, 72 (1), 014308.
 Kaburaki, H., Li, J., Yip, S., & Kimizuka, H. (2007). Dynamical thermal conductivity of argon crystal. Journal of Applied Physics, 102 (4), 043514.
 Turney, J.E., Landry, E.S., McGaughey, A.J.H., & Amon, C.H. (2009). Predicting phonon properties and thermal conductivity from anharmonic lattice dynamics calculations and molecular dynamics simulations. Physical Review B, 79 (6), 064301.
 McGaughey, A.J., & Kaviany, M. (2004). Quantitative validation of the Boltzmann transport equation phonon thermal conductivity model under the single-mode relaxation time approximation. Physical Review B, 69 (9), 094303.
 Li, W., Carrete, J., Katcho, N.A., & Mingo, N. (2014). ShengBTE: A solver of the Boltzmann transport equation for phonons. Computer Physics Communications, 185 (6), 1747-1758.
 Anthony, T. R., Banholzer, W. F., Fleischer, J. F., Wei, L., Kuo, P. K., Thomas, R. L., & Pryor, R. W. (1990). Thermal diffusivity of isotopically enriched C 12 diamond. Physical Review B, 42 (2), 1104.
 Wei, L., Kuo, P. K., Thomas, R. L., Anthony, T. R., & Banholzer, W. F. (1993). Thermal conductivity of isotopically modified single crystal diamond. Physical Review Letters, 70 (24), 3764.
 Onn, D. G., Witek, A., Qiu, Y. Z., Anthony, T. R., & Banholzer, W. F. (1992). Some aspects of the thermal conductivity of isotopically enriched diamond single crystals. Physical Review Letters, 68 (18), 2806.
 Olson, J. R., Pohl, R. O., Vandersande, J. W., Zoltan, A., Anthony, T. R., & Banholzer, W. F. (1993). Thermal conductivity of diamond between 170 and 1200 K and the isotope effect. Physical Review B, 47 (22), 14850.
 Asen-Palmer, M., Bartkowski, K., Gmelin, E., Cardona, M., Zhernov, A. P., Inyushkin, A. V., Taldenkov, A., Ozhogin, V. I., Kohei M. Itoh, & Haller, E. E. (1997). Thermal conductivity of germanium crystals with different isotopic compositions. Physical Review B, 56 (15), 9431.
 Omini, M., & Sparavigna, A. (1995). An iterative approach to the phonon Boltzmann equation in the theory of thermal conductivity. Physica B, 212 (2), 101-112.
 Omini, M., & Sparavigna, A. (1997). Effect of phonon scattering by isotope impurities on the thermal conductivity of dielectric solids. Physica B, 233 (2-3), 230-240.
 Holland, M. G., & Neuringer, L. J. (1962). The effect of impurities on the lattice thermal conductivity of silicon. In Proc. Int. Conf. on Physics of Semiconductors (pp. 474-481).
 Srivastava, G. P. (1976). Calculation of lattice thermal conductivity of Ge from 4 to 900 K. Philosophical Magazine, 34 (5), 795-809.
 Sparavigna, A. (2002). Influence of isotope scattering on the thermal conductivity of diamond. Physical Review B, 65 (6), 064305.
 Holland, M. G. (1964). Phonon scattering in semiconductors from thermal conductivity studies. Physical Review, 134 (2A), A471.
 Capinski, W.S., Maris, H.J., Bauser, E., Silier, I., Asen-Palmer, M., Ruf, T., Cardona, M., & Gmelin, E., (1997). Thermal conductivity of isotopically enriched Si. Applied Physics Letters, 71 (15), 2109-2111.
 Holland, M. G. (1963). Analysis of lattice thermal conductivity. Physical Review, 132 (6), 2461.
 Ruf, T., Henn, R. W., Asen-Palmer, M., Gmelin, E., Cardona, M., Pohl, H-J., Devyatych, G. G., & Sennikov, P. G. (2000). Thermal conductivity of isotopically enriched silicon. Solid State Communications, 115 (5), 243-247.
 Touloukian, Y. S., Powell, R. W., Ho, C. Y., & Klemens, P. G. (1970). TPRC data Series. Thermal Conductivity Nonmetallic Solids, IFI Plenum Press, Vol.1.
 Capinski, W. S., Maris, H. J., & Tamura, S. (1999). Analysis of the effect of isotope scattering on the thermal conductivity of crystalline silicon. Physical Review B, 59 (15), 10105.
 Inyushkin, A. V. (2002). Thermal Conductivity of Isotopically Modified Silicon: Current Status of Research Inorganic Materials, 38 (5), 427–433. Translated from Neorganicheskie Materialy, 38(5), 527–534. Original Russian Text.
 Ruf, T., Henn, R. W., Asen-Palmer, M., Gmelin, E., Cardona, M., Pohl, H-J., Devyatych, G. G., & Sennikov, P. G. (2000). Erratum to ''Thermal conductivity of isotopically enriched silicon''. Solid State Communications, 115, 243–247.
 Kremer, R. K., Graf, K., Cardona, M., Devyatykh, G. G., Gusev, A. V., Gibin, A. M., Inyushkin, A. V., Taldenkov, A. N., & Pohl, H. J. (2004). Thermal conductivity of isotopically enriched 28 Si: revisited. Solid state communications, 131 (8), 499-503.
 Gusev, A. V., Gibin, A. M., Morozkin, O. N., Gavva, V. A., & Mitin, A. V. (2002). Thermal Conductivity of 28 Si from 80 to 300 K. Inorganic materials, 38 (11), 1100-1102.
Information about this Article
This Article has not yet been peer-reviewed
This Article was published on 30th October, 2017 at 11:12:07 and has been viewed 390 times.
This work is licensed under a Creative Commons Attribution 2.5 License.
The full citation for this Article is:|
Sparavigna, A. (2017). Thermal Conductivity of the Crystalline Silicon. PHILICA.COM Article number 1143.