Analytical Gradients of the Second-Order Moller-Plesset Energy Using Cholesky Decompositions

An algorithm for computing analytical gradients of the second-order MOller-Plesset (MP2) energy using density fitting (DF) is presented. The algorithm assumes that the underlying canonical Hartree-Fock reference is obtained with the same auxiliary basis set, which we obtain by Cholesky decomposition (CD) of atomic electron repulsion integrals. CD is also used for the negative semidefinite MP2 amplitude matrix. Test calculations on the weakly interacting dimers of the S22 test set (Jureka et al., Phys. Chem. Chem. Phys. 2006, 8, 1985) show that the geometry errors due to the auxiliary basis set are negligible. With double-zeta basis sets, the error due to the DF approximation in intermolecular bond lengths is better than 0.1 pm. The computational time is typically reduced by a factor of 6-7. (c) 2013 Wiley Periodicals, Inc. (Less)

[1]  Florian Weigend,et al.  Hartree–Fock exchange fitting basis sets for H to Rn † , 2008, J. Comput. Chem..

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

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

[4]  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.

[5]  Martin Head-Gordon,et al.  Scaled opposite-spin second order Møller-Plesset correlation energy: an economical electronic structure method. , 2004, The Journal of chemical physics.

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

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

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

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

[10]  D. Bernholdt,et al.  Large-scale correlated electronic structure calculations: the RI-MP2 method on parallel computers , 1996 .

[11]  Mark S. Gordon,et al.  A derivation of the frozen-orbital unrestricted open-shell and restricted closed-shell second-order perturbation theory analytic gradient expressions , 2003 .

[12]  Frederick R Manby,et al.  Analytical energy gradients for local second-order Møller-Plesset perturbation theory using density fitting approximations. , 2004, The Journal of chemical physics.

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

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

[15]  Michio Katouda,et al.  Efficient parallel algorithm of second‐order Møller–Plesset perturbation theory with resolution‐of‐identity approximation (RI‐MP2) , 2009 .

[16]  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.

[17]  Jirí Cerný,et al.  Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. , 2006, Physical chemistry chemical physics : PCCP.

[18]  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.

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

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

[21]  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.

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

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

[24]  Yihan Shao,et al.  An improved algorithm for analytical gradient evaluation in resolution‐of‐the‐identity second‐order Møller‐Plesset perturbation theory: Application to alanine tetrapeptide conformational analysis , 2007, J. Comput. Chem..

[25]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[26]  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.

[27]  Jonas Boström,et al.  Analytical Gradients of Hartree-Fock Exchange with Density Fitting Approximations. , 2013, Journal of chemical theory and computation.

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

[29]  Holger Patzelt,et al.  RI-MP2: optimized auxiliary basis sets and demonstration of efficiency , 1998 .

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

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

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

[33]  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.

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