Dataset Information

Tinker-HP: Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields Using GPUs and Multi-GPU Systems.

ABSTRACT: We present the extension of the Tinker-HP package (Lagardère, Chem. Sci. 2018, 9, 956-972) to the use of Graphics Processing Unit (GPU) cards to accelerate molecular dynamics simulations using polarizable many-body force fields. The new high-performance module allows for an efficient use of single- and multiple-GPU architectures ranging from research laboratories to modern supercomputer centers. After detailing an analysis of our general scalable strategy that relies on OpenACC and CUDA, we discuss the various capabilities of the package. Among them, the multiprecision possibilities of the code are discussed. If an efficient double precision implementation is provided to preserve the possibility of fast reference computations, we show that a lower precision arithmetic is preferred providing a similar accuracy for molecular dynamics while exhibiting superior performances. As Tinker-HP is mainly dedicated to accelerate simulations using new generation point dipole polarizable force field, we focus our study on the implementation of the AMOEBA model. Testing various NVIDIA platforms including 2080Ti, 3090, V100, and A100 cards, we provide illustrative benchmarks of the code for single- and multicards simulations on large biosystems encompassing up to millions of atoms. The new code strongly reduces time to solution and offers the best performances to date obtained using the AMOEBA polarizable force field. Perspectives toward the strong-scaling performance of our multinode massive parallelization strategy, unsupervised adaptive sampling and large scale applicability of the Tinker-HP code in biophysics are discussed. The present software has been released in phase advance on GitHub in link with the High Performance Computing community COVID-19 research efforts and is free for Academics (see https://github.com/TinkerTools/tinker-hp).

SUBMITTER: Adjoua O

PROVIDER: S-EPMC8047816 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:We present Tinker-HP, a massively MPI parallel package dedicated to classical molecular dynamics (MD) and to multiscale simulations, using advanced polarizable force fields (PFF) encompassing distributed multipoles electrostatics. Tinker-HP is an evolution of the popular Tinker package code that conserves its simplicity of use and its reference double precision implementation for CPUs. Grounded on interdisciplinary efforts with applied mathematics, Tinker-HP allows for long polarizable MD simulations on large systems up to millions of atoms. We detail in the paper the newly developed extension of massively parallel 3D spatial decomposition to point dipole polarizable models as well as their coupling to efficient Krylov iterative and non-iterative polarization solvers. The design of the code allows the use of various computer systems ranging from laboratory workstations to modern petascale supercomputers with thousands of cores. Tinker-HP proposes therefore the first high-performance scalable CPU computing environment for the development of next generation point dipole PFFs and for production simulations. Strategies linking Tinker-HP to Quantum Mechanics (QM) in the framework of multiscale polarizable self-consistent QM/MD simulations are also provided. The possibilities, performances and scalability of the software are demonstrated via benchmarks calculations using the polarizable AMOEBA force field on systems ranging from large water boxes of increasing size and ionic liquids to (very) large biosystems encompassing several proteins as well as the complete satellite tobacco mosaic virus and ribosome structures. For small systems, Tinker-HP appears to be competitive with the Tinker-OpenMM GPU implementation of Tinker. As the system size grows, Tinker-HP remains operational thanks to its access to distributed memory and takes advantage of its new algorithmic enabling for stable long timescale polarizable simulations. Overall, a several thousand-fold acceleration over a single-core computation is observed for the largest systems. The extension of the present CPU implementation of Tinker-HP to other computational platforms is discussed.

Project description:A coupled polarization-matrix inversion and iteration (CPII) method is described to achieve and accelerate the convergence of induced dipoles for condensed phase systems employing polarizable intermolecular potential functions (PIPF). The present PIPF is based on the Thole interaction dipole model in which all atomic pair interactions are considered, including those that are directly bonded covalently. Although induced dipoles can be obtained both by inverting a 3N x 3N polarization-matrix where N is the number of polarizable sites, or by a direct iterative approach, the latter approach is more efficient computationally for large systems in molecular dynamics simulations. It was found that induced dipole moments failed to converge in the direct iterative approach if 1-2, 1-3, and 1-4 intramolecular interactions are included in the Thole model. However, it is necessary to include all intramolecular interactions in the Thole model to yield the correct molecular anisotropic polarizability tensor. To solve this numerical stability problem, we reformulated the Thole interaction dipole model in terms of molecular block matrices, which naturally leads to a coupled, preconditioning algorithm that involves a polarization-matrix inversion term to account for intramolecular interactions, and an iterative procedure to incorporate the mutual polarization effects between different molecules. The CPII method is illustrated by applying to cubic boxes of water and NMA molecules as well as an alanine pentapeptide configuration, and it was shown that the CPII method can achieve convergence for the dipole induction polarization rapidly in all cases, whereas the direct iterative approach failed to reach convergence in these cases. In addition, the CPII reduces the overall computational costs by decreasing the number of iteration steps in comparison with the direct iteration approach in which intramolecular bonded interactions are excluded to ensure that induced dipole convergence is obtained.

Dataset Information

Tinker-HP: Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields Using GPUs and Multi-GPU Systems.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets