View source for Rep:MOD:hly3614

== '''
== '''Introduction''' ==
=== Potential Energy Surface (PES) ===
A potential energy surface (PES) is an energy 'landscape' that interpretates the total energy of a molecule as a function of molecular geometries based on Born-Oppenheimer Approximation. According to the fact that a polyatomic molecule often has 3N-6 ( N= Number of atoms and N>2 ) geometric degrees of freedom, the PES is likely to possess 3N-6 dimensions but is usually reduced to 3-D representation.

=== A Minimum And A Transition State ===
A minimum on the PES corresponds to a stationary point where the gradient of the energy curve (i.e the first derivative) is zero and the curvature (i.e the second derivative) reflects that any change being made to this point leads to an increase of energy, whereas a transition state is a saddle point at which the change of total energy with respect to the change of geometry is zero but in terms of curvature, moving away from the point, the energy drops in one direction. 
=== Relevant Gaussian Calculations ===
A geometry optimization in Gaussian helps find out the stationary points (i.e gradient=0) by calculating the first derivatives of the energy curve and the frequency calculations are used to obtain information about the curvature of the points (i.e the second derivative) by calculating harmonic frequencies of the polyatomic molecule, confirming whether a minimum or a transition state is reached, based on the fact that all the vibrational frequencies for a minimum are real numbers whereas an imaginary frequency occurs for a transition structure.

[[User:Nf710|Nf710]] ([[User talk:Nf710|talk]]) 23:19, 14 December 2016 (UTC) That is a nice way of interpret it, but the correct way is to say that you taking the second derivatives 

== '''Exercise 1: Reaction of Butadiene with Ethylene''' ==
=== Reaction Scheme ===
[[File: EX1 REACTION SCHEME.png|x250px|400px|]]
=== Frontier Molecular Orbital Analysis ===
[[File:Finallllll ex1 mo.png|thumb|left|x600px|600px|Fig. 1: MO Diagram of Butadiene/Ethylene TS]] 
[[File:Cycloaddition_lecture_handout.PNG|thumb|left|x450px|450px|Fig. 2: MO Diagram for an example of cycloaddtion with normal electron demand <sup>[1]</sup>]] 
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Molecule and orbital:
! style="background: #0D4F8B; color: white;" | MOs:
! style="background: #0D4F8B; color: white;" | Symmetry [gerade/ungerade 
(g/u for short)]:
|-
| Butadiene, HOMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 16; mo 11; rotate y 120; mo nodots nomesh fill translucent; mo titleformat ""; </script>
 <uploadedFileContents>NOV29 HLY EX1 PREOPTIMIZED FRAGMENT.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| g
|-
| Butadiene, LUMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 16; mo 12; rotate y 120; mo nodots nomesh fill translucent; mo titleformat ""; </script>
 <uploadedFileContents>NOV29 HLY EX1 PREOPTIMIZED FRAGMENT.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| u
|-
| Ethene, HOMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 14; mo 6; rotate y 90; mo nodots nomesh fill translucent; mo titleformat ""; </script>
 <uploadedFileContents>NOV_29_ETHENE_MINIMUM_OPT_PM6.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| u
|-
| Ethene, LUMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 14; mo 7; rotate y 90; mo nodots nomesh fill translucent; mo titleformat ""; </script>
 <uploadedFileContents>NOV_29_ETHENE_MINIMUM_OPT_PM6.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| g
|-
| Butadiene/ethylene Transition state, HOMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 16; mo 17; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.01 </script>
 <uploadedFileContents>NOV30 HLY EX1 TS OPT PM6.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| u
|-
| Butadiene/ethylene Transition state,  LUMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 16; mo 18; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.02 </script>
 <uploadedFileContents>NOV30 HLY EX1 TS OPT PM6.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| u
|-
| Butadiene/ethylene Transition state, 2nd HOMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 16; mo 16; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.02 </script>
 <uploadedFileContents>NOV30 HLY EX1 TS OPT PM6.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| g
|-
| Butadiene/ethylene Transition state, 2nd LUMO
| <jmol>
 <jmolApplet>
 <color>white</color>
 <size>200</size>
 <script>frame 16; mo 19; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.02 </script>
 <uploadedFileContents>NOV30 HLY EX1 TS OPT PM6.LOG</uploadedFileContents>
 </jmolApplet>
</jmol>
| g
|}

According to the results of Gaussian calculations, the HOMO of the transition state is formed from the LUMO of the butadiene and the HOMO of the ethylene, both of which have the π<sub>u</sub> symmetry. Theoretically, however, as illustrated in Fig 2, a Diels-Alder Reaction with normal electron demand  favours the interaction between the HOMO of diene and the LUMO of the dienophile due to a smaller energy gap and better orbital overlap. In the case of Fig.2, the small energy gap results from attaching an electron withdrawing group to the dienophile substrate that lowers its LUMO energy. This can be further enhanced by making the diene substrate more electron-rich, leading to a higher energy HOMO. The inconsistency of the computed results is likely to arise from the fact that the unsubstituted butadiene and ethylene are not a good diene and not a good dienophile respectively to ensure the normal energy match. As for symmetry, interacting orbitals also need to have the same symmetry to allow a non-zero orbital overlap, so the g-g and u-u interactions are symmetry allowed whereas a g-u interaction is symmetry disallowed because the orbital overlap integral is equal to zero.

=== C-C Bond Length Analysis ===
'''Changes in Bond Lengths:'''
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Molecule:
! style="background: #0D4F8B; color: white;" | Illustration:
! style="background: #0D4F8B; color: white;" | C-C Bond Lengths/Å:
|-
| Butadiene
|[[File:NUMBERER_BUTADIENE.png|x100px|100px|]]
| C<sub>1</sub>-C<sub>2</sub>=1.33531

C<sub>2</sub>-C<sub>3</sub>=1.46835

C<sub>3</sub>-C<sub>4</sub>=1.33531
|-
| Ethene
|[[File:Ethene.png|x80px|80px|]]
|C<sub>5</sub>-C<sub>6</sub>=1.32731
|-
| Transition State
|[[File:NUMBERED TRANSITION STATE.png|x100px|100px|]]
| C<sub>1</sub>-C<sub>2</sub>=1.37977
C<sub>2</sub>-C<sub>3</sub>=1.41110

C<sub>3</sub>-C<sub>4</sub>=1.37979

C<sub>4</sub>-C<sub>6</sub>=2.11469

C<sub>6</sub>-C<sub>5</sub>=1.38177

C<sub>5</sub>-C<sub>1</sub>=2.11479
|-
| Product
|[[File:Product 1.png|x100px|100px|]]
| C<sub>1</sub>-C<sub>2</sub>=1.50034

C<sub>2</sub>-C<sub>3</sub>=1.33766

C<sub>3</sub>-C<sub>4</sub>=1.50034

C<sub>4</sub>-C<sub>6</sub>=1.54004

C<sub>6</sub>-C<sub>5</sub>=1.54076

C<sub>5</sub>-C<sub>1</sub>=1.54003
|-
|}
<br />
[The bond lengths are measured by extracting the bond lengths of corresponding optimized structures in GaussView.]
<br />
'''Representative C-C bond lengths with respect to hybridization/Å:'''
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | sp3-sp3 
! style="background: #0D4F8B; color: white;" | sp3-sp2
! style="background: #0D4F8B; color: white;" | sp2-sp2
|-
|literature value<sup>[2]</sup>: 1.54Å
computed value: 1.54Å
|literature value<sup>[2]</sup>: 1.50Å
computed value: 1.50Å
|literature value<sup>[2]</sup>: 1.47Å
computed value: 1.47Å
(except C<sub>2</sub>=C<sub>3</sub> in the product = 1.34Å)
|-
|}
<br />
By comparing the computed bond lengths with the literature values, the C=C double bonds (C<sub>1</sub>-C<sub>2</sub>, C<sub>3</sub>-C<sub>4</sub>, C<sub>5</sub>-C<sub>6</sub>) in both of the reactants start to elongate until eventually become single bonds in the product during the reaction, while the previous C-C single bond in butadiene (C<sub>2</sub>-C<sub>3</sub>) is shortend to form a double bond although the bond length is slightly shorter than expected. This might be due to the ring strain leading to a distortion of the bond. Also, based on an appropriate trajectory of the reactants, two new single bonds (C<sub>4</sub>-C<sub>6</sub> and C<sub>5</sub>-C<sub>1</sub>) are formed with typical sp3 C-C bond lengths of 1.54Å via a transition state where the distances between the butadiene terminal carbons and the ethylene carbons are appox. 2.11Å. A typical Van der Waals' radius of a carbon is ~1.70Å <sup>[3]</sup>, and by looking at the IRC profile, the two reactants start to approach each other from a distance of ~3.41Å, which is approximately twice of the Van der Waals' radius. The computed distance of 2.11Å at transition state lying in between the Van der Waals' distance of ~3.41Å and the final bond lengths of 1.54Å indicates a partially formed C-C single bond at the transition state.

=== Vibration at Transition State Corresponding to Reaction Path ===
[[File: Gif_vibration_at_transition_state_ex1.gif|1500x400px|gif 1: Imaginary vibrational frequency at transition state corresponding to reaction path: -948.57cm<sup>-1</sup>]]
<br />
''Gif 1'': Vibration at transition state with an imaginary vibrational frequency: -948.57cm<sup>-1</sup>
<br />
<br />
[[File: 145.10.gif|1500x400px|gif 2: first positive vibrational frequency at transition state: 145.10cm<sup>-1</sup>]]
<br />
''Gif 2'': vibration with first positive vibrational frequency: 145.10cm<sup>-1</sup>
<br />
<br />
<br />
The vibration occurring at -948.57cm<sup>-1</sup> in the transition state is illustrated in ''Gif.1''. The negative sign actually represents an imaginary vibrational frequency, which results from solving the square root of a negative force constant, i.e the second derivative (curvature) of the energy curve related to this particular reaction pathway. Since mathematically, a negative second derivative indicates a maximum point, in this case, the negative value occurring at -948.57cm<sup>-1</sup> suggests the vibration of bonds corresponding to an energy maximum in the transition state. 
<br />
<br />
Comparing the vibrations shown in Gif.1 and Gif.2, the vibration occurring at -948.57cm<sup>-1</sup> shows a trend towards bond formation between the butadiene terminal carbons and ethylene carbons, whereas the lowest positive vibration mode indicates the vibrations of the two reactant molecules are rather in two parallel planes without any interactions to allow a bond formation.
<br />
<br />
In addition, the illustration of the vibration at -948.57cm<sup>-1</sup> suggests that the formation of the 2 bonds is synchronous, but an IRC animation makes the idea of this concerted reaction more evident:
<br />
[[File: DEC1 HLY ex1 irc.gif|1500x375px|gif 3: IRC Animation of the Diels-Alder reaction of butadiene/ethylene]]
<br />
''gif 3'': IRC Animation of the Diels-Alder reaction of butadiene/ethylene
<br />

CLever use of IRC to identity this. Good section

== '''Exercise 2: Reaction of Cyclohexadiene and 1,3-Dioxole''' ==
=== Reaction Scheme ===
[[File:EX2 REACTION SCHEME.png|600px|Reaction Scheme of cycloaddtion of cyclohexadiene and 1,3-Dioxole]]

<br />

=== Frotier Molecular Orbital ( FMO ) Analysis ===
<br />
'''FMOs of the reactants:'''
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Molecules and Orbitals
! style="background: #0D4F8B; color: white;" | Illustration
! style="background: #0D4F8B; color: white;" | Symmetry (g/u)
|-
|1,3-Dioxole HOMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 12; mo 19; rotate y 30; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.02</script>
 <uploadedFileContents>DEC1 HLY ex2 dioxole OPT 631G.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|u
|-
|1,3-Dioxole LUMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 12; mo 20; rotate y 30; mo nodots nomesh fill translucent; mo titleformat "";  mo cutoff 0.02</script>
 <uploadedFileContents>DEC1 HLY ex2 dioxole OPT 631G.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|g
|-
|Cyclohexadiene HOMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 18; mo 22; rotate y 30; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.02</script>
 <uploadedFileContents>DEC1 HLY EX2 CYCLOHEXADIENE 631G OPT.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|(distorted)g
|-
|Cyclohexadiene LUMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 18; mo 23; rotate y 30; mo nodots nomesh fill translucent; mo titleformat "";  mo cutoff 0.02</script>
 <uploadedFileContents>DEC1 HLY EX2 CYCLOHEXADIENE 631G OPT.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|(distorted)u
|}
<br />
<br />
'''Illustration of FMOs during the cyclohexadiene/1,3-Dioxole transition state:'''
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Types of Transition States and FMOs
! style="background: #0D4F8B; color: white;" | Illustration [optimized at B3LYP/6-31G(d) level]
! style="background: #0D4F8B; color: white;" | Symmetry (g/u)
|-
|ENDO HOMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 34; mo 41; rotate y 30; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.01</script>
 <uploadedFileContents>NOV 30 HLY EX2 ENDO TS OPT 631G.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|u
|-
|ENDO LUMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 34; mo 42; rotate y 30; mo nodots nomesh fill translucent; mo titleformat "";  mo cutoff 0.01</script>
 <uploadedFileContents>NOV 30 HLY EX2 ENDO TS OPT 631G.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|u
|-
|EXO HOMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 18; mo 41; rotate y 30; mo nodots nomesh fill translucent; mo titleformat ""; mo cutoff 0.01</script>
 <uploadedFileContents>NOV 30 HLY EX2 EXO TS OPT 631G.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|u
|-
|EXO LUMO
|<jmol>
 <jmolApplet>
 <color>white</color>
 <size>250</size>
 <script>frame 18; mo 42; rotate y 30; mo nodots nomesh fill translucent; mo titleformat "";  mo cutoff 0.01</script>
 <uploadedFileContents>NOV 30 HLY EX2 EXO TS OPT 631G.log</uploadedFileContents>
 </jmolApplet>
</jmol>
|u
|}
<br />
According to the FMOs illustrated above, the HOMOs and LUMOs of both ENDO and EXO transition states are dominated by the interactions between the HOMO of 1,3-Dioxole and the LUMO of cyclohexadiene, both of which have a π<sub>u</sub> symmetry. Therefore ENDO and EXO transition states also have π<sub>u</sub> symmetries. This computed result refers to an inverse electron demand, which is likely due to the fact that being attached to two electron donating groups (i.e the O atoms ) makes 1,3-dioxole a very poor dienophile with a high energy LUMO and thus making the LUMO of dienophile and HOMO of diene too far in energy to interact. 
<br />

=== Analysis on The Reaction Barrier ΔG† and The Reaction Energy ΔG ===
'''Sum of Electronic and Thermal Free Energies at each stage during the reaction'''
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Molecules
! style="background: #0D4F8B; color: white;" | Sum of Electronic and Thermal Free Energies/ Hartrees/Particle
! style="background: #0D4F8B; color: white;" | Sum of Electronic and Thermal Free Energies ×627.5095×4.184 / KJ/mol
|-
|1,3-Dioxole
| -267.068132
| -701187.3133
|-
|Cyclohexadiene
| -233.324375
| -612593.0878
|-
|ENDO TS
| -500.332149
| -1313621.931
|-
|EXO TS
| -500.329168
| -1313614.105
|-
|ENDO PROD
| -500.418691
| -1313849.147
|-
|EXO PROD
| -500.417322
| -1313845.553
|-
|}
*'''Unit conversion:'''
The way to convert unit from ''Hartree/particle'' to ''KJ/mol'' is by multiplying constants 627.5095 and 4.184 <sup>[4]</sup>
<br />
*'''Calculating Reaction Barrier ΔG†:'''
ΔG†= [G<sub>TS</sub> - (G<sub>cyclohexadiene</sub>+ G<sub>dioxole</sub>)]    '''<eq.a>'''
<br />
*'''Calculating Reaction Energy ΔG:'''
ΔG= [G<sub>Product</sub> - (G<sub>cyclohexadiene</sub>+ G<sub>dioxole</sub>)]  '''<eq.b>'''
<br />
<br />
Therefore, based on the computed results and the calculation methods listed above:
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | 
! style="background: #0D4F8B; color: white;" | For ENDO
! style="background: #0D4F8B; color: white;" | For EXO
|-
|ΔG†/KJ/mol
| +158.47
| +166.30
|-
|ΔG/KJ/mol
| -68.75
| -65.15
|-
|}
<br />
The ENDO product turns out to be both kinetically and thermodynamically more favored as it can be reached via the lowest transition barrier among the three, and can lead to the lowest energy product as well. 
<br />
'''Secondary Orbital Interactions'''
<br />
Having all of the pi systems aggregated close in space, the ENDO transition state is expected to be more sterically clashed than the EXO TS and high in energy. However, the lower reaction barrier via an ENDO transition state can be explained by the secondary orbital interactions since a strong, in-phase orbital overlap is easily observed from the HOMO of the ENDO transition state. The dominant secondary orbital interaction between the p orbitals on O atoms and C atoms stablizes the HOMO of the ENDO TS, lowering the activation barrier and yielding the kinetically more stable ENDO product formation. 
<br />

[[User:Nf710|Nf710]] ([[User talk:Nf710|talk]]) 23:58, 14 December 2016 (UTC) Your energies are correct with the correct reasoning, try not to us u and g for symmetry. as the group theory behind those symmetry labels requires a center of inversion. 

== Exercise 3: Diels-Alder vs Cheletropic ==
=== Reaction Scheme ===
[[File: REACTION SCHEME HLY EX3.png|x500px|600px]]

=== Optimized Transition Structure at PM6 Level ===
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Reaction Pathways
! style="background: #0D4F8B; color: white;" | Optimized TS 
|-
|EXO Cycloaddition
|[[File: Ex3_exo_ts_opt.gif|1200x400px|thumb|left|Vibrational frequency at TS = -351.77cm<sup>-1</sup>]]
|-
|ENDO Cycloaddition
|[[File: Ex3_endo_ts_opt.gif|1200x400px|thumb|left|Vibrational frequency at TS = -334.39cm<sup>-1</sup>]]
|-
|Chelatropic Reaction
|[[File: Chelatropic ex3 ts opt.gif|1000x400px|thumb|left|Vibrational frequency at TS = -486.71 cm<sup>-1</sup>]]
|-
|}

<br />

=== Visualising reactin coordinates: IRC of The Three Possible Pathways ===
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Reaction Pathways
! style="background: #0D4F8B; color: white;" | IRC 
|-
|EXO Cycloaddition
|[[File: DEC1 HLY ex3 EXO IRC.gif|1200x400px]]
|-
|ENDO Cycloaddition
|[[File: DEC1 HLY EX3 ENDO IRC.gif|1200x400px]]
|-
|Chelatropic Reaction
|[[File: Chelatropic hly dec1 IRC.gif|1000x400px]]
|}
<br />
<br />
The animation of IRCs illustrates that during all pathways, as the SO<sub>2</sub> approaches close enough to the xylylene, the pi electrons in the 6-membered ring in xylylene starts to conjugate across the system and eventually the 6-membered ring becomes aromatic.
<br />
Besides, the IRC also indicates a concerted reaction for the chelatropic pathway, whereas the formation of C-S and C-O in both of the cycloaddition pathways are asynchronous. The reason why the formation of C-S follows the C-O bond formation is that since S atom is likely to have a larger Van der Waals' radius than O, it approaches the xylylene from a further distance than O.

=== Calculating the Reaction Barriers ΔG† and the Reaction Energies ΔG ===
'''Sum of Electronic and Thermal Free Energies at each stage during the reaction'''
<br />
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | Molecules
! style="background: #0D4F8B; color: white;" | Sum of Electronic and Thermal Free Energies/ Hartrees/Particle
|-
|SO<sub>2</sub>
| -0.119267
|-
|xylylene
| 0.178545
|-
|ENDO TS
| 0.090562
|-
|EXO TS
| 0.092078
|-
|Chelatropic TS
| 0.099062
|-
|ENDO PROD
| 0.021703
|-
|EXO PROD
| 0.021452
|-
|Chelatropic PROD
| 0.000001
|-
|}
<br />
Therefore, based on the computed results and the calculation methods listed in Exercise 2, eq.a and eq.b:
<br />
'''Calculating Reaction Barrier ΔG†:'''
'''ΔG†= [G<sub>TS</sub> - (G<sub>cyclohexadiene</sub>+ G<sub>dioxole</sub>)]×627.5095×4.184 / KJ/mol'''    '''<eq.a>'''
<br />
'''Calculating Reaction Energy ΔG:'''
'''ΔG= [G<sub>Product</sub> - (G<sub>cyclohexadiene</sub>+ G<sub>dioxole</sub>)]×627.5095×4.184 / KJ/mol'''  '''<eq.b>'''
<br />
the ΔG† and ΔG are calculated:
{| class="wikitable"
! style="background: #0D4F8B; color: white;" | 
! style="background: #0D4F8B; color: white;" | For ENDO
! style="background: #0D4F8B; color: white;" | For EXO
! style="background: #0D4F8B; color: white;" | For Chelatropic
|-
|ΔG†/KJ/mol
| +82.14
| +86.12
| +104.45
|-
|ΔG/KJ/mol
| -98.65
| -99.31
| -155.63
|-
|}

Based on the calculations, the ENDO Diels-Alder Reaction is kinetically preferred as it has the lowest reaction barrier among the three pathways, while the Chelatropic product is thermodynamically most stable due to a much lower energy level of product than the others. 
<br />
<br />
=== Reaction Profiles for the three possible pathways ===
[[File:Ex3 reaction profile hly dec2.png|600px|Reaction Profiles for the three possible pathways]]

(Don't use curved arrows that go backwards [[User:Tam10|Tam10]] ([[User talk:Tam10|talk]]) 14:22, 6 December 2016 (UTC))

=== Reaction beween SO<sub>2</sub> and the second s-cis diene fragment ===
<br />
IRC For ENDO Cycloaddition
[[File: DEC1_EXTENTION_ENDO_IRC_animation.gif|1200x400px]]
<br />
<br />
IRC For EXO Cycloaddition
[[File: DEC2_EXTENTION_EXO_IRC_animation.gif|1200x400px]]


<br />
After calculating the reaction barrier and reaction energies, it turns out that both of the ENDO and EXO pathways result in high reaction barriers and higher energies of products comparing with the reactions with the outside diene substrates. Therefore the reactions with the second diene are both thermodinamically and kinetically disfavored.

=== Conclusion ===
To sum up, different calculation methods of choice using Gaussian allows us to locate (by geometry optimization and frequency calculations), visualize (by extracting the IRC) and extract various information about particular reaction pathways (by checking the MOs for HOMO-LUMO interactions) on a PES. However, the program can only make predictions based on various theoretical approximations and models in quantum mechanics, the results of which are possible to disagree with experimental observations that are affected by a lot of other factors such as solvent, temperature, pH and so on. Also, the Gaussian calculations can be time consuming because complicated calculations can be established according to the calculation methods of choice.


=== Reference ===
[1] Lecture Handout 2; Sue Gibson; 2nd year Pericyclic Reactions; 
[2] [[http://www.kshitij-iitjee.com/Orbital-Hybridization-Bond-Angles]]
[3] [[http://www.kshitij-iitjee.com/Orbital-Hybridization-Bond-Angles]]