Open AccessCommunication

Phonon Transport through Nanoscale Contact in Tip-Based Thermal Analysis of Nanomaterials

Jay Dulhani

and

Bong Jae Lee

Department of Mechanical Engineering, Korea Advanced Institute of Science and Technology, Daejeon 34141, Korea

Author to whom correspondence should be addressed.

Nanomaterials 2017, 7(8), 200; https://doi.org/10.3390/nano7080200

Submission received: 18 June 2017 / Revised: 17 July 2017 / Accepted: 21 July 2017 / Published: 28 July 2017

(This article belongs to the Special Issue Nanomaterials for Renewable and Sustainable Energy)

Download

Browse Figures

Versions Notes

Abstract

Nanomaterials have been actively employed in various applications for energy and sustainability, such as biosensing, gas sensing, solar thermal energy conversion, passive radiative cooling, etc. Understanding thermal transports inside such nanomaterials is crucial for optimizing their performance for different applications. In order to probe the thermal transport inside nanomaterials or nanostructures, tip-based nanoscale thermometry has often been employed. It has been well known that phonon transport in nanometer scale is fundamentally different from that occurred in macroscale. Therefore, Fourier’s law that relies on the diffusion approximation is not ideally suitable for describing the phonon transport occurred in nanostructures and/or through nanoscale contact. In the present study, the gray Boltzmann transport equation (BTE) is numerically solved using finite volume method. Based on the gray BTE, phonon transport through the constriction formed by a probe itself as well as the nanoscale contact between the probe tip and the specimen is investigated. The interaction of a probe and a specimen (i.e., treated as a substrate) is explored qualitatively by analyzing the temperature variation in the tip-substrate configuration. Besides, each contribution of a probe tip, tip-substrate interface, and a substrate to the thermal resistance are analyzed for wide ranges of the constriction ratio of the probe.

Keywords:

nanoscale constriction and contact; Boltzmann transport equation; phonon transport

Graphical Abstract

1. Introduction

Recently, we are witnessing extreme miniaturization of thermal devices and other components [1,2,3]. Besides, rapid development in nanotechnology enables that nanostructures or nanomaterials have been widely employed in various applications for energy and sustainability, such as biosensing, gas sensing, solar thermal energy conversion, and passive radiative cooling [4,5,6,7]. Broad application of nanomaterials or nanostructures is posing new challenges for the thermal science community. One such challenge is in fundamental understanding of the geometry-induced effects on phonon transport in nanostructures and/or nanoscale contacts, which is crucial for optimizing the thermal performance of nanostructures.

Development and continuous evolution of a tip-based thermal microscopy are revolutionizing our perception about thermal transport at sub-micron scale. In a thermal microscope system, a nanometer-sized probe is often employed to interact with a specimen for measuring local temperature. Such tip-substrate configuration also has enormous practical applications, such as in data writing [8], chemical and biomolecular applications [9], nanofabrication [10] and nanoscale heat transfer analysis [11]. These applications can be broadly classified into two major area: nanomanufacturing and nanometrology. For the former application, the temperature of the probe is raised; thus, the probe acts as a localized heat source while interacting with a specimen (or a substrate). We are particularly interested in understanding the tip-substrate interaction in such systems.

It is known that the thermal conductivity of nanostructures is significantly lower than its bulk counterpart because of the phonon boundary scattering. Recent research also suggested that it is possible to alter the transport regime and thus the transport property by tailoring the geometry of nanostructures [12,13]. Besides, Jean et al. [14] concluded that the shape of the probe tip has a reasonable effect on probe’s measurement capacity. Thus, fundamental understanding of temperature distribution in the probe-substrate system is essential for the accurate use of tip-based thermal analysis.

Another area of our interest is in the quantification of thermal resistance for a tip-based thermal analysis system in which a probe with a nanoscale tip is in contact with a substrate. For macroscale thermal transport, the thermal contact resistance across the interface has been subjects of multiple studies [15,16]. As far as nanoscale thermal transport is concerned, the difference in the vibrational properties of phonon and partial transmission of phonon through the interface give rise to the thermal boundary resistance (TBR). Initial research has shown that thermal transport between the nanostructures can be limited by TBR [17]. Attempts are being made to model the thermal boundary resistance theoretically. Cooper et al. [15] derived the constriction resistance of a point contact between two semi-infinite objects of the same material using diffusion equation. This expression was strictly derived for the diffusive regime and will lead to significant error if transport is in quasi-ballistic regime. Taking it forward, Prasher [18] developed an analytical model to calculate the thermal resistance of nanoscale constrictions.

In light of the preceding discussion, the objective of present work is twofold: (i) fundamental understanding of phonon transport between a probe with nanoscale tip and a specimen (treated as a substrate made of the same material with the probe for simplicity) and (ii) quantification of thermal resistance of such geometry. These will be achieved by solving the phonon Boltzmann transport equation. In a recent study [13], the authors have used finite volume method to investigate the effect of constriction on phonon transport in thin films and nanowires. A similar approach is used here but for the different purpose.

2. Results and Discussion

Being aware of the complexity of the configuration under consideration, following assumptions are made. First, we assumed that perfect contact exists between the probe tip and the substrate; that is, there is no boundary scattering and/or resampling at the tip-substrate interface. Secondly, it is assumed that experiments occur in the vacuum; thus, convection loss through the air is absent.

The schematic of geometric configuration considered in the present study is shown in Figure 1. Axisymmetric model is used here as it closely mimics the probe used in atomic force microscope. In Figure 1, length of the tip is

L_{p}

, and the radius of broad end of tip is

r_{m a x}

, whereas the radius of narrow end of tip is

r_{m i n}

. The length of substrate is

L_{s u b}

and the radius of substrate is

r_{s u b}

. Broad end of the tip is maintained at temperature

T_{H}

, whereas, the bottom of the substrate is maintained at temperature

T_{L}

At steady state, the gray Boltzmann transport equation (BTE) can be written for an axisymmetric system with the relaxation time approximation as [19]:

\frac{μ}{r} \frac{\partial (r e^{″})}{\partial r} - \frac{1}{r} \frac{\partial (η e^{″})}{\partial ϕ} + ξ \frac{\partial e^{″}}{\partial z} = \frac{e^{0} - e^{″}}{Λ}

(1)

where

e^{″}

is the net phonon energy density,

e^{0}

is the equilibrium energy density,

μ

η

ξ

are direction cosines, and

Λ

is the phonon mean free path. Details about formulation and boundary conditions can be found elsewhere [13,19,20,21]; thus, they are not repeated here.

All the results presented here are for Silicon at 300 K. Thermophysical properties required for the BTE simulations are calculated by phonon branch averaging procedure [22]. Silicon’s properties are group velocity

v_{g} = 1804

m/s, mean free path

Λ = 260.4

nm, and specific heat

c_{p} = 0.93 \times 10^{6}

J/m

^{3} \cdot

K. For all results presented here, angular domain of

2 π

(i.e.,

0 \leq θ \leq π

and

0 \leq ϕ \leq π

due to symmetry) is discretized into

128 \times 128

non-overlapping control angles. Further refinement in it does not bring any reasonable change to the computational accuracy. The geometric dimensions are

L_{s u b} = 3 Λ

and

r_{s u b} = 4 Λ

for the substrate (refer to Figure 1). As far as the probe tip is concerned, its length is fixed to be

L_{p} = 1500

nm and radius of the broad end varies as

r_{m a x} = 100

, 150, or 200 nm. Radius of the narrow end (i.e., size of the tip-substrate contact) is determined according to the constriction ratio,

C_{r} = r_{m i n} / r_{m a x}

Qualitative features of the tip-substrate interaction are explored first by analyzing the temperature distribution. Figure 2 represents the non-dimensional temperature distribution along the z-axis of the tip-substrate system when

r_{m a x} = 100

nm. In the figure, temperature is normalized as

Θ = (T - T_{L}) / (T_{H} - T_{L})

, and the z-axis is normalized as

z^{*} = z / (L_{p} + L_{s u b})

. For the case of

C_{r} = 1

, a constant temperature gradient is observed on the probe side, which implies that the phonon transport is mainly in diffusive regime. Substrate’s response is similar to the case when the substrate is subjected to a localized surface heating, which is not surprising as in the considered configuration the tip-substrate interface acts as a localized heat source for the substrate. With increase in the degree of constriction (i.e., lower

C_{r}

values), phonon transport in the probe side becomes increasingly ballistic due to the constriction effect (i.e., phonon boundary scattering), as evident from nearly constant temperature profile in the probe side for

C_{r} = 0.15

(i.e.,

r_{m i n} = 15

nm). Moreover, irrespective of the constriction ratio, the temperature of the substrate is found to be roughly constant. This is due to high thermal resistance across the tip-substrate interface. Another interesting observation is a drastic change in the temperature gradient across the tip-substrate interface. Such trend points to the non-equilibrium nature of phonon transport across the nanoscale tip-substrate contact even though perfect contact is assumed.

Next, we want to explore the effect of absolute dimensions on the qualitative features of thermal transport. To achieve this, constriction ratio of the probe tip is kept constant while the absolute dimension (

r_{m a x}

and

r_{m i n}

) is varied. Figure 3 shows the temperature distribution along the z-axis of the tip-substrate system for different

r_{m a x}

values. As can be seen from Figure 3, if

C_{r}

is kept constant, variation in

r_{m a x}

(or in

r_{m i n}

) does not affect the temperature distribution significantly. Thus, it would be inferred that for the range of length scale considered, the constriction ratio and not the absolute dimension is the critical parameter for determining thermal transport inside the tip-substrate system. Similar observations were also reported by Cheney [23]. It should be also noted that size confinement could alter the mean free path of the low-frequency phonons more than that of the high-frequency phonons [24]. To properly model the spectrally dependent phonon mean free path, the phonon dispersion relation as well as the three phonon scattering relaxation time should be considered, which is beyond the scope of this study.

The study of the thermal resistance is critical as it can significantly affect the phonon transport in nanostructures. The heat transfer rate Q from

T_{H}

T_{L}

is obtained directly from the BTE simulations. Various components of the thermal resistance are defined as: total thermal resistance of the tip-substrate system

R_{t o t a l} = (T_{H} - T_{L}) / Q

, thermal resistance in the substrate

R_{s u b} = Δ T_{s u b} / Q

, thermal resistance at the interface between the tip and the substrate

R_{i n t e r} = Δ T_{i n t e r} / Q

, and thermal resistance in the probe tip

R_{t i p} = R_{t o t a l} - R_{s u b} - R_{i n t e r}

. Here, the temperature profile in the substrate is extrapolated linearly to calculate

Δ T_{s u b}

, and

Δ T_{i n t e r}

is calculated by averaging the temperature of control volumes adjacent to the interface and then taking their difference. Care has been taken to ensure that the cell adjacent to the interface are at equal distance from the tip-substrate interface. The similar procedure was also used by Jean et al. [14].

Figure 4 shows variation of the thermal resistance with constriction ratio for

r_{m a x} = 100

nm. As we observed from the temperature profile, the effect of the substrate is broadly restricted by high thermal resistance at the interface. In fact,

R_{s u b}

is less than 2% of the total resistance and thus can be neglected. Besides,

(R_{t i p} + R_{i n t e r})

is much higher than the total resistance predicted by Fourier’s law, and this difference increases as

C_{r}

decreases (i.e., the degree of constriction increases). Fourier’s law is found to substantially underpredict the thermal resistance as compared to that predicted by the BTE. In Figure 4,

R_{t i p}

and

R_{i n t e r}

increase with increase in the degree of constriction. Interestingly, when observed from non-dimensional basis,

C_{r}

has more profound impact on

R_{i n t e r}

than on

R_{t i p}

. Please note that in reality, the probe tip cannot make a perfect contact with the substrate. In other words, there must be additional boundary scattering at the tip-substrate interface due to lattice mismatch. Therefore,

R_{i n t e r}

estimated in this work indicates the minimum possible value.

In general, phonon boundary scattering depends on the surface area to volume ratio. In this regard, a two-dimensional (2-D) model of the BTE may differently predict when compared with the axisymmetric model. Here, thermal transport in the tip-substrate system has been analyzed using both the axisymmetric model and the 2-D model. The Gray BTE code utilized for the 2-D model is same as used previously by the authors [13]. Figure 5 shows the temperature distribution along the z-axis of the tip-substrate system for

r_{m a x} = 100

nm. The qualitative comparison indicates that constriction can play a dominating role in governing the thermal transport especially when the degree of constriction is high (i.e., when

C_{r}

value is low). For a given heat flux, the thermal resistance of different components is proportional to the temperature drop across the corresponding component. It can be noticed that the interface provides significantly higher thermal resistance to phonons in the axisymmetric model. A possible reason behind such trend is that the area to volume ratio is greater for the axisymmetric case.

3. Summary

In summary, we have investigated the phonon transport through the constriction formed by a probe itself as well as the nanoscale contact between the probe tip and the specimen based on the gray Boltzmann transport equation. We observed that constriction of the probe tip could play a dominating role in governing the phonon transport especially when the degree of constriction is high. It was also found that the combined thermal resistance of the probe tip and the tip-substrate interface predicted by the BTE is much greater than the total resistance predicted by Fourier’s law and this difference increases as the degree of constriction increases. Lastly, the tip-substrate interface is found to offer significantly higher resistance to phonons in the axisymmetric model as compared to the 2-D model. The insights obtained from the present study may be applied to improve our understanding of a tip-based device.

Acknowledgments

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and future Planning (NRF-2015R1A2A1A10055060).

Author Contributions

Jay Dulhani developed numerical code and conducted simulation; Jay Dulhani and Bong Jae Lee analyzed the data and wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

Cahill, D.G.; Ford, W.K.; Goodson, K.E.; Mahan, G.D.; Majumdar, A.; Maris, H.J.; Merlin, R.; Phillpot, S.R. Nanoscale Thermal Transport. J. Appl. Phys. 2003, 93, 793–818. [Google Scholar] [CrossRef]
Cahill, D.G.; Braun, P.V.; Chen, G.; Clarke, D.R.; Fan, S.; Goodson, K.E.; Keblinski, P.; King, W.P.; Mahan, G.D.; Majumdar, A.; et al. Nanoscale Thermal Transport. II. 2003–2012. Appl. Phys. Rev. 2014, 1, 11305. [Google Scholar] [CrossRef]
Shi, L.; Dames, C.; Lukes, J.R.; Reddy, P.; Duda, J.; Cahill, D.G.; Lee, J.; Marconnet, A.; Goodson, K.E.; Bahk, J.H.; et al. Evaluating Broader Impacts of Nanoscale Thermal Transport Research. Nanoscale Microscale Thermophys. Eng. 2015, 19, 127–165. [Google Scholar] [CrossRef]
Yun, J.; Ahn, J.H.; Lee, B.J.; Moon, D.I.; Choi, Y.K.; Park, I. Temperature Measurement of Joule Heated Silicon Micro/Nanowires Using Selectively Decorated Quantum Dots. Nanotechnology 2016, 27, 505705. [Google Scholar] [CrossRef] [PubMed]
Jeon, J.; Park, S.; Lee, B.J. Optical Property of Blended Plasmonic Nanofluid Based on Gold Nanorods. Opt. Express 2014, 22, A1101–A1111. [Google Scholar] [CrossRef] [PubMed]
Han, S.; Shin, J.H.; Jung, P.H.; Lee, H.; Lee, B.J. Broadband Solar Thermal Absorber Based on Optical Metamaterials for High-Temperature Applications. Adv. Opt. Mater. 2016, 4, 1265–1273. [Google Scholar] [CrossRef]
Zhai, Y.; Ma, Y.; David, S.N.; Zhao, D.; Lou, R.; Tan, G.; Yang, R.; Yin, X. Scalable-Manufactured Randomized Glass-Polymer Hybrid Metamaterial for Daytime Radiative Cooling. Science 2017, 355, 1062–1066. [Google Scholar] [CrossRef] [PubMed]
Binnig, G.; Despont, M.; Drechsler, U.; Haeberle, W.; Lutwyche, M.; Vettiger, P.; Mamin, H.; Chui, B.; Kenny, T.W. Ultrahigh-Density Atomic Force Microscopy Data Storage with Erase Capability. Appl. Phys. Lett. 1999, 74, 1329–1331. [Google Scholar] [CrossRef]
Poggi, M.A.; Lillehei, P.T.; Bottomley, L.A. Chemical Force Microscopy on Single-Walled Carbon Nanotube Paper. Chem. Mater. 2005, 17, 4289–4295. [Google Scholar] [CrossRef]
King, W.P.; Bhatia, B.; Felts, J.R.; Kim, H.J.; Kwon, B.; Lee, B.; Somnath, S.; Rosenberger, M. Heated Atomic Force Microscope Cantilevers and Their Applications. Ann. Rev. Heat Transf. 2013, 16, 287–326. [Google Scholar] [CrossRef]
Kim, K.; Song, B.; Fernández-Hurtado, V.; Lee, W.; Jeong, W.; Cui, L.; Thompson, D.; Feist, J.; Reid, M.H.; García-Vidal, F.J.; et al. Radiative Heat Transfer in the Extreme Near Field. Nature 2015, 528, 387–391. [Google Scholar] [CrossRef] [PubMed]
Jean, V.; Fumeron, S.; Termentzidis, K.; Zianni, X.; Lacroix, D. Monte Carlo Simulations of Phonon Transport in Si Nanowires with Constrictions. Int. J. Heat Mass Transf. 2015, 86, 648–655. [Google Scholar] [CrossRef]
Dulhani, J.; Lee, B.J. Effect of Constriction on Phonon Transport in Silicon Thin Films and Nanowires. Smart Sci. 2016, 4, 173–179. [Google Scholar] [CrossRef]
Jean, V.; Randrianalisoa, J.; Trannoy, N. Monte Carlo Simulation of Phonon Transport Across Si-Si and SiO₂ Interfaces. In Proceedings of the 21st International Workshop on Thermal Investigations of ICs and Systems (THERMINIC), Paris, France, 30 September–2 October 2015; pp. 1–6. [Google Scholar]
Cooper, M.; Mikić, B.; Yovanovich, M. Thermal Contact Conductance. Int. J. Heat Mass Transf. 1969, 12, 279–300. [Google Scholar] [CrossRef]
Mikić, B. Thermal Contact Conductance; Theoretical Considerations. Int. J. Heat Mass Transf. 1974, 17, 205–214. [Google Scholar] [CrossRef]
Singh, D.; Murthy, J.Y.; Fisher, T.S. Phonon Transport Across Mesoscopic Constrictions. J. Heat Transf. 2011, 133, 42402. [Google Scholar] [CrossRef]
Prasher, R. Predicting the Thermal Resistance of Nanosized Constrictions. Nano Lett. 2005, 5, 2155–2159. [Google Scholar] [CrossRef] [PubMed]
Lemonnier, D. Solution of the Boltzmann Equation for Phonon Transport. In Microscale and Nanoscale Heat Transfer; Springer: Berlin, Germany, 2007; pp. 77–106. [Google Scholar]
Chai, J.C.; Parthasarathy, G.; Lee, H.S.; Patankar, S.V. Finite Volume Radiative Heat Transfer Procedure for Irregular Geometries. J. Thermophys. Heat Transf. 1995, 9, 410–415. [Google Scholar] [CrossRef]
Kim, M.Y.; Baek, S.W. Analysis of Radiative Transfer in Cylindrical Enclosures using the Finite Volume Method. J. Thermophys. Heat Transf. 1997, 11, 246–252. [Google Scholar] [CrossRef]
Chen, G. Thermal Conductivity and Ballistic-Phonon Transport in the Cross-Plane Direction of Superlattices. Phys. Rev. B 1998, 57, 14958. [Google Scholar] [CrossRef]
Cheney, D.A. Computational Modeling of Geometry Dependent Phonon Transport in Silicon Nanostructures. Ph.D. Thesis, University of Pennsylvania, Philadelphia, PA, USA, 2013. [Google Scholar]
Muñoz, R.M.; Abad, B.; Manzano, C.; Torres, P.; Cartoixà, X.; Alvarez, F.; Martín, G.M. Thermal Conductivity of Bi₂Te₃ Nanowires: How Size Affects Phonon Scattering. Nanoscale 2017, 9, 6741–6747. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Geometric configuration of the probe tip and the substrate.

Figure 2. Temperature distribution along the z-axis when

r_{m a x} = 100

nm for various

C_{r}

values.

Figure 2. Temperature distribution along the z-axis when

r_{m a x} = 100

nm for various

C_{r}

values.

Figure 3. Temperature distribution along the z-axis for constant constriction ratio of (a)

C_{r} = 0.4

and (b)

C_{r} = 0.7

, while varying

r_{m a x}

values.

Figure 3. Temperature distribution along the z-axis for constant constriction ratio of (a)

C_{r} = 0.4

and (b)

C_{r} = 0.7

, while varying

r_{m a x}

values.

Figure 4. Variation of the thermal resistance with respect to the constriction ratio when

r_{m a x} = 100

nm.

Figure 4. Variation of the thermal resistance with respect to the constriction ratio when

r_{m a x} = 100

nm.

Figure 5. Temperature distribution along the z-axis as predicted by the axisymmetric model and the 2-D model when

r_{m a x} = 100

nm. In the 2-D model, the width of broad end of tip is

2 \times r_{m a x}

, whereas the width of narrow end of tip is

2 \times r_{m i n}

Figure 5. Temperature distribution along the z-axis as predicted by the axisymmetric model and the 2-D model when

r_{m a x} = 100

nm. In the 2-D model, the width of broad end of tip is

2 \times r_{m a x}

, whereas the width of narrow end of tip is

2 \times r_{m i n}

© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Dulhani, J.; Lee, B.J. Phonon Transport through Nanoscale Contact in Tip-Based Thermal Analysis of Nanomaterials. Nanomaterials 2017, 7, 200. https://doi.org/10.3390/nano7080200

AMA Style

Dulhani J, Lee BJ. Phonon Transport through Nanoscale Contact in Tip-Based Thermal Analysis of Nanomaterials. Nanomaterials. 2017; 7(8):200. https://doi.org/10.3390/nano7080200

Chicago/Turabian Style

Dulhani, Jay, and Bong Jae Lee. 2017. "Phonon Transport through Nanoscale Contact in Tip-Based Thermal Analysis of Nanomaterials" Nanomaterials 7, no. 8: 200. https://doi.org/10.3390/nano7080200

APA Style

Dulhani, J., & Lee, B. J. (2017). Phonon Transport through Nanoscale Contact in Tip-Based Thermal Analysis of Nanomaterials. Nanomaterials, 7(8), 200. https://doi.org/10.3390/nano7080200

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Phonon Transport through Nanoscale Contact in Tip-Based Thermal Analysis of Nanomaterials

Abstract

1. Introduction

2. Results and Discussion

3. Summary

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI