Module 1

Introduction to Computational Chemistry

Introduction

Computational chemistry is a specialized branch of chemical science that uses computers,

mathematical models, and simulations to understand, explain, and predict chemical
behaviour. It acts as a bridge between theoretical chemistry and experimental chemistry.

In traditional chemistry, experiments are carried out in laboratories using chemicals and
instruments. However, many chemical problems are:

 Too complex to solve analytically

 Time-consuming and costly to investigate experimentally
 Difficult or unsafe to study in the laboratory

Computational chemistry addresses these challenges by performing virtual experiments on

computers, allowing scientists and engineers to explore molecular systems efficiently

Computational chemistry is important because it:

 Reduces dependence on trial-and-error laboratory experiments

 Saves time, cost, and chemical resources
 Helps visualize molecules at the atomic and molecular level
 Supports experimental chemistry by providing theoretical insight

Real-world relevance

 Drug discovery and development

 Battery and energy-storage materials
 Catalysts and industrial chemicals
 Environmental and climate-related materials
 Biochemical and pharmaceutical research

Note : Computational chemistry does not replace experiments, but guides and
complements them.

Scope and Key Applications

Computational chemistry is applied across various domains of chemical science and

engineering.

[Link] Geometry
Molecular geometry refers to the three-dimensional arrangement of atoms in a molecule.

Computational chemistry helps to:

 Determine bond lengths and bond angles

 Predict molecular shape
 Identify stable molecular conformations

Example:
Why is the water molecule bent instead of linear?
Computational calculations provide an explanation based on energy minimization.

2. Molecular Spectroscopy

Spectroscopy studies the interaction of molecules with electromagnetic radiation.

Computational methods are used to:

 Predict IR, UV–Visible, and NMR spectra

 Interpret experimental spectral data
 Assign vibrational and electronic transitions

It helps confirm molecular structures and understand experimental observations.

3. Potential Energy Surfaces (PES)

A Potential Energy Surface describes how the energy of a molecular system changes with
atomic positions.

 Identifies stable structures (energy minima)

 Helps locate transition states
 Explains reaction pathways and reaction barriers

A central concept in modelling is that the energy of a molecule is a function of its

coordinates. The Potential Energy Surface is a multidimensional "map" representing how
the potential energy changes as the molecular geometry varies.
• Stationary Points: These occur where the first derivative of the energy with respect to all
coordinates is zero.
• Minima: Points on the PES that correspond to stable structures. A "local minimum" is the
lowest energy point in a specific region, while the "global minimum" is the lowest energy
state for the entire surface.
• Saddle Points: These represent transition states between two minima, characterized by a
maximum in one direction and minima in all others
 

4. Drug Design

One of the most impactful applications of computational chemistry is drug design.

It is used to:

 Design and screen drug molecules

 Study drug–protein interactions
 Reduce time and cost in pharmaceutical development

Example:
Many modern drugs are first tested computationally before laboratory synthesis.

Computational Chemistry and Computer Science

Computational chemistry lies at the intersection of chemistry and computer science. It

involves:

 Developing mathematical models of molecular systems

 Implementing efficient algorithms to solve chemical equations
  Running simulations that would be impossible to perform manually

Computer science provides the algorithms and software, while chemistry provides
the physical and chemical principles.

Models Used in Computational Chemistry

Since real molecules are highly complex, computational chemistry uses theoretical models
to approximate their behaviour.

What is a model?

A model, on the other hand, typically involves the deliberate introduction of simplifying
approximations into a more general theory so as to extend its practical utility. For example
one would certainly not want to design the pendulum for a mechanical for a clock using the
complicated mathematics of quantum mechanics. It could be done to fairly the same
amount of accuracy using the simpler equations of the more restricted classical theories.
Another very useful model is the Lewis dot structure description of chemical bonding.
In a way, a model is a simplified way of describing and predicting the scientific results, which
might sometimes be an incomplete description. Models are very useful and they allow us to
predict and understand phenomena without performing complicated mathematical
manipulations dictated by a rigorous theory. Experienced researchers continue to use
models, however realizing that there will always be exceptions to the rules of these models.
Valance Shell Electron Pair Repulsion (VSEPR) theory is one such model used to predict the
geometry of molecules from the number of electron pairs surrounding a central atom. The
VSEPR model does not give a complete description of a molecule and does not predict many
properties of the molecule. However, it does a good job in predicting molecular geometry,
based on a very little information (the number of electron pairs).
Another example of a very useful model is the Lewis dot structure description of chemical
bonding. Lewis dot structures are the diagram that show the bonding between atoms in a
molecule and the lone pairs of electrons present in a molecule. It is not a complete
description of the molecule, since it does not contain the columbic interactions between the
electrons and the nuclei or other interactions. It also does not predict the structure of the
molecule. However, Lewis model takes into account the pairing of electrons. This results in
Lewis model being able to predict chemical bonding patterns and give an indications of the
strength of the bonds (single bonds, double bonds, etc.).
The Lewis dot structures of methane

Methods/Models used in Computational Chemistry

1. Quantum Mechanics (QM) Methods

Quantum mechanical methods describe electrons and nuclei explicitly.

i. Ab initio methods (e.g., Hartree–Fock, MP2, CCSD(T))

→ High accuracy, high computational cost
ii. Density Functional Theory (DFT)
→ Balance between accuracy and computational efficiency
iii. Semi-empirical methods
→ Faster calculations with additional approximations for large systems

These methods are essential for studying electronic structure and chemical bonding.

2. Molecular Mechanics and Molecular Dynamics (MM/MD)

These methods use classical physics to model molecular systems.

Molecular Mechanics (MM)

 Treats atoms as balls and bonds as springs

 Uses force fields (e.g., CHARMM, AMBER, OPLS)
 Suitable for large molecules like proteins and polymers

Molecular mechanics models molecules as collections of spheres (atoms) held together by

springs (bonds). The total energy of a system is calculated using a force field, which
typically includes the following terms:
 Bond Stretching: Energy associated with deviations from ideal bond lengths.
 Angle Bending: Energy associated with deviations from ideal bond angles.
 Torsional Terms: Energy changes resulting from rotation around single bonds (eg:
dihedral angle)
 Non-bonded Interactions: These include van der Waals forces and electrostatic
interactions (modelled using Coulomb's law).
Unlike quantum mechanics, molecular mechanics ignores electronic motion and treats atoms
as the smallest fundamental units.

Molecular Dynamics (MD)

 Molecular Dynamics (MD) is a computer simulation that acts like a virtual movie,
showing how atoms and molecules move and interact over time. Studies conformational
changes.

o MD treats atoms like small moving objects that follow Newton's Laws of Motion
(Force=Mass × Acceleration) . If you know where an atom starts and how fast it’s
going, you can predict where it will be next .
o The push or pull on each atom is calculated from the system's potential energy,
which models chemical bonds like springs .
o Because atoms move so fast, the computer updates their positions every femtosecond
(one quadrillionth of a second) . By connecting millions of these "snapshots," the
computer creates a trajectory, or a movie of the molecule's life .
o Unlike a static 3D picture, MD shows us the vibrations and shape-changes that
molecules undergo in the real world, allowing us to calculate averages for properties
like pressure or heat .

3. Monte Carlo (MC) Methods

 Use statistical sampling techniques

 Explore molecular configurations randomly
 Useful for studying equilibrium properties

4. Molecular Modelling and Simulations

Molecular Modelling

Molecular modelling is described as a set of techniques used to mimic the behaviour of

molecules. It ranges from simple graphical representations to complex dynamical
simulations of systems containing thousands of atoms. While once the domain of specialists,
it is now an integral tool in chemistry and biology for rationalizing experimental data and
making predictions

Molecular Simulations

 Observing how molecules move and interact over time

 Studying dynamic behaviour rather than static structures

Modelling is a photograph, simulation is a video.

Typical Workflow in Computational Chemistry
A standard computational chemistry study follows these steps:

1. Define the problem and select an appropriate method

2. Build the molecular model (structure, charge, spin state)
3. Optimize geometry to obtain stable structures
4. Perform frequency analysis to confirm minima or transition states
5. Calculate required properties (energies, spectra, barriers, etc.)
6. Validate results by comparison with experimental data
7. Document and report methods, parameters, and results

Software Tools Used in Computational Chemistry

Calculation Engines

 Gaussian
 GAMESS
 ORCA

Used for electronic structure and energy calculations.

Visualization and Modelling Tools

 PyMOL – molecular visualization

 Avogadro – molecular building and modelling

Chemical bonds

Ionic Bonds
Ionic bonding is a type of chemical bond that occurs when electrons are transferred from
one atom to another due to the difference in electronegativity. In this process, one atom
loses an electron, which is then gained by a different atom. When an atom loses an electron,
it becomes positively charged and is known as a cation. Conversely, the atom that gains the
electron becomes negatively charged and is called an anion. The electrostatic interaction
between these two ions keeps them together. Normally, one of the atoms which is involved
in the bond formation is a metal and the other one is a non-metal. The strength of an ionic
bond comes from the difference in charge between the two atoms. The greater the charge
difference, the stronger the ionic bond will be.
For example, in sodium chloride (NaCl), sodium (Na) loses one electron to become a sodium
ion (Na⁺), while chlorine (Cl) gains that electron to become a chloride ion (Cl⁻). This transfer
creates a strong ionic bond between the two ions. Another example is magnesium oxide
(MgO), where magnesium (Mg) loses two electrons to become Mg²⁺, and oxygen (O) gains
those electrons to become O²⁻, resulting in a strong ionic bond due to the higher charge
difference.

Fig 1.3: Schematic representation of the formation of an ionic bond between sodium and
chlorine.
Covalent bonds

In a covalent bond, one atom shares one or more pairs of electrons with another atom. This
sharing occurs because the atoms want to satisfy the octet rule, which means they aim to
have a full outer electron shell like noble gases. Covalent bonding usually happens between
two nonmetals and is the strongest and most common type of chemical bond found in living
organisms. Along with ionic bonds, covalent bonds are among the two most important types
of chemical bonds.

Covalent bonds can be classified into two types:

1. Nonpolar Covalent Bonds: Electrons are shared equally between the two atoms.
2. Polar Covalent Bonds: Electrons are shared unequally, resulting in a difference in
charge between the atoms.
 In polar covalent chemical bonding, electrons are shared unequally since the more
electronegative atom pulls the electron pair closer to itself and away from the less
electronegative atom. Water is an example of such a polar molecule. A difference in
charge arises in different areas of the atom due to the uneven spacing of the electrons
between the atoms. One end of the molecule tends to be partially positively charged,
and the other end tends to be partially negatively charged.

Examples:

 Two iodine atoms (I) combine to form iodine gas (I₂).

 One carbon atom (C) bonds with two oxygen atoms (O) to create a double covalent
bond in carbon dioxide (CO₂).

Module 1
Introduction to
Computational
Chemistry
Computational
Chemistry is the branch
of chemistry that uses
computer simulations to
solve chemical
problems. In
computational
chemistry, we use the
results of theoretical
chemistry incorporated
into computer
programs to calculate
the structure, energy
and other
properties of atoms and
molecules.
The term computational
chemistry is generally
used when a
mathematical model is
well
developed, that it can
be well automated for
implementation on a
computer. However, it
should
be remembered that
very few aspects of
chemistry can be
computed exactly,
rather, almost
every aspect of
chemistry has been
described in a
qualitative or
approximately
quantitative
computational scheme.
While computational
results normally
complement the
information obtained
from the
chemical experiments, it
can sometimes predict
some unobserved
chemical phenomenon.
It
can predict the reaction
mechanisms of new
reactions. This branch of
chemistry is also widely
used in the
pharmaceutical industry
today. Computational
chemistry is involved in
many stages
of drug development. It
allows scientists to
predict how specific
molecules will interact
with
each other These types of bonds normally occur between two
atoms of a non-metal.
Other examples of covalent compounds: Hydrogen (H2), Methane (CH4), Nitrogen (N2),
Water (H2O)
Fig 1.4: Schematic representation of sharing of electron in a covalent bond

Metallic Bond

Metallic bonds are the chemical bonds that hold atoms together in metals. They differ from
covalent and ionic bonds because the electrons in metallic bonding are delocalized, that is,
they are not shared between only two atoms. Instead, the electrons in metallic bonds float
freely through the lattice of metal nuclei. Metallic bonding accounts for many of the
distinctive properties of metals, such as their shiny appearance, malleability, and good
conductivity for heat and electricity.

For instance, while gallium primarily exhibits covalent bonding, it can form crystal structures
that are also stabilized by metallic bonds. Additionally, the mercurous ion demonstrates both
metallic and covalent bonding characteristics.

Several factors influence the strength of metallic bonds, including:

 The total number of delocalized electrons.

 The magnitude of the positive charge on the metal cation.
 The ionic radius of the cation.

Bond Parameters
Every bond has some characteristic feature associated with it, known as bond parameters.
Chemical bonds can be characterized with the help of these parameters. Some of the bond
parameters are:
1. Bond length
2. Bond angle
3. Dihedral angle
4. Bond enthalpy
5. Bond order
In this discussion, we will be looking at some of these parameters-bond length, bond
angle and the dihedral angle.
Bond length- bond length is the average distance between the nuclei of two bonded
atoms in a molecule.
Bond angle – bond angle is the angle between three atoms in a molecule.

Fig 1.5: Bond angles

Dihedral angle- the dihedral angle in a molecule is an angle between four atoms in a
molecule. Mathematically, it is the angle between two intersecting planes. In chemistry,
dihedral angle is the angle between two planes formed by two sets containing three
atoms each.

Fig.1.6: Dihedral angles

The dihedral angles in n-butane can be visualized easily with the help of sawhorse and
Newman projections.

Fig.1.7: Two methyl substitutes of syn n-butane are at a dihedral angle of 60 0.

 Representing molecules using Cartesian Co-ordinates.

For the purpose of carrying out calculations involving molecules in computational

chemistry, we need to represent the geometry of molecules in mathematical terms. One
of the simple ways to do this is to represent them by mentioning the cartesian
coordinates of each atom of the molecule.
In this method, three variables (x, y and z ) are used to describe the positions of each
atom in space. consider a hydrogen atom located at the centre of the Cartesian systems
we are able to show the position of the hydrogen atom in space and the simple way to
do this is to represent them by mentioning the Cartesian coordinates of each atom of
the molecule. In this method, three variables (x, y, z) are used to describe the position of
each atom in space. Consider a hydrogen atom located at the centre of the Cartesian
system we are using. We can represent as:

H 0.00000 0.00000 0.00000

For hydrogen molecule, the cartesian coordinates of both the hydrogen atom need to be
mentioned and can be represented as:

H 0.00000 0.00000 0.00000

H
H 0.00000 0.74000 0.00000

In the above case, one of the hydrogen atoms is located at the centre and the other one is
located at a distance of 0.74Å along the y-axis. Here, we have used a distance of 0.74 Å
between the hydrogen, as it is the typical bond length in a H2 molecule.

For a water molecule, we represent the coordinates of one oxygen atom and two
hydrogen atoms

O 0.00000 0.00000 0.00000

H 0.00000 0.27017 0.00000
H 0.00000 0.81192 0.73679

From the above table it may observed that there is no representation of bond length of
OH (0.975 Å ) or the H-O-H bond angle (104.5) in the coordinate table. So, It is not easy
to see the bond length or bond angle in the coordinate table even after much complex
mathematical calculations.
For more complex molecules we will need to mention the three cartesian variables of
each of the atoms. So, in a molecule with N atoms, a total of 3N variables or coordinates
will be used.

Representation using Z-matrix

Z matrix uses the approach of tracking the relative positions of atoms in a molecule
rather than keeping track of the absolute positions. In the Cartesian system, we mention
the absolute positions of atoms in space. However, a H2 molecule centred around (0,0,0)
is no different from the one centred around (0,0,1). Both of them will have the same
energy, dipole moment and other properties, provided they have equal bond length.
However the properties changes with the change in bond length, which is an internal
coordinate. In Z-matrix representation, we use these internal coordinates to represent a
molecule. So the Z-matrix is also known as the internal coordinate representation.
In the Z-matrix approach, we use the bond parameters such as the bond length, bond
angle and the dihedral angle to represent a molecule.
The hydrogen atom can be represented by just writing the symbol for hydrogen, H, as
there are no bond lengths or bond angles involved. So the Z-matrix representation of a
hydrogen atom is given below as:

Now consider the hydrogen molecule. The Z-matrix of H2 is given below

H
H 1 r2
Variables:
r2 =1.4352
Here there are two atoms of hydrogen involved. The first hydrogen atom is again
represented by the symbol H. The representation of the second hydrogen atom is given
in the second line of the Z-matrix, which read as,

H 1 r2

This should be interpreted as, the atom given in the second line (H) as at a distance of r2
from the atom given in the first line (the first hydrogen atom). The value of the r2
variable (in Angstrom, Å) is given below in the Z-matrix. In this case, r2=1.4352 Å. The
advantage of this representation is that the only parameter which has any effect on the
properties of the molecule, including energy, is the distance between the atoms (the
 bond length). We can change the value of the bond length by changing the value of r2
and calculate the value of r2 and calculate the energy of different bond lengths.
Now consider a bigger molecule, H2O molecule
Its Z-matrix can be given as :

O
H 1 r2
H 1 r3 2 a3
Variables:
r2 = 0.9700
r3 = 0.9700
a3 = 104.5
The first letters of the first three lines in the Z matrix represent the atoms present in the
molecule. From these, we understand that the molecule contains one oxygen atom and
two hydrogen atoms. The second line reads as,
H 1 r2
Just like the previous example, this line should be interpreted as the atom in the given
line (H) is at a distance of r2 from the atom in the first line (Oxygen).The value of r2 is
given as 0.9700A, which mean the O-H bond length is 0.97 A.

Now the third line is read as,

H 1 r3 2 a3

This should be interpreted as the atom in the given line (the second hydrogen atom) is at
a distance of r3 from the atom in the first line (O) and makes an angle a3 with the atom
in the first line (O) and the one in the second line (H). the values of r2 and a3 are given
as 0.97 Å and 104.5 degree, respectively. The values of r1, r2 and a3 represent the two
bond lengths and a bond angle in a water molecule.
Z matrix formulation in general can be :
The first three atom entries of a Z-matrix are different from the subsequent entries. The
first Z-matrix line declares a single atom. The second line of the Z-matrix input declares a
second atom, refers to the first atom and gives the distance between them. The third line
declares the third atom, refers to either the first or second atom, gives the distance between
them, refers to the remaining atom and gives the angle between them. All subsequent entries
begin with an atom declaration, a reference atom and a distance, a second reference atom and
an angle, a third reference atom and a dihedral angle. This can be summarized as:

1. First atom.
2. Second atom, reference atom, distance.
3. Third atom, reference atom A, distance between A and the third atom, reference atom
B, angle defined by atoms A, B and the third atom.
4. Fourth atom, reference atom A, distance, reference atom B, angle, reference atom C,
dihedral angle (A, B, C and the fourth atom).
5. All subsequent atoms follow the same basic form as (4)

Consider hydrogen peroxide. Z-matrix for this structure would be:

H
O 1 0.9
O 2 1.4 1 105.0
H 3 0.9 2 105.0 1 120.0

The first line of the Z-matrix simply specifies a hydrogen. The next line lists an oxygen atom
and specifies the internuclear distance between it and the hydrogen as 0.9 Angstroms. The
third line defines another oxygen with an O-O distance of 1.4 Angstroms (i.e., from atom 2,
the other oxygen) and having an O-O-H angle (with atoms 2 and 1) of 105 degrees. The
fourth and final line is the only one for which all three internal coordinates need be given. It
defines the other hydrogen as bonded to the second oxygen with an H-O distance of 0.9
Angstroms, an H-O-O angle of 105 degrees and a H-O-O-H dihedral angle of 120 degrees.

Z MATRIX of Formaldehyde

C
O 1 1.2
H 1 1.0 2 120
H 1 1.0 2 120 3 180
Z Matrix of Carbon dioxide -Dummy atoms

 Bond angles equal to 180 degrees are not accepted by Z matrix

 Therefore, we use dummy atoms when we are dealing with linear molecules
 When the input is read into a computer code, these dummy atoms are just used as
reference points, and do not enter into any calculation

O
C 1 1.2
X 1 1.0 2 90
O 1 1.2. 2 90 3 180

Z MATRIX of Ammonia

N1
H2 1 r2
H3 1 r3 2 a3
H4 1 r4 2 a4 3 d4

r2=1.0
r3=1.0
r4=1.0
a3=110.
a4=110.
d4=120.