Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method
Abstract
We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.
- Authors:
-
- Univ. of California, Berkeley, CA (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Univ. of California, Berkeley, CA (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Kavli Energy NanoScience Inst. of Berkeley, Berkeley, CA (United States)
- SLAC National Accelerator Lab., Menlo Park, CA (United States). Photon Ultrafast Laser Science and Engineering Inst. (PULSE); Stanford Univ., CA (United States)
- Publication Date:
- Research Org.:
- SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States); University of California, Berkeley, CA (United States); Stanford Univ., CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Basic Energy Sciences (BES); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC)
- OSTI Identifier:
- 1767927
- Alternate Identifier(s):
- OSTI ID: 1755455
- Grant/Contract Number:
- AC02-76SF00515; AC02-05CH11231; SciDAC
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Chemical Physics
- Additional Journal Information:
- Journal Volume: 154; Journal Issue: 1; Journal ID: ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
Citation Formats
Song, Chenchen, Neaton, Jeffrey B., and Martínez, Todd J. Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method. United States: N. p., 2021.
Web. doi:10.1063/5.0035233.
Song, Chenchen, Neaton, Jeffrey B., & Martínez, Todd J. Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method. United States. https://doi.org/10.1063/5.0035233
Song, Chenchen, Neaton, Jeffrey B., and Martínez, Todd J. Thu .
"Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method". United States. https://doi.org/10.1063/5.0035233. https://www.osti.gov/servlets/purl/1767927.
@article{osti_1767927,
title = {Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method},
author = {Song, Chenchen and Neaton, Jeffrey B. and Martínez, Todd J.},
abstractNote = {We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.},
doi = {10.1063/5.0035233},
journal = {Journal of Chemical Physics},
number = 1,
volume = 154,
place = {United States},
year = {Thu Jan 07 00:00:00 EST 2021},
month = {Thu Jan 07 00:00:00 EST 2021}
}
Works referenced in this record:
Analytic first derivatives of floating occupation molecular orbital-complete active space configuration interaction on graphical processing units
journal, July 2015
- Hohenstein, Edward G.; Bouduban, Marine E. F.; Song, Chenchen
- The Journal of Chemical Physics, Vol. 143, Issue 1
Classification of configurations and the determination of interacting and noninteracting spaces in configuration interaction
journal, February 1973
- McLean, A. D.; Liu, B.
- The Journal of Chemical Physics, Vol. 58, Issue 3
Analytic gradients for the single-reference driven similarity renormalization group second-order perturbation theory
journal, July 2019
- Wang, Shuhe; Li, Chenyang; Evangelista, Francesco A.
- The Journal of Chemical Physics, Vol. 151, Issue 4
Multiconfigurational perturbation theory with level shift — the Cr2 potential revisited
journal, October 1995
- Roos, Björn O.; Andersson, Kerstin
- Chemical Physics Letters, Vol. 245, Issue 2-3
Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory
journal, July 2012
- Hohenstein, Edward G.; Parrish, Robert M.; Martínez, Todd J.
- The Journal of Chemical Physics, Vol. 137, Issue 4
A simple method for the evaluation of the second-order-perturbation energy from external double-excitations with a CASSCF reference wavefunction
journal, April 1982
- Roos, Björn O.; Linse, Per; Siegbahn, Per E. M.
- Chemical Physics, Vol. 66, Issue 1-2
An atomic orbital-based formulation of the complete active space self-consistent field method on graphical processing units
journal, June 2015
- Hohenstein, Edward G.; Luehr, Nathan; Ufimtsev, Ivan S.
- The Journal of Chemical Physics, Vol. 142, Issue 22
Tensor Hypercontraction Second-Order Møller–Plesset Perturbation Theory: Grid Optimization and Reaction Energies
journal, June 2015
- Kokkila Schumacher, Sara I. L.; Hohenstein, Edward G.; Parrish, Robert M.
- Journal of Chemical Theory and Computation, Vol. 11, Issue 7
Tensor hypercontraction. II. Least-squares renormalization
journal, December 2012
- Parrish, Robert M.; Hohenstein, Edward G.; Martínez, Todd J.
- The Journal of Chemical Physics, Vol. 137, Issue 22
Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation
journal, January 2008
- Ufimtsev, Ivan S.; Martínez, Todd J.
- Journal of Chemical Theory and Computation, Vol. 4, Issue 2
Direct configuration interaction with a reference state composed of many reference configurations
journal, November 1980
- Siegbahn, Per E. M.
- International Journal of Quantum Chemistry, Vol. 18, Issue 5
Analytical Derivative Coupling for Multistate CASPT2 Theory
journal, May 2017
- Park, Jae Woo; Shiozaki, Toru
- Journal of Chemical Theory and Computation, Vol. 13, Issue 6
Elimination of energy denominators in Møller—Plesset perturbation theory by a Laplace transform approach
journal, June 1991
- Almlöf, Jan
- Chemical Physics Letters, Vol. 181, Issue 4
A perspective on the CASPT2 method
journal, April 2011
- Pulay, Peter
- International Journal of Quantum Chemistry, Vol. 111, Issue 13
Reduced scaling CASPT2 using supporting subspaces and tensor hyper-contraction
journal, July 2018
- Song, Chenchen; Martínez, Todd J.
- The Journal of Chemical Physics, Vol. 149, Issue 4
Linear scaling geometry optimisation and transition state search in hybrid delocalised internal coordinates
journal, January 2000
- Billeter, Salomon R.; Turner, Alexander J.; Thiel, Walter
- Physical Chemistry Chemical Physics, Vol. 2, Issue 10
Use of approximate integrals in ab initio theory. An application in MP2 energy calculations
journal, June 1993
- Feyereisen, Martin; Fitzgerald, George; Komornicki, Andrew
- Chemical Physics Letters, Vol. 208, Issue 5-6
Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen
journal, January 1989
- Dunning, Thom H.
- The Journal of Chemical Physics, Vol. 90, Issue 2
Second‐order perturbation theory with a complete active space self‐consistent field reference function
journal, January 1992
- Andersson, Kerstin; Malmqvist, Per‐Åke; Roos, Björn O.
- The Journal of Chemical Physics, Vol. 96, Issue 2
An efficient internally contracted multiconfiguration–reference configuration interaction method
journal, November 1988
- Werner, Hans‐Joachim; Knowles, Peter J.
- The Journal of Chemical Physics, Vol. 89, Issue 9
Coulombic potential energy integrals and approximations
journal, May 1973
- Whitten, J. L.
- The Journal of Chemical Physics, Vol. 58, Issue 10
Molecular gradients for the second-order generalized Van Vleck variant of multireference perturbation theory
journal, July 2003
- Dudley, Timothy J.; Khait, Yuriy G.; Hoffmann, Mark R.
- The Journal of Chemical Physics, Vol. 119, Issue 2
A modified definition of the zeroth-order Hamiltonian in multiconfigurational perturbation theory (CASPT2)
journal, September 2004
- Ghigo, Giovanni; Roos, Björn O.; Malmqvist, Per-Åke
- Chemical Physics Letters, Vol. 396, Issue 1-3
A direct MP2 gradient method
journal, February 1990
- Frisch, Michael J.; Head-Gordon, Martin; Pople, John A.
- Chemical Physics Letters, Vol. 166, Issue 3
GVVPT2 energy gradient using a Lagrangian formulation
journal, July 2011
- Theis, Daniel; Khait, Yuriy G.; Hoffmann, Mark R.
- The Journal of Chemical Physics, Vol. 135, Issue 4
A direct-compatible formulation of the coupled perturbed complete active space self-consistent field equations on graphical processing units
journal, May 2017
- Snyder, James W.; Fales, B. Scott; Hohenstein, Edward G.
- The Journal of Chemical Physics, Vol. 146, Issue 17
Laplace transform techniques in Mo/ller–Plesset perturbation theory
journal, January 1992
- Häser, Marco; Almlöf, Jan
- The Journal of Chemical Physics, Vol. 96, Issue 1
CASPT2 study of inverse sandwich-type dinuclear 3d transition metal complexes of ethylene and dinitrogen molecules: similarities and differences in geometry, electronic structure, and spin multiplicity
journal, January 2015
- Nakagaki, Masayuki; Sakaki, Shigeyoshi
- Physical Chemistry Chemical Physics, Vol. 17, Issue 25
Construction of open shell perturbation theory invariant with respect to orbital degeneracy
journal, August 1994
- Kozlowski, P. M.; Davidson, E. R.
- Chemical Physics Letters, Vol. 226, Issue 5-6
A climbing image nudged elastic band method for finding saddle points and minimum energy paths
journal, December 2000
- Henkelman, Graeme; Uberuaga, Blas P.; Jónsson, Hannes
- The Journal of Chemical Physics, Vol. 113, Issue 22, p. 9901-9904
Geometry optimization of atomic microclusters using inverse‐power distance coordinates
journal, December 1996
- Baker, Jon; Pulay, Peter
- The Journal of Chemical Physics, Vol. 105, Issue 24
Generalized Møller—Plesset perturbation theory applied to general MCSCF reference wave functions
journal, September 1991
- Murphy, Robert B.; Messmer, Richard P.
- Chemical Physics Letters, Vol. 183, Issue 5
Discovering chemistry with an ab initio nanoreactor
journal, November 2014
- Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert
- Nature Chemistry, Vol. 6, Issue 12
Nuclear Energy Gradients for Internally Contracted Complete Active Space Second-Order Perturbation Theory: Multistate Extensions
journal, July 2016
- Vlaisavljevich, Bess; Shiozaki, Toru
- Journal of Chemical Theory and Computation, Vol. 12, Issue 8
Pressure-Dependent Kinetics of Initial Reactions in Iso-octane Pyrolysis
journal, April 2015
- Ning, HongBo; Gong, ChunMing; Li, ZeRong
- The Journal of Physical Chemistry A, Vol. 119, Issue 18
An efficient method for the evaluation of coupling coefficients in configuration interaction calculations
journal, April 1988
- Knowles, Peter J.; Werner, Hans-Joachim
- Chemical Physics Letters, Vol. 145, Issue 6
Ab initio molecular dynamics: basic concepts, current trends and novel applications
journal, December 2002
- Tuckerman, Mark E.
- Journal of Physics: Condensed Matter, Vol. 14, Issue 50
Geometry optimization made simple with translation and rotation coordinates
journal, June 2016
- Wang, Lee-Ping; Song, Chenchen
- The Journal of Chemical Physics, Vol. 144, Issue 21
Extended multi-configuration quasi-degenerate perturbation theory: The new approach to multi-state multi-reference perturbation theory
journal, June 2011
- Granovsky, Alexander A.
- The Journal of Chemical Physics, Vol. 134, Issue 21
The Complete Active Space Self-Consistent Field Method and its Applications in Electronic Structure Calculations
book, January 1987
- Roos, Björn O.
- Advances in Chemical Physics
Communication: Extended multi-state complete active space second-order perturbation theory: Energy and nuclear gradients
journal, August 2011
- Shiozaki, Toru; Győrffy, Werner; Celani, Paolo
- The Journal of Chemical Physics, Vol. 135, Issue 8
A complete active space SCF method (CASSCF) using a density matrix formulated super-CI approach
journal, May 1980
- Roos, Björn O.; Taylor, Peter R.; Sigbahn, Per E. M.
- Chemical Physics, Vol. 48, Issue 2
Partitioning technique, perturbation theory, and rational approximations: PARTITIONING TECHNIQUE
journal, January 1982
- Löwdin, Per-Olov
- International Journal of Quantum Chemistry, Vol. 21, Issue 1
Multireference perturbation theory for large restricted and selected active space reference wave functions
journal, April 2000
- Celani, Paolo; Werner, Hans-Joachim
- The Journal of Chemical Physics, Vol. 112, Issue 13
A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency
journal, August 2002
- Weigend, Florian
- Physical Chemistry Chemical Physics, Vol. 4, Issue 18
Quantum Chemistry on Graphical Processing Units. 3. Analytical Energy Gradients, Geometry Optimization, and First Principles Molecular Dynamics
journal, August 2009
- Ufimtsev, Ivan S.; Martinez, Todd J.
- Journal of Chemical Theory and Computation, Vol. 5, Issue 10
Analytic energy gradients for multiconfigurational self-consistent field second-order quasidegenerate perturbation theory (MC-QDPT)
journal, April 1998
- Nakano, Haruyuki; Hirao, Kimihiko; Gordon, Mark S.
- The Journal of Chemical Physics, Vol. 108, Issue 14
Semi-direct algorithms for the MP2 energy and gradient
journal, February 1990
- Frisch, Michael J.; Head-Gordon, Martin; Pople, John A.
- Chemical Physics Letters, Vol. 166, Issue 3
Analytical energy gradients for internally contracted second-order multireference perturbation theory
journal, September 2003
- Celani, Paolo; Werner, Hans-Joachim
- The Journal of Chemical Physics, Vol. 119, Issue 10
State-specific multireference perturbation theory: development and present status: State-specific multireference perturbation theory
journal, February 2016
- Chattopadhyay, Sudip; Chaudhuri, Rajat K.; Mahapatra, Uttam Sinha
- Wiley Interdisciplinary Reviews: Computational Molecular Science, Vol. 6, Issue 3
Imaginary Shift in CASPT2 Nuclear Gradient and Derivative Coupling Theory
journal, May 2019
- Park, Jae Woo; Al-Saadon, Rachael; Strand, Nils E.
- Journal of Chemical Theory and Computation, Vol. 15, Issue 7
Analytical gradients for tensor hyper-contracted MP2 and SOS-MP2 on graphical processing units
journal, October 2017
- Song, Chenchen; Martínez, Todd J.
- The Journal of Chemical Physics, Vol. 147, Issue 16
Communication: Automatic code generation enables nuclear gradient computations for fully internally contracted multireference theory
journal, February 2015
- MacLeod, Matthew K.; Shiozaki, Toru
- The Journal of Chemical Physics, Vol. 142, Issue 5
Analytic energy derivatives for coupled-cluster methods describing excited states: General formulas and comparison of computational costs
journal, July 1995
- Szalay, P�ter G.
- International Journal of Quantum Chemistry, Vol. 55, Issue 2
An atomic orbital-based formulation of analytical gradients and nonadiabatic coupling vector elements for the state-averaged complete active space self-consistent field method on graphical processing units
journal, October 2015
- Snyder, James W.; Hohenstein, Edward G.; Luehr, Nathan
- The Journal of Chemical Physics, Vol. 143, Issue 15
Second-order perturbation theory with a CASSCF reference function
journal, July 1990
- Andersson, Kerstin.; Malmqvist, Per Aake.; Roos, Bjoern O.
- The Journal of Physical Chemistry, Vol. 94, Issue 14
Analytical energy gradients for second-order multireference perturbation theory using density fitting
journal, March 2013
- Győrffy, Werner; Shiozaki, Toru; Knizia, Gerald
- The Journal of Chemical Physics, Vol. 138, Issue 10
Analytical Gradient Theory for Strongly Contracted (SC) and Partially Contracted (PC) N -Electron Valence State Perturbation Theory (NEVPT2)
journal, September 2019
- Park, Jae Woo
- Journal of Chemical Theory and Computation, Vol. 15, Issue 10
Quantum Chemistry on Graphical Processing Units. 2. Direct Self-Consistent-Field Implementation
journal, March 2009
- Ufimtsev, Ivan S.; Martinez, Todd J.
- Journal of Chemical Theory and Computation, Vol. 5, Issue 4
Multireference Møller—Plesset method
journal, March 1992
- Hirao, K.
- Chemical Physics Letters, Vol. 190, Issue 3-4
Analytical Gradient Theory for Quasidegenerate N -Electron Valence State Perturbation Theory (QD-NEVPT2)
journal, November 2019
- Park, Jae Woo
- Journal of Chemical Theory and Computation, Vol. 16, Issue 1
Atomic orbital-based SOS-MP2 with tensor hypercontraction. I. GPU-based tensor construction and exploiting sparsity
journal, May 2016
- Song, Chenchen; Martínez, Todd J.
- The Journal of Chemical Physics, Vol. 144, Issue 17
Toward an Understanding of Diamond sp 2 -Defects with Unsaturated Diamondoid Oligomer Models
journal, May 2015
- Zhuk, Tatyana S.; Koso, Tatyana; Pashenko, Alexander E.
- Journal of the American Chemical Society, Vol. 137, Issue 20