1. |
- Coriani, Sonia, et al.
(author)
-
Linear-scaling implementation of molecular response theory in self-consistent field electronic-structure theory
- 2007
-
In: Journal of Chemical Physics. - : AIP Publishing. - 0021-9606 .- 1089-7690. ; 126:15, s. 11930-11935
-
Journal article (peer-reviewed)abstract
- A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field theories for the calculation of frequency-dependent molecular response properties and excitation energies is presented, based on a nonredundant exponential parametrization of the one-electron density matrix in the atomic-orbital basis, avoiding the use of canonical orbitals. The response equations are solved iteratively, by an atomic-orbital subspace method equivalent to that of molecular-orbital theory. Important features of the subspace method are the use of paired trial vectors (to preserve the algebraic structure of the response equations), a nondiagonal preconditioner (for rapid convergence), and the generation of good initial guesses (for robust solution). As a result, the performance of the iterative method is the same as in canonical molecular-orbital theory, with five to ten iterations needed for convergence. As in traditional direct Hartree-Fock and Kohn-Sham theories, the calculations are dominated by the construction of the effective Fock/Kohn-Sham matrix, once in each iteration. Linear complexity is achieved by using sparse-matrix algebra, as illustrated in calculations of excitation energies and frequency-dependent polarizabilities of polyalanine peptides containing up to 1400 atoms.
|
|
2. |
- Salek, Pawel, et al.
(author)
-
Linear-scaling implementation of molecular electronic self-consistent field theory
- 2007
-
In: Journal of Chemical Physics. - : AIP Publishing. - 0021-9606 .- 1089-7690. ; 126:11, s. 85-98
-
Journal article (peer-reviewed)abstract
- A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field (SCF) theories is presented and illustrated with applications to molecules consisting of more than 1000 atoms. The diagonalization bottleneck of traditional SCF methods is avoided by carrying out a minimization of the Roothaan-Hall (RH) energy function and solving the Newton equations using the preconditioned conjugate-gradient (PCG) method. For rapid PCG convergence, the Lowdin orthogonal atomic orbital basis is used. The resulting linear-scaling trust-region Roothaan-Hall (LS-TRRH) method works by the introduction of a level-shift parameter in the RH Newton equations. A great advantage of the LS-TRRH method is that the optimal level shift can be determined at no extra cost, ensuring fast and robust convergence of both the SCF iterations and the level-shifted Newton equations. For density averaging, the authors use the trust-region density-subspace minimization (TRDSM) method, which, unlike the traditional direct inversion in the iterative subspace (DIIS) scheme, is firmly based on the principle of energy minimization. When combined with a linear-scaling evaluation of the Fock/Kohn-Sham matrix (including a boxed fitting of the electron density), LS-TRRH and TRDSM methods constitute the linear-scaling trust-region SCF (LS-TRSCF) method. The LS-TRSCF method compares favorably with the traditional SCF/DIIS scheme, converging smoothly and reliably in cases where the latter method fails. In one case where the LS-TRSCF method converges smoothly to a minimum, the SCF/DIIS method converges to a saddle point.
|
|