Unit - I

Introduction to Computational Chemistry

Department of Chemistry & Biochemistry

1
SYLLABUS
Overview of Computational Chemistry, its scope and applications in [Link]. Introduction to modern
methods of Computational Chemistry; applications include Molecular geometry, Molecular
spectroscopy, Potential energy surfaces and Drug design. Models used in Computational Chemistry.
Principles of Molecular mechanics, Basics of Molecular Modelling and Simulations.

2
OVERVIEW
Computational Chemistry
Computational chemistry is a branch of chemistry that uses computer simulations to solve chemical
problems.

In computational chemistry, we use results of theoretical chemistry incorporated into computer programs
to calculate the structure, energy and other properties of atoms and molecules.
4
SCOPE
Scope of Computational Chemistry in
Computer Science
Computational chemistry sits at the intersection of chemistry and computer science, leveraging diverse
computational techniques and methodologies to solve complex molecular problems.

● Machine Learning & AI: Neural networks for property prediction, pattern recognition in molecular data,
and automated discovery of chemical relationships. Deep learning models for drug design and materials
science.
● Algorithms & Data Structures: Sophisticated algorithms for molecular dynamics, quantum mechanics
calculations, and optimization problems. Efficient data structures for handling large molecular datasets and
simulation trajectories.
● Data Science & Visualization: Big data analytics for chemical databases, statistical analysis of molecular
properties, and interactive 3D visualization of molecular structures and simulation results.
● Software Engineering: Development of specialized computational chemistry software, user interfaces for
molecular modeling, and integration of various computational tools into comprehensive platforms.
● High-Performance Computing: Parallel processing, distributed computing, and GPU acceleration for
computationally intensive molecular simulations. Cluster computing and cloud-based solutions for
large-scale calculations.
Scope of Computational Chemistry in various
fields
Computational chemistry is utilized across various
fields:
● Biochemistry: Understanding enzyme mechanisms
and protein folding.
● Materials Science: Designing new materials with
specific properties.
● Environmental Chemistry: Modeling pollutant
behavior and degradation.
Importance: These applications enable researchers to
predict outcomes and design experiments more
efficiently, reducing the need for extensive laboratory
work.
BASICS
Finding Molecular Geometry using Computational Chemistry

Molecular geometry refers to the 3D arrangement of atoms in a molecule, including:

● Bond lengths (distance between atoms)
● Bond angles (angle between two bonds)
● Dihedral/torsional angles (angle between two planes of atoms)
 Molecular Geometry using Geometry
Optimization
Geometry Optimization is the process of finding the most stable 3D arrangement of atoms in a
molecule where the potential energy is at a minimum. Its working can be explained in the
following steps.

● Start with an initial guess of molecular geometry.

● Compute energy and forces on atoms.
● Adjust atomic positions to reduce energy.
● The above two steps will be repeated many number
of times. With each step, there will be a change in
energy. The process will be repeated until we get an
energy minimum.
● The output will be optimized bond lengths, angles,
and dihedral angles corresponding minimum
potential energy.
Potential Energy Diagram
A Potential Energy Surface (PES) is a multidimensional surface that shows how the potential energy of a
molecular system changes with different arrangements of its atoms.
● Typically, the x- and y-axes represent nuclear coordinates
(bond lengths, angles, etc.), and the z-axis represents the
potential energy.
● The lowest point (valley) on the PES represents the most
stable structure or global minimum energy conformation.
● Saddle points on the PES correspond to transition states
— high-energy intermediates in chemical reactions

● PES helps in understanding reaction mechanisms,

transition states, and stable molecular conformations.
● Computational chemistry methods are used to calculate
energies for different molecular geometries and construct
the PES.
APPLICATIONS
Applications of Computational Chemistry
Computational chemistry is utilized across various fields:
● Predict Molecular Properties: Enables accurate prediction of physical and chemical properties
of molecules, eliminating the need for costly and time-consuming lab experiments.
● Elucidate Reaction Mechanisms: Provides detailed insights into reaction pathways, transition
states, and intermediate species, deepening our understanding of chemical transformations and
drug interactions.
● Design Novel Materials & Drugs: Facilitates the rational design of new molecules with desired
properties, from advanced materials to life-saving pharmaceuticals, by simulating their
behaviour before synthesis.
● Accelerate Research & Development: Significantly reduces the time and cost associated with
discovering and optimising new drugs and materials by simulating various scenarios.
1. Molecular Spectroscopy
Molecular spectroscopy is the study of the interaction of electromagnetic radiation with molecules to provide
information about their structure, energy levels, and dynamics.

Types of Spectroscopy
● Ultraviolet-Visible (UV-Vis) Spectroscopy –
studies the absorption of ultraviolet and visible
light, promoting electrons to higher electronic
energy levels.
● Infrared (IR) Spectroscopy – studies the
vibrational transitions of molecules by measuring
the absorption of infrared light.
● Microwave Spectroscopy – measures the
absorption of microwave radiation, which
induces rotational transitions in molecules.
Molecular Spectroscopy
A spectra shows a signal (a band or a peak) corresponding to the transitions occurring in a molecule.
Role of Computational Chemistry in
Molecular Spectroscopy

● Predicting Spectral Lines: Computational methods calculate energy levels of molecules and predict
frequencies or wavelengths at which absorption or emission occurs.
● Spectral Assignment: Helps interpret complex spectra by assigning peaks to specific molecular motions
or electronic transitions.
● Structural Verification: Comparing computed spectra with experimental spectra aids in verifying
molecular structure or identifying unknown compounds.
2. Computational Chemistry in Drug
Designing
● Understanding Molecular Interactions
○ Computational tools model how drug molecules interact with biological targets like proteins and DNA.
○ This helps predict binding affinity, interaction sites, and molecular conformation.
● Virtual Screening of Drug Candidates
○ Thousands of potential drug molecules can be screened virtually using docking and scoring
algorithms.
○ Reduces the need for costly and time-consuming lab experiments.
● Structure-Based Drug Design (SBDD)
○ Based on the 3D structure of a target (e.g., enzyme or receptor), computational models suggest
molecules that fit well into the active site.
○ Uses molecular docking and molecular dynamics simulations to refine drug candidates.
● Optimization of Drug-Like Properties
○ Simulations help predict and optimize properties such as solubility, permeability, stability, toxicity.
● Reduces Cost and Time
○ Significantly lowers the cost of early-stage drug discovery.
○ Speeds up identification of promising leads by reducing the number of physical experiments
required.
Molecular Modelling
Molecular Modelling is the use of computers to build and visualize molecules in 2D or 3D. It helps scientists
predict how molecules look, move, and react without without carrying out the experiments in lab.

Molecular Modelling involves:

● Building molecules: Atoms and bonds are arranged to form a molecule.
● Energy minimization: Molecules are adjusted to their most stable (lowest energy) shape.
● Simulations: Software can show how molecules move and change over time.

Applications of Molecular Modelling

● Drug Design: Finding the right molecule to fit a disease
target
● Material Science: Studying polymers, batteries, etc.
● Biology: Understanding protein structures
● Chemistry: Predicting reactions and properties
METHODS
Methods in Computational Chemistry
Methods in computational chemistry refer to the mathematical models and algorithms used to
simulate, calculate, and predict the structure, properties, reactions, and behavior of molecules and
materials using computers.
Methods can be broadly classified as:
● Quantum Mechanical Methods: Solve the Schrödinger equation to describe electrons and
atomic interactions. Examples: Ab initio methods, Density Functional Theory (DFT),
Semi-empirical methods.
● Classical Methods: Use Newtonian mechanics and force fields to model molecular structures.
Examples: AMBER, CHARMM, OPLS force fields.
● Hybrid Methods: Combine quantum mechanics for reactive sites with classical mechanics for
the rest of the system.
Examples of Methods in Computational
Chemistry
Some modern methods used in computational chemistry are:
● Ab Initio Methods
● Density Functional Theory (DFT)
● Semi-Empirical Methods
● Molecular Mechanics (MM)
Ab Initio Methods
● Based purely on quantum mechanics without using any empirical data.
● Solve the Schrödinger equation numerically for molecular systems.
● Most accurate but computationally expensive.

Density Functional Theory (DFT)

● Based on electron density rather than the many-body wavefunction.
● Offers a good balance between accuracy and speed.
● Widely used for molecular geometry, reaction energies, and materials modeling.
● Less computational cost than ab initio methods, but more accurate than semi-empirical.
Semi-Empirical Methods
● Based on quantum mechanical foundations but use experimental data to simplify calculations.
● Faster but less accurate than ab initio or DFT.
● Suitable for large molecules where full quantum treatment is impractical.

Molecular Mechanics (MM)

● Treats atoms as balls and springs based on classical physics (Newton's laws).
● Uses force fields to describe bonded and non-bonded interactions (e.g., bond stretching, van der
Waals forces).
● Ideal for biomolecules, proteins, polymers, and long timescale simulations.
Molecular Mechanics
Molecular Mechanics is a classical approach used to model the structure, energy, and dynamics of molecules
by treating atoms as balls and bonds as springs. Instead of solving quantum equations, it applies Newton’s
laws of motion to simulate molecular systems.

● Molecular mechanics focuses on bond lengths, bond angles, torsional angles and non-bonded
interactions (like van der Waals and electrostatic forces).

● The total potential energy of a molecule is calculated using a force field formula.
Etotal=Ebonds+Eangles+Etorsions+Enon-bonded
Where, Ebonds , Eangles , Etorsions and Enon-bonded represent the bond stretching energy, angle bending energy,
torsional (dihedral) angle energy and the non-bonded interaction energy, respectively.

● Force fields are predefined sets of parameters and equations used to describe molecules.
Advantages of Molecular Mechanics
● Very fast and efficient
● Can handle large systems like proteins, polymers, and DNA
● Suitable for geometry optimization, molecular dynamics, and docking studies

Disadvantages of Molecular Mechanics

● Does not account for electrons
● Cannot model bond breaking/forming (no reaction mechanisms)
● Limited accuracy for transition states and charged species

Applications
● Protein structure modeling
● Drug design
● Material property prediction
MODELLING
Representing a molecule using Z-matrix
● A Z-matrix provides a description of each atom in a molecule in terms of its atomic number,
bond length, bond angle and the dihedral angle.
● Z-matrix representation is also called the internal coordinate representation.

Hydrogen atom
Z-matrix:

The hydrogen atom can be represented by just writing

the symbol for hydrogen, ‘H’, as there are no bond
lengths or bond angles involved (no internal
coordinates).
Hydrogen molecule

Z-matrix:

H
H 1 r2
Variables:
r2= 1.4352

● The first hydrogen atom is represented by the symbol ‘H’ in the first line.
● The representation of the second hydrogen atom is given in the second line of the Z-matrix. This
should be interpreted as, the atom given in the second line (H) is at a distance of r2 (r2= 1.4352
Å) from the atom given in the first line (the first hydrogen atom).
 Z-matrix:
Water molecule
O
H 1 r2
H 1 r3 2 a3
Variables:
r2= 0.9700
r3= 0.9700
a3= 109.47
● The oxygen atom is represented by the symbol ‘O’ in the first line.
● The representation of the first hydrogen atom is given in the second line of the Z-matrix. This should be
interpreted as, the atom given in the second line (H) is at a distance of r2 (r2= 0.9700 Å) from the atom given in
the first line (the oxygen atom).

● The second hydrogen atom is represented in the third line. This should be interpreted as, the atom in the given
line (the second hydrogen atom) is at a distance of r3 (r3= 0.9700 Å) from the atom in the first line (O) and makes
and angle a3 (a3= 109.47o) with the atom in the first line (O) and the one in the second line(H).
 Hydrogen peroxide (H2O2) molecule

● The first oxygen atom is represented by the symbol ‘O’ in the first line.

● The representation of the second oxygen atom is given in the second line of the Z-matrix. This
should be interpreted as, the atom given in the second line (O) is at a distance of r2 (r2=
Z-matrix: 1.4740 Å) from the atom given in the first line (the first oxygen atom).
O
O 1 r2 ● The first hydrogen atom is represented in the third line. This should be interpreted as, the atom
in the given line (the first hydrogen atom) is at a distance of r3 (r3= 0.9500 Å) from the atom in
H 1 r3 2 a3
the first line (O) and makes and angle a3 (a3= 94.8o) with the atom in the first line (O) and the
H 2 r4 1 a4 3 d4
one in the second line(H).
Variables:
r2= 1.4740 ● The H in the fourth line represents the second hydrogen atom. This line tells us that the
r3= 0.9500 second hydrogen atom is at a distance of r4 (r4= 0.9500 Å) from the atom in the second line
a3= 94.8 (the second oxygen atom), make an angle of a4 (a4= 94.8o) with the atoms in the second and

r4= 0.9500 the first lines and finally makes a dihedral angle, d4 (d4= 115.5o) with the atoms in the first line
(O), the second line (O) and the third line (H).
a4= 94.8
d4= 115.5
Z-matrix in computational chemistry
● Z-matrix is used to represent the internal coordinates of molecules and used as input in
calculations.

A sample input file used for calculations. The Z-matrix of the

molecule can be seen.
Constructing Z-matrix of molecule using Avogadro
Avogadro is a free and open source software used for designing and creating molecular structures.

Avogadro can be downloaded and installed for

free from [Link]
Drawing molecules and constructing Z-matrix of simple
molecules

In a similar manner, we’ll construct the molecular structures and the z-matrixes of H2O2, HCHO,
NH3 and CO2.
Drawing molecules and constructing Z-matrix of simple
molecules

Once the Z-matrix is saved, we can open it using any text editor like ‘Notepad’ or ‘Gedit’. Further, we
can alter the Z-matrix by changing the bond lengths or bond angles (the values of r1, a1, d1 etc.) to
our liking.
 REFERENCES
● Jensen, Frank. (2017). Introduction to computational chemistry. (2nd ed., Vol. 4). Wiley
● Dawber, J.G., Moore, A.T. (1980). Atomic structure and chemical bonding. In: Chemistry for the
Life Sciences. Palgrave, London.
● A. K. Haghi, L. P. (2021). Applied Chemistry and Chemical Engineering (2nd ed., Vol. 2). CRC
Press.

35