SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Reduced temperature-dependent thermal conductivity of magnetite thin films by controlling film thickness

No-Won Park1, Won-Yong Lee1, Jin-A Kim2, Kyungjun Song3, Hyuneui Lim3, Wan-Doo Kim3, Soon-Gil Yoon2* and Sang-Kwon Lee1*

Author Affiliations

1 Department of Physics, Chung-Ang University, Seoul, 156–756, Republic of Korea

2 Department of Materials Engineering, Chungnam National University, Daejeon 305-764, Republic of Korea

3 Department of Nature-Inspired Nanoconvergence Systems, Korean Institute of Machinery and Materials (KIMM), Daejeon 305-343, Republic of Korea

For all author emails, please log on.

Nanoscale Research Letters 2014, 9:96  doi:10.1186/1556-276X-9-96

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/9/1/96

Received:15 January 2014
Accepted:19 February 2014
Published:26 February 2014

© 2014 Park et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.


We report on the out-of-plane thermal conductivities of epitaxial Fe3O4 thin films with thicknesses of 100, 300, and 400 nm, prepared using pulsed laser deposition (PLD) on SiO2/Si substrates. The four-point probe three-omega (3-ω) method was used for thermal conductivity measurements of the Fe3O4 thin films in the temperature range of 20 to 300 K. By measuring the temperature-dependent thermal characteristics of the Fe3O4 thin films, we realized that their thermal conductivities significantly decreased with decreasing grain size and thickness of the films. The out-of-plane thermal conductivities of the Fe3O4 films were found to be in the range of 0.52 to 3.51 W/m · K at 300 K. For 100-nm film, we found that the thermal conductivity was as low as approximately 0.52 W/m · K, which was 1.7 to 11.5 order of magnitude lower than the thermal conductivity of bulk material at 300 K. Furthermore, we calculated the temperature dependence of the thermal conductivity of these Fe3O4 films using a simple theoretical Callaway model for comparison with the experimental data. We found that the Callaway model predictions agree reasonably with the experimental data. We then noticed that the thin film-based oxide materials could be efficient thermoelectric materials to achieve high performance in thermoelectric devices.

Iron oxide (Fe3O4); Thermal conductivity; 2D thin films; 3-ω technique; Callaway model; In-plane and out-of-plane


In recent decades, there has been a great interest in the application of thermoelectric (TE) effects in alternative clean energy sources [1-6]. For the evaluation of the thermoelectric performances of TE devices, their efficiencies can usually be quantified by a dimensionless figure of merit (ZT), S2σT/κ or a power factor S2σ, where S is the Seebeck coefficient, σ is the electrical conductivity, κ is the thermal conductivity, and T is the absolute temperature. High-performance thermoelectric materials with high ZT values should have a large Seebeck coefficient, high electrical conductivity, and low thermal conductivity [2,7,8]. To obtain an efficiently comparable to a household refrigerator, a ZT value at least 3 is desired for more widespread applications [6]. Recently, several researchers have alternatively studied two-dimensional (2D) thin films [9,10] to overcome the limitations of 1D nanostructured materials whose thermal properties are highly dependent on their dimensionality and morphology [3,11-13]. In 2010, Tang et al. reported that the thermal conductivity of holey Si thin film consistently reduces by around 2 orders of magnitude with a reduction in the pitch of the hexagonal holey pattern down to approximately 55 nm with approximately 35% porosity [9]. Similarly, Yu et al. reported that a Si nanomesh structure exhibits a substantially lower thermal conductivity than an equivalently prepared array of Si nanowires [10]. Hence, we believe that the 2D materials (i.e., thin film formation) could be highly promising candidates as TE materials for scalable and practical TE device applications.

Magnetite (Fe3O4) is a well-known half-metallic material, whose electronic density of states is 100% spin polarized at the Fermi level [14,15]. These properties allow Fe3O4 to be a promising candidate for spintronic devices [16]. However, the thermal property of this metal compound has not been widely studied. In 1962, Slack extensively studied and analyzed the thermal conductivity of a single crystal of paramagnetic bulk Fe3O4 materials at temperatures of 3 to 300 K [17]. He found that the thermal conductivity of Fe3O4 falls sharply with increasing temperature at the approximately 121 ± 2 K transition and reported a notable effect of vacancy and impurities on Fe3O4, particularly below 30 K. The thermal conductivity of pure Fe3O4 was as low as approximately 6 W/m · K at 300 K, owing to phonon scattering by local disorder in the materials, thus implying that pure Fe3O4 is a promising TE material. To the best of our knowledge, there have been no studies on the thermal properties of Fe3O4 thin films.

In this work, we present the out-of-plane thermal conductivities of epitaxial Fe3O4 thin films with thicknesses of 100 to 400 nm having different grain sizes and surface roughness. The films were grown at a deposition temperature of 300°C using pulsed laser deposition (PLD). We successfully demonstrated the temperature-dependent thermal conductivities of epitaxial Fe3O4 thin films via four-point probe 3-ω method in the temperature range of 20 to 300 K. The measured out-of-plane thermal conductivities of the Fe3O4 thin films (0.52 to 3.51 W/m · K) at 300 K are considerably reduced compared to those of the bulk materials (approximately 6 W/m · K) [17] because of strongly enhanced phonon-boundary scattering, as expected in the Callaway model [18]. Furthermore, we clearly realized that the thermal conductivity increased with an increase in film thickness and grain size, which agreed well with the theoretical predictions of the Callaway model.


The epitaxial magnetite thin films were synthesized on SiO2/Si (100) substrates at a temperature of 300°C using PLD. The detailed growth processes can be found in our previous publication [19]. In brief, a krypton fluoride (KrF, 248 nm in wavelength) excimer laser whose energy density was approximately 2.1 J/cm2 at repetition rate of 4 Hz at a pressure of 10-3 Pa was used along with a ceramic target (pure, homogeneous, and highly dense α-Fe2O3 ceramic). Our previous results confirmed that the surface roughness of the films increased with increasing temperature. Consequently, the deposition temperature was maintained at 300°C to obtain a uniform quality in the grown films. The deposition rate of the films was maintained at approximately 1.2 nm/min. To measure the thermal conductivity, we prepared three Fe3O4 thin films with thicknesses of 100, 300, and 400 nm using PLD. X-ray diffraction confirmed that the films were grown with a (111) preferred orientation with high-quality epitaxial growth, as detected from the in-plane phi-scans of the films [19]. Figure 1a,b,c shows the cross-sectional scanning electron microscope (SEM) images of the as-grown Fe3O4 thin films, confirming that the thicknesses of the films were in the range of 100 to 400 nm. Atomic force microscope (AFM) images (insets of Figure 1ab,c) showed that the grown films exhibit smooth grain morphologies with a root-mean-square (rms) roughness of 1.4 to 6.0 nm, as summarized in Figure 1d. We also found that the grain size of the films increased from approximately 13.2 ± 5.2 nm to approximately 230 ± 23.10 nm when the film thickness was increased from 100 to 400 nm, indicating that thicker films have much rougher surface morphology and larger grain size.

thumbnailFigure 1. SEM cross-sectional images of Fe3O4 thin films grown on a SiO2/Si substrate at 300°C using PLD. (a) 100 nm, (b) 300 nm, and (c) 400 nm. The insets show the AFM images of each thin film. (d) A summary of the prepared Fe3O4 thin film, including rms roughness, film thickness, deposition time, and grain size information.

Prior to the thermal conductivity measurements, a 100-nm-thick SiO2 layer was deposited on the thin film through plasma-enhanced chemical vapor deposition for electrical insulation of the films. Finally, a narrow metal strip (Ti/Au = 10/300 nm) consisting of four-point probe electrodes acting as a heater wire and probe pads was patterned onto the specimen through a conventional photolithography process. The thermal transport measurements were performed in closed cycle refrigerator (CCR) system with a shielding box, as shown in Figure 2a, which helped maintain the temperature in the range of 20 to 300 K and provided a high-vacuum (approximately 10-6 Torr) environment to avoid heat loss. In the current study, we utilized a four-point probe 3-ω method based on the application of an alternating current (AC) with angular modulation frequency (1-ω), which was first developed by Cahill in 1990 [20] to measure the temperature-dependent thermal conductivities of as-grown Fe3O4 thin films. It has been proved the most promising technique to extract thermal conductivities of 1D nanostructures such as nanowires [21,22] and carbon nanotubes [23,24] and thin films [25-27]. We have also proved this technique to be one of the powerful methods to extract the thermal conductivity of most low-dimensional materials [21]. Our experimental setup reported previously [21] is similar to the original design by Cahill [20] and adheres to the experimental design by Feser et al.[25].

thumbnailFigure 2. Experimental setup including the circuit connections with thermal management and electrical measurement systems. Experimental setup and circuit (a) and the corresponding circuit (right side) (b), equipped with thermal management and electrical measurement systems for thermal conductivity measurements via the 3-ω method at temperature ranges of 20 to 300 K.

Figure 2a,b shows the experimental setup including the circuit connections with thermal management and electrical measurement systems for out-of-plane thermal conductivity measurements via the 3-ω method. In brief, the sample was first attached to a printed circuit board substrate with vacuum grease for mounting inside a CCR with a shielding box. The source meter (Keithley 6221, Cleveland, OH, USA) was connected to both metallic pads to generate an AC (I0), as shown in Figure 2b. I0 with an angular modulation frequency of 1-ω was applied to generate Joule heat and temperature fluctuations at a frequency of 2-ω. The resistance of the narrow metal strip is proportional to the temperature that leads to a voltage fluctuation (V = IR) of 3-ω across the specimen. A lock-in amplifier (A-B mode, SR-850, Stanford Research System, Sunnyvale, CA, USA) connected to the two electrodes in the middle received the 3-ω voltage fluctuation along the narrow metal strip; this gives the information on the thermal conductivity of the films (as indicated in Figure 2b). To measure the thermal conductivity of the thin films, we then plotted the third-harmonic voltage (V3ω) against the natural logarithm of the applied frequencies (ln ω), which showed a linear relationship. Consequently, we determined the out-of-plane thermal conductivity from the slope in the linear region. The difference between two V3ω values (i.e., V3ω1 and V3ω2) is equated to the temperature drop across the Fe3O4 film and is used to calculate the cross-plane thermal conductivity, which is defined by the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M1">View MathML</a>


Here, V0 and R0 are the applied voltage and electrical resistance, respectively, along the heater wire of length l. <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M2">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M3">View MathML</a> are the third-harmonic voltages at input current frequencies of ω1 and ω2, respectively, and dR/dT (temperature coefficient resistance, TCR) is the rate of the resistance change of the heater at temperatures of 20 to 300 K. Figure 3a shows a schematic of the four-point probe electrodes patterned onto SiOx/Fe3O4/SiO2/Si substrate for thermal conductivity measurements using the 3-ω method. To confirm our results of thermal conductivity measured using the four-point probe 3-ω method, we used bismuth (Bi) films (50 nm in thickness) whose thermal conductivity is well known, as a reference sample. We determined its thermal conductivity to be 2.7 to 2.9 W/m · K, which is in good agreement with the previous reported results by Völklein and Kessler [28] and Völklein et al. [29] who reported that the thermal conductivity of 60-nm Bi thin films was approximately 3.6 W/m · K at 300 K. Thus, our experimental setup and the associated analysis via the four-point probe 3-ω method were clearly validated through a comparison with the results for reference sample. Figure 3b shows temperature-dependent resistances of the three Fe3O4 thin films (100, 300, 400 nm in thickness) in the temperature range of 20 to 300 K. The relationship between the resistance changes in the heater wire and the temperature is linear. Figure 3b shows that the TCR for the 100-, 300-, and 400-nm Fe3O4 thin films is approximately 0.104 Ω/K, approximately 0.041 Ω/K, and approximately 0.026 Ω/K, respectively. These values can be used for estimating thermal conductivity as defined in Equation 1.

thumbnailFigure 3. Four-point probe 3-ω method and temperature-dependent resistances. (a) Schematic view of the four-point probe 3-ω method where the out-of-plane thermal conductivity can be measured. (b) The temperature-dependent resistances of three Fe3O4 thin films (100, 300, 400 nm in thickness) at temperature ranges of 20 to 300 K.

Results and discussion

To ensure that the measured V3ω signal is generated by the Fe3O4 thin film, we investigated the variation in the signal with the applied frequency (ln ω) from the 3-ω measurements. This applied frequency usually provides a suitable current range for an estimation of the V3ω signal from the sample. As discussed previously by Cahill [20], the linear relationship of ln ω with V3ω should be satisfied as shown in Figure 4a. Figure 4a presents the V3ω distribution of the 100-nm Fe3O4 thin film for different applied frequencies. The V3ω signal follows ln ω dependence very well in the frequency range of 205 to 495 Hz, which agrees well with Equation 1. To validate the measured V3ω signal and the thermal conductivity (κ) from the 3-ω measurements, we studied the applied current dependence on the thermal conductivity by applying an AC of 5 to 10 μA. As shown in Figure 4b, the measured thermal conductivities of the films with thicknesses of 100, 300, and 400 nm were approximately 0.52 ± 0.05, approximately 1.92 ± 0.06, and approximately 3.51 ± 0.12 W/m · K, respectively, in the applied current range, indicating that κ is independent of the applied current (I0). We found that the errors in the thermal conductivity measurements are less than approximately 3% to 9%, depending on the film thickness.

thumbnailFigure 4. V3ωdistribution and thermal conductivities of the Fe3O4 film. (a) Linear regions of the third-harmonic voltage versus the applied frequency at various applied alternating currents (AC) ranging from 5 to 10 μA. (b) Thermal conductivities of Fe3O4 film with different film thicknesses (100, 300, and 400 nm) with respect to the applied AC (5 to 10 μA). Variation in the thermal conductivity with modulation of the input AC current could be assumed as measurement errors in thermal conductivity.

Figure 5a shows the temperature dependence of out-of-plane thermal conductivity of three Fe3O4 films at temperatures of 20 to 300 K and a simple theoretical calculation based on the Callaway model (solid lines in the figure) to compare with the experimental results (discussed in the next section). For the 400-nm-thick films, the thermal conductivity increased with increasing temperature up to approximately 40 K, then decreased with increasing temperature up to 300 K. Similar behaviors were observed for the other thin films (100 and 300 nm), as shown in Figure 5a. The phonon-phonon Umklapp and phonon-boundary scattering play an important role in phonon transport, and thus, the thermal conductivity decreases with increasing temperature [30,31]. Thus, we characterized the peaks of thermal conductivity (Umklapp peak) for the thin films whose thicknesses were 100, 300, and 400 nm, respectively. Our results presented in Figure 5a show that with the decrease in the film thickness from 400 to 100 nm, the corresponding Umklapp peaks shifted by approximately 20 K. According to the previous work in bulk F3O4, the Umklapp peak was generally observed at approximately 30 K [17], which is much lower than that for the thin films (approximately 40 to 60 K as shown in Figure 5a). From the shift in the Umklapp peaks, we can also confirm that phonon-boundary scattering is clearly dominant in the films in the temperature range of 40 to 60 K as a result of the grain size and film thicknesses [32,33]. In addition, when the temperature is above 50 K, the phonon-phonon Umklapp scattering becomes more pronounced. Our observation was in good agreement with a previous report on the thermal conductivity of 1D Bi nanowires [21]. Figure 5b presents the dependence of out-of-plane thermal conductivity on the film thickness at temperatures of 100, 200, and 300 K. The values of κ for the corresponding film thicknesses 100, 300, and 400 nm at 300 K increased gradually to approximately 0.52, approximately 1.85, and approximately 3.51 W/m · K, respectively. We also found that the thermal conductivities of the films were 1.7 to 11.5 times lower than that of bulk Fe3O4 (approximately 6 W/m · K) [17]. It has been well understood that the significant reduction in the thermal conductivity of the thin films (100 to 400 nm in thickness) compared to the bulk materials could be due to the enhanced phonon-boundary scattering in thin films predicted previously by Callaway [18]. In addition, we added the theoretical calculation results of Callaway's model in the same figure (solid line in Figure 5a,b). The results predicted by the Callaway model agree reasonably well with the experimental data, including the results for bulk Fe3O4. We can thus confirm that the significant reduction in the thermal conductivity for nanoscale thin films is principally a result of phonon-boundary scattering. In the following section, the calculation model is discussed in detail.

thumbnailFigure 5. Temperature-dependent conductivities of three Fe3O4 films and a simple theoretical calculation based on the Callaway model. (a, b) Measured thermal conductivities of 100-, 300-, and 400-nm-thick Fe3O4 thin films at temperatures of 20 to 300 K using the 3-ω method, including the thermal conductivity of bulk materials. The solid line denotes thermal conductivity of bulk materials from the theoretical Callaway model, which includes the effect of the impurity, Umklapp process, boundary scattering with film grain size, and film thickness.

To determine the temperature dependence of the thermal conductivity, κ(T), in Fe3O4 thin films quantitatively, we performed a theoretical calculation (i.e., fitting) based on the relaxation time model using the following expression predicted by Callaway in 1959 [18]:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M4">View MathML</a>


where ω is the phonon frequency, kB is the Boltzman constant, is the reduced Planck constant, x denotes the dimensionless parameter, x = ℏω/kBT, θD is the Debye temperature, T is the absolute temperature, and c is the velocity of sound. The total combined phonon scattering rate (relaxation time, τc) is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M5">View MathML</a>


where d1 is the grain size of the thin films (approximately 13.2, approximately 86, approximately 230 nm for the 100-, 300-, and 400-nm-thick films, respectively, from the AFM measurements shown in Figure 1), A and B are independent parameters of temperature and fitting, respectively, and c is the sound velocity, which is highly dependent on the direction of movement of phonons (average c = 2,500 m/s) [17]. To add the film thickness in Equation 3, we modified the phonon scattering rate given as

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M6">View MathML</a>


where d2 is the corresponding film thickness. For the Fe3O4 films, we estimated that the values of A and B in Equation 4 were numerically optimized as approximately 8.46 × 10-43S3 and approximately 7.89 × 10-18S/K, respectively, from the fitting to the bulk material values [17]. According to the Callaway model in Equations 3 and 4, the first term <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M7">View MathML</a> represents the boundary scattering; the second term Aω4 represents the scattering by point impurities or isotopes, and the third term <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M8">View MathML</a> represents the Umklapp process. Theoretical fits of the temperature dependence of the out-of-plane thermal conductivities of the Fe3O4 films from 20 to 300 K of Equations 2 and 4, which were obtained using the commercial application Mathematica (http://www.wolfram.com webcite), are compared with the experimental results in Figure 5a,b. From the numerical calculation of the temperature dependence of thermal conductivity, it was noted that the κ values indisputably decreased when the grain size was reduced, indicating that the effect of the nano-grained thin films on the thermal conductivity is essentially due to the relaxation time model based on phonon-boundary scattering. As shown in Figure 5a,b, the theoretical modeling based on the Callaway model agrees well quantitatively with the experimental data even though there is a difference in the κ values between the theoretical and experimental results for the 100-nm Fe3O4 film. The measured thermal conductivity results in the 100-nm films were approximately five times lower than the Callaway model prediction. This deviation can be explained by two arguments. First, the deviation in the thermal conductivity for the 100-nm thick film could be explained by the boundary effect, i.e., surface boundary scattering of the thinner films, in which the surface boundary scattering is more dominant compared to that of bulk and bulk-like thicker films, providing more phonon-boundary effect in thermal conductivity. However, in our theoretical model, no size and surface boundary scattering effects were considered. Thus, the measured temperature dependence of the thermal conductivity (0.52 W/m · K at 300 K) was relatively lower than the results expected from the theoretical calculation (1.9 to 2.4 W/m · K at 300 K), as shown in Figure 5b [2,34,35]. Previously, Li et al. also reported a similar observation for the thermal conductivity of Bi2Se3 nanoribbon [36]. Second, to numerically calculate the thermal conductivity using the Callaway model, we used the fitting parameters of A and B in the relaxation rate from the bulk materials. Thus, the theoretical calculation could be closer to the bulk material values. To clearly understand this inconsistency between the theoretical and experimental results, especially in nanoscale thin films (100-nm thin film in our case), the size and surface boundary effects in the Callaway model should be studied in detail for 1D and 2D nanostructures.

To measure the in-plane thermal conductivity of the nanoscale thin films from the out-of-plane thermal conductivity, we re-analyzed the theoretical Callaway model with the longitudinal component of sound velocity and Debye temperature, which were highly dependent on the direction of movement of phonons. In this calculation, we assumed that our measured results are close to the theoretical prediction, as shown in Figure 5a,b. The average sound velocity can be decomposed into the transverse and longitudinal components as defined in [37]:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M9">View MathML</a>


where cT (approximately 2,305.4 m/s) and cL (approximately 3,263.3 m/s) are the transverse and longitudinal velocities, respectively. In addition, the Debye temperature depends on sound velocity. Thus, using the calculated cT in Equation 5, we could calculate the Debye temperature for transverse component (<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M10">View MathML</a>) given as

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M11">View MathML</a>


where V is the volume per atom (approximately 10.54 × 10-30 m3). The Debye temperature with transverse sound velocity is then determined to be approximately 313.1 K. Finally, we calculated the in-plane thermal conductivity of the Fe3O4 films with transverse components of sound velocity (cT) and Debye temperature (<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/96/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/96/mathml/M12">View MathML</a>) using Equation 2. Figure 6a,b shows calculated both in-plane and out-of-plane thermal conductivities of 100-, 300-, and 400-nm-thick Fe3O4 thin films at temperatures of 20 to 300 K obtained using the simple Callaway phonon scattering model. As shown in Figure 6a, the deviation in thermal conductivity between the out-of-plane and in-plane thermal conductivities decreased with increasing temperature. At room temperature, the out-of-plane and in-plane thermal conductivities were determined to be 1.7 to 3.0 and 1.6 to 2.8 W/m · K, respectively. It was also noticed that the calculated out-of-plane thermal conductivity values are slightly higher than the in-plane thermal conductivity values in the Fe3O4 thin film as shown in Figure 6. This behavior could be due to the columnar structures of the grains (see Figure 1), where the phonons moving transversally in the Fe3O4 films are scattered by the columnar grains in the films. Similar results can be seen in diamond thin film grown by chemical vapor deposition (CVD), where the measured out-of-plane thermal conductivity consistently show a higher thermal conductivity along the columnar grains than the in-plane thermal conductivity [38].

thumbnailFigure 6. Thermal conductivities of 100-, 300-, and 400-nm-thick Fe3O4 thin films. (a, b) Calculated thermal conductivities of 100-, 300-, and 400-nm-thick Fe3O4 thin films at temperatures of 20 to 300 K obtained using the simple Callaway phonon scattering model. The temperature-dependent in-plane thermal conductivity was calculated by modifying the Debye temperature and sound velocity in the Callaway phonon scattering model.


In summary, we first present the thermal conductivity of epitaxial Fe3O4 thin films with thicknesses of 100 to 400 nm prepared on SiO2/Si (100) substrates using PLD. By measuring the temperature-dependent thermal characteristics of three Fe3O4 thin films using the effective four-point probe 3-ω method, we found that the thermal conductivities of the films are greatly reduced when compared with those of the corresponding bulk materials and that the thermal conductivity decreases with decreasing film thickness from 400 to 100 nm. Both theoretical and experimental results indicate that the Umklapp peaks of the thermal conductivity of Fe3O4 films move toward higher temperatures with decreasing film thickness, owing to the phonon-boundary scattering. The thermal conductivity was found to be in the range of 0.52 to 3.51 W/m · K at 300 K, which was 1.7 to 11.5 orders of magnitude lower than that of bulk materials at 300 K. We found that the modified Callaway theoretical model including the film thickness effect agreed well with the results obtained using the 3-ω method. Furthermore, we indirectly measured the in-plane thermal conductivity by re-analyzing the Callaway model using the measured out-of-plane thermal conductivity. We then suggested that the thin film-based oxide materials could be promising candidates as thermoelectric materials to achieve high-performance TE devices.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NWP and WYL, and JAK carried out all the experiments and analysis including the sample growth. KS, HEL, SGY, and WDK helped discuss the sample analysis and provided part of the financial support. SKL organized the final manuscript. All authors read and approved the final manuscript.


This study was supported by the Priority Research Centers Program and by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2013R1A12012685, NRF-2013R1A4A1069528). This study was also supported by a grant from the Global Excellent Technology Innovation R&D Program funded by the Ministry of Knowledge Economy, Republic of Korea (10038702-2010-01).


  1. Majumdar A: Thermoelectricity in semiconductor nanostructures.

    Science 2004, 303:777-778. PubMed Abstract | Publisher Full Text OpenURL

  2. Hochbaum AI, Chen RK, Delgado RD, Liang WJ, Garnett EC, Najarian M, Majumdar A, Yang PD: Enhanced thermoelectric performance of rough silicon nanowires.

    Nature 2008, 451:163-U5. PubMed Abstract | Publisher Full Text OpenURL

  3. Li DY, Wu YY, Kim P, Shi L, Yang PD, Majumdar A: Thermal conductivity of individual silicon nanowires.

    Appl Phys Lett 2003, 83:2934-2936. Publisher Full Text OpenURL

  4. Lim JW, Hippalgaonkar K, Andrews SC, Majumdar A, Yang PD: Quantifying surface roughness effects on phonon transport in silicon nanowires.

    Nano Lett 2012, 12:2475-2482. PubMed Abstract | Publisher Full Text OpenURL

  5. Kim DH, Kim C, Ha DW, Kim H: Fabrication and thermoelectric properties of crystal-aligned nano-structured Bi2Te3.

    J Alloys Comp 2011, 509:5211-5215. Publisher Full Text OpenURL

  6. DiSalvo FJ: Thermoelectric cooling and power generation.

    Science 1999, 285:703-706. PubMed Abstract | Publisher Full Text OpenURL

  7. Kim W, Wang R, Majumdar A: Nanostructuring expands thermal limits.

    Nano Today 2007, 2:40-47. OpenURL

  8. Kim W, Singer SL, Majumdar A, Zide JMO, Klenov D, Gossard AC, Stemmer S: Reducing thermal conductivity of crystalline solids at high temperature using embedded nanostructures.

    Nano Lett 2008, 8:2097-2099. PubMed Abstract | Publisher Full Text OpenURL

  9. Tang JY, Wang HT, Lee DH, Fardy M, Huo ZY, Russell TP, Yang PD: Holey silicon as an efficient thermoelectric material.

    Nano Lett 2010, 10:4279-4283. PubMed Abstract | Publisher Full Text OpenURL

  10. Yu JK, Mitrovic S, Tham D, Varghese J, Heath JR: Reduction of thermal conductivity in phononic nanomesh structures.

    Nat Nanotechnol 2010, 5:718-721. PubMed Abstract | Publisher Full Text OpenURL

  11. Jeon SJ, Jeon H, Na S, Kang SD, Lyeo HK, Hyun S, Lee HJ: Microstructure evolution of sputtered BiSb-Te thermoelectric films during post-annealing and its effects on the thermoelectric properties.

    J Alloy Compd 2013, 553:343-349. OpenURL

  12. Shi L, Hao Q, Yu CH, Mingo N, Kong XY, Wang ZL: Thermal conductivities of individual tin dioxide nanobelts.

    Appl Phys Lett 2004, 84:2638-2640. Publisher Full Text OpenURL

  13. Wang JA, Wang JS: Carbon nanotube thermal transport: ballistic to diffusive.

    Appl Phys Lett 2006, 88:111909. Publisher Full Text OpenURL

  14. Wolf SA, Awschalom DD, Buhrman RA, Daughton JM, von Molnar S, Roukes ML, Chtchelkanova AY, Treger DM: Spintronics: a spin-based electronics vision for the future.

    Science 2001, 294:1488-1495. PubMed Abstract | Publisher Full Text OpenURL

  15. Versluijs JJ, Bari MA, Coey JMD: Magnetoresistance of half-metallic oxide nanocontacts.

    Phys Rev Lett 2001, 87:026601. OpenURL

  16. Zutic I, Fabian J, Das Sarma S: Spintronics: fundamentals and applications.

    Rev Mod Phys 2004, 76:323-410. Publisher Full Text OpenURL

  17. Slack G: Thermal conductivity of MgO, Al2O3, MgAl2O4 and Fe3O4 crystals from 3 to 300 K.

    Phys Rev 1962, 126:427-441. Publisher Full Text OpenURL

  18. Callaway J: Model for lattice thermal conductivity at low temperatures.

    Phys Rev 1959, 113:1046-1051. Publisher Full Text OpenURL

  19. Yun JG, Lee YM, Lee WJ, Kim CS, Yoon SG: Selective growth of pure magnetite thin films and/or nanowires grown in situ at a low temperature by pulsed laser deposition.

    J Mater Chem C 2013, 1:1977-1982. Publisher Full Text OpenURL

  20. Cahill DG: Thermal-conductivity measurement from 30-K to 750-K- the 3-omega method.

    Rev Sci Instrum 1990, 61:802-808. Publisher Full Text OpenURL

  21. Lee SY, Kim GS, Lee MR, Lim H, Kim WD, Lee SK: Thermal conductivity measurements of single-crystalline bismuth nanowires by the four-point-probe 3-omega technique at low temperatures.

    Nanotechnology 2013, 24:185401. PubMed Abstract | Publisher Full Text OpenURL

  22. Lee KM, Choi TY, Lee SK, Poulikakos D: Focused ion beam-assisted manipulation of single and double beta-SiC nanowires and their thermal conductivity measurements by the four-point-probe 3-omega method.

    Nanotechnology 2010, 21:125301. PubMed Abstract | Publisher Full Text OpenURL

  23. Choi TY, Poulikakos D, Tharian J, Sennhauser U: Measurement of the thermal conductivity of individual carbon nanotubes by the four-point three-omega method.

    Nano Lett 2006, 6:1589-1593. PubMed Abstract | Publisher Full Text OpenURL

  24. Choi TY, Poulikakos D, Tharian J, Sennhauser U: Measurement of thermal conductivity of individual multiwalled carbon nanotubes by the 3-omega method.

    Appl Phys Lett 2005, 87:013108. Publisher Full Text OpenURL

  25. Feser JP, Chan EM, Majumdar A, Segalman RA, Urban JJ: Ultralow thermal conductivity in polycrystalline CdSe thin films with controlled grain size.

    Nano Lett 2013, 13:2122-2127. PubMed Abstract | Publisher Full Text OpenURL

  26. Feser JP, Sadhu JS, Azeredo BP, Hsu KH, Ma J, Kim J, Seong M, Fang NX, Li XL, Ferreira PM, Sinha S, Cahill DG: Thermal conductivity of silicon nanowire arrays with controlled roughness.

    J Appl Phys 2012, 112:114306. Publisher Full Text OpenURL

  27. Wang ZJ, Alaniz JE, Jang WY, Garay JE, Dames C: Thermal conductivity of nanocrystalline silicon: importance of grain size and frequency-dependent mean free paths.

    Nano Lett 2011, 11:2206-2213. PubMed Abstract | Publisher Full Text OpenURL

  28. Völklein F, Kessler E: A method for the measurement of thermal-conductivity, thermal-diffusivity, and other transport-coefficients of thin-films.

    Phys Status Solidi A 1984, 81:585-596. Publisher Full Text OpenURL

  29. Völklein F, Reith H, Cornelius TW, Rauber M, Neumann R: The experimental investigation of thermal conductivity and the Wiedemann-Franz law for single metallic nanowires.

    Nanotechnology 2009, 20:325706. PubMed Abstract | Publisher Full Text OpenURL

  30. Bui CT, Xie R, Zheng M, Zhang Q, Sow CH, Li B, Thong JT: Diameter-dependent thermal transport in individual ZnO nanowires and its correlation with surface coating and defects.

    Small 2012, 8:738-745. PubMed Abstract | Publisher Full Text OpenURL

  31. Guthy C, Nam CY, Fischer JE: Unusually low thermal conductivity of gallium nitride nanowires.

    J Appl Phys 2008, 103:064319. Publisher Full Text OpenURL

  32. Jezowski A, Danilchenko BA, Bockowski M, Grzegory I, Krukowski S, Suski T, Paszkiewicz T: Thermal conductivity of GaN crystals in 4.2-300 K range.

    Solid State Commun 2003, 128:69-73. Publisher Full Text OpenURL

  33. Mamand SM, Omar MS, Muhammad AJ: Nanoscale size dependence parameters on lattice thermal conductivity of Wurtzite GaN nanowires.

    Mater Res Bull 2012, 47:1264-1272. Publisher Full Text OpenURL

  34. Boukai AI, Bunimovich Y, Tahir-Kheli J, Yu JK, Goddard WA, Heath JR: Silicon nanowires at efficient thermoelectric materials.

    Nature 2008, 451:168-171. PubMed Abstract | Publisher Full Text OpenURL

  35. Sansoz F: Surface faceting dependence of thermal transport in silicon nanowires.

    Nano Lett 2011, 11:5378-5382. PubMed Abstract | Publisher Full Text OpenURL

  36. Li GD, Liang D, Qiu RLJ, Gao XPA: Thermal conductivity measurement of individual Bi2Se3 nano-ribbon by self-heating three-omega method.

    Appl Phys Lett 2013, 102:033106. Publisher Full Text OpenURL

  37. Alvarez-Quintana J, Martinez E, Perez-Tijerina E, Perez-Garcia SA, Rodriguez-Viejo J: Temperature dependent thermal conductivity of polycrystalline ZnO films.

    Appl Phys Lett 2010, 107:063713. OpenURL

  38. Garebner JE, Reiss ME, Seibles L: Phonon scattering in chemical-vapor-deposited diamond.

    Phys Rev B 1994, 50:3702-3713. Publisher Full Text OpenURL