Cholesky Decomposition Techniques in Electronic Structure Theory

We review recently developed methods to efficiently utilize the Cholesky decomposition technique in electronic structure calculations. The review starts with a brief introduction to the basics of the Cholesky decomposition technique. Subsequently, examples of applications of the technique to ab inito procedures are presented. The technique is demonstrated to be a special type of a resolution-of-identity or density-fitting scheme. This is followed by explicit examples of the Cholesky techniques used in orbital localization, computation of the exchange contribution to the Fock matrix, in MP2, gradient calculations, and so-called method specific Cholesky decomposition. Subsequently, examples of calibration of the method with respect to computed total energies, excitation energies, and auxiliary basis set pruning are presented. In particular, it is demonstrated that the Cholesky method is an unbiased method to derive auxiliary basis sets. Furthermore, details of the implementational considerations are put forward and examples from a parallel Cholesky decomposition scheme is presented. Finally, an outlook and perspectives are presented, followed by a summary and conclusions section. We are of the opinion that the Cholesky decomposition method is a technique that has been overlooked for too long. We have just recently started to understand how to efficiently incorporate the method in existing ab initio programs. The full potential of the Cholesky technique has not yet been fully explored.

[1]  T. Helgaker,et al.  Variational and robust density fitting of four-center two-electron integrals in local metrics. , 2008, The Journal of chemical physics.

[2]  Donald G Truhlar,et al.  Density functionals for inorganometallic and organometallic chemistry. , 2005, The journal of physical chemistry. A.

[3]  Jonas Boström,et al.  Ab Initio Density Fitting: Accuracy Assessment of Auxiliary Basis Sets from Cholesky Decompositions. , 2009, Journal of chemical theory and computation.

[4]  Thomas Bondo Pedersen,et al.  Reduced scaling in electronic structure calculations using Cholesky decompositions , 2003 .

[5]  I. Røeggen Analytic functions for the three-body potential of the helium trimer. , 2007, The Journal of chemical physics.

[6]  U. Ryde,et al.  Reaction mechanism of manganese superoxide dismutase studied by combined quantum and molecular mechanical calculations and multiconfigurational methods. , 2009, The journal of physical chemistry. B.

[7]  Gene H. Golub,et al.  Matrix computations , 1983 .

[8]  Paul G. Mezey,et al.  A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions , 1989 .

[9]  F. Weigend,et al.  RI-MP2: first derivatives and global consistency , 1997 .

[10]  Walter Thiel,et al.  Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3. , 2008, The Journal of chemical physics.

[11]  H. Koch,et al.  Cholesky decomposition-based definition of atomic subsystems in electronic structure calculations. , 2010, The Journal of chemical physics.

[12]  J. Simons,et al.  Application of cholesky-like matrix decomposition methods to the evaluation of atomic orbital integrals and integral derivatives , 1989 .

[13]  Branislav Jansík,et al.  Maximum locality in occupied and virtual orbital spaces using a least-change strategy. , 2009, The Journal of chemical physics.

[14]  J. Whitten,et al.  Electron repulsion integral approximations and error bounds: Molecular applications , 1974 .

[15]  Martin W. Feyereisen,et al.  Use of approximate integrals in ab initio theory. An application in MP2 energy calculations , 1993 .

[16]  C. Cramer,et al.  What Active Space Adequately Describes Oxygen Activation by a Late Transition Metal? CASPT2 and RASPT2 Applied to Intermediates from the Reaction of O2 with a Cu(I)-α-Ketocarboxylate. , 2009, Journal of chemical theory and computation.

[17]  J. L. Whitten,et al.  Coulombic potential energy integrals and approximations , 1973 .

[18]  Martin Head-Gordon,et al.  Auxiliary basis expansions for large-scale electronic structure calculations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[19]  A. St.-Amant,et al.  Linear scaling for the charge density fitting procedure of the linear combination of Gaussian-type orbitals density functional method , 1996 .

[20]  Henrik Koch,et al.  Method specific Cholesky decomposition: coulomb and exchange energies. , 2008, The Journal of chemical physics.

[21]  F. Weigend Accurate Coulomb-fitting basis sets for H to Rn. , 2006, Physical chemistry chemical physics : PCCP.

[22]  Francesco Aquilante,et al.  Fast noniterative orbital localization for large molecules. , 2006, The Journal of chemical physics.

[23]  Andreas Heßelmann,et al.  Accurate Intermolecular Interaction Energies from a Combination of MP2 and TDDFT Response Theory. , 2010, Journal of chemical theory and computation.

[24]  B. I. Dunlap,et al.  Robust variational fitting: Gáspár's variational exchange can accurately be treated analytically , 2000 .

[25]  Evert Jan Baerends,et al.  Self-consistent molecular Hartree—Fock—Slater calculations I. The computational procedure , 1973 .

[26]  T. Crawford,et al.  Gas phase optical rotation calculated from coupled cluster theory with zero-point vibrational corrections from density funcional theory. , 2009, Chirality.

[27]  Ronald H. Felton,et al.  A new computational approach to Slater’s SCF–Xα equation , 1975 .

[28]  Valera Veryazov,et al.  Not innocent: verdict from ab initio multiconfigurational second-order perturbation theory on the electronic structure of chloroiron corrole. , 2008, The journal of physical chemistry. B.

[29]  Patrizia Calaminici,et al.  Robust and efficient density fitting. , 2009, The Journal of chemical physics.

[30]  H. Koch,et al.  Carbon nanorings: A challenge to theoretical chemistry. , 2006, Chemphyschem : a European journal of chemical physics and physical chemistry.

[31]  C. de Graaf,et al.  Analysis of the magnetic coupling in M3(dpa)4Cl2 systems (M = Ni, Pd, Cu, Ag) by ab initio calculations. , 2010, The journal of physical chemistry. A.

[32]  I. Røeggen An ab initio study of the fcc and hcp structures of helium. , 2006, The Journal of chemical physics.

[33]  Roland Lindh,et al.  Accurate ab initio density fitting for multiconfigurational self-consistent field methods. , 2008, The Journal of chemical physics.

[34]  John R. Sabin,et al.  On some approximations in applications of Xα theory , 1979 .

[35]  S. F. Boys Construction of Some Molecular Orbitals to Be Approximately Invariant for Changes from One Molecule to Another , 1960 .

[36]  T. Crawford,et al.  On the importance of vibrational contributions to small-angle optical rotation: Fluoro-oxirane in gas phase and solution. , 2009, The Journal of chemical physics.

[37]  Florian Janetzko,et al.  A MinMax self-consistent-field approach for auxiliary density functional theory. , 2009, The Journal of chemical physics.

[38]  Christof Hättig,et al.  Analytic gradients for excited states in the coupled-cluster model CC2 employing the resolution-of-the-identity approximation , 2003 .

[39]  I Røeggen,et al.  Cholesky decomposition of the two-electron integral matrix in electronic structure calculations. , 2008, The Journal of chemical physics.

[40]  Joseph E. Subotnik,et al.  Linear scaling density fitting. , 2006, The Journal of chemical physics.

[41]  N. H. Beebe,et al.  Simplifications in the generation and transformation of two‐electron integrals in molecular calculations , 1977 .

[42]  H. Koch,et al.  Benzene-argon S1 intermolecular potential energy surface , 1999 .

[43]  H. Koch,et al.  Variation of polarizability in the [4n+2] annulene series: from [22]- to [66]-annulene. , 2008, Physical chemistry chemical physics : PCCP.

[44]  H. Koch,et al.  Coupled cluster calculations of interaction energies in benzene–fluorobenzene van der Waals complexes , 2007 .

[45]  Klaus Ruedenberg,et al.  Localized Atomic and Molecular Orbitals , 1963 .

[46]  I. Røeggen An Extended Group Function Model for Intermolecular Interactions , 2006 .

[47]  Roland Lindh,et al.  Density fitting with auxiliary basis sets from Cholesky decompositions , 2009 .

[48]  F. Manby,et al.  Explicitly correlated second-order perturbation theory using density fitting and local approximations. , 2006, The Journal of chemical physics.

[49]  Roland Lindh,et al.  Unbiased auxiliary basis sets for accurate two-electron integral approximations. , 2007, The Journal of chemical physics.

[50]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials , 1995 .

[51]  Francesco Aquilante,et al.  Cholesky Decomposition-Based Multiconfiguration Second-Order Perturbation Theory (CD-CASPT2): Application to the Spin-State Energetics of Co(III)(diiminato)(NPh). , 2008, Journal of chemical theory and computation.

[52]  R. Lindh,et al.  Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals. , 2007, The Journal of chemical physics.

[53]  Evert Jan Baerends,et al.  Towards an order , 1998 .

[54]  FRANCESCO AQUILANTE,et al.  MOLCAS 7: The Next Generation , 2010, J. Comput. Chem..

[55]  Frederick R. Manby,et al.  Fast Hartree–Fock theory using local density fitting approximations , 2004 .

[56]  F. Weigend,et al.  Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations , 2002 .

[57]  C. Cramer,et al.  The Ru-Hbpp water oxidation catalyst. , 2009, Journal of the American Chemical Society.

[58]  C. Cramer,et al.  Generating Cu(II)-oxyl/Cu(III)-oxo species from Cu(I)-alpha-ketocarboxylate complexes and O2: in silico studies on ligand effects and C-H-activation reactivity. , 2009, Chemistry.

[59]  Pavel Hobza,et al.  On the Structure and Geometry of Biomolecular Binding Motifs (Hydrogen-Bonding, Stacking, X-H···π): WFT and DFT Calculations. , 2010, Journal of chemical theory and computation.

[60]  Dmitrij Rappoport,et al.  Analytical time-dependent density functional derivative methods within the RI-J approximation, an approach to excited states of large molecules. , 2005, The Journal of chemical physics.

[61]  Per-Olov Löwdin,et al.  Studies in perturbation theory. XI - Lower bounds to energy eigenvalues, ground state, and excited states. , 1965 .

[62]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[63]  J. G. Snijders,et al.  Towards an order-N DFT method , 1998 .

[64]  Roland Lindh,et al.  Atomic Cholesky decompositions: a route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency. , 2009, The Journal of chemical physics.

[65]  L. Cederbaum,et al.  On the Cholesky decomposition for electron propagator methods: General aspects and application on C(60). , 2009, The Journal of chemical physics.

[66]  P. Neogrády,et al.  Three- and four-body nonadditivities in nucleic acid tetramers: a CCSD(T) study. , 2010, Physical chemistry chemical physics : PCCP.

[67]  Roland Lindh,et al.  Analytic derivatives for the Cholesky representation of the two-electron integrals. , 2008, The Journal of chemical physics.

[68]  Henrik Koch,et al.  Size-intensive decomposition of orbital energy denominators , 2000 .

[69]  Jörg Kussmann,et al.  Linear‐scaling Cholesky decomposition , 2008, J. Comput. Chem..

[70]  Jon Baker,et al.  The location of transition states: A comparison of Cartesian, Z‐matrix, and natural internal coordinates , 1996 .

[71]  Christof Hättig,et al.  Geometry optimizations with the coupled-cluster model CC2 using the resolution-of-the-identity approximation , 2003 .

[72]  Filipp Furche,et al.  Nuclear second analytical derivative calculations using auxiliary basis set expansions , 2004 .

[73]  K. Pierloot,et al.  Relative energy of the high-(5T2g) and low-(1A1g) spin states of the ferrous complexes [Fe(L)(NHS4)]: CASPT2 versus density functional theory. , 2008, The Journal of chemical physics.

[74]  Paweł Sałek,et al.  Linear-scaling implementation of molecular electronic self-consistent field theory. , 2007, The Journal of chemical physics.

[75]  Florian Weigend,et al.  A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency , 2002 .

[76]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[77]  H. Koch,et al.  Origin invariant calculation of optical rotation without recourse to London orbitals , 2004 .

[78]  Jonas Boström,et al.  Calibration of Cholesky Auxiliary Basis Sets for Multiconfigurational Perturbation Theory Calculations of Excitation Energies. , 2010, Journal of chemical theory and computation.

[79]  Poul Jørgensen,et al.  The second-order approximate coupled cluster singles and doubles model CC2 , 1995 .

[80]  J. Connolly,et al.  On first‐row diatomic molecules and local density models , 1979 .

[81]  Philippe Y. Ayala,et al.  Linear scaling coupled cluster and perturbation theories in the atomic orbital basis , 1999 .

[82]  L. Curtiss,et al.  Assessment of Gaussian-2 and density functional theories for the computation of enthalpies of formation , 1997 .

[83]  Henrik Koch,et al.  Polarizabilities of small annulenes from Cholesky CC2 linear response theory , 2004 .

[84]  Asbjörn M. Burow,et al.  Resolution of identity approximation for the Coulomb term in molecular and periodic systems. , 2009, The Journal of chemical physics.

[85]  Francesco Aquilante,et al.  Calculation of protein-ligand interaction energies by a fragmentation approach combining high-level quantum chemistry with classical many-body effects. , 2009, The journal of physical chemistry. B.

[86]  Thomas Bondo Pedersen,et al.  Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions. , 2004, The Journal of chemical physics.

[87]  Martin Head-Gordon,et al.  Hartree-Fock exchange computed using the atomic resolution of the identity approximation. , 2008, The Journal of chemical physics.

[88]  Emily A Carter,et al.  Cholesky decomposition within local multireference singles and doubles configuration interaction. , 2010, The Journal of chemical physics.

[89]  I. Røeggen,et al.  On the Beebe-Linderberg two-electron integral approximation , 1986 .

[90]  Marek Sierka,et al.  Fast evaluation of the Coulomb potential for electron densities using multipole accelerated resolution of identity approximation , 2003 .

[91]  J. Almlöf,et al.  Integral approximations for LCAO-SCF calculations , 1993 .

[92]  Francesco Aquilante,et al.  Quartic scaling evaluation of canonical scaled opposite spin second-order Møller Plesset correlation energy using Cholesky decompositions , 2007 .

[93]  F. Aquilante,et al.  p-Benzoquinone in aqueous solution: Stark shifts in spectra, asymmetry in solvent structure. , 2007, Physical chemistry chemical physics : PCCP.

[94]  I. Røeggen,et al.  Interatomic potential for the X1Σ +g state of Be2, revisited , 2005 .

[95]  K. Pierloot,et al.  Binding of CO, NO, and O2 to heme by density functional and multireference ab initio calculations. , 2008, The journal of physical chemistry. A.