# qml.qchem.gaussian_overlap¶

gaussian_overlap(la, lb, ra, rb, alpha, beta)[source]

Compute overlap integral for two primitive Gaussian functions.

The overlap integral between two Gaussian functions denoted by $$a$$ and $$b$$ can be computed as [Helgaker (1995) p803]:

$S_{ab} = E^{ij} E^{kl} E^{mn} \left (\frac{\pi}{p} \right )^{3/2},$

where $$E$$ is a coefficient that can be computed recursively, $$i-n$$ are the angular momentum quantum numbers corresponding to different Cartesian components and $$p$$ is computed from the exponents of the two Gaussian functions as $$p = \alpha + \beta$$.

Parameters
• la (integer) – angular momentum for the first Gaussian function

• lb (integer) – angular momentum for the second Gaussian function

• ra (float) – position vector of the first Gaussian function

• rb (float) – position vector of the second Gaussian function

• alpha (array[float]) – exponent of the first Gaussian function

• beta (array[float]) – exponent of the second Gaussian function

Returns

overlap integral between primitive Gaussian functions

Return type

array[float]

Example

>>> la, lb = (0, 0, 0), (0, 0, 0)
>>> ra, rb = np.array([0., 0., 0.]), np.array([0., 0., 0.])
>>> alpha = np.array([np.pi/2])
>>> beta = np.array([np.pi/2])
>>> o = gaussian_overlap(la, lb, ra, rb, alpha, beta)
>>> o
array([1.])


Using PennyLane

Development

API