Abstract
We investigate the influence of optical phonon coupling across interfaces comprised of different materials with varying crystallographic orientations on the overall thermal boundary conductance. We show that for interfaces formed between a fcc solid and a L1_{0} solid (where L1_{0} solids exhibit alternating atomic layers in certain orientations), coupling between acoustic phonons in the fcc crystal and optical phonons on the L1_{0}side of the interface leads to a highly anisotropic thermal boundary conductance, where optical phonons can considerably enhance the conductance in a preferred crystallographic orientation of the layered solid. We attribute this in part to directionally dependent group velocities of optical phonons in the different crystallographic directions. For interfaces comprised of materials exhibiting diamond cubic crystal structures, higher conductances are observed for interfaces where there is a better overlap of acoustic phonons on either side of the interface, whereas, acoustic phonons directly coupling with high frequency optical phonons is shown to lower the overall conductance, especially at high temperatures where anharmonic interactions become important. Unique to the interfaces formed between the materials with diamond cubic crystal structures studied in this work, the presence of localized interfacial optical modes mediate thermal conductance across these interfaces.
Introduction
In most modernday nanoscale devices, thermal transport is limited by the high densities of interfaces rather than the materials that comprise the device^{1, 2}. The relation between the nanoscale temperature discontinuity, ΔT, at the interface between the materials and the impinging heat flux, Q, at the interface quantifies the thermal boundary conductance (h _{K} = Q/ΔT)^{3}, which is a measure of the efficacy with which heat propagates across the interface. Over the past few decades, h _{K} has been shown to be mostly limited by the acoustic impedance and mismatch of phonon spectra between the materials comprising the interface, to a first approximation^{4,5,6,7,8,9,10,11}. Moreover, prior works have also demonstrated the influence of extrinsic factors such as defect concentration and roughness around the interface^{2, 6, 12,13,14,15,16}, strength of crossspecies interaction and chemistry around the interface^{17,18,19,20,21,22} and the relative crystallographic orientations of the two materials comprising the interface^{23,24,25} in dictating interfacial heat flow. However, despite these advances in understanding interfacial heat flow, relatively few studies have focused on understanding the role of optical phonons on thermal boundary conductance at solidsolid interfaces, even though optical phonons can make up to 90% of the available vibrational modes in some materials.
Historically, theoretical calculations have suggested that optical phonons do not contribute to dominant mean free paths and resulting thermal conductivity of crystalline solids due to their relatively small group velocities and short lifetimes relative to acoustic phonons^{26, 27}. The advent of rigorous computational approaches such as first principles calculations, molecular dynamics simulations and solutions to the phonon Boltzmann transport equation, however, have highlighted the important contributions of optical modes on the thermal conductivity of bulk materials and nanostructures^{28,29,30,31,32,33,34}. In this regard, surely, the understanding of the role of optical phonons on thermal boundary conductance across solidsolid interfaces will be crucial in designing and optimizing thermal transport across nanoscale devices with high densities of interfaces that are composed of semiconducting materials with optical branches making up a significant proportion of the dispersion relations. As such, theoretical models based on the diffuse mismatch theory have been reformulated to account for optical phonons, which can contribute significantly to the total thermal boundary conductance across interfaces of several representative material systems due to elastic and inelastic scattering^{35, 36}.
In this work, we use molecular dynamics simulations and lattice dynamics calculations to investigate the role of optical phonons on thermal boundary conductance across interfaces comprised of materials with different crystallographic configurations. In particular, we study interfaces formed: 1) between a homogeneous fcc crystal and a layered L1_{0} crystal described by 6–12 LennardJones (LJ) potentials; and 2) between a homogeneous diamond cubic crystal and a layered zincblende crystal described by the StillingerWeber (SW) potentials. We show that for interfaces formed between a fcc solid and a layered L1_{0} solid, optical phonons from the L1_{0} crystal preferentially couple to acoustic phonons of the monatomic crystal when the interface is comprised of the layered crystal oriented in the direction which accommodates for the higher group velocities of the optical modes (which is the inplane direction for the layered crystal). A better overlap between these optical modes (in the inplane direction) and the acoustic modes in the homogeneous crystal is shown to dramatically enhance the thermal boundary conductance. In contrast, interfaces comprised between diamond cubic and zincblende crystal structures demonstrate higher conductances for interfaces formed with a better overlap of acoustic phonons on either side of the interface. Moreover, at higher temperatures where anharmonic interactions become important, acoustic phonons directly coupling with high frequency optical phonons is shown to lower the overall conductance for these structures. Furthermore, the presence of localized optical modes at the interface mediate thermal conductance across these interfaces, which are unique to the diamond cubic crystal structures studied in this work.
Methods and computational details
In a crystalline solid, a multiatom basis facilitates the presence of optical branches in the available vibrational modes; there are only three acoustic branches, while 3 p3 branches of the optical variety populate the vibrational spectrum, where p denotes the number of atoms forming the basis^{37}. In order to accommodate for optical phonons in the dispersion relation, we create a L1_{0} crystal with a fouratom basis (with two atom types occupying the basis sites that form the unit cell with a lattice parameter, a _{0} as shown in the schematic in Fig. 1a. The two atom types, which we refer to as material A and material B, are differentiated by mass only. The computational domains with a two point basis forms an ordered crystal that is arranged in a 1 × 1 superlatticetype structure in the[001] direction, which we refer to as the crossplane direction (cp), whereas the[010] direction is referred to as the inplane direction (ip); we refer to the L1_{0} crystal formed by the A and B atoms as an AB alloy. If the four basis atoms forming the unit cell in the L1_{0} structure are identical, a homogeneous fcc crystal is formed.
For the cubic fcc and L1_{0} crystals, the interactomic potential is defined by the widely used LJ potential, \(U(r)=4\varepsilon [{(\sigma /r)}^{12}{(\sigma /r)}^{6}]\), where U is the interatomic potential, r is the interatomic separation, and σ and ε are the LJ length and energy parameters, respectively. For computational efficiency, the cutoff distance is set to 2.5 σ for all the simulations. The length and energy parameters are modeled for LJ argon (σ = 3.405 Å and ε = 0.0103 eV, repectively) with the lattice constant a _{0} = 1.56 σ. The sizes of the computational domains are 10 a _{0} × 10 a _{0} × 80 a _{0} with periodic boundary conditions applied in the x and ydirections, whereas, fixed boundaries with 4 monolayers of atoms at each end are placed in the zdirection. Figure 1b and c show the schematics of interfaces formed between the AB alloy and materials A and B in the inplane and crossplane directions, respectively.
Similarly, for the diamond cubic crystal structures, a layered domain (accommodating optical phonons and a band gap between the optical and acoustic phonons) is constructed by replacing the four basis atoms that form the unit cell with atoms of different mass (which forms a zincblende structure as shown in Fig. 1d). For the zincblende and diamond structures, all interatomic interactions are defined by the Stillinger Weber (SW) potential^{38} for Si and only the mass of the atoms are altered for simplicity. We note that as the main purpose of this work is to understand the general effect of optical modes on the thermal boundary conductance as opposed to predicting material specific properties, the use of the LJ potential for the fcc and L1_{0} crystals and SW potential for the diamond cubic and zincblende crystals will provide sufficient qualitative insights into the role of optical phonons on thermal boundary conductance across these structures.
To begin, we prescribe the masses of materials A and B to 40 g mol^{−1} (m _{A} = m _{Ar}) and 160 g mol^{−1} (m _{B} = 4 m _{Ar}), respectively, for the LJbased crystals. For the thermal boundary conductance (h _{K}) calculations, the computational domains are separated midway to create atomically smooth and perfectly matched interfaces with the difference in mass between the materials creating the acoustic mismatch. The computational domains for the AB alloy oriented in the crossplane or inplane directions in contact with the cubic A or B materials are shown in Fig. 2a. After the domains are setup, an equilibration scheme is imposed to generate the relaxed structures. Starting with the NoseHoover thermostat^{39}, the number of atoms, volume and temperature of the simulation is held constant followed by the NPT integration (which is the isothermalisobaric ensemble with the number of particles, pressure and temperature of the system held constant) for another 1 ns at 0 bar pressure. The time step for all simulations is set to 1 fs throughout the simulations.
After equilibration, a fixed amount of energy is added per time step to a warm bath at one end and removing the equal amount of energy from a cool bath at the other end (with the lengths of the baths at 10 a _{0} in the zdirection) under a microcanonical ensemble with the number of particles, volume and energy held constant. This procedure establishes a steadystate temperature gradient and the temperature of atoms in each atomic monolayer in the direction of the applied heat flux are averaged for a total of 20 ns. The temperature profiles for the first 5 ns are discarded in order for the systems to reach steady state after which the timeaveraged temperature profiles are used to calculate the thermal boundary conductances. As NEMD simulations are known to exhibit size effects due to scattering of phonons with the heat baths^{40, 41}, we have simulated all systems with two different sizes along the direction of the applied heat flux to make sure that our MDpredicted thermal boundary conductances are not artifacts of the simulated system sizes. Figure 2b shows the temperature profiles obtained for AB(cp)/A structures considering two different system sizes (with energy added and removed at a rate of 12.5 meV ps^{−1} from the hot and cold baths, respectively). Along with these temperature profiles, we also include the temperature profile obtained by applying a lower applied heat flux (with addition and removal of energy at a rate of 6.25 meV ps^{−1} from the heat baths). As is clear from the figure, neither the computational domain size nor the applied heat flux alter h _{K}, which suggests that the results are independent of the system size as well the applied heat flux for our simulations. Note, to reduce the uncertainty in determining the temperature drop at the interface, linear fits to the temperature gradients are performed for the two materials on either side of the interface (excluding the nonlinear regimes near the boundaries and that comprising the heat baths).
The phonon dispersion relations of material A, B and the AB alloy (chosen in the direction of high symmetry) are determined via harmonic lattice dynamics calculations carried out with the General Utility Lattice Program (GULP)^{42}. For the isotropic materials A and B, the Brillouin zone (represented by the truncated octahedron) is analyzed in the[001] direction, which corresponds to Γ → X. Note, the dispersions in the [001], [010] and [100] directions are the same for the isotropic fcc solid. For the AB alloy, we use the conventional unit cell with four atom basis and a simple cubic lattice to specify the unit cell. This description results in a simple cubic Brillouin zone with a length of π/a _{0}, where a _{0} is the lattice constant. As the AB alloy are oriented in the [001] and[010] directions to form interfaces with the isotropic solid (as shown in Fig. 1b and c), the dispersion relations are analyzed in these two directions of high symmetry. The bulk phonon density of states are calculated by taking the Fourier transform of the velocity autocorrelation function and using the standard procedure for power spectral density calculations^{43}. For each atom under consideration, a velocity fluctuation time series is obtained at 10 fs intervals for a total of 500 ps under the NVE integration at 30 K after the equilibration process. The phonon density of states for material A, material B and the AB (L1_{0}) alloy in the crossplane and inplane directions are also plotted in Fig. 3 along side their respective phonon dispersion relations.
Results and Discussions
LennardJones based structures
The wavevectors are normalized by 2 π/a _{0} (k ^{*} = k/(2π/a _{0})) for all phonon dispersion relations shown in Fig. 3. For the LJbased monatomic crystals, although the Brillouin zone extends to k ^{*} value of unity, the branches have been folded at k ^{*} = 0.5 to be consistent with the dispersions for the L1_{0} AB alloy. For comparison, we have also plotted the dispersion relation for the material A with the four atom conventional unit cell (similar to the calculations for the AB alloy) in the [001] direction. The additional dispersion curves (denoted by the blue diamond symbols in Fig. 3a) correspond to dispersion curves along other directions in the octahedron Brillouin zone (defined by the primitive unit cell) that arise due to the change in the shape of the Brillouin zone and are therefore mapped in the[001] direction of the conventional unit cell^{44}. The phonon dispersion and DOS calculations for materials A and B display characteristics typical to fcc solids, whereas, the dispersions of the L1_{0} AB alloy in the inplane and crossplane directions are markedly different from the monatomic (isotropic) structures. The DOS for the AB alloy is also distinct from that of materials A and B with the formation of stop bands arising between frequencies of ~0.5 and ~1.3 THz, which is attributed to the two atom basis in the unit cell of the AB alloy (see Fig. 1). The acoustic phonons in the AB alloy extend up to 0.8 THz, and therefore the vibrational spectrum of material B with a cutoff frequency of f = 1 THz demonstrates a better spectral overlap with the acoustic phonons of the AB alloy. The vibrational spectra of Material A (which extends till ~2 THz) demonstrates an overall better spectral overlap with the AB alloy in terms of their cutoff frequencies since the optical phonons in the AB alloy extend from ~1 to 2 THz. Also, the group velocities, \({v}_{{\rm{g}}}=\partial \omega /\partial k\) of the optical phonons in the inplane directions are higher than that in the crossplane direction, the consequences of which will be discussed in terms of the thermal boundary conductance in the following paragraphs.
Figure 4 shows the thermal boundary conductances as a function of temperature for interfaces comprised of the AB alloy in the inplane and crossplane directions with materials A and B. For comparison, h _{K} across the material A/material B interface is also plotted. The uncertainties in the MDpredicted h _{K} are calculated from three independent simulations and taking into consideration the uncorrelated fluctuations in the data sampling. For the A/B interface, we have considered the effect of different crystallographic orientations on h _{K} and found no correlation between the different crystal orientations within uncertainties of the MDpredictions. This result is in agreement with previous simulations and experimental results on cubic crystals where h _{K} does not depend on the relative orientation^{23, 24, 45}.
For the interfaces comprised of the L1_{0} AB alloys, h _{K} demonstrates clear anisotropy with the inplane direction showing considerably higher conductances as compared to the crossplane direction for both AB/A and AB/B interfaces. In the crossplane direction, the conductance across the AB/B interface is higher than the AB/A interface, even though material A has a better overlap in the DOS with the AB alloy as noted earlier. However, as material B has a better overlap with the acoustic phonons in the AB alloy, h _{K} across the interface comprising of AB alloy oriented in the crossplane direction and material B is greater in comparison to that across the AB(cp)/A interface. This points to the fact that a better relative match between the cutoff frequencies (or Debye temperatures) of two materials does not necessarily result in a high conductance when the two materials have different crystallographic configurations and corresponding phonon dispersions. While the crossplane conductances are higher for the case with the better match in the acoustic modes (across the AB(cp)/B interface), the conductances for the AB(ip)/A interface with better overall spectral overlap are relatively higher at lower temperatures. At higher temperatures (>20 K), the linear increase in h _{K} for the AB(ip)/B interface leads to higher conductances as compared to that for the AB(ip)/A interface (where the conductance at higher temperatures studied in this work are similar within uncertainties).
The temperature dependence for the AB(ip)/A interface is unique as all other interfaces considered show a linear increase in h _{K} (c.f. Fig. 4). Along these lines, the A/B interface demonstrates a linear increase in h _{K}, which is in agreement with the results from prior studies that have considered LJ argon/“heavyargon” interfaces (refs 46 and 47). These prior studies have attributed the linear increase in h _{K} with temperature to anharmonic effects in the bulk and at the interface, which enhances transmission of phonon frequencies that exceed the cutoff frequency of the heavier solid. Until recently, it was largely assumed that the bulk properties of the two materials adjoining the interface (such as the spectral overlap between the two materials) could sufficiently describe phonon transport across interfaces. However, Gordiz and Henry have shown that the vibrational properties of the interface (which are dictated by the two materials that form the interface) need to be considered in order to properly describe the energy transport mechanisms occurring at the interface^{48, 49}.
Generally, the linear increase in h _{K} with temperature is related to multiplephonon inelastic scattering events at the interface^{3, 50, 51}. This would suggest that the ability of inelastic channels to increase interfacial conductance for our AB(ip)/A interface saturates at ~20 K for the temperature range studied in this work. Although, the optical modes in the AB alloy are “elastically accessible” to material A due to the DOS overlap in both cases when the AB alloy is oriented inplane or in the crossplane directions, one of the reasons for the difference in the large anisotropy in h _{K} might potentially be due to higher group velocities of the optical modes in the inplane direction as compared to that in the crossplane direction.
To understand the differing temperature trends in more detail, we calculate the DOS of atoms in the bilayer of the AB alloy that are immediately adjacent to the interface for the two orientations (c.f., Fig. 5a and b). Even though the DOS are calculated for bilayers directly adjacent to different materials in two different orientations, a drastic deviation from the bulk DOS is not observed for all cases. However, for the AB(ip)/A case, a mode depletion is observed at the frequency interval of 1.25 to 1.75 THz. This frequency range corresponds to the optical phonon modes that demonstrate higher group velocities in the inplane direction as compared to that in the crossplane direction as noted earlier. This is quantitatively depicted in Fig. 5c, which shows the calculations of group velocities as a function of frequency in the inplane and crossplane directions. The mode depletion that is observed in only one direction (Fig. 5b) suggests that the acoustic phonons in material A are coupled to the optical phonons in the AB alloy oriented in the inplane direction more so than the crossplane direction. This is because the phonon energy of these modes are distributed between the two solids, which leads to the diminished DOS of these modes in the bilayer of the AB alloy adjacent to the interface.
To quantitatively support the observations presented in the above paragraph, we turn to the spectral decomposition method at the interface (the details of which are presented in our previous work)^{20} to investigate which modes carry the significant amount of heat across the interface. Briefly, the heat flux is spectrally resolved by the relation^{46}, \(Q={\int }_{0}^{\infty }\frac{d\omega }{2\pi }q(\omega )\), where ω is the angular frequency and q(ω) is the spectral heat current. For pairwise interactions between an atom i and j, the heat current is proportional to the correlation between the interatomic force \({\vec{F}}_{ij}\) between the atoms and the velocities, \({q}_{i\to j}(\omega )\propto \langle {\vec{F}}_{ij}\cdot ({\vec{v}}_{i}+{\vec{v}}_{j})\rangle \), where the brackets denote steadystate nonequilibrium ensemble average^{52,53,54}. To this effect, under the NVE integration at 30 K, velocities of atoms near the interface for the AB alloy are tabulated for a total of 10 ns with 10 fs time intervals along with the forces due to atoms on the other side of the interface. The resulting heat current accumulations from the AB alloy to the other side of the interfaces are shown in Fig. 5d. As is clear, the high frequency optical modes in the inplane orientation of the AB alloy contribute more towards thermal boundary conductance than in the crossplane orientation. Moreover, for the AB(ip)/A interface, the modes >1.5 THz contribute ~30% of the heat flux, whereas for the other three cases, the relative contribution to the heat current from these modes is considerably lower (<10%). This is in line with the observation of mode depletion observed for the AB(ip)/A interface (while for the other three cases, the interfacial DOS is similar to that of the bulk) as discussed in the above paragraph and shown in Fig. 5a and b.
To better understand the role of relative match of the cutoff frequencies (or Debye temperatures) on the interfacial heat flow, we performed additional simulations by varying the mass of the isotropic fcc solid while the AB alloy remains unchanged; note, increasing the mass decreases the cutoff frequency via the relation, \(\omega \propto \mathrm{1/}\sqrt{m}\). Figure 6a shows h _{K} as a function of mass of the fcc solid for interfaces comprised of inplane and crossplane configurations of the AB alloy. For the crossplane configuration, h _{K} increases monotonically up to a mass of 130 g mol^{−1} and decreases thereafter. In Fig. 6b the DOS for the solid with m = 130 g mol^{−1} is plotted along with the DOS for the AB alloy. For comparison, we also plot the DOS for the solid with m = 55 g mol^{−1}, which represents a lighter solid with a broader spectrum of frequencies that can elastically accommodate the optical phonons in the AB alloy. The acoustic modes in the heavier solid (with m = 130 g mol^{−1}) better overlap the acoustic modes in the AB alloy, which results in an increase in h _{K} when the AB alloy is oriented in the crossplane direction. This suggests that for interfaces comprised of AB alloys in the crossplane configuration, interfacial heat flow is mainly mediated through acoustic phonons. However, a further increase in the mass leads to a mismatch in the heat carrying acoustic phonons between the solids, which could potentially be one of the reasons for the decrease in h _{K}. In contrast, for the inplane configuration, the maximum in h _{K} is observed for m = 55 g mol^{−1}, which has a better overall spectral overlap. In this case, the high frequency acoustic modes in the fcc solid coincide with the high frequency optical modes in the AB alloy (in the 1.25 to 1.75 THz range). Going back to the phonon dispersion relations and group velocities shown in Figs 3b,c and 5c, the optical modes in the inplane configuration have a higher group velocity as compared to that in the crossplane direction and therefore carry a more significant proportion of heat. This manifests in a higher h _{K} when the optical modes are “elastically accessible” to the acoustic modes in the fcc solid across the interface, as noted earlier. It is also worth noting that the anisotropy in h _{K} is pronounced for comparatively lighter masses with broader frequency spectrums that can elastically couple with the optical phonons in the AB alloy.
StillingerWeber based structures
As most of the semiconductor devices are based on Si, we also investigate the role of optical phonons across interfaces comprised of diamond cubic crystal structures defined by the SW potential. The AB alloy in this case is formed from a combination of Si and massheavy Si with m = 72.6 g mol^{−1} (which we refer to as Ge) arranged in a 1 × 1 layered structure with the atoms placed in a fcc lattice with a 2 point basis with a_{0} = 5.44 Å, which we have referred to as the zincblende structure. Similar to the LJ structures, we vary the mass of the homogeneous crystal (forming the interface with the SiGe alloy) to gauge the effect of DOS overlap (of the optical phonons) between the two materials. Figure 7a shows the results of the NEMD calculations and Fig. 7b shows representative DOS for the SiGe alloy and the homogeneous diamond structures with m = 50 and 112 g mol^{−1}. Similar to the dependence of h _{K} on the mass of material B for the AB(cp)/B interfaces in the LJ systems, h _{K} increases linearly up to 50 g mol^{−1} and decreases thereafter. As is clear from Fig. 7b, the acoustic phonons in the 50 g mol^{−1} perfectly coincide with the acoustic phonons in the SiGe alloy. This suggests that acoustic phonons are the primary heat carriers mediating interfacial transport across these diamond structured materials. This is in line with the modal decomposition of thermal boundary conductance across Si/Ge interfaces carried out in ref. 55 where it is shown that ~80% of the heat is carried by frequencies that fall in the acoustic spectrum of Si.
Figure 8a shows the temperature dependence of h _{K} for interfaces comprised of solids with m = 112 and 20 g mol^{−1} and the SiGe alloy. Within uncertainties, h _{K} values for the two interfaces are similar for almost the entire temperature range. However, at 1400 K, h _{K} is greater for the interface comprised of m = 112 g mol^{−1} due to the relatively stronger anharmonic interactions of the acoustic phonons that are better overlapped with the acoustic phonons of the SiGe alloy as compared to that of the solid with m = 20 g mol^{−1}.
To investigate the effect of acoustic phonons in the homogeneous solid coupling with the high frequency optical phonons in the AB alloy, we run additional simulations on interfaces comprised of SiGe alloy and solid with m = 12 g mol^{−1} and an AB alloy with the two atoms defined by m = 120 g mol^{−1} and m = 160 g mol^{−1}. The band gap for the alloy with the lower mass ratio between the atoms is relatively smaller as compared to that of the SiGe (c.f., Fig. 8b and c for the DOS comparison between the alloys), allowing for the overlap of the acoustic phonons in the solid with m = 12 g mol^{−1} with the high frequency phonons in the AB alloy. The results for these interfaces are also plotted in Fig. 8a, which show that h _{K} at low temperatures (<700 K) are similar for the two cases, whereas, the increase in anharmonic interactions at the interface comprised of the SiGe alloy leads to greater conductances as compared to the case where the acoustic phonons in the homogeneous solid are coincidental with the high frequency optical phonons in the AB alloy at higher temperatures. It is also worth noting that the interfaces composed of diamond and zincblende structures demonstrate linear temperature dependencies for the temperature range studied, even when the acoustic phonons are “elastically accessible” to the high frequency optical phonons, which is unlike in the interfaces formed between the fcc and L1_{0} structures where the conductance of the AB(ip)/A interface deviates from the linear dependence (c.f., Fig. 4). Similar linear temperature trends of h _{K} for interfaces comprised of an fcc solid and GaN (with wurtzite crystal structure) were observed in ref. 22. For these structures, the better overlap of acoustic phonons of the fcc solid to the acoustic modes in the GaN was shown to enhance the conductance^{22}.
We now study the local phonon DOS immediately adjacent to the interfaces for our diamond cubic structures to better understand the results discussed in the previous paragraphs. To this end, we calculate the DOS of bilayers for the SiGe alloy and diamond structures with m = 20 g mol^{−1} and m = 112 g mol^{−1} that comprise the interface, the results of which are plotted in Fig. 9. Unlike the LJ systems, the local DOS reveals the existence of interfacial modes in two materials that are absent in the bulk DOS for each of the corresponding material. These interfacial modes are highlighted in Fig. 9. It is interesting to note that these modes are above the cutoff frequency for either the SiGe alloy (for the interface comprised of SiGe and diamond structure with m = 20 g mol^{−1}) or the homogeneous diamond structure (for the interface comprised of SiGe and diamond structure with m = 112 g mol^{−1}). This suggests that the contributions to h _{K} from these modes are a direct result of inelastic energy exchange due to the anharmonicity at the interface. Similar results for interfaces comprised of Si and Ge have also been reported, where considerable amount of heat is carried by interfacial modes at 12–13 THz region that are beyond the cutoff frequency of Ge^{55,56,57}. The large DOS for the interfacial modes is also not an artifact of the SW potential used in this work as similar results (for DOS calculations as shown Fig. 9) are obtained when the atomic interactions are defined by the Tersoff potential^{58, 59}. These results suggest that the description of heat transfer across these interfaces cannot be accurately depicted by only considering the bulk phonon properties of the two materials and it becomes necessary to consider the localized and nondispersive interfacial modes to accurately describe interfacial heat transfer across these diamond cubic structures. Moreover, Gordiz and Henry postulate that these interfacial modes facilitate interfacial energy transport among other phonon modes since the conductance decreases considerably when these interfacial modes are excluded (see Fig. 5 of ref. 55). For a detailed investigation of the role of these interfacial modes on heat transfer across Si/Ge interfaces, the reader is referred to ref. 55. We note that the fcc structures defined by the LJ potential did not show similar existence of pronounced interfacial modes that mediate thermal transport. Moreover, in ref. 60, an interface comprised of InGaAs/InP (note, InGaAs and InP both exhibit zincblende crystal structures) was also shown to exhibit little to no interfacial modes that mediate thermal transport across these materials. This suggests that the existence of large interfacial modes could be unique to interfaces composed of diamond cubic crystal structures of Si and Ge.
Conclusion
We assessed the role of optical phonons on the thermal boundary conductance across materials with different crystal structures and orientations. We show that for interfaces formed between a fcc solid and a L1_{0} solid (where the L1_{0} solid exhibits alternating atomic layers in certain orientations), optical phonons are preferentially better coupled for the case when the interface is comprised of the layered crystal oriented in the direction which possesses higher group velocities of these modes, thus increasing the overall thermal boundary conductance. In contrast, interfaces comprised of diamond cubic crystal structures demonstrate higher conductances for interfaces formed with a better overlap of acoustic phonons on either side. At high temperatures where anharmonic interactions become important, acoustic phonons directly coupling with high frequency optical phonons is shown to lower the overall conductance for these structures, which is in contrast to the results obtained for the LJbased solids. Furthermore, the presence of localized interfacial optical modes, which in this work are unique to the diamond cubic crystal structures are shown to mediate thermal conductance across these interfaces.
References
 1.
Pop, E. Energy dissipation and transport in nanoscale devices. Nano Research 3, 147–169 (2010).
 2.
Hopkins, P. E. Thermal transport across solid interfaces with nanoscale imperfections: effects of roughness, disorder, dislocations, and bonding on thermal boundary conductance. ISRN Mechanical Engineering 682586 (2013).
 3.
Swartz, E. T. & Pohl, R. O. Thermal boundary resistance. Reviews of Modern Physics 61, 605–668 (1989).
 4.
Stoner, R. J. & Maris, H. J. Kapitza conductance and heat flow between solids at temperatures from 50 to 300 K. Phys. Rev. B 48, 16373–16387 (1993).
 5.
Stoner, R. J., Maris, H. J., Anthony, T. R. & Banholzer, W. F. Measurements of the kapitza conductance between diamond and several metals. Phys. Rev. Lett. 68, 1563–1566 (1992).
 6.
Stevens, R. J., Smith, A. N. & Norris, P. M. Measurement of thermal boundary conductance of a series of metaldielectric interfaces by the transient thermoreflectance technique. Journal of Heat Transfer 127, 315–322 (2005).
 7.
Hopkins, P. E., Stevens, R. J. & Norris, P. M. Influence of inelastic scattering at metaldielectric interfaces. Journal of Heat Transfer 130, 022401–022401 (2008).
 8.
Lyeo, H.K. & Cahill, D. G. Thermal conductance of interfaces between highly dissimilar materials. Phys. Rev. B 73, 144301 (2006).
 9.
Freedman, J. P., Yu, X., Davis, R. F., Gellman, A. J. & Malen, J. A. Thermal interface conductance across metal alloy–dielectric interfaces. Phys. Rev. B 93, 035309 (2016).
 10.
Cheaito, R. et al. Thermal boundary conductance accumulation and interfacial phonon transmission: Measurements and theory. Phys. Rev. B 91, 035432 (2015).
 11.
Dechaumphai, E. et al. Ultralow thermal conductivity of multilayers with highly dissimilar debye temperatures. Nano Letters 14, 2448–2455 (2014).
 12.
Stevens, R. J., Zhigilei, L. V. & Norris, P. M. Effects of temperature and disorder on thermal boundary conductance at solid–solid interfaces: Nonequilibrium molecular dynamics simulations. International Journal of Heat and Mass Transfer 50, 3977–3989 (2007).
 13.
Giri, A., Hopkins, P. E., Wessel, J. G. & Duda, J. C. Kapitza resistance and the thermal conductivity of amorphous superlattices. Journal of Applied Physics 118, 165303 (2015).
 14.
Giri, A., Braun, J. L. & Hopkins, P. E. Effect of crystalline/amorphous interfaces on thermal transport across confined thin films and superlattices. Journal of Applied Physics 119 (2016).
 15.
Merabia, S. & Termentzidis, K. Thermal boundary conductance across rough interfaces probed by molecular dynamics. Phys. Rev. B 89, 054309 (2014).
 16.
Zhou, X. W., Jones, R. E., Kimmer, C. J., Duda, J. C. & Hopkins, P. E. Relationship of thermal boundary conductance to structure from an analytical model plus molecular dynamics simulations. Phys. Rev. B 87, 094303 (2013).
 17.
Duda, J. C. et al. Influence of interfacial properties on thermal transport at gold:silicon contacts. Applied Physics Letters 102, 081902 (2013).
 18.
Losego, M. D., Grady, M. E., Sottos, N. R., Cahill, D. G. & Braun, P. V. Effects of chemical bonding on heat transport across interfaces. Nat Mater 11, 502–506 (2012).
 19.
Giri, A. et al. Mechanisms of nonequilibrium electronphonon coupling and thermal conductance at interfaces. Journal of Applied Physics 117, 105105 (2015).
 20.
Giri, A., Braun, J. L. & Hopkins, P. E. Implications of interfacial bond strength on the spectral contributions to thermal boundary conductance across solid, liquid, and gas interfaces: A molecular dynamics study. The Journal of Physical Chemistry C 120, 24847–24856 (2016).
 21.
Duda, J. C. et al. Implications of crossspecies interactions on the temperature dependence of kapitza conductance. Phys. Rev. B 84, 193301 (2011).
 22.
Zhou, X. W., Jones, R. E., Duda, J. C. & Hopkins, P. E. Molecular dynamics studies of material property effects on thermal boundary conductance. Phys. Chem. Chem. Phys. 15, 11078–11087 (2013).
 23.
Duda, J. C. et al. Influence of crystallographic orientation and anisotropy on kapitza conductance via classical molecular dynamics simulations. Journal of Applied Physics 112 (2012).
 24.
Hopkins, P. E. et al. Influence of anisotropy on thermal boundary conductance at solid interfaces. Phys. Rev. B 84, 125408 (2011).
 25.
Duda, J. C., Smoyer, J. L., Norris, P. M. & Hopkins, P. E. Extension of the diffuse mismatch model for thermal boundary conductance between isotropic and anisotropic materials. Applied Physics Letters 95 (2009).
 26.
Srivastava, G. Phonon conductivity of insulators and semiconductors. Journal of Physics and Chemistry of Solids 41, 357–368 (1980).
 27.
Chen, G. Thermal conductivity and ballisticphonon transport in the crossplane direction of superlattices. Phys. Rev. B 57, 14958–14973 (1998).
 28.
Tian, Z., Esfarjani, K., Shiomi, J., Henry, A. S. & Chen, G. On the importance of optical phonons to thermal conductivity in nanostructures. Applied Physics Letters 99, 053122 (2011).
 29.
Ward, A. & Broido, D. A. Intrinsic phonon relaxation times from firstprinciples studies of the thermal conductivities of si and ge. Phys. Rev. B 81, 085205 (2010).
 30.
Henry, A. S. & Chen, G. Spectral phonon transport properties of silicon based on molecular dynamics simulations and lattice dynamics. Journal of Computational and Theoretical Nanoscience 5, 141–152 (2008).
 31.
Sellan, D. P., Turney, J. E., McGaughey, A. J. H. & Amon, C. H. Crossplane phonon transport in thin films. Journal of Applied Physics 108, 113524 (2010).
 32.
Jain, A. & McGaughey, A. J. H. Thermal conductivity of compound semiconductors: Interplay of mass density and acousticoptical phonon frequency gap. Journal of Applied Physics 116 (2014).
 33.
Shulumba, N., Hellman, O. & Minnich, A. J. Intrinsic localized mode and low thermal conductivity of pbse. Phys. Rev. B 95, 014302 (2017).
 34.
McGaughey, A. J. H., Hussein, M. I., Landry, E. S., Kaviany, M. & Hulbert, G. M. Phonon band structure and thermal transport correlation in a layered diatomic crystal. Phys. Rev. B 74, 104304 (2006).
 35.
Duda, J. C., Beechem, T. E., Smoyer, J. L., Norris, P. M. & Hopkins, P. E. Role of dispersion on phononic thermal boundary conductance. Journal of Applied Physics 108, 073515 (2010).
 36.
Beechem, T., Duda, J. C., Hopkins, P. E. & Norris, P. M. Contribution of optical phonons to thermal boundary conductance. Applied Physics Letters 97, 061907 (2010).
 37.
Kittel, C. Introduction to Solid State Physics. 6th edn, (John Wiley & Sons, Inc, New York, 1986).
 38.
Stillinger, F. H. & Weber, T. A. Computer simulation of local order in condensed phases of silicon. Phys. Rev. B 31, 5262–5271 (1985).
 39.
Hoover, W. G. Canonical dynamics: Equilibrium phasespace distributions. Phys. Rev. A 31, 1695–1697 (1985).
 40.
Sellan, D. P., Landry, E. S., Turney, J. E., McGaughey, A. J. H. & Amon, C. H. Size effects in molecular dynamics thermal conductivity predictions. Phys. Rev. B 81, 214305 (2010).
 41.
Evans, W. J., Hu, L. & Keblinski, P. Thermal conductivity of graphene ribbons from equilibrium molecular dynamics: Effect of ribbon width, edge roughness, and hydrogen termination. Applied Physics Letters 96, 203112 (2010).
 42.
Gale, J. D. & Rohl, A. L. The general utility lattice program (gulp). Molecular Simulation 29, 291–341 (2003).
 43.
Allen, M. P. & Tildesley, D. J. Computer simulation of liquids (Oxford Science Publications). Oxford science publications (Oxford University Press, 1989), reprint edn.
 44.
Turney, J. E. Predicting phonon properties and thermal conductivity from anharmonic lattice dynamics calculations. Ph.D. thesis, Carnegie Mellon University, Pittsburgh, Pennsylvania (2009).
 45.
Costescu, R. M., Wall, M. A. & Cahill, D. G. Thermal conductance of epitaxial interfaces. Phys. Rev. B 67, 054302 (2003).
 46.
Sääskilahti, K., Oksanen, J., Tulkki, J. & Volz, S. Role of anharmonic phonon scattering in the spectrally decomposed thermal conductance at planar interfaces. Phys. Rev. B 90, 134312 (2014).
 47.
Le, N. Q. et al. Effects of bulk and interfacial anharmonicity on thermal conductance at solid/solid interfaces. ArXiv eprints (2016).
 48.
Gordiz, K. & Henry, A. Phonon transport at interfaces: Determining the correct modes of vibration. Journal of Applied Physics 119 (2016).
 49.
Gordiz, K. & Henry, A. A formalism for calculating the modal contributions to thermal interface conductance. New Journal of Physics 17, 103002 (2015).
 50.
Duda, J. C., Norris, P. M. & Hopkins, P. E. On the linear temperature dependence of phonon thermal boundary conductance in the classical limit. Journal of Heat Transfer 133, 074501 (2011).
 51.
Landry, E. S. & McGaughey, A. J. H. Thermal boundary resistance predictions from molecular dynamics simulations and theoretical calculations. Phys. Rev. B 80, 165304 (2009).
 52.
Sääskilahti, K., Oksanen, J., Linna, R. P. & Tulkki, J. Thermal conduction and interface effects in nanoscale fermipastaulam conductors. Phys. Rev. E 86, 031107 (2012).
 53.
Ong, Z.Y. & Pop, E. Frequency and polarization dependence of thermal coupling between carbon nanotubes and sio_{2}. Journal of Applied Physics 108, 103502 (2010).
 54.
Domingues, G., Volz, S., Joulain, K. & Greffet, J.J. Heat transfer between two nanoparticles through near field interaction. Phys. Rev. Lett. 94, 085901 (2005).
 55.
Gordiz, K. & Henry, A. Phonon transport at crystalline si/ge interfaces: The role of interfacial modes of vibration. Scientific Reports 6, 23139 (2016).
 56.
Chalopin, Y. & Volz, S. A microscopic formulation of the phonon transmission at the nanoscale. Applied Physics Letters 103 (2013).
 57.
Murakami, T., Hori, T., Shiga, T. & Shiomi, J. Probing and tuning inelastic phonon conductance across finitethickness interface. Applied Physics Express 7, 121801 (2014).
 58.
Tersoff, J. Modeling solidstate chemistry: Interatomic potentials for multicomponent systems. Phys. Rev. B 39, 5566–5568 (1989).
 59.
Tersoff, J. New empirical approach for the structure and energy of covalent systems. Phys. Rev. B 37, 6991–7000 (1988).
 60.
Gordiz, K. & Henry, A. Interface conductance modal analysis of lattice matched ingaas/inp. Applied Physics Letters 108 (2016).
Acknowledgements
We would like to thank the support from the Office of Naval research (N000141512769).
Author information
Affiliations
Contributions
P.E.H. and A.G. conceived the project and A.G. performed the simulations. All authors wrote and reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Giri, A., Hopkins, P.E. Role of interfacial mode coupling of optical phonons on thermal boundary conductance. Sci Rep 7, 11011 (2017). https://doi.org/10.1038/s4159801710482z
Received:
Accepted:
Published:
Further reading

High thermal conductive copper/diamond composites: state of the art
Journal of Materials Science (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.