An efficient approach for self-consistent-field energy and energy second derivatives in the atomic-orbital basis.

Based on self-consistent-field (SCF) perturbation theory, we recast the SCF and the coupled-perturbed SCF (CPSCF) equations for time-independent molecular properties into the atomic-orbital basis. The density matrix and the perturbed density matrix are obtained iteratively by solving linear equations. Only matrix multiplications and additions are required, and this approach can exploit sparse matrix multiplications and thereby offer the possibility of evaluating second-order properties in computational effort that scales linearly with system size. Convergence properties are similar to conventional molecular-orbital-based CPSCF procedures, in terms of the number of derivative Fock matrices that must be constructed. We also carefully address the issue of the numerical accuracy of the calculated second derivatives of the energy, in order to specify the minimum precision necessary in the CPSCF procedure. It is found that much looser tolerances for the perturbed density matrices are adequate when using an expression for the second derivatives that is correct through second order in the CPSCF error.

[1]  Michel Dupuis,et al.  Molecular symmetry and closed‐shell SCF calculations. I , 1977 .

[2]  S. Goedecker Linear scaling electronic structure methods , 1999 .

[3]  K. Wolinski,et al.  Self-consistent perturbation theory: Open-shell states in perturbation-dependent non-orthogonal basis sets This work was partly supported by the Institute of Low Temperatures and Structure Research of the Polish Academy of Sciences under contract No. MR-I.9.4.3/2. , 1980 .

[4]  Christian Ochsenfeld,et al.  Linear scaling exchange gradients for Hartree–Fock and hybrid density functional theory , 2000 .

[5]  Matt Challacombe,et al.  Density matrix perturbation theory. , 2003, Physical review letters.

[6]  A. J. Sadlej,et al.  Molecular electric polarizabilities. Electronic-field-variant (EFV) gaussian basis set for polarizability calculations , 1977 .

[7]  Jörg Kussmann,et al.  Ab initio NMR spectra for molecular systems with a thousand and more atoms: a linear-scaling method. , 2004, Angewandte Chemie.

[8]  Valéry Weber,et al.  Ab initio linear scaling response theory: electric polarizability by perturbed projection. , 2004, Physical review letters.

[9]  Peter M. W. Gill,et al.  Molecular integrals Over Gaussian Basis Functions , 1994 .

[10]  Harry F. King,et al.  Analytic computation of energy derivatives. Relationships among partial derivatives of a variationally determined function , 1986 .

[11]  S. Yokojima,et al.  Generalized linear-scaling localized-density-matrix method , 1999 .

[12]  R. Parr Density-functional theory of atoms and molecules , 1989 .

[13]  Michael J. Frisch,et al.  Direct analytic SCF second derivatives and electric field properties , 1990 .

[14]  C. Daul,et al.  Improved coupled perturbed Hartree–Fock and Kohn–Sham convergence acceleration , 2003 .

[15]  Michel Dupuis,et al.  Numerical integration using rys polynomials , 1976 .

[16]  Yihan Shao,et al.  Efficient evaluation of the Coulomb force in density-functional theory calculations , 2001 .

[17]  J. Simons,et al.  Ab initio analytical molecular gradients and Hessians , 1983 .

[18]  Yihan Shao,et al.  Improved Fermi operator expansion methods for fast electronic structure calculations , 2003 .

[19]  P. Pulay Improved SCF convergence acceleration , 1982 .

[20]  David R. Bowler,et al.  Recent progress in linear scaling ab initio electronic structure techniques , 2002 .

[21]  W. Lipscomb,et al.  Perturbed Hartree—Fock Calculations. I. Magnetic Susceptibility and Shielding in the LiH Molecule , 1963 .

[22]  Matt Challacombe,et al.  A simplified density matrix minimization for linear scaling self-consistent field theory , 1999 .

[23]  Jon Baker,et al.  Q‐Chem 2.0: a high‐performance ab initio electronic structure program package , 2000, J. Comput. Chem..

[24]  Martin Head-Gordon,et al.  Approaching the Basis Set Limit in Density Functional Theory Calculations Using Dual Basis Sets without Diagonalization , 2004 .

[25]  Michel Dupuis,et al.  Computation of electron repulsion integrals using the rys quadrature method , 1983 .

[26]  Peter Pulay,et al.  Second and third derivatives of variational energy expressions: Application to multiconfigurational self‐consistent field wave functions , 1983 .

[27]  David E. Manolopoulos,et al.  Canonical purification of the density matrix in electronic-structure theory , 1998 .

[28]  Gustavo E. Scuseria,et al.  Linear Scaling Density Functional Calculations with Gaussian Orbitals , 1999 .

[29]  Yihan Shao,et al.  Curvy steps for density matrix-based energy minimization: Application to large-scale self-consistent-field calculations , 2003 .

[30]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .

[31]  A. J. Sadlej,et al.  Analytic techniques for the coupled multiconfiguration SCF perturbation theory , 1977 .

[32]  Paweł Sałek,et al.  Linear-scaling formation of Kohn-Sham Hamiltonian: application to the calculation of excitation energies and polarizabilities of large molecular systems. , 2004, The Journal of chemical physics.

[33]  Christian Ochsenfeld,et al.  A reformulation of the coupled perturbed self-consistent field equations entirely within a local atomic orbital density matrix-based scheme , 1997 .

[34]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[35]  Guanhua Chen,et al.  Linear-scaling time-dependent density-functional theory , 2003 .

[36]  R. Mcweeny Perturbation Theory for the Fock-Dirac Density Matrix , 1962 .

[37]  Trygve Helgaker,et al.  Direct optimization of the AO density matrix in Hartree-Fock and Kohn-Sham theories , 2000 .

[38]  H. Schaefer,et al.  Analytic second derivatives in restricted Hartree–Fock theory. A method for high‐spin open‐shell molecular wave functions , 1982 .

[39]  Michel Dupuis,et al.  Molecular symmetry. II. Gradient of electronic energy with respect to nuclear coordinates , 1978 .

[40]  P. Schleyer Encyclopedia of computational chemistry , 1998 .

[41]  G. Diercksen,et al.  Self‐Consistent Perturbation Theory. I. General Formulation and Some Applications , 1966 .

[42]  Ian M. Mills,et al.  Force Constants and Dipole-Moment Derivatives of Molecules from Perturbed Hartree-Fock Calculations. I , 1968 .

[43]  H. Schaefer,et al.  Analytic evaluation and basis set dependence of intensities of infrared spectra , 1986 .

[44]  Trygve Helgaker,et al.  Geometrical derivatives and magnetic properties in atomic-orbital density-based Hartree-Fock theory , 2001 .

[45]  Christian Ochsenfeld,et al.  Linear and sublinear scaling formation of Hartree-Fock-type exchange matrices , 1998 .

[46]  Li,et al.  Density-matrix electronic-structure method with linear system-size scaling. , 1993, Physical review. B, Condensed matter.

[47]  M. Head‐Gordon,et al.  Fast evaluation of a linear number of local exchange matrices , 2002 .