Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory

Melissa J. Lucero, Anders M.N. Niklasson, Sergei Tretiak, Matt Challacombe

Research output: Contribution to journalArticlepeer-review

17 Citations (Scopus)


A nonlinear conjugate gradient optimization scheme is used to obtain excitation energies within the random phase approximation (RPA). The solutions to the RPA eigenvalue equation are located through a variational characterization using a modified Thouless functional, which is based upon an asymmetric Rayleigh quotient, in an orthogonalized atomic orbital representation. In this way, the computational bottleneck of calculating molecular orbitals is avoided. The variational space is reduced to the physically-relevant transitions by projections. The feasibility of an RPA implementation scaling linearly with system size N is investigated by monitoring convergence behavior with respect to the quality of initial guess and sensitivity to noise under thresholding, both for well- and ill-conditioned problems. The molecular-orbital-free algorithm is found to be robust and computationally efficient, providing a first step toward large-scale, reduced complexity calculations of time-dependent optical properties and linear response. The algorithm is extensible to other forms of time-dependent perturbation theory including, but not limited to, time-dependent density functional theory.

Original languageEnglish
Article number064114
JournalJournal of Chemical Physics
Issue number6
Publication statusPublished - 2008
Externally publishedYes


Dive into the research topics of 'Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory'. Together they form a unique fingerprint.

Cite this