Towards excitation energies and (hyper)polarizability calculations of large molecules. Application of parallelization and linear scaling techniques to time‐dependent density functional response theory

We document recent improvements in the efficiency of our implementation in the Amsterdam Density Functional program (ADF) of the response equations in time‐dependent density functional theory (TDDFT). Applications to quasi one‐dimensional polyene chains and to three‐dimensional water clusters show that, using our all‐electron atomic orbital (AO)‐based implementation, calculations of excitation energies and (hyper)polarizabilities on molecules with several hundred atoms and several thousand basis functions are now feasible, even on (a small cluster of) personal computers. The matrix elements, which are required in TDDFT, are calculated on an AO basis and the same linear scaling techniques as used in ADF for the iterative solution of the Kohn–Sham (KS) equations are applied to the determination of these matrix elements. Near linear scaling is demonstrated for this part of the calculation, which used to be the time‐determining step. Transformations from the AO basis to the KS orbital basis and back exhibit N3 scaling, but due to a very small prefactor this N3 scaling is still of little importance for currently accessible system sizes. The main CPU bottleneck in our current implementation is the multipolar part of the Coulomb potential, scaling quadratically with the system size. It is shown that the parallelization of our code leads to further significant reductions in execution times, with a measured speed‐up of 70 on 90 processors for both the SCF and the TDDFT parts of the code. This brings high‐level calculations on excitation energies and dynamic (hyper)polarizabilities of large molecules within reach. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1511–1523, 2000

[1]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[2]  M. Kappes,et al.  Experiment versus Time Dependent Density Functional Theory Prediction of Fullerene Electronic Absorption , 1998 .

[3]  Jack J. Dongarra,et al.  Automatically Tuned Linear Algebra Software , 1998, Proceedings of the IEEE/ACM SC98 Conference.

[4]  J. G. Snijders,et al.  Improved density functional theory results for frequency‐dependent polarizabilities, by the use of an exchange‐correlation potential with correct asymptotic behavior , 1996 .

[5]  Evert Jan Baerends,et al.  Accurate density functional calculations on frequency-dependent hyperpolarizabilities of small molecules , 1998 .

[6]  George F. Bertsch,et al.  Time-dependent local-density approximation in real time , 1996 .

[7]  M. E. Casida.  Time-Dependent Density Functional Response Theory for Molecules , 1995 .

[8]  Benoît Champagne,et al.  Electric field dependence of the exchange-correlation potential in molecular chains , 1999 .

[9]  E. Baerends,et al.  On the origin of the field-counteracting term of the Kohn-Sham exchange-correlation potential of molecular chains in an electric field , 2000 .

[10]  J. G. Snijders,et al.  Time-dependent density functional result for the dynamic hyperpolarizabilities of C60. , 1997 .

[11]  J. G. Snijders,et al.  Implementation of time-dependent density functional response equations , 1999 .

[12]  J. G. Snijders,et al.  APPLICATION OF TIME-DEPENDENT DENSITY FUNCTIONAL RESPONSE THEORY TO RAMAN SCATTERING , 1996 .

[13]  Evert Jan Baerends,et al.  Molecular calculations of excitation energies and (hyper)polarizabilities with a statistical average of orbital model exchange-correlation potentials , 2000 .

[14]  J. G. Snijders,et al.  EXCITATION ENERGIES FOR TRANSITION METAL COMPOUNDS FROM TIME-DEPENDENT DENSITY FUNCTIONAL THEORY. APPLICATIONS TO MNO4-, NI(CO)4, AND MN2(CO)10 , 1999 .

[15]  Evert Jan Baerends,et al.  A density-functional theory study of frequency-dependent polarizabilities and Van der Waals dispersion coefficients for polyatomic molecules , 1995 .

[16]  Marco Häser,et al.  CALCULATION OF EXCITATION ENERGIES WITHIN TIME-DEPENDENT DENSITY FUNCTIONAL THEORY USING AUXILIARY BASIS SET EXPANSIONS , 1997 .

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

[18]  N. Handy,et al.  Excitation energies of benzene from Kohn–Sham theory , 1999 .

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

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

[21]  Nicholas C. Handy,et al.  Improving virtual Kohn-Sham orbitals and eigenvalues: Application to excitation energies and static polarizabilities , 1998 .

[22]  Eric Schwegler,et al.  Linear scaling computation of the Fock matrix. IV. Multipole accelerated formation of the exchange matrix , 1999 .

[23]  D. Sundholm.  Density functional theory calculations of the visible spectrum of chlorophyll a , 1999 .

[24]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[25]  Dennis R. Salahub,et al.  Dynamic polarizabilities and excitation spectra from a molecular implementation of time‐dependent density‐functional response theory: N2 as a case study , 1996 .

[26]  Nicholas C. Handy,et al.  Density functional calculations of the hyperpolarisabilities of small molecules , 1999 .

[27]  E. Baerends,et al.  A Density Functional Study of the Optical Spectra and Nonlinear Optical Properties of Heteroleptic Tetrapyrrole Sandwich Complexes: The Porphyrinato−Porphyrazinato−Zirconium(IV) Complex as a Case Study , 2000 .

[28]  Vincenzo Barone,et al.  Accurate excitation energies from time-dependent density functional theory: Assessing the PBE0 model , 1999 .

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

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

[31]  Benoît Champagne,et al.  Assessment of Conventional Density Functional Schemes for Computing the Polarizabilities and Hyperpolarizabilities of Conjugated Oligomers: An Ab Initio Investigation of Polyacetylene Chains , 1998 .

[32]  Gustavo E. Scuseria,et al.  Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations , 1997 .

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

[34]  R. Ahlrichs,et al.  Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory , 1996 .

[35]  E. Baerends,et al.  Time-dependent density functional calculations on the electronic absorption spectrum of free base porphin , 1999 .

[36]  Serdar Ogut,et al.  Ab Initio Excitation Spectra and Collective Electronic Response in Atoms and Clusters , 1999 .

[37]  M. Petersilka,et al.  Excitation energies from time-dependent density-functional theory. , 1996 .