LAB MANUAL

CHY1005: Introduction to Computational Chemistry

NAME : NARESH SUTHAR

REGISTRATION NUMBER : 24BAI10260

FACULTY NAME : DR THAMIM M.

SLOT NUMBER : A24 + C21 + F22

CLASS NUMBER : AB -103

BRANCH : CSE Sp. In AI & ML

- Date:28/11/2024

Index

S. Title of the Experiment Page

No. No.
1. Molecular Visualization 2

2. Single Point Energy Calculation of Molecules 7

3. Energy Minimization of Diatomic Molecules 9

4. Geometry Optimization of Polyatomic Molecules 11

5. Generation and Visualization of Molecular Orbitals 13

6. Potential Energy Scan of Molecules 16

7. Frequency Calculation of Molecules 18

8. Computational Thermochemistry 21

9. Visualization of Molecules In VMD 23

1
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry

Introduction to Computational Chemistry
Name: Naresh Suthar Branch: cse(aiml)

Registration Number: 24BAI10260 Duration: 90 min.

Practical 1: Molecular Visualization

Aim
To build and visualize molecules using Avogadro.

Theory:
There are many tools to visualize the molecule. In this session, we will be using Avogadro, an
opensource software. Generally, computational chemistry allows chemists to analyze properties of
molecules using a computer. The computer can predict quantities such as: molecular geometry (bond
length, bond angle), molecular orbitals, bonding, spectroscopy, and energy. It can also simulate the
motion of molecules (molecular dynamics). The first step for these calculations is building the
molecule. Building a molecule in computer involves describe the molecule to the computer and tell
it what you want it to calculate. If you wanted to calculate the bond length in carbon dioxide you
would specify an initial guess for the (x,y,z) coordinates of the three atoms. Then you would add
some keywords to tell the computer to find the optimum geometry. This can be accomplished using
one of the dedicated software. The generated structures can be saved in either of the following two
formats:
1. Cartesian format
2. Z-matrix
Cartesian Geometry of Water:

2
- Date:28/11/2024

Z-matrix of Geometry of Benzene:

Procedure:
• Molecules can be built using the Draw Tool of Avogadro GUI window.

• Open Avogadro, select the Draw Tool.

• Select the Draw tool (at the left end of the toolbar; looks like a pencil) to start building a molecule.
In the left pane, make sure that the Adjust Hydrogens checkbox is checked.

• Select carbon from the dropdown menu. Click once in the drawing window (the black field). To
create a bond to a second carbon, click the first carbon and drag away from it.

• If you want to delete an atom, you can right-click it while in drawing mode. If you want to add
another type of atom, you can select it from the dropdown menu (select “Other…” and choose any
element from the periodic table)

• Single click puts an isolated atom, eventually with its hydrogen atoms. If you want to create a
bond, drag until the bond is long enough and release the mouse button.

• Bond order can be changed through the Bond Order drop down menu, or by typing the numbers
"1", "2", or "3". Bonds can also be added to a molecule by left clicking on a bond that has already
been created (this process cycles through single, double, and triple bonds). Right clicking on a
bond will delete the bond, and the atom it's bonded to.

• Click "File" --> "Save" to save the molecule for later use.

3
- Date:28/11/2024

Results:
Draw these molecules and paste here the output as a snapshot of the molecule and XYZ file.
1. Ethylene
Grey : carbon
White : hydrogen

2. Acetaldehyde
Red : oxygen , grey : carbon , white : hydrogen

4
- Date:28/11/2024

3. Glucose
Red : oxygen
Grey : carbon
White : hydrogen

4. Alanine
Red : oxygen , grey : carbon , blue : nitrogen , white : hydrogen

5
- Date:28/11/2024

5. 2-heptene
Grey : carbon
White : hydrgen

6
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry

Introduction to Computational Chemistry
Name: Naresh Suthar Branch: cse(aiml)

Registration Number: 24BAI10260 Duration: 90 min.

Practical 2: Single Point Energy Calculation of Molecules

Aim
To calculate single point energy of molecules using a quantum mechanical method.

Theory:
Single point energy is the potential energy of a molecule for a given arrangement of the atoms in the
molecule. Single point energies are the lowest energy solution for the Schrödinger equation and are
the simplest properties one might aim to obtain. It consists in the calculation of the wave function
and energy for a given system with a well-specified geometric structure, i.e. at a single, fixed point
on the potential energy surface. The computed energy is the total energy, sum of the electronic energy
and nuclear repulsion energy. To perform a single point energy calculation a well-defined level of
calculation must be specified. A level of calculation is uniquely defined by the combination of a
theoretical method with a basis set.
Many aspects of the system’s behaviour can be understood when the total energy of the system is
known, or more often, when the energy difference between two or more electronic or nuclear
configurations is known as a function of some parameters.

Procedure:
• Build a water molecule using ChemCompute server or any other molecular builder.
• Prepare the input for single point energy calculation. Choose an appropriate method and basis set.
• Submit the job by clicking on submit button.
• Note down the electronic energy of the water molecule.
• To perform calculations on a new molecule click the Submit Tab at the top.

Results:

7
- Date:28/11/2024

Draw the following molecules and report their single point energy.

S. No. Molecule Single Point Energy (Hartree)

1 H2O -76.27183 a.u.

2 He2 -5.7 a.u.

3 NH3 -56.511 a.u.

4 CH4 -40.48329 a.u.

5 CO -115.64713 a.u.

8
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry Introduction to

Computational Chemistry
Name: Naresh Suthar Branch: cse(aiml)

Registration Number: 24BAI10260 Duration: 90 min.

Practical 3: Energy Minimization of Diatomic Molecules Aim

To perform energy minimization of diatomic molecules using a quantum mechanical method.

Theory:
The energy minimization is a way to get the stable structure of molecule which is the main aim of
computational chemistry. The resultant molecular structure is refereed as “Optimized geometry”.
Here, we mainly do the minimization of energy as the function of nuclear coordinate.
Let us take the Energy (E) is a function of cartesian coordinate (‘x,y,z’) of a given molecular structure
i.e;, E (xyz), and represented as E(Q). As the geometry changes its energy also varies. We change the
geometry of the molecule is such a way that energy gets minimum. This will be the most stable
structure, and known as “Optimized geometry”.
Mathematically,
We do the following operations:
𝜕𝐸(𝑄)
=0
𝜕𝑄
And,
𝜕2𝐸(𝑄)
>0
𝜕𝑄

Procedure:
• Build the oxygen molecule using ChemCompute server or any other molecular builder.
• Prepare the input for energy minimization calculation. Choose an appropriate method and basis
set.
• Submit the job by clicking on submit button.

9
- Date:28/11/2024

• Note down the optimized energy of the molecule.

• Report the cartesian optimized geometry of oxygen molecule.

Results:
Draw these molecules and optimize it. Report the optimized energy.

S. No. Molecule Optimized Energy Bond Distance (Å)

(Hartree) Initial Optimized
1 O2 -150.14543 a.u. 0.144 nm 0.144 nm
2 He2 -5.7 a.u. 0.3 nm 0.3 nm
3 Cl2 -920.21861 a.u. 0.202 nm 0.202 nm
4 Br2 -5148.05345 a.u. 0.23 nm 0.23 nm
5 I2 -13840.65655 a.u. 0.266 nm 0.266 nm

10
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry Introduction to

Computational Chemistry
Name : Naresh Suthar Branch:cse(aiml)
Registration Number: 24BAI10260 Duration: 90 min.

Practical 4: Geometry Optimization of Polyatomic Molecules

Aim
To perform energy minimization of diatomic molecules using a quantum mechanical method.

Theory:
The energy minimization is a way to get the stable structure of molecule which is the main aim of
computational chemistry. The resultant molecular structure is refereed as “Optimized geometry”.
Here, we mainly do the minimization of energy as the function of nuclear coordinate.
Let us take the Energy (E) is a function of cartesian coordinate (‘x,y,z’) of a given molecular structure
i.e;, E (xyz), and represented as E(Q). As the geometry changes its energy also varies. We change the
geometry of the molecule is such a way that energy gets minimum. This will be the most stable
structure, and known as “Optimized geometry”.
Mathematically,
We do the following operations:
𝜕𝐸(𝑄)
=0
𝜕𝑄
And,
𝜕2𝐸(𝑄)
>0
𝜕𝑄

Procedure:
• Build a H2O water molecule using ChemCompute server or any other molecular builder.
• Prepare the input for energy minimization calculation. Choose an appropriate method and basis
set.
• Submit the job by clicking on submit button.
• Note down the optimized energy of the molecule.

11
- Date:28/11/2024

• Report the cartesian optimized geometry of molecule.

•
Results:
Draw these molecules and optimize it. Report the optimized energy.

S. No. Molecule Optimized Energy (Hartree)

1 H2O -76.37183 a.u.

2 Ethyelene -78.52925 a.u

3 Glucose -686.76526 a.u.

4 Methane -40.4824 a.u.

5 NH3 -56.51141 a.u.

12
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry

Introduction to Computational Chemistry
Name : Naresh Suthar Branch: cse(aiml)
Registration Number: 24BAI10260 Duration: 90 min.

Practical 5: Generation And Visualization of Molecular Orbitals

Aim
To generate and visualize the molecular orbitals using computational chemistry tools.

Theory:
The molecular orbital is a fundamental concept in Quantum chemistry. In the Schrodinger picture,
the MO’s represents the wavefunction of molecules which stores all information pertaining to the
molecule. It is important to visualize the MO’s to get the bonding feature of the molecule. Generally,
these MO’s are linear combination of the Hydrogenic wavefinction “s, p, d, f” orbitals, known as
atomic orbitals. In this session, we will generate the atomic orbitals of Ne atom. In Quantum
Chemistry,

Ψ = ∑ 𝐶𝑖Φi

Where Φ is molecular orbitals. In case of atoms these are atomic orbitals.

In this lab, we will generate and visualize the Mos of the Neon atom. The electronic configuration of Ne
atom is 1s2, 2s2, 2px2, 2py2, 2pz2. All these orbitals will be visualized here.

Procedure: • Build a “Carbon” atom in Chemcompute or any other molecular

builder.

• Prepare the input for single point energy calculation. Choose an appropriate method and basis
set.

• Visualize and report the energy of its different orbitals.

Results:

13
- Date:28/11/2024

Draw the following molecules and generate their atomic/molecular orbitals.

• C

Figures: 1s, 2s, 2px, 2py, 2pz orbitals

1s,2s.

2px,2py,2pz.

14
- Date:28/11/2024

• O2

Figures: Highest Occupied Molecular Orbital (HOMO), Lowest Unoccupied Molecular Orbital (LUMO)

• Ethylene

Figures: Highest Occupied Molecular Orbital (HOMO), Lowest Unoccupied Molecular Orbital (LUMO)

15
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry

Introduction to Computational Chemistry
Name : Naresh Suthar Branch: cse(aiml)
Registration Number: 24BAI10260 Duration: 90 min.

Practical 6: Potential Energy Scan of Molecules Aim

To scan the potential energy scan of a molecule using a quantum mechanical method.

Theory:
The potential energy curve (PEC) or surface is the plot of electronic energy as a variation of nuclear
geometry. The origin of PEC is due to heavy mass of nuclei, the electron move much faster than
nuclei, as a result we may separate the electronic and nuclei motion – Born-Oppenheimer
Approximation. This is the fundamental theory in the area of molecular dynamics. The PEC is also
considered as a tool for visualization of chemical reaction. It is assumed that molecule undergoes to
chemical reaction by moving over the PEC.

Figure: General Potential energy surface and Potential energy curve

In this session, we will generate the PEC of O2 molecule by varying its bond distance.

16
- Date:28/11/2024

Procedure:
• Build O2 the atom in ChemCompute or any other tool.
• Prepare the input for single point energy calculation. Choose an appropriate method and basis set.
• Change the bond distance of O2 molecule at intervals of 0.1 Å and repeat the calculation.
• Repeat the procedure for 15 steps by changing the distance by 0.2 Å in each step.
• Record the Electronic Energy (above the molecule box) and the bond distance

Results:
Draw O2 molecule and generate its PECs.

S. No. Bond Distance (Å) Electronic Energy (Hartree) Relative Energy (Hartree)
1 0.80 -0.3244 0000
2 1.00 -0.1792 0.1452
3 1.20 -0.1480 0.1764
4 1.40 -0.1604 0.1640
5 1.60 -0.1857 0.1387
6 1.80 -0.2113 0.1131
7 2.00 -0.2332 0.0912
8 2.20 -0.2504 0.0740
9 2.40 -0.2633 0.0611
10 2.60 -0.2727 0.0517
11 2.80 -0.2796 0.0448
12 3.00 -0.2844 0.0400
13 3.20 -0.2879 0.0365
14 3.40 -0.2903 0.0341
15 3.60 -0.2920 0.0324

Figure: Relative energy vs bond distance

17
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry

Introduction to Computational Chemistry
Name : Naresh Suthar Branch: cse(aiml)
Registration Number: 24BAI10260 Duration: 90 min.

Practical 7: Frequency calculation of Molecules Aim

To perform frequency calculation of a molecule by computational chemistry tools.

Theory:
Frequency calculation is a crucial step in molecular modeling, which helps in predicting the
vibrational spectrum of a molecule. This involves calculating the frequencies of various modes of
vibrations of a molecule and their corresponding energies. The vibrations of a molecule can be
divided into three types: stretching, bending, and torsional. The stretching mode is associated with
the stretching of bonds between atoms in the molecule, the bending mode with the bending of bonds,
and the torsional mode with the rotation of atoms around a bond.

To get the vibrational structure of the molecule, Hessian calculation is performed. Hessian is the
second order energy derivative with respect to the nuclear coordinate. It represents the curvature of
the energy surfaces. The eigen value of Hessian matrix leads to the vibrational frequencies and
eigenvector leads to the vibrational mode of the molecules. The Hessian gives important information
of the molecules whether it is minimum or not? If the all eigen-values is positive, then the molecule
is present at one of the minima. Therefore, to confirm the minima structure of molecule, Hessian
calculation is important. Additionally, after getting the vibrations of the molecule, Infrared spectra
can also be simulated.

In this session, we generate the vibrational frequencies, vibrational mode and IR spectra of the water
molecule.

18
- Date:28/11/2024

Procedure:
• Build the water molecule in Chemcompute or any other software.

• Perform energy optimization of the water molecule at B3LYP/6-31G level of theory.

• Set up the calculation for the Frequency calculation using the optimized geometry.
• Submit the job.
• Under the Vibrationals visualizations, choose the vibrational frequencies to animate the three
vibrational mode.

Results:
Draw the following molecules and generate its Vibrational frequencies. Report the four vibrational
frequencies.
• CO2

19
- Date:28/11/2024

• CH4

20
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry Introduction to

Computational Chemistry
Name : Naresh Suthar Branch: cse(aiml)
Registration Number: 24BAI10260 Duration: 90 min.

Practical 8: Computational Thermochemistry Aim

To calculate reaction enthalpy for a reaction computational chemistry tools.

Theory:
Chemical reactions involve changes in energy as well as matter. Enthalpy change of a reaction is one
of the most useful thermodynamic quantities. The experimental determination of the reaction using
a standard bomb calorimeter requires great attention to detail to achieve the desired precision and
accuracy. The tools of quantum chemistry can also be used to calculate thermodynamic properties,
such as enthalpy (ΔH), entropy (ΔS), and Gibbs free energy (ΔG) of reactions.

In this laboratory, you will determine the ΔH associated with the Haber-Bosch process:
N2 + 3 H2 ⇌ 2 NH3

Procedure: • Build the N2 molecule in Chemcompute or any other

software.

• Perform energy optimization of the water molecule at B3LYP/6-31G level of theory.

• Set up the calculation for the Frequency and thermochemistry calculation using the optimized
geometry.
• Repeat the optimization for the molecules involved in the Haber-Bosch process at the same level
of theory.
• Calculate the reaction enthalpy (ΔH).

Results:
Calculate (ΔH) for the following reactions:

21
- Date:28/11/2024

• N2 + 3 H2 ⇌ 2 NH3
N2 = -109.3426403817 = -109.3426403817

H2 = 3 (-1.1673385838)= - 3.5020157514

NH3 = 2x(-56.5113596852)= -113.0227193704

ΔH= ΔHp -ΔHr

ΔHrxn = (-113.0227193704) – [(-109.3426403817)+ (- 3.5020157514)

= (-113.0227193704)- (-112.8446561331)

• CH4 + H2O ⇌ CO + 3 H2

ΔH of CH4 = -40.14400691

ΔH of H2O = -75.97965784

ΔH of CO = -112.72637063

ΔH of 3 H2 = 3x(-1.11178571) = -3.33535713

ΔH= ΔHp -ΔHr

ΔHrxn -40.14400691)+( -75.97965784)] – [(-112.72637063)+( -3.33535713)]

= [(

= (-116.12366475)-( -116.06172776)

22
- Date:28/11/2024

School of Advanced Sciences and Languages, Chemistry Introduction to

Computational Chemistry
Name: Naresh Suthar Branch: cse(aiml)

Registration Number: 24BAI10260 Duration: 90 min.

Practical 9: Visualization of Molecules in VMD

Aim
The aim of this lab is to learn how to use VMD (Visual Molecular Dynamics) software to visualize
molecules and analyze their properties.
Theory:
VMD is a powerful visualization tool for studying and analyzing molecules, biomolecules, and materials.
It provides a wide range of capabilities for visualizing structures, trajectories, and data from molecular
simulations. VMD is designed to be highly extensible and customizable, with support for a wide range of
molecular file formats and data types.

In this lab, you will learn how to use VMD to visualize and analyze molecules. You should now be
able to load, manipulate, and analyze various types of molecules and generate visualizations and
plots of their properties.

23
- Date:

Procedure:
• Installing VMD
o Download and install VMD from the official website:
https://www.ks.uiuc.edu/Research/vmd/
o Follow the installation instructions for your operating system.
• Loading a molecule o Open VMD and select "New Molecule" from the File menu.
o Choose the appropriate file format for your molecule (e.g. PDB, XYZ, etc.) and select the
corresponding file.
o VMD will automatically load the molecule and display it in the graphics window. o
Manipulating the molecule o Use the mouse to rotate, translate, and zoom in/out of the
molecule in the graphics window.
o Use the "Display" menu to change the representation of the molecule (e.g. ball-and-stick,
space-filling, etc.). o Use the "Color" menu to change the color of the molecule based on
various properties (e.g. atom type, charge, etc.).
o Use the "Graphics" menu to adjust various display settings (e.g. lighting, shadows, etc.).

• Analyzing the molecule o Use the "Measure" menu to calculate various properties of the
molecule (e.g. bond length, angle, dihedral, etc.).
o Use the "Graph" menu to generate plots of various properties (e.g. RMSD, energy, etc.)
over time.

24
- Date:28/11/2024

o Use the "Extensions" menu to access additional plugins and tools for analyzing and
visualizing molecular data.
o Saving the molecule
o Select "Save Molecule" from the File menu to save the current state of the molecule
(including all display and analysis settings) to a VMD file.
o Select "Write PDB" from the File menu to save the current state of the molecule in PDB
format (with optional additional data fields).

Figure: Steps to create high quality rendering.

25
- Date:

Results:
Load “Ubiquitin” protein (PDB ID: 1UBQ) in VMD. Submit rendered geometry of the molecule.

26