Chemistry:STO-nG basis sets

From HandWiki
Short description: Basis sets used in quantum chemistry


STO-nG basis sets are minimal basis sets, where [math]\displaystyle{ n }[/math] primitive Gaussian orbitals are fitted to a single Slater-type orbital (STO). [math]\displaystyle{ n }[/math] originally took the values 2 – 6. They were first proposed by John Pople. A minimum basis set is where only sufficient orbitals are used to contain all the electrons in the neutral atom. Thus for the hydrogen atom, only a single 1s orbital is needed, while for a carbon atom, 1s, 2s and three 2p orbitals are needed. The core and valence orbitals are represented by the same number of primitive Gaussian functions [math]\displaystyle{ \mathbf \phi_i }[/math]. For example, an STO-3G basis set for the 1s, 2s and 2p orbital of the carbon atom are all linear combination of 3 primitive Gaussian functions. For example, a STO-3G s orbital is given by:

[math]\displaystyle{ \mathbf \psi_{\mathrm{STO}-3\mathrm{G}}(s)=c_1\phi_1 + c_2\phi_2 + c_3\phi_3 }[/math]

where

[math]\displaystyle{ \mathbf \phi_1 = \left (\frac{2\alpha_1}{\pi} \right ) ^{3/4}e^{-\alpha_1 r^2} }[/math]
[math]\displaystyle{ \mathbf \phi_2 = \left (\frac{2\alpha_2}{\pi} \right ) ^{3/4}e^{-\alpha_2 r^2} }[/math]
[math]\displaystyle{ \mathbf \phi_3 = \left (\frac{2\alpha_3}{\pi} \right ) ^{3/4}e^{-\alpha_3 r^2} }[/math]

The values of c1, c2, c3, α1, α2 and α3 have to be determined. For the STO-nG basis sets, they are obtained by making a least squares fit of the three Gaussian orbitals to the single Slater-type orbitals. (Extensive tables of parameters have been calculated for STO-1G through STO-5G for s orbitals through g orbitals.)[1] This differs from the more common procedure where the criteria often used is to choose the coefficients (c's) and exponents (α's) to give the lowest energy with some appropriate method for some appropriate molecule. A special feature of this basis set is that common exponents are used for orbitals in the same shell (e.g. 2s and 2p) as this allows more efficient computation.[2]

The fit between the Gaussian orbitals and the Slater orbital is good for all values of r, except for very small values near to the nucleus. The Slater orbital has a cusp at the nucleus, while Gaussian orbitals are flat at the nucleus.[3][4]

Use of STO-nG basis sets

The most widely used basis set of this group is STO-3G, which is used for large systems and for preliminary geometry determinations. This basis set is available for all atoms from hydrogen up to xenon.[5]

STO-2G basis set

The STO-2G basis set is a linear combination of 2 primitive Gaussian functions. The original coefficients and exponents for first-row and second-row atoms are given as follows.[2]

STO-2G α1 c1 α2 c2
1s 0.151623 0.678914 0.851819 0.430129
2s 0.0974545 0.963782 0.384244 0.0494718
2p 0.0974545 0.61282 0.384244 0.511541

Accuracy

The exact energy of the 1s electron of H atom is −0.5 hartree, given by a single Slater-type orbital with exponent 1.0. The following table illustrates the increase in accuracy as the number of primitive Gaussian functions increases from 3 to 6 in the basis set.[2]

Basis set Energy [hartree]
STO-3G −0.49491
STO-4G −0.49848
STO-5G −0.49951
STO-6G −0.49983

See also

  • List of quantum chemistry and solid state physics software

References

  1. Stewart, Robert F. (1 January 1970). "Small Gaussian Expansions of Slater‐Type Orbitals". The Journal of Chemical Physics 52 (1): 431–438. doi:10.1063/1.1672702. 
  2. 2.0 2.1 2.2 Hehre, W. J.; R. F. Stewart; J. A. Pople (1969). "Self-Consistent Molecular-Orbital Methods. I. Use of Gaussian Expansions of Slater-Type Atomic Orbitals". Journal of Chemical Physics 51 (6): 2657–2664. doi:10.1063/1.1672392. Bibcode1969JChPh..51.2657H. 
  3. Chemical Modeling From Atoms to Liquids, Alan Hinchliffe, John Wiley & Sons, Ltd., 1999. pg 294.
  4. Molecular Modelling, Andrew R. Leach, Longman, 1996. pg 68 – 73.
  5. Computational Chemistry, David Young, Wiley-Interscience, 2001. pg 86.