Extended Hückel method

The extended Hückel method is a semiempirical quantum chemistry method, developed by Roald Hoffmann since 1963.[1] It is based on the Hückel method but, while the original Hückel method only considers pi orbitals, the extended method also includes the sigma orbitals.

The extended Hückel method can be used for determining the molecular orbitals, but it is not very successful in determining the structural geometry of an organic molecule. It can however determine the relative energy of different geometrical configurations. It involves calculations of the electronic interactions in a rather simple way for which the electron-electron repulsions are not explicitly included and the total energy is just a sum of terms for each electron in the molecule. The off-diagonal Hamiltonian matrix elements are given by an approximation due to Wolfsberg and Helmholz that relates them to the diagonal elements and the overlap matrix element.[2]

K is the Wolfsberg–Helmholz constant, and is usually given a value of 1.75. In the extended Hückel method, only valence electrons are considered; the core electron energies and functions are supposed to be more or less constant between atoms of the same type. The method uses a series of parametrized energies calculated from atomic ionization potentials or theoretical methods to fill the diagonal of the Fock matrix. After filling the non-diagonal elements and diagonalizing the resulting Fock matrix, the energies (eigenvalues) and wavefunctions (eigenvectors) of the valence orbitals are found.

It is common in many theoretical studies to use the extended Hückel molecular orbitals as a preliminary step to determining the molecular orbitals by a more sophisticated method such as the CNDO/2 method and ab initio quantum chemistry methods. Since the extended Hückel basis set is fixed, the monoparticle calculated wavefunctions must be projected to the basis set where the accurate calculation is to be done. One usually does this by adjusting the orbitals in the new basis to the old ones by least squares method. As only valence electron wavefunctions are found by this method, one must fill the core electron functions by orthonormalizing the rest of the basis set with the calculated orbitals and then selecting the ones with less energy. This leads to the determination of more accurate structures and electronic properties, or in the case of ab initio methods, to somewhat faster convergence.

The method was first used by Roald Hoffmann who developed, with Robert Burns Woodward, rules for elucidating reaction mechanisms (the Woodward–Hoffmann rules). He used pictures of the molecular orbitals from extended Hückel theory to work out the orbital interactions in these cycloaddition reactions.

A closely similar method was used earlier by Hoffmann and William Lipscomb for studies of boron hydrides.[3][4][5] The off-diagonal Hamiltonian matrix elements were given as proportional to the overlap integral.[6]

This simplification of the Wolfsberg and Helmholz approximation is reasonable for boron hydrides as the diagonal elements are reasonably similar due to the small difference in electronegativity between boron and hydrogen.

The method works poorly for molecules that contain atoms of very different electronegativity. To overcome this weakness, several groups have suggested iterative schemes that depend on the atomic charge. One such method, that is still widely used in inorganic and organometallic chemistry is the Fenske-Hall method.[7][8][9]

A program for the extended Hückel method is YAeHMOP which stands for "yet another extended Hückel molecular orbital package".[10] YAeHMOP has also been merged with the Avogadro open-source molecular editor and visualizer to enable calculations directly from the Avogadro graphical user interface for materials that are periodic in one, two, or three dimensions.[11] This integration also enables visualization of band structures, total and projected density of states, and crystal orbital overlap/Hamilton populations (COOPs/COHPs).

See also

edit

References

edit
  1. ^ Hoffmann, R. (1963). "An Extended Hückel Theory. I. Hydrocarbons". J. Chem. Phys. 39 (6): 1397–1412. Bibcode:1963JChPh..39.1397H. doi:10.1063/1.1734456.
  2. ^ M. Wolfsberg; L. J. Helmholz (1952). "The Spectra and Electronic Structure of the Tetrahedral Ions MnO4−, CrO4−−, and ClO4−". J. Chem. Phys. 20 (5): 837. Bibcode:1952JChPh..20..837W. doi:10.1063/1.1700580.
  3. ^ R. Hoffmann; W. N. Lipscomb (1962). "Theory of Polyhedral Molecules. I. Physical Factorizations of the Secular Equation". J. Chem. Phys. 36 (8): 2179. Bibcode:1962JChPh..36.2179H. doi:10.1063/1.1732849.
  4. ^ R. Hoffmann; W. N. Lipscomb (1962). "Boron Hydrides: LCAO—MO and Resonance Studies". J. Chem. Phys. 37 (12): 2872. Bibcode:1962JChPh..37.2872H. doi:10.1063/1.1733113.
  5. ^ W. N. Lipscomb Boron Hydrides, W. A. Benjamin Inc., New York, 1963, Chapter 3
  6. ^ Avery, Patrick. "Avogadro with YAeHMOP - Installation". avogadro-yaehmop.github.io. Retrieved 2023-03-15.
  7. ^ Charles Edwin Webster; Michael B. Hall (2005). "Chapter 40. Forty years of Fenske-Hall molecular orbital theory". Theory and Applications of Computational Chemistry: The First Forty Years. pp. 1143–1165. doi:10.1016/B978-044451719-7/50083-4. ISBN 978-0-444-51719-7.
  8. ^ Hall, M. B.; Fenske, R. F. (1972). "Electronic structure and bonding in methyl- and perfluoromethyl(pentacarbonyl)manganese". Inorg. Chem. 11 (4): 768. doi:10.1021/ic50110a022.
  9. ^ jimp2 program
  10. ^ Computational Chemistry, David Young, Wiley-Interscience, 2001. Appendix A. A.3.3 pg 343, YAeHMOP
  11. ^ Avery, Patrick; Ludowieg, Herbert; Autschbach, Jochen; Zurek, Eva (2018-02-13). "Extended Hückel Calculations on Solids Using the Avogadro Molecular Editor and Visualizer". Journal of Chemical Education. 95 (2): 331–337. doi:10.1021/acs.jchemed.7b00698. ISSN 0021-9584.
edit
  • [1] Online Extended Hückel Calculator.