Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Chemical Theory and Computation

Quantum Chemistry on Graphical Processing Units. 2. Direct Self-Consistent-Field Implementation

Ivan S. Ufimtsev, Todd J. Martı́nez

Published Mar 6, 2009• 451 citations
Publisher

Abstract

We demonstrate the use of graphical processing units (GPUs) to carry out complete self-consistent-field calculations for molecules with as many as 453 atoms (2131 basis functions). Speedups ranging from 28× to 650× are achieved as compared to a mature third-party quantum chemistry program (GAMESS) running on a traditional CPU. The computational organization used to construct the Coulomb and exchange operators is discussed. We also present results using three GPUs in parallel, combining coarse and fine-grained parallelism.

DOI
10.1021/ct800526s

Related papers

Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation

Ivan S. Ufimtsev, Todd J. Martı́nez

Modern videogames place increasing demands on the computational and graphical hardware, leading to novel architectures that have great potential in the context of high performance computing and molecular simulation. We demonstrate that Graphical Processing Units (GPUs) can be used very efficiently to calculate two-electron repulsion integrals over Gaussian basis functions [Formula: see text] the first step in most quantum chemistry calculations. A benchmark test performed for the evaluation of approximately 10(6) (ss|ss) integrals over contracted s-orbitals showed that a naïve algorithm implemented on the GPU achieves up to 130-fold speedup over a traditional CPU implementation on an AMD Opteron. Subsequent calculations of the Coulomb operator for a 256-atom DNA strand show that the GPU advantage is maintained for basis sets including higher angular momentum functions.

Journal of Chemical Theory and Computation•Jan 25, 2008•606 citations

Many-body quantum chemistry on graphics processing units

A. Eugene DePrince, Jeff R. Hammond, Stephen K. Gray

Heterogeneous nodes composed of a multicore CPU and at least one graphics processing unit (GPU) are increasingly common in high-performance scientific computing, and significant programming effort is currently being undertaken to port existing scientific algorithms to these unique architectures. We present implementations for two many-body quantum chemistry methods on heterogeneous nodes: the coupled-cluster with single and double excitations (CCSD) and time-dependent configuration interaction with single and double excitations (TD-CISD) methods. Both methods can be implemented on a computer as a series of dense matrix-matrix multiplications, operations that GPUs are particularly adept at performing. The GPU-accelerated CCSD algorithm is as much as 4.3 times faster than the corresponding CPU algorithm and 9.7 times faster than the algorithm in the Molpro package. The TD-CISD algorithm is accelerated by as much as a factor of 3.9 by GPUs. Enhanced performance is achieved by overlapping CPU and GPU computations.

•Jan 1, 2011•1 citation

Extending GPU-accelerated Gaussian integrals in the TeraChem software package to f type orbitals: Implementation and applications

Yuanheng Wang, Diptarka Hait, Katherine Grace Johnson +4 more

The increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH3) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.

arXiv:2406.14920•The Journal of Chemical Physics•Nov 6, 2024•19 citations

Parallel Implementation of the Electronic Structure Code GAMESS

Theresa L. Windus, Michael W. I. Schmidt, Mark S. Gordon

This paper outlines various tools and techniques for the parallelization of quantum chemistry codes; in particular, for the electronic structure code GAMESS. A general overview of the parallel capabilities of GAMESS are also presented. The parallelization of quantum chemistry codes has become a very active area of research over the last decade( 1, 2, 3, 4 ). Until recently, most of this research has dealt with self-consistent field (SCF) theory( 1 ). However, in the last few years parallel implementations of post-SCF methods have been presented ( 2 ). Most of the post-SCF methods and analytic Hessians for SCF wavefunctions face the substantial problem of parallelizing the atomic orbital (AO) integral to molecular orbital (MO) integral transformation ( 3 ). The objective of this paper is to provide general information about the parallel implementation of GAMESS. The following sections are presented in this paper: (A) a brief overview of the functionality of the ab initio code GAMESS (General Atomic and Molecular Electronic Structure System); (B) a short discussion of the model, software, and general ideas used to parallelize GAMESS; (C) spécifics concerning the parallelization of the SCF; (D) discussion concerning the AO to MO integral transformation; (E) the transformation as applied to multi-configuration SCF (MCSCF); (F) the transformation as applied to analytic Hessians; (G) a brief overview of the parallel MP2 code; and (H) conclusions and future areas of research will be discussed.

ACS symposium series•May 17, 1995•7 citations

GPU-Accelerated Graph-Based Semiempirical Quantum Chemistry

Maksim Kulichenko, Robert M. Stanton, Cheng-Han Li +6 more

Graph-based electronic structure theory offers a scalable approach to study large, complex atomistic systems using distributed and hybrid computational platforms. We demonstrate the coupling of graph-based linear scaling electronic structure theory, as implemented in the Scalable Ecosystem, Driver, and Analyzer for Complex Chemistry Simulations (SEDACS), with semiempirical quantum chemistry methods as implemented in the PySEQM code, with Graphics Processing Unit (GPU) acceleration. This powerful combination enables efficient, scalable electronic structure calculations over many nodes, significantly reducing computational cost while naturally harnessing parallelism. Detailed analyses of parallelization efficiency, computational accuracy, and communication overheads are provided, highlighting an order-of-magnitude speedup for systems of up to 10,000 atoms.

Journal of Chemical Theory and Computation•Oct 27, 2025•1 citation

Quantum supercharger library: Hyper-parallelism of the Hartree-Fock method

Kyle D. Fernandes, C. Alicia Renison, Kevin J. Naidoo

We present here a set of algorithms that completely rewrites the Hartree-Fock (HF) computations common to many legacy electronic structure packages (such as GAMESS-US, GAMESS-UK, and NWChem) into a massively parallel compute scheme that takes advantage of hardware accelerators such as Graphical Processing Units (GPUs). The HF compute algorithm is core to a library of routines that we name the Quantum Supercharger Library (QSL). We briefly evaluate the QSL's performance and report that it accelerates a HF 6-31G Self-Consistent Field (SCF) computation by up to 20 times for medium sized molecules (such as a buckyball) when compared with mature Central Processing Unit algorithms available in the legacy codes in regular use by researchers. It achieves this acceleration by massive parallelization of the one- and two-electron integrals and optimization of the SCF and Direct Inversion in the Iterative Subspace routines through the use of GPU linear algebra libraries. © 2015 Wiley Periodicals, Inc.

Journal of Computational Chemistry•May 14, 2015•10 citations
We demonstrate the use of graphical processing units (GPUs) to carry out complete self-consistent-field calculations for molecules with as many as 453 atoms (2131 basis functions). Speedups ranging from 28× to 650× are achieved as compared to a mature third-party quantum chemistry program (GAMESS) running on a traditional CPU. The computational organization used to construct the Coulomb and exchange operators is discussed. We also present results using three GPUs in parallel, combining coarse and fine-grained parallelism.
Modern videogames place increasing demands on the computational and graphical hardware, leading to novel architectures that have great potential in the context of high performance computing and molecular simulation. We demonstrate that Graphical Processing Units (GPUs) can be used very efficiently to calculate two-electron repulsion integrals over Gaussian basis functions [Formula: see text] the first step in most quantum chemistry calculations. A benchmark test performed for the evaluation of approximately 10(6) (ss|ss) integrals over contracted s-orbitals showed that a naïve algorithm implemented on the GPU achieves up to 130-fold speedup over a traditional CPU implementation on an AMD Opteron. Subsequent calculations of the Coulomb operator for a 256-atom DNA strand show that the GPU advantage is maintained for basis sets including higher angular momentum functions.
Heterogeneous nodes composed of a multicore CPU and at least one graphics processing unit (GPU) are increasingly common in high-performance scientific computing, and significant programming effort is currently being undertaken to port existing scientific algorithms to these unique architectures. We present implementations for two many-body quantum chemistry methods on heterogeneous nodes: the coupled-cluster with single and double excitations (CCSD) and time-dependent configuration interaction with single and double excitations (TD-CISD) methods. Both methods can be implemented on a computer as a series of dense matrix-matrix multiplications, operations that GPUs are particularly adept at performing. The GPU-accelerated CCSD algorithm is as much as 4.3 times faster than the corresponding CPU algorithm and 9.7 times faster than the algorithm in the Molpro package. The TD-CISD algorithm is accelerated by as much as a factor of 3.9 by GPUs. Enhanced performance is achieved by overlapping CPU and GPU computations.
The increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH3) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.
This paper outlines various tools and techniques for the parallelization of quantum chemistry codes; in particular, for the electronic structure code GAMESS. A general overview of the parallel capabilities of GAMESS are also presented. The parallelization of quantum chemistry codes has become a very active area of research over the last decade( 1, 2, 3, 4 ). Until recently, most of this research has dealt with self-consistent field (SCF) theory( 1 ). However, in the last few years parallel implementations of post-SCF methods have been presented ( 2 ). Most of the post-SCF methods and analytic Hessians for SCF wavefunctions face the substantial problem of parallelizing the atomic orbital (AO) integral to molecular orbital (MO) integral transformation ( 3 ). The objective of this paper is to provide general information about the parallel implementation of GAMESS. The following sections are presented in this paper: (A) a brief overview of the functionality of the ab initio code GAMESS (General Atomic and Molecular Electronic Structure System); (B) a short discussion of the model, software, and general ideas used to parallelize GAMESS; (C) spécifics concerning the parallelization of the SCF; (D) discussion concerning the AO to MO integral transformation; (E) the transformation as applied to multi-configuration SCF (MCSCF); (F) the transformation as applied to analytic Hessians; (G) a brief overview of the parallel MP2 code; and (H) conclusions and future areas of research will be discussed.
Graph-based electronic structure theory offers a scalable approach to study large, complex atomistic systems using distributed and hybrid computational platforms. We demonstrate the coupling of graph-based linear scaling electronic structure theory, as implemented in the Scalable Ecosystem, Driver, and Analyzer for Complex Chemistry Simulations (SEDACS), with semiempirical quantum chemistry methods as implemented in the PySEQM code, with Graphics Processing Unit (GPU) acceleration. This powerful combination enables efficient, scalable electronic structure calculations over many nodes, significantly reducing computational cost while naturally harnessing parallelism. Detailed analyses of parallelization efficiency, computational accuracy, and communication overheads are provided, highlighting an order-of-magnitude speedup for systems of up to 10,000 atoms.
We present here a set of algorithms that completely rewrites the Hartree-Fock (HF) computations common to many legacy electronic structure packages (such as GAMESS-US, GAMESS-UK, and NWChem) into a massively parallel compute scheme that takes advantage of hardware accelerators such as Graphical Processing Units (GPUs). The HF compute algorithm is core to a library of routines that we name the Quantum Supercharger Library (QSL). We briefly evaluate the QSL's performance and report that it accelerates a HF 6-31G Self-Consistent Field (SCF) computation by up to 20 times for medium sized molecules (such as a buckyball) when compared with mature Central Processing Unit algorithms available in the legacy codes in regular use by researchers. It achieves this acceleration by massive parallelization of the one- and two-electron integrals and optimization of the SCF and Direct Inversion in the Iterative Subspace routines through the use of GPU linear algebra libraries. © 2015 Wiley Periodicals, Inc.