Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory
Abstract
Within this paper, we present two efficient iterative algorithms for solving the linear response eigenvalue problem arising from the time dependent density functional theory. Although the matrix to be diagonalized is nonsymmetric, it has a special structure that can be exploited to save both memory and floating point operations. In particular, the nonsymmetric eigenvalue problem can be transformed into an eigenvalue problem that involves the product of two matrices M and K. We show that, because MK is self-adjoint with respect to the inner product induced by the matrix K, this product eigenvalue problem can be solved efficiently by a modified Davidson algorithm and a modified locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm that make use of the K-inner product. Additionally, the solution of the product eigenvalue problem yields one component of the eigenvector associated with the original eigenvalue problem. We show that the other component of the eigenvector can be easily recovered in an inexpensive postprocessing procedure. As a result, the algorithms we present here become more efficient than existing methods that try to approximate both components of the eigenvectors simultaneously. In particular, our numerical experiments demonstrate that the new algorithms presented here consistently outperform the existing state-of-the-artmore »
- Authors:
-
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Division
- Pacific Northwest National Lab. (PNNL), Richland, WA (United States). Environmental Molecular Sciences Lab.
- Publication Date:
- Research Org.:
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Basic Energy Sciences (BES); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); Ministry of Education, Youth and Sports (Cambodia); USDOE Office of Science (SC), Biological and Environmental Research (BER)
- OSTI Identifier:
- 1425428
- Alternate Identifier(s):
- OSTI ID: 1395271; OSTI ID: 1576835
- Report Number(s):
- PNNL-SA-114405
Journal ID: ISSN 0010-4655; ark:/13030/qt8rf6x4kg; TRN: US1802105
- Grant/Contract Number:
- AC02-05CH11231; AC02-05CH1123; KC-030106062653; AC05-76RL01830; AC05-76RL1830
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Computer Physics Communications
- Additional Journal Information:
- Journal Volume: 221; Journal Issue: C; Journal ID: ISSN 0010-4655
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; Time dependent density functional theory; Linear response eigenvalue problem; Preconditioned eigensolvers
Citation Formats
Vecharynski, Eugene, Brabec, Jiri, Shao, Meiyue, Govind, Niranjan, and Yang, Chao. Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory. United States: N. p., 2017.
Web. doi:10.1016/j.cpc.2017.07.017.
Vecharynski, Eugene, Brabec, Jiri, Shao, Meiyue, Govind, Niranjan, & Yang, Chao. Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory. United States. https://doi.org/10.1016/j.cpc.2017.07.017
Vecharynski, Eugene, Brabec, Jiri, Shao, Meiyue, Govind, Niranjan, and Yang, Chao. Fri .
"Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory". United States. https://doi.org/10.1016/j.cpc.2017.07.017. https://www.osti.gov/servlets/purl/1425428.
@article{osti_1425428,
title = {Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory},
author = {Vecharynski, Eugene and Brabec, Jiri and Shao, Meiyue and Govind, Niranjan and Yang, Chao},
abstractNote = {Within this paper, we present two efficient iterative algorithms for solving the linear response eigenvalue problem arising from the time dependent density functional theory. Although the matrix to be diagonalized is nonsymmetric, it has a special structure that can be exploited to save both memory and floating point operations. In particular, the nonsymmetric eigenvalue problem can be transformed into an eigenvalue problem that involves the product of two matrices M and K. We show that, because MK is self-adjoint with respect to the inner product induced by the matrix K, this product eigenvalue problem can be solved efficiently by a modified Davidson algorithm and a modified locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm that make use of the K-inner product. Additionally, the solution of the product eigenvalue problem yields one component of the eigenvector associated with the original eigenvalue problem. We show that the other component of the eigenvector can be easily recovered in an inexpensive postprocessing procedure. As a result, the algorithms we present here become more efficient than existing methods that try to approximate both components of the eigenvectors simultaneously. In particular, our numerical experiments demonstrate that the new algorithms presented here consistently outperform the existing state-of-the-art Davidson type solvers by a factor of two in both solution time and storage.},
doi = {10.1016/j.cpc.2017.07.017},
journal = {Computer Physics Communications},
number = C,
volume = 221,
place = {United States},
year = {Fri Dec 01 00:00:00 EST 2017},
month = {Fri Dec 01 00:00:00 EST 2017}
}
Web of Science
Figures / Tables:
Works referenced in this record:
Stability analysis for solutions of the closed shell Kohn–Sham equation
journal, June 1996
- Bauernschmitt, Rüdiger; Ahlrichs, Reinhart
- The Journal of Chemical Physics, Vol. 104, Issue 22
Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT
journal, October 2015
- Brabec, Jiri; Lin, Lin; Shao, Meiyue
- Journal of Chemical Theory and Computation, Vol. 11, Issue 11
Minimization Principles for the Linear Response Eigenvalue Problem II: Computation
journal, January 2013
- Bai, Zhaojun; Li, Ren-Cang
- SIAM Journal on Matrix Analysis and Applications, Vol. 34, Issue 2
Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations
journal, March 2014
- Challacombe, Matt
- Computation, Vol. 2, Issue 1
On the Efficiency of Algorithms for Solving Hartree–Fock and Kohn–Sham Response Equations
journal, April 2011
- Kauczor, Joanna; Jørgensen, Poul; Norman, Patrick
- Journal of Chemical Theory and Computation, Vol. 7, Issue 6
An indefinite variant of LOBPCG for definite matrix pencils
journal, August 2013
- Kressner, Daniel; Pandur, Marija Miloloža; Shao, Meiyue
- Numerical Algorithms, Vol. 66, Issue 4
Solution of the large matrix equations which occur in response theory
journal, February 1988
- Olsen, Jeppe; Jensen, Hans Jørgen Aa; Jørgensen, Poul
- Journal of Computational Physics, Vol. 74, Issue 2
Passing the one-billion limit in full configuration-interaction (FCI) calculations
journal, June 1990
- Olsen, Jeppe; Jørgensen, Poul; Simons, Jack
- Chemical Physics Letters, Vol. 169, Issue 6
An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules
journal, November 1998
- Stratmann, R. Eric; Scuseria, Gustavo E.; Frisch, Michael J.
- The Journal of Chemical Physics, Vol. 109, Issue 19
Advances in molecular quantum chemistry contained in the Q-Chem 4 program package
journal, September 2014
- Shao, Yihan; Gan, Zhengting; Epifanovsky, Evgeny
- Molecular Physics, Vol. 113, Issue 2
NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations
journal, September 2010
- Valiev, M.; Bylaska, E. J.; Govind, N.
- Computer Physics Communications, Vol. 181, Issue 9, p. 1477-1489
Minimization Principles for the Linear Response Eigenvalue Problem I: Theory
journal, January 2012
- Bai, Zhaojun; Li, Ren-Cang
- SIAM Journal on Matrix Analysis and Applications, Vol. 33, Issue 4
A Trace Minimization Algorithm for the Generalized Eigenvalue Problem
journal, December 1982
- Sameh, Ahmed H.; Wisniewski, John A.
- SIAM Journal on Numerical Analysis, Vol. 19, Issue 6
The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices
journal, January 1975
- Davidson, Ernest R.
- Journal of Computational Physics, Vol. 17, Issue 1
A Jacobi–Davidson Iteration Method for Linear Eigenvalue Problems
journal, April 1996
- G. Sleijpen, Gerard L.; Van der Vorst, Henk A.
- SIAM Journal on Matrix Analysis and Applications, Vol. 17, Issue 2
An iterative method for calculating several of the extreme eigensolutions of large real non-symmetric matrices
journal, January 1982
- Rettrup, Sten
- Journal of Computational Physics, Vol. 45, Issue 1
Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method
journal, January 2001
- Knyazev, Andrew V.
- SIAM Journal on Scientific Computing, Vol. 23, Issue 2
Basis selection in LOBPCG
journal, October 2006
- Hetmaniuk, U.; Lehoucq, R.
- Journal of Computational Physics, Vol. 218, Issue 1
A Block Orthogonalization Procedure with Constant Synchronization Requirements
journal, January 2002
- Stathopoulos, Andreas; Wu, Kesheng
- SIAM Journal on Scientific Computing, Vol. 23, Issue 6
A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix
journal, June 2015
- Vecharynski, Eugene; Yang, Chao; Pask, John E.
- Journal of Computational Physics, Vol. 290
Density‐functional thermochemistry. III. The role of exact exchange
journal, April 1993
- Becke, Axel D.
- The Journal of Chemical Physics, Vol. 98, Issue 7, p. 5648-5652
Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules
journal, March 1972
- Hehre, W. J.; Ditchfield, R.; Pople, J. A.
- The Journal of Chemical Physics, Vol. 56, Issue 5, p. 2257-2261
The influence of polarization functions on molecular orbital hydrogenation energies
journal, January 1973
- Hariharan, P. C.; Pople, J. A.
- Theoretica Chimica Acta, Vol. 28, Issue 3
Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation
journal, February 1992
- Godbout, Nathalie; Salahub, Dennis R.; Andzelm, Jan
- Canadian Journal of Chemistry, Vol. 70, Issue 2
Large-scaleRPA calculations of chiroptical properties of organic molecules: Program RPAC
journal, February 1983
- Bouman, Thomas D.; Hansen, Aage E.; Voigt, Bj�rn
- International Journal of Quantum Chemistry, Vol. 23, Issue 2
A classical mechanics technique for quantum linear response
journal, June 2001
- Tsiper, E. V.
- Journal of Physics B: Atomic, Molecular and Optical Physics, Vol. 34, Issue 12
Vibrational states of nuclei in the random phase approximation
journal, January 1961
- Thouless, D. J.
- Nuclear Physics, Vol. 22, Issue 1
Variational procedure and generalized Lanczos recursion for small-amplitude classical oscillations
journal, December 1999
- Tsiper, E. V.
- Journal of Experimental and Theoretical Physics Letters, Vol. 70, Issue 11
Works referencing / citing this record:
S erenity : A subsystem quantum chemistry program
journal, January 2018
- Unsleber, Jan P.; Dresselhaus, Thomas; Klahr, Kevin
- Journal of Computational Chemistry, Vol. 39, Issue 13
A FEAST Algorithm for the Linear Response Eigenvalue Problem
journal, August 2019
- Teng, Zhongming; Lu, Linzhang
- Algorithms, Vol. 12, Issue 9