ISSN: 2637-4609
Adeboye Omolara*
Received: February 02, 2018; Published: February 21, 2018
Corresponding author: Adeboye Omolara, Chemistry Department, Emmanuel Alayande College of Education, P. M. B 1010, Nigeria
DOI: 10.32474/AOICS.2018.01.000123
Computational chemistry is a technique used only by people that have experience in using tools that are difficult to understand and apply. This reviewed article give insight to simple procedure that can be used to optimize geometry of molecules to obtain geometric parameters such as bond length, bond angles and dihedral angles and also to carry out kinetic and thermo dynamic calculations on uni- molecular reactions involving a first order and bimolecular reactions involving second order reactions.
Keywords: Procedures; Optimization; Kinetics; Spartan; Thermodynamics
Today, advances in software have produced programs that are easily used by chemists. Along with new software comes new literature on the subject. Different books are available to describe the fundamental principles of computational chemistry at almost every level of detail. A number of books also exist that explain how to apply computational chemistry techniques to simple calculations appropriate for student assignments [1]. However, there is dearth of information on detailed and simple progressive procedure for calculations using Spartan software, hence this review.
Conformational search was done using molecular mechanics force field (MMFF) to obtain the most stable conformer. The most stable structure is the structure that has the lowest energy. Draw the structure. From Calculation: Calculate click Conformer distribution from the menu list with Molecular Mechanics MMFF. Click Global calculations and Click submit
Reaction path study is important to obtain a clear insight into the true and correct transition state so as to obtain correct kinetic and thermodynamic values. This was carried out by selecting reaction coordinates for the molecule or the reaction, here the internal coordinates would be varied from its initial bond distance in the stable reactant form to its bond distance in the product form values in the product form and set calculation into desired iterations (2030 steps of calculations). Set the transition State Structure select the reaction coordinate. Set the bond distance in the reactant form and set the bond distance in the product form. Select Energy Profile in calculate. Click Submit. Open the result and click to animate the different Structures and the corresponding Energies of formation. The first structure has the Lowest Energy and it is the most Stable Structure.
Geometry optimization is very important in Computational Chemistry. It involves taking the geometric approximations and making them as exact as possible. It involves a series of computing process or repeated mathematical calculations performed on the molecule until the energy of the molecule reach a minimum. Three types of calculations can be performed:
a) Molecular Mechanics
b) Ab-Initio
c) Semi-Empirical Calculation
These Calculations are used to optimize the geometry and this depends basically on the type of level of calculations and the basis set used. The geometry of a molecule determines the physical and chemical properties so to have accurate molecular geometry is very important. Three measurements that are very important in calculating the molecular geometry are:
i Bond Angles (degrees)
ii. Bond Distance (Angstrom)
iii. Dihedral Angles (degrees)
Molecular geometry is simply the molecule's arrangement of atoms in space. This arrangement and the electrons around the atom determine the energy level of that molecule. The energy of a molecular system varies even when there is little change in its structure. The main aim of geometry optimization is to determine the atomic arrangement that makes a molecule stable. Molecules are most stable and most likely to be found in nature at very low energy. In order to optimize a molecular geometry, one would want to test various possibilities to see which will give the lowest energy value. This is achieved by creating a Potential Energy Surface (PES). PES is the mathematical expression that shows the relationship between the different molecular geometries and their corresponding single - point energies.
This is performed to obtain the correct molecular structure with the lowest energy value.
Click new from file on the worksheet of spartan. build a molecule of ethyl acetate for example using appropriate hybridization in the model kit. choose geometry and click constrain dihedral
Click On the C, O, C, C Atoms (Figure 1)
Figure 1:
Click on the pink constrain marker that appears on the C-O Bond
Choose display click properties
Click on the dynamic box in the constrain properties window. Change the values to read 180.000 to 0.000 and change steps to 20. Hit the enter key
Choose Set up, click calculation select energy profile at the ground state with semi - empirical pm3 or any other levels of calculations.
Click submit and save as ethylacetate spartan.
Click ok to started calculation and click ok to completed calculation.
Click yes to open all the possible molecular structures.
Choose display. Click spreadsheet click on the box with grey shade click add select rel. energy.
Click ok.
To add angles to the spreadsheet,
Choose geometry click measure dihedral.
Click on the constrain marker on the c - o bond.
Click on the yellow p box on the right of the screen Choose display click plot
Choose X axis click dihedral and Y axis click e (kj/mol).
Click Ok.
Press play to view the animation of the bond rotation.
Save and Close The File.
Click New from file to display the worksheet
Use the model kit to build a molecule using appropriate hybridizations to build ethyl acetate (Figure 2).
Figure 2:
Choose the structure with the lowest minimization energy value.
To minimize the energy click on the icon E the energy value is shown at the lower part of the workspace.
Choose set up to calculate the equilibrium geometry at the ground state with semi empirical PM3 or other calculation methods.
Start from current geometry click orbitals and energies and charges and bond order on print.
Set Options OPTCYCLE=200-5000 in case the calculation ran out of cycle.
Click Converge
Click Submit
Save as ethylacetate. Spartan in a folder.
Click ok to start the calculation Click ok to completed calculation Click display
Select output to read the energy of formation and other parameters
Choose display to display the output, properties, orbiral energies (homo-lumo), spectra, formulas.
Click on the icon ? Measure distance to determine the bond length.
Click on the icon <? Measure angleto determine the bond angle
Click on the icon'?'Measure dihedral
Record the values for the optimized structure.
Click New from file.
Build a molecule using the appropriate hybridization from the model kit.
Click on E to minimize the Energy of the structure
Record the minimization energy.
Click View to display the Structure on the Screen.
Choose geometry and click constrain bond distance from the predicted transition state bond distance to the normal bond distance (path studies).
Click the dynamic box to set the values.
Hit the enter key
Click set up to choose calculation
Select energy profile at the ground state with semi empirical pm3
Click Submit and Save.
Click ok to start calculation
Click ok to completed calculation
Click yes to open all the possible molecular structures.
Click ok
Press Play to view the animation of the bond breaking and bond formation.
Build the structure with proper orientation and mechanism using the model kit.
Minimize the energy using the icon E.
Click on the icon to guess transition state and show the suggested mechanism of the reaction.
Click on a double arrow head at the lower right corner of the worksheet to display the Transition State Structure and possibly to confirm the mechanism.
Click set up to calculate the transition state geometry at the ground state with semi empirical PM3 and other calculation methods.
Select compute ir, orbital energies, thermodynamics, vibrational modes, charges and bond order on print.
Click options optcycle =250
Click converge
Click submit
Click ok to started calculation
Click ok to completed calculation
Select output to read the heat of formation and other parameters [2-4]
The kinetic and thermodynamic parameters were obtained using the sum of the ground state energy (GSE) and the statistical mechanically calculated enthalpy to arrive at a closer approximation of the true energy of the molecule. The enthalpy of a species will be defined as:
H_{i} = GSE_{i} + H_{i}^{sm}
Where the superscript 'sm' is the statistical mechanically calculated enthalpy. Substitute this into the initial definition of the heat of reaction we have:
ΔH_{rxn} = (GSE_{product} + H^{sm} _{product})- (GSE_{reactant} + H^{sm} _{reactnt})
Activation energy (Ea) was calculated according to the transition state theory for a unimolecular reaction at 698K [2,5].
E_{a} = ΔH * - RT
The entropy of the reaction was calculated by taking the difference of product and reactant entropies that is
ΔS_{reaction} = S_{Product} - S_{reactants}
And ΔS_{activated}=S_{transition}- S_{reactant}
The Gibbs free energy and Gibbs free energy of activation was calculated using the modified version of the heat of reaction equation. Knowing that
G = H - TS and ΔG = ΔH * - TΔS
The first order coefficient k(T) was calculated using transition state theory (TST) [6] assuming that the transmission coefficient is unity as shown in the following equation
Where ΔG^* is the Gibbs free energy change between the reactant and the transition state and K^'and h are the Boltzmann and Planck constants respectively.
Pre -exponential factor is given as
Arrhenius rate was obtained using the rate equation [7]
The bond lengthening character was used to describe the extent of bond lengthening to show the rate determining step and to explain the concertedness and synchronicity of the reaction.
This detailed and simple progressive procedure for calculations described using Spartan software has been used for computational modeling of kinetic and thermodynamic [9-16] studies of various reactions.
Bio chemistry
University of Texas Medical Branch, USADepartment of Criminal Justice
Liberty University, USADepartment of Psychiatry
University of Kentucky, USADepartment of Medicine
Gally International Biomedical Research & Consulting LLC, USADepartment of Urbanisation and Agricultural
Montreal university, USAOral & Maxillofacial Pathology
New York University, USAGastroenterology and Hepatology
University of Alabama, UKDepartment of Medicine
Universities of Bradford, UKOncology
Circulogene Theranostics, EnglandRadiation Chemistry
National University of Mexico, USAAnalytical Chemistry
Wentworth Institute of Technology, USAMinimally Invasive Surgery
Mercer University school of Medicine, USAPediatric Dentistry
University of Athens , GreeceThe annual scholar awards from Lupine Publishers honor a selected number Read More...
The annual scholar awards from Lupine Publishers honor a selected number read more...