View source for Cbl08:module 1

= Module 1: Organic =
==The Hydrogenation of Cyclopentadiene Dimer==
===Diels-Alder Cycloaddition:  Dimerisation of Cyclopentadiene===

[[Image:dimerisation_reaction_scheme.png|thumb|centre|500px|Fig.1.1a - Dimerisation of Cyclopentadiene]]



Cyclopentadiene can undergo self '''4π<sub>s</sub> +2π<sub>s</sub>''' Diels-Alder cycloaddition to form product with two possible conformations: the '''endo''' and the '''exo''' product. From the molecular calculation ( MM2 force field in ChemBio 3D Ultra) as shown in table 1.1i, the endo product is thermodynamically less stable than that of the exo product and the differece is largely contributed by the torsion strain in the endo dimer. This increase in torsion strain is most likely due to the steric interaction between the 6-membered ring with carbon bridge and the 5-membered ring.


{| align=center class="wikitable" border="1"
|+ '''Table 1.1i – ChemBio 3D Energy Minimisation data'''
! Energy/kcal mol<sup>-1</sup>!!Exo Product (1) !! Endo Product (2)!! Hydrogenated Product (3)!!Hydrogenated Product (4)
|-
|'''Bend''' ||align="center"|20.5593||align="center"|20.8436||align="center"|18.9520||align="center"|14.5538
|-
|'''Torsion'''||align="center"|7.6657||align="center"|9.5132||align="center"|12.1125||align="center"|12.5029
|-
|'''1,4 VDW'''||align="center"| 4.2424||align="center"|4.3237||align="center"|5.7159||align="center"|4.4976
|-
|'''Total Energy''' ||align="center"|31.8814 ||align="center"|33.9987||align="center"|35.9410||align="center"|31.1658
|}


[[Image:CpD orbitals new.png|thumb|right|400px|Fig.1.1b - Orbital interactions during dimerisation]]


Though the exo product is thermodynamically favoured, it is shown experimentally that the reaction yields exclusively the endo product.<ref name="ja01494a060">Stanley J. Cristol, Wolfgang K. Seifert, S. Barney Soloway, ''Bridged Polycyclic Compounds'','''1960''',''82(9)'',''2351-2356''.{{DOI|10.1021/ja01494a060}}</ref> By considering the frontier molecular orbitals of the reactants, the endo product is kinetically favoured due to secondary orbital stabilisation as shown in Fig.1.1b. In the endo transition state shown in fig. 1.1b, the HOMO of the dienophile will interact more favourably with the LUMO of the diene due to the secondary orbital stabilisation arise from better orbital overlap. This lowers its activation pathway and therefore the endo product is kinetically favoured.

























===Hydrogenation of ''Endo'' Cyclopentadiene Dimer ===

The unsaturated cyclopentadiene dimer can undergo further hydrogenation through H2, Pd/C, with possible attacks on two different sites: the C=C bond in the 5-membered ring  or the C=C bond in the norborene ring as shown in Fig. 1.1c.


[[Image:hydrogenation_cbl.png|thumb|centre|350px|Fig.1.1c - Hydrogenation of ''Endo'' Cyclopentadiene Dimer]]

From the thermodynamics calculations in Tabla 1.1, structure 3 raises the total energy of the molecule relative to the unhydrogenated endo dimer, whereas structure 4 is lower in energy and therefore it is thermodynamically more favourable.


[[Image:Fig1.1d.png|thumb|right|150px|Fig.1.1d - Completely reduced product]]

The main contribution to the difference in energy between structure 2 and structure 3/4 can be attributed to torsion and bending. Hydrogenation of C=C bond reliefs bending strain significantly but increases torsion strain as a result of carbon centre going from sp2 to sp3 with additional hydrogen atoms. In structure 3, hydrogenation of C=C bond in the 5 membered ring will force the added hydrogen atoms to be eclipsing each other due to ring conformational restriction and hence increase steric strains. Whereas for structure 4, the bending relief is much more significant in structure 4 than 3 because the ring strain is relieved with respect to the bridging carbon and the reduced C-C bond can rotate more freely.   With reference to the literature <ref name="Petroleum">D. Skala, J. Hanika, ''Petroleum and Coal'', '''2003''', Vol. 45, '''105-108'''.</ref>, initial hydrogenation of C=C bond in noborene ring is favoured over that in the 5 membered ring. The double bond in the 5 membered ring will be reduced upon further reduction as shown in Fig. 1.1d.

== Stereochemistry of Nucleophilic additions to a pyridinium ring (NAD<sup>+</sup> analogue)==

===Nucleophilic Addition of Prolinol with Grignard Reagent, MeMgI===


[[Image:updated_prolinol.png|thumb|center|700px|Fig.1.2a - Nucleophilic Addition of Prolinol with Grignard Reagent, MgI
]]

Methylation of pyridine rign can be highly stereo- and regioselective through nucleophilic addition with Grignard reagents. The methyl group is expected to add specifically at the C(4) position of the pyridine ring due to chelate effect from the magnesium and the carbonyl group. Selectivities as high as 99:1 is observed experimentally with bulkier halide substituent such as iodide in the Grignard reagent.<ref name="jo00356a016"> A. G. Shultz, L. Flood and J. P. Springer, J. Org. Chemistry, '''1986''', 51, ''838'' {{DOI|10.1021/jo00356a016}}</ref> From the product shown in structure 6, there are two chiral centres: one at the alkylating centre and one at the hydrogen anti to the alkylated carbon. 


In order to explain the stereoselectivity, the position of the carbonyl group relative to the alkylating carbon centre is crucial and hence the dihedral angle between the C=O bond and C-H bond, and the conformation minima of the reactant has first to be determined.


To find potential different dihedral angles between the carbonyl bond and the aromatic ring under minimised energy of Prolinol structure 5, the MM2 energy minimisation was first carried out on automatic geometry after inputted the structure. A first set of data was achieved with dihedral angle of 9.0°. The molecule was then manually distorted with both oxygen in the molecule pulled 90° below the π-system and the result had no significant differences between the optimised energy and the dihedral angle. Hence the same optimisation was done with MMFF94 force field calculation and the results were tabulated as follow:


{|align=center class="wikitable" border="1"
|+ '''Table 1.2i – ChemBio 3D Energy Minimisation data using MM2 and MMFF94 Force Field of Structure 5'''
!   Calculation Type    !! Pre-Opt Structure!! Dihedral Angle/<sup>o</sup>!!Energy/kcal mol<sup>-1</sup>
|-
|'''MM2'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MM2_9degrees.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol>
 ||Automatic input structure||align="center"|9.0||align="center"|43.1301
|-
|'''MM2'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MM2_11.5degrees.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol>||Both Oxygen 90<sup>o</sup> below the pi-system||align="center"|11.5||align="center"|43.1436
|-
|'''MMFF94'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MMFF94_2degrees.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol>|| Automatic input structure||align="center"|-1.0||align="center"|57.4874
|-
|'''MMFF94'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MMFF94_13degrees.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol> ||Both Oxygen 90<sup>o</sup> below the pi-system ||align="center"|13.2||align="center"|57.4256
|-
|'''MMFF94'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MMFF94_19.7degrees.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol>|| Carbonyl O 90<sup>o</sup> below the pi ring and enol O 90<sup>o</sup> above the ring||align="center"|19.7||align="center"|58.7154
|-
|'''MMFF94'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MMFF94_36degrees.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol> ||System widely stretched and strained ||align="center"|36.4||align="center"|59.1591
|}

*The magnesium atom is excluded from modelling as the molecular mechanics/energy minisation was unable to calculate with metal atoms.


From the calculations of MMFF94, two conformation minima were found with dihedral angle of -1.0° and 13.2°. This is in good agreement with the literature<ref name="jo00356a016"> A. G. Shultz, L. Flood and J. P. Springer, J. Org. Chemistry, '''1986''', 51, ''838'' {{DOI:10.1021/jo00356a016}}</ref> in which one conformer has the carbonyl group and the pyridinium ring coplanar, and the other conformer has the carbonyl group above the plane of the pyridinium ring and anti to the hydrogen at the other chiral centre. Taking either of these two conformers as the starting molecule, the incoming magnesium in the Grignard reagent will chelate onto the carbonyl oxygen from above the plane. This chelate effect restricts the position of the methyl group, allowing it only to attack from the top of the plane which is anti to the chiral hydrogen.

===Nucleophilic Addition to Pyridinium Ring with Aniline===


[[Image:pyridinium_7_to_8.png|thumb|center|700px|Fig.1.2b - Nucleophilic Addition of PAD+ analogue with aniline
]]


Methylation onto molecule 7 displays a total opposite stereochemistry with the aniline attacks anti to the carbonyl. In this reaction scheme, there is no Grignard reagent present to chelate onto the carbonyl group and hence steric factor dominates in this reaction.

Regardless of the distortion extent of the starting configurations, the molecule optimized into similar structures, with the carbonyl group syn to the chiral centre as marked in *.Without the presence of chelate effect, the incoming aniline favours to attack on the face anti to the carbonyl group for both steric and electronic considerations. The aniline has a bulky phenyl group attached and hence it attack from anti face minimised steric clashes. Also, the presence of phenyl group makes the aniline a highly electron rich nucleophile and hence anti attack is favoured to reduce electronic repulsion.

{| align=center class="wikitable" border="1"
|+ '''Table 1.2ii – ChemBio 3D Energy Minimisation data using MMFF94 Force Field of Structure 7'''
! Calculation Type !! Dihedral Angle/<sup>o</sup>!!Energy/kcal mol<sup>-1</sup>
|-
| '''MMFF94'''<jmol>
<jmolAppletButton>
<uploadedFileContents>MMFF94_35d.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol>||align="center"|-35.0||align="center"|99.6127
|-
| '''MMFF94'''<jmol>
<jmolAppletButton>
<uploadedFileContents>Structure7_37d.mol</uploadedFileContents>
<text>3D Jmol</text>
</jmolAppletButton>
</jmol>||align="center"|-37.2||align="center"|99.6799
|}


The simulation ran by ChemBio 3D program is a thermodynamic approach that explains the relative energies of the molecules. However, this is insufficient to have a full understanding of a reaction as reaction depends on both kinetics and thermodynamics. To improve the model, the transition state of the molecule can be studied and the relative energy and orbital alignment of respective HOMOs and LUMOs can also be taken into consideration to determine the preferred outcome.




== Stereochemistry and Reactivity of an Intermediate in the Synthesis of Taxol==


===Introduction===
'''Atropisomers''' are stereoisomers that have high kinetic barrier for rotating its single bond during isomer interconversion. This type of isomerism is commonly found in molecules with large rings in which its steric strains impose by the ring also play a part in rising its kinetic barrier for interconversion.

===Results and Analysis===
In this example, two atropisomers are found in the intermediate of taxol synthesis with restricted bond rotation in C-C=O bond as shown in Fig. 1.3a. To determine the most thermodynamically stable isomers, both structures were optimised with MM2 and MMFF94 force field.

Upon minimising structure 9 and 10, the optimised structures are both isomers with cyclohexane ring in twistboat conformations. This example is hence further investigated with the cyclohexane ring manually manipulated into chair conformation and the relative energy data were tablulated as shown in Table 1.3i:

[[Image:Taxol_atrop.png|thumb|right|443px|Fig.1.3a - Taxol Intermediate Atropisomers]]
{| align=center class="wikitable" border="1"
|+ '''Table 1.3i – Energy Minimisation data of Atropisomers using MM2 and MMFF94 Force Field'''
!   Structure  !! MM2 Energy/kcal mol<sup>-1</sup>!!MMFF94 Energy/kcal mol<sup>-1</sup>
|-
|'''9-Twistboat'''<jmol>
<jmolAppletButton>
<uploadedFileContents>9_twistboat.mol</uploadedFileContents>
<text>9-Twistboat</text>
</jmolAppletButton>
</jmol>||align="center"|54.3865||align="center"|79.3064
|-
|'''9-Chair'''<jmol>
<jmolAppletButton>
<uploadedFileContents>9_chair.mol</uploadedFileContents>
<text>9-Chair</text>
</jmolAppletButton>
</jmol>||align="center"|48.8929||align="center"|70.5486
|-
|'''10-Twistboat'''<jmol>
<jmolAppletButton>
<uploadedFileContents>10_twistboat.mol</uploadedFileContents>
<text>10-Twistboat</text>
</jmolAppletButton>
</jmol>||align="center"|58.7616||align="center"|81.5824
|-
|'''10-Chair'''<jmol>
<jmolAppletButton>
<uploadedFileContents>10_chair.mol</uploadedFileContents>
<text>10-Chair</text>
</jmolAppletButton>
</jmol>||align="center"|53.1676||align="center"|74.8324
|}

From the 2D structure, structure 10 seems to be thermodynamically more stable than 9 due to lesser steric clashes between the bridging carbon and the carbonyl oxygen. However when the structure is modelled in 3D, the other end of the norborane ring tends to fold down and inwards, making the two methyl groups on the bridge pointing outwards and hence relieves the steric strains. This observation is supported by the modelling results where structure 9 in chair conformation is the most thermodynamically stable isomer whereas isomer 10 with twistboat conformation is the least stable isomer with highest energy. Isomers in structure 9 experience instability that arises from the steric clash between the oxygen and the bridging carbons. Similar steric clashes are also present in structure 10 isomers where the downward facing carbonyl oxygen clashes with the carbons in norborane ring. On top of that, structure 10 isomers also experience steric clash with the β–carbon in the cyclohexane ring which is closer and hence stronger steric interactions. This interaction is more significant in twistboat conformation than the chair conformation. Therefore if given enough energy to rotate the bond, structure 10 will isomerise to the more stable structure 9 as it is thermodynamically favourable.

For further addition to the carbonyl bond, the incoming group will attrack from Burgi- Dunitz angle of 109° to the carbon.  From the 3D view of the molecules, addition will occur more favourably for structure 9 than 10 as the Burgi-Dunitz angle for 9 is spatially more available whereas that for 10 is more hindered by the hydrogens on the cyclohexane ring and bridging carbons.

===Further Discussion: Reaction at Bridgehead C=C double Bond===
[[Image:Hydrogenated_taxol_9.png|thumb|right|150px|Fig.1.3b - Hydrogenated twistboat of Structure 9<br clear="all" />
<jmol>
<jmolAppletButton>
<uploadedFileContents>Hydrogenated_9_twistboat.mol</uploadedFileContents>
<text>Hydrogenated 9</text>
</jmolAppletButton>
</jmol>]]

Both atropisomers studied includes a bridgehead double bond and a medium size polycyclic ring system. This double bond display unusually high stability and is regarded as a hyperstable alkene.<ref name="ja000398a003">W.F. Maier, P.V.R.Schleyer, J.Am.Chem. Rev., 1976,7, 31{{DOI|10.1021/ja00398a003}}</ref> A hyperstable alkene is thermodynamically more stable than its hydrogenated derivative and hence reluctant to undergo hydrogenation. The reluctance of hydrogenation can be attributed to as suggested in two literatures <ref name="ja000398a003">W.F. Maier, P.V.R.Schleyer, J.Am.Chem. Rev., 1976,7, 31{{DOI|10.1021/ja00398a003}}</ref> <ref name="ja000274a016">Alan B. McEwen, Paul v. R. Schleyer,J. Am. Chem. Soc., '''1986''', 108 ('''14'''), pp 3951–3960{{DOI: 10.1021/ja00274a016}}</ref> that the cage structure adds extra stability to the alkene bond and the carbon centre will carbon centre will experience more strains when it goes from sp2 to sp3 in this conformationally restricted ring system.

A simple way to justify is to compare the optimised energy of the hydrogenated specie and its parent compound. The twistboat conformation of structure 9 is hydrogenated and it has a higher energy of 65.9243kcal mol<sup>-1</sup>.




==Modelling Using Semi-empirical Molecular Orbital Theory==
===Part 1: Regioselective Addition of Dichlorocarbene===
Predicting the major product in a reaction with classic mechanic parameters such as thermodynamics and sterics can sometimes proof to be limited, for example in a kinetically controlled reaction. As such, understanding energies of the molecular orbital that is built on quantum mechanics will be useful to learn about the electron donating /accepting power of the specie and the kinetically favoured pathway.

[[Image:Updated_carbene.png|thumb|center|500px|Fig.1.4a - Electrophilic addition of Structure 12 and dichlorocarbene]]

This example investigates the electrophilic addition of dichlorocarbene with structure 12 as shown above. In structure 12, two nucleophilic sites of C=C bond are susceptible to electrophilic addition from potential electrophile such as dichlorocarbene. This electrophilic addition reaction involves a flow of electron density from the nucleophile to the electrophile during bond formation and hence it is interesting to investigate the HOMO-LUMO interaction in this process.

Using the MOPAC Interface followed by PM6 calculation method, the wavefunction of structure 12 is calculated and the relevant molecular orbitals were modelled as shown in the figures below.

{| align=center class="wikitable" border="1"
|+ '''Table 1.4i – Calculated Molecular Orbitals of Structure 12'''
|-
| ||align="center"|[[Image:cbl_homo-1.png|thumb|150px|Fig.1.4b - HOMO-1 of Structure 12]]||align="center"|[[Image:cbl_homo.png|thumb|150px|Fig.1.4c - HOMO of Structure 12]]||align="center"|[[Image:cbl_lumo.png|thumb|centre|150px|Fig.1.4d - LUMO of Structure 12]]||align="center"|[[Image:cbl_lumo+1.png|thumb|150px|Fig.1.4e - LUMO+1 of Structure 12]]||align="center"|
[[Image:cbl_lumo+2.png|thumb|right|150px|Fig.1.4f - LUMO+2 of Structure 12]]
|-
| Comments||align="center"|Electron density localised on the exo alkene bond because this π bond received favourable stabilising interaction from C-Cl σ* orbital and hence lower in energy than the endo π orbitals||align="center"|Electron density concentrated on endo alkene π orbitals, increases the nucleophilicity in this bond and hence susceptible to electrophilic addition||align="center"|resembles the shape of π* orbital of exo alkene||align="center"|C-Cl σ* orbital that is aligned with π* orbital of exo alkene||align="center"|Electron density localised on the π* orbital of the endo alkene
|}

[[Image:Lumo_carbenee.png|thumb|right|150px|Fig.1.4g - LUMO of dichlorocarbene]]

With the HOMO of the nucleophile has its largest electron density localised in the endo alkene bond and the LUMO of the dichlorocarbene eletrophile in the largely p-based antibonding orbital of C-Cl bond, the HOMO of the nucleophile will donate its electron density to the LUMO of the dichlorocarbene when the orbitals are correctly aligned (i.e. the LUMO came in perpendicular to the HOMO). This addition will be highly favoured as there is no steric concern in here.


With the HOMO of the nucleophile has its largest electron density localised in the endo alkene bond and the LUMO of the dichlorocarbene eletrophile in the largely p-based antibonding orbital of C-Cl bond, the HOMO of the nucleophile will donate its electron density to the LUMO of the dichlorocarbene when the orbitals are correctly aligned (i.e. the LUMO came in perpendicular to the HOMO). This addition will be highly favoured as there is no steric concern in here.




===Part 2: Influence of C-Cl bond on the Vibrational Frequency===
[[Image:Cbl_mono_alkene.png|thumb|right|80px|Fig.1.4h - Structure of 12i]]

In this section, to study the influence of the C-Cl bond on the vibrational energies of the original structure 12 and 12i as shown in Fig. 1.4h , IR spectra for both of the molecules are modelled.

 


[[Image:Cbl_carbene_IR.png|thumb|350px|Fig.1.4i - IR of Structure 12]]
[[Image:Cbl_hydrogenated_IR.png|thumb|350px|Fig.1.4j - IR of Structure 12i]]


The key bond vibrations are highlighted as follow:

 

{| align=center class="wikitable" border="1"
|+ '''Table 1.2ii – ChemBio 3D Energy Minimisation data using MMFF94 Force Field of Structure 7'''
! Molecule !! C-Cl stretch/ cm<sup>-1</sup>!!''syn''C=C stretch/cm<sup>-1</sup>!!''anti''C=C stretch/cm<sup>-1</sup>
|-
| Structure 12||align="center"|770.65||align="center"|1757.48||align="center"|1737.59
|-
| Structure 12i||align="center"|775.03||align="center"|1753.66||align="center"|-
|} 




The C-Cl bond increases in its vibrational energy as it goes from 12 to 12i and this indicates that the C-Cl bond is stronger in 12i than 12. In structure 12, the C-Cl bond is aligned in anti-periplanar (app) fashion to the πC=C,anti orbital, allowing electron density from π double bond to donate into σ*C-Cl orbital and hence weakening the C-Cl bond due to this app stabilisation. However in structure 12i, this app alignment is lost and electron density is now more concentrated on the C-Cl bond. This increase in localised electron density in turn donates into the π orbital of syn C=C bond and stabilised it, hence leading to a weaker bond and a decrease in syn C=C stretching frequency.











.

== Mini Project - Investigation on ''Endo'' and ''Exo'' Camphor using DFT based Molecular Orbital Method==
===Introduction===

'''''Camphor''''' is a waxy, white and aromatic natural product that can be found in the wood of  Camphor Laurel, a large evergreen tree. Synthetically, camphor can also be produced from the oil of turpentine. Its aromatic scent is commonly used as an ingredient in cooking, as an embalming fluid for medicinal purposes and in modern days, it is used as a plsticizer for nitrocellulose as a moth repellent and even in fireworks. Further reactions such as oxidation, bromination and Mannich reaction can be carried out on camphor. In this mini project, the Mannich reaction of 3-(hydroxymethylene)-d-camphor is studied with the following reaction scheme:

[[Image:Rxn_scheme_camphor200.png|thumb|center|600px|Fig.2a - Mannich Reaction of 3-(hydroxymethylene)-d-camphor]]



Hine and co-workers first treated 3-(hydroxymethylene)-d-camphor with dimethylamine followed byhydrogenation in the presence of platinum oxide and hydrochloric acid<ref name="hine">Hine. J.,W.-S.,J.Am.Chem.Soc.9 75, 97,3550-3551</ref>. This results in a very good yield of 90% with 80% exo-isomer and 20% endo isomer. The presence of this diastereomers was proven with NMR results and could be isolated with repeated crystallisation.

To identify and differentiate between isomers, <sup>13</sup>C NMR and <sup>1</sup>H NMR will be suitable to indentify as the protons and the carbons couples to each other differently. Also, optical rotation will be good to identify stereoisomers as it will rotate the plane polarized light in opposite directions. It is suggested by Natalie and co-workers that the difference in the dihedral angles between H<sub>a</sub> and H<sub>b</sub> does influence the H coupling in <sup>1</sup>NMR<ref name="jo00246a042">Natalie L. McClure, Gui Yuan Dai, Harry S.Mosher,J.Org.Chem., 1988, 53(11), pp2617-2620. {{DOI|10.1021/jo00246a042}}</ref>. In this example, camphor derivatives were modelled and optimised to predict <sup>13</sup>C NMR spectra.

===Predict <sup>13</sup>C NMR spectra of Camphor derivative isomers===

The optimised energies of both diastereomers were first run and tabulated in '''Table 2i'''. It is in good agreement with theoretical expectation that the Endo product is thermodynamically stable based on steric considerations similar to that as discussed in 1.1.

{| align=center class="wikitable" border="1"
|+ '''Table 2i – ChemBio 3D Energy Minimisation data of camphor product derivatives'''
! Energy/kcal mol<sup>-1</sup>!!Exo Product <jmol>
<jmolAppletButton>
<uploadedFileContents>camphor_exo.mol</uploadedFileContents>
<text>Exo</text>
</jmolAppletButton>
</jmol> !! Endo Product <jmol>
<jmolAppletButton>
<uploadedFileContents>camphor_endo.mol</uploadedFileContents>
<text>Endo</text>
</jmolAppletButton>
</jmol>
|-
|'''Stretch'''||align="center"|2.0363||align="center"|1.9525
|-
|'''Bend''' ||align="center"|14.8838||align="center"|14.1463
|-
|'''Stretch-Bend'''||align="center"|-0.2020||align="center"|-0.1154
|-
|'''Torsion'''||align="center"|15.5082||align="center"|16.0331
|-
|'''Non-1,4 VDW'''||align="center"| -0.3975||align="center"|-1.9941
|-
|'''1,4 VDW'''||align="center"|8.8985||align="center"|9.4383
|-
|'''Dipole-Dipole'''||align="center"|0.6232||align="center"|-0.1182
|-
|'''Total Energy''' ||align="center"|41.9506 ||align="center"|39.3424
|}

The computed camphor derivative is fully optimised at the DFT-mpw1pw91/6-31(d,p) level and a comparison between the literature and the predicted values are tabulated as follow.

[[Image:nmr_assigned.png|thumb|center|350px|Fig.2b -Labelled carbons to corresponding NMR values]]



{| align=center class="wikitable" border="1"
|+ '''Table 2ii – NMR data of ''endo'' Camphor derivative'''
! Corresponding Carbon!!Predicted Value/ppm !! Lit. Value/ppm !!Difference/ppm
|-
|align="center"|'''1''' ||align="center"|37.3194||align="center"|31.3||align="center"|-6.0194
|-
|align="center"|'''2'''||align="center"|13.7768||align="center"|49.9||align="center"|36.1232
|-
|align="center"|'''3'''||align="center"|16.9504||align="center"|20.6||align="center"|3.6496
|-
|align="center"|'''4''' ||align="center"|50.6142||align="center"|56.3||align="center"|5.6858
|-
|align="center"|'''5''' ||align="center"|205.378||align="center"|222.6||align="center"|17.222
|-
|align="center"|'''6'''||align="center"|41.1653||align="center"|46.6||align="center"|5.4347
|-
|align="center"|'''7'''||align="center"|37.0318||align="center"|46.4||align="center"|9.3682
|-
|align="center"|'''8''' ||align="center"|13.1555||align="center"|19.4||align="center"|6.2445
|-
|align="center"|'''9'''||align="center"|12.8694||align="center"|18.9||align="center"|6.0306
|-
|align="center"|'''10''' ||align="center"|11.9053||align="center"|9.3||align="center"|-2.6053
|-
|align="center"|'''11''' ||align="center"|46.6986||align="center"|59.8||align="center"|7.0143
|-
|align="center"|'''12'''||align="center"|37.1857||align="center"|44.2||align="center"|7.0143
|-
|align="center"|'''13'''||align="center"|31.4493||align="center"|43.8||align="center"|12.3507
|}





{| align=center class="wikitable" border="1"
|+ '''Table 2ii – NMR data of ''Exo'' Camphor derivative'''
! Corresponding Carbon!!Predicted Value/ppm !! Lit. Value/ppm !!Difference/ppm
|-
|align="center"|'''1''' ||align="center"|21.6696||align="center"|29.2||align="center"|7.5304
|-
|align="center"|'''2'''||align="center"|40.2072||align="center"|50.2||align="center"|9.9929
|-
|align="center"|'''3'''||align="center"|16.2684||align="center"|28.4||align="center"|12.1316
|-
|align="center"|'''4''' ||align="center"|50.9022||align="center"|58.4||align="center"|7.4978
|-
|align="center"|'''5''' ||align="center"|204.854||align="center"|222.9||align="center"|18.046
|-
|align="center"|'''6'''||align="center"|45.0371||align="center"|48.6||align="center"|3.5629
|-
|align="center"|'''7'''||align="center"|37.9855||align="center"|47.7||align="center"|9.7145
|-
|align="center"|'''8''' ||align="center"|15.3405||align="center"|21.1||align="center"|5.7595
|-
|align="center"|'''9'''||align="center"|13.5321||align="center"|19.8||align="center"|6.2679
|-
|align="center"|'''10''' ||align="center"|12.3594||align="center"|9.3||align="center"|-3.0593
|-
|align="center"|'''11''' ||align="center"|53.372||align="center"|56.3||align="center"|2.9328
|-
|align="center"|'''12'''||align="center"|37.5169||align="center"|43.8||align="center"|6.2831
|-
|align="center"|'''13'''||align="center"|33.6605||align="center"|43.5||align="center"|9.8395
|}



[[Image:Endo_nmr.png|thumb|left|400px|Fig.2c - Endo Camphor NMR]]
[[Image:EXO_nmr.png|thumb|right|400px|Fig.2d - Exo Camphor NMR]]
























The results were not in very good agreement with that of the literatures but it is generally out of literature by a consistent range of ~10ppm. This error may be attributed to wrong setting in setting the reference for the NMR peaks and the limitations of this calculation method itself. The largest deviation between the literature and the predicted carbon shift corresponds to the carbonyl carbon. This shift is abnormally high due to its high deshielding nature caused by oxygen and its rigid structure. Also, the larger the chemical shift, it is likely to incur a greater inaccuracy in the calculation.

=Reference=
William L. Jorgensen, Dongchul Lim, James F. Blake,''J. Am. Chem. Soc.'', '''1993''', ''115'' (7), pp ''2936–2942''.{{DOI|10.1021/ja00060a048}}

Brian Halton, Sarah G. G. Russell,''J. Org. Chem.'', '''1991,''' 56 (19), pp 5553–5556.{{DOI|10.1021/jo00019a015}}
<references/>