Open main menu
Home
Random
Recent changes
Special pages
Community portal
Preferences
About Wikipedia
Disclaimers
Incubator escapee wiki
Search
User menu
Talk
Dark mode
Contributions
Create account
Log in
Editing
Gaussian orbital
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
In [[computational chemistry]] and [[molecular physics]], '''Gaussian orbitals''' (also known as '''Gaussian type orbitals''', '''GTOs''' or '''Gaussians''') are [[function (mathematics)|function]]s used as [[atomic orbital]]s in the [[linear combination of atomic orbitals molecular orbital method|LCAO method]] for the representation of [[Molecular orbital|electron orbitals]] in [[molecule]]s and numerous properties that depend on these.<ref>{{cite journal|last=Gill|first=Peter M.W.|title=Molecular integrals Over Gaussian Basis Functions|journal=Advances in Quantum Chemistry|year=1994|volume=25|pages=141–205|doi=10.1016/S0065-3276(08)60019-2|url=http://cs.anu.edu.au/~Alistair.Rendell/hons09/qc/GillReview.pdf|accessdate=17 June 2011|bibcode = 1994AdQC...25..141G |isbn=9780120348251}}</ref> == Rationale == The use of Gaussian orbitals in electronic structure theory (instead of the more physical [[Slater-type orbital]]s) was first proposed by [[S. Francis Boys|Boys]]<ref>{{cite journal|first1=S. F.|last1=Boys|title=Electronic Wave Functions. I. A General Method of Calculation for the Stationary States of Any Molecular System|jstor=98423 | journal= Proc. R. Soc. Lond. A |volume=200 |issue=1063|pages=542–554|year=1950|bibcode=1950RSPSA.200..542B|doi=10.1098/rspa.1950.0036|s2cid=122709395|url=http://elib.bsu.by/handle/123456789/154387}}</ref> in 1950. The principal reason for the use of Gaussian [[Basis set (chemistry)|basis functions]] in molecular quantum chemical calculations is the 'Gaussian Product Theorem', which guarantees that the product of two GTOs centered on two different atoms is a finite sum of Gaussians centered on a point along the axis connecting them. In this manner, four-center integrals can be reduced to finite sums of two-center integrals, and in a next step to finite sums of one-center integrals. The speedup by 4-5 orders of magnitude compared to [[Slater orbitals]] outweighs the extra cost entailed by the larger number of basis functions generally required in a Gaussian calculation. For reasons of convenience, many quantum chemistry programs work in a basis of Cartesian Gaussians even when spherical Gaussians are requested, as integral evaluation is much easier in the Cartesian basis, and the spherical functions can be simply expressed using the Cartesian functions.<ref>{{cite journal|last=Schlegel|first=H.|author2=Frisch, M.|title=Transformation between Cartesian and pure spherical harmonic Gaussians|journal=International Journal of Quantum Chemistry|year=1990|volume=54|issue=2|pages=83–87|doi=10.1002/qua.560540202 |s2cid=94417974 }}</ref> <ref> {{cite journal|last=Mathar|first=Richard J.|title=Mutual Conversion of Three Flavors of Gaussian Type Orbitals|journal=International Journal of Quantum Chemistry|year=2002|volume=90|issue=1|pages=227–243|doi=10.1002/qua.10085|arxiv=physics/9907051 |bibcode=2002IJQC...90..227M|s2cid=119100125 }} </ref> == Mathematical form == The Gaussian basis functions obey the usual radial-angular decomposition : <math>\ \Phi(\mathbf{r}) = R_l(r) Y_{lm}(\theta,\phi)</math>, where <math>Y_{lm}(\theta,\phi)</math> is a [[Spherical harmonics|spherical harmonic]], <math>l</math> and <math>m</math> are the angular momentum and its <math>z</math> component, and <math>r,\theta,\phi</math> are spherical coordinates. While for Slater orbitals the radial part is : <math>\ R_l(r) = A(l,\alpha) r^l e^{-\alpha r}, </math> <math>A(l,\alpha)</math> being a normalization constant, for Gaussian primitives the radial part is : <math>\ R_l(r) = B(l,\alpha) r^l e^{-\alpha r^2},</math> where <math>B(l,\alpha)</math> is the normalization constant corresponding to the Gaussian. The normalization condition which determines <math>A(l,\alpha)</math> or <math>B(l,\alpha)</math> is :<math>\int _0 ^\infty \mathrm{d}r \, r^2 \left| R_l (r) \right|^2 = 1</math> which in general does not impose orthogonality in <math>l</math>. Because an individual primitive Gaussian function gives a rather poor description for the electronic wave function near the nucleus, Gaussian basis sets are almost always contracted: :<math>\ R_l(r) = r^l \sum_{p=1}^P c_p B(l,\alpha_p) \exp(-\alpha_p r^2)</math>, where <math>c_p</math> is the contraction coefficient for the primitive with exponent <math>\alpha_p</math>. The coefficients are given with respect to normalized primitives, because coefficients for unnormalized primitives would differ by many orders of magnitude. The exponents are reported in [[atomic units]]. There is a large library of published Gaussian basis sets optimized for a variety of criteria available at the [http://basissetexchange.org/ Basis Set Exchange portal]. ===Cartesian coordinates=== In Cartesian coordinates, Gaussian-type orbitals can be written in terms of exponential factors in the <math>x</math>, <math>y</math>, and <math>z</math> directions as well as an exponential factor <math>\alpha</math> controlling the width of the orbital. The expression for a Cartesian Gaussian-type orbital, with the appropriate normalization coefficient is :<math>\Phi(x,y,z;\alpha,i,j,k)=\left(\frac{2\alpha}{\pi}\right)^{3/4}\left[\frac{(8\alpha)^{i+j+k}i!j!k!}{(2i)!(2j)!(2k)!}\right]^{1/2}x^i y^j z^k e^{-\alpha(x^2+y^2+z^2)}</math> In the above expression, <math>i</math>, <math>j</math>, and <math>k</math> must be integers. If <math>i+j+k=0</math>, then the orbital has spherical symmetry and is considered an s-type GTO. If <math>i+j+k=1</math>, the GTO possesses axial symmetry along one axis and is considered a p-type GTO. When <math>i+j+k=2</math>, there are six possible GTOs that may be constructed; this is one more than the five canonical d orbital functions for a given angular quantum number. To address this, a linear combination of two d-type GTOs can be used to reproduce a canonical d function. Similarly, there exist 10 f-type GTOs, but only 7 canonical f orbital functions; this pattern continues for higher angular quantum numbers.<ref>{{cite book |last1=Cramer |first1=Christopher J. |title=Essentials of computational chemistry : theories and models |date=2004 |publisher=Wiley |location=Chichester, West Sussex, England |isbn=9780470091821 |pages=167 |edition=2nd}}</ref> ==Molecular integrals== Taketa et al. (1966) presented the necessary mathematical equations for obtaining matrix elements in the Gaussian basis.<ref>{{cite journal|last=Taketa|first=Hiroshi|author2=Huzinaga, Sigeru |author3=O-ohata, Kiyosi |title=Gaussian-Expansion Methods for Molecular Integrals|journal=Journal of the Physical Society of Japan|year=1966|volume=21|issue=11|pages=2313–2324|doi=10.1143/JPSJ.21.2313|bibcode = 1966JPSJ...21.2313T }}</ref> Since then much work has been done to speed up the evaluation of these integrals which are the slowest part of many quantum chemical calculations. Živković and Maksić (1968) suggested using [[Hermite]] Gaussian functions,<ref>{{cite journal|last=Živković|first=T.|author2=Maksić, Z. B.|title=Explicit Formulas for Molecular Integrals over Hermite-Gaussian Functions|journal=Journal of Chemical Physics|year=1968|volume=49|issue=7|pages=3083–3087|doi=10.1063/1.1670551 |bibcode = 1968JChPh..49.3083Z }}</ref> as this simplifies the equations. McMurchie and Davidson (1978) introduced recursion relations,<ref>{{cite journal|last=McMurchie|first=Larry E.|author2=Davidson, Ernest R.|title=One- and two-electron integrals over Cartesian Gaussian functions|journal=Journal of Computational Physics|year=1978|volume=26|issue=2|pages=218–31|doi=10.1016/0021-9991(78)90092-X|bibcode = 1978JCoPh..26..218M |url=https://digital.library.unt.edu/ark:/67531/metadc1057442/}}</ref> which greatly reduces the amount of calculations. [[John Pople|Pople]] and Hehre (1978) developed a local coordinate method.<ref>{{cite journal|last=Pople|first=J. A.|author2=Hehre, W. J.|title=Computation of electron repulsion integrals involving contracted Gaussian basis functions.|journal=J. Comput. Phys.|year=1978|volume=27|issue=2|pages=161–168|doi=10.1016/0021-9991(78)90001-3|bibcode = 1978JCoPh..27..161P }}</ref> Obara and Saika introduced efficient recursion relations in 1985,<ref>{{cite journal|last=Obara|first=S.|author2=Saika, A.|title=Efficient recursive computation of molecular integrals over Cartesian Gaussian functions|journal=J. Chem. Phys.|year=1986|volume=84|issue=7|pages=3963–74|doi=10.1063/1.450106|bibcode = 1986JChPh..84.3963O }}</ref> which was followed by the development of other important recurrence relations. Gill and Pople (1990) introduced a 'PRISM' algorithm which allowed efficient use of 20 different calculation paths.<ref>{{cite journal|last=Gill|first=Peter M. W.|author2=Pople, John A.|title=The Prism Algorithm for Two-Electron Integrals|journal=International Journal of Quantum Chemistry|date=December 1991|volume=40|issue=6|pages=753–772|doi=10.1002/qua.560400605|url=http://rscweb.anu.edu.au/~pgill/papers/026PRISM.pdf|accessdate=17 June 2011}}</ref> ==The POLYATOM System== The POLYATOM System<ref>{{cite journal|first1=I.G.|last1=Csizmadia|first2=M.C.|last2=Harrison|first3=J.W.|last3=Moskowitz|first4=B.T.|last4=Sutcliffe|title=Nonempirical LCAO-MO-SCF-CI calculations on organic molecules with gaussian-type functions. Introductory review and mathematical formalism|journal=Theoretica Chimica Acta|volume=6|number=3|page=191|year=1966|doi=10.1007/BF02394698|s2cid=198176437 }}</ref> was the first package for ''ab initio'' calculations using Gaussian orbitals that was applied to a wide variety of molecules.<ref>A.C. Wahl, ''Chemistry by computer'', Scientific American, pages 54-70, April, 1970.</ref> It was developed in [[John Clarke Slater|Slater's]] Solid State and Molecular Theory Group (SSMTG) at [[MIT]] using the resources of the Cooperative Computing Laboratory.<!-- verification in acknowledgements in TCA paper--> The mathematical infrastructure and operational software were developed by Imre Csizmadia,<ref>[https://books.google.com/books?id=3xCzIQDYOvoC&dq=Imre+Csizmadia+quantum&pg=PA248 ''Imre Csizmadia'', Professor Emeritus of Chemistry, University of Toronto, in Reviews in Computational Chemistry, vol.15, p.248]</ref> Malcolm Harrison,<ref>[http://cs.nyu.edu/cs/review95/node12.html ''Malcolm C. Harrison'', Professor of Computer Science, New York University]</ref> Jules Moskowitz<ref>{{dead link|date=April 2025}}[https://www.google.com/#hl=en&sugexp=ldymls&xhr=t&q=Jules+Moskowitz&cp=10&pf=p&sclient=psy&aq=0&aqi=&aql=&oq=Jules+Mosk&pbx=1&bav=on.1,or.&fp=5e2b21bd614e0a97 ''Jules W. Moskowitz'', Professor Emeritus of Chemistry, New York University]</ref> and Brian Sutcliffe.<ref>[http://www.chilton-computing.org.uk/acl/associates/users/sutcliffe.htm ''Brian T. Sutcliffe'', Professor of Chemistry, York University]</ref> ==See also== *[[Quantum chemistry computer programs]] ==References== {{reflist}} ==External links== * [http://www.shef.ac.uk/chemistry/orbitron/ A visualization of all common and uncommon atomic orbitals, from 1s to 7g] (''Note that the radial part of the expressions given corresponds to [[Slater orbitals]] rather than Gaussians. The angular parts, and hence their shapes as displayed in figures, are the same as those of spherical Gaussians.'') * [http://www.ccl.net/cca/documents/basis-sets/basis.html Explanation of Gaussian basis set] * [https://bse.pnl.gov/bse/portal Basis set exchange] * {{cite journal|url=http://lnu.diva-portal.org/smash/get/diva2:282089/FULLTEXT01 |title=Detailed derivation of Gaussian orbital based matrix elements in electron structure calculations |doi=10.1088/0143-0807/31/1/004 |year=2010|volume=31 |issue=1 | page=37 | first1=T | last1=Petersson | first2=B. |last2=Hellsing | journal= Eur. J. Phys.|bibcode=2010EJPh...31...37P|s2cid=122528581 }} [[Category:Molecular physics]] [[Category:Quantum chemistry]] [[Category:Computational chemistry]]
Edit summary
(Briefly describe your changes)
By publishing changes, you agree to the
Terms of Use
, and you irrevocably agree to release your contribution under the
CC BY-SA 4.0 License
and the
GFDL
. You agree that a hyperlink or URL is sufficient attribution under the Creative Commons license.
Cancel
Editing help
(opens in new window)
Pages transcluded onto the current version of this page
(
help
)
:
Template:Cite book
(
edit
)
Template:Cite journal
(
edit
)
Template:Dead link
(
edit
)
Template:Reflist
(
edit
)