{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,7,30]],"date-time":"2025-07-30T11:15:21Z","timestamp":1753874121816,"version":"3.41.2"},"reference-count":43,"publisher":"AIP Publishing","issue":"23","license":[{"start":{"date-parts":[[2020,6,17]],"date-time":"2020-06-17T00:00:00Z","timestamp":1592352000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"},{"start":{"date-parts":[[2020,6,17]],"date-time":"2020-06-17T00:00:00Z","timestamp":1592352000000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["pubs.aip.org"],"crossmark-restriction":true},"short-container-title":[],"published-print":{"date-parts":[[2020,6,21]]},"abstract":"<jats:p>We report a closed formula expressing the McMurchie\u2013Davidson (MD) key intermediates {[r](0); rx + ry + rz \u2264 L} directly in terms of the set of basic integrals {[0](m); m \u2264 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\u03ba](m+m\u2032)}\u21a6{[l\u03ba](m)} (\u03ba \u2208 {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\u2014one 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.<\/jats:p>","DOI":"10.1063\/5.0002880","type":"journal-article","created":{"date-parts":[[2020,6,17]],"date-time":"2020-06-17T10:28:46Z","timestamp":1592389726000},"update-policy":"https:\/\/doi.org\/10.1063\/aip-crossmark-policy-page","source":"Crossref","is-referenced-by-count":2,"title":["Molecular integrals from Fast Fourier Transforms (FFT) instead of recurrences: The McMurchie\u2013Davidson case"],"prefix":"10.1063","volume":"152","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-3389-9393","authenticated-orcid":false,"given":"Mieke","family":"Peels","sequence":"first","affiliation":[{"name":"Department of Chemistry, Pennsylvania State University , University Park, Pennsylvania 16802, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7163-4823","authenticated-orcid":false,"given":"Gerald","family":"Knizia","sequence":"additional","affiliation":[{"name":"Department of Chemistry, Pennsylvania State University , University Park, Pennsylvania 16802, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"317","published-online":{"date-parts":[[2020,6,17]]},"reference":[{"key":"2023080302020032300_c1","doi-asserted-by":"publisher","first-page":"218","DOI":"10.1016\/0021-9991(78)90092-x","article-title":"One-and two-electron integrals over Cartesian Gaussian functions","volume":"26","year":"1978","journal-title":"J. Comput. Phys."},{"key":"2023080302020032300_c2","doi-asserted-by":"publisher","first-page":"3963","DOI":"10.1063\/1.450106","article-title":"Efficient recursive computation of molecular integrals over Cartesian Gaussian functions","volume":"84","year":"1986","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c3","doi-asserted-by":"publisher","first-page":"1540","DOI":"10.1063\/1.455717","article-title":"General recurrence formulas for molecular integrals over Cartesian Gaussian functions","volume":"89","year":"1988","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c4","doi-asserted-by":"publisher","first-page":"3072","DOI":"10.1039\/b605188j","article-title":"A simple algebraic derivation of the Obara\u2013Saika scheme for general two-electron interaction potentials","volume":"8","year":"2006","journal-title":"Phys. Chem. Chem. Phys."},{"key":"2023080302020032300_c5","doi-asserted-by":"publisher","first-page":"5119","DOI":"10.1039\/b413539c","article-title":"Efficient evaluation of three-center two-electron integrals over Gaussian functions","volume":"6","year":"2004","journal-title":"Phys. Chem. Chem. Phys."},{"key":"2023080302020032300_c6","doi-asserted-by":"publisher","first-page":"5777","DOI":"10.1063\/1.455553","article-title":"A method for two-electron Gaussian integral and integral derivative evaluation using recurrence relations","volume":"89","year":"1988","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c7","doi-asserted-by":"publisher","first-page":"5564","DOI":"10.1021\/j100377a031","article-title":"Efficient computation of two-electron-repulsion integrals and their nth-order derivatives using contracted Gaussian basis sets","volume":"94","year":"1990","journal-title":"J. Phys. Chem."},{"key":"2023080302020032300_c8","doi-asserted-by":"publisher","first-page":"1505","DOI":"10.1002\/1096-987x(200012)21:16&lt;1505::aid-jcc7&gt;3.0.co;2-4","article-title":"Rapid evaluation of two-center two-electron integrals","volume":"21","year":"2000","journal-title":"J. Comput. Chem."},{"key":"2023080302020032300_c9","doi-asserted-by":"publisher","first-page":"111","DOI":"10.1063\/1.432807","article-title":"Evaluation of molecular integrals over Gaussian basis functions","volume":"65","year":"1976","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c10","doi-asserted-by":"publisher","first-page":"154","DOI":"10.1002\/jcc.540040206","article-title":"Computation of electron repulsion integrals using the Rys quadrature method","volume":"4","year":"1983","journal-title":"J. Comput. Chem."},{"key":"2023080302020032300_c11","doi-asserted-by":"publisher","first-page":"1664","DOI":"10.1002\/jcc.23981","article-title":"Libcint: An efficient general integral library for Gaussian basis functions","volume":"36","year":"2015","journal-title":"J. Comput. Chem."},{"key":"2023080302020032300_c12","doi-asserted-by":"publisher","first-page":"5889","DOI":"10.1063\/1.461610","article-title":"The reduced multiplication scheme of the Rys quadrature and new recurrence relations for auxiliary function based two-electron integral evaluation","volume":"95","year":"1991","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c13","doi-asserted-by":"publisher","first-page":"2067","DOI":"10.1063\/1.1336541","article-title":"The Rys quadrature revisited: A novel formulation for the efficient computation of electron repulsion integrals over Gaussian functions","volume":"114","year":"2001","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c14","doi-asserted-by":"publisher","first-page":"334","DOI":"10.1002\/jcc.20779","article-title":"Two-electron integral evaluation on the graphics processor unit","volume":"29","year":"2008","journal-title":"J. Comput. Chem."},{"key":"2023080302020032300_c15","doi-asserted-by":"publisher","first-page":"696","DOI":"10.1021\/ct9005079","article-title":"Uncontracted Rys quadrature implementation of up to g functions on graphical processing units","volume":"6","year":"2010","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c16","doi-asserted-by":"publisher","first-page":"2387","DOI":"10.1021\/ct200239p","article-title":"Determining the numerical stability of quantum chemistry algorithms","volume":"7","year":"2011","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c17","doi-asserted-by":"publisher","first-page":"949","DOI":"10.1021\/ct100701w","article-title":"Dynamic precision for electron repulsion integral evaluation on graphical processing units (GPUs)","volume":"7","year":"2011","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c18","doi-asserted-by":"publisher","first-page":"4088","DOI":"10.1021\/acs.jctc.8b00321","article-title":"Double precision is not needed for many-body calculations: Emergent conventional wisdom","volume":"14","year":"2018","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c19","doi-asserted-by":"publisher","first-page":"213","DOI":"10.1021\/ct300321a","article-title":"Generating efficient quantum chemistry codes for novel architectures","volume":"9","year":"2012","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c20","doi-asserted-by":"publisher","first-page":"5319","DOI":"10.1021\/acs.jctc.9b00560","article-title":"Calculation of quantum chemical two-electron integrals by applying compiler technology on GPU","volume":"15","year":"2019","journal-title":"J. Chem. Theory Comput."},{"volume-title":"Molecular Electronic Structure Theory","year":"2000","key":"2023080302020032300_c21"},{"key":"2023080302020032300_c22","doi-asserted-by":"publisher","first-page":"2570","DOI":"10.1021\/acs.jctc.9b01296","article-title":"Fast evaluation of two-center integrals over Gaussian charge distributions and Gaussian orbitals with general interaction kernels","volume":"16","year":"2020","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c23","doi-asserted-by":"publisher","first-page":"269","DOI":"10.1002\/qua.560360831","article-title":"An efficient algorithm for the generation of two-electron repulsion integrals over Gaussian basis functions","volume":"36","year":"1989","journal-title":"Int. J. Quantum Chem."},{"key":"2023080302020032300_c24","doi-asserted-by":"publisher","first-page":"809","DOI":"10.1002\/qua.560400610","article-title":"Exact and approximate solutions to the one-center McMurchie\u2013Davidson tree-search problem","volume":"40","year":"1991","journal-title":"Int. J. Quantum Chem."},{"key":"2023080302020032300_c25","doi-asserted-by":"publisher","first-page":"572","DOI":"10.1021\/acs.jctc.7b00788","article-title":"LIBRETA: Computerized optimization and code synthesis for electron repulsion integral evaluation","volume":"14","year":"2018","journal-title":"J. Chem. Theory Comput."},{"edition":"2nd ed.","volume-title":"Circulant Matrices","year":"2012","key":"2023080302020032300_c27"},{"key":"2023080302020032300_c28","doi-asserted-by":"publisher","first-page":"216","DOI":"10.1109\/jproc.2004.840301","article-title":"The design and implementation of FFTW3","volume":"93","year":"2005","journal-title":"Proc. IEEE"},{"key":"2023080302020032300_c29","article-title":"Implementing FFTs in practice","volume-title":"Fast Fourier Transforms","author":"Burrus","year":"2008"},{"first-page":"285","article-title":"Cache-oblivious algorithms","year":"1999","key":"2023080302020032300_c30"},{"key":"2023080302020032300_c32","doi-asserted-by":"publisher","first-page":"e103","DOI":"10.7717\/peerj-cs.103","article-title":"SymPy: Symbolic computing in Python","volume":"3","year":"2017","journal-title":"PeerJ Comput. Sci."},{"first-page":"1","article-title":"Generation of kernels for calculating electron repulsion integrals of high angular momentum functions on GPUs\u2013preliminary results","year":"2010","key":"2023080302020032300_c33"},{"key":"2023080302020032300_c34","doi-asserted-by":"publisher","first-page":"1449","DOI":"10.1021\/ct500984t","article-title":"Acceleration of high angular momentum electron repulsion integrals and integral derivatives on graphics processing units","volume":"11","year":"2015","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c35","doi-asserted-by":"publisher","first-page":"3160","DOI":"10.1021\/acs.jctc.7b00030","article-title":"Arbitrary angular momentum electron repulsion integrals with graphical processing units: Application to the resolution of identity Hartree\u2013Fock method","volume":"13","year":"2017","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c36","doi-asserted-by":"publisher","first-page":"044103","DOI":"10.1063\/1.4732310","article-title":"Tensor hypercontraction density fitting. I. Quartic scaling second-and third-order M\u00f8ller\u2013Plesset perturbation theory","volume":"137","year":"2012","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c37","doi-asserted-by":"publisher","first-page":"4915","DOI":"10.1021\/acs.jctc.6b00770","article-title":"Two-electron integrals over Gaussian geminals","volume":"12","year":"2016","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c38","doi-asserted-by":"publisher","first-page":"024103","DOI":"10.1063\/1.4991733","article-title":"Three- and four-electron integrals involving Gaussian geminals: Fundamental integrals, upper bounds, and recurrence relations","volume":"147","year":"2017","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c39","doi-asserted-by":"publisher","first-page":"221106","DOI":"10.1063\/1.2817618","article-title":"A simple and efficient CCSD(T)-F12 approximation","volume":"127","year":"2007","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c40","doi-asserted-by":"publisher","first-page":"054104","DOI":"10.1063\/1.3054300","article-title":"Simplified CCSD(T)-F12 methods: Theory and benchmarks","volume":"130","year":"2009","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c41","doi-asserted-by":"publisher","first-page":"154103","DOI":"10.1063\/1.2889388","article-title":"Explicitly correlated RMP2 for high-spin open-shell reference states","volume":"128","year":"2008","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c42","doi-asserted-by":"publisher","first-page":"034113","DOI":"10.1063\/1.3528720","article-title":"Explicitly correlated multireference configuration interaction: MRCI-F12","volume":"134","year":"2011","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c43","doi-asserted-by":"publisher","first-page":"1735","DOI":"10.1021\/acs.jctc.6b00130","article-title":"Many-electron integrals over Gaussian basis functions. I. Recurrence relations for three-electron integrals","volume":"12","year":"2016","journal-title":"J. Chem. Theory Comput."},{"key":"2023080302020032300_c44","doi-asserted-by":"publisher","first-page":"6397","DOI":"10.1063\/1.1461814","article-title":"Explicitly correlated second-order M\u00f8ller\u2013Plesset methods with auxiliary basis sets","volume":"116","year":"2002","journal-title":"J. Chem. Phys."},{"key":"2023080302020032300_c45","doi-asserted-by":"publisher","first-page":"190","DOI":"10.1016\/j.cplett.2004.07.061","article-title":"Improving on the resolution of the identity in linear R12 ab initio theories","volume":"395","year":"2004","journal-title":"Chem. Phys. Lett."}],"container-title":["The Journal of Chemical Physics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/pubs.aip.org\/aip\/jcp\/article-pdf\/doi\/10.1063\/5.0002880\/15574757\/231103_1_online.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/pubs.aip.org\/aip\/jcp\/article-pdf\/doi\/10.1063\/5.0002880\/15574757\/231103_1_online.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/pubs.aip.org\/aip\/jcp\/article-pdf\/doi\/10.1063\/5.0002880\/15574757\/231103_1_online.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,5,27]],"date-time":"2025-05-27T20:20:29Z","timestamp":1748377229000},"score":1,"resource":{"primary":{"URL":"https:\/\/pubs.aip.org\/jcp\/article\/152\/23\/231103\/946316\/Molecular-integrals-from-Fast-Fourier-Transforms"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2020,6,17]]},"references-count":43,"journal-issue":{"issue":"23","published-print":{"date-parts":[[2020,6,21]]}},"URL":"https:\/\/doi.org\/10.1063\/5.0002880","relation":{},"ISSN":["0021-9606","1089-7690"],"issn-type":[{"type":"print","value":"0021-9606"},{"type":"electronic","value":"1089-7690"}],"subject":[],"published-other":{"date-parts":[[2020,6,21]]},"published":{"date-parts":[[2020,6,17]]},"article-number":"231103"}}