We report on the reduction of the thermal conductivity in ultra-thin suspended Si membranes with high crystalline quality. A series of membranes with thicknesses ranging from 9 nm to 1.5 μm was investigated using Raman thermometry, a novel contactless technique for thermal conductivity determination. A systematic decrease in the thermal conductivity was observed as reducing the thickness, which is explained using the Fuchs-Sondheimer model through the influence of phonon boundary scattering at the surfaces. The thermal conductivity of the thinnest membrane with d = 9 nm resulted in (9 ± 2) W/mK, thus approaching the amorphous limit but still maintaining a high crystalline quality.
The long quest for thermoelectric materials with an improved figure of merit (ZT = S2σ/κL + κe, where S is the Seebeck coefficient and σ is the electrical conductivity) has mainly followed the roadmap of spoiling the lattice component of the thermal conductivity, κL, due to its larger contribution as compared to its electrical counterpart, κe, in undoped semiconductors. So far, this concept was based on producing materials with an inherent nanostructuration or simply including defects and/or impurities with the purpose of scattering the acoustic phonons responsible for most of the heat transport, e.g., for Si bulk λ ≈ 1 μm at 300 K.1,2 This approach was shown to be extremely efficient in a variety of material systems such as, e.g., Bi2Te3, PbTe, ZnO, Si1-xGex, and Si.3,4 Although the case of Si results technologically relevant due to its leading role in the semiconductor industry, its large κL in its bulk form of 148 W/mK5 at room temperature leads to a poor figure of merit of ZT300K ≈ 0.01. This naturally prevents the use of Si in thermoelectric modules designed for room temperature operation. Many recent publications have demonstrated that by smart engineering of Si in its one- or two-dimensional forms (i.e., nanowires (NWs) or thin films, respectively) it is possible to reduce κL by almost two orders of magnitude with respect to its bulk value. The reduction of the thermal conductivity in these systems has been associated to two principal factors: (i) the modification of the acoustic dispersion relations due to the additional periodicity (superlattices and phononic crystal structures)6–10 or spatial confinement of the phonon modes (thin films and nanowires);12–14 and (ii) the shortening of the phonon mean free path due to diffuse scattering of phonons at the boundaries.15–18
Patterning of suspended Si membranes has shown to be an effective method to decreases the κL. Values as low as ∼2 W/mK at room temperature were found in Refs. 10 and 19, thus approaching the amorphous limit of ∼1.8 W/mK.20 The case of Si nanowires also results exciting due to their large surface-to-volume ratio, i.e., κL ≈ 7 W/mK21 was obtained for NWs with 22 nm in diameter, and κL ≈ 5 W/mK was obtained for NWs with 70 nm in diameter and large surface roughness.22 Although 1D systems such as NWs are ideal regarding κL reduction, their applicability to produce devices is by no means straightforward due to the obvious technical complications arising from their manipulation and precise positioning. In fact, Si 2D systems patterned with holes10,19 led to similar reductions in κL as compared to Si NWs.21–25 However, the necessary patterning stage to introduce efficient phonon scatters (structural holes) introduces a large density of defect centres which might spoil the electrical conductivity of the devices, thus, leading to a poor ZT.
Several previous works have investigated the thermal properties of Si thin films using SOI-based (Silicon-on-Insulator) electrical metrology techniques.26–29 The in-plane and out-of-plane components of the thermal conductivity were shown to depend on the thickness of the thin films. Although these measurements provided the correct thickness dependence, in most cases several error sources arising from thermal boundary resistances, heat spreading within the heater/thermometers, and heat leakage to the underlying SiO2 had to be considered and modelled. These problems can be easily overcome using free-standing Si membranes in combination with a contactless method for thermal conductivity determination such as Raman thermometry,30,31 time-domain thermoreflectance (TDTR),32,33 frequency-domain thermoreflectance (FDTR),34 thermal transient grating (TTG),35 etc. Free-standing ultra-thin Si membranes have recently attracted considerable attention, since they offer model system to study the influence of low dimensionality on the thermal properties as recently reported.36–38 For example, size effects on the acoustic dispersion relations were investigated by inelastic Brillouin light scattering spectroscopy. A dramatic decrease of at least one order of magnitude of the phase velocity of zero order flexural modes was recently observed,38 which could have a direct impact on the group velocity and, thus, on the thermal conductivity. Furthermore, reduction of the lifetime of the fundamental dilatational mode (D1) of about two orders of magnitude was observed by investigating the dynamics of the reflectivity of the membranes using ultrafast methods.37 Although the flexural and dilatational modes are not the only responsible for the heat transport in the membranes, it is expected that they reflect the influence of size effects in the acoustic branches and, thus, on the thermal properties.
In this work an extra step is taken to account for the influence of low dimensionality in the thermal properties of free-standing ultrathin Si membranes by measuring in a contactless fashion the thermal conductivity, namely, by Raman thermometry. The origin of the size-dependent reduction of the thermal conductivity is discussed using the well-established Fuchs-Sondheimer model. This model assumes that the main physical mechanism behind this reduction is the shortening of the phonon mean free path due to diffuse scattering of phonons at the boundaries. For the thinnest membrane with d = 9 nm we have observed a 16-fold reduction of the thermal conductivity with respect to the bulk value. Finally, we compare our data with previous experiments in supported Si thin films.
Figure 1(a) shows a schematic of the membranes and the measurements configuration. A SOI substrate was gradually oxidized to thin down the Si film to the desired thickness, followed by removal of the grown oxide. The thermal processes tend to create compressive stress in the Si film, thus leading to buckling of the membranes after release as shown in a representative optical micrograph in Fig. 1(b). We point out that full back etching of the substrate is essential to perform the Raman thermometry measurements, since the absorption coefficient of the membranes must be accurately determined in order to extract the thermal conductivity.
(a) Schematics of the Si membranes supporting structure and the measurement configuration showing the 514.5 nm incident laser, and the power meters (PI, PR, and PT) used to determine the power absorption coefficient of each membrane. (b) Optical micrographs of a representative Si membrane showing their surface corrugation arising from relaxation of built-in compressive strain. (c) Representative high-resolution transmission electron microscopy (HR-TEM) image of a representative high-stress Si membrane.11
(a) Schematics of the Si membranes supporting structure and the measurement configuration showing the 514.5 nm incident laser, and the power meters (PI, PR, and PT) used to determine the power absorption coefficient of each membrane. (b) Optical micrographs of a representative Si membrane showing their surface corrugation arising from relaxation of built-in compressive strain. (c) Representative high-resolution transmission electron microscopy (HR-TEM) image of a representative high-stress Si membrane.11
The Raman thermometry technique is based on focusing a laser beam onto the centre of the membranes and monitoring the Raman shift (redshift) of the longitudinal optical (LO) phonon of Si as a function of the absorbed power. The lattice temperature can be directly determined given a previous calibration of these quantities.39 The thermal conductivity is extracted from the knowledge of the absorbed power and the temperature rise at the spot by solving the 2D heat equation, ∇2T = PA/(κπa2d)e(− (r/a)2) where a is the spot radius and d the thickness of the membranes.40–45 To obtain the absorbed power in each membrane we first determine their power absorption coefficient. The power absorption (A), reflection (R), and transmission (T) coefficients of the Si membranes as a function of their thickness are shown in Fig. 2. The solid points represent the experimental data. The reflection (R = PR/PI) and transmission (T = PT/PI) coefficients were obtained measuring the incident (PI), reflected (PR), and transmitted (PT) powers after focusing the laser spot on the surface of the membranes using a 50× (NA = 0.55) microscope objective. The power absorption (A) coefficient of each membrane was computed considering that A = 1 − R − T = 1 − PR/PI − PT/PI. The absolute error in the determination of A was below 1%. The solid lines in Fig. 2 are the result from theoretical calculations at 300 K obtained applying the laws of reflection and refraction in a plane-parallel film and using the data for the dielectric function of Si reported in Ref. 46. We note that since the thickness of the membranes are well below, or comparable to the wavelength of the incident laser (λ = 514.5 nm), the membranes behave as Frabry-Pérot optical cavities. For this reason A, R, and T do not exhibit a monotonous response as in the case of bulk Si39 but rather an oscillatory behaviour. The agreement between the experimental data and the calculated curves results exceptionally good for the power absorption coefficient. We note that the data shown in Fig. 2 were obtained for low excitation powers to ensure no local heating of the membranes. Moreover, the absorption coefficient was also determined as a function of the incident power, since it systematically increases with the local temperature at the laser spot.
Power absorption (A), reflection (R), and transmission (T) coefficients as a function of the thickness of the membranes. Solid lines are calculations obtained from the laws of reflection and refraction in a plane-parallel geometry as shown in the inset.
Power absorption (A), reflection (R), and transmission (T) coefficients as a function of the thickness of the membranes. Solid lines are calculations obtained from the laws of reflection and refraction in a plane-parallel geometry as shown in the inset.
Figure 3 displays the temperature of each membrane at the laser spot as a function of the absorbed power, which is computed as PA = PI.A. The Si LO phonon Raman shift was measured with an accuracy of ±0.05 cm−1 after fitting the spectra with a Lorentzian lineshape. Using the temperature dependence of the Si LO phonon the Raman shift can be converted into local temperature with a corresponding accuracy of about ±2 K.39 The solid lines are fits to the data points using a linear relation. For any given absorbed power the temperature rise at the laser spot increases as the thickness of the membranes decreases. The 2D heat equation is solved using finite element (FE) simulations for a given power. We varied the thermal conductivity of each membrane until the temperature rise corresponding to a given power was obtained using the data from Fig. 3. We have chosen a wide range of absorbed powers for each membrane, since in this way the method provides the larger accuracy. Finally, we also show in Fig. 3 a representative solution of the 2D heat equation. The temperature distribution on the membranes exhibits a maximum at the laser spot position and gradually decays to room temperature towards the frame of the membranes, which is in contact with the substrate.
Raman shift (right axis) of the longitudinal optical (LO) Si phonon of the membranes as a function of the absorbed power and membrane's thickness. The right axis represents the temperature obtained from the temperature dependence of the LO mode extracted from Ref. 35. The full circular symbols represent the low-stress membranes, whereas stars symbols account for the stressed membranes. The inset shows a schematic of the experiments. The incident laser (λ = 514.5 nm) is focus at the centre of the membranes.
Raman shift (right axis) of the longitudinal optical (LO) Si phonon of the membranes as a function of the absorbed power and membrane's thickness. The right axis represents the temperature obtained from the temperature dependence of the LO mode extracted from Ref. 35. The full circular symbols represent the low-stress membranes, whereas stars symbols account for the stressed membranes. The inset shows a schematic of the experiments. The incident laser (λ = 514.5 nm) is focus at the centre of the membranes.
Figure 4 shows the normalized thermal conductivity (κ/κBulk) for the Si membranes as a function of their thickness (d). A systematic decrease is observed with decreasing the thickness, d. The dashed line shows the theoretical prediction of the dependence of κ on d using the Fuchs-Sondheimer model. The details of the calculations are presented in Ref. 47. The agreement between the theoretical model and the experimental data is satisfactory for the lower thicknesses. For the thicker membranes an underestimation is observed, which arises from an artifact of the model, since in this regime boundary scattering becomes gradually negligible. Regarding the influence of surface roughness we surmise that this may also contribute to the reduction of the thermal conductivity. Nevertheless, we expect this effect to be small at room temperature given the <0.5 nm RMS values obtained for all the membranes. A 16-fold reduction of the thermal conductivity (κ9nm ≈ 9 W/mK) with respect to the bulk value was observed in the thinnest membrane with d = 9 nm. In Fig. 4 the experimental data on SOI thin films from Refs. 26–29 as well as in thick Si membranes30 are also shown for comparison. The thermal conductivities of the SOI thin films are systematically larger than those of the free-standing membranes. This deviation may be associated with the detailed structure of the Si/SiO2 interface or with a different impurity concentration, but these are not yet experimentally confirmed. A more logical explanation is to consider that in the case of the SOI samples, the substrate acts as an extra heat sink, thus, leading to a larger effective thermal conductivity.
Thermal conductivities of the membranes, κ/κBulk, normalized to the bulk Si value as a function of the thickness are shown in solid red symbols. The theoretical description of the data using the Fuchs-Sonheimer model is shown in dashed lines. Comparison with previous works on SOI and membranes is also shown. The thermal conductivity of bulk Si was taken as κBulk = 148 W/mK.5 The amorphous Si (a-Si) limit is also shown in dotted-dashed line.
Thermal conductivities of the membranes, κ/κBulk, normalized to the bulk Si value as a function of the thickness are shown in solid red symbols. The theoretical description of the data using the Fuchs-Sonheimer model is shown in dashed lines. Comparison with previous works on SOI and membranes is also shown. The thermal conductivity of bulk Si was taken as κBulk = 148 W/mK.5 The amorphous Si (a-Si) limit is also shown in dotted-dashed line.
In conclusion, we have investigated the thickness dependence of the thermal conductivity in ultrathin free-standing Si membranes of high crystalline quality using Raman thermometry. The power absorption coefficient of the membranes was determined experimentally and theoretically calculated. We have found that the thermal conductivity of the membranes systematically reduces as their thickness decreases. This was successfully modelled considering the shortening of the phonon mean free path due to the diffuse scattering at the boundaries. The thermal conductivity of the thinnest membrane with d = 9 nm resulted in (9 ± 2) W/mK, which approaches the amorphous limit while still maintaining a high crystalline quality. The dependence of the thermal conductivity on dimensionality opens new possibilities for building thermoelectric modules based on suspended two-dimensional systems.
The authors acknowledge the financial support from the FP7 project MERGING (Grant No. 309150), the Spanish MICINN projects nanoTHERM (Grant No. CSD2010–0044), TAPHOR (MAT2012–31392), and the Academy of Finland (Grant No. 252598). E.C.-A. gratefully acknowledges a Becas Chile 2010 CONICYT fellowship from Chilean government.