Thermal Transport across Solid Interfaces with Nanoscale Imperfections: Effects of Roughness, Disorder, Dislocations, and Bonding on Thermal Boundary Conductance

e efficiency in modern technologies and green energy solutions has boiled down to a thermal engineering problem on the nanoscale. Due to themagnitudes of the thermalmean free paths approaching or overpassing typical length scales in nanomaterials (i.e., materials with length scales less than one micrometer), the thermal transport across interfaces can dictate the overall thermal resistance in nanosystems. However, the fundamental mechanisms driving these electron and phonon interactions at nanoscale interfaces are difficult to predict and control since the thermal boundary conductance across interfaces is intimately related to the characteristics of the interface (structure, bonding, geometry, etc.) in addition to the fundamental atomistic properties of the materials comprising the interface itself. In this paper, I review the recent experimental progress in understanding the interplay between interfacial properties on the atomic scale and thermal transport across solid interfaces. I focus this discussion speci�cally on the role of interfacial nanoscale “imperfections,” such as surface roughness, compositional disorder, atomic dislocations, or interfacial bonding. Each type of interfacial imperfection leads to different scattering mechanisms that can be used to control the thermal boundary conductance. is offers a unique avenue for controlling scattering and thermal transport in nanotechnology.


Introduction
As the world's population and power demands increase, our technological solutions continue to rely on power-hungry applications. e shear population increase in the world along with skyrocketing electricity, power, and transportation needs in emerging markets has led to necessary minimum levels of energy for sustainable growth and livelihood of human kind. e trends in energy use are marked by our society's continued advancement of technology, communication, transportation, and quality of life. For example, in 2010, the United States used an estimated 28.7 trillion kW-h of energy [1]. ese massive energy demands both in the USA and abroad are a necessary requirement for continued capabilities and qualities of life that have been made possible with the advancement of technologies. For example, increased computing needs (internet, cloud computing, and wireless communications), although more energy efficient in the recent past, have led to an increase in data center energy usage by roughly a factor of six over the past decade [2]; although technologies have become more energy efficient, due to the shear increase in the numbers of users and information, energy consumption and needs are skyrocketing, leading to these astounding levels of energy usage. e primary fuel for this electricity generation comes from coal, a nonrenewable resource that leads to 2,250 lbs/ MW-h of carbon dioxide, 13 lbs/MW-h of sulfer dioxide, and 6 lbs/MW-h of nitrogen oxides [2]. is leads to 12.6 billion lbs of CO 2 emissions, along with 73 and 34 million lbs of sulfer dioxide and nitrogen oxides, respectively. In addition, more than half (57%) of the used energy in the USA is rejected as heat where only 43% of the consumed energy is actually used (work). ese sti�ing statistics create two major foci for energy research: (1) improving the efficiency of the energy usage to reduce the net wasted heat and (2) reducing the harmful emissions of the technologies generating the energy. Improving energy efficiency of devices has been an active area of material research for the last few decades, stemming from the drive of the computing industry to maintain Moore's Law [3] and recently focused on mitigating the thermal transport in the nanoscale devices [4,5]. Recent novel approaches in reducing the net wasted heat have focused on waste heat recycling via thermoelectric power generation [6][7][8]. ermoelectric solutions also show promise in reducing harmful emissions in the grid, as the entirely solid state makeup of thermoelectrics requires no moving parts, cycles, or input fuel other than heat. ese thermoelectric solutions combined with "renewable source" technologies (technologies that do not rely on coal as a fuel-such as photovoltaics [9,10]) could substantially reduce the carbon footprint of the energy grid.
e issue with these technologies lies in their operation efficiencies, which at this stage of research have boiled down to a thermal engineering problem on the nanoscale. For example, power dissipation limits the performances of electronic systems from individual microprocessors to data centers. In silicon-based microprocessors (by far the most commonly implemented technology), operation beyond a few GHz is not possible due to on-chip power densities exceeding 100 W cm −2 , higher than a typical hotplate and far greater than typical cooling capabilities [4]. Even greater power densities have been observed in RF and ampli�er devices using III-V or group IV compound semiconductors (e.g., GaAs or SiC) [11]. ese huge power densities lead to hot-spots and device failure, warranting the need for increased thermal conduction to mitigate the heat load. As another example, current thermoelectric solutions for waste heat recovery are far outperformed by traditional solutions such as heat exchangers with heat pumps or heat engine cycles such as the Organic Rankine or Kalina cycles [6], making major implementations of the more-carbon-footprintfriendly thermoelectric devices impractical. Improvement in the thermoelectric efficiency can be achieved by controlled reduction in thermal transport of the thermoelectric materials, since the thermoelectric efficiency (quanti�ed by the Figure of Merit, ZT) is inversely related to the thermal conductivity. Recent massive improvements in ZT in thermoelectric materials have been discovered by nanostructuring [8,[12][13][14][15], leaving the research problem to improve energy efficiency centered around understanding thermal transport on the nanoscale. As both increases and decreases in thermal conductivity are required for different applications, the ability to control and tune the thermal properties of materials is therefore the ultimate goal to improve device efficiency and improve energy consumption trends.
is desired thermal control in nanosystems, however, is not a trivial task. Due to the magnitudes of the thermal mean free paths approaching or overpassing typical length scales in nanomaterials (i.e., materials with length scales less than 1 micron), the thermal transport across interfaces can dictate the overall resistance in nanosystems. is leads to a regime in which unique yet relatively unknown mechanisms such as interfacial scattering, ballistic transport, and waveeffects offer potential new degrees of freedom in the thermal engineering of material systems [5,16]. However, the fundamental mechanisms driving these electron and phonon Perfect interface Substrate Film Imperfect interface Film Substrate F 1: Schematic of a perfect and an atomically imperfect interface. A perfect interface is an atomically abrupt, coherent interface between two solids in intimate contact. e atomically imperfect interface depicted here occurs when the interface is not well de�ned due to atomic disorder (both structural and compositional), resulting in roughness and impurities around the interface.
interactions at nanoscale interfaces are difficult to predict and control since the thermal boundary conductance across interfaces (oen called the Kapitza conductance), ℎ K [17,18], is intimately related to the characteristics of the interface (structure, bonding, geometry, etc.) in addition to the fundamental atomistic properties of the materials comprising the interface itself. is interplay between interfacial properties and thermal boundary conductance can lead to additional carrier scattering mechanisms and phenomena that will not only change the thermal processes, but can also be used to systematically manipulate and control thermal transport in nanosystems.
As an example, consider a "perfect" interface and an atomically disordered, "imperfect" interface, depicted in Figure 1. e imperfect interface introduces roughness along with atomic imperfections that cause additional thermal scattering which can change ℎ K , as I have discussed in previous works [19][20][21][22][23][24] and will detail throughout this paper. e properties of this interface dictate the additional scattering mechanisms. erefore, an understanding of the effects of interfacial imperfections and other "non-idealities" (such as roughness, disorder, and bonding) on thermal transport across complex nanoscale interfaces is warranted to fully develop the potential to control electron and phonon thermal processes in nanosystems.
To this end, in this paper I will review recent advancements in the underlying physics of thermal boundary conductance at solid interfaces, focusing speci�cally on the role of interfacial nanoscale "imperfections, " such as surface roughness, compositional disorder, atomic dislocations, or interfacial bonding. In the next section, I present an overview of thermal boundary conductance and discuss pertinent measurements of thermal boundary conductance between two solids to put this thermophysical quantity into perspective. Following this background, I discuss time domain thermore-�ectance in Section 3. TDTR has become a common, effective experimental technique to accurately measure the thermal boundary conductance between two solids. In Section 4, I will discuss the mathematical semiclassical formalisms for predicting thermal boundary conductance, which are simple, analytical expressions that can be used to understand the electron or phonon transmission across interfaces. Following this, the remainder of the sections will focus on examining the effects of various interfacial imperfections on either the electron-or phonon-dominated thermal boundary conductance across interfaces. ese sections will focus on the effects of interfacial interspecies mixing on electron and phonon thermal boundary conductance (Sections 5 and 6, resp.) and roughness, dislocations, and bonding effects on phonon thermal boundary conductance (Sections 7, 8, and 9, resp.).

Background: Thermal
Boundary Conductance e thermal transport across solid interfaces is quanti�ed by the thermal boundary conductance, ℎ K , or Kapitza conductance [18], which quanti�es the e�cacy of heat �ow across an interface. Mathematically, the thermal boundary conductance is the proportionality constant that relates the heat �ux across an interface to the temperature drop associated with that interfacial region via where int is the thermal �ux across the solid interfacial region and Δ is the temperature drop across the interfacial region. It is also oen intuitive to conceptualize the inverse of this quantity, the thermal boundary resistance, K = 1/ℎ K , which can be paralleled to an electrical resistance network. Note that thermal boundary conductance is discussed here in terms of two materials in intimate contact, thereby making this inherently different than a "contact conductance" which is related to the voids or lack of contact between two solids. Measurements of the thermal boundary conductance between two solids typical occur in a thin �lm geometry or at ultralow temperatures. In both of these cases, the resistance due to the interface is a major resistance in the measured system of interest, thereby allowing for accurate measurement sensitivity. One of the �rst comprehensive reviews of thermal boundary conductance was detailed by Swartz and Pohl in 1989 [17], and their review mainly focused on low temperature measurements and phenomena. Recently, with the advances in sensitive measurement techniques for interrogating thermal transport on the nanoscale [25], accurate measurements of thermal boundary conductance in nanosystems and at elevated temperatures have been achieved [4,5].
As the goal of this paper is to discuss the effects of interfacial "non-idealities" on thermal boundary conductance, only a subset of the measurements of thermal transport across various interfaces will be discussed here. e reader is referred to [4,5,26] for overviews of thermal boundary conductance at "smooth" or seemingly "perfect" interfaces. However, it is informative to discuss the limits of thermal 75 100 200 300 400 500  [27], Al/Cu [28], TiN/MgO [29], Al/Si [19,21], Al/single layer graphene (SLG) [24], GaSb/GaAs [23], and Bi/diamond [30]). Typical values span over three orders of magnitude. For comparison, the equivalent conductances of various thicknesses of SiO 2 are shown as the solid lines. boundary conductance at "perfect" interfaces and the difference between electron-dominated and phonon-dominated thermal boundary conductances. Figure 2 shows the thermal boundary conductance as a function of temperature for various interfaces in which ℎ K is dominated by either electron or phonon scattering [21,23,24,[27][28][29][30]. e TiN/MgO and Bi/diamond data represent the highest and lowest phonon-dominated thermal boundary conductance ever measured, respectively [29,30]. Note that the electron-dominated thermal boundary conductances (Al/Cu and Pd/Ir) are roughly one order of magnitude higher than the phonon-dominated thermal boundary conductances and have a steeper increase with temperature [27,28]. Along with these data, the equivalent conductances of various thicknesses of SiO 2 are indicated by the solid lines; this equivalent conductance is calculated by ℎ SiO 2 = SiO 2 / , where SiO 2 is the thermal conductivity of SiO 2 and is the thickness, as indicated in [31]. erefore, at room temperature, the thermal boundary conductance between a Pd/Ir interface offers the equivalent resistance as roughly 1 Å of SiO 2 . e highest phonon-dominated thermal boundary conductance offers the same resistance as ∼2 nm of SiO 2 , nearly 20 times greater resistance than the electron-electrondominated interface. At the lower end of the spectrum, 4 ISRN Mechanical Engineering the phonon transport across a Bi/diamond interface is equivalent to the resistance of ∼200 nm of SiO 2 . Clearly, in nanomaterials and composites, the thermal boundary conductance across solid interfaces can offer a major impediment to heat �ow that can dominate the thermal transport in the overall system. e additional data in Figure 2 (Al/Si, Al/single layer graphene (SLG), and GaSb/GaAs interfaces) represent our recent measurements to examine the effects of some interfacial imperfections on thermal transport across interfaces [21,23,24]. e speci�c interfacial imperfection that was investigated in each material case is indicated by labels in Figure 2. ese, along with other measurements, will be discussed in detail throughout this paper. Brie�y, these three example interfaces show the effects of either roughness (Al/Si), bonding (Al/SLG), or dislocations (GaSb/GaAs) on the thermal boundary conductance. In the Al/Si case, a rougher substrate leads to a decrease in thermal boundary conductance [21]. At the Al/SLG interfaces, by increasing the bonding between otherwise nonreacting Al �lms and SLG substrates, Al can bond covalently to the SLG, thereby increasing the thermal transport [24]. At the GaSb/GaAs interfaces, the large lattice mismatch between the GaSb and GaAs leads to strain dislocations. e number of dislocations can be controlled via the growth technique, and we observed that making a GaSb/GaAs interface with a higher number of dislocations per unit area results in a lower ℎ K . ese effects-roughness, dislocations, and bonding-will be discussed in more detail in Sections 7, 8, and 9, respectively. However, it is apparent in this preliminary overview that the thermal boundary conductance between two solids is not only a function of the materials comprising the interface, but has a strong dependency on the atomic properties of the interface. erefore, with careful nanoscale engineering of the interfacial properties, the ability to tune the thermal boundary conductance between two solids can be achieved. In the remainder of this paper, I will discuss speci�c experimental results and accompanying theory on the thermal transport across imperfect interfaces between two solids, including more detail on the speci�c imperfect interfaces discussed above. However, before delving into the transport physics driving thermal boundary conductance across imperfect interfaces, measurements of ℎ K with time domain thermore�ectance will be detailed in the following section.

�� Time �omain Thermore�ectance�
An Effective Tool for Measuring Thermal Boundary Conductance between Two Solids Over the past few decades, contact-based techniques or experimental platforms based on resistive thermometry apparati have made substantial progress in measuring the thermal transport in nanosystems due to advances in nanofabrication and increasingly sensitive detection schemes [25]. Common resistive thermometry techniques that have been used to determine temperature rises and heat transport properties in nanosystems are steady-state resistance-based techniques [32][33][34][35][36][37] and the 3 technique [31,[38][39][40][41][42][43]. ese methods are based on the fact that the electrical resistivity of a sample or test strip is related to the temperature chang; therefore, by monitoring the change in electrical resistivity of a reference, the change in temperature can be determined. With precise knowledge of geometry of the nanosystem of interest, the thermal conductivity can be determined from the measured temperatures. e advantages of these techniques are that they are direct measurement of the temperature rise and can be directly related to the thermal conductivity of the sample of interest through only a carefully calibrated electrical resistivity reference and precise knowledge of the sample geometry. e disadvantages, however, in addition to the potential uncertainty in sample geometries, are that these techniques require intimate contact between the electrical leads and the sample platform of interest. Inherent in these contacts are both electrical and thermal contact resistances that must be accounted for [44]. In addition to these measurement difficulties, direct measurement of the thermal boundary conductance between two materials of interest is very difficult without the use of multiple samples. For example, if one were to use a steady-state or frequency domain (such as 3 ) electrical resistivity technique to measure the thermal properties of a bilayer sample (say, thin �lm on thin �lm on substrate), the measured data would represent the total thermal resistance of the sample. Although the individual �lm properties could be inferred from multilayer thermal modeling [45], a "heat capacity-less" thin interface could not be differentiated from the individual layers. is is not the case with time domain techniques, however. Time domain techniques for measuring the thermal transport of materials have proven invaluable for monitoring the thermal conductivity and thermal boundary conductance of multilayer structures due to the fact that the thermal effusivity or thermal diffusivity of a material responds differently than the thermal transport across an interfacial region in the time domain. Several time domain techniques have proven quite useful in measurement of these interfacial thermal properties, including laser �ash diffusivity [46], transient electrothermal techniques [47], and thermore�ectance techniques [48,49]. With recent progress in pulsed laser systems, thermore�ectance techniques employing short pulses have become a robust metrology to accurately determine the thermal transport in nanosystems, as I outline in the following. ese thermore�ectance techniques, which rely on the fundamentals of thermore�ectance spectroscopy of materials, can offer ultrafast characterization of the temperature change of materials in the time domain, which alleviates many of the radiative loss, optical effusivity, and contact issues that plague other techniques.
In thermore�ectance spectroscopy, an oscillating perturbation is applied to the sample of interest and the AC component of the re�ected light is monitored [50][51][52][53][54][55][56][57]. For example, in a standard pump-probe geometry, the pump is a modulated laser beam or train of laser pulses that is used to optically and thermally excite some sample at frequency, , while the probe is used to interrogate the small changes in re�ectivity of the sample caused by the pump excitation at frequency . To back out heat transport characteristics using thermore�ectance, the majority of optical perturbation that causes a change in re�ectivity must be directly due to a thermal perturbation. For this reason, metallic systems are ideal to transduce the change in re�ectivity into the change in temperature due to the large electronic number density around the Fermi level in metals. is resulting thermore�ectivity change in a metal is related to an electron population that is thermalizing or thermalized to a new distribution due to the pump pulse. Hence, the thermore�ectance of the metal, de�ned as the change in re�ectance due to a change in temperature, is given by where is the re�ectivity, 0 is the baseline temperature of the sample, and excited is the re�ectivity of the metal in the state during and aer pump pulse absorption. e thermore�ectance of a metal is dictated by its thermore-�ectance coe�cient, , which is typically small, on the order of 10 −6 -10 −4 K −1 depending on the temperature and metallic band structure at a given wavelength [58][59][60]. erefore lock-in techniques are typically used to measure the thermore�ectance of metals in pump-probe geometries. Using ultrashort laser pulses in pump-probe geometries, the time evolution of the thermore�ectance can be measured [49,61].
e time evolution of the thermore�ectance signal is the basis behind time domain thermore�ectance (TDTR). TDTR has been used in several different geometries with various different pulsed lasers, ranging from nanoseconds to femtoseconds [5,25,48,49,[62][63][64][65][66][67]. ese experiments can be separated into two classes: (1) experiments in which the laser pulses are separated far enough in time in which the interrogated sample is relaxed back down to its equilibrium state before the next pulse arrives and (2) experiments in which the pulses are so closely spaced in time that the pulse arrives at the sample surface before full thermal relaxation of the previous pulse. Recently [68], the experiments in the former class have been characterized as "transient thermore�ectance� where experiments utilizing laser systems producing the later class of excitation have been commonly referred to as "time domain thermore�ectance� (TDTR). In TDTR experiments, the buildup of laser pulse energy leads to a steady-state temperature rise in the sample [69] that can provide additional information for thermal analysis. is paper will limit the review of measurement techniques to TDTR, as the majority of the data presented in the later sections were measured with TDTR techniques.
Due to the recent decrease in cost of femtosecond Ti : Sapphire oscillators, most TDTR systems are centered around these high frequency (∼80 MHz) laser systems. A typical cost of a Ti : Sapphire oscillator from a commercial vendor is ∼$100,000-$200,000, depending on the output power and pulse width. e TDTR experiment that I have in my lab at University of Virginia is based off of a Spectra Physics Tsunami with ∼3 W output of 90 fs laser pulses centered at 800 nm. is means that 90 fs laser pulses are emanating from the oscillator with about 37.5 J of energy every 12.5 ns. e pulse energy is split with a variable beam splitter into a pump path and a probe path. e pump path is directed through an electro-optic modulator (EOM) that can modulate the pulse train upwards of = MHz (I use the electro-optic modulators sold commercially through ConOptics). In my particular setup, I then focus the modulated pump pulses through a bismuth triborate (BiBO) crystal to frequency double the near infrared pulses (800 nm = 1.55 eV) to visible wavelengths (400 nm = 3.11 eV). I then direct the pump pulses through an objective to focus the 400 nm pump pulses onto the sample surface. e modulated pump pulses produce a temperature rise on the sample surface at the frequency of the EOM. e pulse energy that is not directed down the pump path aer the variable beam splitter is directed down a mechanical delay stage and then focused collinearly with the pump pulses through the objective; these probe pulses are focused on the sample surface in the middle of the pump pulses. I typically make the pump spot size larger than the probe spot size on the sample surface. is reduces error in our measurements due to probe beam dri and movement. is also creates a onedimensional heating event in the sample, cross-plane to the sample surface (note, if the pump beam waist is decreased to the same spot size as the pump, in-plane thermal transport properties can also be measured; however, this is not of interest for this paper since I will only focus on measurements of thermal boundary conductance across a thin interfacial region where heat transfer is nearly entirely cross-plane) [65,66]. I typically use in�nity corrected objectives with anywhere from 5×-20× magni�cation, which leads to ranges of pump and probe 1 radii anywhere from 25-6 m and 6-4 m, respectively, depending on the objective. Note that this is very dependent on the particular setup and the spot size of the pump and probe going into the back face of the objective. e re�ected probe beam is then re�ected along the same incident path, passed through polarizing and color �ltering optics, and directed to a photodiode. e voltage generated from the photodiode contains both a DC signal and a small AC component which represents the change in the probe beam re�ectivity due to the pump heating event at frequency . is AC component is �ltered from the DC component with lock-in ampli�cation, and the in-phase and out-of-phase components of the AC signal are monitored as a function of pump-probe delay time, which we can monitor over ∼7 ns of delay time (although typically I �nd that 4-5 ns is su�cient to exact the thermal boundary conductance of the majority of the samples of interest). A schematic of the TDTR system in my lab at the University of Virginia is shown in Figure 3.
TDTR data taken on a 117 nm Al �lm evaporated onto a single crystalline Si substrate are shown in Figure 4. Characteristic in the TDTR data is a rapid rise in signal which is related to the increase in electronic temperature due to the pulse absorption followed by an exponential decay with different time constants related to the various cooling processes in the Al �lm. In the �rst ∼10 ps, the signal decay is related to the hot electron equilibration with the lattice, indicative of short pulsed laser heating [  thermal diffusion through the Al �lm occurs, which typically occurs on a time scale of 10-100 ps [48,80]. From 100 ps to several nanoseconds, the TDTR signal is related to the thermal boundary conductance across the Al/Si interface and the thermal effusivity or thermal diffusivity into the Si, depending on the modulation frequency of the pump [81]. For purposes of measuring the thermal boundary conductance between two solids, the TDTR data between ∼100 ps to a few nanoseconds contains the pertinent information.
To quantify the thermal boundary conductance from the TDTR data, an accurate thermal model that relates the measured lock-in output to the thermal excitation and propagation in the structures of interest must be used. For most TDTR experiments, the use of Gaussian pump and probe shapes along with a metal transducer that is thicker than the optical penetration depths facilitates the solution to the cylindrical heat equation in a multilayer structure (a metal �lm that is ∼80-120 nm should suffice) [69,82,83]. e speci�cs of the thermal model to predict the temperature change on the surface of the metal transducer have been outlined in several works previously [65-67, 69, 84]. It is important to note that due to the laser pulse accumulation in the sample from a Ti : Sapphire oscillator, the temperature rise on the sample surface consists of two parts: (1) a real component that represents the temperature rise due to the excitation of a single pulse and subsequent thermal diffusion from the surface of the metal into the underlying material and across the material interfaces; this real response is directly related to the time domain response of the surface temperature following the optical pump pulse; and (2) an imaginary component that represents the temperature rise from the pulse accumulation; this imaginary response is predominantly controlled by the change in surface temperature subjected to the periodic heat source at frequency . Both the real and imaginary components of the heating event at the sample surface contain information about the thermal properties of the sample and can be used to analyze TDTR data for robust data analysis. e temperature rise predicted from the thermal model is related to the measured output of the lock-in ampli�er during TDTR experiments via where in is the in-phase component of the lock-in transfer function and related to the real component of the temperature rise and out is the out-of-phase component of the lockin transfer function and related to the imaginary component of the temperature rise.
, which is the lock-in transfer function, is a function of pump modulation frequency, , and surface temperature rise, , and is given by where is a constant related to the electronics, thermore-�ectance coefficient, laser powers, and frequency of the laser source, laser , and is the pump-probe delay time. By nature of (4), pulse accumulation is taken into account. Solving (3) for various pump-probe delay times, , yields the output of the lock-in ampli�er that is related to the measured data through the various thermal properties of the system. To analyze the data, the negative ratio of in to out − in / out yields a quantity that is directly related to the real temperature rise on the surface of the sample and inversely related to the imaginary temperature rise on the surface of the sample. Using this ratio is much more robust than using either in and out alone since it minimizes error due to probe beam dri and dephasing, decreases the noise due to the electronics, and eliminates the need to scale the model to the experimental data or need to know the thermore�ectance coefficient of the metal transducer [29,85]. In Figure 4, I show the best �t of the model to the TDTR data along with the values of Al/Si ℎ K and the Si that achieve this best �t. I oen use an Si substrate as a calibration of TDTR to ensure our experimental measurement can accurately reproduce the thermal conductivity of single crystalline silicon at room temperature [86,87] (note, substantially below room temperature, the thermal conductivity of single crystalline silicon cannot be accurately measured due to ballistic phonon effects, which will not be discussed in this paper; however, the reader is referred to [88,89] that discussed this mean free path effect in more detail).

Semiclassical Formalism for Predicting Thermal Boundary Conductance
Semiclassical modeling based on analytical expressions of the electron or phonon �ux from a particle basis serves as a robust, relatively simple approach to gain fundamental insight into the underlying process driving thermal boundary conductance across solid interfaces. Although there are more rigorous approaches to model the thermal boundary conductance between two solids, such as molecular dynamics simulations [21,[90][91][92][93][94][95] or nonequilibrium Green's Function methods [96][97][98][99][100], due to their simplicity, semiclassical approaches can be very powerful for quick predictions that require very minimal computational expenses. erefore, I will restrict the discussion of thermal boundary conductance modeling in this paper to semiclassical approaches and draw on speci�c variations of these approaches for modeling imperfect interfaces in the sections that follow. e two most common approaches to modeling thermal boundary conductance are rooted in the assumptions of either a perfectly specular interface or a perfectly diffusive interface, resulting in specular or diffuse scattering of the electrons or phonons. Typical approaches for modeling the specular re�ection of electron or phonon waves at interfaces use transmission calculations rooted in the Fresnel equations [101]. An example of this is the Acoustic Mismatch Model for predicting the thermal boundary conductance of phonons due to specular phonon wave interactions and transmission at interfaces [102,103]. ese approaches are rooted in the assumptions that the electrons or phonons will see the interface as a perfectly specular wall. From simple analytical arguments, it is straight forward to show that this assumption breaks down at temperatures above ∼30 K for most materials [104], even in the event of perfect interfaces as was demonstrated experimentally [29]. However, in the event that imperfections or disorder are present at the interface, this specular assumption breaks down, even at temperatures as low as 1 K for certain interfaces [105]. In this case, the assumption of diffusive particle scattering must be employed to model the thermal transport across interfaces.
e most common and straight forward approach for modeling diffusive phonon scattering across solid interfaces and the resulting thermal boundary conductance is implementation of the Diffuse Mismatch Model (DMM) [17,105]. e DMM is rooted in the fact that electrons or phonons scatter diffusively at the interface between two materials, and therefore the probability of the energy of these carriers to transmit energy from material 1 to material 2 is the probability of energy re�ection from material 2 to material 1. Mathematically, this is expressed as where and are the probability of energy transmission and re�ection, respectively. is diffusive transmission is then multiplied by the phonon �ux in side 2 incident on the interface between material 1 and material 2 to give the transmitted phonon �ux from side 1 to side 2, which, in the semiclassical formalism, is approximated by [ where is an index corresponding to the energy level (or phonon mode/branch), is the energy of the carriers, is the statistical distribution function of the carriers, is the velocity of the carriers, and min and max are the minimum and maximum energy of the carrier in level . From (1), in the limit of an in�nitesimally small temperature rise, the semiclassical, "particle" formalism for the thermal boundary conductance across interfaces is given by which is applicable for both electrons and phonons (and photons, for that matter). More rigorous derivations of ℎ K,1,2 are discussed elsewhere [95]. In addition, several works have examined the assumptions going into calculations of (7), including the use of phonon dispersion [106], realistic density of states [107], and isotropy [108][109][110]. e calculation of 1,2 relies on the de�nition of diffusive scattering (5) and the fact that a diffusive scattering event at an interface will lead to a black, thermalizing boundary. In this case, the principle of detailed balance can be evoked [111] and the DMM transmission is given by Note that this transmission expression is valid for both electrons and phonons (and, again, photons) scattering diffusively at interfaces. I will not discuss the speci�cs of the transmission calculations in this paper, since the goal of this paper is to review recent advances in experimental observations of thermal boundary conductance at rough interfaces. However, several additional works have described the derivation of the DMM transmission coefficient in great detail, and the reader is directed to the following references to understand the fundamental assumptions of detailed balance and how it is applied to the DMM transmission [17,104], the role of inelastic scattering on this energy transmission [112][113][114], and the role of �uxes and transport around these interfaces [95,115,116]. I will discuss the modi�cations to (8) that have been made to describe the additional processes affecting transmission across imperfect interfaces in the sections that follow.

Thermal Transport between Two Metals: The Electron-Dominated Thermal Boundary Conductance and the Effects of Compositional Mixing
Typically, the thermal boundary conductance in a system is only a major concern when the interfacial transport is low compared to the transport in either of the materials comprising the interface. However, at metal/metal interfaces, due to the relatively high thermal conductivity of pure metals, the transport across these interfaces can be a limiting resistance, even though the thermal boundary conductance can reach upwards of 10,000 MW m −2 K −1 at room temperature. ere have only been limited experimental measurements of the thermal boundary conductance across metal/ metal interfaces compared to those at interfaces comprised of at least one nonmetal. One of the �rst studies by Clemens et al. [117] reported on the thermal boundary conductance across a series of metal/metal bilayer and multilayer �lms. In the bilayer �lms, they were able to directly measure the thermal boundary conductance across Ni/Ti and Ni/Zr interfaces. ese values, shown in Figure 5, are similar to the room temperature values of the TiN/MgO interface [29], which as discussed previously, is the highest phonon-dominated thermal boundary conductance ever measured. Gundrum et al. [28] measured the thermal boundary conductance across an interface between two free electron metals, Al and Cu. e thermal boundary conductance was ∼4,000 MW m −2 K −1 at room temperature, half-an-order of magnitude larger than the highest phonon-dominated thermal boundary conductance ever measured. is was theorized to be due to purely free-electron dominated transport across the Al/Cu interfaces. ey derived a variation of the DMM for free electrons, given by where is the linear coefficient to the electron heat capacity (oen referred to as the Sommerfeld coefficient) [118] and F is the Fermi velocity (recently, we have derived a more general form of this electron DMM for nonfree electron metals [119]). Calculations of (9) for an Al/Cu interface are shown as the line in Figure 5, which show acceptable agreement with the experimental data, alluding to the fact that this interface is dominated by free electron transport.
Recently, Wilson and Cahill [27] followed up to this Al/Cu measurement by measuring the thermal boundary conductance deduced from thermal conductivity measurements of Pd/Ir multilayers. ese multilayers had very low electrical resistivity, and the deduced thermal boundary conductance is, to date, the highest thermal boundary conductance ever measured. eir work also demonstrated that the Wiedemann-Franz Law [120] is valid at solid interfaces [121], further proving that the transport at metal/metal interfaces is an electron-dominated phenomenon. e origin of the low thermal boundary conductance of the Ni/Ti and Ni/Zr interfaces compared to the Pd/Ir and Al/Cu interfaces is currently unknown; however, strong electron-phonon F 5: ermal boundary conductance at metal-metal interfaces (Pd/Ir [27], Al/Cu [28], Ni/Ti [117], and Ni/Zr [117]) compared to the highest phonon-dominated thermal boundary conductance ever measured (TiN/MgO [29]). Inset: thermal boundary conductance at the Al/Cu interfaces as a function of ion dose.
coupling effects in the Ni, Ti, and Zr transition metals, which would lead to high electrical resistivities compared to Al, Cu, or the Pd/Ir multilayers, could be playing a role. However, with the recent result that the Wiedemann-Franz Law is valid at metal/metal interfaces [27], these underlying mechanisms could be deduced from electrical resistivity data. Due to this validity of the Wiedemann-Franz law at metalmetal interfaces, disorder, which would lead to an enhancement in electron-phonon and electron-defect scattering and thereby decrease the electrical conductivity (increase the electrical resistivity), should also lead to a decrease in electron thermal boundary conductance. is was veri�ed by Gundrum et al. [28] by ion irradiating the Al/Cu interfaces with 1 MeV Kr ions to induce an AlCu intermixed layer at the interface, where the thickness of this layer would increase with increasing ion dose. As seen in the inset of Figure 5, the thermal boundary conductance across the Al/Cu interface systematically decreases with an increase in ion dose. is is consistent with the increase in thickness of the AlCu intermixed layer, which leads to increased electron scattering around the interface and therefore a decrease in the thermal boundary conductance from the Al to the Cu.
In this section, I reviewed the limited experimental data on thermal boundary conductance at metal/metal interfaces and the effects of AlCu mixing at the interface leading to a decrease in thermal boundary conductance. In general, metal/metal interfaces do not pose a limiting thermal resistance in most technology due to their relatively high thermal boundary conductances compared to metal/nonmetal or nonmetal/nonmetal interfaces where phonon scattering dominates transport. erefore, for the remainder of this paper, I will focus on the effects of interfacial imperfections on phonon transport by considering different types of imperfections.

Effects of Compositional Intermixing on Phonon Thermal Boundary Conductance
In the previous section, I reviewed the effects of intermixing of Al and Cu on electron scattering and resulting thermal boundary conductance. At the Al/Cu interface, the thermal boundary conductance is dominated by electron transport. However, at interfaces in which one material is a nonmetal, the thermal boundary conductance is dominated by phonon transport [30].
We have studied the problem of compositional intermixing on phonon thermal boundary conductance by measuring ℎ K across Cr/Si interfaces with varying degrees of CrSi intermixed regions [22]. We accomplished this via in situ processing of the Si surface prior to Cr �lm sputtering followed by in situ heat treatments. is led to compositionally mixed regions of CrSi around the Cr/Si interface that were measured with Auger Electron Spectroscopy and Transmission Electron Microscopy. We measured the thermal boundary conductance from the Cr to Si with a variation of TDTR with a regeneratively ampli�ed laser system (repetition rate 250 kHz) [22,122]. e measured thermal boundary conductance across these Cr/Si interfaces is shown in Figure 6. e measured thermal boundary conductance decreases as the compositionally mixed layer increases. For comparison, I show the thermal boundary conductance from the Al/Cu interfaces reported by Gundrum et al. [28] and discussed in Section 5. e dependency of the Al/Cu electron dominated thermal boundary conductance as a function of layer thickness is similar to that at the Cr/Si phonon dominated interface-that is, as the thickness of the compositionally mixed interfacial layer is increased, the thermal boundary conductance decreases.
To understand the underlying phonon mechanisms driving thermal boundary conductance at compositionally mixed interfacial regions, we extended the DMM for phonon transport outlined in Section 4 to a three-layer system. We have discussed this approach in detail in our previous works [123,124]. In short, in this approach, we described the interfacial, compositionally mixed region as a thin layer of thickness with averaged properties of material 1 and material 2 comprising the interface. erefore, for the Cr/Si interface, the interfacial layer had averaged properties of Cr and Si. e properties of this compositionally mixed interfacial layer were calculated via the virtual crystal approximation [125]. e thermal boundary conductance was then calculated by F 6: ermal boundary conductance of a series of Cr/Si interfaces as a function of CrSi mixed layer thickness [22]. For comparison, I also show the thermal boundary conductance from the Al/Cu interfaces reported by Gundrum et al. [28] and discussed in Section 5. Compositional mixing leads to a decrease in both electron-and phonon-dominated thermal boundary conductance.
where "int" refers to the interfacial mixture, the properties which are described by a virtual crystal averaged property of material 1 and material 2, and ℎ K is calculated via (7). e transmission coefficient was then calculated via (8), where the phonon �uxes in each of the materials and the interfacial layer are calculated via (6). As this model is based on the phonon DMM, we referred to this model as the Virtual Crystal DMM (VCDMM). Calculations of our VCDMM as a function of of mixed layer thickness, , for a Cr/Si interface including a varying width of a CrSi compositionally mixed region are shown in Figure 6. Details of these particular calculations are discussed in our previous work [22]. e model shows very good agreement with our experimentally measured data. e calculations discussed above for the VCDMM for the Cr/Si interface were speci�c for phonon interactions and scattering in a compositionally disordered CrSi layer of varying thickness . Inherent to these calculations were the virtual crystal assumption of averaged layer properties based on a model derived for crystalline alloys [125]. e properties of an interfacial mixed layer can take on a range of structures and morphologies, and depending on the materials comprising the interface, the region may or may not form an alloy compound. We developed an extension to the VCDMM to account for complete disorder in an intermixed region (i.e., amorphous) [124]; however, in general, these approaches must be further validated. With the limited experimental data on the effects of compositional mixing on thermal boundary conductance, it is difficult to draw any concrete conclusions on either the fundamental electron or phonon mechanisms driving transport in an interfacial mixed region. However, it is apparent that the presence of this compositionally mixed interfacial region will decrease the interfacial thermal transport from one material to another.

Effects of Roughness on Phonon Thermal
Boundary Conductance e previous section discussed the effects of compositional mixing on phonon-mediated thermal boundary conductance. In this case, the interface is not compositionally abrupt but consists of some region of graded compositional properties that was referred to as a mixed layer. is begs the question, what type of thermal perturbations and phonon scattering mechanisms can exist at compositional abrupt, yet nonperfect interfaces? We studied this problem by considering thermal transport across interfaces between thin Al �lms and geometrically roughened Si substrates. With pre-Al-deposition Si processing, the silicon surface created a geometrically rough boundary between the Al and Si. However, by ensuring that the native oxide layer reforms on the Si aer processing, the Al will not migrate or spike into the Si as the native oxide layer in the Si acts as a diffusion barrier. is creates a compositionally abrupt, yet geometrically rough Al/Si interface. In a recent series of papers [19][20][21], we explored the extent of manipulation of phonon thermal boundary conductance through either chemical roughening of the Si [20,21] or creating roughness by synthesizing Si 1− Ge quantum dots on an Si surface. In both situations, the samples were exposed to ambient atmosphere to allow the formation of a silica or germania native oxide layer aer surface modi�cation to serve as a diffusion barrier to the Al transducer, as mentioned above; this was con�rmed with cross-sectional transmission electron microscopy [21]. We quanti�ed the roughness of the various surfaces with atomic force microscopy and used the root mean square (RMS) roughness of the silicon surface ( ) as the metric for comparison. We measured the thermal boundary conductance across the Al/Si interfaces from 100 to 300 K with TDTR. e thermal boundary conductance as a function of temperature for the Al/chemically roughened Si interfaces [21] is shown in Figure 7 (�lled symbols). In addition, I have included the Kapitza conductance at a nominally �at and oxide-free Al/Si interface as reported in [89] (open squares). As the data indicate, even a thin native oxide layer at the interface substantially reduces the effective Kapitza conductance (>50% reduction at room temperature). In addition, these two datasets demonstrate signi�cantly different temperature dependencies, suggesting that the oxide layer inhibits multiple-phonon scattering events which would otherwise contribute to Kapitza conductance [112,113]. Similarly, comparing the four datasets of the present study,  [89], and the �lled symbols are the data measured in our previous work [21]. It is evident that both the presence of a native oxide layer and interface roughness can have a signi�cant effect on Kapitza conductance. Not only does roughness decrease Kapitza conductance, but also it suppresses the temperature dependence. e agreement between the dash/dot lines and the data suggests that the DMM can be adjusted to take into account both the presence of an oxide layer and interface roughness.
increased interface roughness both reduces the magnitude of Kapitza conductance as well as suppresses its temperature dependence-that is, Kapitza conductance is less temperature dependent as interface roughness increases. In addition to the data, I have plotted several different predictive models in Figure 7. All models are calculated assuming that elastic phonon-phonon interactions dominate Kapitza conductance-that is, phonons in silicon at frequencies higher than the maximum phonon frequency of aluminum do not participate in transport. e diffuse mismatch model [17] is calculated using assumptions that we have previously outlined in [106], where the vibrational properties of �lm and substrate are approximated by �tting polynomials to the phonon dispersion curves of aluminum [126] and silicon [127] along the [100] crystallographic direction; assumed spherical Brillouin zones are then constructed via an isotropic revolution of these polynomial �ts in wavevector space. As seen in Figure 7, the prediction of the DMM falls between the data of [89] and that at the smoothest interface presently considered. is suggests that without an oxide layer, inelastic phonon-phonon scattering could play a role in thermal transport across aluminum-silicon interfaces [112,113]. On the other hand, we attribute the difference between the predicted and measured values at the smoothest interface considered (black squares) to the native oxide layer. e conductance of this oxide layer is described by its thermal conductivity divided by its thickness, When evaluating (12), we use the temperature-dependent bulk thermal conductivity of :SiO 2 , as it has been shown that the thermal conductivity of thin �lm :SiO 2 does not substantially differ from that of bulk [29]. A series resistor approach then yields is prediction is represented by the solid black line in Figure  7 and agrees well with our experimental data. In order to take interfacial roughness into account, we introduce a spectral attenuation coefficient (discussed in detail in our previous work [19]) and insert this coefficient into the integral expression of the DMM. is coefficient, , is unity when the phonon wavelength, , is greater than the RMS roughness, . On the other hand, = − when . at is, phonons with wavelengths greater than are unaffected by the roughness of the interface, whereas those with wavelengths less than are affected in a fashion similar to that of photons in an absorptive media, for example, the Beer-Lambert law. Qualitatively speaking, this approach suggests that as the "absorptivity" of the interface increases, so too does the temperature drop across it. With the spectral attenuation coefficient implemented, the DMM is once again plotted in Figure 7 for roughnesses of 6.5 nm and 11.4 nm. As is evident in the plot, this approach not only captures the reduction in Kapitza conductance due to interface roughening, but also captures the reduction in temperature dependence as well.
Finally, I plot room temperature Kapitza conductance as a function of RMS roughness in Figure 8 (squares), comparing the chemically roughened data from our paper [21], the aforementioned roughness DMM model calculated at 300 K, and two prior sets of experimental data from our previous works in which we explored other methods of roughening the silicon surface (different chemical etchants (le pointing triangles) [20] and quantum dot roughening (circles) [19]), as mentioned at the beginning of this section. Regardless of the silicon roughening procedure, the thermal boundary conductances all follows a similar relationship with RMS roughness.
To summarize, we have extensively characterized roughness effects on thermal boundary conductance across Al/Si interfaces. In general, we found that an increase in RMS roughness leads to a systematic decrease in thermal boundary conductance. We developed an extension to the DMM that is based around phonon wavelength considerations in DMM + oxide + roughness Chemical etch [21] Quantum dots [19] Chemical etch [20] Thermal boundary conductance, F 8: Room-temperature predicted (line) and measured thermal boundary conductance at Al/Si interfaces plotted as a function of interface roughness. Regardless of the silicon roughening procedure (indicated in the plot), the thermal boundary conductances all follows a similar relationship with RMS roughness.
which only certain phonon wavelengths see the interfacial roughness. is model captures both the magnitude of the experimental data and the temperature and roughness trends well. However, this RMS roughness phenomenon, to my knowledge, has only been intricately explored with wellcharacterized interfaces in our three recent works [19][20][21]. Different materials with different phonon spectra and oxide reactivities (e.g., gold, which does not wet a Si native oxide layer) should be explored to generalize the phonon physics at geometrically rough interfaces.

Effects of Lattice Dislocations on Phonon
Thermal Boundary Conductance us far, I have discussed the effects of compositional mixing and geometric roughness on thermal transport across solid interfaces. ese interfacial characteristics can commonly be produced by nanoscale processing or modifying the growth/synthesis conditions. However, combinations of materials can lead to additional interfacial imperfections based on the crystal structure and atomic properties. e most common imperfections that naturally form at the interface between two solids are dislocations produced from strain relaxation that stem from the lattice mismatch between two materials. Dislocations, which can readily scatter phonons and add to thermal resistance in solids [128], have been shown to cause drastic reductions in the thermal conductivity of thin �lms and nanostructures. For example, at large periodic thicknesses (∼10 nm), the thermal conductivity of Si/Ge superlattices has been shown to suddenly drop by a factor of two due to dislocation formation from strain relaxation of the Si or Ge layers on Ge or Si, respectively, leading to poor crystalline quality [129]. Similarly, threading dislocation densities of × 10 10 cm −1 have led to AlN �lms on Si� substrates to exhibit thermal conductivities less than the theoretical minimum limit [130,131]. As these dislocations stem from interfacial growth and lattice relaxation, understanding how dislocations affect thermal boundary conductance adds an important piece to the study of the phonon mechanisms affecting thermal transport across imperfect interfaces.
is issue was addressed in part by Costescu et al. [29], who measured the thermal conductivity across TiN/Al 2 O 3 interfaces in addition to TiN/MgO interfaces (one of the TiN/MgO samples that was measured is shown in Figure 2). e measured TiN/MgO (001), which is essentially free of mis�t dislocations, TiN/MgO (111), which contains a high density of stacking faults but few mis�t dislocations, and TiN/Al 2 O 3 , which should contain a large density of both stacking faults and mis�t dislocations, show no appreciable difference in the thermal boundary conductance. is indicates that in these well acoustically matched materials, the dislocation densities did not affect the phonon transport. However, without quantifying the dislocation density, it is difficult to understand the effects of these dislocations on thermal boundary conductance. In addition, without quan-ti�cation, this does not add any direct insight into previous studies observing drastically reduced thermal conductivities of thin �lms due to dislocation formation (such as the AlN thin �lm example above [130]).
We quantitatively studied the effects of dislocation density on thermal boundary conductance by measuring the thermal boundary conductance across well-characterized GaSb/GaAs interfaces, as is outlined in detail in our previous work [23]. e lattice mismatch between GaSb and GaAs is 7.78%, and the critical thickness for strain relaxation is formed within the �rst monolayer. Different epitaxial growth modes were utilized to change the dislocation densities at the GaSb/GaAs interfaces, as characterized with plan view and cross-sectional transmission electron microscopy [132,133]. ese growth techniques were able to produce GaSb/GaAs interfaces with two ranges of dislocations densities: 5 × 10 6 -5 × 10 8 dislocations cm −2 and 10 9 -10 11 dislocations cm −2 . We measured the thermal boundary conductances across the GaSb/GaAs interfaces from 100-500 K with TDTR.
e thermal boundary conductances as a function of temperature for the two GaSb/GaAs interfaces with different dislocation densities are shown in Figure 9 (�lled and open triangles for the relatively low and high dislocation dense samples, resp.). It is interesting to note that roughly a two order of magnitude change in dislocation density only leads to a factor of two change in thermal boundary conductance. To put this into perspective, at similar dislocation levels, for two orders of magnitude change in dislocation density, the thermal conductivity of GaN drops by an order of magnitude [130]. is indicates the seemingly small effect that changes in dislocation densities have on thermal boundary conductance across heavily dislocated interfaces. For comparison, I also plot the thermal boundary conductance of TiN/MgO [29] and Bi/diamond [30] interfaces (as discussed in Section 2, TiN/MgO and Bi/diamond represent the highest and lowest phonon-mediated thermal boundary conductances ever measured). e GaSb/GaAs dislocation F 9: ermal boundary conductances at the two different dislocation dense GaSb/GaAs interfaces from our previous work [23]. ese interfaces have dislocation densities ranging from 5 × 10 6 -5 × 10 8 dislocations cm −2 and 10 9 -10 11 dislocations cm −2 . For comparison, the plot includes previously measured data at TiN/MgO [29] and Bi/diamond [30] interfaces. DMM calculations are shown for the GaSb/GaAs interfaces, which drastically overpredict the measured data since the traditional formulation of the DMM does not account for dislocations interfaces. It is also interesting to note that roughly a two order of magnitude change in dislocation density only leads to a factor of two change in thermal boundary conductance. dense interfaces are approaching the measured values of Bi/diamond. However, Bi and diamond have drastically different phonon spectra and are extremely acoustically mismatched (over an order of magnitude different acoustic velocities and Debye Temperatures) [134]. Based on DMM theory, acoustically mismatched interfaces will have lower phonon transmission probability than acoustically matched interfaces [17]. GaSb and GaAs have similar phonon spectra and phonon velocities and are well acoustically matched. e TiN/MgO interfaces are well acoustically matched and were epitaxially grown (like the GaSb). e differences in thermal boundary conductances between the GaSb/GaAs and TiN/MgO well acoustically matched interfaces are over an order of magnitude, indicating the potential strong effect that dislocations can have on phonon transmission. is is further exempli�ed by comparing our measured thermal boundary conductances across the GaSb/GaAs interfaces to DMM calculations for GaSb/GaAs. For these calculations, we used the approach that we outlined previously in which we �t a polynomial to a measured or simulated phonon dispersion in one crystallographic direction and assume an isotropic medium to calculate the DMM [106]. For DMM calculations, we used the dispersion in the Γ → X direction of the Brillouin zone from [135] for GaSb and [136] for GaAs (we include the optical branches of GaSb and GaAs in our calculations) [137]. I plot these calculations of the DMM for GaSb/GaAs as the solid line in Figure 9. Clearly, the DMM greatly overpredicts the measured values. As the DMM assumes a perfect interface and a single-phonon scattering event, this approach is not valid for the dislocation dense GaSb/GaAs interfaces.
e discussion in this section demonstrates the detrimental effects that dislocations can have on thermal boundary conductance. However, due to the limited data, the fundamental mechanisms driving phonon dislocation interactions around interfaces are unknown. One outstanding question is the functional form of the dependency of thermal boundary conductance on dislocation density. In the thermal conductivity of a homogeneous material, the thermal conductivity is inversely proportional to the square root of the dislocation density [130]. is implies that changes in dislocation densities at a dislocation light interface will be greater than at a dislocation dense interface. Assuming the GaSb/GaAs interface is in the dislocation dense regime, this could explain the relatively small change in thermal boundary conductance (factor of two) with a two order of magnitude change in dislocation density. However, as discussed by Costescu et al. [29], there is essentially no difference in thermal boundary conductance between TiN/MgO (001), TiN/MgO (111), and TiN/Al 2 O 3 , which has varying levels of dislocation densities. is implies that either the dislocation densities, although varying, were extremely low at these TiNbased interfaces; however, without exact characterization of the defect densities at these interfaces, quantitative analyses are not possible. erefore, to elucidate the underlying role of dislocations on phonon thermal boundary conductance, controlled experiments of thin �lms on substrates with varying degrees of dislocation formation at the �lm substrate boundaries need to be pursued.

Effect of Bonding on Phonon Thermal
Boundary Conductance us far, I have discussed the roles that imperfections in the lattice structure or order of the masses around solid interfaces have on thermal boundary conductance. However, phonon transport is also affected by the strength of the interatomic bonding. e question then arises: what effect does interfacial bond strength have on phonon thermal boundary conductance? Very recently, this has been an active area of research in nanoscale heat transport. Several works have investigated this problem computationally or theoretically [92,93,[138][139][140][141][142]. Experimentally, the effects of bond strength on thermal boundary conductance have been pin-pointed in a few spe-ci�c systems using self-assembled monolayers, diamond anvil technologies, or surface termination to tune the interfacial interactions. Hsieh et al. [143] used a diamond anvil cell to vary the pressure at Al/SiC interfaces and demonstrated an increase in thermal boundary conductance with pressure when an acoustically so layer is placed between the Al and SiC (such as graphene, which interacts with the Al via van der Waal interactions, or SiO , which has as low acoustic impedance). e stiffness of these "weak" interfaces increased with pressure, and a nearly order of magnitude increase in thermal boundary conductance between the Al and SiC was demonstrated at pressures of ∼10 GPa compared to no applied pressure. Losego et al. [144] tuned the bonding between Au and quartz by using self-assembled monolayers (SAMs) sandwiched at the interface. e end group chemistries of the SAMs covalently bonded to the quartz gave a method of tuning the SAM/Au interaction strength. is led to tunability in the thermal transport from the Au to the quartz by increasing the density of covalent bonds that bridge the Au to the substrate. Collins and Chen [145] and Monachon and Weber [146] demonstrated the use of surface termination of diamond as a means to increasing the thermal boundary conductance between various metals and diamond substrates. Most notably, they found that oxygen termination can increase the thermal boundary conductance between Al and diamond by a factor of two (Collins and Chen [145]) and Cr and diamond by a factor of 11 (Monachon and Weber [146]) at room temperature as compared to untreated diamond. We investigated this problem by utilizing adsorbed elemental functional groups at graphene contacts [24]. As previously mentioned, graphene interfaces are typically described as weakly interacting due to the van der Waal forces that are characteristic of cross plane graphene/graphite. Using low pressure electron beam generated plasmas [147,148], a series of graphene �lms on SiO 2 /Si substrates were functionalized with O 2 . Aer plasma functionalization, the graphene surfaces were coated with Al and we measured the thermal boundary conductance between the Al and graphene. I show the thermal boundary conductance between Al/single layer graphene (SLG) and Al/oxygen functionalized SLG (Al/O-SLG) in Figure 10. We observed a factor of two increase in ℎ K across the Al/O-SLG interface as compared with the Al/SLG interface. We attributed this behavior to the enhanced reactivity of the oxygenated surface, which is consistent with the increase observed across the Al/O-diamond interface as compared to the Al/diamond interface [145,149], also shown in Figure 10. e presence of oxygen-containing functional groups increases the surface energy and thus makes graphene more reactive, allowing a lower energy bonding state with the evaporated Al �lm as compared to the untreated SLG which is chemically inert. is lower energy bonding state promotes Al-O formation, leading to a stronger Al-O-SLG junction and thereby an increase in ℎ K .
To quantify this, we turn to the DMM, discussed in Section 4. However, in this case, the transmission coefficient must be rederived to account for the two-dimensional nature of the graphene. To calculate the transmission coefficient, we used the assumption that we discussed in a previous work [110], in which we treat the graphene sheet as a twodimensional solid, thereby performing detailed balance on F 10: ermal boundary conductance across our Al/SLG (�lled triangles) and Al/O-SLG (un�lled triangles) interfaces [24]. e thermal transport across the Al/SLG interface increases by a factor of two with oxygen functionalization of the graphene. Our measured values are in good agreement with conductances measured across Al/diamond and Al/O-diamond interfaces, which exhibit an increase in thermal boundary conductance with oxygen termination [145,149]. We model the thermal conductance across the Al/SLG interface with the DMM, as shown by the solid line. We adjust the velocity of the SLG in our DMM calculations to model the Al/O-SLG thermal boundary conductance (dashed line), and �nd that the resultant velocity in the oxygenated SLG is about a factor of two higher than the non-functionalized sample. is is indicative of the increase in bond strength between the Al and the SLG via the oxygen adsorbates. the �uxes in the Al and SLG. For the �ux in the Al, we can make the usual three-dimensional isotropic approximation of phonon �ux discussed in Section 4. However, given the two-dimensional nature of the SLG, the phonon �ux in this material is given by [110] SLG = 1 2 0 ℏ 2D SLG SLG (13) where is the interlayer spacing of graphite (i.e., = 3.35 Å) [150]. For calculations of this anisotropic SLG-modi�ed DMM (which we refer to in Figure 10 as "DMM-2D variation"), we �t polynomials to the phonon dispersions of Al in the Γ → X direction [126] and of graphene in the Γ → K direction [151] to give greater accuracy in the DMM calculations as compared to the Debye approximation [106].
We note, however, that the modes in the Γ → K direction in graphene are extremely Debye-like (i.e., non dispersive) in the frequency regime that is elastically accessible to the Al (frequencies below 10 THz). We ignore the ZA mode in the graphene dispersion since the ZA mode is heavily suppressed in encased and supported SLG [43,152]. Given that our measurements represented the cross plane conductance of Al/SLG, the cross plane velocity of SLG is meaningless. erefore, we de�ne SLG as a �tting parameter that we assume is constant with frequency. is is a similar approach as performed in the analysis of Koh et al. [153] in which they adjust the transmission coe�cient to �t the DMM to their Au/Ti/SLG data and Schmidt et al. [154] in which they adjust the phonon transmission to �t the DMM to their metal/graphite data. We note that this approach of adjusting only SLG gives us more direct insight into how the bonding at the Al/SLG interface changes due to the functionalization since we are not adjusting any aspect of the graphene dispersion, only the transport velocity. e �ts of the DMM are shown in Figure 10. e solid line is the DMM calculation assuming the graphene cross plane velocities are 2,455 and 1,480 m s −1 ) for the longitudinal and transverse modes of the SLG, respectively (labeled "DMM-2D variation"). e dashed line is the DMM calculation assuming 4,687 and 2,825 m s −1 for the longitudinal and transverse modes of the O-SLG, respectively. e velocities that result in the best �t of the DMM to the SLG data are in good agreement with the cross plane velocity of bulk graphite [150], which could be indicative of similar bonding between the Al and SLG as the van der Waals bonds cross plane in graphite. e velocities resulting in the best �t in the Al/O-SLG data are nearly a factor of two higher, which we attribute to the increased bond strength between the Al and SLG due to the presence of oxygen leading to covalent bonding between the Al and O-SLG. It is important to note that we use the velocity parameter and �tting only to lend insight into the nature of the bonding at the two different interfaces. Since the velocity is proportional to the square root of the spring constant, the increase by a factor of two in the best-�t velocity with oxygen functionalization implies that increase in thermal transport is due to an increase in interfacial bond strength by roughly 40%. In summary, in this section I have discussed the role of interfacial bonding on thermal boundary conductance at solid interfaces. is has been a topic of several recent works, as outlined in detailed in this section. e strength of the interfacial bond can be manipulated to directly tune the interfacial conductance, offering a unique "knob" for thermal transport. In our speci�c work outlined in [24], we showed that the conductance of a single layer graphene interface can be manipulated by elemental adsorbates from low energy plasma functionalization. It is interesting to note that the temperature dependency of the thermal boundary conductance also changes, which could be evidence of higher order phonon interactions, such as 3-or 4-phonon conversion processes at the interface [138]. However, the exact phonon conversion mechanisms occurring at variably bonded interfaces are unclear, and more research must be pursued to elucidate the fundamental phonon interaction problem.

Summary
e bottle neck of many modern technologies and energy solutions has boiled down to thermal engineering problem on the nanoscale. As both increases and decreases in thermal conductivity are required for different applications, the ability to control and tune the thermal properties of materials is therefore the ultimate goal to improve both device efficiency and energy consumption trends. is desired thermal control in nanosystems, however, is not a trivial task. Due to the magnitudes of the thermal mean free paths approaching or overpassing typical length scales in nanomaterials (i.e., materials with length scales less than 1 micron), the thermal transport across interfaces can dictate the overall resistance in nanosystems. However, the fundamental mechanisms driving these electron and phonon interactions at nanoscale interfaces are difficult to predict and control since the thermal boundary conductance across interfaces (oen called the Kapitza conductance), ℎ K , is intimately related to the characteristics of the interface (structure, bonding, geometry, etc.) in addition to the fundamental atomistic properties of the materials comprising the interface itself. is interplay between interfacial properties and thermal boundary conductance can lead to additional carrier scattering mechanisms and phenomena that will not only change the thermal processes but can be used to systematically manipulate and control thermal transport in nanosystems.
In this paper, I have reviewed recent experimental progress in understanding this interplay between interfacial properties on the atomic scale and thermal transport across solid interfaces. I focused this discussion speci�cally on the role of interfacial nanoscale "imperfections, " such as surface roughness, compositional disorder, atomic dislocations, or interfacial bonding. In general, these interfacial imperfections lead to decreases in thermal boundary conductance, where each type of imperfection leads to different scattering mechanisms that can be used to control the thermal boundary conductance. is offers a unique avenue for controlling scattering and thermal transport in nanotechnology.