skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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 » Davidson type solvers by a factor of two in both solution time and storage.« less

Authors:
 [1];  [1];  [1];  [2];  [1]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Division
  2. 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:
Journal Article: 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. 2017. "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},
url = {https://www.osti.gov/biblio/1425428}, journal = {Computer Physics Communications},
issn = {0010-4655},
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}
}

Journal Article:

Citation Metrics:
Cited by: 8 works
Citation information provided by
Web of Science

Figures / Tables:

Algorithm 1 Algorithm 1: The Davidson algorithm8 for solving the linear response eigenvalue problem.

Save / Share:

Works referenced in this record:

Stability analysis for solutions of the closed shell Kohn–Sham equation
journal, June 1996


Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT
journal, October 2015


Minimization Principles for the Linear Response Eigenvalue Problem II: Computation
journal, January 2013


On the Efficiency of Algorithms for Solving Hartree–Fock and Kohn–Sham Response Equations
journal, April 2011


An indefinite variant of LOBPCG for definite matrix pencils
journal, August 2013


Solution of the large matrix equations which occur in response theory
journal, February 1988


Passing the one-billion limit in full configuration-interaction (FCI) calculations
journal, June 1990


An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules
journal, November 1998


Advances in molecular quantum chemistry contained in the Q-Chem 4 program package
journal, September 2014


NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations
journal, September 2010


Minimization Principles for the Linear Response Eigenvalue Problem I: Theory
journal, January 2012


A Trace Minimization Algorithm for the Generalized Eigenvalue Problem
journal, December 1982


A Jacobi–Davidson Iteration Method for Linear Eigenvalue Problems
journal, April 1996


Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method
journal, January 2001


Basis selection in LOBPCG
journal, October 2006


A Block Orthogonalization Procedure with Constant Synchronization Requirements
journal, January 2002


A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix
journal, June 2015


Density‐functional thermochemistry. III. The role of exact exchange
journal, April 1993


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


The influence of polarization functions on molecular orbital hydrogenation energies
journal, January 1973


Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation
journal, February 1992


Large-scaleRPA calculations of chiroptical properties of organic molecules: Program RPAC
journal, February 1983


A classical mechanics technique for quantum linear response
journal, June 2001


Vibrational states of nuclei in the random phase approximation
journal, January 1961


Variational procedure and generalized Lanczos recursion for small-amplitude classical oscillations
journal, December 1999


Works referencing / citing this record:

S erenity : A subsystem quantum chemistry program
journal, January 2018


A FEAST Algorithm for the Linear Response Eigenvalue Problem
journal, August 2019