ChemWiki - User contributions [en]

FD915 Masters Proposal

2018-06-07T08:39:16Z

Fd915:

<jmol>
<jmolApplet>
<title>12 Epoxyobtusallene IV</title>
<color>white</color>
<size>800</size>
<script>frame 2; mo 6; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>12-EPOXYOBTUSALLENEIV_B3LYP_MINIMISATION.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

12-EPOXYOBTUSALLENEIV_B3LYP_MINIMISATION.LOG

FD915 Masters Proposal

2018-06-05T11:17:46Z

Fd915:

<jmol>
<jmolApplet>
<title>12 Epoxyobtusallene IV</title>
<color>white</color>
<size>500</size>
<script>frame 2; mo 6; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>12-EPOXYOBTUSALLENEIV_B3LYP_MINIMISATION.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

12-EPOXYOBTUSALLENEIV_B3LYP_MINIMISATION.LOG

FD915 Masters Proposal

2018-06-05T11:17:17Z

Fd915: Created page with " <jmol> <jmolApplet> <title>12 Epoxyobtusallene IV</title> <color>white</color> <size>300</size> <script>frame 2; mo 6; mo nodots nomesh fill translucent; mo titleformat "..."

<jmol>
<jmolApplet>
<title>12 Epoxyobtusallene IV</title>
<color>white</color>
<size>300</size>
<script>frame 2; mo 6; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>12-EPOXYOBTUSALLENEIV_B3LYP_MINIMISATION.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

12-EPOXYOBTUSALLENEIV_B3LYP_MINIMISATION.LOG

2018-03-14T10:26:42Z

Fd915:

= Transition Structure and Reactivity =

Felix de Courcy-Ireland
01062960
Transition States and Reactivity Computational Lab

==Introduction ==

Stationary points on potential energy surfaces (or ‘hypersurfaces’) are defined by a value of zero for the first derivative at that point. However, it is through the use of second derivatives that the distinction can be made between the various types of stationary points, including reactants, transition states, intermediates and products.

Products, reactants and intermediates are all minima on a potential energy surface. They are characterised by the fact that all paths through these points have a positive second derivative.

Transition states are characterised by having a single normal coordinate for which the second derivative is negative - this path is defined as a maximum – with all other paths through the transition state giving a positive second derivative. They are known as ‘saddle-points’. If a point on a potential energy surface is a maximum in more than one normal coordinate, it is diagnostic that there is a lower energy path for the reaction to proceed in the vicinity of this point.

The calculations carried out were effectively an interrogation of the potential energy surface of a simple reaction. The aims of the calculations were to locate the stationary points of the potential energy surfaces. This was done through using the Gaussian software, which traced paths across the potential energy surface, iteratively calculating the second derivative of each point on the path until a point with a zero first derivative was reached.

In this laboratory, a successfully located transition state was denoted by there being a single negative frequency in the frequency calculation, which when animated, showed the desired reaction coordinate for the reaction.

In the relatively simple reactions investigated in this laboratory, where the reactants pass through only one transition state, it is simple to calculate important quantities such as the activation energy and the reaction energy, from knowledge of the free energies of the reactants, products and transition states.

In this laboratory, two computational methods were used; the semi-empirical PM6 method and the DFT method B3LYP/6-31(d).

The semi-empirical PM6 is a less accurate method but not as computationally expensive. Therefore, it proved useful in optimising structures, be they reactants, products or transition states, in preparation for the use of the more accurate and expensive B3LYP method.
B3LYP/6-31(d) is based in density functional theory (DFT). DFT finds approximate solutions to unsolvable many-electron wave functions by analysis of functionals of the one electron density, ρ(r), of the molecule.

In some instances, use of chemical intuition in the location of transition states was most useful. There were cases where it proved fruitful to optimise a product or reactant(s) at the PM6 level, and then freeze the atoms immediately involved in the transition and further minimise the molecule around this frozen ‘transition state’ at the B3LYP level. A final calculation to find the transition state of at the B3LYP level was often more successful when using this ‘bond freezing’ technique.

==Exercise 1 ==

=== Reaction Scheme ===

[[File:Reaction_scheme.png|550px|center]]

=== Log Files ===

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - Butadiene + Ethene
! colspan = "2" style="background: #0D4F8B; color: white;" | Log Files - Butadiene + Ethene
|-
| Minimisation of butadiene (PM6)
| [[Media:CIS_BUTADIENE_ATTEMPT3.LOG| Minimisation of butadiene (PM6)]]
|-
| Minimisation of ethene (PM6)
| [[Media:ETHENE_MINIMISATION_ATTEMPT1.LOG|Minimisation of ethene (PM6)]]
|-
| Minimisation of cyclohexene (PM6)
| [[Media:CYCLOHEXENE_OPTIMISATION_ATTEMPT4.LOG|Minimisation of cyclohexene (PM6)]]
|-
| Transition State (PM6)
| [[Media:BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG|Transition State (PM6)]]
|-
| IRC (PM6)
| [[Media:IRC_FROM_CYCLOHEXENE_ATTEMPT1_BOTH_WAYS.LOG|IRC (PM6)]]
|}

=== MO diagram ===

[[File:Ex1_MO_diagram.png|500px|center| MO diagram for reaction of butadiene and ethene]]

The above molecular orbital diagram shows the interaction between the HOMOs and LUMOs of the ethene and butadiene. Each MO is assigned a symmetry value, from these assignments it can be seen that only orbitals of like symmetry interact. The orbital overlap integral is zero for the case of a symmetric-(anti-symmetric) interaction and non-zero for either a symmetric-symmetric or (anti-symmetric)-(anti-symmetric) interaction. The orbital energies are relative to one another rather than being absolute. An attempt has been made to draw the MOs of the transition state - the calculated MOs can be found in the table below.

=== Jmols of the HOMOs and LUMOS of Butadiene and Ethene ===

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Table 1 showing the HOMO and LUMO of ethene with relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Ethene
|-
! style="background: #0D4F8B; color: white;" | HOMO
| <jmol>
<jmolApplet>
<title>Ethene HOMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 6; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>ETHENE_MINIMISATION_ATTEMPT1.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.39228 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | LUMO
| <jmol>
<jmolApplet>
<title>Ethene LUMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 7; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>ETHENE_MINIMISATION_ATTEMPT1.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.04256 E<sub>h</sub>
|}

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Table 2 showing the HOMO and LUMO of butadiene, with relative energies in Hartrees.
! colspan = "3" style="background: #0D4F8B; color: white;" | Butadiene
|-
! style="background: #0D4F8B; color: white;" | HOMO
| <jmol>
<jmolApplet>
<title>Butadiene HOMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 11; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>CIS_BUTADIENE_ATTEMPT3.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.35899 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | LUMO
| <jmol>
<jmolApplet>
<title>Butadiene LUMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 12; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>CIS_BUTADIENE_ATTEMPT3.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01943 E<sub>h</sub>
|}

=== Jmols of the MOs in the Transition State ===

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Table 3 showing the four MOs produced from the interaction of the ethene and butadiene MOs with associated relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | MO 16
| <jmol>
<jmolApplet>
<title> MO 16 </title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 16; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.32754 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 17
| <jmol>
<jmolApplet>
<title> MO 17 </title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 17; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.32532 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 18
| <jmol>
<jmolApplet>
<title> MO 18 </title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 18; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.03</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01732 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 19
| <jmol>
<jmolApplet>
<title> MO 19 </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 19; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.03 </script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.03066 E<sub>h</sub>
|}

=== C-C Bond Lengths of Reactants, Transition State and Products ===

[[File:Reaction_scheme.png|550px|center]]

Reaction Scheme repeated for clarity when analysing the bond lengths below.

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 4 showing bond distances in the reactants.
! colspan = "2" style="background: #0D4F8B; color: white;" | Reactants
|-
! style="background: #0D4F8B; color: white;" | Bond
! style="background: #0D4F8B; color: white;" | C-C Bond Distance (Å)
|-
| C1-C2
| 1.33344
|-
| C2-C3
| 1.47077
|-
| C3-C4
| 1.33344
|-
| C5-C6
| 1.32731
|}

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 5 showing bond distances in the transition state
! colspan = "2" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | Bond
! style="background: #0D4F8B; color: white;" | C-C Bond Distance (Å)
|-
| C1-C2
| 1.37978
|-
| C2-C3
| 1.41104
|-
| C3-C4
| 1.37981
|-
| C4-C5
| 2.11452
|-
| C5-C6
| 1.38177
|-
| C6-C1
| 2.11469
|}

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 6 showing bond distances in the product
! colspan = "2" style="background: #0D4F8B; color: white;" | Cyclohexene
|-
! style="background: #0D4F8B; color: white;" | Bond
! style="background: #0D4F8B; color: white;" | C-C Bond Distance (Å)
|-
| C1-C2
| 1.49118
|-
| C2-C3
| 1.36309
|-
| C3-C4
| 1.49118
|-
| C4-C5
| 1.58345
|-
| C5-C6
| 1.56027
|-
| C6-C1
| 1.58345
|}

==== Discussion of carbon-carbon (C-C) bond lengths ====

The C-C single bond length of typical n-hydrocarbon is 1.533 Å, and the bond length of a typical C-C double bond is 1.33 Å. These lengths will be used as a point of comparison in the following discussion.

The two pairs of carbon atoms which make up the two alkene bonds in the butadiene (C1-C2 and C3-C4) lengthen from a separation of 1.33 Å to a separation of 1.49 Å. This can be rationalised by observing that these two double bonds in the butadiene reactant will be present in the cyclohexene product as single bonds in the cyclohexene ring.

The internal C-C bond in butadiene (C2-C3) shortens from a distance of 1.47 Å to a distance of 1.36 Å. Again, although this bond is shorter than a typical C-C single bond due to the conjugation of the butadiene, it shortens to a distance more typical of a C=C double bond which it adopts in the cyclohexene product.

The ethene C-C double bond lengthens from 1.327 Å to 1.560 Å. This is due to the carbons moving from the sp2 hybridised environment of the ethene molecule to the sp3 hybridised environment of the cyclohexene ring.

The two new C-C bonds which are formed upon the Diels-Alder [4+2] cycloaddition between the butadiene and ethene have a value of 2.115 Å in the transition state, which shortens to a value of 1.583 Å in the cyclohexene product. This change in bond distance shows a shortening of bond lengths from an effectively infinite value in the reactants, through to a bond distance of 1.583 Å in the cyclohexene product.

Van der Waals radius of a carbon atom is 1.7 Å, and this can be used in discussion of the partially formed bond between the ethene and butadiene in the transition state.

The length of the partly formed C-C bonds between butadiene and ethene is 2.115 Å. This is less than two times the van der Waals radius of carbon. Therefore, although there is not formally a bond formed between the atoms, the atoms are close enough to interact through van der Waals forces.

==== Vibration corresponding to the transition state ====

<jmol>
<jmolApplet>
<title> Reaction Path Vibration </title>
<color>white</color>
<size>500</size>
<script>frame 15; vibration 2;rotate x -20;</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

The Jmol above shows the vibration that corresponds to the reaction path at the transition state. It shows that the formation of the two bonds is synchronous - both bonds form simultaneously.

==Exercise 2 ==

=== Reaction Scheme ===

[[File:Ex2_reaction_scheme.png|500px|center]]

=== Log files ===

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - 1,3-Dioxole + Cyclohexadiene
! colspan = "2" style="background: #0D4F8B; color: white;" | Log Files - 1,3-Dioxole + Cyclohexadiene
|-
| Minimised Cyclohexadiene (B3LYP/6-31G(d))
| [[Media:CYCLOHEXADIENE_REPEAT_AGAIN_B3LYP_MINIMISATION.LOG|Minimised Cyclohexadiene (B3LYP/6-31G(d))]]
|-
| Minimised 1,3-Dioxole (B3LYP/6-31G(d))
| [[Media:13DIOXOLE_PROPER_MINIMISATION_B3LYP_ATTEMPT3.LOG|Minimised 1,3-Dioxole (B3LYP/6-31G(d))]]
|-
| Minimised exo-product (B3LYP/6-31G(d))
| [[Media:EXO_FULL_MOLECULE_B3LYP_MINIMISATION.LOG|Minimised exo-product (B3LYP/6-31G(d))]]
|-
| Minimised endo-product (B3LYP/6-31G(d))
| [[Media:ENDO_FULL_MOLECULE_INITIAL_MINIMISATION_B3LYP.LOG|Minimised endo-product (B3LYP/6-31G(d))]]
|-
| Minimised exo-transition state (B3LYP/6-31G(d))
| [[Media:EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG|Minimised exo-transition state (B3LYP/6-31G(d))]]
|-
| Minimised endo-transition state (B3LYP/6-31G(d))
| [[Media:ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG|Minimised endo-transition state (B3LYP/6-31G(d))]]
|-
| IRC on Exo Transition State (PM6)
| [[Media:EXO_FULL_MOLECULE_IRC_ATTEMPT1.LOG|IRC on Exo Transition State (PM6)]]
|-
| IRC on Endo Transition State (PM6)
| [[Media:ENDO_MOLECULE_PM6_IRC_ATTEMPT1.LOG|IRC on Endo Transition State (PM6)]]
|-
| Single Point Energy on Exo Transition State (B3LYP/6-31G(d))
| [[Media:EXO_FULL_MOLECULE_B3LYP_SINGLEPOINTENERGY.LOG|Single Point Energy on Exo Transition State (B3LYP/6-31G(d))]]
|-
| Single Point Energy on Endo Transition State (B3LYP/6-31G(d))
| [[Media:ENDO_MOLECULE_B3LYP_SINGLEPOINTENERGY.LOG|Single Point Energy on Endo Transition State (B3LYP/6-31G(d))]]
|}

=== Jmols and Molecular Orbital Diagram for the Endo Reaction of Cyclohexadiene and 1,3-Dioxole ===

[[File:Ex2_MO_diagram_endo.png|500px|right| MO diagram for reaction with endo arrangement between 1,3-dioxole and cyclohexadiene]]

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 7 showing the MOs of the endo transition state with associated relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | MO 40 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 40 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 40; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.19648 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 41 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 41 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 41; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.19052 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 42 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 42 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 42; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.00462 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 43 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 43 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 43; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01543 E<sub>h</sub>
|}

Above is the new MO diagram from the endo reaction of cyclohexadiene and 1,3-dioxole. The endo reaction involves the two oxygen atoms of the 1,3-dioxole tucking under butadiene fragment. This leads to a stabilising secondary orbital interaction which shall be discussed later. The MO's produced by the interaction of the cyclohexadiene and 1,3-dioxole orbitals are shown in Table 7.

=== Jmols and Molecular Orbital Diagram for the Exo Reaction of Cyclohexadiene and 1,3-Dioxole ===

[[File:Ex2_MO_diagram_exo.png|500px|right| MO diagram for reaction with exo arrangement between 1,3-dioxole and cyclohexadiene]]

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 8 showing the MOs of the exo transition state with associated relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | MO 40 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 40 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 40; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.19801 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 41 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 41 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 41; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.18560 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 42 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 42 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 42; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.00699 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 43 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 43 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 43; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01019 E<sub>h</sub>
|}

Above is the MO diagram from the exo transition state, with the MOs produced from interaction of the reactants HOMOs and LUMOs showing in Table 8. It should be noted that, similarly to the previous Diels-Alder reaction of butadiene and ethene, that only orbitals of like symmetry interact to give MOs in the transition state.

=== Inverse and Normal Electron Demand Diels-Alder Reactions ===

The reaction between 1,3-dioxole and cyclohexadiene is an inverse electron demand Diels-Alder reaction. Inverse electron demand Diels-Alder reactions are characterised by there being a smaller energy gap between the HOMO of the dienophile (1,3-dioxole) and the LUMO of the diene (cyclohexadiene) than the energy between the HOMO of the diene and the LUMO of the dienophile. (Clayden).

[[File:Inverse_normal_DA.png|400px|center|Figure highlighting the relative energy differences that lead to definition of normal and inverse electron demand Diels-Alder reactions.]]

=== Thermochemistry ===

If the formation of more than one product is possible from a given set of reactants, knowledge of reaction energies and activation energies can lead to assignment of the kinetic and thermodynamic products. The kinetically favourable product is that from the reaction which has the lowest energy barrier to conversion between reactants and products. The thermodynamically favourable product is the lower energy product relative to the reactant(s).

In the case that there are two possible products of a reaction, as is the situation in this Diels-Alder reaction between 1,3-dioxole and cyclohexadiene, it is not necessary for one product to be the thermodynamic product and the other product to be the kinetic product. It may be the case that one of the products is both the thermodynamic and kinetic product.

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Energies Endo Reaction
|-
! style="background: #0D4F8B; color: white;" | Molecule
! style="background: #0D4F8B; color: white;" | Gibbs Free Energy (kJ/mol)
|-
| 1,3dioxole
| -701,188.414
|-
| cyclohexadiene
| -612,592.877
|-
| 1,3dioxole and cyclohexadiene
| -1,313,781.291
|-
| Endo Transition State
| -1,313,621.479
|-
| Endo Product
| -1,313,848.695
|}

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Endo Reaction Parameters
|-
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
| 159.811
|-
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
| -67.404
|}

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Energies Exo Reaction
|-
! style="background: #0D4F8B; color: white;" | Molecule
! style="background: #0D4F8B; color: white;" | Gibbs Free Energy (kJ/mol)
|-
| 1,3dioxole
| -701,188.414
|-
| cyclohexadiene
| -612,592.877
|-
| 1,3dioxole and cyclohexadiene
| -1,313,781.291
|-
| Exo Transition State
| -1,313,613.639
|-
| Exo Product
| -1,313,845.098
|}

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Exo Reaction Parameters
|-
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
| 167.651
|-
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
| -63.807
|}

As shown in the above tables, the reaction barrier for the exo reaction is higher at a value of 167.651 kJ/mol, due to the transition state through which the reaction proceeds being of higher energy. Therefore, due to the lower reaction barrier of 159.811 kJ/mol, the endo product is kinetic product.

The endo product proceeds through a lower energy transition state due to the secondary orbital interactions between the oxygen p orbitals and the p orbitals of the internal carbons of the cis-butadiene fragment.

The Jmol of the HOMO of the endo transition state is repeated here, where it can be seen that there is a stabilising interaction between the orbitals which lie on the oxygens of 1,3-dioxole and the two carbons which eventually form the double bond in the product, for the portion of the orbital located in these areas is the same phase.

<jmol>
<jmolApplet>
<title> MO 41 - Endo Transition State </title>
<color>white</color>
<size>350</size>
<script> frame 2; mo 41; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.01 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

The endo product is again the thermodynamic product, as the energy of the endo product (-67.404 kJ/mol) is lower than that of the exo product (-63.807) - although the difference in energy is less in the products than in the reactants. This reason for the greater stabilisation of the endo product is possibly due to this same orbital interaction between the p orbitals on the oxygens with the alkene pi system. The fact that the difference in energy between the products is less than the difference in energy of the transition states is a consequence of this orbital interaction being poorer in the products.

==Exercise 3 ==

[[File:Ex3_reaction_scheme.png|500px|center]]

=== Log Files ===

==== Reaction of Sulfur Dioxide with Exo-Cyclic Cis-Butadiene Fragment ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - Dioxide + Exo-Cyclic Butadiene
! colspan = "2" style="background: #0D4F8B; color: white;" | Reaction of Sulfur Dioxide with Exo-Cyclic Cis-Butadiene Fragment
|-
| Xylylene Minimisation (PM6)
|[[File:XYLYLENE_MINIMISATION_PM6.LOG| Xylylene Minimisation (PM6)]]
|-
| Sulfur Dioxide Minimisation (PM6)
| [[File:SULFUR_DIOXIDE.LOG| Sulfur Dioxide Minimisation (PM6)]]
|-
| Cheletropic Transition State (PM6)]
| [[File:CHELETROPIC_PM6_PRODUCT_TRANSITIONSTATE_NOTFROZEN.LOG| Cheletropic Transition State (PM6)]]
|-
| Cheletropic Product (PM6)]
| [[File:CHELETROPIC_PRODUCT_INITIAL_MINIMISATION_2.LOG| Cheletropic Product (PM6)]]
|-
| Endo (Xylylene-SO2) Transition State (PM6)
| [[File:DA_ENDO_PM6_TRANSITIONSTATE_ATTEMPT1.LOG| Endo (Xylylene-SO2) Transition State (PM6)]]
|-
| Endo (Xylylene-SO2) Product (PM6)]
| [[File:DA_ENDO_PM6_PRODUCT_ATTEMPT1_NOTFROZEN_MINIMISATION.LOG| Endo (Xylylene-SO2) Product (PM6)]]
|-
| Exo (Xylylene-SO2) Transition State (PM6)]
| [[File:DA_EXO_PM6_TRANSITIONSTATE_SECOND.LOG| Exo (Xylylene-SO2) Transition State (PM6)]]
|-
| Cheletropic IRC (PM6)]
| [[File:CHELETROPIC_PM6_IRC_ATTEMPT1.LOG | Cheletropic IRC (PM6)]]
|-
| Endo (Xylylene-SO2) IRC (PM6)
| [[File:DA_ENDO_PM6_IRC_ATTEMPT1.LOG| Endo (Xylylene-SO2) IRC (PM6)]]
|-
| Exo (Xylylene-SO2) IRC (PM6)]
| [[File:DA_EXO_PM6_IRC_ATTEMPT1.LOG| Exo (Xylylene-SO2) IRC (PM6)]]
|}

==== Reaction of Sulfur with Alternate, Endo-Cyclic Cis-Butadiene Fragment of Xylylene ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - Dioxide + Endo-Cyclic Butadiene
! colspan = "2" style="background: #0D4F8B; color: white;" | Reaction of Sulfur Dioxide with Exo-Cyclic Cis-Butadiene Fragment
|-
| Endo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]
|[[File:ENDO EXTRA NEW TRANSITIONSTATE ATTEMPT2.LOG| Endo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Endo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:ENDO_EXTRA_NEW_MINIMISATION_ATTEMPT1.LOG| Endo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Exo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:EXO_EXTRA_TRANSITIONSTATE_ATTEMPT1.LOG| Exo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Exo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:EXO_EXTRA_MINIMISATION_ATTEMPT1.LOG| Exo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Cheletropic Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:CHELETROPIC_EXTRA_TRANSITIONSTATE_NOTFROZEN_ATTEMPT2.LOG| Cheletropic Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Cheletropic Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:CHELETROPIC_PRODUC_MINIMISATION_HOPEFULLY_ATTEMPT2.LOG| Cheletropic Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|}

=== GIFs of IRC for Reaction of Sulfur Dioxide and Xylylene ===

It can be seen in all of the GIFs of the reactions below that the 6-membered ring of the highly unstable xylylene molecule becomes aromatic on reaction with sulfur dioxide.

==== GIF for Endo Reaction ====

[[File:Endo_IRC_GIF.gif]]

==== GIF for Exo Reaction ====

[[File:Exo_IRC_GIF.gif]]

This reaction goes from products to reactants.

==== GIF for Cheletropic Reaction ====

[[File:Cheletropic_IRC_GIF.gif]]

=== Thermochemistry ===

{| class="wikitable"
! colspan = "3" style="background: #0D4F8B; color: white;" | Reaction Parameters For Reaction of Sulfur Dioxide With Xylylene
|-
! style="background: #0D4F8B; color: white;" | Reaction Type
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
|-
| style="background: #0D4F8B; color: white;" | Endo
| 81.724
| -99.050
|-
| style="background: #0D4F8B; color: white;" | Exo
| 85.709
| -99.714
|-
| style="background: #0D4F8B; color: white;" | Cheletropic
| 104.049
| -156.041
|}

The data above shows that the kinetic product is the endo product, due to it having the lowest reaction barrier of 111.946 kJ/mol. This is probably a consequence of the secondary orbitals interactions that were previously discussed, between the second sulfur oxygen and the carbon p orbitals that eventually form the alkene. The cheletropic reaction has the highest reaction barrier of 140.666 kJ/mol. This is potentially due to the 5-membered transition states that forms being of higher energy than that of the 6-membered transition states that form in the endo and exo forms of the reaction.

The thermodynamic product is also that from the endo reaction, due to it having the lowest energy (-156.041 kJ/mol). This is best explained in terms of bond enthalpies. In the endo and exo reactions, a sulfur-oxygen bond is broken and a carbon-oxygen and sulfur-carbon bond are formed, whereas in the cheletropic reaction, no bonds are broken, and two bonds are formed.

Although the cheletropic reaction proceeds through a higher energy transition states, the formation of two bonds with the loss of none leads to a lower energy product.

=== Reaction Profile for the Reaction of Sulfur Dioxide and the Exo-Cyclic Cis-Butadiene Fragment of Xylylene ===

[[File:Reaction_Profile_for_Exercise_3.png|500px|center]]

=== Further Comment ===

==== Discussion of Stability of Xylylene ====

Xylylene is an 8 electron system, and therefore anti-aromatic (Huckel's 4n rule, with n=2). These systems are often bent to avoid planarity of the p orbitals of the molecule, however xylylene is somewhat constrained to a planar configuration due to the sp<sup>2</sup) nature of the 6 ring carbons. This makes for a very strained molecule with a relative energy of 467.460 kJ/mol.

=== Reaction at the Endo-Cyclic, Cis-Butadiene Fragment ===

[[File:Ex3 reaction scheme extension.png|500px|center]]

Reaction scheme above shows the same reactions are possible at the endo-cyclic cis-butadiene fragment of xylylene.

{| class="wikitable"
! colspan = "3" style="background: #0D4F8B; color: white;" | Reaction Parameters For Reaction at the Alternate Cis-Butadiene Site of Xylylene
|-
! style="background: #0D4F8B; color: white;" | Reaction Type
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
|-
| style="background: #0D4F8B; color: white;" | Endo
| 111.946
| 16.197
|-
| style="background: #0D4F8B; color: white;" | Exo
| 119.783
| 20.668
|-
| style="background: #0D4F8B; color: white;" | Cheletropic
| 140.666
| 47.238
|}

The values for activation energy (reaction barrier) and reaction energy above show these reactions are far less favoured. All the reactions are endothermic with larger reactions barriers than the reactions at the original cis-butadiene fragment.

The reason for the positive values for the reaction energy (showing these are endothermic reactions) is due to the fact that there is no gain of aromaticity in these reactions, like there was previously.

The reason for the higher activation barrier is potentially a consequence of the endo-cyclic cis-butadiene fragment being harder to access, therefore the molecules need to adopt higher energy conformations in the transition states to yield the products

== Further Work ==

=== Ring Opening of the Dimethyl Ester Cyclobutene ===

[[File:Cyclobutene.png|500px|center]]

==== Log Files ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Ring Opening of the Dimethyl Ester Cyclobutene
! colspan = "2" style="background: #0D4F8B; color: white;" | Ring Opening of the Dimethyl Ester Cyclobutene
|-
| Dimethyl Ester of Cyclobutene
| [[Media:CYCLOBUTENE_REMINIMISATION_PM6_ATTEMPT1.LOG|Dimethyl Ester of Cyclobutene]]
|-
| Transition State
| [[Media:FRESHCYCLOBUTENE_PM6_TRANSITIONSTATE_ATTEMPT3.LOG|Transition State]]
|-
| Diene Product
| [[Media:DIENE_PRODUCT_MINIMISED_ATTEMPT1.LOG|Diene Product]]
|-
| IRC Calculation
| [[Media:BIS_CYCLOHEXENE_IRC_ATTEMPT1.LOG| IRC Calculations]]
|}

==== Jmols of MOs ====

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Jmols of MOs for Ring Opening of the Dimethyl Ester Cyclobutene
! style="background: #0D4F8B; color: white;" |
! style="background: #0D4F8B; color: white;" | HOMO
! style="background: #0D4F8B; color: white;" | LUMO
|-
| Reactant (ring closed)
|
|
|-
| Transition State
|
|
|-
| Product (ring open)
|
|
|}

==== Disrotatory or Conrotatory? ====

=== Cyclisation of [1,1']bicyclohexyl-1,1'-diene ===

[[File:Diene_one.png|500px|center]]

==== Log Files ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Cyclisation of [1,1']bicyclohexyl-1,1'-diene
! colspan = "2" style="background: #0D4F8B; color: white;" | Cyclisation of [1,1']bicyclohexyl-1,1'-diene
|-
| [1,1']bicyclohexyl-1,1'-diene
| [[Media:REACTANT_EXTENSION_MINIMISATION.LOG|[1,1']bicyclohexyl-1,1'-diene]]
|-
| Transition State
| [[Media:PRODUCT_TRANSITIONSTATE_ATTEMPT2.LOG|Transition State]]
|-
| Cyclised Product
| [[Media:PRODUCT_INITIAL_MINIMISATION.LOG|Cyclised Product]]
|-
| IRC Calculation
| [[Media:FRESHCYCLOBUTENE_PM6_IRC_ATTEMPT1.LOG|IRC Calculation]]
|}

==== Jmols of MOs ====

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Jmols of MOs for Cyclisation of [1,1']bicyclohexyl-1,1'-diene
! style="background: #0D4F8B; color: white;" |
! style="background: #0D4F8B; color: white;" | HOMO
! style="background: #0D4F8B; color: white;" | LUMO
|-
| Reactant (ring closed)
|
|
|-
| Transition State
|
|
|-
| Product (ring open)
|
|
|}

==== Disrotatory or Conrotatory? ====

==Conclusion==

A computational analysis of three different reactions of increasingly complexity were studied. Through the introduction to Gaussian, two computational methods were used (PM6 and B3LYP/6-31(d)) to optimise both the reactants and products of a reaction. From these initial minimisations, and using chemical intuition, the transition states for these reactions were located.

From using the relative energies of the transition states, reactants and products, thermodynamic properties of the systems were calculated (the activation energy (reaction barrier) and reaction energy). These values were used in discussion of the thermodynamic and kinetic products of the reaction.

Further work was carried out on two electrocyclic reactions, however the student ran out of time to conduct a full analysis of these reactions, for said was enjoying carrying out the calculations too much.

The tutorial and further reading provided facilitated an introduction to potential energy surfaces and the underlying quantum mechanics of the Gaussian code. A computational chemist must marry a working knowledge of the computational methods (and underlying approximations) with the savoir-faire of an experimental chemist.

== References ==

Rep:FD915 TRANSITION

2018-03-14T10:20:02Z

Fd915: /* Further Work */

= Transition Structure and Reactivity =

Felix de Courcy-Ireland
01062960
Transition States and Reactivity Computational Lab

==Introduction ==

Stationary points on potential energy surfaces (or ‘hypersurfaces’) are defined by a value of zero for the first derivative at that point. However, it is through the use of second derivatives that the distinction can be made between the various types of stationary points, including reactants, transition states, intermediates and products.

Products, reactants and intermediates are all minima on a potential energy surface. They are characterised by the fact that all paths through these points have a positive second derivative.

Transition states are characterised by having a single normal coordinate for which the second derivative is negative - this path is defined as a maximum – with all other paths through the transition state giving a positive second derivative. They are known as ‘saddle-points’. If a point on a potential energy surface is a maximum in more than one normal coordinate, it is diagnostic that there is a lower energy path for the reaction to proceed in the vicinity of this point.

The calculations carried out were effectively an interrogation of the potential energy surface of a simple reaction. The aims of the calculations were to locate the stationary points of the potential energy surfaces. This was done through using the Gaussian software, which traced paths across the potential energy surface, iteratively calculating the second derivative of each point on the path until a point with a zero first derivative was reached.

In this laboratory, a successfully located transition state was denoted by there being a single negative frequency in the frequency calculation, which when animated, showed the desired reaction coordinate for the reaction.

In the relatively simple reactions investigated in this laboratory, where the reactants pass through only one transition state, it is simple to calculate important quantities such as the activation energy and the reaction energy, from knowledge of the free energies of the reactants, products and transition states.

In this laboratory, two computational methods were used; the semi-empirical PM6 method and the DFT method B3LYP/6-31(d).

The semi-empirical PM6 is a less accurate method but not as computationally expensive. Therefore, it proved useful in optimising structures, be they reactants, products or transition states, in preparation for the use of the more accurate and expensive B3LYP method.
B3LYP/6-31(d) is based in density functional theory (DFT). DFT finds approximate solutions to unsolvable many-electron wave functions by analysis of functionals of the one electron density, ρ(r), of the molecule.

In some instances, use of chemical intuition in the location of transition states was most useful. There were cases where it proved fruitful to optimise a product or reactant(s) at the PM6 level, and then freeze the atoms immediately involved in the transition and further minimise the molecule around this frozen ‘transition state’ at the B3LYP level. A final calculation to find the transition state of at the B3LYP level was often more successful when using this ‘bond freezing’ technique.

==Exercise 1 ==

=== Reaction Scheme ===

[[File:Reaction_scheme.png|550px|center]]

=== Log Files ===

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - Butadiene + Ethene
! colspan = "2" style="background: #0D4F8B; color: white;" | Log Files - Butadiene + Ethene
|-
| Minimisation of butadiene (PM6)
| [[Media:CIS_BUTADIENE_ATTEMPT3.LOG| Minimisation of butadiene (PM6)]]
|-
| Minimisation of ethene (PM6)
| [[Media:ETHENE_MINIMISATION_ATTEMPT1.LOG|Minimisation of ethene (PM6)]]
|-
| Minimisation of cyclohexene (PM6)
| [[Media:CYCLOHEXENE_OPTIMISATION_ATTEMPT4.LOG|Minimisation of cyclohexene (PM6)]]
|-
| Transition State (PM6)
| [[Media:BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG|Transition State (PM6)]]
|-
| IRC (PM6)
| [[Media:IRC_FROM_CYCLOHEXENE_ATTEMPT1_BOTH_WAYS.LOG|IRC (PM6)]]
|}

=== MO diagram ===

[[File:Ex1_MO_diagram.png|500px|center| MO diagram for reaction of butadiene and ethene]]

The above molecular orbital diagram shows the interaction between the HOMOs and LUMOs of the ethene and butadiene. Each MO is assigned a symmetry value, from these assignments it can be seen that only orbitals of like symmetry interact. The orbital overlap integral is zero for the case of a symmetric-(anti-symmetric) interaction and non-zero for either a symmetric-symmetric or (anti-symmetric)-(anti-symmetric) interaction. The orbital energies are relative to one another rather than being absolute. An attempt has been made to draw the MOs of the transition state - the calculated MOs can be found in the table below.

=== Jmols of the HOMOs and LUMOS of Butadiene and Ethene ===

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Table 1 showing the HOMO and LUMO of ethene with relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Ethene
|-
! style="background: #0D4F8B; color: white;" | HOMO
| <jmol>
<jmolApplet>
<title>Ethene HOMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 6; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>ETHENE_MINIMISATION_ATTEMPT1.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.39228 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | LUMO
| <jmol>
<jmolApplet>
<title>Ethene LUMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 7; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>ETHENE_MINIMISATION_ATTEMPT1.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.04256 E<sub>h</sub>
|}

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Table 2 showing the HOMO and LUMO of butadiene, with relative energies in Hartrees.
! colspan = "3" style="background: #0D4F8B; color: white;" | Butadiene
|-
! style="background: #0D4F8B; color: white;" | HOMO
| <jmol>
<jmolApplet>
<title>Butadiene HOMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 11; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>CIS_BUTADIENE_ATTEMPT3.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.35899 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | LUMO
| <jmol>
<jmolApplet>
<title>Butadiene LUMO</title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 12; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>CIS_BUTADIENE_ATTEMPT3.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01943 E<sub>h</sub>
|}

=== Jmols of the MOs in the Transition State ===

{| class="wikitable" style="margin-left: auto; margin-right: auto; border: none;"
|+ align="bottom" style="caption-side: bottom" | Table 3 showing the four MOs produced from the interaction of the ethene and butadiene MOs with associated relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | MO 16
| <jmol>
<jmolApplet>
<title> MO 16 </title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 16; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.32754 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 17
| <jmol>
<jmolApplet>
<title> MO 17 </title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 17; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.32532 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 18
| <jmol>
<jmolApplet>
<title> MO 18 </title>
<color>white</color>
<size>250</size>
<script>frame 2; mo 18; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.03</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01732 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 19
| <jmol>
<jmolApplet>
<title> MO 19 </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 19; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.03 </script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.03066 E<sub>h</sub>
|}

=== C-C Bond Lengths of Reactants, Transition State and Products ===

[[File:Reaction_scheme.png|550px|center]]

Reaction Scheme repeated for clarity when analysing the bond lengths below.

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 4 showing bond distances in the reactants.
! colspan = "2" style="background: #0D4F8B; color: white;" | Reactants
|-
! style="background: #0D4F8B; color: white;" | Bond
! style="background: #0D4F8B; color: white;" | C-C Bond Distance (Å)
|-
| C1-C2
| 1.33344
|-
| C2-C3
| 1.47077
|-
| C3-C4
| 1.33344
|-
| C5-C6
| 1.32731
|}

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 5 showing bond distances in the transition state
! colspan = "2" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | Bond
! style="background: #0D4F8B; color: white;" | C-C Bond Distance (Å)
|-
| C1-C2
| 1.37978
|-
| C2-C3
| 1.41104
|-
| C3-C4
| 1.37981
|-
| C4-C5
| 2.11452
|-
| C5-C6
| 1.38177
|-
| C6-C1
| 2.11469
|}

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 6 showing bond distances in the product
! colspan = "2" style="background: #0D4F8B; color: white;" | Cyclohexene
|-
! style="background: #0D4F8B; color: white;" | Bond
! style="background: #0D4F8B; color: white;" | C-C Bond Distance (Å)
|-
| C1-C2
| 1.49118
|-
| C2-C3
| 1.36309
|-
| C3-C4
| 1.49118
|-
| C4-C5
| 1.58345
|-
| C5-C6
| 1.56027
|-
| C6-C1
| 1.58345
|}

==== Discussion of carbon-carbon (C-C) bond lengths ====

The C-C single bond length of typical n-hydrocarbon is 1.533 Å, and the bond length of a typical C-C double bond is 1.33 Å. These lengths will be used as a point of comparison in the following discussion.

The two pairs of carbon atoms which make up the two alkene bonds in the butadiene (C1-C2 and C3-C4) lengthen from a separation of 1.33 Å to a separation of 1.49 Å. This can be rationalised by observing that these two double bonds in the butadiene reactant will be present in the cyclohexene product as single bonds in the cyclohexene ring.

The internal C-C bond in butadiene (C2-C3) shortens from a distance of 1.47 Å to a distance of 1.36 Å. Again, although this bond is shorter than a typical C-C single bond due to the conjugation of the butadiene, it shortens to a distance more typical of a C=C double bond which it adopts in the cyclohexene product.

The ethene C-C double bond lengthens from 1.327 Å to 1.560 Å. This is due to the carbons moving from the sp2 hybridised environment of the ethene molecule to the sp3 hybridised environment of the cyclohexene ring.

The two new C-C bonds which are formed upon the Diels-Alder [4+2] cycloaddition between the butadiene and ethene have a value of 2.115 Å in the transition state, which shortens to a value of 1.583 Å in the cyclohexene product. This change in bond distance shows a shortening of bond lengths from an effectively infinite value in the reactants, through to a bond distance of 1.583 Å in the cyclohexene product.

Van der Waals radius of a carbon atom is 1.7 Å, and this can be used in discussion of the partially formed bond between the ethene and butadiene in the transition state.

The length of the partly formed C-C bonds between butadiene and ethene is 2.115 Å. This is less than two times the van der Waals radius of carbon. Therefore, although there is not formally a bond formed between the atoms, the atoms are close enough to interact through van der Waals forces.

==== Vibration corresponding to the transition state ====

<jmol>
<jmolApplet>
<title> Reaction Path Vibration </title>
<color>white</color>
<size>500</size>
<script>frame 15; vibration 2;rotate x -20;</script>
<uploadedFileContents>BUTADIENE_ETHENE_TRANSITIONSTATE_PM6.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

The Jmol above shows the vibration that corresponds to the reaction path at the transition state. It shows that the formation of the two bonds is synchronous - both bonds form simultaneously.

==Exercise 2 ==

=== Reaction Scheme ===

[[File:Ex2_reaction_scheme.png|500px|center]]

=== Log files ===

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - 1,3-Dioxole + Cyclohexadiene
! colspan = "2" style="background: #0D4F8B; color: white;" | Log Files - 1,3-Dioxole + Cyclohexadiene
|-
| Minimised Cyclohexadiene (B3LYP/6-31G(d))
| [[Media:CYCLOHEXADIENE_REPEAT_AGAIN_B3LYP_MINIMISATION.LOG|Minimised Cyclohexadiene (B3LYP/6-31G(d))]]
|-
| Minimised 1,3-Dioxole (B3LYP/6-31G(d))
| [[Media:13DIOXOLE_PROPER_MINIMISATION_B3LYP_ATTEMPT3.LOG|Minimised 1,3-Dioxole (B3LYP/6-31G(d))]]
|-
| Minimised exo-product (B3LYP/6-31G(d))
| [[Media:EXO_FULL_MOLECULE_B3LYP_MINIMISATION.LOG|Minimised exo-product (B3LYP/6-31G(d))]]
|-
| Minimised endo-product (B3LYP/6-31G(d))
| [[Media:ENDO_FULL_MOLECULE_INITIAL_MINIMISATION_B3LYP.LOG|Minimised endo-product (B3LYP/6-31G(d))]]
|-
| Minimised exo-transition state (B3LYP/6-31G(d))
| [[Media:EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG|Minimised exo-transition state (B3LYP/6-31G(d))]]
|-
| Minimised endo-transition state (B3LYP/6-31G(d))
| [[Media:ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG|Minimised endo-transition state (B3LYP/6-31G(d))]]
|-
| IRC on Exo Transition State (PM6)
| [[Media:EXO_FULL_MOLECULE_IRC_ATTEMPT1.LOG|IRC on Exo Transition State (PM6)]]
|-
| IRC on Endo Transition State (PM6)
| [[Media:ENDO_MOLECULE_PM6_IRC_ATTEMPT1.LOG|IRC on Endo Transition State (PM6)]]
|-
| Single Point Energy on Exo Transition State (B3LYP/6-31G(d))
| [[Media:EXO_FULL_MOLECULE_B3LYP_SINGLEPOINTENERGY.LOG|Single Point Energy on Exo Transition State (B3LYP/6-31G(d))]]
|-
| Single Point Energy on Endo Transition State (B3LYP/6-31G(d))
| [[Media:ENDO_MOLECULE_B3LYP_SINGLEPOINTENERGY.LOG|Single Point Energy on Endo Transition State (B3LYP/6-31G(d))]]
|}

=== Jmols and Molecular Orbital Diagram for the Endo Reaction of Cyclohexadiene and 1,3-Dioxole ===

[[File:Ex2_MO_diagram_endo.png|500px|right| MO diagram for reaction with endo arrangement between 1,3-dioxole and cyclohexadiene]]

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 7 showing the MOs of the endo transition state with associated relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | MO 40 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 40 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 40; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.19648 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 41 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 41 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 41; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.19052 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 42 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 42 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 42; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.00462 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 43 - Endo Transition State
| <jmol>
<jmolApplet>
<title> MO 43 - Endo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 43; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01543 E<sub>h</sub>
|}

Above is the new MO diagram from the endo reaction of cyclohexadiene and 1,3-dioxole. The endo reaction involves the two oxygen atoms of the 1,3-dioxole tucking under butadiene fragment. This leads to a stabilising secondary orbital interaction which shall be discussed later. The MO's produced by the interaction of the cyclohexadiene and 1,3-dioxole orbitals are shown in Table 7.

=== Jmols and Molecular Orbital Diagram for the Exo Reaction of Cyclohexadiene and 1,3-Dioxole ===

[[File:Ex2_MO_diagram_exo.png|500px|right| MO diagram for reaction with exo arrangement between 1,3-dioxole and cyclohexadiene]]

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Table 8 showing the MOs of the exo transition state with associated relative energies in Hartrees
! colspan = "3" style="background: #0D4F8B; color: white;" | Transition State
|-
! style="background: #0D4F8B; color: white;" | MO 40 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 40 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 40; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.19801 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 41 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 41 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 41; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.18560 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 42 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 42 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 42; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| -0.00699 E<sub>h</sub>
|-
! style="background: #0D4F8B; color: white;" | MO 43 - Exo Transition State
| <jmol>
<jmolApplet>
<title> MO 43 - Exo Transition State </title>
<color>white</color>
<size>250</size>
<script> frame 2; mo 43; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.02 </script>
<uploadedFileContents>EXO_FULL_MOLECULE_B3LYP_TRANSITIONSTATE.LOG</uploadedFileContents>
</jmolApplet>
</jmol>
| 0.01019 E<sub>h</sub>
|}

Above is the MO diagram from the exo transition state, with the MOs produced from interaction of the reactants HOMOs and LUMOs showing in Table 8. It should be noted that, similarly to the previous Diels-Alder reaction of butadiene and ethene, that only orbitals of like symmetry interact to give MOs in the transition state.

=== Inverse and Normal Electron Demand Diels-Alder Reactions ===

The reaction between 1,3-dioxole and cyclohexadiene is an inverse electron demand Diels-Alder reaction. Inverse electron demand Diels-Alder reactions are characterised by there being a smaller energy gap between the HOMO of the dienophile (1,3-dioxole) and the LUMO of the diene (cyclohexadiene) than the energy between the HOMO of the diene and the LUMO of the dienophile. (Clayden).

[[File:Inverse_normal_DA.png|400px|center|Figure highlighting the relative energy differences that lead to definition of normal and inverse electron demand Diels-Alder reactions.]]

=== Thermochemistry ===

If the formation of more than one product is possible from a given set of reactants, knowledge of reaction energies and activation energies can lead to assignment of the kinetic and thermodynamic products. The kinetically favourable product is that from the reaction which has the lowest energy barrier to conversion between reactants and products. The thermodynamically favourable product is the lower energy product relative to the reactant(s).

In the case that there are two possible products of a reaction, as is the situation in this Diels-Alder reaction between 1,3-dioxole and cyclohexadiene, it is not necessary for one product to be the thermodynamic product and the other product to be the kinetic product. It may be the case that one of the products is both the thermodynamic and kinetic product.

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Energies Endo Reaction
|-
! style="background: #0D4F8B; color: white;" | Molecule
! style="background: #0D4F8B; color: white;" | Gibbs Free Energy (kJ/mol)
|-
| 1,3dioxole
| -701,188.414
|-
| cyclohexadiene
| -612,592.877
|-
| 1,3dioxole and cyclohexadiene
| -1,313,781.291
|-
| Endo Transition State
| -1,313,621.479
|-
| Endo Product
| -1,313,848.695
|}

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Endo Reaction Parameters
|-
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
| 159.811
|-
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
| -67.404
|}

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Energies Exo Reaction
|-
! style="background: #0D4F8B; color: white;" | Molecule
! style="background: #0D4F8B; color: white;" | Gibbs Free Energy (kJ/mol)
|-
| 1,3dioxole
| -701,188.414
|-
| cyclohexadiene
| -612,592.877
|-
| 1,3dioxole and cyclohexadiene
| -1,313,781.291
|-
| Exo Transition State
| -1,313,613.639
|-
| Exo Product
| -1,313,845.098
|}

{| class="wikitable"
! colspan = "2" style="background: #0D4F8B; color: white;" | Exo Reaction Parameters
|-
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
| 167.651
|-
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
| -63.807
|}

As shown in the above tables, the reaction barrier for the exo reaction is higher at a value of 167.651 kJ/mol, due to the transition state through which the reaction proceeds being of higher energy. Therefore, due to the lower reaction barrier of 159.811 kJ/mol, the endo product is kinetic product.

The endo product proceeds through a lower energy transition state due to the secondary orbital interactions between the oxygen p orbitals and the p orbitals of the internal carbons of the cis-butadiene fragment.

The Jmol of the HOMO of the endo transition state is repeated here, where it can be seen that there is a stabilising interaction between the orbitals which lie on the oxygens of 1,3-dioxole and the two carbons which eventually form the double bond in the product, for the portion of the orbital located in these areas is the same phase.

<jmol>
<jmolApplet>
<title> MO 41 - Endo Transition State </title>
<color>white</color>
<size>350</size>
<script> frame 2; mo 41; mo nodots nomesh fill translucent; mo titleformat ""; set antialiasdisplay on; mo cutoff 0.01 </script>
<uploadedFileContents>ENDO_MOLECULE_B3LYP_FRESH_TRANSITIONSTATE2.LOG</uploadedFileContents>
</jmolApplet>
</jmol>

The endo product is again the thermodynamic product, as the energy of the endo product (-67.404 kJ/mol) is lower than that of the exo product (-63.807) - although the difference in energy is less in the products than in the reactants. This reason for the greater stabilisation of the endo product is possibly due to this same orbital interaction between the p orbitals on the oxygens with the alkene pi system. The fact that the difference in energy between the products is less than the difference in energy of the transition states is a consequence of this orbital interaction being poorer in the products.

==Exercise 3 ==

[[File:Ex3_reaction_scheme.png|500px|center]]

=== Log Files ===

==== Reaction of Sulfur Dioxide with Exo-Cyclic Cis-Butadiene Fragment ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - Dioxide + Exo-Cyclic Butadiene
! colspan = "2" style="background: #0D4F8B; color: white;" | Reaction of Sulfur Dioxide with Exo-Cyclic Cis-Butadiene Fragment
|-
| Xylylene Minimisation (PM6)
|[[File:XYLYLENE_MINIMISATION_PM6.LOG| Xylylene Minimisation (PM6)]]
|-
| Sulfur Dioxide Minimisation (PM6)
| [[File:SULFUR_DIOXIDE.LOG| Sulfur Dioxide Minimisation (PM6)]]
|-
| Cheletropic Transition State (PM6)]
| [[File:CHELETROPIC_PM6_PRODUCT_TRANSITIONSTATE_NOTFROZEN.LOG| Cheletropic Transition State (PM6)]]
|-
| Cheletropic Product (PM6)]
| [[File:CHELETROPIC_PRODUCT_INITIAL_MINIMISATION_2.LOG| Cheletropic Product (PM6)]]
|-
| Endo (Xylylene-SO2) Transition State (PM6)
| [[File:DA_ENDO_PM6_TRANSITIONSTATE_ATTEMPT1.LOG| Endo (Xylylene-SO2) Transition State (PM6)]]
|-
| Endo (Xylylene-SO2) Product (PM6)]
| [[File:DA_ENDO_PM6_PRODUCT_ATTEMPT1_NOTFROZEN_MINIMISATION.LOG| Endo (Xylylene-SO2) Product (PM6)]]
|-
| Exo (Xylylene-SO2) Transition State (PM6)]
| [[File:DA_EXO_PM6_TRANSITIONSTATE_SECOND.LOG| Exo (Xylylene-SO2) Transition State (PM6)]]
|-
| Cheletropic IRC (PM6)]
| [[File:CHELETROPIC_PM6_IRC_ATTEMPT1.LOG | Cheletropic IRC (PM6)]]
|-
| Endo (Xylylene-SO2) IRC (PM6)
| [[File:DA_ENDO_PM6_IRC_ATTEMPT1.LOG| Endo (Xylylene-SO2) IRC (PM6)]]
|-
| Exo (Xylylene-SO2) IRC (PM6)]
| [[File:DA_EXO_PM6_IRC_ATTEMPT1.LOG| Exo (Xylylene-SO2) IRC (PM6)]]
|}

==== Reaction of Sulfur with Alternate, Endo-Cyclic Cis-Butadiene Fragment of Xylylene ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Log Files - Dioxide + Endo-Cyclic Butadiene
! colspan = "2" style="background: #0D4F8B; color: white;" | Reaction of Sulfur Dioxide with Exo-Cyclic Cis-Butadiene Fragment
|-
| Endo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]
|[[File:ENDO EXTRA NEW TRANSITIONSTATE ATTEMPT2.LOG| Endo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Endo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:ENDO_EXTRA_NEW_MINIMISATION_ATTEMPT1.LOG| Endo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Exo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:EXO_EXTRA_TRANSITIONSTATE_ATTEMPT1.LOG| Exo Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Exo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:EXO_EXTRA_MINIMISATION_ATTEMPT1.LOG| Exo Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Cheletropic Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:CHELETROPIC_EXTRA_TRANSITIONSTATE_NOTFROZEN_ATTEMPT2.LOG| Cheletropic Transition State - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|-
| Cheletropic Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE
| [[File:CHELETROPIC_PRODUC_MINIMISATION_HOPEFULLY_ATTEMPT2.LOG| Cheletropic Product Minimisation - OTHER CIS-BUTADIENE FRAGMENT OF XYLYLENE]]
|}

=== GIFs of IRC for Reaction of Sulfur Dioxide and Xylylene ===

It can be seen in all of the GIFs of the reactions below that the 6-membered ring of the highly unstable xylylene molecule becomes aromatic on reaction with sulfur dioxide.

==== GIF for Endo Reaction ====

[[File:Endo_IRC_GIF.gif]]

==== GIF for Exo Reaction ====

[[File:Exo_IRC_GIF.gif]]

This reaction goes from products to reactants.

==== GIF for Cheletropic Reaction ====

[[File:Cheletropic_IRC_GIF.gif]]

=== Thermochemistry ===

{| class="wikitable"
! colspan = "3" style="background: #0D4F8B; color: white;" | Reaction Parameters For Reaction of Sulfur Dioxide With Xylylene
|-
! style="background: #0D4F8B; color: white;" | Reaction Type
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
|-
| style="background: #0D4F8B; color: white;" | Endo
| 81.724
| -99.050
|-
| style="background: #0D4F8B; color: white;" | Exo
| 85.709
| -99.714
|-
| style="background: #0D4F8B; color: white;" | Cheletropic
| 104.049
| -156.041
|}

The data above shows that the kinetic product is the endo product, due to it having the lowest reaction barrier of 111.946 kJ/mol. This is probably a consequence of the secondary orbitals interactions that were previously discussed, between the second sulfur oxygen and the carbon p orbitals that eventually form the alkene. The cheletropic reaction has the highest reaction barrier of 140.666 kJ/mol. This is potentially due to the 5-membered transition states that forms being of higher energy than that of the 6-membered transition states that form in the endo and exo forms of the reaction.

The thermodynamic product is also that from the endo reaction, due to it having the lowest energy (-156.041 kJ/mol). This is best explained in terms of bond enthalpies. In the endo and exo reactions, a sulfur-oxygen bond is broken and a carbon-oxygen and sulfur-carbon bond are formed, whereas in the cheletropic reaction, no bonds are broken, and two bonds are formed.

Although the cheletropic reaction proceeds through a higher energy transition states, the formation of two bonds with the loss of none leads to a lower energy product.

=== Reaction Profile for the Reaction of Sulfur Dioxide and the Exo-Cyclic Cis-Butadiene Fragment of Xylylene ===

[[File:Reaction_Profile_for_Exercise_3.png|500px|center]]

=== Further Comment ===

==== Discussion of Stability of Xylylene ====

Xylylene is an 8 electron system, and therefore anti-aromatic (Huckel's 4n rule, with n=2). These systems are often bent to avoid planarity of the p orbitals of the molecule, however xylylene is somewhat constrained to a planar configuration due to the sp<sup>2</sup) nature of the 6 ring carbons. This makes for a very strained molecule with a relative energy of 467.460 kJ/mol.

=== Reaction at the Endo-Cyclic, Cis-Butadiene Fragment ===

[[File:Ex3 reaction scheme extension.png|500px|center]]

Reaction scheme above shows the same reactions are possible at the endo-cyclic cis-butadiene fragment of xylylene.

{| class="wikitable"
! colspan = "3" style="background: #0D4F8B; color: white;" | Reaction Parameters For Reaction at the Alternate Cis-Butadiene Site of Xylylene
|-
! style="background: #0D4F8B; color: white;" | Reaction Type
! style="background: #0D4F8B; color: white;" | Reaction Barrier (kJ/mol)
! style="background: #0D4F8B; color: white;" | Reaction Energy (kJ/mol)
|-
| style="background: #0D4F8B; color: white;" | Endo
| 111.946
| 16.197
|-
| style="background: #0D4F8B; color: white;" | Exo
| 119.783
| 20.668
|-
| style="background: #0D4F8B; color: white;" | Cheletropic
| 140.666
| 47.238
|}

The values for activation energy (reaction barrier) and reaction energy above show these reactions are far less favoured. All the reactions are endothermic with larger reactions barriers than the reactions at the original cis-butadiene fragment.

The reason for the positive values for the reaction energy (showing these are endothermic reactions) is due to the fact that there is no gain of aromaticity in these reactions, like there was previously.

The reason for the higher activation barrier is potentially a consequence of the endo-cyclic cis-butadiene fragment being harder to access, therefore the molecules need to adopt higher energy conformations in the transition states to yield the products

== Further Work ==

=== Ring Opening of the Dimethyl Ester Cyclobutene ===

[[File:Cyclobutene.png|500px|center]]

==== Log Files ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Ring Opening of the Dimethyl Ester Cyclobutene
! colspan = "2" style="background: #0D4F8B; color: white;" | Ring Opening of the Dimethyl Ester Cyclobutene
|-
| Dimethyl Ester of Cyclobutene
| [[Media:CYCLOBUTENE_REMINIMISATION_PM6_ATTEMPT1.LOG|Dimethyl Ester of Cyclobutene]]
|-
| Transition State
| [[Media:FRESHCYCLOBUTENE_PM6_TRANSITIONSTATE_ATTEMPT3.LOG|Transition State]]
|-
| Diene Product
| [[Media:DIENE_PRODUCT_MINIMISED_ATTEMPT1.LOG|Diene Product]]
|-
| IRC Calculation
| [[Media:BIS_CYCLOHEXENE_IRC_ATTEMPT1.LOG| IRC Calculations]]
|}

=== Cyclisation of [1,1']bicyclohexyl-1,1'-diene ===

[[File:Diene_one.png|500px|center]]

==== Log Files ====

{| class="wikitable"
|+ align="bottom" style="caption-side: bottom" | Cyclisation of [1,1']bicyclohexyl-1,1'-diene
! colspan = "2" style="background: #0D4F8B; color: white;" | Cyclisation of [1,1']bicyclohexyl-1,1'-diene
|-
| [1,1']bicyclohexyl-1,1'-diene
| [[Media:REACTANT_EXTENSION_MINIMISATION.LOG|[1,1']bicyclohexyl-1,1'-diene]]
|-
| Transition State
| [[Media:PRODUCT_TRANSITIONSTATE_ATTEMPT2.LOG|Transition State]]
|-
| Cyclised Product
| [[Media:PRODUCT_INITIAL_MINIMISATION.LOG|Cyclised Product]]
|-
| IRC Calculation
| [[Media:FRESHCYCLOBUTENE_PM6_IRC_ATTEMPT1.LOG|IRC Calculation]]
|}

==Conclusion==

A computational analysis of three different reactions of increasingly complexity were studied. Through the introduction to Gaussian, two computational methods were used (PM6 and B3LYP/6-31(d)) to optimise both the reactants and products of a reaction. From these initial minimisations, and using chemical intuition, the transition states for these reactions were located.

From using the relative energies of the transition states, reactants and products, thermodynamic properties of the systems were calculated (the activation energy (reaction barrier) and reaction energy). These values were used in discussion of the thermodynamic and kinetic products of the reaction.

Further work was carried out on two electrocyclic reactions, however the student ran out of time to conduct a full analysis of these reactions, for said was enjoying carrying out the calculations too much.

The tutorial and further reading provided facilitated an introduction to potential energy surfaces and the underlying quantum mechanics of the Gaussian code. A computational chemist must marry a working knowledge of the computational methods (and underlying approximations) with the savoir-faire of an experimental chemist.

== References ==

File:FRESHCYCLOBUTENE PM6 IRC ATTEMPT1.LOG

2018-03-14T10:19:42Z

Fd915:

2018-03-14T08:20:35Z

Fd915:

Rep:FD915 TRANSITION

2018-03-14T08:18:35Z

Fd915: