Abstract:We
develop a new implementation of coupled-cluster singles and
doubles (CCSD) optimized for the most recent graphical processing
unit (GPU) hardware. We find that a single node with 8 NVIDIA V100
GPUs is capable of performing CCSD computations on roughly 100 atoms
and 1300 basis functions in less than 1 day. Comparisons against massively
parallel implementations of CCSD suggest that more than 64 CPU-based
nodes (each with 16 cores) are required to match this performance.
“…The SOS‐MP2 method 86 neglects exchange‐like terms to arrive at a formal scaling of O (N 4 ), and THC‐SOS‐MP2 further reduces this scaling to O (N 3 ). Recently, a GPU‐accelerated coupled‐cluster code has also been implemented in TeraChem, enabling coupled‐cluster singles and doubles (CCSD) as well as any method that can be written as a subset of CCSD diagrams 87 …”
Section: Single Reference Post‐scf Methodsmentioning
TeraChem was born in 2008 with the goal of providing fast on‐the‐fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low‐cost parallel computer architecture that became more accessible for general‐purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large‐scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large‐scale nonadiabatic dynamics simulations.
This article is categorized under:
Electronic Structure Theory > Ab Initio Electronic Structure Methods
Software > Quantum Chemistry
Structure and Mechanism > Computational Biochemistry and Biophysics
“…The SOS‐MP2 method 86 neglects exchange‐like terms to arrive at a formal scaling of O (N 4 ), and THC‐SOS‐MP2 further reduces this scaling to O (N 3 ). Recently, a GPU‐accelerated coupled‐cluster code has also been implemented in TeraChem, enabling coupled‐cluster singles and doubles (CCSD) as well as any method that can be written as a subset of CCSD diagrams 87 …”
Section: Single Reference Post‐scf Methodsmentioning
TeraChem was born in 2008 with the goal of providing fast on‐the‐fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low‐cost parallel computer architecture that became more accessible for general‐purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large‐scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large‐scale nonadiabatic dynamics simulations.
This article is categorized under:
Electronic Structure Theory > Ab Initio Electronic Structure Methods
Software > Quantum Chemistry
Structure and Mechanism > Computational Biochemistry and Biophysics
“…Dr. Hohenstein then discussed an implementation of CCSD and EOM-CCSD algorithms on graphics cards, which makes the calculation of 100 atom systems computationally feasible. The CCSD calculations on eight graphics cards matched the performance of sixty-four 16-core CPU-based nodes …”
“…14,15 They have, for instance, demonstrated that NCIs are strongly attenuated in a solvated environment. 16 Despite the intense efforts associated with developing a posteriori dispersion corrections [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34] or with reducing the cost of wavefunction-based methods, [35][36][37][38][39][40][41] the finite-temperature description of competing NCIs remains coarse if the solvent is represented as a continuum. In particular, interaction energies estimated by dispersion-corrected computations in implicitly modeled solvent can be one magnitude larger than experimental energies in solutions.…”
Non-covalent bonding patterns are commonly harvested as a design principle in the field of catalysis, supramolecular chemistry and functional materials to name a few. Yet, their computational description generally neglects finite temperature and environment effects, which promote competing interactions and alter their static gas-phase properties. Recently, neural network potentials (NNPs) trained on Density Functional Theory (DFT) data have become increasingly popular to simulate molecular phenomena in condensed phase with an accuracy comparable to ab initio methods. To date, most applications have centered on solid-state materials or fairly simple molecules made of a limited number of elements. Herein, we focus on the persistence and strength of chalcogen bonds involving a benzotelluradiazole in condensed phase. While the tellurium-containing heteroaromatic molecules are known to exhibit pronounced interactions with anions and lone pairs of different atoms, the relevance of competing intermolecular interactions, notably with the solvent, are complicated to monitor experimentally but also challenging to model at an accurate electronic structure level. Here, we train a direct and baselined NNPs to reproduce hybrid DFT energies and forces in order to identify what are the most prevalent non-covalent interactions occurring in a solute-Cl --THF mixture. The simulations in explicit solvent highlight the clear competition with chalcogen bonds formed with the solvent and the short-range directionality of the interaction with direct consequences for the molecular properties in the solution. The comparison with other potentials (e.g., AMOEBA, direct NNP and continuum solvent model) also demonstrates that baselined NNPs offer a reliable picture of the non-covalent interaction interplay occurring in solution.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.