1
|
Hicks CB, Martinez TJ. Massively scalable workflows for quantum chemistry: BigChem and ChemCloud. J Chem Phys 2024; 160:142501. [PMID: 38591672 DOI: 10.1063/5.0190834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 03/14/2024] [Indexed: 04/10/2024] Open
Abstract
Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. We present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet.
Collapse
Affiliation(s)
- Colton B Hicks
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J Martinez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
2
|
Ugandi M, Roemelt M. A configuration-based heatbath-CI for spin-adapted multireference electronic structure calculations with large active spaces. J Comput Chem 2023; 44:2374-2390. [PMID: 37589287 DOI: 10.1002/jcc.27203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Revised: 07/24/2023] [Accepted: 07/27/2023] [Indexed: 08/18/2023]
Abstract
This work reports on a spin-pure configuration-based implementation of the heatbath configuration interaction (HCI) algorithm for selective configuration interaction. Besides the obvious advantage of being spin-pure, the presented method combines the compactness of the configurational ansatz with the known efficiency of the HCI algorithm and a variety of algorithmic and conceptual ideas to achieve a high level of performance. In particular, through pruning of the selected configurational space after HCI selection by means of a more strict criterion, a more compact wavefunction representation is obtained. Moreover, the underlying logic of the method allows us to minimize the number of redundant matrix-matrix multiplications while making use of just-in-time compilation to achieve fast diagonalization of the Hamiltonian. The critical search for 2-electron connections within the configurational space is facilitated by a tree-based representation thereof as suggested previously by Gopal et al. Usage of a prefix-based parallelization and batching during the calculation of the PT2-correction leads to a good load balancing and significantly reduced memory requirements for these critical steps of the calculation. In this way, the need for a semistochastic approach to the PT2 correction is avoided even for large configurational spaces. Finally, several test-cases are discussed to demonstrate the strengths and weaknesses of the presented method.
Collapse
Affiliation(s)
- Mihkel Ugandi
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Michael Roemelt
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| |
Collapse
|
3
|
Williams-Young DB, Asadchev A, Popovici DT, Clark D, Waldrop J, Windus TL, Valeev EF, de Jong WA. Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory. J Chem Phys 2023; 158:234104. [PMID: 37326157 DOI: 10.1063/5.0151070] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2023] [Accepted: 05/26/2023] [Indexed: 06/17/2023] Open
Abstract
With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. In the present work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn-Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.
Collapse
Affiliation(s)
- David B Williams-Young
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Andrey Asadchev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Doru Thom Popovici
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - David Clark
- NVIDIA Corporation, Santa Clara, California 95051, USA
| | - Jonathan Waldrop
- Chemical and Biological Sciences Division, Ames National Laboratory, Ames, Iowa 50011, USA
| | - Theresa L Windus
- Chemical and Biological Sciences Division, Ames National Laboratory, Ames, Iowa 50011, USA
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Wibe A de Jong
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
4
|
Whitten JL. Correction of Residual Errors in Configuration Interaction Electronic Structure Calculations. J Chem Phys 2022; 157:124101. [DOI: 10.1063/5.0098793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Methods for correcting residual energy errors of configuration interaction (CI) calculations of molecules and other electronic systems are discussed based on the assumption that the energy defect can be mapped onto atomic regions. The methods do not consider the detailed nature of excitations, but instead define a defect energy per electron that that is unique to a specific atom. Defect energy contributions are determined from calculations on diatomic and hydride molecules and then applied to other systems. Calculated energies are compared with experimental thermodynamic and spectroscopic data for a set of forty-one mainly organic molecules representing a wide range of bonding environments. The most stringent test is based on a severely truncated virtual space in which higher spherical harmonic basis functions are removed. The errors of the initial CI calculations are large, but in each case, including defect corrections brings calculated CI energies into agreement with experimental values. The method is also applied to a NIST compilation of coupled-cluster calculations that employ a larger basis set and no truncation of the virtual space. The corrections show excellent consistency with total energies in very good agreement with experimental values. An extension of the method is applied to dmsn states of Sc, Ti, V, Mn, Cr, Fe, Co, Ni and Cu, significantly improving the agreement of calculated transition energies with spectroscopic values.
Collapse
Affiliation(s)
- Jerry L Whitten
- North Carolina State University, North Carolina State University, United States of America
| |
Collapse
|
5
|
Nandy A, Duan C, Kulik HJ. Audacity of huge: overcoming challenges of data scarcity and data quality for machine learning in computational materials discovery. Curr Opin Chem Eng 2022. [DOI: 10.1016/j.coche.2021.100778] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
6
|
Jeong W, Gaggioli CA, Gagliardi L. Active Learning Configuration Interaction for Excited-State Calculations of Polycyclic Aromatic Hydrocarbons. J Chem Theory Comput 2021; 17:7518-7530. [PMID: 34787422 PMCID: PMC8675132 DOI: 10.1021/acs.jctc.1c00769] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Indexed: 11/30/2022]
Abstract
We present the active learning configuration interaction (ALCI) method for multiconfigurational calculations based on large active spaces. ALCI leverages the use of an active learning procedure to find important electronic configurations among the full configurational space generated within an active space. We tested it for the calculation of singlet-singlet excited states of acenes and pyrene using different machine learning algorithms. The ALCI method yields excitation energies within 0.2-0.3 eV from those obtained by traditional complete active-space configuration interaction (CASCI) calculations (affordable for active spaces up to 16 electrons in 16 orbitals) by including only a small fraction of the CASCI configuration space in the calculations. For larger active spaces (we tested up to 26 electrons in 26 orbitals), not affordable with traditional CI methods, ALCI captures the trends of experimental excitation energies. Overall, ALCI provides satisfactory approximations to large active-space wave functions with up to 10 orders of magnitude fewer determinants for the systems presented here. These ALCI wave functions are promising and affordable starting points for the subsequent second-order perturbation theory or pair-density functional theory calculations.
Collapse
Affiliation(s)
- WooSeok Jeong
- Department
of Chemistry, Nanoporous Materials Genome Center, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Carlo Alberto Gaggioli
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Argonne
National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
7
|
Khedkar A, Roemelt M. Modern multireference methods and their application in transition metal chemistry. Phys Chem Chem Phys 2021; 23:17097-17112. [PMID: 34355719 DOI: 10.1039/d1cp02640b] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.
Collapse
Affiliation(s)
- Abhishek Khedkar
- Lehrstuhl für theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany.
| | | |
Collapse
|
8
|
Cruzeiro VWD, Manathunga M, Merz KM, Götz AW. Open-Source Multi-GPU-Accelerated QM/MM Simulations with AMBER and QUICK. J Chem Inf Model 2021; 61:2109-2115. [PMID: 33913331 DOI: 10.1021/acs.jcim.1c00169] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The quantum mechanics/molecular mechanics (QM/MM) approach is an essential and well-established tool in computational chemistry that has been widely applied in a myriad of biomolecular problems in the literature. In this publication, we report the integration of the QUantum Interaction Computational Kernel (QUICK) program as an engine to perform electronic structure calculations in QM/MM simulations with AMBER. This integration is available through either a file-based interface (FBI) or an application programming interface (API). Since QUICK is an open-source GPU-accelerated code with multi-GPU parallelization, users can take advantage of "free of charge" GPU-acceleration in their QM/MM simulations. In this work, we discuss implementation details and give usage examples. We also investigate energy conservation in typical QM/MM simulations performed at the microcanonical ensemble. Finally, benchmark results for two representative systems in bulk water, the N-methylacetamide (NMA) molecule and the photoactive yellow protein (PYP), show the performance of QM/MM simulations with QUICK and AMBER using a varying number of CPU cores and GPUs. Our results highlight the acceleration obtained from a single or multiple GPUs; we observed speedups of up to 53× between a single GPU vs a single CPU core and of up to 2.6× when comparing four GPUs to a single GPU. Results also reveal speedups of up to 3.5× when the API is used instead of FBI.
Collapse
Affiliation(s)
- Vinícius Wilian D Cruzeiro
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States.,Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Madushanka Manathunga
- Department of Chemistry, Department of Biochemistry and Molecular Biology, Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| | - Kenneth M Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
9
|
Levine BG, Durden AS, Esch MP, Liang F, Shu Y. CAS without SCF-Why to use CASCI and where to get the orbitals. J Chem Phys 2021; 154:090902. [PMID: 33685182 DOI: 10.1063/5.0042147] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.
Collapse
Affiliation(s)
- Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Michael P Esch
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
10
|
Whitten JL. Estimates of electron correlation based on density expansions. J Chem Phys 2020; 153:244103. [PMID: 33380085 DOI: 10.1063/5.0031279] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Two methods for estimating the correlation energy of molecules and other electronic systems are discussed based on the assumption that the correlation energy can be partitioned between atomic regions. In the first method, the electron density is expanded in terms of atomic contributions using rigorous electron repulsion bounds, and in the second method, correlation contributions are associated with basis function pairs. These methods do not consider the detailed nature of localized excitations but instead define a correlation energy per electron factor that is unique to a specific atom. The correlation factors are basis function dependent and are determined by configuration interaction (CI) calculations on diatomic and hydride molecules. The correlation energy estimates are compared with the results of high-level CI calculations for a test set of 27 molecules representing a wide range of bonding environments (average error of 2.6%). An extension based on truncated CI calculations in which d-type and hydrogen p-type functions are eliminated from the virtual space combined with estimates of dynamical correlation contributions using atomic correlation factors is discussed and applied to the dissociation of several molecules.
Collapse
Affiliation(s)
- Jerry L Whitten
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
11
|
Seritan S, Bannwarth C, Fales BS, Hohenstein EG, Isborn CM, Kokkila‐Schumacher SIL, Li X, Liu F, Luehr N, Snyder JW, Song C, Titov AV, Ufimtsev IS, Wang L, Martínez TJ. TeraChem
: A graphical processing unit
‐accelerated
electronic structure package for
large‐scale
ab initio molecular dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1494] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Stefan Seritan
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christoph Bannwarth
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Bryan S. Fales
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Edward G. Hohenstein
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christine M. Isborn
- Department of Chemistry University of California Merced Merced California USA
| | | | - Xin Li
- Division of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Fang Liu
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts USA
| | | | | | - Chenchen Song
- Department of Physics University of California Berkeley Berkeley California USA
- Molecular Foundry Lawrence Berkeley National Laboratory Berkeley California USA
| | | | - Ivan S. Ufimtsev
- Department of Structural Biology Stanford University School of Medicine Stanford California USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California Davis Davis California USA
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| |
Collapse
|
12
|
Seritan S, Bannwarth C, Fales BS, Hohenstein EG, Kokkila-Schumacher SIL, Luehr N, Snyder JW, Song C, Titov AV, Ufimtsev IS, Martínez TJ. TeraChem: Accelerating electronic structure and ab initio molecular dynamics with graphical processing units. J Chem Phys 2020; 152:224110. [PMID: 32534542 PMCID: PMC7928072 DOI: 10.1063/5.0007615] [Citation(s) in RCA: 76] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Accepted: 05/19/2020] [Indexed: 11/15/2022] Open
Abstract
Developed over the past decade, TeraChem is an electronic structure and ab initio molecular dynamics software package designed from the ground up to leverage graphics processing units (GPUs) to perform large-scale ground and excited state quantum chemistry calculations in the gas and the condensed phase. TeraChem's speed stems from the reformulation of conventional electronic structure theories in terms of a set of individually optimized high-performance electronic structure operations (e.g., Coulomb and exchange matrix builds, one- and two-particle density matrix builds) and rank-reduction techniques (e.g., tensor hypercontraction). Recent efforts have encapsulated these core operations and provided language-agnostic interfaces. This greatly increases the accessibility and flexibility of TeraChem as a platform to develop new electronic structure methods on GPUs and provides clear optimization targets for emerging parallel computing architectures.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | | | - Ivan S. Ufimtsev
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California 94305, USA
| | | |
Collapse
|
13
|
Manathunga M, Miao Y, Mu D, Götz AW, Merz KM. Parallel Implementation of Density Functional Theory Methods in the Quantum Interaction Computational Kernel Program. J Chem Theory Comput 2020; 16:4315-4326. [DOI: 10.1021/acs.jctc.0c00290] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Madushanka Manathunga
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Yipu Miao
- Facebook, 1 Hacker Way, Menlo Park, California 94025, United States
| | - Dawei Mu
- National Center for Supercomputing Applications, University of Illinois at Urbana−Champaign, 1205 W Clark St, Urbana, Illinois 61801, United States
| | - Andreas W. Götz
- San Diego Supercomputer Center, University of California, San Diego, 9500 Gilman Drive, La Jolla, California 92093-0505, United States
| | - Kenneth M. Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| |
Collapse
|
14
|
Fales BS, Martínez TJ. Fast transformations between configuration state function and Slater determinant bases for direct configuration interaction. J Chem Phys 2020; 152:164111. [PMID: 32357800 DOI: 10.1063/5.0005155] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A hybrid configuration state function (CSF) and Slater determinant (SD) basis full configuration interaction (CI) program was developed to simultaneously take advantage of fast SD basis algorithms for σ = Hc formation and the smaller CI vector length and more robust convergence offered by a CSF basis. Graphical processing unit acceleration of the direct CSF-SD and SD-CSF basis transformation algorithms ensures that the combined transformation time per iteration relative to σ formation is small (∼15%). In addition to the obvious benefits of reducing the memory footprint of the CI vector, additional computational savings are demonstrated that rely directly on the size of the CI basis, in one particular case reducing the CI time-to-solution of a HF-CAS-(16,16)-CI/6-31G calculation of ethylene from 1954.79 s to 956 s by using a CSF basis, a 2.0× speedup.
Collapse
Affiliation(s)
- B Scott Fales
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|