Molecular integrals from Fast Fourier Transforms (FFT) instead of recurrences: The McMurchie-Davidson case

Mieke Peels, Gerald Knizia

Research output: Contribution to journalArticlepeer-review

Abstract

We report a closed formula expressing the McMurchie-Davidson (MD) key intermediates {[r](0); rx + ry + rz ≤ L} directly in terms of the set of basic integrals {[0](m); m ≤ L}, without any recurrences. This formula can be evaluated at O(L) cost per output [r](0) with dense matrix multiplications and Fast Fourier Transforms (FFT). Key to this is the fact that the transformation that builds Cartesian angular momentum from the basic integrals, {[0κ](m+m′)}→{[lκ](m)} (κ ϵ {x, y, z}), can be phrased as a circulant-matrix/vector product, which is susceptible to FFTs. After simplification, a simple formula yields the final [r](0) in one step, as contraction of four auxiliary vectors over a common Fourier index k-one vector for the [0](m) and one for each Cartesian axis. Similar transformations occur in many integral approaches beside MD, making this idea potentially broadly applicable. The simple resulting code and data structures may make it attractive for novel hardware platforms.

Original languageEnglish (US)
Article number231103
JournalJournal of Chemical Physics
Volume152
Issue number23
DOIs
StatePublished - Jun 21 2020

All Science Journal Classification (ASJC) codes

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

Fingerprint Dive into the research topics of 'Molecular integrals from Fast Fourier Transforms (FFT) instead of recurrences: The McMurchie-Davidson case'. Together they form a unique fingerprint.

Cite this