Electrides (aka solvated electrons).

Peter Edwards has just given the 2015 Hofmann lecture here at Imperial on the topic of solvated electrons. An organic chemist knows this species as “e” and it occurs in ionic compounds known as electrides; chloride = the negative anion of a chlorine atom, hence electride = the negative anion of an electron. It struck me how very odd these molecules are and so I thought I might share here some properties I computed after the lecture for a specific electride known as GAVKIS.[1] If you really want to learn (almost) everything about these strange species, go read the wonderful review by Zurek, Edwards and Hoffmann,[2] including a lesson in the history of chemistry stretching back almost 200 years.

GAVKIS consists of a tricyclic aza-ether ligand or cryptand wrapping a potassium atom in the centre, the overall unit having no charge. The oxygen and nitrogen heteroatoms coordinate to the metal, in the process evicting its single electron. The question that struck me is “where does that electron go?”. You see in all normal molecules that electrons are associated with either one, two (or rarely) three nuclei, to form one-centred monosynaptic basins (lone pairs), two-centre or disynaptic basins (i.e. bonds) and more rarely three-centre bonds. The shared-electron two-centre manifestation was of course famously introduced by Gilbert N. Lewis in 1916 (note the centenary coming up!). Knowing where the electron (pairs) are has enabled the technique popular with organic chemists known as arrow pushing, or the VSEPR analysis of inorganic compounds. But an electride has no nucleus associated with it! So how can one describe its location?

Click for 3D structure of GAVFIS

The crystal structure of GAVFIS shows the potassium to be 8-coordinate. Remember, x-rays are diffracted not by a nucleus but by electrons in the molecule. The highest densities are of course associated with electrons in inner shells centered on nuclei and the much lower densities found in conventional bonds are not normally located by this technique (but see here). So it is no surprise to find that this x-ray analysis[1] did not succeed in answering the question posed above; where is the single electron liberated from the potassium atom? They did look for it, but surmised only that would be found in the “noise level electron density in the spaces between them (molecules)“. For GAVFIS, that empty space is actually dumb-bell shaped, and so perhaps an answer is that the electron occupies the dumb-bell shaped spaces between the ligand-potassium complex.

X-ray analysis was defeated by noise; it is an experimental technique after all. But the noise in a quantum mechanical calculation is much smaller; can this reveal where the evicted electron is? Here is the spin density (unpaired electron) distribution for one molecule of GAVFIS computed using the UωB97XD/6-31++(G) DFT method.[3].[4] It is a stratocumulus-like cloud that enshrouds the molecule (click on the diagram below and you can rotate the function to view it from your own point of interest) but interestingly avoiding the regions along the N….N axis. There are also tiny amounts of (negative) spin density on the ligand atoms. So even when the “empty space” is infinitely large, the shape of the electride anion is nevertheless quite specific, but a holistic function of the shape of the entire molecule rather than its component atoms.

Click for 3D

Another way of describing where electrons are is using functions known as molecular orbitals. Below is the SOMO (singly occupied MO) and its shape in this case coincides with that of the spin density.

Click for 3D interaction

The molecular electrostatic potential is rather wackier (red = attractive to protons).

Click for 3D interaction

Odder still is the ELF (electron localisation function) and the identification of the centroids of its basins. These centroids normally coincide with the two-centre basins (bonds) and one-centre basins (lone pairs, inner shell electrons) in normal molecules, both being close to nuclear centres (atoms). For GAVFIS, two unexpected one-centre basins are found close to the two nitrogen atoms in the molecule, each with a population of 0.48 electrons, along with regular one-centre “lone pair” basins pointing inwards to the potassium (2.38 electrons each). The odd-looking pair of locations identified for the electride anion may have little physical reality, except for reminding us that the electride can indeed be in more than one location simultaneously!

I often also use the NCI (non-covalent-interaction) property of the electron density in these blogs. It tells us about regions of non-covalent electron density which represent attractive weak interactions between or within molecules. Here, it again shows us the weak non-covalent density (as the reduced density gradients) wrapping the molecule (green=weakly stabilizing).

Click for 3D interaction

The obvious next question is that if each molecule is surrounded by weak spin density arising from an unpaired electron, would two such species form a dimer in which the spins are paired in an manner analogous to the conventional single bond? The overlap is not going to be fantastic if the spin distribution has the shape shown above, but what the hell. Here is the HOMO of such a species.[5] It appears the shape of the electride is very pliable indeed; they have been squeezed out of the contact region between the two molecules (which form a close contact pair) into wrapping the dimer rather than the monomer! The spin-coupled singlet by the way is about 4.6 kcal/mol more stable in free energy ΔG298 than two isolated monomer doublets, and 5.5 kcal/mol lower than the triplet species[6] which retains two unpaired electrons. A sort of weak molecule-pair bond rather than an atom-pair bond.

Click for 3D interaction

This has hardly started to scratch the surface of the strange properties of electrides. If your appetite has been whetted, go read the article I noted at the beginning.[2]


For normal molecules, a Mulliken or other population analysis reduces the charge and spin density down to an atom-centered distribution. If this is done for GAVFIS, the spin density collapses down to the molecular centroid, in this case the potassium (spin density 1.15). This of course is horribly misleading, and serves to remind us that such atom-centered distributions can sometimes be far from realistic.


References

  1. D.L. Ward, R.H. Huang, and J.L. Dye, "Structures of alkalides and electrides. I. Structure of potassium cryptand[2.2.2] electride", Acta Crystallographica Section C Crystal Structure Communications, vol. 44, pp. 1374-1376, 1988. http://dx.doi.org/10.1107/S0108270188002847
  2. E. Zurek, P. Edwards, and R. Hoffmann, "A Molecular Perspective on Lithium–Ammonia Solutions", Angewandte Chemie International Edition, vol. 48, pp. 8198-8232, 2009. http://dx.doi.org/10.1002/anie.200900373
  3. Henry S Rzepa., "C 18 H 36 K 1 N 2 O 6", 2015. http://dx.doi.org/10.14469/ch/191347
  4. Henry S Rzepa., "C 18 H 36 K 1 N 2 O 6", 2015. http://dx.doi.org/10.14469/ch/191354
  5. Henry S Rzepa., "C 36 H 72 K 2 N 4 O 12", 2015. http://dx.doi.org/10.14469/ch/191348
  6. Henry S Rzepa., "C 36 H 72 K 2 N 4 O 12", 2015. http://dx.doi.org/10.14469/ch/191350
Henry Rzepa

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

View Comments

  • There have been a couple of tweets asking why I am still using Java to display surfaces.

    1. The alternative, JSmol is still rather slow at displaying JVXL surfaces. It is possible if one uses low resolution surfaces, but medium and high res is not really viable.

    2. I did try some time ago to get the JSmol Wordpress extension working on this instance of the blog. Despite revisiting several times, it still does not work. I suspect only a complete reinstall of the Wordpress software and rebuilding from scratch would solve this issue.

    3. The syntax used for WP-JSMol is quite different from that used to invoke Java directly on this blog. So I would have to edit the syntax for >200 posts if I were to replace Java by JSmol.

    These issues combine to make the migration quite difficult. I have just updated the official Java Jmol applet to the latest release, using the official signed version. It seems to be working.

  • As an experiment, I have cloned the above post as http://www.rzepa.net/blog/?p=14272

    This version uses JSmol to view the interactive objects, and unlike the post here, it works (at least for me). I have been baffled by why it works on the other version but not here for more than a year, and have not yet found the cause. It appears to be some sort of conflict with the other extensions I have installed here, but which one eludes me. Go try the other version, but it might be slow to load.

  • A nice example of electrides. I used your wavefunction for QTAIM analysis (for the moment just topological analysis) to see where is the electron. It has been suggested that non-nuclear attractors of AIM theory coincide with free electrons. For instance please see:
    http://pubs.rsc.org/En/content/articlelanding/2015/cc/c5cc00215j#!divAbstract (doi: 10.1039/C5CC00215J)
    or an older paper:
    http://www.sciencedirect.com/science/article/pii/S0166128009000931 (doi: 10.1016/j.theochem.2009.02.003)
    So, I was curious to see which picture is consistent with AIM picture. You may see the molecular graph here:
    https://www.academia.edu/13810211/GAVFIS
    AIM picture is in line with the picture obtained from MESP, showing three basins for the free electron.

    It is quite different from ELF. I can provide some more information about delocalization of the electronic basins to see which atom is interacting more with them or how localized they are in their positions.

  • Very nice Cina! Yes, do provide more information. Also, I have edited in your QTAIM analysis into your comment above. And I would add that your topological analysis also matches the NCI function.

    For completeness, the doi for the wfn (wavefunction) file I used for the ELF, and which can also be used for QTAIM can be found here: 10.14469/ch/191354

    And finally a comment on non-nuclear attractors. My favourite "NNA" is in the H3+ cation, with an NNA at its centre. So its not restricted to free electrons! F2 is another!

  • You are absolutely right in case of H3+ and the fact that NNAs are not restricted to free electrons. There are more examples of molecules having NNA but no free electron. I had looked at three of them in this old paper in early years of my PhD.
    http://www.sciencedirect.com/science/article/pii/S0166128009000931
    Interestingly, I noticed that NNAs behave differently! For instance NNAs in the center of H3+ and B3- are delocalized and interact with atoms. I would say they are floppy. It seems to me they are a manifestation of highly delocalized but shared electrons, which promote covalency in the systems that they are associated with.
    On the other hand, the NNA in the center of Li3+ is essentially different. Electrons are localized in the center. It reminds me of a definition of metalic bonding I learned in high school. A sea of free electrons loosely bonded to positive atoms. The NNA in the Li3+ interacts with three Li+ atoms as if it is an ionic system with free electrons as anions. I have not seen a detailed analysis of NNA-bonding as in my paper!

  • I have obtained the QTAIM analysis of the singlet weak dimer of GAVKIS, but no NNAs are detected (the Poincare-Hopf theorem is satisfied).

  • I checked it with AIM2000 last night. Now, I am analyzing the wavefunction by AIMAll of Todd Keith. I am curious to see if it can find NNAs or not. AIM2000 is a very powerful tool for topological analysis. Would you please let me know which software you employed for AIM analysis.

  • I used AIMALL (latest version). It takes a few hours to hunt everything down to the satisfaction of Poincare-Hopf, it tends to be more thorough than AIM2000.

    The archive for the singlet dimer can be found at 10.14469/ch/191360.

  • You are right. AIMAll does not recognize NNAs. My computations finished with no NNA few minutes ago.
    However, if you completely remove NNAs still Poincare-Hopf satisfies since the set of critical points that are associated with the NNA are all very flat. So, if AIMAll cannot find NNA, it certainly cannot find other CPs connecting the NNA to the other attractors. I will add NNAs manually into the wavefunction and repeat the analysis by AIMAll to see what I can get.
    AIMAll is fast for integration of basin properties but unfortunately is not as precise as AIM2000 for topological analysis. Besides, it is too much user-friendly! It is almost impossible to manipulate it. This case seems to be more interesting than it appeared at first glance.

1 2 3 9

Recent Posts

Detecting anomeric effects in tetrahedral boron bearing four oxygen substituents.

In an earlier post, I discussed a phenomenon known as the "anomeric effect" exhibited by…

2 days ago

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