Anomalous heat transport in one-dimensional nanostructures, such as nanotubes and nanowires, is a widely debated problem in condensed matter and statistical physics, with contradicting pieces of evidence from experiments and simulations. Using a comprehensive modeling approach, comprised of lattice dynamics and molecular dynamics simulations, we proved that the infinite length limit of the thermal conductivity of a (10,0) single-wall carbon nanotube is finite but this limit is reached only for macroscopic lengths due to a thermal phonon mean free path of several millimeters. Our calculations showed that the extremely high thermal conductivity of this system at room temperature is dictated by quantum effects. Modal analysis showed that the divergent nature of thermal conductivity, observed in one-dimensional model systems, is suppressed in carbon nanotubes by anharmonic scattering channels provided by the flexural and optical modes with polarization in the plane orthogonal to the transport direction.

Anomalous heat transport in one-dimensional nanostructures, such as nanotubes and nanowires, is a widely debated problem in condensed matter and statistical physics, with contradicting pieces of evidence from experiments and simulations. Using a comprehensive modeling approach, comprised of lattice dynamics and molecular dynamics simulations, we proved that the infinite length limit of the thermal conductivity of a (10,0) single-wall carbon nanotube is finite but this limit is reached only for macroscopic lengths due to thermal phonon mean free path of several millimeters. Our calculations showed that the extremely high thermal conductivity of this system at room temperature is dictated by quantum effects. Modal analysis showed that the divergent nature of thermal conductivity, observed in one-dimensional model systems, is suppressed in carbon nanotubes by anharmonic scattering channels provided by the flexural and optical modes with polarization in the plane orthogonal to the transport direction.

We derive expressions of interatomic force and heat current for many-body potentials such as the Tersoff, the Brenner, and the Stillinger-Weber potential used extensively in molecular dynamics simulations of covalently bonded materials. Although these potentials have a many-body nature, a pairwise force expression that follows Newton's third law can be found without referring to any partition of the potential. Based on this force formula, a stress applicable for periodic systems can be unambiguously defined. The force formula can then be used to derive the heat current formulas using a natural potential partitioning. Our heat current formulation is found to be equivalent to most of the seemingly different heat current formulas used in the literature, but to deviate from the stress-based formula derived from two-body potential. We validate our formulation numerically on various systems described by the Tersoff potential, namely three-dimensional silicon and diamond, two-dimensional graphene, and quasi-one-dimensional carbon nanotube. The effects of cell size and production time used in the simulation are examined.

Grain boundaries in graphene are inherent in wafer-scale samples prepared by chemical vapor deposition. They can strongly influence the mechanical properties and electronic and heat transport in graphene. In this work, we employ extensive molecular dynamics simulations to study thermal transport in large suspended polycrystalline graphene samples. Samples of different controlled grain sizes are prepared by a recently developed efficient multiscale approach based on the phase field crystal model. In contrast to previous works, our results show that the scaling of the thermal conductivity with the grain size implies bimodal behavior with two effective Kapitza lengths. The scaling is dominated by the out-of-plane (flexural) phonons with a Kapitza length that is an order of magnitude larger than that of the in-plane phonons. We also show that, to get quantitative agreement with the most recent experiments, quantum corrections need to be applied to both the Kapitza conductance of grain boundaries and the thermal conductivity of pristine graphene, and the corresponding Kapitza lengths must be renormalized accordingly.

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here, we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green's function method in the ballistic regime and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it into the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification.

Two-dimensional materials have unusual phonon spectra due to the presence of flexural (out-of-plane) modes. Although molecular dynamics simulations have been extensively used to study heat transport in such materials, conventional formalisms treat the phonon dynamics isotropically. Here, we decompose the microscopic heat current in atomistic simulations into in-plane and out-of-plane components, corresponding to in-plane and out-of-plane phonon dynamics, respectively. This decomposition allows for direct computation of the corresponding thermal conductivity components in two-dimensional materials. We apply this decomposition to study heat transport in suspended graphene, using both equilibrium and nonequilibrium molecular dynamics simulations. We show that the flexural component is responsible for about two-thirds of the total thermal conductivity in unstrained graphene, and the acoustic flexural component is responsible for the logarithmic divergence of the conductivity when a sufficiently large tensile strain is applied.