Efficient molecular dynamics simulations with many-body potentials on graphics processing units

Graphics processing units have been extensively used to accelerate classical molecular dynamics simulations. However, there is much less progress on the acceleration of force evaluations for many-body potentials compared to pairwise ones. In the conventional force evaluation algorithm for many-body potentials, the force, virial stress, and heat current for a given atom are accumulated within different loops, which could result in write conflict between different threads in a CUDA kernel. In this work, we provide a new force evaluation algorithm, which is based on an explicit pairwise force expression for many-body potentials derived recently (Fan et al., 2015). In our algorithm, the force, virial stress, and heat current for a given atom can be accumulated within a single thread and is free of write conflicts. We discuss the formulations and algorithms and evaluate their performance. A new open-source code, GPUMD, is developed based on the proposed formulations. For the Tersoff many-body potential, the double precision performance of GPUMD using a Tesla K40 card is equivalent to that of the LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) molecular dynamics code running with about 100 CPU cores (Intel Xeon CPU X5670 @ 2.93 GHz). (C) 2017 Elsevier B.V. All rights reserved.Peer reviewe

Magic angle in thermal conductivity of twisted bilayer graphene

We report a local minimum in thermal conductivity in twisted bilayer graphene (TBG) at the angle of 1.08$^\circ$, which corresponds to the 'magic angle' in the transition of several other reported properties. Within the supercell of a moir\'e lattice, different stacking modes generate phonon scattering sites which reduce the thermal conductivity of TBG. The thermal magic angle arises from the competition between the delocalization of atomic vibrational amplitudes and stresses on one hand, and the increased AA stacking density on the other hand. The former effect weakens the scattering strength of a single scatterer while the latter one increases the density of scatterers. The combination of these two effects eventually leads to the apparition of the highlighted irregularity in heat conduction. The manifestation of a magic angle, disclosing new thermal mechanisms at nanoscale, further uncovers the unique physics of two-dimensional materials.Comment: 15 pages, 5 figure

Energetics and structure of grain boundary triple junctions in graphene

© 2017 The Author(s). Grain boundary triple junctions are a key structural element in polycrystalline materials. They are involved in the formation of microstructures and can influence the mechanical and electronic properties of materials. In this work we study the structure and energetics of triple junctions in graphene using a multiscale modelling approach based on combining the phase field crystal approach with classical molecular dynamics simulations and quantum-mechanical density functional theory calculations. We focus on the atomic structure and formation energy of the triple junctions as a function of the misorientation between the adjacent grains. We find that the triple junctions in graphene consist mostly of five-fold and seven-fold carbon rings. Most importantly, in addition to positive triple junction formation energies we also find a significant number of orientations for which the formation energy is negative

Extremely anisotropic van der Waals thermal conductors

The densification of integrated circuits requires thermal management strategies and high thermal conductivity materials1–3. Recent innovations include the development of materials with thermal conduction anisotropy, which can remove hotspots along the fast-axis direction and provide thermal insulation along the slow axis4,5. However, most artificially engineered thermal conductors have anisotropy ratios much smaller than those seen in naturally anisotropic materials. Here we report extremely anisotropic thermal conductors based on large-area van der Waals thin films with random interlayer rotations, which produce a room-temperature thermal anisotropy ratio close to 900 in MoS2, one of the highest ever reported. This is enabled by the interlayer rotations that impede the through-plane thermal transport, while the\ua0long-range intralayer crystallinity maintains high in-plane thermal conductivity. We measure ultralow thermal conductivities in the through-plane direction for MoS2 (57 \ub1 3 mW m−1 K−1) and WS2 (41 \ub1 3 mW m−1 K−1) films, and we quantitatively explain these values using molecular dynamics simulations that reveal one-dimensional glass-like thermal transport. Conversely, the in-plane thermal conductivity in these MoS2 films is close to the single-crystal value. Covering nanofabricated gold electrodes with our anisotropic films prevents overheating of the electrodes and blocks heat from reaching the device surface. Our work establishes interlayer rotation in crystalline layered materials as a new degree of freedom for engineering-directed heat transport in solid-state systems