A statistical theory of interfacial thermal conductivity
Interfacial thermal conductivity model based on statistical features of two faces in contact, discussing surface roughness and heat flow
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Interfacial thermal conductivity model based on statistical features of two faces in contact, discussing surface roughness and heat flow
Dynamical control of thermal transport at the nanoscale provides a time-domain strategy for optimizing thermal management in nanoelectronics, magnetic devices, and thermoelectric devices. However, the rate of change available for thermal switches and regulators is limited to millisecond time scales, calling for a faster modulation speed. Here, in this work, time-resolved X-ray diffraction measurements and thermal transport modeling reveal an ultrafast modulation of the interfacial thermal conductance of an FeRh/MgO heterostructure as a result of a structural phase transition driven by optical excitation. Within 90 ps after optical excitation, the interfacial thermal conductance is reduced by a factor of 5 and lasts for a few nanoseconds, in comparison to the value at the equilibrium FeRh/MgO interface. The experimental results combined with thermal transport calculations suggest that the reduced interfacial thermal conductance results from enhanced phonon scattering at the interface where the lattice experiences transient in-plane biaxial stress due to the structural phase transition of FeRh. Our results suggest that optically driven phase transitions can be utilized for ultrafast nanoscale thermal switches for device application.
Hot isostatic pressing (HIPing) was found to increase the thermal diffusivity/conductivity of uniaxial silicon carbide fiber-reinforced reaction-bonded silicon nitride (RBSN) matrix composites, as the result of the densification of the matrix, the increase in the grain size of the silicon carbide and the improved thermal contact between the fibers and the matrix. Transverse to the fiber direction the thermal diffusivity/conductivity was found to be a function of the surrounding gaseous atmosphere due to the access of the gas phase to the fiber-matrix interface, which was facilitated by the existence of an interfacial gap due to the thermal expansion mismatch between the fibers and the matrix. The interfacial conductance was found to exhibit a strong positive temperature dependence as the result of the closure of the interfacial gap with increasing temperature.
Graphene continues to demonstrate promise as a highly effective barrier coating, even at only one atom thick. The thermal properties of this coating are also promising to allow diffusion of heat across the surface, as the isolated graphene is an intrinsically good thermal conductor, yet this and its behavior under mechanical deformation has been less extensively studied. This report demonstrates that the in-plane thermal conductivity and interfacial thermal conductance of graphene coatings on copper are affected by mechanical strain. By inducing a strain in the copper substrate, the Raman-active 2D peak exhibits a change in position as well as a change in laser power dependence as the copper substrate is uniaxially elongated to a maximum of 0.5%. Non-linear trends in thermal conductivity are observed with tensile strain in samples with differing strain transfer rates from the substrate, indicating the close correlation between intrinsic thermal conduction and interfacial properties in atomically thin coatings transferred onto metals.
Thermal and electric conductivities for interface separating two parts of system with different temperature and chemical and electrical potential
Thermal conductivity in nanoparticle solids has been previously reported in the range of 0.1-1 W m -1 K -1 , which is a much smaller variation than the orders of magnitude differences achievable in electrical conductivity of similar systems. Both the low absolute magnitude of thermal conductivity and the relative insensitivity compared to electrical conductivity may be largely attributed to the poor interfacial thermal conductance of the many interfaces of the nanocrystal solid, but a direct experimental study of these interfaces is challenging. Here, in this work, we overcome this challenge via direct spectroscopic observation of heat flow within the components of a nanocrystal solid. These thermal transfers are studied by mixing two distinct types of particles: one that serves as a selectively excited antenna to inject heat and the other as the thermal acceptor to report the time-dependent change in temperature. Using transient spectroscopy, the equilibration between the heat donor and heat acceptor is observed to require similar to ~300 ps at room temperature, speeds up at reduced temperature, and has only weak sensitivity to the relative stoichiometry of the components or the intervening ligands. These results contrast strongly with the 10-20 ps time-scale of through-bond heat transfer at the ligand-particle surface and highlight the substantially lower interfacial thermal conductance of particle-to-particle transport without covalent bonding. It is also found, serendipitously, that the mixed composite films show an unexpected, substantially enhanced nonlinear absorption at the resonant wavelength of the plasmonic particles, which is tentatively attributed to a local field enhancement.
A theoretical analysis of the principles of interfacial thermal conductance between two heated metal surfaces in mutual contact is presented. It was determined that when heat flows from one metal to the other, the flow lines converge toward the actual contact points since the thermal conductivities of the metals are much greater than any intervening medium such as air or fluids. The converging lines cause contact resistance, which is usually high compared to the resistances offered to heat flow by metals. The resultant variation in resistance can produce unreliable equipment for aerospace and industrial applications unless the conductance can be determined accurately.
Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.
Here, the interfacial thermal conductance (h lig–NC ) between a cadmium selenide (CdSe) nanocrystal (NC) and three related organic ligands—olealdehyde, oleyl alcohol, and oleic acid—was investigated computationally. These ligands have the same carbon backbone but differ in the number and type of oxygen-coordinated headgroups (carbonyl and/or hydroxyl), leading to distinct bonding geometries involving monodentate and bidentate bonds. For a fully encapsulated NC, h lig–NC increases in the order of olealdehyde, oleic acid, and oleyl alcohol ligands. To isolate the contributions of h lig–NC from each headgroup type, the distinct bonding geometries were analyzed. Aldehyde and alcohol ligands, each featuring a single oxygen headgroup (carbonyl or hydroxyl), exhibit similar O–Cd separations and nearly identical h lig–NC per ligand at full surface coverage. However, the hydroxyl group in the alcohol ligand enables a higher ligand grafting density on the NC surface, resulting in a greater overall h lig–NC than the aldehyde-grafted NCs. In contrast, the oleic acid ligand forms multidentate bonds with the NC, leading to a shorter average O–Cd separation and a higher h lig–NC per ligand compared to monodentate bonds. Nevertheless, steric hindrance from the acid ligand's larger headgroup reduces its grafting density relative to alcohol ligands, ultimately resulting in a lower overall h lig–NC .
Statistical theory of interfacial thermal conductivity and crystal growth under weightlessness
We report optimization of thermal transport across the interface of two different materials is critical to micro-/nanoscale electronic, photonic, and phononic devices. Although several examples of compositional intermixing at the interfaces having a positive effect on interfacial thermal conductance (ITC) have been reported, an optimum arrangement has not yet been determined because of the large number of potential atomic configurations and the significant computational cost of evaluation. On the other hand, computation-driven materials design efforts are rising in popularity and importance. Yet, the scalability and transferability of machine learning models remain as challenges in creating a complete pipeline for the simulation and analysis of large molecular systems. In this work we present a scalable Bayesian optimization framework, which leverages dynamic spawning of jobs through the Message Passing Interface (MPI) to run multiple parallel molecular dynamics simulations within a parent MPI job to optimize heat transfer at the silicon and aluminum (Si/Al) interface. We found a maximum of 50% increase in the ITC when introducing a two-layer intermixed region that consists of a higher percentage of Si. Because of the random nature of the intermixing, the magnitude of increase in the ITC varies. We observed that both homogeneity/heterogeneity of the intermixing and the intrinsic stochastic nature of molecular dynamics simulations account for the variance in ITC.
Colloidal quantum dots (QDs) and their thin films are increasingly used in electronic and photonic devices replacing traditional bulk semiconductors. However, thermal properties of QDs remain underexplored relative to device development efforts. This study shows the use of time-resolved X-ray diffraction as a contact-free method to probe the thermal response of QDs in environments representative of the active layer in QD optoelectronic devices, providing in situ insights for future thermal management strategies. Through the extraction of Debye–Waller factors on a subnanosecond time scale, we directly capture the heating and cooling of core/shell CdSe/CdS QDs following pulsed optical excitation. In a QD thin film that actively provides optical gain, the thermal conductivity is found to be as low as 0.13 W m –1 K –1 , because of the poor heat flow between close-packed QD solids. For QDs dispersed in liquids, interfacial thermal conductance dominates thermal relaxation, with a conductance of 15 MW m –2 K –1 .
Time-resolved x-ray diffraction has been used to measure the low-temperature thermal transport properties of a Pt/Gd3Fe5O12//Gd3Ga5O12 metal/oxide heterostructure relevant to applications in spin caloritronics. A pulsed femtosecond optical signal produces a rapid temperature rise in the Pt layer, followed by heat transport into the Gd3Fe5O12 (GdIG) thin film and the Gd3Ga5O12 (GGG) substrate. The time dependence of x-ray diffraction from the GdIG layer was tracked using an accelerator-based femtosecond x-ray source. The ultrafast diffraction measurements probed the intensity of the GdIG (1 −1 2) x-ray reflection in a grazing-incidence x-ray diffraction geometry. The comparison of the variation of the diffracted x-ray intensity with a model including heat transport and the temperature dependence of the GdIG lattice parameter allows the thermal conductance of the Pt/GdIG and GdIG//GGG interfaces to be determined. Complementary synchrotron x-ray diffraction studies of the low-temperature thermal expansion properties of the GdIG layer provide a precise calibration of the temperature dependence of the GdIG lattice parameter. The interfacial thermal conductance of the Pt/GdIG and GdIG//GGG interfaces determined from the time-resolved diffraction study is of the same order of magnitude as previous reports for metal/oxide and epitaxial dielectric interfaces. The thermal parameters of the Pt/GdIG//GGG heterostructure will aid in the design and implementation of thermal transport devices and nanostructures.
Experimental thermal diffusivity data transverse to the fiber direction for composites composed of a reaction bonded silicon nitride matrix reinforced with uniaxially aligned carbon-coated silicon carbide fibers indicate the existence of a significant thermal barrier at the matrix-fiber interface. Calculations of the interfacial thermal conductances indicate that at 300 C and 1-atm N2, more than 90 percent of the heat conduction across the interface occurs by gaseous conduction. Good agreement is obtained between thermal conductance values for the oxidized composite at 1 atm calculated from the thermal conductivity of the N2 gas and those inferred from the data for the effective composite thermal conductivity.
Thermal conductivity of two‐dimensional (2D) materials is important to characterize as this will govern thermal transport physics in the various thermal and energy applications benefitting from the unique properties of these materials. Here, analytical and numerical techniques are presented to determine the thermal conductivity and interfacial conductance of a 2D material coating using confocal micro‐Raman spectroscopy as a noncontact diagnostic. While several methods have been proposed for calculating the values of thermal conductivity, interfacial conductance, and spatial temperature profile, these models often do not consider effects of radiation, convection, and substrate resistance on the temperature and heat flux profiles. We present a model to calculate the thermal conductivity and interfacial conductance which accounts for convection, radiation, and substrate effects to characterize a variety of 2D material coatings, which is demonstrated using large area graphene transferred onto copper and nickel substrates. Convection and radiation effects are found to have a negligible effect on the temperature profile of supported coatings, while the substrate effects have a considerable effect and therefore impact the quantification of thermal conductivity and interfacial conductance in addition to limiting the spatial resolution of the technique. These findings will allow for more accurate and reproducible extraction of the thermal conductivity of graphene and other 2D material coatings when employing noncontact optothermal Raman spectroscopy methods.
The effect of functionalization of carbon nanotubes on the thermal conductivity of nanocomposites has been studied using a multi-scale modeling approach. These results predict that grafting linear hydrocarbon chains to the surface of a single wall carbon nanotube with covalent chemical bonds should result in a significant increase in the thermal conductivity of these nanocomposites. This is due to the decrease in the interfacial thermal (Kapitza) resistance between the single wall carbon nanotube and the surrounding polymer matrix upon chemical functionalization. The nanocomposites studied here consist of single wall carbon nanotubes in a bulk poly(ethylene vinyl acetate) matrix. The nanotubes are functionalized by end-grafting linear hydrocarbon chains of varying length to the surface of the nanotube. The effect which this functionalization has on the interfacial thermal resistance is studied by molecular dynamics simulation. Interfacial thermal resistance values are calculated for a range of chemical grafting densities and with several chain lengths. These results are subsequently used in an analytical model to predict the resulting effect on the bulk thermal conductivity of the nanocomposite.
The transverse thermal conductivity of reaction-bonded Si3N4 is significantly affected by an interfacial barrier at the interface formed with SiC reinforcing fibers. A comparative study of composites with and without reinforcing-fiber carbon coatings found the coating to reduce effective thermal conductivity by a factor of about 2; this, however, is partially due to a thermal expansion-mismatch gap between fiber and matrix. HIPing of composites with coated fibers led to an enhancement of thermal conductivity via improved interfacial thermal contact and greater grain size and crystallinity of the fibers.
Gen3, Embedded Cooling, promises to revolutionize thermal management of advanced microelectronic systems by eliminating the sequential conductive and interfacial thermal resistances which dominate the present 'remote cooling' paradigm. Single-phase interchip microfluidic flow with high thermal conductivity chips and substrates has been used successfully to cool single transistors dissipating more than 40kW/sq cm, but efficient heat removal from transistor arrays, larger chips, and chip stacks operating at these prodigious heat fluxes would require the use of high vapor fraction (quality), two-phase cooling in intra- and inter-chip microgap channels. The motivation, as well as the challenges and opportunities associated with evaporative embedded cooling in realistic form factors, is the focus of this paper. The paper will begin with a brief review of the history of thermal packaging, reflecting the 70-year 'inward migration' of cooling technology from the computer-room, to the rack, and then to the single chip and multichip module with 'remote' or attached air- and liquid-cooled coldplates. Discussion of the limitations of this approach and recent results from single-phase embedded cooling will follow. This will set the stage for discussion of the development challenges associated with application of this Gen3 thermal management paradigm to commercial semiconductor hardware, including dealing with the effects of channel length, orientation, and manifold-driven centrifugal acceleration on the governing behavior.