Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations

Abstract

In this paper, we develop a highly efficient molecular dynamics code fully implemented on graphics processing units for thermal conductivity calculations using the Green-Kubo formula. We compare two different schemes for force evaluation, a previously used thread-scheme where a single thread is used for one particle and each thread calculates the total force for the corresponding particle, and a new block-scheme where a whole block is used for one particle and each thread in the block calculates one or several pair forces between the particle associated with the given block and its neighbor particle(s) associated with the given thread. For both schemes, two different classical potentials, namely, the Lennard-Jones potential and the rigid-ion potential are implemented. While the thread-scheme performs a little better for relatively large systems, the block-scheme performs much better for relatively small systems. The relative performance of the block-scheme over the thread-scheme also increases with the increasing of the cutoff radius. We validate the implementation by calculating lattice thermal conductivities of solid argon and lead telluride.

References in corpus (6)

• High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA
• Harvesting graphics power for MD simulations
• Thermal conduction of carbon nanotubes using molecular dynamics
• Enhanced molecular dynamics performance with a programmable graphics processor
• Exact diagonalization of the Hubbard model on graphics processing units
• Computation of electron quantum transport in graphene nanoribbons using GPU

Cited by in corpus (21)

• Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids
• Thermal conductivity decomposition in two-dimensional materials: Application to graphene
• Improving the accuracy of the neuroevolution machine learning potential for multi-component systems
• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles
• Multiscale modeling of polycrystalline graphene: A comparison of structure and defect energies of realistic samples from phase field crystal models
• Thermal Transport in MoS${}_2$ from Molecular Dynamics using Different Empirical Potentials
• Kapitza thermal resistance across individual grain boundaries in graphene
• Reduced Thermal Conductivity of Supported and Encased Monolayer and Bilayer MoS${}_2$
• Heat transport in pristine and polycrystalline single-layer hexagonal boron nitride
• Bimodal grain-size scaling of thermal transport in polycrystalline graphene from large-scale molecular dynamics simulations
• Efficient linear-scaling quantum transport calculations on graphics processing units and applications on electron transport in graphene
• Thermal transport properties of single-layer black phosphorous from extensive molecular dynamics simulations
• Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations
• Efficient Cross-layer Thermal Transport with Atypical Glassy-like Phenomena in Crystalline CsCu${}_4$Se${}_3$
• GPUQT: An efficient linear-scaling quantum transport code fully implemented on graphics processing units
• Exactly equivalent thermal conductivity in finite systems from equilibrium and nonequilibrium molecular dynamics simulations