Abstract:The conductor-like polarization model (C-PCM) with switching/Gaussian smooth discretization is a widely used implicit solvation model in chemical simulations. However, its application in quantum mechanical calculations of large-scale biomolecular systems can be limited by computational expense of both the gas phase electronic structure and the solvation interaction. We have previously used graphical processing units (GPUs) to accelerate the first of these steps. Here, we extend the use of GPUs to accelerate electronic structure calculations including C-PCM solvation. Implementation on the GPU leads to significant acceleration of the generation of the required integrals for C-PCM. We further propose two strategies to improve the solution of the required linear equations: a dynamic convergence threshold and a randomized block-Jacobi preconditioner. These strategies are not specific to GPUs and are expected to be beneficial for both CPU and GPU implementations. We benchmark the performance of the new implementation using over 20 small proteins in solvent environment. Using a single GPU, our method evaluates the C-PCM related integrals and their derivatives more than 10X faster than a conventional CPU based implementation. Our improvements to the linear solver provide a further 3X acceleration. The overall calculations including C-PCM solvation require typically 20-40% more effort than their gas phase counterparts for moderate basis set and molecule surface discretization level. The relative cost of the C-PCM solvation correction decreases as the basis sets and/or cavity radii increase. Therefore description of solvation with this model should be routine. We also discuss applications to the study of the conformational landscape of an amyloid fibril.

Introducing GPU-acceleration into the Python-based Simulations of Chemistry Framework

Enhancing GPU-acceleration in the Python-based Simulations of Chemistry Framework

The Implementation of the Three-Dimensional Unified Gas-Kinetic Wave-Particle Method on Multiple Graphics Processing Units

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

GPU-Accelerated On-the-fly Nonadiabatic Semiclassical Dynamics

Accelerating Relativistic Exact-Two-Component Density Functional Theory Calculations with Graphical Processing Units

The Python-based Simulations of Chemistry Framework (PySCF)

Fast and Scalable GPU-Accelerated Quantum Chemistry for Periodic Systems with Gaussian Orbitals: Implementation and Hybrid Density Functional Theory Calculations

Extending GPU-Accelerated Gaussian Integrals in the TeraChem Software Package to f Type Orbitals: Implementation and Applications

Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units

Multi-GPU RI-HF Energies and Analytic Gradients $-$ Towards High Throughput Ab Initio Molecular Dynamics

Quantum Chemistry for Solvated Molecules on Graphical Processing Units (GPUs)using Polarizable Continuum Models

Accurate QM/MM Molecular Dynamics for Periodic Systems in \textsc{GPU4PySCF} with Applications to Enzyme Catalysis

GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions

GPU-accelerated Effective Hamiltonian Calculator

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

Robust and effective ab initio molecular dynamics simulations on the GPU cloud infrastructure using the Schrödinger Materials Science Suite

GPU Acceleration of Numerical Atomic Orbitals-Based Density Functional Theory Algorithms within the ABACUS package

OpenCL-accelerated first-principles calculations of all-electron quantum perturbations on HPC resources

High-performance GPU-accelerated evaluation of electron repulsion integrals

Accelerating Pythonic coupled cluster implementations: a comparison between CPUs and GPUs