A differential derivation of the Obara-Saika relation for Gaussian electron repulsion integrals

Charles C. Forgy; David A. Mazziotti

arxiv: 2607.00102 · v1 · pith:5S5Z4X7Jnew · submitted 2026-06-30 · ⚛️ physics.chem-ph · quant-ph

A differential derivation of the Obara-Saika relation for Gaussian electron repulsion integrals

Charles C. Forgy , David A. Mazziotti This is my paper

Pith reviewed 2026-07-02 17:09 UTC · model grok-4.3

classification ⚛️ physics.chem-ph quant-ph

keywords Obara-Saika relationelectron repulsion integralsGaussian basis functionsvertical recurrencedifferential derivationprimitive termsquantum chemistryrecurrence relations

0 comments

The pith

A derivation of the Obara-Saika recurrence uses only differential relations between Gaussian basis functions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes a compact derivation of the Obara-Saika vertical recurrence relation for electron repulsion integrals that starts from differential relations among Gaussian basis functions rather than from direct relations between integrals of different angular momentum. By working out the explicit derivative expressions and cataloging every non-zero primitive term, the authors obtain the standard recurrence in hierarchical form. A sympathetic reader would care because the Obara-Saika method is a core computational step in quantum chemistry, and a differential route supplies an independent, self-contained route to the same relations while exposing the structure of the contributing terms.

Core claim

The central claim is that the Obara-Saika vertical recurrence follows directly once the differential action of angular-momentum operators on Gaussian products is written out and every non-zero primitive contribution to the full electron repulsion integral is identified, yielding the familiar recurrence relations without additional integral identities.

What carries the argument

Differential relations between Gaussian basis functions, which generate the vertical recurrence by isolating non-zero primitive derivative terms.

If this is right

The recursion relations admit a hierarchical formulation ordered by the independent primitive derivative quantities.
The derivation stands as a self-contained alternative to the original Obara-Saika argument.
The separation into independent primitive quantities supports modular code generation.
The same organization lends itself to parallel evaluation on GPU hardware.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The explicit primitive terms could be precomputed once and reused across multiple integral batches.
A similar differential starting point might be applied to other families of molecular integrals that obey recurrence relations.
The hierarchical view could guide automatic differentiation or symbolic simplification tools for integral code.

Load-bearing premise

Explicit derivative expressions alone suffice to capture every non-zero primitive term in the electron repulsion integral.

What would settle it

Compute a low-angular-momentum electron repulsion integral both by the derived recurrence and by direct six-dimensional numerical quadrature of the Gaussian product; any numerical mismatch falsifies the claim that the derivatives recover the full set of terms.

read the original abstract

The Obara-Saika (OS) method is one of the most widely used techniques in quantum chemistry for evaluating electron repulsion integrals (ERIs) via a set of recurrence relations that build higher angular momentum integrals from lower-order ones. The original derivation by Obara and Saika proceeded by directly relating integrals of differing angular momentum. In this work, we present a compact novel derivation of the OS vertical recurrence relation based solely on differential relations between Gaussian basis functions, expanding on a method suggested in earlier work. By explicitly deriving the required derivative expressions we identify all non-zero primitive terms contributing to the full ERI to develop a hierarchical formulation of the OS recursion relations. This approach has pedagogical value as a rigorous and self-contained derivation. Additionally, the resulting organization exposes independent primitive derivative quantities and may be useful for code generation and parallel implementations on modern GPU architectures.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

A re-derivation of the Obara-Saika recurrence via Gaussian derivatives that checks out algebraically but adds little new capability.

read the letter

The main thing here is a derivation of the OS vertical recurrence that starts from differential identities on the Gaussian primitives rather than from direct angular-momentum relations. The authors apply the product rule and the known Gaussian derivative formulas repeatedly, collect the non-zero terms, and arrive at the familiar recurrence without extra pruning steps. The stress-test confirms the algebra terminates cleanly at the lower-momentum integrals and that no terms are missing or added by hand.

What the paper actually supplies is therefore an alternative route to an established result, plus a hierarchical organization of the primitive derivative quantities. That organization is presented as potentially helpful for code generation and GPU parallelism, which is plausible on its face. The pedagogical framing is also reasonable; anyone who has taught the OS scheme might find the differential path cleaner for some students.

The limitation is straightforward: the recurrence itself is not new, and the paper does not report timing comparisons, memory savings, or new molecular calculations that would show a practical edge. The abstract and the stress-test both indicate the derivation is self-contained, so there is no obvious circularity or hidden fitting. Still, the contribution sits at the level of a methods note rather than a result that changes what codes can compute.

This is the sort of paper that belongs in a quantum-chemistry methods journal if the derivation is fully spelled out and the hierarchical structure is shown to be useful for implementation. It is worth a referee's time because the algebra is reproducible and the claim is modest. I would not bring it to a general reading group, and I would not cite it unless I were writing about ERI code organization.

Referee Report

0 major / 2 minor

Summary. The manuscript presents a compact derivation of the Obara-Saika vertical recurrence relation for Gaussian electron repulsion integrals that relies exclusively on repeated application of the product rule and standard Gaussian differentiation identities. It explicitly constructs the required derivative expressions, enumerates all non-zero primitive terms contributing to the ERI, and organizes the resulting relations into a hierarchical formulation that terminates at lower-angular-momentum integrals without external recurrence relations or case-by-case pruning.

Significance. If the derivation holds, the work supplies a self-contained algebraic route to the OS relations that is pedagogically clearer than the original angular-momentum-based proof and exposes independent primitive derivative quantities. This organization is noted as potentially useful for automated code generation and GPU-parallel implementations. The approach is parameter-free and reproduces the known recurrence structure without hidden dependencies or post-hoc selection rules.

minor comments (2)

The introduction states that the method expands on 'earlier work' but does not give an explicit citation in the first paragraph; adding the reference at that point would improve traceability.
In the hierarchical formulation (around the discussion of primitive terms), the notation for the independent derivative quantities could be introduced with a small table or explicit list to make the independence immediately visible to readers implementing the scheme.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript and the recommendation to accept. The comments correctly identify the self-contained differential derivation and its potential utility.

Circularity Check

0 steps flagged

Derivation is algebraically self-contained; no circular reductions identified

full rationale

The manuscript derives the Obara-Saika vertical recurrence solely by repeated application of the product rule and standard Gaussian differentiation identities (e.g., ∂/∂A_x (exp(-α|r-A|^2)) = -2α(x-A_x)exp(...)) directly to the ERI kernel. Each algebraic step terminates at lower-angular-momentum primitive integrals without introducing fitted parameters, selection rules, or external recurrence relations. The abstract's reference to 'expanding on a method suggested in earlier work' is non-load-bearing; the central chain relies only on the explicit derivative expressions and the product rule, which are independent of any prior result by the same authors. No self-definitional loops, fitted-input predictions, or ansatz smuggling appear in the hierarchical organization of non-zero primitive terms.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities can be extracted or verified.

pith-pipeline@v0.9.1-grok · 5677 in / 1064 out tokens · 25462 ms · 2026-07-02T17:09:55.594610+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

35 extracted references · 22 canonical work pages

[1]

Obara, S., Saika, A.: Efficient recursive computation of molecular integrals over Cartesian Gaussian functions. J. Chem. Phys.84(7), 3963–3974 (1986)

1986
[2]

McMurchie, L.E., Davidson, E.R.: One- and two-electron integrals over cartesian Gaussian functions. J. Comput. Phys.26(2), 218–231 (1978) https://doi.org/10. 1016/0021-9991(78)90092-X

1978
[3]

Proceedings of the Royal Society of London A: Mathe- matical, Physical and Engineering Sciences258(1294), 402–411 (1960) https: //doi.org/10.1098/rspa.1960.0195

Boys, S.F.: The integral formulae for the variational solution of the molecu- lar many-electron wave equations in terms of Gaussian functions with direct electronic correlation. Proceedings of the Royal Society of London A: Mathe- matical, Physical and Engineering Sciences258(1294), 402–411 (1960) https: //doi.org/10.1098/rspa.1960.0195

work page doi:10.1098/rspa.1960.0195 1960
[4]

PhD thesis, Budapest University of Technology and Economics (2018)

Samu, G.: Effcient evaluation of three-center coulomb-integrals and their geo- metrical first derivatives. PhD thesis, Budapest University of Technology and Economics (2018)

2018
[5]

WIREs Comput

Lehtola, S., Karttunen, A.J.: Free and open source software for computational chemistry education. WIREs Comput. Mol. Sci.12(5), 1610 (2022) https://doi. org/10.1002/wcms.1610

work page doi:10.1002/wcms.1610 2022
[6]

Asadchev, A., Valeev, E.F.: Memory-efficient recursive evaluation of 3-center Gaussian integrals. J. Chem. Theory Comput.19(6), 1698–1710 (2023) https: //doi.org/10.1021/acs.jctc.2c00995

work page doi:10.1021/acs.jctc.2c00995 2023
[7]

Asadchev, A., Valeev, E.F.: High-performance evaluation of high angular momen- tum 4-center Gaussian integrals on modern accelerated processors. J. Phys. Chem. A127(51), 10889–10895 (2023) https://doi.org/10.1021/acs.jpca.3c04574

work page doi:10.1021/acs.jpca.3c04574 2023
[8]

Asadchev, A., Valeev, E.F.: 3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors. J. Chem. Phys.160(24), 244109 (2024) https: //doi.org/10.1063/5.0217001

work page doi:10.1063/5.0217001 2024
[9]

Kim, T.D., Pujal, L., Richer, M., Zyl, M., Mart´ ınez-Gonz´ alez, M., Tehrani, A., Chuiko, V., S´ anchez-D´ ıaz, G., Sanchez, W., Adams, W., Huang, X., Kelly, B.D., V¨ ohringer-Martinez, E., Verstraelen, T., Heidar-Zadeh, F., Ayers, P.W.: GBasis: A Python library for evaluating functions, functionals, and integrals expressed with Gaussian basis functions....

work page doi:10.1063/5.0216776 2024
[10]

Pritchard, B.P., Chow, E.: Horizontal vectorization of electron repulsion integrals. J. Comput. Chem.37(28), 2537–2546 (2016) https://doi.org/10.1002/jcc.24483

work page doi:10.1002/jcc.24483 2016
[12]

Gordon, M.S., Barca, G., Leang, S.S., Poole, D., Rendell, A.P., Galvez Vallejo, J.L., Westheimer, B.: Novel computer architectures and quantum chemistry. J. Phys. Chem. A124(23), 4557–4582 (2020) https://doi.org/10.1021/acs.jpca. 0c02249

work page doi:10.1021/acs.jpca 2020
[13]

Li, R., Sun, Q., Zhang, X., Chan, G.K.-L.: Introducing GPU acceleration into the Python-based simulations of chemistry framework. J. Chem. Phys. A129(5), 1459–1468 (2025) https://doi.org/10.1021/acs.jpca.4c05876

work page doi:10.1021/acs.jpca.4c05876 2025
[14]

WIREs Comput

Sun, Q., Berkelbach, T.C., Blunt, N.S., Booth, G.H., Guo, S., Li, Z., Liu, J., McClain, J.D., Sayfutyarova, E.R., Sharma, S., Wouters, S., Chan, G.K.-L.: PySCF: the Python-based simulations of chemistry framework. WIREs Comput. Mol. Sci.8(1), 1340 (2018) https://doi.org/10.1002/wcms.1340

work page doi:10.1002/wcms.1340 2018
[15]

Dupuis, M., Rys, J., King, H.F.: Evaluation of molecular integrals over Gaussian basis functions. J. Chem. Phys.65(1), 111–116 (1976) https://doi.org/10.1063/ 1.432807

1976
[16]

Palethorpe, E., Barca, G.M.J.: High-performance, high-angular-momentum J engine on graphics processing units. J. Chem. Theory Comput.21(19), 9388–9403 (2025) https://doi.org/10.1021/acs.jctc.5c00775

work page doi:10.1021/acs.jctc.5c00775 2025
[17]

Fujii, H., Ito, Y., Yokogawa, N., Suzuki, K., Tsuji, S., Nakano, K., Parque, V., Kasagi, A.: Efficient GPU implementation of the McMurchie–Davidson method for shell-based ERI computations. Appl. Sci.15(5) (2025)

2025
[18]

Ahlrichs, R.: A simple algebraic derivation of the Obara–Saika scheme for general two-electron interaction potentials. Phys. Chem. Chem. Phys.8(26), 3072–3077 (2006)

2006
[19]

Limpanuparb, T., Hollett, J.W., Gill, P.M.W.: Resolutions of the Coulomb oper- ator. VI. computation of auxiliary integrals. J. Chem. Phys.136(10), 104102 (2012) https://doi.org/10.1063/1.3691829

work page doi:10.1063/1.3691829 2012
[20]

Head-Gordon, M., Pople, J.A.: A method for two-electron Gaussian integral and integral derivative evaluation using recurrence relations. J. Chem. Phys.89(9), 5777–5786 (1988) https://doi.org/10.1063/1.455553

work page doi:10.1063/1.455553 1988
[21]

Pulay, P.: Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules. Mol. Phys.17(2), 197–204 (1969)

1969
[22]

Gerratt, J., Mills, I.M.: Force constants and dipole-moment derivatives of molecules from perturbed Hartree–Fock calculations. I. J. Chem. Phys.49(4), 1719–1729 (1968) https://doi.org/10.1063/1.1670299 19

work page doi:10.1063/1.1670299 1968
[23]

Park, J.W., Al-Saadon, R., MacLeod, M.K., Shiozaki, T., Vlaisavljevich, B.: Multireference electron correlation methods: Journeys along potential energy surfaces. Chem. Rev.120(13), 5878–5909 (2020) https://doi.org/10.1021/acs. chemrev.9b00496

work page doi:10.1021/acs 2020
[24]

Park, J.W.: Analytical nuclear gradient and derivative coupling theories for multireference perturbation methods. Phys. Chem. Chem. Phys.27, 3531–3551 (2025) https://doi.org/10.1039/D4CP03671A

work page doi:10.1039/d4cp03671a 2025
[25]

Hennefarth, M.R., Kim, Y., Jangid, B., Wardzala, J., Hermes, M.R., Truhlar, D.G., Gagliardi, L.: MC-PDFT nuclear gradients and L-PDFT energies with meta and hybrid meta on-top functionals for ground- and excited-state geometry optimization and vertical excitation energies. J. Chem. Theory Comput.21(16), 7890–7902 (2025) https://doi.org/10.1021/acs.jctc.5c00899

work page doi:10.1021/acs.jctc.5c00899 2025
[26]

Dover Publications, Inc., Mineola, New York (1989)

Szabo, A., Ostlund, N.S.: Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, 1st edn. Dover Publications, Inc., Mineola, New York (1989)

1989
[27]

Pearson, London, England (2008)

Levine, I.N.: Quantum Chemistry, 6th edn. Pearson, London, England (2008)

2008
[28]

Schlegel, H.B., Frisch, M.J.: Transformation between Cartesian and pure spherical harmonic Gaussians. Int. J. Quantum Chem.54(2), 83–87 (1995) https://doi. org/10.1002/qua.560540202

work page doi:10.1002/qua.560540202 1995
[29]

Reine, S., Tellgren, E., Helgaker, T.: A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaus- sians. Phys. Chem. Chem. Phys.9, 4771–4779 (2007) https://doi.org/10.1039/ B705594C

2007
[30]

Advances in Quantum Chemistry, vol

Gill, P.M.W.: Molecular integrals over Gaussian basis functions. Advances in Quantum Chemistry, vol. 25, pp. 141–205. Academic Press (1994). https://doi. org/10.1016/S0065-3276(08)60019-2 . http://www.sciencedirect.com/science/ article/pii/S0065327608600192

work page doi:10.1016/s0065-3276(08)60019-2 1994
[31]

Wiley, West Sussex, England (2000)

Helgaker, T., Jørgensen, P., Olsen, J.: Molecular Electronic-Structure Theory, 1st edn. Wiley, West Sussex, England (2000)

2000
[32]

Boys, S.F.: Electronic wave functions. I. A general method of calculation for the stationary states of any molecular system. Proc. Roy. Soc. London Ser. A 200(1063), 542–554 (1950) https://doi.org/10.1098/rspa.1950.0036

work page doi:10.1098/rspa.1950.0036 1950
[33]

Electron

Hardy, M.: Combinatorics of partial derivatives. Electron. J. Combinatorics13(1), 1–13113 (2006)

2006
[34]

Schlegel, H.B.: An efficient algorithm for calculatingab initioenergy gradients usings,pCartesian Gaussians. J. Chem. Phys.77(7), 3676–3681 (1982) https: 20 //doi.org/10.1063/1.444270

work page doi:10.1063/1.444270 1982
[35]

Ryu, U., Lee, Y.S., Lindh, R.: An efficient method of implementing the horizontal recurrence relation in the evaluation of electron repulsion integrals using cartesian Gaussian functios. Chem. Phys. Lett.185(5), 562–568 (1991) https://doi.org/10. 1016/0009-2614(91)80260-5

1991
[36]

Chem, Phys.150(2), 163–171 (1991) https://doi.org/10.1016/0301-0104(91)80126-3 21

Hamilton, T.P., Schaefer, H.F.: New variations in two-electron integral evaluation in the context of direct scf procedures. Chem, Phys.150(2), 163–171 (1991) https://doi.org/10.1016/0301-0104(91)80126-3 21

work page doi:10.1016/0301-0104(91)80126-3 1991

[1] [1]

Obara, S., Saika, A.: Efficient recursive computation of molecular integrals over Cartesian Gaussian functions. J. Chem. Phys.84(7), 3963–3974 (1986)

1986

[2] [2]

McMurchie, L.E., Davidson, E.R.: One- and two-electron integrals over cartesian Gaussian functions. J. Comput. Phys.26(2), 218–231 (1978) https://doi.org/10. 1016/0021-9991(78)90092-X

1978

[3] [3]

Proceedings of the Royal Society of London A: Mathe- matical, Physical and Engineering Sciences258(1294), 402–411 (1960) https: //doi.org/10.1098/rspa.1960.0195

Boys, S.F.: The integral formulae for the variational solution of the molecu- lar many-electron wave equations in terms of Gaussian functions with direct electronic correlation. Proceedings of the Royal Society of London A: Mathe- matical, Physical and Engineering Sciences258(1294), 402–411 (1960) https: //doi.org/10.1098/rspa.1960.0195

work page doi:10.1098/rspa.1960.0195 1960

[4] [4]

PhD thesis, Budapest University of Technology and Economics (2018)

Samu, G.: Effcient evaluation of three-center coulomb-integrals and their geo- metrical first derivatives. PhD thesis, Budapest University of Technology and Economics (2018)

2018

[5] [5]

WIREs Comput

Lehtola, S., Karttunen, A.J.: Free and open source software for computational chemistry education. WIREs Comput. Mol. Sci.12(5), 1610 (2022) https://doi. org/10.1002/wcms.1610

work page doi:10.1002/wcms.1610 2022

[6] [6]

Asadchev, A., Valeev, E.F.: Memory-efficient recursive evaluation of 3-center Gaussian integrals. J. Chem. Theory Comput.19(6), 1698–1710 (2023) https: //doi.org/10.1021/acs.jctc.2c00995

work page doi:10.1021/acs.jctc.2c00995 2023

[7] [7]

Asadchev, A., Valeev, E.F.: High-performance evaluation of high angular momen- tum 4-center Gaussian integrals on modern accelerated processors. J. Phys. Chem. A127(51), 10889–10895 (2023) https://doi.org/10.1021/acs.jpca.3c04574

work page doi:10.1021/acs.jpca.3c04574 2023

[8] [8]

Asadchev, A., Valeev, E.F.: 3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors. J. Chem. Phys.160(24), 244109 (2024) https: //doi.org/10.1063/5.0217001

work page doi:10.1063/5.0217001 2024

[9] [9]

Kim, T.D., Pujal, L., Richer, M., Zyl, M., Mart´ ınez-Gonz´ alez, M., Tehrani, A., Chuiko, V., S´ anchez-D´ ıaz, G., Sanchez, W., Adams, W., Huang, X., Kelly, B.D., V¨ ohringer-Martinez, E., Verstraelen, T., Heidar-Zadeh, F., Ayers, P.W.: GBasis: A Python library for evaluating functions, functionals, and integrals expressed with Gaussian basis functions....

work page doi:10.1063/5.0216776 2024

[10] [10]

Pritchard, B.P., Chow, E.: Horizontal vectorization of electron repulsion integrals. J. Comput. Chem.37(28), 2537–2546 (2016) https://doi.org/10.1002/jcc.24483

work page doi:10.1002/jcc.24483 2016

[11] [12]

Gordon, M.S., Barca, G., Leang, S.S., Poole, D., Rendell, A.P., Galvez Vallejo, J.L., Westheimer, B.: Novel computer architectures and quantum chemistry. J. Phys. Chem. A124(23), 4557–4582 (2020) https://doi.org/10.1021/acs.jpca. 0c02249

work page doi:10.1021/acs.jpca 2020

[12] [13]

Li, R., Sun, Q., Zhang, X., Chan, G.K.-L.: Introducing GPU acceleration into the Python-based simulations of chemistry framework. J. Chem. Phys. A129(5), 1459–1468 (2025) https://doi.org/10.1021/acs.jpca.4c05876

work page doi:10.1021/acs.jpca.4c05876 2025

[13] [14]

WIREs Comput

Sun, Q., Berkelbach, T.C., Blunt, N.S., Booth, G.H., Guo, S., Li, Z., Liu, J., McClain, J.D., Sayfutyarova, E.R., Sharma, S., Wouters, S., Chan, G.K.-L.: PySCF: the Python-based simulations of chemistry framework. WIREs Comput. Mol. Sci.8(1), 1340 (2018) https://doi.org/10.1002/wcms.1340

work page doi:10.1002/wcms.1340 2018

[14] [15]

Dupuis, M., Rys, J., King, H.F.: Evaluation of molecular integrals over Gaussian basis functions. J. Chem. Phys.65(1), 111–116 (1976) https://doi.org/10.1063/ 1.432807

1976

[15] [16]

Palethorpe, E., Barca, G.M.J.: High-performance, high-angular-momentum J engine on graphics processing units. J. Chem. Theory Comput.21(19), 9388–9403 (2025) https://doi.org/10.1021/acs.jctc.5c00775

work page doi:10.1021/acs.jctc.5c00775 2025

[16] [17]

Fujii, H., Ito, Y., Yokogawa, N., Suzuki, K., Tsuji, S., Nakano, K., Parque, V., Kasagi, A.: Efficient GPU implementation of the McMurchie–Davidson method for shell-based ERI computations. Appl. Sci.15(5) (2025)

2025

[17] [18]

Ahlrichs, R.: A simple algebraic derivation of the Obara–Saika scheme for general two-electron interaction potentials. Phys. Chem. Chem. Phys.8(26), 3072–3077 (2006)

2006

[18] [19]

Limpanuparb, T., Hollett, J.W., Gill, P.M.W.: Resolutions of the Coulomb oper- ator. VI. computation of auxiliary integrals. J. Chem. Phys.136(10), 104102 (2012) https://doi.org/10.1063/1.3691829

work page doi:10.1063/1.3691829 2012

[19] [20]

Head-Gordon, M., Pople, J.A.: A method for two-electron Gaussian integral and integral derivative evaluation using recurrence relations. J. Chem. Phys.89(9), 5777–5786 (1988) https://doi.org/10.1063/1.455553

work page doi:10.1063/1.455553 1988

[20] [21]

Pulay, P.: Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules. Mol. Phys.17(2), 197–204 (1969)

1969

[21] [22]

Gerratt, J., Mills, I.M.: Force constants and dipole-moment derivatives of molecules from perturbed Hartree–Fock calculations. I. J. Chem. Phys.49(4), 1719–1729 (1968) https://doi.org/10.1063/1.1670299 19

work page doi:10.1063/1.1670299 1968

[22] [23]

Park, J.W., Al-Saadon, R., MacLeod, M.K., Shiozaki, T., Vlaisavljevich, B.: Multireference electron correlation methods: Journeys along potential energy surfaces. Chem. Rev.120(13), 5878–5909 (2020) https://doi.org/10.1021/acs. chemrev.9b00496

work page doi:10.1021/acs 2020

[23] [24]

Park, J.W.: Analytical nuclear gradient and derivative coupling theories for multireference perturbation methods. Phys. Chem. Chem. Phys.27, 3531–3551 (2025) https://doi.org/10.1039/D4CP03671A

work page doi:10.1039/d4cp03671a 2025

[24] [25]

Hennefarth, M.R., Kim, Y., Jangid, B., Wardzala, J., Hermes, M.R., Truhlar, D.G., Gagliardi, L.: MC-PDFT nuclear gradients and L-PDFT energies with meta and hybrid meta on-top functionals for ground- and excited-state geometry optimization and vertical excitation energies. J. Chem. Theory Comput.21(16), 7890–7902 (2025) https://doi.org/10.1021/acs.jctc.5c00899

work page doi:10.1021/acs.jctc.5c00899 2025

[25] [26]

Dover Publications, Inc., Mineola, New York (1989)

Szabo, A., Ostlund, N.S.: Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, 1st edn. Dover Publications, Inc., Mineola, New York (1989)

1989

[26] [27]

Pearson, London, England (2008)

Levine, I.N.: Quantum Chemistry, 6th edn. Pearson, London, England (2008)

2008

[27] [28]

Schlegel, H.B., Frisch, M.J.: Transformation between Cartesian and pure spherical harmonic Gaussians. Int. J. Quantum Chem.54(2), 83–87 (1995) https://doi. org/10.1002/qua.560540202

work page doi:10.1002/qua.560540202 1995

[28] [29]

Reine, S., Tellgren, E., Helgaker, T.: A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaus- sians. Phys. Chem. Chem. Phys.9, 4771–4779 (2007) https://doi.org/10.1039/ B705594C

2007

[29] [30]

Advances in Quantum Chemistry, vol

Gill, P.M.W.: Molecular integrals over Gaussian basis functions. Advances in Quantum Chemistry, vol. 25, pp. 141–205. Academic Press (1994). https://doi. org/10.1016/S0065-3276(08)60019-2 . http://www.sciencedirect.com/science/ article/pii/S0065327608600192

work page doi:10.1016/s0065-3276(08)60019-2 1994

[30] [31]

Wiley, West Sussex, England (2000)

Helgaker, T., Jørgensen, P., Olsen, J.: Molecular Electronic-Structure Theory, 1st edn. Wiley, West Sussex, England (2000)

2000

[31] [32]

Boys, S.F.: Electronic wave functions. I. A general method of calculation for the stationary states of any molecular system. Proc. Roy. Soc. London Ser. A 200(1063), 542–554 (1950) https://doi.org/10.1098/rspa.1950.0036

work page doi:10.1098/rspa.1950.0036 1950

[32] [33]

Electron

Hardy, M.: Combinatorics of partial derivatives. Electron. J. Combinatorics13(1), 1–13113 (2006)

2006

[33] [34]

Schlegel, H.B.: An efficient algorithm for calculatingab initioenergy gradients usings,pCartesian Gaussians. J. Chem. Phys.77(7), 3676–3681 (1982) https: 20 //doi.org/10.1063/1.444270

work page doi:10.1063/1.444270 1982

[34] [35]

Ryu, U., Lee, Y.S., Lindh, R.: An efficient method of implementing the horizontal recurrence relation in the evaluation of electron repulsion integrals using cartesian Gaussian functios. Chem. Phys. Lett.185(5), 562–568 (1991) https://doi.org/10. 1016/0009-2614(91)80260-5

1991

[35] [36]

Chem, Phys.150(2), 163–171 (1991) https://doi.org/10.1016/0301-0104(91)80126-3 21

Hamilton, T.P., Schaefer, H.F.: New variations in two-electron integral evaluation in the context of direct scf procedures. Chem, Phys.150(2), 163–171 (1991) https://doi.org/10.1016/0301-0104(91)80126-3 21

work page doi:10.1016/0301-0104(91)80126-3 1991