Posts Tagged ‘crystalline calixarene network’

Solid carbon dioxide: hexacoordinate carbon?

Friday, September 17th, 2010

Carbon dioxide is much in the news, not least because its atmospheric concentration is on the increase. How to sequester it and save the planet is a hot topic. Here I ponder its solid state structure, as a hint to its possible reactivity, and hence perhaps for clues as to how it might be captured. The structure was determined (DOI 10.1103/PhysRevB.65.104103) as shown below.

The structure of solid carbon dioxide. Click for 3D

The two nominal double bond distances are 1.33Å, whilst a further four O…C contacts in the shape of a square complete the coordination (2.38Å each). All would probably agree that the central carbon is best described as hexa-coordinated. This is also a hot topic. For example, note the claim made recently to have created a hexa-coordinated carbon species by design (Synthesis and Structure of a Hexacoordinate Carbon Compound, DOI: 10.1021/ja710423d) based on a motif derived from an allene:

Designed hexacoordinate carbon. Click for 3D

This claim was supported by an unusual measured property, the electron density ρ(r) and its Laplacian in the putative O…C region. These two properties are one of those (relatively rare) meetings between experiment and quantum mechanics, and their usefulness has been noted in this blog on previous occasions. However, note that in this designed structure, the O…C distances are merely 2.65-2.7Å, significantly longer than in solid carbon dioxide! So carbon dioxide, in a form many of us are familiar with (solid), can certainly be justified as being described as having a hexacoordinate carbon (although we might draw the line at describing it as having hexavalent carbon).

If oxygen atoms can approach the carbon in CO2 to within ~2.4Å, an interesting question can be posed. How close can another carbon get to CO2 without actually reacting and forming a new molecule? C-C bonds, even weak ones, are so much more interesting than C-O bonds! It would have to be a particularly nucleophilic carbon, of course. A search of the August 2010 version of the Cambridge structural database (CSD) reveals no really close approaches of another carbon to CO2. Only about 8 weak examples are found, and here the C-C distances are ~3.0-3.2Å, with the O=C=O angle in the CO2 never less than 170°. In this context, there is an intriguing and very recent report (which has not yet made it into the searchable CSD) of the structure of CO2 trapped in a cavity next to what was claimed to be a molecule of 1,3-dimethyl cyclobutadiene, or CBD (see 10.1126/science.1188002 and the discussion of this article in my earlier blog post). The focus in that report was on the “Mona Lisa of organic chemistry”, namely the CBD unit. One feels that the structure of the adjacent CO2 was of lesser interest to the authors. According to a visual image of this system, the CBD and CO2 pair show quite an intimate approach via their carbon atoms (a ghostly C-C bond is clearly represented). This raises the interesting question of whether the description of this pair should be of two intimate but nevertheless separate and relatively unperturbed molecules not connected by a covalent bond (“more indicative of a strong van der Waals contact than of covalent bonding“) or of a pair fully bound by a covalent C-C bond between them?

The issue of what is an interaction, and what is a bond continues to raise its often controversial head. And quantum theory continues to provide a multitude of interpretations as well.

Reactions in supramolecular cavities – trapping a cyclobutadiene: ! or ?

Sunday, August 8th, 2010

Cavities promote reactions, and they can also trap the products of reactions. Such (supramolecular) chemistry is used to provide models for how enzymes work, but it also allows un-natural reactions to be undertaken. A famous example is the preparation of P4 (see blog post here), an otherwise highly reactive species which, when trapped in the cavity is now sufficiently protected from the ravages of oxygen for its X-ray structure to be determined. A colleague recently alerted me to a just-published article by Legrand, van der Lee and Barboiu (DOI: 10.1126/science.1188002) who report the use of cavities to trap and stabilize the notoriously (self)reactive 1,3-dimethylcyclobutadiene (3/4 in the scheme below). Again sequestration by the host allowed an x-ray determination of  the captured species!

Scheme for production of 1,3-dimethylcyclobutadiene 3 and CO2.

The colleague tells me he has himself already penned an article on the topic and submitted this to a conventional journal. However, their rules decree that whilst it is being refereed, I could not discuss the article here, or indeed even name its author. Assuming his article is published, I will indeed reveal his identity, so that he gets the credit he deserves! Meanwhile, I will concentrate in this blog purely on two other aspects of this reaction which caught my own eye after he brought the article to my attention.

The reaction involves imobilising a precursor 1 in a crystalline calixarene network as shown above, and then in situ photolysis to form the Dewar lactone 2. Further photolysis then results in extrusion of carbon dioxide and the formation of 1,3-dimethyl cyclobutadiene 3 and CO2, both still trapped in the host crystals. Thus imobilised, here they both apparently remain (at 175K) for long enough for their X-ray structure to be determined. What attracted me to this chemistry was the potential of this reaction as a nice example of a Diels Alder reaction occuring in a cavity. The first example of such catalysis was reported by Endo, Koike, Sawaki, Hayashida, Masuda, and Aoyama (DOI: 10.1021/ja964198s) and I have used this in my lectures for many years. This latter example however illustrates the promotion of a cycloaddition, which inside a cavity is accelerated by a factor of ~105, rather than of the reverse cycloelimination. I explain this to students by invoking entropy. Normally, when two molecules react together, there is an entropic penalty, which can add 8 or more kcal/mol to the free energy of activation of a bimolecular reaction in the absence of the cavity.

Structure of entrapped 1,3-dimethylcyclobutadiene, obtained from the CIF file provided via DOI: 10.1126/science.1188002

By a strange coincidence, my name is also on a recently published article (DOI: 10.1021/ol9024259) with other colleagues on the use of (Lewis) acid catalysts to accelerate a type of reaction known as the Prins. This involves the addition of an alkene to a carbonyl group. Now as it happens, the reaction in the scheme above showing 42 happens to combine these features; it is both a Diels-Alder cycloaddition and also involves an alkene adding to a carbonyl compound! It is therefore noteworthy that the claimed reaction 123 + CO2 is done in the presence of a strong acid catalyst, the guanidinium cation 5, which is itself part of the structure of the calixarene-based host. It is represented as X in the scheme above, and can also be identified in the above 3D model via the light blue atoms.

There are however crucial differences between these two earlier examples I quoted and the reaction of 23; the latter is in fact a cycloelimination and not a (cyclo)addition. In other words, according to literature precedent, the guanidinium cation-based cavity should act to accelerate the reverse cycloaddition 42 rather than the forward cycloelimination. Since the isomerisation 34 is thought to be fast, the question arises: how rapid is the reverse reaction 42? In particular, is it slow enough to allow X-ray diffraction data to be collected for 3/4 over the necessary period of 24 hours or more? Legrand, van der Lee and Barboiu do not address this point in their article. Nor is there discussion there of how the cavity and the acid catalyst might influence the position of the equilbrium 23 + CO2.

This is where calculations can help. At the B3LYP/6-311G(d,p) level four different models were selected.

  1. Model A is a simple gas phase calculation for the singlet state, which reveals the free energy barrier for 42 is already quite modest for a Diels-Alder reaction (more typical values are ~22 kcal/mol), due no doubt to the instability/reactivity of the cyclobutadiene. However, at 175K, that would still be quite sufficient to prevent the reverse reaction from occurring to any extent over the period of X-ray data collection.
  2. Model B adds a condensed phase (water) to the model. This serves in part to simulate the condensed crystal environment (which is pretty ionic being a tetra ion-pair). The barrier drops to 12.1 kcal/mol.
  3. Adding one guanidinium cation to both these models (C and D) which simulate the Prins conditions, drops the barrier to 8.3 kcal/mol (model 4).
  4. You can inspect details of any of the calculations by clicking on the digital repository entry (shown as dr in the table), where full data is available.

None of these models includes the entropic effects of full constraint in a cavity (which I estimated above as capable of reducing the free energy barrier for reaction by ~8 or more kcal/mol). If this correction is applied to model D, it would reduce the barrier to ~0 kcal/mol! The calculations also reveal that the reverse reaction 42 is exothermic, and this exothermicity is enhanced by the acid catalyst 5. It would be further enhanced by reducing the entropy of reaction by pre-organizing the reactants in the cavity. The tendency must therefore be for 3/4 to revert to 2 on purely thermodyamic grounds, and only the presence of a significant kinetic barrier would allow them to exist as separate species. This barrier, as one might infer from the calculations shown in the table below, may not be a large one. Even a barrier of 8 kcal/mol might require cooling to significantly lower than 175K to render the reaction slow on a ~24 hour timescale.

Model ΔG4 → 2
kcal/mol
ΔGreac 4 → 2 Singlet-triplet
separation
A. Gas phase,X=none dr ts 16.8 dr -3.5dr +5.7 dr
B. Continuum solvent (water),X=none dr ts 12.1dr -6.0 dr +7.7 dr
C. Gas phase,X=guanidinium+ dr ts 6.1 dr -19.5dr +2.1dr
D. Continuum solvent (water),X=guanidinium+ dr ts 8.3 dr -10.1 dr +7.7 dr

So I end my own speculations here on the nature of the reaction reported by Legrand, van der Lee and Barboiu by asking: are the products they claim (1,3-dimethylcyclobutadiene and carbon dioxide) capable of existing as separate species for long enough inside their cavity, even at 175K, to allow for the collection of X-ray data for a structure determination?

I tend to think probably not (? rather than !). But do decide for yourselves.


Archived as http://www.webcitation.org/5rpkn2Z5S on 08/08/2010. See also this post.