Faber and Newton polynomial integrators for open-system density matrix propagation

Wilhelm Huisinga*, Lorenzo Pesce, Ronnie Kosloff, Peter Saalfrank

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

64 Scopus citations


Two polynomial expansions of the time-evolution superoperator to directly integrate Markovian Liouville-von Neumann (LvN) equations for quantum open systems, namely the Newton interpolation and the Faber approximation, are presented and critically compared. Details on the numerical implementation including error control, and on the performance of either method are given. In a first physical application, a damped harmonic oscillator is considered. Then, the Faber approximation is applied to compute a condensed phase absorption spectrum, for which a semianalytical expression is derived. Finally, even more general applications are discussed. In all applications considered here it is found that both the Newton and Faber integrators are fast, general, stable, and accurate.

Original languageEnglish (US)
Pages (from-to)5538-5547
Number of pages10
JournalJournal of Chemical Physics
Issue number12
StatePublished - Mar 22 1999

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Fingerprint Dive into the research topics of 'Faber and Newton polynomial integrators for open-system density matrix propagation'. Together they form a unique fingerprint.

Cite this