OPEN

SUBJECT AREAS:

CHARACTERIZATION AND ANALYTICAL TECHNIQUES

Length Dependent Thermal Conductivity Measurements Yield Phonon Mean Free Path Spectra in Nanostructures

NANOWIRES Hang Zhang, Chengyun Hua, Ding Ding & Austin J. Minnich

Received 29 October 2014

Accepted 19 February 2015

Published 13 March 2015

Correspondence and requests for materials should be addressed to A.J.M. (aminnich@ caltech.edu)

Division of Engineering and Applied Science California Institute of Technology Pasadena, CA 91125.

Thermal conductivity measurements over variable lengths on nanostructures such as nanowires provide important information about the mean free paths (MFPs) of the phonons responsible for heat conduction. However, nearly all of these measurements have been interpreted using an average MFP even though phonons in many crystals possess a broad MFP spectrum. Here, we present a reconstruction method to obtain MFP spectra of nanostructures from variable-length thermal conductivity measurements. Using this method, we investigate recently reported length-dependent thermal conductivity measurements on SiGe alloy nanowires and suspended graphene ribbons. We find that the recent measurements on graphene imply that 70% of the heat in graphene is carried by phonons with MFPs longer than 1 micron.

Thermal transport in nanostructures has been a topic of intense interest in recent years1-3. When the characteristic dimensions of nanostructures such as the diameter of a nanowire approach phonon mean free paths (MFPs), the thermal conductivity can be substantially smaller than the bulk value due to scattering from sample boundaries. Significant thermal conductivity reductions have been observed in a number of nanos-cale systems, including nanowires4-6, nanotubes7, thin Si membranes8, and micron size beams at cryogenic temperatures9. This concept has been widely adopted in thermoelectrics applications10-13.

Understanding and engineering the thermal conductivity reduction in nanostructures requires knowledge of phonon scattering mechanisms in the form of the phonon MFPs. The MFP accumulation function, which we term the MFP spectrum in this work, has been demonstrated to be a particularly useful quantity to describe the values of the MFPs relevant for heat conduction14. In several works, information about MFPs was obtained by measuring the thermal conductivity over variable lengths of nanostructures such as nanotubes7, graphene rib-bons15 and SiGe nanowires16. If phonons have MFPs exceeding the distance between the heat source and sink their contribution to thermal conductivity is reduced compared to that in the bulk material, and thus the deviations of the measured thermal conductivity from the bulk value provide information on the phonon MFPs. However, prior studies extracted only an average MFP despite the fact that recent works have demonstrated that in many solids phonon MFPs vary over orders of magnitude, making the approximation of an average MFP for all phonons quite poor17,18.

In principle, information about the full MFP spectrum should be contained in these variable-length thermal conductivity measurements, just as the MFP spectrum can be obtained from thermal conductivity measurements performed over variable thermal length scales in MFP spectroscopy19. In particular, the method proposed by Minnich based on convex optimization20 should be applicable to the present situation provided that the suppression function that describes the effect of the finite length on the thermal conductivity can be identified. Li et al obtained the phonon MFP spectrum of graphite along the c-axis from thickness dependent thermal conductivities obtained with molecular dynamics simulations21, but their suppression function was not rigorously obtained from the Boltzmann Transport Equation (BTE).

In this report, we present a reconstruction approach to obtain MFP spectra from variable-length thermal conductivity measurements. We use a recently reported analytical solution of the BTE, along with efficient numerical simulations, to identify a suppression function that describes the discrepancy between the actual heat flux and that predicted by Fourier's law for a finite length domain. The MFP spectrum is then obtained using the convex optimization method described in Ref. 20. We apply this approach to SiGe nanowires and graphene ribbons. The measurements on graphene ribbons imply that MFPs are exceedingly long, with 70% of the heat being carried by phonons with MFPs longer than 1 micron.

Theory

Our goal is to relate experimentally measured thermal conductivities to the MFP spectrum, or the accumulated thermal conductivity as a fUnction of MFP14. Following the approach of Ref. 20, we therefore seek an equation of the form:

S(Kn„)f (Av)dAv ~

L-1K (Kn„)F(A„)dA„ (1)

where Kn, = A,/L is the Knudsen Number, A, is the MFP, L is the sample length along the direction of the temperature gradient, k denotes thermal conductivity as a function of length L, fA,) and F(A,) are differential and accumulative MFP spectra related by

F(Av)~

f (A)dA, and S is the heat flux suppression function

that equals the ratio of actual heat flux to the Fourier's law prediction. The kernel K is defined as K(Knv) = — dS/dKn,.

The inputs to this equation are a finite number of measured thermal conductivities k as a function of lengths L. To close the problem, we must identify the suppression function by solving the BTE, given by22:

+Wr Lt t,.

where ev is the desired distribution function, v is the angular frequency, eV is the distribution function at the equilibrium state, v is the group velocity of phonons, and tv is the relaxation time of phonons at certain frequency.

We obtain this function using two distinct approaches: a semi-analytic method and a numerical Monte Carlo (MC) method. First, we use a recently reported semi-analytical solution for steady heat

conduction through a crystal of thickness L with two blackbody boundaries23. In this solution, the BTE is linearized and solved using a series expansion method. The full details are given in Ref. 24. The final result for the suppression function and kernel are:

S(Kn„)~1 + 3Kn„ £s(Kn-1)- 1 (3)

K(Kn")~-dKnv~4-3£5(Kn-1)-Knv£4(Kn--1) (4) where En(x) is the exponential integral function, given by:

En(x) ~ mn-2exp[ - x )dm25.

Jo V mJ

This equation was derived by neglecting temperature slip between the black walls. Physically, this assumption is similar to the weakly quasiballistic regime described in Ref. 23 and implies that the ballistic phonons are low frequency modes with a small heat capacity. The assumption has been shown to be quite accurate for experimentally accessible length scales24.

We plot this result in Fig. 1a. For extremely short MFPs compared to the sample length L, the suppression function equals unity, indicating these phonons are diffusive and their heat flux contribution equals the Fourier's law prediction. As the Knudsen number increases, the suppression function decreases and eventually approaches zero, indicating that phonons contribute a smaller amount to the heat flux than predicted by Fourier's law. Physically, this suppression occurs because phonons cannot travel a full MFP before being absorbed by the blackbody boundary.

Figure 1 | (a) Suppression functions obtained from analytical (black solid line) and MC methods (open circles and open squares), and the Kernel (green solid line) from the analytic method. The two suppression functions are in excellent agreement. (b) MFP spectra for pure Si slabs of various thicknesses obtained from MC. Larger suppression for long MFP phonons occurs as the length of the Si slab decreases. (c) Thermal conductivities of pure Si slabs as a function of thickness calculated by MC. (d) Reconstructed (red circles) and the actual (blue solid line) MFP spectrum of bulk Si. All the values are normalized to bulk thermal conductivity of pure Si. The reconstructed result is in excellent agreement with the actual MFP spectrum.

We additionally solve the BTE numerically to validate the calculations above as well as to consider more complex situations such as when boundary scattering occurs. For this calculation, we use a linearized deviational Monte Carlo method to solve the adjoint BTE as described by Peraud et al26. This technique solves BTE by simulating advection and scattering of particles that represent phonons traveling inside the simulation domain. Substantial reductions in computational cost are achieved through a number of simple changes to the original MC algorithm. First, the deviational algorithm simulates only the deviation from a known equilibrium Bose-Einstein distribution, thereby incorporating deterministic information and reducing the variance. Further, for small temperature differences, the collision term in the BTE can be linearized, allowing particles to be simulated completely independently and without spatial and temporal discretization27. Next, we use a variable local equilibrium temperature method that closely matches the steady-state temperature profile.

Finally, we solve the adjoint BTE rather than the traditional BTE28,29. In the original algorithms of Refs. 26, 27, the probability for a certain phonon mode to be sampled is proportional to the density of states. Therefore, low frequency phonons are rarely sampled even though they contribute substantially to thermal conductivity, leading to large stochastic noise. The adjoint method overcomes this limitation by drawing particles with equal probability among all pho-non modes and correcting the bias introduced by this sampling when thermal properties are calculated. With these advances in numerical approach, we are able to solve the BTE in a 100 micron long domain in minutes on a desktop computer. Further, this numerical approach can incorporate boundary scattering mechanisms for arbitrary geometries exactly, unlike the analytical treatment.

To validate the code, we calculate the MFP spectrum for an infinite planar slab with two blackbody boundary conditions, the same problem solved by the semi-analytical method. For this calculation, we use an isotropically averaged dispersion for Si. The original dispersion and relaxation times were calculated by density functional theory (DFT) by Jesus Carrete and N. Mingo with ShengBTE30,31 and Phonopy32, from interatomic force constants obtained with VASP33-36. We reduce computational cost by taking advantage of the cubic symmetry of Si and computing an isotropic equivalent dispersion as described in Ref. 37. Using this dispersion, we calculate the MFP spectrum for variable lengths, as in Fig. 1b. We observe that decreasing the length of the domain results in the suppression of long MFP phonons to thermal conductivity compared to the bulk spectrum. The ratio of the differential MFP spectrum for a finite length to that for an infinite length yields the suppression function and is plotted for two lengths in Fig. 1a, demonstrating that the function obtained from our numerical approach exactly agrees with the analytical result.

With these tools, we now demonstrate the principal result of this work by using the suppression function to reconstruct the MFP spectrum from variable-length thermal conductivity measurements on a Si slab. We synthesized thermal conductivities as a function of length as shown in Fig. 1c. With these length dependent thermal conductivities and the suppression function, we used the same convex optimization method introduced in Ref. 20 to reconstruct the MFP spectrum. As in Figure 1d, the reconstructed accumulative MFP distribution is in excellent agreement with the actual accumulative MFP distribution, which is obtained from the DFT calculation, demonstrating that our approach can accurately reconstruct the MFP spectrum from length-dependent thermal conductivities of Si slabs.

We now numerically demonstrate that our approach can be applied to more general problems than Si slabs with a single scattering mechanism. We consider a SiGe nanowire in which point defects and the nanowire boundaries scatter phonons in addition to the intrinsic phonon-phonon scattering mechanism. The phonon-phonon relaxation times are taken to be the same as those of pure

Si30 36, while the mass defect scattering rate is given by t- 1 ~ x(1-x)A®438, where A is a constant of 3.01 x 10241 s3 for Si12xGex, which is obtained from Ref. 38. This model predicts a thermal conductivity of 14 W/mK for bulk Si09Ge01, which is consistent with other models by DFT calculation and experimental result38,39. These scattering rates are combined using Matthiessen's rule. We incorporate boundary scattering by explicitly simulating phonon trajectories inside a nanowire with a square cross-section of size 100 nm by 100 nm. We use Ziman's specularity parameter, p = exp( — 16p2ff2/l2), to determine the probability of specular or diffuse scattering, where s is surface roughness and l is the phonon wavelength40.

With this framework, we use our MC simulations to calculate the length-dependent thermal conductivities for a Si0.9Ge0.1 nanowire with surface roughness of s = 0.1 nm, a 100 nm by 100 nm square cross-section, over lengths from L = 5 nm to 16 mm as would be obtained in an experiment. Then, using only our knowledge ofthese thermal conductivities and the suppression function, we perform the reconstruction procedure to obtain the MFP spectrum of the nanowire. This result is shown in Fig. 2a. The reconstructed spectrum is in good agreement with the actual one without requiring any knowledge of the scattering mechanisms in the nanowire. Thus, the success of the reconstruction of MFP spectrum of the nanowire demonstrates the self-consistency of our approach.

Discussion

We now use our approach to examine two recent reports of length-dependent thermal conductivities in nanostructures. First, we consider SiGe alloy nanowires as investigated by Hsiao et al16. These nanowires were reported to have ballistic heat conduction persisting over approximately 8 microns. To investigate this experimental report, we calculate the length-dependent thermal conductivities of Si0.9Ge0.1 nanowires, which has the approximately the same cross-sectional area, 100 nm by 100 nm cross section, as that of the nanowires in Ref. 16, whose diameters range from 50 nm to 180 nm. Simulated data are shown in Fig. 2b for nanowires with extremely rough boundaries, s r ', and smooth boundaries, with s = 0.1 nm.

RMS=0.1 nm —RMS=<

MFP (nm)

Figure 2 | (a) The actual (MC data) and reconstructed MFP spectra of a simulated Si09Ge01 nanowire with surface roughness of RMS = 0.1 nm and the same square cross-section with a side length of 100 nm. The reconstructed distribution is in good agreement with the actual distribution, even though it is extracted from merely a series of discrete thermal conductivities without any boundary scattering information. (b)Thermal conductivities of Si09Ge01 nanowire from both MC simulations (lines) and experimental measurements16 (open squares) as a function of nanowire lengths. All simulated nanowires have the same square cross section as that in (a), and were simulated with very smooth surface (s = 0.1 nm, green dashed line) and extremely rough surface (s = infinity, red solid line), respectively. The experimental data does not follow the trend predicted by the simulations. All the thermal conductivities are normalized to their ''bulk'' value, which is thermal conductivity of the infinitely long nanowire.

www.nature.com/scientificreports \ > _

Knudsen Number A/L MFP or Length (nm)

Figure 3 | (a) Suppression function in both 3D (dashed black line) and 2D (solid blue line) spaces and their corresponding kernel functions. The 2D suppression function is very similar to the 3D one. (b) Experimentally measured length dependent thermal conductivity15 (blue open squares) and the corresponding reconstructed accumulative thermal conductivity as a function of phonon MFP (red open circles) in suspended graphene samples. All these thermal conductivities are normalized to thermal conductivities of''bulk'' graphene flakes, which is calculated using the extrapolating method in Ref. 21. Phonons with MFPs longer than 1 mm carry the majority of the heat in suspended graphene.

We observe a discrepancy between the trends of experimental data and our simulations. The experimental data suggest that the thermal conductivities of these SiGe nanowires are mainly due to phonons with a narrow MFP spectrum around ~8.3 mm, but our simulations indicate that even for very smooth nanowires (a = 0.1 nm), phonons within this range only contribute ~15% of the total thermal conductivity (Fig. 2a). Most of the heat is carried by MFPs less than 1 micron, with some heat being carried by longer MFPs in the partially specular case. Additionally, most actual nanowires have surface roughness higher than 0.1 nm41,42, and in this case long MFP phonons contribute only a small amount to heat conduction. Due to the relatively large diameter of the nanowire, changes to the phonon dispersion for thermal phonons due to phonon confinement are unlikely43. The experimental measurements thus do not agree with our self-consistent calculations. Due to the lack of appropriate experimental reports, we are unable to apply our approach to other nanowire data sets. Further experimental investigation is necessary to address this discrepancy.

Next, we consider recent measurements on graphene15. In this work, Xu et al15 performed thermal conductivity measurements over variable lengths on suspended single-layer graphene ribbons to infer an average MFP of 240 nm at room temperature. Using our approach, we can use these same measurements to obtain the MFP spectrum of graphene. Due to the difficulty of fabrication for very long suspended graphene devices and vulnerability of these devices during measurement, the authors of Ref. 15 didn't obtain the saturated thermal conductivity from an extremely long, or "bulk", gra-phene sample. Therefore, we evaluated the saturation value as ~2,000 W/mK with the same extrapolating method in Ref. 21, which is in the range of previous reported experimental results44-46.

There are two subtleties that require discussion before applying our reconstruction approach to graphene. First, our derivation is based on the relaxation time approximation (RTA) of the BTE. It is well known that the RTA with the computed relaxation times from DFT under-predicts the thermal conductivity of graphene due to the importance of normal processes47. However, while our derivation is based on the RTA, we do not make any assumption of the values of the relaxation times, but rather only that an effective relaxation time for each phonon mode can be identified. Our approach can be applied to graphene provided that we regard the MFP variable as an effective MFP that represents the average propagation length for a particular phonon frequency as determined by both normal and Umklapp processes.

Second, we have derived our suppression function for an isotropic material with a three-dimensional phase space. While graphene can be reasonably modeled as isotropic for the in-plane directions, the

phase space is two-dimensional. This dimensionality change requires a modification of the form of the BTE for 2D materials. Repeating the derivation in Ref. 24 for a 2D phase space yields the suppression function S and kernel K as:

S2D(Kn„) = 1 + - Kn, p

cos3 (h)e K-»~(h>d0--

K2D(Kn„)~-

dS2D dKn,

cos3(h)e-

_ 8 ~ 3p

cos2(h)e-

Figure 3a plots the 2D and 3D suppression functions and kernels, demonstrating that the two are quite similar. Using the 2D kernel, we apply our method to obtain the MFP spectrum of a graphene ribbon as in Fig. 3b. This result shows that MFPs in graphene span a large range from ~ 100 nm to ~ 10 mm. In addition, we observe that a large portion of heat in graphene is carried by long MFP phonons: ~70% of thermal conductivity are from phonons with MFPs greater than 1 mm. Further, as reported in Ref. 15, the widths of these graphene ribbons are only between 2 to 4 microns. From Fig. 3b, phonons with MFPs longer than 4 mm still carry 13% of heat, which means that some of these long MFP phonons must be specularly reflected at the edges of the graphene ribbons. This observation further confirms the report in Ref. 15 of weak width-dependent thermal conductivities of suspended graphene ribbons when widths are larger than 1.5 mm. Our MFP reconstruction approach has thus provided valuable insights into the intrinsic and edge scattering mechanisms in graphene ribbons that are difficult to obtain from knowledge of only the average MFP.

Summary

We have presented a reconstruction method that allows MFP spectra of nanostructures to be obtained from length-dependent thermal conductivity measurements. Our approach requires no prior knowledge of the scattering mechanisms in the nanostructure. By applying our approach to recent measurements on graphene ribbons, we find that more than half of the heat in graphene is contributed by phonons with MFPs exceeding 1 micron.

1. Zebarjadi, M., Esfarjani, K., Dresselhaus, M. S., Ren, Z. F. & Chen, G. Perspectives on thermoelectrics: from fundamentals to device applications. Energy Environ. Sci. 5, 5147-5162 (2012).

Kn, cos

Kn, cos

2. Cahill, D. G. etal. Nanoscale thermal transport. ii. 2003-2012. Appl. Phys. Rev. 1, 011305 (2014).

3. Zhao, L.-D., Dravid, V. P. & Kanatzidis, M. G. The panoscopic approach to high performance thermoelectrics. Energy Environ. Sci. 7, 251-268 (2014).

4. Li, D. etal. Thermal conductivity of individual silicon nanowires. Appl. Phys. Lett. 83, 2934-2936 (2003).

5. Hochbaum, A. I. et al. Enhanced thermoelectric performance of rough silicon nanowires. Nature 451, 163-167 (2008).

6. Boukai, A. I. et al. Silicon nanowires as efficient thermoelectric materials. Nature 451, 168-171 (2008).

7. Chang, C. W., Okawa, D., Garcia, H., Majumdar, A. & Zettl, A. Breakdown of fourier's lawin nanotube thermal conductors. Phys. Rev. Lett. 101,075903 (2008).

8. Johnson, J. A. et al. Direct measurement of room-temperature nondiffusive thermal transport over micron distances in a silicon membrane. Phys. Rev. Lett. 110, 025901 (2013).

9. Tighe, T. S., Worlock, J. M. & Roukes, M. L. Direct thermal conductance measurements on suspended monocrystalline nanostructures. Appl. Phys. Lett. 70, 2687-2689 (1997).

10. Poudel, B. et al. High-thermoelectric performance of nanostructured bismuth antimony telluride bulk alloys. Science 320, 634-638 (2008).

11. Biswas, K. etal. High-performance bulk thermoelectrics with all-scale hierarchical architectures. Nature 489, 414-418 (2012).

12. Chowdhury, I. et al. On-chip cooling by superlattice-based thin-film thermoelectrics. Nat. Nanotechnol 4, 235-238 (2009).

13. Tian, Z., Lee, S. & Chen, G. Heat transfer in thermoelectric materials and devices. J. Heat Transfer 135, 061605-061605 (2013).

14. Chen, G. & Dames, C. Thermal conductivity of nanostructured thermoelectric materials. In Thermoelectrics Handbook: Macro to Nano, 42-1-42-16- (CRC Press, 2005).

15. Xu, X. et al. Length-dependent thermal conductivity in suspended single-layer graphene. Nat. Commun. 5, 1-6 (2014).

16. Hsiao, T.-K. et al. Observation of room-temperature ballistic thermal conduction persisting over 8.3 mm in sige nanowires. Nat. Nanotechnol. 8, 534-538 (2013).

17. Esfarjani, K., Chen, G. & Stokes, H. T. Heat transport in silicon from first-principles calculations. Phys. Rev. B 84, 085204 (2011).

18. Henry, A. S. & Chen, G. Spectral phonon transport properties ofsilicon based on molecular dynamics simulations and lattice dynamics. J. Comput. Theor. Nanos. 5, 141-152 (2008).

19. Minnich, A. J. et al. Thermal conductivity spectroscopy technique to measure phonon mean free paths. Phys. Rev. Lett. 107, 095901 (2011).

20. Minnich, A. J. Determining phonon mean free paths from observations of quasiballistic thermal transport. Phys. Rev. Lett. 109, 205901 (2012).

21. Wei, Z. etal. Phonon mean free path of graphite along the c-axis. Appl. Phys. Lett. 104, 081903 (2014).

22. Majumdar, A. Microscale heat conduction in dielectric thin films. J. Heat Transfer 115, 7-16 (1993).

23. Hua, C. & Minnich, A. J. Transport regimes in quasiballistic heat conduction. Phys. Rev. B 89, 094302 (2014).

24. Hua, C. & Minnich, A. J. Cross-plane heat conduction in thin solid films. ArXiv e-prints (2014)Http://arxiv.org/abs/1410.2845, 1410.2845

25. Chen, G. Nanoscale energy transport and conversion: a parallel treatment of electrons, molecules, phonons, and photons (Oxford University Press, Oxford; New York, 2005).

26. Peraud, J.-P. M. & Hadjiconstantinou, N. G. Efficient simulation of multidimensional phonon transport using energy-based variance-reduced monte carlo formulations. Phys. Rev. B 84, 205331 (2011).

27. Peraud, J.-P. M. & Hadjiconstantinou, N. G. An alternative approach to efficient simulation of micro/nanoscale phonon transport. Appl. Phys. Lett. 101, 153114 (2012).

28. Peraud, J.-P. M., Landon, C. D. & Hadjiconstantinou, N. G. Monte carlo methods for solving the boltzmann transport equation. Annual Review ofHeat Transfer 17,(2014).

29. Peraud, J.-P. M., Landon, C. D. & Hadjiconstantinou, N. G. Deviational methods for small-scale phonon transport. Mech. Eng. Rev. 1, FE0013-FE0013 (2014).

30. Li, W., Carrete, J. C., Katcho, N. A. & Mingo, N. Shengbte: A solver of the boltzmann transport equation for phonons. Comput. Phys. Commun. 185, 1747-1758 (2014).

31. Li, W., Carrete, J., Katcho, N. A. & Mingo, N. Shengbte. www.shengbte.org, Date of

access: 01/07/2014

32. Togo, A. Phonopy. Http://phonopy.sourceforge.net/, Date of access: 01/07/2014.

33. Kresse, G. &Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev.B

47, 558-561 (1993).

34. Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metalamorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251-14269 (1994).

35. Kresse, G. & Furthmller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15-50 (1996).

36. Kresse, G. & Furthmiiller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169-11186 (1996).

37. Hua, C. & Minnich, A. J. Analytical green's function of the multidimensional frequency-dependent phonon boltzmann equation. Phys. Rev. B 90, 214306 (2014).

38. Wang, Z. & Mingo, N. Diameter dependence of sige nanowire thermal conductivity. Appl. Phys. Lett. 97, 101903 (2010).

39. Garg, J., Bonini, N., Kozinsky, B. & Marzari, N. Role of disorder and anharmonicity in the thermal conductivity of silicon-germanium alloys: A first-principles study. Phys. Rev. Lett. 106, 045901 (2011).

40. Ziman, J. Electrons and Phonons: The Theory of Transport Phenomena in Solids International series of monographs on physics (Clarendon Press, 1962).

41. Azeredo, B. P. et al. Silicon nanowires with controlled sidewall profile and roughness fabricated by thin-film dewetting and metal-assisted chemical etching. Nanotechnology 24, 225305 (2013).

42. Feser, J. P. et al. Thermal conductivity of silicon nanowire arrays with controlled roughness. J. Appl. Phys. 112, 114306 (2012).

43. Turney, J. E., McGaughey, A. J. H. & Amon, C. H. In-plane phonon transport in thin films. J. Appl. Phys. 107, 024317 (2010).

44. Cai, W. etal. Thermal transport in suspended and supported monolayer graphene grown by chemical vapor deposition. Nano Lett. 10, 1645-1651 (2010).

45. Chen, S. et al. Raman measurements of thermal transport in suspended monolayer graphene of variable sizes in vacuum and gaseous environments. ACS Nano 5, 321-328 (2011).

46. Chen, S. et al. Thermal conductivity of isotopically modified graphene. Nat. Mater. 11,203-207 (2012).

47. Lindsay, L., Broido, D. A. & Mingo, N. Flexural phonons and thermal transport in graphene. Phys. Rev. B 82, 115427 (2010).

Acknowledgments

This work was supported by a start-up fund from the California Institute of Technology and by the National Science Foundation under CAREER Grant CBET 1254213.

Author contributions

H.Z., C.H. and D.D. performed simulations and calculations. H.Z. and A.M. analyzed data. H.Z., A.M., C.H. and D.D. discussed the result. H.Z. and A.M. wrote the main manuscript text.

Additional information

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Zhang, H., Hua, C., Ding, D. & Minnich, A.J. Length Dependent Thermal Conductivity Measurements Yield Phonon Mean Free Path Spectra in Nanostructures. Sci. Rep. 5, 9121; D0I:10.1038/srep09121 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder in order to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/