The importance of being complete.

To (mis)quote Oscar Wilde again, ““To lose one methyl group may be regarded as a misfortune; to lose both looks like carelessness.” Here, I refer to the (past) tendency of molecular modellers to simplify molecular structures. Thus in 1977, quantum molecular modelling, even at the semi-empirical level, was beset by lost groups. One of my early efforts (DOI: 10.1021/ja00465a005) was selected for study because it had nothing left to lose; the mass spectrometric fragmentation of the radical cations of methane and ethane. Methyl, phenyl and other “large” groups were routinely replaced by hydrogen in order to enable the study. Cations indeed were always of interest to modellers; the relative lack of electrons almost always meant unusual or interesting structures and reactions (including this controversial species, DOI: 10.1021/ja00444a012). Inured to such functional loss, we modellers forgot that (unless in a mass spectrometer), cations have to have a counter anion. Here I explore one example of the model being complete(d).

The ion-pair complex of cyclobutadiene.

In the earlier post on this topic I had explored the possibility of a new isomer of cyclobutadiene, induced by the presence nearby of a strong acid, in the form of guanidinium cation. You might note there was no mention of any counterion! Well, here I add it in to complete the model, using perchlorate. I was following in a sense my own advice on Steve Bachrach’s blog, where the NMR spectrum of the adamantly cation was discussed. I had argued there that the anion (I chose SnCl5) might actually have an effect on the NMR. For the cyclobutadiene complex above without a counter-ion, this non-planar form of the cyclobutadiene was calculated earlier to be ~8.5 kcal/mol in free energy higher than the rectangular conventional geometry. Add the perchlorate as above, and this energy difference drops to 4.1 kcal/mol (modelled in water as a solvent). So the counter-ion CAN make a difference!

What are the implications to a modeller of adding counter ions? Well, when you start doing such calculations, you find that the practical matter of optimising the geometry is not quite as straightforward as it is found to be for what I would call covalently bonded systems. These latter have pretty predictable geometries, and these geometries are pretty rigid. Ion-pairs on the other hand are less predictable. Note for example in the above diagram that the perchlorate counterion sits to one side of the molecule, and is not symmetrical. The potential energy surface can be very flat indeed, which means that locating the optimal geometry can be quite a struggle. And unlike a covalent structure, where once the location of the covalent bonds is decided, there is little further ambiguity, ion-pairs may have many different possible relative orientations. Thus the above one may not be unique!

But the last word to this post should be: do not forget counter ions if you a looking at ionic species, and always strive to be complete!

Henry Rzepa

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

View Comments

  • So what is your recipe to solve problems with unsymmetrical multiatomic ion pairs? One has to try as many initial configurations as he can make up and also include various orientations of ions relative to each other? What about the initial separation between ions in the starting configuration?

    P.C. I'm talking about the situations where it is unaffordable to run a complete MD/MC simulation.

  • I do not have a comprehensive recipe, and its clearly a need for someone to develop methodologies in this area.

    What I did have to do to optimise the ion-pair in the (relatively flat) potential is to calculate 2nd derivatives at EACH step of the optimisation! Wow!

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