Infra-red spectroscopy of molecules was introduced 110 years ago by Coblentz[1] as the first functional group spectroscopic method (” The structure of the compound has a great influence on the absorption spectra. In many cases it seems as though certain bonds are due to certain groups.“). It hangs on in laboratories to this day as a rapid and occasionally valuable diagnostic tool, taking just minutes to measure. Its modern utility rests on detecting common functional groups, mostly based around identifying the nature of double or triple bonds, and to a lesser extent in differentiating between different kinds of C-H stretches[2] (and of course OH and NH). One common use is to identify the environment of carbonyl groups, C=O. These tend to come in the form of aldehydes and ketones, esters, amides, acyl halides, anhydrides and carbonyls which are part of small rings. The analysis is performed by assigning the value of the C=O stretching wavenumber to a particular range characteristic of each type of compound. Thus ketones are said to inhabit the range of ~1715-1740 cm-1 and simple esters come at ~1740-1760 cm-1, some 20-30 cm-1 higher. Here I try to analyse how this difference arises.
The analysis is based on trying to understand how the components of an ester interact with each other, and in particular how the alkyl oxygen interacts with the carbonyl group. Three electronic interactions in particular can be focused on (below). The first two of these weaken the C=O bond; the last strengthens it. So which effect wins out?
I will start with computational models, which have the advantage that one can dissect how the vibrations arise. The first two rows show a comparison of the experimental gas phase values[3] with a standard “medium level” ωB97XD/6-311G(d,p) calculation. The discrepancy amounts to ~100-114 cm-1.
The carbonyl stretch in esters and ketones | ||
---|---|---|
Method: | Ester | Ketone |
Expt (gas phase)[3] | 1761‡ | 1737‡ |
Harmonic ωB97XD/6-311G(d,p) | 1860 | 1851 |
Anharmonic ωB97XD/6-311G(d,p) | 1832† | 1828† |
Harmonic ωB97XD/aug-cc-pvQZ | 1836 | 1831 |
Harmonic CCSD(T)/6-311G(d,p) | 1826 | 1792 |
Corrected CCSD(T)/6-311G(d,p) | ~1774 | ~1749 |
Expt (gas phase) | 1761 | 1737 |
Reduced CCSD(T)/6-311G(d,p) | 1764 | 1743 |
There are several possible causes for such errors:
So to really get to the root of why an observed ester carbonyl stretch is higher than that of the equivalent ketone, we have to get a handle on these effects above.
Now that we can assess the accuracy of our computational methods, we need to try to relate the results to the C=O bond itself. Does turning a ketone into an ester really make it stronger? To directly compare the C=O bond of two different molecules, we need to eliminate the effects of mixing the C=O normal stretching mode with similar energy modes arising from other parts of the molecule. A simple way of estimating this is to set the mass of all but two of the atoms to a very small value (0.00001), leaving only the masses of the C and O as normal; this is shown as a reduced frequency in the table above. The harmonic CCSD(T)/6-311G(d,p) C=O “pure” mode reduces to 1764 for methyl ethanoate and 1743 cm-1 for propanone. So after all of this, at least we now know that the force constant for the C=O stretch really is stronger for an ester. The green arrows seem to win out over the blue/red ones.
One calculation too many? The (Wiberg) bond order for the C=O bond can be derived from the wavefunctions. Its value is 1.635 for ester, and 1.681 for ketone (CCSD/6-311G(d,p)) or 1.766/1.848 (ωB97XD/aug-cc-pvQZ). This is the opposite to that inferred from the carbonyl stretch, and hence favours the blue/red arrows over the green arrows. I set out in this post to try to bring clarity to how an adjacent oxygen influences how we think of the properties of the C=O functional group, but as happens quite often, the answer you get depends on the measurement you make.
‡ The solution values in e.g. acetonitrile are reduced by ~20 cm-1, reaching the values often quoted in text books for these functional groups. † The effect on C-H values is greater, e.g. a reduction from 3186 to 2967 cm-1.
This post has been cross-posted in PDF format at Authorea.
In the mid to late 1990s as the Web developed, it was becoming more obvious…
I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…
The recent release of the DataCite Data Citation corpus, which has the stated aim of…
Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…
In the late 1980s, as I recollected here the equipment needed for real time molecular…
On 24th January 1984, the Macintosh computer was released, as all the media are informing…
View Comments
"One can calculate cubic and quartic force constants to get an estimate of the effect of anharmonicity on the (harmonic/quadratic) values, which emerges as 23-28 cm-1"
Do you happen to have a good reference to explain how this is done?
The Gaussian manual states the keyword Anharmonic will "Do numerical differentiation along modes to compute zero-point energies, anharmonic frequencies... This option is only available for methods with analytic second derivatives". One can infer from the latter that this numerical procedure is very slow, and only really possible for small molecules.
Exactly, the computational procedure becomes complicated with larger molecules. The same problem that was experienced while trying to construct Zlatan determinant of larger molecules; any molecule Larger than Hydrogen molecule becomes " Many-Body Problem" - One of the few limitation of Quantum Mechanics.
Hello,
Thanks for the wonderful post. I wanted to bring to your attention the following two papers if you had not already seen them:
Filgueras, C. A. L.; Huheey, J. E. J. Org. Chem. 1976, 41, 49
Drago, R. S.; Vogel, G. C.; Needham, T.E. J. Am. Chem. Soc. 1971, 93, 6014
These two discuss the basicity of carbonyl oxygen. It has been shown that the oxygen in ketone is slightly more basic than ester. This is inconsistent with a picture of resonance donation from the alkyl oxygen in esters.
I stumbled across your post when I was looking for explanations on why an ester is less electrophilic than ketones. I find the standard explanation of resonance donation from oxygen less than convincing precisely because of the IR-stretching frequency issue. I believe the reduced electrophilicity of esters is due to the repulsion experienced by the incoming nucleophile from the increased electron density in the C=O bond. Probably, repulsion by the lone pair on the alkyl oxygen adds to this. While reading your post on esters, I came across your posts on Sharpless Epoxidation and Burgi-Dunitz trajectory. These posts are now reading assignments in my course.