Introduction to Computational Chemistry

Table of Contents

• Introduction
• Ab Initio
• Semiempirical
• Modeling the solid state
• Molecular Mechanics
• Molecular Dynamics
• Statistical Mechanics
• Thermodynamics
• Structure-Property Relationships
• Symbolic Calculations
• Artifical Intelligence
• How to do a computational research project
• Visualization
• Further information

Introduction

Many universities are now offering classes, which are an overview of various aspects of computational chemistry. Since we have had many people wanting to start doing computations before they have had even an introductory course, this document has been written as step one in understanding what computational chemistry is about. Note that this is not intended to teach the fundamentals of chemistry, quantum mechanics or mathematics, only most basic description of how chemical computations are done.

The term theoretical chemistry may be defined as the mathematical description of chemistry. The term computational chemistry is usually used when a mathematical method is sufficiently well developed that it can be automated for implementation on a computer. Note that the words exact and perfect do not appear in these definitions. Very few aspects of chemistry can be computed exactly, but almost every aspect of chemistry has been described in a qualitative or approximate quantitative computational scheme. The biggest mistake that a computational chemists can make is to assume that any computed number is exact. However, just as not all spectra are perfectly resolved, often a qualitative or approximate computation can give useful insight into chemistry if you understand what it tells you and what it doesn't.

Although most chemists avoid the true paper & pencil type of theoretical chemistry, keep in mind that this is what many Nobel prizes have been awarded for.

Ab Initio

The most common type of ab initio calculation is called a Hartree Fock calculation (abbreviated HF), in which the primary approximation is called the central field approximation. This means that the Coulombic electron-electron repulsion is not specifically taken into account. However, it's net effect is included in the calculation. This is a variational calculation, meaning that the approximate energies calculated are all equal to or greater than the exact energy. The energies calculated are usually in units called Hartrees (1 H = 27.2114 eV). Because of the central field approximation, the energies from HF calculations are always greater than the exact energy and tend to a limiting value called the Hartree Fock limit.

The second approximation in HF calculations is that the wave function must be described by some functional form, which is only known exactly for a few one electron systems. The functions used most often are linear combinations of Slater type orbitals exp(-ax) or Gaussian type orbitals exp(-ax^2), abbreviated STO and GTO. The wave function is formed from linear combinations of atomic orbitals or more often from linear combinations of basis functions. Because of this approximation, most HF calculations give a computed energy greater than the Hartree Fock limit. The exact set of basis functions used is often specified by an abbreviation, such as STO-3G or 6-311++g**.

A number of types of calculations begin with a HF calculation then correct for the explicit electron-electron repulsion, referred to as correlation. Some of these methods are Mohlar-Plesset perturbation theory (MPn, where n is the order of correction), the Generalized Valence Bond (GVB) method, Multi-Configurations Self Consistent Field (MCSCF), Configuration Interaction (CI) and Coupled Cluster theory (CC). As a group, these methods are referred to as correlated calculations.

A method, which avoids making the HF mistakes in the first place is called Quantum Monte Carlo (QMC). There are several flavors of QMC .. variational, diffusion and Green's functions. These methods work with an explicitly correlated wave function and evaluate integrals numerically using a Monte Carlo integration. These calculations can be very time consuming, but they are probably the most accurate methods known today.

An alternative ab initio method is Density Functional Theory (DFT), in which the total energy is expressed in terms of the total electron density, rather than the wavefunction. In this type of calculation, there is an approximate Hamiltonian and an approximate expression for the total electron density.

The good side of ab initio methods is that they eventually converge to the exact solution, once all of the approximations are made sufficiently small in magnitude. However, this convergence is not montonic. Sometimes, the smallest calculation gives the best result for a given property.

The bad side of ab initio methods is that they are expensive. These methods often take enormous amounts of computer cpu time, memory and disk space. The HF method scales as N^4, where N is the number of basis functions, so a calculation twice as big takes 16 times as long to complete. Correlated calculations often scale much worse than this. In practice, extremely accurate solutions are only obtainable when the molecule contains half a dozen electrons or less.

In general, ab initio calculations give very good qualitative results and can give increasingly accurate quantitative results as the molecules in question become smaller.

Semiempirical

The good side of semiempirical calculations is that they are much faster than the ab initio calculations.

The bad side of semiempirical calculations is that the results can be eratic. If the molecule being computed is similar to molecules in the data base used to parameterize the method, then the results may be very good. If the molecule being computed is significantly different from anything in the parameterization set, the answers may be very poor.

Semiempirical calculations have been very successful in the description of organic chemistry, where there are only a few elements used extensively and the molecules are of moderate size. However, semiempirical methods have been devised specifically for the description of inorganic chemistry as well.

Modeling the solid state

Band structure calculations have been done for very complicated systems, however the software is not yet automated enough or sufficiently fast that anyone does band structures casually. If you want to do band structure calculations, you had better expect to put a lot of time into your efforts.

Molecular Mechanics

In a molecular mechanics method, the data base of compounds used to parameterize the method (a set of parameters and functions is called a force field) is crucial to it's success. Where as a semiempirical method may be parameterized against a set of organic molecules, a molecular mechanics method may be parameterized against a specific class of molecules, such as proteins. Such a force field would only be expected to have any relevance to describing other proteins.

The good side of molecular mechanics is that it allows the modeling the enormous molecules, such as proteins and segments of DNA, making it the primary tool of computational biochemists.

The bad side of molecular mechanics is that there are many chemical properties that are not even defined within the method, such as electronic excited states. In order to work with extremely large and complicated systems, often molecular mechanics software packages have the most powerful and easiest to use graphical interfaces. Because of this, mechanics is sometimes used because it is easy, but not necessarily a good way to describe a system.

Molecular Dynamics

The application of molecular dynamics to solvent/solute systems allows the computation of properties such as diffusion coeficients or radial distribution functions for use in statistical mechanical treatments. Usually the scheme of a solvent/solute calculation is that a number of molecules (perhaps 1000) are given some initial position and velocity. New positions are calculated a small time later based on this movement and this process is itterated for thousands of steps in order to bring the system to equilibrium and give a good statistical description of the radial distribution function.

In order to analyze the vibrations of a single molecule, many dynamics steps are done, then the data is Fourier transformed into the frequency domain. A given peak can be chosen and transformed back to the time domain, in order to see what the motion at that frequency looks like.

Statistical Mechanics

Thermodynamics

Structure-Property Relationships

The simplest case of structure-property relationships are qualitative thumb rules. For example, an experienced polymer chemist may be able to predict whether a polymer will be soft or brittle based on the geometry and bonding of the monomers.

When structure-property relationships are mentioned in current literature, it usually implies a quantitative mathematical relationship. These relationships are most often derived by using curve fitting software to find the linear combination of molecular properties, which best reproduces the desired property. The molecular properties are usually obtained from molecular modeling computations. Other molecular descriptors such as molecular weight or topological descriptions are also used.

When the property being described is a physical property, such as the boiling point, this is refered to as a Quantitative Structure-Property Relationship (QSPR). When the property being described is a type of biological activity (such as drug activity), this is refered to as a Quantitative Structure-Activity Relationship (QSAR).

Symbolic Calculations

Artifical Intelligence

How to do a computational research project

What do you want to know? How accurately? Why? If you can't answer these questions, then you don't even have a research project yet.

How accurate do you predict the answer will be? In analytical chemistry, you do a number of identical measurements then work out the error from a standard deviation. With computational experiments, doing the same thing should always give exactly the same result. The way that you estimate your error is to compare a number of similar computations to the experimental answers. There are articles and compilations of these studies. If none exist, you will have to guess which method should be reasonable, based on it's assumptions then do a study yourself, before you can apply it to you unknown and have any idea how good the calculation is. When someone just tells you off the top of their head what method to use, they either have a fair amount of this type of information memorized, or they don't know what they are talking about. Beware of someone who tells you a given program is good just because it is the only one they know how to use, rather than the basing their answer on the quality of the results.

How long do you expect it to take? If the world were perfect, you would tell your PC (voice input of course) to give you the exact solution to the Schrödinger equation and go on with your life. However, often ab initio calculations would be so time consuming that it would take a decade to do a single calculation, if you even had a machine with enough memory and disk space. However, a number of methods exist because each is best for some situation. The trick is to determine which one is best for your project. Again, the answer is to look into the literature and see how long each takes. If the only thing you know is how a calculation scales, do the simplest possible calculation then use the scaling equation to estimate how long it will take to do the sort of calculation that you have predicted will give the desired accuracy.

What approximations are being made? Which are significant? This is how you avoid looking like a complete fool, when you successfully perform a calculation that is complete garbage. An example would be trying to find out about vibrational motions that are very anharmonic, when the calculation uses a harmonic oscillator approximation.

Once you have finally answered all of these questions, you are ready to actually do a calculation. Now you must determine what software is available, what it costs and how to use it. Note that two programs of the same type (i.e. ab initio) may calculate different properties, so you have to make sure the program does exactly what you want.

When you are learning how to use a program, you will probably try to do dozens of calculations that will fail because you constructed the input incorrectly. Do not use your project molecule to do this. Make all your mistakes with something really easy, like a water molecule. That way you don't waste enormous amounts of time.

Visualization

Further information

There are many books on the principles of quantum mechanics and every physical chemistry text has an introductory treatment. The work which I am listing here is a two volume set with each chapter broken into a basic and advanced sections making it excellent for both intermediate and advanced users.
C. Cohen-Tannoudji, B. Diu, F. Laloe "Quantum Mechanics Volumes I & II" Wiley-Interscience (1977)

For an introduction to quantum chemistry see
D. A. McQuarrie "Quantum Chemistry" University Science Books (1983)

A graduate level text on quantum chemistry is
I. N. Levine "Quantum Chemistry" Prentice Hall (1991)

For quantum Monte Carlo methods, order the following book using ISBN 981-02-0322-5 because the title is listed incorrectly in 'Books in Print'.
B. L. Hammond, W. A. Lester, Jr., P. J. Reynolds "Monte Carlo Methods in Ab Initio Quantum Chemistry" World Scientific (1994)

For density functional theory see
R. G. Parr, W. Yang "Density-Functional Theory of Atoms and Molecules" Oxford (1989)

For a basic understanding of solid state modeling see
R. Hoffmann "Solids and Surfaces : A Chemist's View of Bonding in Extended Structures", VCH (1988)

For a graduate level description of statistical mechanics see
D. A. McQuarrie "Statistical Mechanics" Harper Collins (1976)

Any physical chemistry text will have a description of thermodynamics but I will recommend
I. N. Levine "Physical Chemistry" McGraw Hill (1995)

There is a comprehensive listing of all available molecular modeling software and structural databanks, free or not, in appendix 2 of
"Reviews in Computational Chemistry Volume 6" Ed. K. B. Lipkowitz and D. B. Boyd, VCH (1995)

There is a write up on computer aided drug design at
gopher://ccl.osc.edu/00/documents/drug.design.guide

Mathematical challenges from theoretical/computational chemistry
http://www.nap.edu/readingroom/books/mctcc/index.html

An online text on molecular modeling using molecular mechanics
http://www.awod.com/netsci/Science/Compchem/feature01.html

A Computational Chemistry Primer
http://www.sdsc.edu/GatherScatter/GSwinter96/taylor1.html

An online text on computational chemistry
http://www.cryst.bbk.ac.uk/~ubcg8ab/course/os_molf.html

Another online text on quantum chemistry
http://zopyros.ccqc.uga.edu/Docs/Knowledge/Fundamental_Theory/quantrev/node1.html

An online introduction to quantum mechanics is at
http://cmcind.far.ruu.nl/webcmc/qm/home.html