Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: A. Kınacı, J.B. Haskins, C. Sevik, and T. Çağın (2012), "Thermal conductivity of BN-C nanostructures", Physical Review B, 86(11), 115410. DOI: 10.1103/physrevb.86.115410.
Abstract: Chemical and structural diversity present in hexagonal boron nitride (h-BN) and graphene hybrid nanostructures provide avenues for tuning various properties for their technological applications. In this paper we investigate the variation of thermal conductivity (κ) of hybrid graphene/h-BN nanostructures: stripe superlattices and BN (graphene) dots embedded in graphene (BN) are investigated using equilibrium molecular dynamics. To simulate these systems, we have parametrized a Tersoff type interaction potential to reproduce the ab initio energetics of the B-C and N-C bonds for studying the various interfaces that emerge in these hybrid nanostructures. We demonstrate that both the details of the interface, including energetic stability and shape, as well as the spacing of the interfaces in the material, exert strong control on the thermal conductivity of these systems. For stripe superlattices, we find that zigzag configured interfaces produce a higher κ in the direction parallel to the interface than the armchair configuration, while the perpendicular conductivity is less prone to the details of the interface and is limited by the κ of h-BN. Additionally, the embedded dot structures, having mixed zigzag and armchair interfaces, affect the thermal transport properties more strongly than superlattices. The largest reduction in thermal conductivity is observed at 50% dot concentration, but the dot radius appears to have little effect on the magnitude of reduction around this concentration.