Updating a worked problem in conformational analysis. Part 2: an answer.

The previous post set out a problem in conformational analysis. Here is my take, which includes an NCI (non-covalent interaction) display as discussed in another post.

The lowest energies of the four diastereomers A-D, each in two conformations (1/2) were calculated at the ωB97D/6-311G(d,p)/SCRF=ethanol level, and are shown here relative to A1 (kcal/mol) as free energies. The values of ΔΔG of each pair (relative to the lower energy conformation) are converted to relative concentrations using ΔΔG = -RT Ln K and these are shown in blue. Those conformations where a C-H bond is aligned anti-periplanar (app) with the C-NMe3 bond are highlighted in red. Various aspects of the questions posed include:

  1. Those conformations with two app combinations are likely to result in a mixture of two alkenes E and F (shown in the previous post), and likewise those with only one app alignment will give a single alkene.
  2. The concentration of the neomenthyl system (C) favours the conformer C1, which is also the one with app allignments. This will react readily. The menthyl system A has only a low relative concentration of A2, and so will react slowly. The experiments show the rate is about 10 times less than C, whereas the ratio of relative concentrations is larger. This may be due to difference in the transition state for the elimination, or possibly effects due to the presence of a counterion.
  3. The concentration of D1 is the lowest of the four pairs of conformational isomers. This might be due to buttressing, in which the relative proximity of the NMe3 and iPr groups forces the former into closer proximity with the axial methyl than is found in A2, the other conformer with an unfavourable 1,3-diaxial interaction.
  4. The relative energies of most of the pairs of combinations of conformers/diasteromers are interesting, and I leave it to you to play with them and draw conclusions (or detect mysteries).

I thought I would finish by adding an NCI analysis using a wavefunction computed for ωB97D/6-31G(d,p)/SCRF=ethanol. The green surfaces are indicative of regions where non-covalent interactions are occurring (which include electrostatic and van der Waals effects). The colour scheme is set to blue(ish) to indicate more attractive NCIs (such as strong hydrogen bonds), green to indicate weak(ish) NCIs, with a red or orange tinge indicating a repulsive  NCI.

A2 diaxial. Click for 3D.

B2 Axial-equatorial. click for 3D.

A1. Di-equatorial. Click for 3D.

B1.Di-equatorial. click for 3D.

C2. Equatorial-axial. Click for 3D.

D2. Equatorial-axial. Click for 3D.

C1 SRR. Axial-equatorial. Click for 3D.

D1 SSS. Axial-equatorial. Click for 3D.

Amongst the many features visible in these plots is the strong 1,3-diaxial NCIs seen for A2 and particularly D1 and the conformational variability in the NCI between the iPr and NMe3 groups.

Henry Rzepa

Henry Rzepa is Emeritus Professor of Computational Chemistry at Imperial College London.

View Comments

Recent Posts

Internet Archeology: reviving a 2001 article published in the Internet Journal of Chemistry.

In the mid to late 1990s as the Web developed, it was becoming more obvious…

1 month ago

Detecting anomeric effects in tetrahedral carbon bearing four oxygen substituents.

I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…

1 month ago

Data Citation – a snapshot of the chemical landscape.

The recent release of the DataCite Data Citation corpus, which has the stated aim of…

2 months ago

Mechanistic templates computed for the Grubbs alkene-metathesis reaction.

Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…

2 months ago

3D Molecular model visualisation: 3 Million atoms +

In the late 1980s, as I recollected here the equipment needed for real time molecular…

3 months ago

The Macintosh computer at 40.

On 24th January 1984, the Macintosh computer was released, as all the media are informing…

3 months ago