Chemistry - Gradient of integrals with Obara Saika recurrence scheme

Solution 1:

I see you have not gotten any answers for some time. I have not myself become familiar enough with the notation used in S. Obara and A. Saika, J. Chem. Phys. 84 (7), 1986. But I have a suggestion for how you can help your self proceed. If you use your derived formula to calculate some derivative, and you try to calculate the same derivatives with finite difference, then you can check and see if you get the same result. If the results are identical to some nummerical error, your equation are proberly correct.

UPDATE:

I have tried to make a derivation of the derivative myself. Here is what I have come up with: The overlap integral is given as:

$$S_{\mu,\nu}=\left(\mu|\nu\right)$$

Here $\mu$ and $\nu$ represent Gaussian functions of the type:

$$\phi_{l,m,n}=N_{l,m,n}x_{a}^{l}y_{a}^{m}z_{a}^{n}\exp\left(-\alpha\left|r-R_{A}\right|^{2}\right)$$

Now the derivative with repsect to the movement of nucleus $A$ can be calculated by using the chain rule:

$$\left.\frac{\partial S_{\mu,\nu}}{\partial X_{A}}\right|_{X_{\mu}\neq X_{\nu}}=\left.\left(\left.\frac{\partial\mu}{\partial X_{A}}\right|\nu\right)\right|_{X_{\mu}=X_{A}}+\left.\left(\mu\left|\frac{\partial\nu}{\partial X_{A}}\right.\right)\right|_{X_{\nu}=X_{A}}$$

Note that $X_{\mu}\neq X_{\nu}$, this is imposed by the invariance of moving all centers [1]:

$$\frac{\partial S_{\mu,\nu}}{\partial X_{A}}+\frac{\partial S_{\mu,\nu}}{\partial X_{B}}=0$$

Here $A=B$. Also a wavefunction must have a depends on $A$ for a derivative to give a value, thus for atom centered wavefunctions, the wavefunction must be placed on the atom for which the derivative is taken. Now the the derivative of the Gaussian functions is given as [2]:

$$\frac{\partial\phi_{l,m,n}}{\partial X_{A}}=\sqrt{(2l+1)\alpha}\phi_{l+1,m,n}-\left.2l\sqrt{\left(\frac{\alpha}{2l-1}\right)}\phi_{l-1,m,n}\right|_{l>0}$$

It should be noted that the factors in front of the Gaussians is because of normalization. Thus now giving:

$$\left.\frac{\partial S_{\mu,\nu}}{\partial X_{A}}\right|_{X_{A}=X_{\mu}\neq X_{\nu}}=\left(\left.\frac{\partial\mu_{l,m,n}}{\partial X_{A}}\right|\nu_{l,m,n}\right)$$

$$=\left(\left.\sqrt{(2l+1)\alpha}\mu_{l+1,m,n}-\left.2l\sqrt{\left(\frac{\alpha}{2l-1}\right)}\mu_{l-1,m,n}\right|_{l>0}\right|\nu_{l,m,n}\right)$$

$$=\sqrt{(2l+1)\alpha}\left(\left.\mu_{l+1,m,n}\right|\nu_{l,m,n}\right)-2l\sqrt{\left(\frac{\alpha}{2l-1}\right)}\left(\left.\left.\mu_{l-1,m,n}\right|_{l>0}\right|\nu_{l,m,n}\right)$$

Now $S_{i,,j}=\left(\mu_{l=i,m,n}|\nu_{l=j,m,n}\right)$and by insertion into the found expression for the derivative:

$$\left.\frac{\partial S_{\mu,\nu}}{\partial X_{A}}\right|_{X_{A}=X_{\mu}\neq X_{\nu}}=\sqrt{(2l+1)\alpha}S_{i+1,j}-\left.2l\sqrt{\left(\frac{\alpha}{2l-1}\right)}S_{i-1.j}\right|_{i>0}$$

Now from here on I was unable to find the equation you found. Also my notation is abit different. I thought you might want it as an alternative way to find the derivative (or maybe you will be able see from this if your equation is right or wrong). Again if you do the calculations with numbers you can see if the two different equations give the same results. Now the Obara-Saika recurrence relation can be considered [3]:

$$S_{i+1,j}=X_{PA}S_{i,j}+\frac{1}{2p}(iS_{i-1,j}+jS_{i,j-1})$$

This can be insterted into the above derivative to eliminate the $i+1$ term:

$$\left.\frac{\partial S_{\mu,\nu}}{\partial X_{A}}\right|_{X_{A}=X_{\mu}\neq X_{\nu}}=\sqrt{(2l+1)\alpha}\left( X_{PA}S_{i,j}+\frac{1}{2p}(iS_{i-1,j}+jS_{i,j-1}) \right)-\left.2l\sqrt{\left(\frac{\alpha}{2l-1}\right)}S_{i-1.j}\right|_{i>0}$$

[1] Trygve Helgaker, Poul Jorgensen, Jeppe Olsen; Molecular Electronic-Structure Theory. Equation: 9.3.4

[2] Attila Szabo, Neil S. Ostlund;Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. page. 442

[3] Trygve Helgaker, Poul Jorgensen, Jeppe Olsen; Molecular Electronic-Structure Theory. Equation: 9.3.8

Solution 2:

I moved my extended explanation from my question to a real answer. The given problem was to calculate the gradient of overlap integrals. Thanks to Erik Kjellgren I noticed that I am looking at the false problem.

So far I have calculated derivatives for cartesian coordinates. This is possible but not very helpful in case I want to have gradient for geometry optimizations. I can use a bit from my previous work with

$$ \frac{\partial}{\partial R_i}\varphi ({\bf r}; \zeta, {\bf n}, {\bf R}) = - \frac{\partial}{\partial r_i}\varphi ({\bf r}; \zeta, {\bf n}, {\bf R}) $$

and get

$$ \frac{\partial}{\partial A_i} (a || b) = 2\zeta_a ( a + {\bf 1}_i || b ) - \sum^{l_a}_{j=1} \delta_{ij}( a - {\bf 1}_j || b ). $$

This is the exactly same as in Erik Kjellgren’s answer. If I use the Obara–Saika recurrsion formula on the expression above I get the following formula.

$$ \frac{\partial}{\partial A_i} (a || b) = 2\zeta_a \Bigl( (P_i-A_i)(a||b) + \sum^{l_a}_{j=1} \delta_{ij}\frac 1{2\zeta}( a - {\bf 1}_j || b ) + \sum^{l_b}_{k=1} \delta_{ik}\frac 1{2\zeta}( a || b - {\bf 1}_k ) \Bigr) - \sum^{l_a}_{j=1} \delta_{ij}( a - {\bf 1}_j || b ) $$

From here on I cannot go any further. The equation I have derived now does not look easy anymore, but I think it is fine, all intermediates are available from the calculation of the overlap integral itself, so I don’t need any additional expensive calculations.