A MATLAB code for Hartree Fock calculation of H-H ground state bondlength and energy using STO-4G
Licensed under Creative Commons.
Category
Published on
Abstract
Hartree Fock (HF) theory is one of the basic theories underlying the current understanding of the electronic structure of materials. It is a simple non-relativistic treatment of many electron system that accounts for the antisymmetric (fermion) nature of electronic wavefunction but does not account for correlation effects. The matrix representation of this theory constitutes the Roothan Equations. This simple MATLAB Code uses the Roothan Equations, within Born Oppenheimer approximation and non-relativistic limit, to calculate the energy of H-H (two Hydrogen atoms) system with different bond lengths using Slater Type Orbital with 4 Gaussian basis functions (STO-4G) for each spatial orbital. It can be used as an example for a class on electronic structure theory or as a starting skeleton for developing more complex ab initio codes, that use a similar basis set.
Credits
Prof. Jorge Rodriguez @ Physics, Purdue University for his lectures on Hartree Fock theory (Phys 570A - Computational Biomolecular Physics), J.M. Thijssen for his book "Computational Physics"
Cite this work
Researchers should cite this work as follows: