Molecular modeling on GPUs

Last updated
Ionic liquid simulation on GPU (Abalone) Hardware-accelerated-molecular-modeling.png
Ionic liquid simulation on GPU (Abalone)

Molecular modeling on GPU is the technique of using a graphics processing unit (GPU) for molecular simulations. [1]

Contents

In 2007, Nvidia introduced video cards that could be used not only to show graphics but also for scientific calculations. These cards include many arithmetic units (as of 2016, up to 3,584 in Tesla P100) working in parallel. Long before this event, the computational power of video cards was purely used to accelerate graphics calculations. What was new is that Nvidia made it possible to develop parallel programs in a high-level application programming interface (API) named CUDA. This technology substantially simplified programming by enabling programs to be written in C/C++. More recently, OpenCL allows cross-platform GPU acceleration.

Quantum chemistry calculations [2] [3] [4] [5] [6] [7] and molecular mechanics simulations [8] [9] [10] (molecular modeling in terms of classical mechanics) are among beneficial applications of this technology. The video cards can accelerate the calculations tens of times, so a PC with such a card has the power similar to that of a cluster of workstations based on common processors.

GPU accelerated molecular modelling software

Programs

API

Distributed computing projects

See also

Related Research Articles

<span class="mw-page-title-main">Computational chemistry</span> Branch of chemistry

Computational chemistry is a branch of chemistry that uses computer simulations to assist in solving chemical problems. It uses methods of theoretical chemistry incorporated into computer programs to calculate the structures and properties of molecules, groups of molecules, and solids. The importance of this subject stems from the fact that, with the exception of some relatively recent findings related to the hydrogen molecular ion, achieving an accurate quantum mechanical depiction of chemical systems analytically, or in a closed form, is not feasible. The complexity inherent in the many-body problem exacerbates the challenge of providing detailed descriptions of quantum mechanical systems. While computational results normally complement information obtained by chemical experiments, it can occasionally predict unobserved chemical phenomena.

<span class="mw-page-title-main">Molecular dynamics</span> Computer simulations to discover and understand chemical properties

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields. The method is applied mostly in chemical physics, materials science, and biophysics.

GROMACS is a molecular dynamics package mainly designed for simulations of proteins, lipids, and nucleic acids. It was originally developed in the Biophysical Chemistry department of University of Groningen, and is now maintained by contributors in universities and research centers worldwide. GROMACS is one of the fastest and most popular software packages available, and can run on central processing units (CPUs) and graphics processing units (GPUs). It is free, open-source software released under the GNU Lesser General Public License (LGPL).

<span class="mw-page-title-main">Molecular mechanics</span> Use of classical mechanics to model molecular systems

Molecular mechanics uses classical mechanics to model molecular systems. The Born–Oppenheimer approximation is assumed valid and the potential energy of all systems is calculated as a function of the nuclear coordinates using force fields. Molecular mechanics can be used to study molecule systems ranging in size and complexity from small to large biological systems or material assemblies with many thousands to millions of atoms.

<span class="mw-page-title-main">Visual Molecular Dynamics</span> Visualization and modelling software

Visual Molecular Dynamics (VMD) is a molecular modelling and visualization computer program. VMD is developed mainly as a tool to view and analyze the results of molecular dynamics simulations. It also includes tools for working with volumetric data, sequence data, and arbitrary graphics objects. Molecular scenes can be exported to external rendering tools such as POV-Ray, RenderMan, Tachyon, Virtual Reality Modeling Language (VRML), and many others. Users can run their own Tcl and Python scripts within VMD as it includes embedded Tcl and Python interpreters. VMD runs on Unix, Apple Mac macOS, and Microsoft Windows. VMD is available to non-commercial users under a distribution-specific license which permits both use of the program and modification of its source code, at no charge.

<span class="mw-page-title-main">Molecular modelling</span> Discovering chemical properties by physical simulations

Molecular modelling encompasses all methods, theoretical and computational, used to model or mimic the behaviour of molecules. The methods are used in the fields of computational chemistry, drug design, computational biology and materials science to study molecular systems ranging from small chemical systems to large biological molecules and material assemblies. The simplest calculations can be performed by hand, but inevitably computers are required to perform molecular modelling of any reasonably sized system. The common feature of molecular modelling methods is the atomistic level description of the molecular systems. This may include treating atoms as the smallest individual unit, or explicitly modelling protons and neutrons with its quarks, anti-quarks and gluons and electrons with its photons.

Q-Chem is a general-purpose electronic structure package featuring a variety of established and new methods implemented using innovative algorithms that enable fast calculations of large systems on various computer architectures, from laptops and regular lab workstations to midsize clusters, HPCC, and cloud computing using density functional and wave-function based approaches. It offers an integrated graphical interface and input generator; a large selection of functionals and correlation methods, including methods for electronically excited states and open-shell systems; solvation models; and wave-function analysis tools. In addition to serving the computational chemistry community, Q-Chem also provides a versatile code development platform.

<span class="mw-page-title-main">PQS (software)</span> Quantum chemistry software program

PQS is a general purpose quantum chemistry program. Its roots go back to the first ab initio gradient program developed in Professor Peter Pulay's group but now it is developed and distributed commercially by Parallel Quantum Solutions. There is a reduction in cost for academic users and a site license. Its strong points are geometry optimization, NMR chemical shift calculations, and large MP2 calculations, and high parallel efficiency on computing clusters. It includes many other capabilities including Density functional theory, the semiempirical methods, MINDO/3, MNDO, AM1 and PM3, Molecular mechanics using the SYBYL 5.0 Force Field, the quantum mechanics/molecular mechanics mixed method using the ONIOM method, natural bond orbital (NBO) analysis and COSMO solvation models. Recently, a highly efficient parallel CCSD(T) code for closed shell systems has been developed. This code includes many other post Hartree–Fock methods: MP2, MP3, MP4, CISD, CEPA, QCISD and so on.

<span class="mw-page-title-main">MOLCAS</span> Computational chemistry software

MOLCAS is an ab initio computational chemistry program, developed as a joint project by a number of international institutes. MOLCAS is developed by scientists to be used by scientists. It is not primarily a commercial product and it is not sold in order to produce a fortune for its owner.

Tinker, previously stylized as TINKER, is a suite of computer software applications for molecular dynamics simulation. The codes provide a complete and general set of tools for molecular mechanics and molecular dynamics, with some special features for biomolecules. The core of the software is a modular set of callable routines which allow manipulating coordinates and evaluating potential energy and derivatives via straightforward means.

Car–Parrinello molecular dynamics or CPMD refers to either a method used in molecular dynamics or the computational chemistry software package used to implement this method.

This is a list of computer programs that are predominantly used for molecular mechanics calculations.

<span class="mw-page-title-main">CP2K</span>

CP2K is a freely available (GPL) quantum chemistry and solid state physics program package, written in Fortran 2008, to perform atomistic simulations of solid state, liquid, molecular, periodic, material, crystal, and biological systems. It provides a general framework for different methods: density functional theory (DFT) using a mixed Gaussian and plane waves approach (GPW) via LDA, GGA, MP2, or RPA levels of theory, classical pair and many-body potentials, semi-empirical and tight-binding Hamiltonians, as well as Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid schemes relying on the Gaussian Expansion of the Electrostatic Potential (GEEP). The Gaussian and Augmented Plane Waves method (GAPW) as an extension of the GPW method allows for all-electron calculations. CP2K can do simulations of molecular dynamics, metadynamics, Monte Carlo, Ehrenfest dynamics, vibrational analysis, core level spectroscopy, energy minimization, and transition state optimization using NEB or dimer method.

<span class="mw-page-title-main">Ascalaph Designer</span>

Ascalaph Designer is a computer program for general purpose molecular modelling for molecular design and simulations. It provides a graphical environment for the common programs of quantum and classical molecular modelling ORCA, NWChem, Firefly, CP2K and MDynaMix . The molecular mechanics calculations cover model building, energy optimizations and molecular dynamics. Firefly covers a wide range of quantum chemistry methods. Ascalaph Designer is free and open-source software, released under the GNU General Public License, version 2 (GPLv2).

TeraChem is a computational chemistry software program designed for CUDA-enabled Nvidia GPUs. The initial development started at the University of Illinois at Urbana-Champaign and was subsequently commercialized. It is currently distributed by PetaChem, LLC, located in Silicon Valley. As of 2020, the software package is still under active development.

BigDFT is a free software package for physicists and chemists, distributed under the GNU General Public License, whose main program allows the total energy, charge density, and electronic structure of systems made of electrons and nuclei to be calculated within density functional theory (DFT), using pseudopotentials, and a wavelet basis.

Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) is a Coarse-grained modeling package, which focuses on the analysis of molecular dynamics data, the development of systematic coarse-graining techniques as well as methods used for simulating microscopic charge transport in disordered semiconductors. It was originally developed at the Max Planck Institute for Polymer Research, and is now maintained by developers at the Max Planck Institute for Polymer Research, Los Alamos National Laboratory, Eindhoven University of Technology and the Beckman Institute for Advanced Science and Technology with contributions from researcher worldwide.

In the context of chemistry and molecular modelling, the Interface force field (IFF) is a force field for classical molecular simulations of atoms, molecules, and assemblies up to the large nanometer scale, covering compounds from across the periodic table. It employs a consistent classical Hamiltonian energy function for metals, oxides, and organic compounds, linking biomolecular and materials simulation platforms into a single platform. The reliability is often higher than that of density functional theory calculations at more than a million times lower computational cost. IFF includes a physical-chemical interpretation for all parameters as well as a surface model database that covers different cleavage planes and surface chemistry of included compounds. The Interface Force Field is compatible with force fields for the simulation of primarily organic compounds and can be used with common molecular dynamics and Monte Carlo codes. Structures and energies of included chemical elements and compounds are rigorously validated and property predictions are up to a factor of 100 more accurate relative to earlier models.

References

  1. Stone JE, Phillips JC, Freddolino PL, Hardy DJ, Trabuco LG, Schulten K (December 2007). "Accelerating molecular modeling applications with graphics processors". Journal of Computational Chemistry. 28 (16): 2618–2640. CiteSeerX   10.1.1.466.3823 . doi:10.1002/jcc.20829. PMID   17894371. S2CID   15313533.
  2. Yasuda K (August 2008). "Accelerating Density Functional Calculations with Graphics Processing Unit". Journal of Chemical Theory and Computation. 4 (8): 1230–1236. doi:10.1021/ct8001046. PMID   26631699.
  3. Yasuda K (February 2008). "Two-electron integral evaluation on the graphics processor unit". Journal of Computational Chemistry. 29 (3): 334–342. CiteSeerX   10.1.1.498.364 . doi:10.1002/jcc.20779. PMID   17614340. S2CID   8078401.
  4. Vogt L, Olivares-Amaya R, Kermes S, Shao Y, Amador-Bedolla C, Aspuru-Guzik A (March 2008). "Accelerating resolution-of-the-identity second-order Møller-Plesset quantum chemistry calculations with graphical processing units". The Journal of Physical Chemistry A. 112 (10): 2049–2057. Bibcode:2008JPCA..112.2049V. doi:10.1021/jp0776762. PMID   18229900. S2CID   4566211.
  5. Ufimtsev IS, Martínez TJ (February 2008). "Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation". Journal of Chemical Theory and Computation. 4 (2): 222–231. doi:10.1021/ct700268q. PMID   26620654.
  6. Ivan S. Ufimtsev & Todd J. Martinez (2008). "Graphical Processing Units for Quantum Chemistry". Computing in Science & Engineering. 10 (6): 26–34. Bibcode:2008CSE....10f..26U. doi:10.1109/MCSE.2008.148. S2CID   10225262.
  7. Tornai GJ, Ladjánszki I, Rák Á, Kis G, Cserey G (October 2019). "Calculation of Quantum Chemical Two-Electron Integrals by Applying Compiler Technology on GPU". Journal of Chemical Theory and Computation. 15 (10): 5319–5331. doi:10.1021/acs.jctc.9b00560. PMID   31503475. S2CID   202555796.
  8. Joshua A. Anderson; Chris D. Lorenz; A. Travesset (2008). "General Purpose Molecular Dynamics Simulations Fully Implemented on Graphics Processing Units". Journal of Computational Physics. 227 (10): 5342–5359. Bibcode:2008JCoPh.227.5342A. CiteSeerX   10.1.1.552.2883 . doi:10.1016/j.jcp.2008.01.047.
  9. Christopher I. Rodrigues; David J. Hardy; John E. Stone; Klaus Schulten & Wen-Mei W. Hwu. (2008). "GPU acceleration of cutoff pair potentials for molecular modeling applications". In CF'08: Proceedings of the 2008 Conference on Computing Frontiers, New York, NY, USA: 273–282.
  10. Peter H. Colberg; Felix Höfling (2011). "Highly accelerated simulations of glassy dynamics using GPUs: Caveats on limited floating-point precision". Comput. Phys. Commun. 182 (5): 1120–1129. arXiv: 0912.3824 . Bibcode:2011CoPhC.182.1120C. doi:10.1016/j.cpc.2011.01.009. S2CID   7173093.
  11. Yousif RH (2020). "Exploring the Molecular Interactions between Neoculin and the Human Sweet Taste Receptors through Computational Approaches" (PDF). Sains Malaysiana. 49 (3): 517–525. doi: 10.17576/jsm-2020-4903-06 .
  12. Bailey N, Ingebrigtsen T, Hansen JS, Veldhorst A, Bøhling L, Lemarchand C, et al. (2017-12-14). "RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles". SciPost Physics. 3 (6): 038. arXiv: 1506.05094 . Bibcode:2017ScPP....3...38B. doi: 10.21468/SciPostPhys.3.6.038 . ISSN   2542-4653. S2CID   43964588.
  13. Harger M, Li D, Wang Z, Dalby K, Lagardère L, Piquemal JP, et al. (September 2017). "Tinker-OpenMM: Absolute and relative alchemical free energies using AMOEBA on GPUs". Journal of Computational Chemistry. 38 (23): 2047–2055. doi:10.1002/jcc.24853. PMC   5539969 . PMID   28600826.