The lattice-based scattering boundary method is applied to compute the phonon mode-resolved transmission coefficients and thermal conductances of in-plane heterostructures built from graphene and hexagonal boron nitride (hBN). The thermal conductance of all structures is dominated by acoustic phonon modes near the Brillouin zone center that have high group velocity, population, and transmission coefficient. Out-of-plane modes make their most significant contributions at low frequencies, whereas in-plane modes contribute across the frequency spectrum. Finite-length superlattice junctions between graphene and hBN leads have a lower thermal conductance than comparable junctions between two graphene leads due to lack of transmission in the hBN phonon bandgap. The thermal conductances of bilayer systems differ by less than 10% from their singlelayer counterparts on a per area basis, in contrast to the strong thermal conductivity reduction when moving from single-to multi-layer graphene. Published by AIP Publishing.

We use molecular dynamics simulations and the lattice-based scattering boundary method to compute the thermal conductance of finite-length Lennard-Jones superlattice junctions confined by bulk crystalline leads. The superlattice junction thermal conductance depends on the properties of the leads. For junctions with a superlattice period of four atomic monolayers at temperatures between 5 and 20 K, those with mass-mismatched leads have a greater thermal conductance than those with mass-matched leads. We attribute this lead effect to interference between and the ballistic transport of emergent junction vibrational modes. The lead effect diminishes when the temperature is increased, when the superlattice period is increased, and when interfacial disorder is introduced, but is reversed in the harmonic limit. C 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.

Three models of the C 60 molecular crystal are studied using molecular dynamics simulations to resolve the roles played by intermolecular and intramolecular degrees of freedom (DOF) in its structural, mechanical, and thermal properties at temperatures between 35 and 400 K. In the full DOF model, all DOF are active. In the rigid body model, the intramolecular DOF are frozen, such that only center of mass (COM) translations and molecular rotations/librations are active. In the point mass model, the molecule is replaced by a point mass, such that only COM translations are active. The zero-pressure lattice constants and bulk moduli predicted from the three models fall within ranges of 0.15 and 20%. The thermal conductivity of the point mass model is the largest across the temperature range, showing a crystal-like temperature dependence (i.e., it decreases with increasing temperature) due to the presence of phonon modes associated with the COM translations. The rigid body model thermal conductivity is the smallest and follows two distinct regimes. It is crystal-like at low temperatures and becomes temperature invariant at high temperatures. The latter is typical of the behavior of an amorphous material. By calculating the rotational diffusion coefficient, the transition between the two regimes is found to occur at the temperature where the molecules begin to rotate freely. Above this temperature, phonons related to COM translations are scattered by the rotational DOF. The full DOF model thermal conductivity is larger than that of the rigid body model, indicating that intramolecular DOF contribute to thermal transport. * S. Kumar and C. Shao contributed equally to this work. † mcgaughey@cmu.edu 1 M. S. Dresselhaus, G. Dresselhaus, and P. C. Eklund, Science of fullerenes and carbon nanotubes: their properties and applications (Academic press, 1996).

No abstract