CONTENTS
Subject |
PAGE |
Chapter one. Introduction and Literature Review …………………………………… |
1 |
1.1 Kinetics of H2S + O3 Reaction ……………………………………………………… |
2 |
1.2 Dynamics of C2H4 + S Reaction ……………………………………………………. |
4 |
Chapter two. Theory and Methodology ……………………………………………… |
7 |
2-1 Potential Energy Surface (PES) ……………………………………………………. |
8 |
2-1-1 The Born-Oppenheimer Approximation ………………………………………… |
9 |
2-1-2 Methods to Calculate the Energy ………………………………………………… |
12 |
2-1-2-1 T1 Diagnostic ………………………………………………………………… |
16 |
2-2 Dynamics …………………………………………………………………………… |
17 |
2-2-1 Interfacing the dynamics calculations to the potential energy surface ………….. |
18 |
2-2-1-1 Analytic Potential Energy Surfaces …………………………………………… |
18 |
2-2-1-2 Straight Direct Dynamics …………………………………………………….. |
19 |
2-2-1-3 Interpolated Variational Transition State Theory (IVTST) …………………… |
19 |
2-2-1-4 Dual-Level Dynamics ………………………………………………………… |
20 |
2-3 Potential Energy Surfaces by Interpolation ………………………………………… |
21 |
2-3-1 Molecular Coordinates …………………………………………………………… |
21 |
2-3-2 Local Shape of the PES ………………………………………………………….. |
23 |
2-4 Minimum Energy Path ……………………………………………………………… |
25 |
2-5 Collision Theory ……………………………………………………………………. |
27 |
2-6 Trajectory Calculation ……………………………………………………………… |
31 |
2-7 The Reaction Probability, Reaction Cross-Section and Rate Constant Calculations in Quasiclassical Trajectory Method (QCT) …………………………………………. |
32 |
2-8 Rate Theory for Simple Barrier Reaction .…………………………………………. |
34 |
2-8-1 Conventional Transition State Theory …………………………………………… |
34 |
2-8-2 Variational Transition State Theory ……………………………………………… |
37 |
2-8-3 Canonical Variational Transition State Theory …………………………………. |
40 |
2-9 Validation Against Accurate Quantum Mechanical Dynamics ……………………. |
42 |
2-9-1 Tunneling and the Transmission coefficient …………………………………….. |
42 |
2-9-2 Multidimensional Tunneling Corrections Based on the Adiabatic Approximation ………………………………………………………………………….. |
45 |
2-10 Unimolecular Reaction Theory …………………………………………………… |
52 |
2-10-1 The Hinshelwood Revision of Lindemann Mechanism ……………………….. |
53 |
2-10-2 The RRKM Theory …………………………………………………………….. |
55 |
2-10-3 Application of the RRKM Theory for Unimolecular Dissociation/EliminationReactions – The Strong Collision Model ……………………………………………… |
57 |
2-10-4 Comparison of RRKM and Transition State Theories …………………………. |
59 |
2-11 Chemically Activated Reactions ………………………………………………….. |
60 |
2-12 Statistical Mechanics ……………………………………………………………… |
63 |
2-13 Molecular Dynamics and Monte Carlo Simulations ……………………………… |
68 |
2-13-1 Molecular Dynamics …………………………………………………………… |
69 |
2-13-2 Monte Carlo Simulations ………………………………………………………. |
72 |
Chapter Three. Results and Discussion: Multichannel RRKM-TST and Direct-Dynamics CVT Study of the Reaction of Hydrogen Sulfide with Ozone ……………. |
74 |
3-1 Introduction ………………………………………………………………………….. |
75 |
3-2 Computational Method ……………………………………………………………… |
76 |
3-3 Geometries and Potential Energy Surfaces …………………………………………. |
84 |
3-4 Calculation Methods of the Rate Constants ………………………………………… |
96 |
3-4-1 RRKM-SSA method for Rw1 to R9 ……………………………………………… |
96 |
3-4-2 Canonical Variational Transition State Theory Calculations …………………….. |
105 |
3-4-2-1 Computational Details …………………………………………………………. |
105 |
3-5 Conclusion …………………………………………………………………………… |
123 |
Chapter Four. Results and Discussion: A Theoretical Study on the Dynamics of the C2H4+ S Reaction on an Interpolated Potential Energy Surface ……………………. |
125 |
4.1 Introduction ………………………………………………………………………….. |
126 |
4.2 Methods and computational details ………………………………………………….. |
128 |
4.2.1 Constructing the Potential Energy Surface ……………………………………….. |
128 |
4.2.2 Confident weights ………………………………………………………………… |
131 |
4.2.3 Iterative development of the PES …………………………………………………. |
133 |
4.2.4 Ab initio calculations ……………………………………………………………… |
135 |
4.3 Classical Dynamics ………………………………………………………………….. |
144 |
4.4 The Iteration Procedure ……………………………………………………………… |
144 |
4.5 Calculation of thermal rate constants ………………………………………………… |
145 |
4.5.1 Rate Constant of Reaction R1 …………………………………………………….. |
146 |
4.5.2 Rate Constant of Reaction R2 …………………………………………………….. |
148 |
4.6 An ab initio investigation of spin-allowed and spin-forbidden pathways of thegas phase reactions of S (3P) + C2H4 …………………………………………………. |
154 |
4.7 Canonical Variational Transition-State Theory Calculations for Dissociation ofC2H4S (Reaction R1‘) ………………………………………………………………….. |
158 |
4.8 Conclusion …………………………………………………………………………… |
162 |
References ……………………………………………………………………………….. |
164 |
LIST OF TABLES
Table |
Page |
3.1 Relative Energies of Various Species at Different Levels of Theory in kJ mol-1 on the Singlet Surface ………………………………………………………………………. |
93 |
3.2 Harmonic Vibrational Wave Numbers (cm-1) and Moments of Inertia (amu Å2) for Various Species at the B3LYP/6-311+G(d,p) Level ……………………………………. |
94 |
3.3 Relative Energies of Various Species at Different Levels of Theory in kJ mol-1 on the Triplet Surface ……………………………………………………………………….. |
95 |
3.4 Microcanonical Variational RRKM Results for Unimolecular Dissociation Reaction R4′ ……………………………………………………………………………………….. |
103 |
3.5 Microcanonical Variational RRKM Results for Unimolecular Dissociation Reaction R8 ……………………………………………………………………………………….. |
104 |
3.6 Microcanonical Variational RRKM Results for Unimolecular Dissociation Reaction R14 ………………………………………………………………………………………. |
109 |
3.7 Microcanonical Variational RRKM Results for Unimolecular Dissociation Reaction R16 ………………………………………………………………………………………. |
110 |
3.8 Rate Constants of the R12 Reaction (in s-1) at mpwb1k method ……………………. |
113 |
3.9 Rate Constants of the R13 Reaction (in s-1) at mpwb1k method ……………………. |
113 |
3.10 Rate Constants of the R15 Reaction (in s-1) at mp2 method ……………………….. |
114 |
4.1 Relative Energies of Various Species on Triplet Surface of Reaction Sulfur and Ethylene at Different Levels of Theory in kJ mol-1 ……………………………………… |
139 |
4.2 Relative Energies of Various Species at Different Levels of Theory in kJ mol-1 on the singlet surface reaction sulfur and ethylene …………………………………………. |
141 |
4.3 Summary of the trajectory calculations for reaction of sulfur and ethylene on triplet surface …………………………………………………………………………………… |
160 |
4.4 Summary of the trajectory calculations for reaction of sulfur and ethylene on singlet surface …………………………………………………………………………………… |
161 |
LIST OF FIGURES
Figure |
Page |
2.1 Illustration of how the steepest descent algorithm follows a path that oscillates around the minimum energy path ……………………………………………………….. |
27 |
2.2 Contour plot of a model bimolecular reaction that indicates the possible tunneling paths ……………………………………………………………………………………… |
51 |
2.3 Rate coefficient of reaction (2.84) at T=400 K as a function of pressure …………… |
54 |
2.4 Schematic diagram of the Lindemann-Hinshelwood mechanism …………………… |
55 |
2.5 Schematic energy diagram of a chemically activated reaction ……………………… |
61 |
2.6 Fraction of molecules that will be found at various energies above the ground- stateenergy for two different temperatures …………………………………………………… |
68 |
3.1 IRC for reaction R1 at the B3LYP level …………………………………………….. |
77 |
3.2 IRC for reaction R2 at the B3LYP level …………………………………………….. |
77 |
3.3 IRC for reaction R3 at the B3LYP level …………………………………………….. |
78 |
3.4 IRC for reaction R4 at the B3LYP level …………………………………………….. |
78 |
3.5 IRC for reaction R5 at the B3LYP level …………………………………………….. |
79 |
3.6 IRC for reaction R6 at the B3LYP level …………………………………………….. |
79 |
3.7 IRC for reaction R7 at the B3LYP level …………………………………………….. |
80 |
3.8 IRC for reaction R9 at the B3LYP level …………………………………………….. |
80 |
3.9 IRC for reaction R10 at the B3LYP level …………………………………………… |
81 |
3.10 IRC for reaction R11 at the B3LYP level ………………………………………….. |
81 |
3.11 IRC for reaction R12 at the B3LYP level ………………………………………….. |
82 |
3.12 IRC for reaction R13 at the B3LYP level ………………………………………….. |
82 |
3.13 IRC for reaction R15 at the MP2 level ……………………………………………… |
83 |
3.14 Optimized geometries of the stationary points at the MP2/Aug-cc-pVTZ level …… |
89 |
3.15 Schematic of the PES of the reaction of H2S + 1O3 at the CCSD(T)/Aug-cc-pVTZ level. The energies are corrected for ZPEs ……………………………………………… |
91 |
3.16 Schematic of the PES of the reaction of H2S + 3O3 at the CCSD(T)/Aug-cc-pVTZ level. The energies are corrected for ZPEs ……………………………………………… |
92 |
3.17 Arrhenius plot of the calculated rate constants for various channels of H2S+O3 reaction at 760 Torr pressure of N2 ……………………………………………………… |
102 |
3.18 (a) Variation of frequencies along the MEP and (b) the vibrationally adiabatic ground-state potential , MEP , and ZPE for reaction R10 at the MPWB1K level ……………………………………………………………………………………… |
115 |
3.19 (a) Variation of frequencies along the MEP and (b) the vibrationally adiabatic ground-state potential , MEP , and ZPE for reaction R12 at the MPWB1K level ……………………………………………………………………………………… |
116 |
3.20 (a) Variation of frequencies along the MEP and (b) the vibrationally adiabatic ground-state potential , MEP , and ZPE for reaction R13 at the MPWB1K level ……………………………………………………………………………………… |
117 |
3.21 (a) Variation of frequencies along the MEP and (b) the vibrationally adiabatic ground-state potential , MEP , and ZPE for reaction R15 at the MPWB1K level ……………………………………………………………………………………… |
118 |
3.22 Thermally averaged transmission probability at 298 K for R10 reaction ………….. |
119 |
3.23 Thermally averaged transmission probability at 298 K for R12 reaction ………….. |
119 |
3.24 Thermally averaged transmission probability at 298 K for R13 reaction ………….. |
120 |
3.25 Thermally averaged transmission probability at 298 K for R15 reaction ………….. |
120 |
3.26 Arrhenius plot for reactions R10 and R11 at the CCSD(T)/Aug-cc-pVTZ level ….. |
121 |
3.27 Arrhenius plot at the CCSD(T)/Aug-cc-pVTZ level for channels R12-R16′ ……… |
122 |
4.1 A schematic representation of the twenty one bond lengths defining the C2H4S system, with atoms labeled as shown ……………………………………………………. |
130 |
4.2 Reaction probablity for reaction R1, as a function of the PES data set size ………… |
134 |
4.3 Reaction probablity for reaction R2, as a function of the PES data set size ………… |
135 |
4.4 Optimized geometries of the stationary points on triplet surface of reaction sulfur and ethylene at the B3LYP/cc-pVDZ level ……………………………………………… |
137 |
4.5 Relative energies of the stationary points in kJ mol-1 on triplet surface of reaction sulfur and ethylene at the CCSD(T)/aug-cc-pVTZ level ………………………………… |
138 |
4.6 IRC for reaction R1 at the B3LYP level …………………………………………….. |
140 |
4.7 IRC for reaction R1‘ at the B3LYP level …………………………………………….. |
140 |
4.8 Optimized geometries of the stationary points at the B3LYP/cc-pVDZ level in Singlet Surface of reaction sulfur and ethylene …………………………………………. |
142 |
4.9 Relative Energies of the Stationary Points in kJ mol-1 on Singlet Surface of ReactionSulfur and Ethylene at the CCSD(T)/cc-pVDZ level …………………………………… |
143 |
4.10 Reactive Cross Section as a function of relative kinetic energy for reaction R1 for2000 trajectory ………………………………………………………………………….. |
150 |
4.11 Dependence on temperature of rate constant for the reaction R1 ………………… |
151 |
4.12 Reactive Cross Section as a function of relative kinetic energy for reaction R2 for2000 trajectory ………………………………………………………………………….. |
152 |
4.13 Dependence on temperature of rate constant for the reaction R2 ………………….. |
153 |
4.14 The singlet and triplet potential energy profile of C2H4 + S system rigorous CCSD(T) level including ZPE …………………………………………………………… |
156 |
4.15 IRC for the singlet and triplet surfaces of C2H4 + S system, Calculated at B3lyp/cc-pVDZ level ……………………………………………………………………. |
157 |
4.16 Arrhenius plot at the CCSD(T)/Aug-cc-pVTZ level for channels R1′ …………….. |
159 |
TABLE OF SCHEMES
Scheme |
Page |
Scheme 1. Compelete Mechanism of The H2S + O3 Reaction ………………………….. |
85 |
Scheme 2. Calculated Rate coefficients from RRKM-SSA method ……………………. |
98 |