Sign inversions in optical rotation as a function of wavelength (ORD spectra)

December 9th, 2019

I have been discussing some historical aspects of the absolute configuration of molecules and how it was connected to their optical rotations. The nomenclature for certain types of molecules such as sugars and less commonly amino acids includes the notation (+) to indicate that the specific optical rotation of the molecule has a positive (rather than a negative) value. What is rarely mentioned is the implicit wavelength at which the rotation is measured. Historically polarimeters operated at the so-called sodium Fraunhofer D-line (588.995nm, hence the name [α]D) and only much more recently at the mercury e-line (546.073nm). The former was used for uncoloured organic molecules, since it was realised early on that colour and optical rotation did not mix well. Here I take a closer look at this aspect by constructing the hypothetical molecule shown below.

The rational behind this choice is that it is (a) based on indigo, which is deep blue in colour and (b) has a bridge of four methylene groups added to make it (axially) asymmetric. The calculated UV/Vis spectrum (ωB97XD/Def2-SVP/SCRF=water, FAIR DOI: 10.14469/hpc/6457) is shown below and you can see the very intense absorption at 535nm (corresponding to a visually blue colour).

The electronic circular dichroism version of this spectrum (simply the difference in absorbance between left and right polarised light instead of absorbance by unpolarised light) is shown below, and this form of chiroptical spectroscopy in large measure replaced the use of specific optical rotations as a means of assigning absolute configurations from the 1960s onwards. Note that the large peak at 535nm is replaced by a much smaller one (the Cotton effect) in the ECD spectrum.

Now I show the original optical rotation as a function of wavelength in 10nm increments. At 589nm ([α]D) it is negative (-1364°), but what on earth is going on at a wavelength of ~535nm, which as you can see above is the value of the first electronic excitation?

An expansion in 0.2nm increments shows more clearly what is happening. The negative value suddently shoots down to -1,200,000°, frankly an absurd value, before discontinuously reversing sign to a positive value of 75,000°. At exactly the value of the electronic absorption it is zero. Most people seeing this happen would conclude that the mathematics derived from the solutions of the quantum mechanical equations is resulting in an unphysical discontinuity. It is in fact the result of the behaviour of the electric and magnetic dipole moment vectors, and it CAN be seen experimentally, albeit never in quite such extreme form![1] The sign of the optical rotation CAN invert, but in this very strange manner whereby if it starts as negative, it first becomes infinitely negative before passing through zero and becoming infinitely positive and finally settling down to a normal positive value. The reason by the way why the “blip” in the ORD spectrum above is +ve, but -ve in the expansion below is “digital resolution”, with the top trace having too coarse a resolution to capture the detail. Now the reason why optical rotation measurement at 589nm becomes clear; it avoids any inversions caused by this effect for the majority of less-coloured molecules. However, if you do have a molecule that were to absorb at 589nm itself, the sodium D-line is the last wavelength you would want to use to measure its optical rotation!

So the notation (+) or (-) used to describe the sign of the specific rotation of a chiral molecule might give the misleading impression that it is a characteristic of the molecule at all wavelengths used to measure it. The rotation can change sign an impressive number of times as the wavelength changes! Does anyone know of any coloured pharmaceutical drugs that are available as pure enantiomers?  It would be fun to repeat the above on such a molecule.

References

  1. M.S. Andrade, V.S. Silva, A.M. Lourenço, A.M. Lobo, and H.S. Rzepa, "Chiroptical Properties of Streptorubin B: The Synergy Between Theory and Experiment", Chirality, vol. 27, pp. 745-751, 2015. https://doi.org/10.1002/chir.22486

What effect do explicit solvent molecules have on calculated optical rotation: D-("+")-Glyceraldehyde.

December 7th, 2019

In this series of posts on optical rotations, I firstly noted Kirkwood’s 1937 attempt to correlate the optical rotation of butan-2-ol with its absolute configuration. He had identified as essential knowing the relative orientation (the term conformation was not yet in common use) of the two methyl groups (the modern terms are gauche vs anti) and also that of the hydroxyl group, noting that anisotropy from this group could influence his result (he had assumed it was linear, or axially symmetric). I then looked at D-(+)-glyceraldehyde, noting that this species itself has a strongly negative rotation and that it is the hydrated diol that results in a positive rotation and hence the (+) designation. Here I take another look at this latter system to see what effect adding explicit water molecules to the unhydrated form of glyceraldehyde might have on its computed rotation, on the premise that strong hydrogen bonds can also contribute anisotropy to the system.

Firstly, here again are the computed results for glyceraldehyde on its own, albeit encased in a continuum solvent field for water (SCRF=water). At 303K and 589nm, the computed rotation is -193° compared to -147° inferred from the population of the aldehyde.[1]

Here are the new values (FAIR data DOI: 10.14469/hpc/6445) obtained by adding an explicit water molecule to the original conformations. This also introduces extra conformations, some of which are included below. This reduces the calculated value by ~25° to improve the agreement with measurement (-147°). The value at 436nm is -411 (calc), -380 (obs).

So a 25° correction is not entirely insignificant, but does not change the overall conclusion that the optical rotation of D-(+) glyceraldehyde is (-). How might the model be improved further?

  1. Adding more water molecules, in theory until a limit is reached where further anisotropy is not added to the model. But the major disadvantage is that each extra water molecule increases the conformational space to be explored. With say seven added water molecules, there are probably 100s of conformations that would have to be searched, now a major undertaking.
  2. Another interesting avenue to explore is the temperature dependence of the optical rotation. The experimental values are shown below. This is due to the change in Boltzmann populations as a function of temperature.

    At 343K, the original calculations without inclusion of a water molecule reduce the calculated rotation changes from -193 to -182.5, or about 10°. The observed value is a change of 39°. However, with the model including an extra water molecule, the value changes from -168.5 to -171°. So it might well be that reproduction of the temperature effects will require more water molecules added to the model.

 

References

  1. M. Fedoroňko, "Optical activity of D-glyceraldehyde in aqueous solutions", Collection of Czechoslovak Chemical Communications, vol. 49, pp. 1167-1172, 1984. https://doi.org/10.1135/cccc19841167

What effect do explicit solvent molecules have on calculated optical rotation: D-(“+”)-Glyceraldehyde.

December 7th, 2019

In this series of posts on optical rotations, I firstly noted Kirkwood’s 1937 attempt to correlate the optical rotation of butan-2-ol with its absolute configuration. He had identified as essential knowing the relative orientation (the term conformation was not yet in common use) of the two methyl groups (the modern terms are gauche vs anti) and also that of the hydroxyl group, noting that anisotropy from this group could influence his result (he had assumed it was linear, or axially symmetric). I then looked at D-(+)-glyceraldehyde, noting that this species itself has a strongly negative rotation and that it is the hydrated diol that results in a positive rotation and hence the (+) designation. Here I take another look at this latter system to see what effect adding explicit water molecules to the unhydrated form of glyceraldehyde might have on its computed rotation, on the premise that strong hydrogen bonds can also contribute anisotropy to the system.

Firstly, here again are the computed results for glyceraldehyde on its own, albeit encased in a continuum solvent field for water (SCRF=water). At 303K and 589nm, the computed rotation is -193° compared to -147° inferred from the population of the aldehyde.[1]

Here are the new values (FAIR data DOI: 10.14469/hpc/6445) obtained by adding an explicit water molecule to the original conformations. This also introduces extra conformations, some of which are included below. This reduces the calculated value by ~25° to improve the agreement with measurement (-147°). The value at 436nm is -411 (calc), -380 (obs).

So a 25° correction is not entirely insignificant, but does not change the overall conclusion that the optical rotation of D-(+) glyceraldehyde is (-). How might the model be improved further?

  1. Adding more water molecules, in theory until a limit is reached where further anisotropy is not added to the model. But the major disadvantage is that each extra water molecule increases the conformational space to be explored. With say seven added water molecules, there are probably 100s of conformations that would have to be searched, now a major undertaking.
  2. Another interesting avenue to explore is the temperature dependence of the optical rotation. The experimental values are shown below. This is due to the change in Boltzmann populations as a function of temperature.

    At 343K, the original calculations without inclusion of a water molecule reduce the calculated rotation changes from -193 to -182.5, or about 10°. The observed value is a change of 39°. However, with the model including an extra water molecule, the value changes from -168.5 to -171°. So it might well be that reproduction of the temperature effects will require more water molecules added to the model.

 

References

  1. M. Fedoroňko, "Optical activity of D-glyceraldehyde in aqueous solutions", Collection of Czechoslovak Chemical Communications, vol. 49, pp. 1167-1172, 1984. https://doi.org/10.1135/cccc19841167

The (+) in D-(+)-glyceraldehyde means it has a positive optical rotation? Wrong!

December 6th, 2019

Text books often show the following diagram, famously consolidated over many years by Emil Fischer from 1891 onwards. At the top sits D-(+)-glyceraldehyde, to which all the monosaccharides below are connected by painstaking chemical transformations.

In this notation, D (for all these structures) indicates the absolute configuration of the series, deriving from placing the last highest priority group before an achiral carbon at the bottom of the representation (the hydroxyl group) to the right of the diagram (= dextro = D). The mirror images of all these species would be designated L (since the hydroxide group would now be on the left). The challenge emerges in connecting this absolute configuration to some measurable property, and during the late 19th century, the only feasible measurement would be the sign of the optical rotation. It was decided that since no method existed at the time to connect absolute configuration to the sign of optical rotation, D-glyceraldehyde (nowadays more commonly referred to by its CIP notatation (R)-glyceraldehyde), would be by convention connected to the enantiomer giving a positive rotation, hence D-(+)-glyceraldehyde. This species has indeed mostly been reported as having a small positive rotation of e.g. +9.2[1] in water. If a future procedure were to confirm that the D-enantiomer had a positive rotation, Fischer’s guess would be proven correct. Possibly.

I quickly mention that this experiment was of course carried out in 1951, as noted in the previous post[2], but on a rubidium tartrate salt and not on glyceraldehyde itself (which is a liquid). Since the two could be connected by chemical synthesis, Fischer’s guess was indeed confirmed as correct. This proof rapidly made D-(+)-glyceraldehyde less relevant as the tentpole of absolute configuration (in sugars) and people stopped worrying about it. But the assertion that the D/(R) enantiomer has a positive rotation by indicating (+) continues in all the diagrams related to the above that I have come across. So is it correct?

Well, an earlier attempt in 1937 to assign optical rotation to absolute configuration had been attempted by Kirkwood[3] using quantum mechanical theory to make that connection. I thought I might apply the modern version of that approach to (R)-glyceraldehyde, using the same procedure I had tried earlier. A conformational analysis of the molecule, followed by dispersion-corrected calculation of free energies and then of the optical rotations gives the following table. Data is at 10.14469/hpc/6436These calculations show that the rotation is strongly negative! Can they be in serious error? Perhaps however there are two different species in solution. Indeed, the equilibrium below is shown to favour the rhs for aldehydes, as shown by the following concentration profile as a function of temperature.[4] This shows that at 30°, only 4% of this molecule is in the aldehyde form.

The article also reports optical rotations, at three wavelengths. I should note here that the notation (+) means only the rotation at 589nm (yes, you can get sign inversions on moving from one wavelength to another![5]). At 30° the observed rotation for D-glyceraldehyde is indeed strongly negative, whilst it is the hydrated form that is moderately (+).

There is a discrepancy of ~47° between these measured values and the calculated one at 589nm. I hope to explore the origins of this error in a separate post.

So we see by both experiment and theory that D-(+)-glyceraldehyde, if the meaning of that (+) is to be upheld, should really be D-(-)-glyceraldehyde. Note that many of the sugars shown in the top diagram have (-) as well as (+), and so there is no reason that glyceraldehyde should be any different.

References

  1. H.J. Lamble, M.J. Danson, D.W. Hough, and S.D. Bull, "Engineering stereocontrol into an aldolase-catalysed reaction", Chemical Communications, pp. 124, 2005. https://doi.org/10.1039/b413255f
  2. J.M. BIJVOET, A.F. PEERDEMAN, and A.J. van BOMMEL, "Determination of the Absolute Configuration of Optically Active Compounds by Means of X-Rays", Nature, vol. 168, pp. 271-272, 1951. https://doi.org/10.1038/168271a0
  3. J.G. Kirkwood, "On the Theory of Optical Rotatory Power", The Journal of Chemical Physics, vol. 5, pp. 479-491, 1937. https://doi.org/10.1063/1.1750060
  4. M. Fedoroňko, "Optical activity of D-glyceraldehyde in aqueous solutions", Collection of Czechoslovak Chemical Communications, vol. 49, pp. 1167-1172, 1984. https://doi.org/10.1135/cccc19841167
  5. M.S. Andrade, V.S. Silva, A.M. Lourenço, A.M. Lobo, and H.S. Rzepa, "Chiroptical Properties of Streptorubin B: The Synergy Between Theory and Experiment", Chirality, vol. 27, pp. 745-751, 2015. https://doi.org/10.1002/chir.22486

Prediction preceding experiment in chemistry – how unlucky was John Kirkwood?

November 30th, 2019

Some areas of science progressed via very famous predictions that were subsequently verified by experiments. Think of Einstein and gravitational waves or of Dirac and the positron. There are fewer well-known examples in chemistry; perhaps Watson and Crick’s prediction of the structure of DNA, albeit based on the interpretation of an existing experimental result. Here I take a look at a what if, that of John Kirkwood’s prediction of the absolute configuration of a small molecule based entirely on matching up the sign of a measured optical rotation with that predicted by (his) theory.

The confirmation that Emil Fischer’s 1891 proposed convention for the absolute configuration of sugars was in fact correct was famously made by Bijvoet in 1951 using crystallography.[1] I first told this story in 2012, noting that Kirkwood apparently made his seminal contribution a year later in 1952[2] using his quantum mechanical theory of optical rotation to independently come up with the same result. Nowadays he rarely gets the credit for solving the problem of absolute configuration. But wait, Kirkwood’s first stab at solving this problem in fact came in 1937,[3] a full 14 years before Bijvoet’s famous result (for which incidentally the Nobel prize was not awarded). 

I have been asked to talk about this story at a Historical meeting of the Royal Society of Chemistry in March 2020, and for this purpose thought I should take a closer look at Kirkwood’s 1937 article. In it, he sets out his quantum mechanical theory of optical rotation. Remember that in that era, there was no recourse to computers and solving the required (heavily approximated) equations had to be done entirely using mechanical calculators. Kirkwood chooses to analyse the following molecule;

I have redrawn it below in more modern form (and name). The difference between the two is in the notations and (R). The former relates to the sign of the optical rotation [α]D where d stands for dextrorotation, or clockwise and is also often represented by (+) and being the sign of the measured rotation. (R) is the modern notation for the absolute configuration shown by Kirkwood in his diagram (Fig. 1).

He is asserting below that the enantiomer of butan-2-ol with a measured rotation of 13.9° has (in modern notation) the absolute configuration (R) because his calculations predict somewhere between 9.5° and 21.9° for this specific three-dimensional geometry. So why is Kirkwood not lauded for solving this problem in 1937? Well, because we now know that (R)-butan-2-ol has a negative rotation of -13.9°![4]

Kirkwood is however very aware of the potential problems with his approach. In a nutshell, conformation! In particular, the conformers resulting from rotation about the central C-C bond and especially the C-O bond, where for the purposes of his theory he assumed axial symmetry about that bond.

and

Now, in 1937 the area of such conformational analysis was hardly known; only in 1948[5] would Barton first put it firmly on the map (and win the Nobel prize for this work). So, in attempting his connection between and (R), Kirkwood was in a sense far too ahead of his time. It worth asking what modern quantum mechanical theory makes of this problem and does it cast any light on why Kirkwood actually got his assignment wrong (in 1937,[3] although he WAS correct in 1952[2]).

  1. Firstly, I carried out a comprehensive search of the rotamers about the C-C and C-O bond using molecular mechanics (a method first introduced by Barton in 1948[5]) and using the MMFF94 forcefield. This identifies seven distinct conformations arising from rotations about these two bonds. Some warning signs area aready present; these seven are bounded by an energy of only 1.2 kcal/mol! All are likely to have a significant Boltzmann population.
  2. Next, to ramp up the level of theory to density functional quantum mechanics, at the B3LYP+GD3+BJ/Def2-TZVPP/SCRF=diethyl ether level (FAIR Data: 10.14469/hpc/6367) and at the minimum energy geometry for each conformation, an optical rotation is calculated (at the ωB97XD/Def2-TZVPP/SCRF=diethyl ether) level. Whilst not the highest practical level possible nowadays, it far exceeds in accuracy what Kirkwood had at his disposal in 1937. The results are at 10.14469/hpc/6367 and available as a spreadsheet if you want to adapt this for your own needs.

What can we conclude?

  1. All seven conformations have a significant population (at 298K). The ordering, with one small exception, is the same for both Molecular and Quantum mechanics. The relative free energies span a slightly larger range than the steric energies obtained from molecular mechanics (1.7 kcal/mol) but all have a population of >6%.
  2. Three conformations have a negative or (-) predicted rotation, and four are positive (+).
  3. When weighted by population, the overall predicted rotation is -18°, which compares well with that observed (~-13°).
  4. Kirkwood, who was not able to include all seven conformations in his analysis, was deeply unlucky that his particular choices/assumptions of conformations happened to have (+) rotations. But he was very much aware that this result could happen (although he does underestimate this in concluding that his tentative result is “probably accurate”. To be fair, if he had been more realistic, the referees might well have rejected his article!).

So the “what if“.

  • If Kirkwood had chosen a conformationally simpler molecule which did not have as many as seven populated conformations, he may well have got his prediction correct (and for mostly the right reasons!). 
  • But he has to be given lots of credit for recognising that optical rotations can be sensitive to conformational analysis. In this regard he could be regarded as one of the early fathers of that entire (Nobel prizing winning) field.
  • He was 14 years ahead of the eventual unambiguous experiment that verified the Fischer convention. Of course he would have needed to correlate the absolute configuration of butan-2-ol with those of both sugars and amino acids using chemical transformations. In 1937, this may well have been quite synthetically challenging (but of course perhaps this correlation may have been actually known at the time. Does anybody reading this know?)
  • But given that two other discoveries, both of which won the Nobel prize (the structure of peptides and the structure of DNA), depended on knowing with certainty the absolute configurations of amino acids and sugars respectively, Kirkwood’s method could be argued directly impacted upon no less than three Nobel prizes within two decades of his initial work.
  • Remember that Watson and Crick “predicted” that the DNA helix is right-handed, as it happens on the basis of a single “short” H…H contact in their model which apparently disfavoured a left-handed helix. Boy was that a lucky guess (since that conclusion cannot be sustained nowadays on the basis of short H…H contacts). And that Pauling, in his own initial structures suggesting an α-helix in some proteins, predicted (wrongly) that the helix was left handed.  

So I think that yes, Kirkwood was pretty unlucky in his 1937 effort. And by 1952 (when he was correct), the opportunity for widespread recognition for this work and perhaps even a Nobel prize, had passed.


Stereochemical notation has suffered from some measure of confusion over the years. Much of that confusion was cleared up with the introduction of the CIP rules, but historical vestiges remain. Thus d was originally used by Fischer himself to indicate configuration using the sense of direction on his diagram, but by others (including Kirkwood in 1937) to indicate the sense of direction of polarised light, an entirely different property. Eventually, the configurational sense became distinguished from the rotational light sense by capitalising the former (which of itself can still lead to confusions). Nowadays, the configuration of an entire molecule tends to be described by specifying the absolute configuration of all the asymmetric units using the CIP formalism, whilst Fischer’s formalism (now rationalised D/L), is only applied to sugars and cannot be used generally for other molecules. If you want to explore the temperature dependence of the Boltzmann populations and hence the predicted change in rotation with temperature, do please download the spreadsheet and try it out for yourself! The elapsed time for these 14 calculations took about 2 hours. The exhaustive molecular mechanics calculations took <20 seconds.

References

  1. J.M. BIJVOET, A.F. PEERDEMAN, and A.J. van BOMMEL, "Determination of the Absolute Configuration of Optically Active Compounds by Means of X-Rays", Nature, vol. 168, pp. 271-272, 1951. https://doi.org/10.1038/168271a0
  2. W.W. Wood, W. Fickett, and J.G. Kirkwood, "The Absolute Configuration of Optically Active Molecules", The Journal of Chemical Physics, vol. 20, pp. 561-568, 1952. https://doi.org/10.1063/1.1700491
  3. J.G. Kirkwood, "On the Theory of Optical Rotatory Power", The Journal of Chemical Physics, vol. 5, pp. 479-491, 1937. https://doi.org/10.1063/1.1750060
  4. A.Z. Gonzalez, J.G. Román, E. Gonzalez, J. Martinez, J.R. Medina, K. Matos, and J.A. Soderquist, "9-Borabicyclo[3.3.2]decanes and the Asymmetric Hydroboration of 1,1-Disubstituted Alkenes", Journal of the American Chemical Society, vol. 130, pp. 9218-9219, 2008. https://doi.org/10.1021/ja803119p
  5. D.H.R. Barton, "83. Interactions between non-bonded atoms, and the structure of cis-decalin", Journal of the Chemical Society (Resumed), pp. 340, 1948. https://doi.org/10.1039/jr9480000340

The Structure of Tetrodotoxin as a free base – with a better solvation model.

November 26th, 2019

In the previous post, I discussed the structure of the free base form of tetrodotoxin, often represented as originally suggested by Woodward[1] below in an ionic form:

Quantum calculations suggested that this form was higher in energy than neutral forms devoid of the zwitterionic charge separation in a relatively non polar solvent such as chloroform. For this, a so-called continuum solvation model was used. But even chloroform is capable of forming rather strong O…H-C hydrogen bonds, and these specific isotropic interactions are not well modelled using a continuum solvent; you need to include the specific hydrogen bonds to do that. So here are two better models, the first including one or two chloroform molecules (in continuum chloroform) and the second three water molecules (in continuum water). The (FAIR) data for these results are available at DOI: 10.14469/hpc/6278 and you can view the 3D models by clicking on the images below.

System ΔΔG298, + 1CHCl3 ΔΔG298, + 2CHCl3 ΔΔG298, + 3H2O
Ionic form 0.0 0.0 0.0
non-ionic form -7.1 -4.9 +1.5

Adding additional solvent molecules in a “non-stochastic” manner is clearly an approximation. For example, three water molecules added to the neutral (non-zwitterionic) form could take up residency in many different ways, given the number of oxygen atoms present in tetrodotoxin. The energies reported above are for the lowest energy forms of the two that I located; I did not investigate more possibilities. However, the ionic form, with the water molecules directly hydrogen bonding with the ionic oxygen are more certain.

Ionic form: Click image to view 3D model

Neutral form: Click image to view 3D model

So we now see that Woodward’s original proposal for a charge-separated ionic form for the free base of tetrodotoxin may indeed be accurate, but only for polar solutions such as water. In rather less polar solutions such as chloroform it is probably not ionic. Clearly the stability also depends on the number of solvent molecules included in the model. The rather large chloroform molecule is less likely to accumulate around the ionic oxygen in numbers >>2 (and including more makes the calculation much slower), but already with just three water molecules the ionic form becomes the more stable. Also, e.g. two chloroforms or merely three waters are insufficient to form any sort of stabilizing “bridge” connecting the two ionic centres. So these results must be treated with a little caution.

It is always rather risky to bet against Woodward’s chemical intuitions and insights (as for example Robinson found when debating the structure of strychnine[2] with him). He may well have been correct with the ionic structure of tetrodotoxin as well (at least in aqueous solutions)!

References

  1. R.B. Woodward, "The structure of tetrodotoxin", Pure and Applied Chemistry, vol. 9, pp. 49-74, 1964. https://doi.org/10.1351/pac196409010049
  2. R.B. Woodward, M.P. Cava, W.D. Ollis, A. Hunger, H.U. Daeniker, and K. Schenker, "THE TOTAL SYNTHESIS OF STRYCHNINE", Journal of the American Chemical Society, vol. 76, pp. 4749-4751, 1954. https://doi.org/10.1021/ja01647a088

The Structure of Tetrodotoxin as a free base.

November 9th, 2019

The notorious neurotoxin Tetrodotoxin is often chemically represented as a zwitterion, shown below as 1. This idea seems to originate from a famous article written in 1964 by the legendary organic chemist, Robert Burns Woodward.[1] This structure has propagated on to Wikipedia and is found in many other sources.
With the elegance and the unique style that is typical Woodward, his article is a tour de force because of the way in which he deploys a large armoury of spectroscopic (X-ray crystal, NMR, IR) as well as physicochemical (pKa) tools to infer this structure; an approach that has been subsequently widely emulated. The article a well worth a read for the elegant logic that slowly builds to a climax on page 73 (sic!) of the article, when he unveils his final structure (XXXVIII, or 38). The lecture(s) from which the article is apparently derived must have been one hell of an occasion.

One technique not available in 1964 to Woodward was quantitative quantum calculation of the molecular free energies. This property is now routinely computable, but is still only rarely used for structural problems even today. Here I add this property to Woodward’s collection (FAIR data DOI: 10.14469/hpc/6278, method=ωB97XD/Def2-TZVPP, SCRF=chloroform).

My initial hypothesis was based on the observation of a rather large separation of the charges in the zwitterion. The surface below is the computed molecular electrostatic potential (value 0.1 au), which indeed shows large charge separation. The computed dipole moment is 18D.

MEP isosurface; orange = -ve,  yellow = +ve.

Charges by and large do not much like being separated, even though of course the actual separation is often misleadingly indicated by simple Lewis structures such as the above. So I started by moving one proton to produce 2 to see if that improved the free energy, by reducing charge separation. The dipole moment was indeed smaller at 14.6D (although the energy was not lower).

I then included three neutral forms (35) where the nominal charge separation is eliminated entirely.

Structure Diagram

ΔΔG

kcal/mol

Dipole moment
1 +10.4 18.0
2 +21.1 14.6
3 0.0 3.9
4 +1.8 6.0
5 +4.2 6.5

This computed model suggests 3 is the predominant species present (in chloroform solutions), although of course all the species are in a prototropic equilibrium. 3 also happens to have the smallest calculated dipole moment. Species 5 was also a front-runner in 1964[2] (as 1c in that article). Certainly the zwitterion 1 as suggested by Woodward, with its large degree of charge separation, is unlikely to feature and 5 is also less likely to be the dominant species.

The computed free energy for such molecules takes < 1 day of elapsed time to produce, and so I here argue that this should be a mandatory reported property for such structural problems.


The crystal structure is recorded[3] for the protonated species, as the HCl or HBr salt. These of course do not indicate where the protons are for the deprotonated neutral base. I only ever attended one lecture by Woodward. It lasted the typical two hours and was indeed hugely memorable in several regards. I also note the parsimony of stereochemical notation (just one dashed bond), presumably on the grounds that omitting such notation does not actually result in ambiguity. Disambiguation does depend of course on perceived “hidden line” removal.

References

  1. R.B. Woodward, "The structure of tetrodotoxin", Pure and Applied Chemistry, vol. 9, pp. 49-74, 1964. https://doi.org/10.1351/pac196409010049
  2. T. Goto, Y. Kishi, S. Takahashi, and Y. Hirata, "Further studies on the structure of tetrodotoxin", Tetrahedron Letters, vol. 5, pp. 779-786, 1964. https://doi.org/10.1016/0040-4039(64)83035-5
  3. A. Furusaki, Y. Tomiie, and I. Nitta, "The Crystal and Molecular Structure of Tetrodotoxin Hydrobromide", Bulletin of the Chemical Society of Japan, vol. 43, pp. 3332-3341, 1970. https://doi.org/10.1246/bcsj.43.3332

Does Kekulene have Kekulé vibrational modes? Yes!

October 19th, 2019

Increasingly, individual small molecules are having their structures imaged using STM, including cyclo[18]carbon that I recently discussed. The latest one receiving such treatment is Kekulene.[1]

As with cyclo[18]carbon, the point of interest was which of the two resonance structures shown below most closely resembled the measured structure. The one on the left has six cyclohexatriene motifs (red), also known as Clar rings (although I have argued they could also be called Armstrong rings), with six largely isolated double bonds. These each follow the 4n+2 π-electron aromaticity rule (n=1). The structure on the right has instead two fully conjugated larger rings, the inner one blue following 4n+2, n=4 (18) and the outer one (magenta) following 4n+2, n=7 (30).

Now, the STM result[1] clearly came down in favour of the left hand structure, with a sextet of six-rings, rather than the larger annulenes implied by the right hand structure. A previously reported crystal structure[2] came to a similar conclusion, and the species is shown as fully planar. There is just one oddity about this structure. The distance between the inner six hydrogens is reported as 1.938Å. In reality, it must be shorter, since the C-H bond lengths of these hydrogens are reported as 1.00Å whereas a more realistic value is ~1.09Å. Applying this well-known correction means that the non-bonded H-H distance in Kekulene is actually closer to 1.83Å, which is unusually close. Recollect a similar issue arose in the imaging of purportedly planar tetraphenylporphrin, in which H…H distances as short as 0.8Å were implied in a published article.[3]

Previously, I have shown that conjugated π-electron rings can exhibit Kekulé vibrational modes for all these ring sizes, and so the question arises as to the nature of any Kekulé vibration in Kekulene itself. This vibration corresponds to a motion which creates two alternative cyclohexatriene structures (for benzene itself) from the symmetrical molecule. For benzene, this mode is real, ie it takes energy to distort benzene into a cyclohexatriene. For larger annulenes, this vibration is imaginary and hence corresponds to a transition state between two bond alternating lower energy forms. So the Kekulé vibration in Kekulene should be distinctly different for the two structures shown above. I am using two different density functional methods, B3LYP+GD3BJ/Def2-TZVPP and ωB97XD/Def2-TZVPP (FAIR DOI: 10.14469/hpc/6226), which probably straddle reality in terms of their propensity for bond length alternation (BLA) in larger annulene rings.

The results show that a fully planar molecule with D6h symmetry is actually a transition state in which three of these six inner hydrogen atoms move up and three down, giving a non-planar molecule with D3d symmetry and increasing the H…H separation. However, both DFT methods also agree that there is just one Kekulé vibration, which is real (ν(A2g) 1389/B3LYP or 1395/ωB97XD cm-1) and corresponds to the below.

B3LYP

wB97XD

All six 6-membered Clar/Armstrong rings are executing their own Kekulé vibration, but in perfect synchrony with each other. There is almost no activity in the six passive double bonds. There is certainly no extended Kekulé vibration similar to the one shown here which involves the entire larger ring, whether in the inner or outer periphery.

So these results match nicely with the initial crystal structure and the recent STM imaging one Curiously, no-one had ever mentioned the Kekulé vibration in this species, and so this omission is rectified here.


Update: Some variations on the theme above.

  1. The inner ring of hydrogens is replaced by carbons
    This has the following Kekule modes, a B2u mode centered on just the inner six carbons (red, ν 1423 cm-1)

    and a degenerate E2g mode (blue, ν 1383 cm-1) for the six outer Clar rings.
  2. The inner ring is replaced by borons.

References

  1. I. Pozo, Z. Majzik, N. Pavliček, M. Melle-Franco, E. Guitián, D. Peña, L. Gross, and D. Pérez, "Revisiting Kekulene: Synthesis and Single-Molecule Imaging", Journal of the American Chemical Society, vol. 141, pp. 15488-15493, 2019. https://doi.org/10.1021/jacs.9b07926
  2. H.A. Staab, F. Diederich, C. Krieger, and D. Schweitzer, "Cycloarenes, a New Class of Aromatic Compounds, II. Molecular Structure and Spectroscopic Properties of Kekulene", Chemische Berichte, vol. 116, pp. 3504-3512, 1983. https://doi.org/10.1002/cber.19831161022
  3. J. Lee, K.T. Crampton, N. Tallarida, and V.A. Apkarian, "Visualizing vibrational normal modes of a single molecule with atomically confined light", Nature, vol. 568, pp. 78-82, 2019. https://doi.org/10.1038/s41586-019-1059-9

Catalytic Mitsunobu reaction.

October 9th, 2019

If, as a synthetic chemist, you want to invert the configuration of an alcohol in which the OH group is at a chiral centre, then the Mitsunobu reaction has been a stalwart for many years. Now a catalytic version has been published, [1] along with a proposed mechanism. Here I apply computation as a reality check to see what the energetics of this mechanism might be.

The chosen computational procedure was B3LYP+GD3BJ dispersion, Def2-TZVPP basis and SCRF=toluene as solvent and for the alcohol, R1=R2=Me. Data is at DOI: 10.14469/hpc/6186

Molecule ΔΔG, kcal/mol
1 0.0
2 25.4
3 15.0
TS 36.8
4 +4.5

In the mechanism above, the proton transfers have not been included in the modelling, with the presumption that the step involving the inversion at carbon is rate limiting.

Overall, the reaction is slightly endoenergic, and this is supported by the need to azeotropically remove water from the reaction. The activation barrier for TS from species 1 however is perhaps somewhat high for a reaction occuring in refluxing toluene. Annoyingly so, since high barriers are sometimes indicative that the mechanism modelled is missing some key aspect. So there must be (a small?) possibility that a more complex mechanism might operate in this catalytic cycle, or that a systemic error in the  DFT (density functional theory) approach might be present.

References

  1. R.H. Beddoe, K.G. Andrews, V. Magné, J.D. Cuthbertson, J. Saska, A.L. Shannon-Little, S.E. Shanahan, H.F. Sneddon, and R.M. Denton, "Redox-neutral organocatalytic Mitsunobu reactions", Science, vol. 365, pp. 910-914, 2019. https://doi.org/10.1126/science.aax3353

Bond length alternation (BLA) in large conjugated rings: an (anti-aromatic) update.

October 3rd, 2019

In the previous post, I looked at a class of molecule known as hexaphyrins, inspecting bond length alternation (BLA) at the so-called meso position, the carbon atom joining two pyrrole rings. A search of the difference in bond lengths at this position had shown two significant clusters of crystal structures.
Molecules in the bottom left of this diagram shows little or no bond length alternation. The right middle shows another cluster with more extreme (and unequal) bond length alternation. I have selected one molecule from this cluster, EGIJEK and it differs from EGIHUY in having four NH units in the ring, whereas the latter has only two.

EGIJEK

EGIHUY

The effect this has is illustrated below. A pyrrole with an NH group contributes two π-electrons to the conjugated periphery whereas a pyrrole with just an N contributes only one. Thus the four NH rings in EGIJEK contribute overall two more π-electrons to the periphery than EGIHUY with just two NH rings. This increases the π-electron count from 26 to 28, driving EGIJEK from a 4n+2 into a 4n π-electron count (n=7) and making it formally anti-aromatic. This in turn induces strong bond length alternation (as for example in cyclobutadiene).

Here is a reprise of the bond length table shown in the previous post, but for EGIJEK.

Meso distances, Å abs(Δr)
EGIJEK crystal, Ci symmetry, DOI: 10.5517/ccrts2d
1,43657 1.37034 0.06623
1.35661 1.44895 0.09234
1.42356 1.37287 0.05069
B3LYP+GD3BJ/Def2-SVPP (FAIR DOI: 10.14469/hpc/6194)
1.44019 1.38259 0.0576
1.37732 1.44481 0.06749
1.43261 1.38568 0.04693
ωB97XD/Def2-SVPP (FAIR DOI: 10.14469/hpc/6194 )
1.44575 1.37537 0.07038
1.36214 1.45508 0.09294
1.44221 1.36958 0.07263

Of the two functionals used in the calculations, the ωB97XD form slightly over-estimates the BLA, with the B3LYP slightly under-estimating it. This seems to tally with the earlier observations made for cyclo[n]carbons.

Perhaps these sorts of molecules might form useful reality checks on calculating bond length alternations in large ring cyclic conjugated molecules.