On the Efficiency of Algorithms for Solving Hartree-Fock and Kohn-Sham Response Equations.

The response equations as occurring in the Hartree-Fock, multiconfigurational self-consistent field, and Kohn-Sham density functional theory have identical matrix structures. The algorithms that are used for solving these equations are discussed, and new algorithms are proposed where trial vectors are split into symmetric and antisymmetric components. Numerical examples are given to compare the performance of the algorithms. The calculations show that the standard response equation for frequencies smaller than the highest occupied molecular orbital-lowest unoccupied molecular orbital gap is best solved using the preconditioned conjugate gradient or conjugate residual algorithms where trial vectors are split into symmetric and antisymmetric components. For larger frequencies in the standard response equation as well as in the damped response equation in general, the preconditioned iterative subspace approach with symmetrized trial vectors should be used. For the response eigenvalue equation, the Davidson algorithm with either paired or symmetrized trial vectors constitutes equally good options.

[1]  Poul Jørgensen,et al.  Time‐dependent Hartree–Fock calculations in the Pariser–Parr–Pople model. Applications to aniline, azulene and pyridine , 1970 .

[2]  Trond Saue,et al.  Linear response at the 4-component relativistic level: Application to the frequency-dependent dipole polarizabilities of the coinage metal dimers , 2003 .

[3]  D. M. Bishop,et al.  Nonlinear response theory with relaxation: the first-order hyperpolarizability. , 2005, The Journal of chemical physics.

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

[5]  Radovan Bast,et al.  Relativistic adiabatic time-dependent density functional theory using hybrid functionals and noncollinear spin magnetization , 2009 .

[6]  Brian J. Orr,et al.  Perturbation theory of the non-linear optical polarization of an isolated system , 1971 .

[7]  Sten Rettrup,et al.  Large-scaleRPA calculations of chiroptical properties of organic molecules: Program RPAC , 1983 .

[8]  J. Olsen,et al.  Passing the one-billion limit in full configuration-interaction (FCI) calculations , 1990 .

[9]  J. Flament,et al.  Equations‐of‐motion method: Calculation of the k lowest or highest solutions , 1979 .

[10]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[11]  J. MacDonald,et al.  Successive Approximations by the Rayleigh-Ritz Variation Method , 1933 .

[12]  P. Wormer,et al.  Conjugate gradient method for the solution of linear equations: Application to molecular electronic structure calculations , 1982 .

[13]  E. Stiefel,et al.  Relaxationsmethoden bester Strategie zur Lösung linearer Gleichungssysteme , 1955 .

[14]  J. Olsen,et al.  Solution of the large matrix equations which occur in response theory , 1988 .

[15]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[16]  Elizabeth D. Getzoff,et al.  Structure at 0.85 Å resolution of an early protein photocycle intermediate , 1998, Nature.

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

[18]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[19]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[20]  Josef Paldus,et al.  Stability Conditions for the Solutions of the Hartree—Fock Equations for Atomic and Molecular Systems. Application to the Pi‐Electron Model of Cyclic Polyenes , 1967 .

[21]  D. M. Bishop,et al.  Near-resonant absorption in the time-dependent self-consistent field and multiconfigurational self-consistent field approximations , 2001 .

[22]  Rodney J. Bartlett,et al.  The reduced linear equation method in coupled cluster theory. , 1981 .

[23]  J. Hopcroft,et al.  Triangular Factorization and Inversion by Fast Matrix Multiplication , 1974 .

[24]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[25]  Trygve Helgaker,et al.  Hartree–Fock and Kohn–Sham atomic-orbital based time-dependent response theory , 2000 .

[26]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[27]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

[28]  Trygve Helgaker,et al.  Hartree-Fock and Kohn-Sham time-dependent response theory in a second-quantization atomic-orbital formalism suitable for linear scaling. , 2008, The Journal of chemical physics.