Posts Tagged ‘Electron’

The H4 (2+) dication and its bonding.

Wednesday, February 15th, 2017

This post arose from a comment attached to the post on Na2He and relating to peculiar and rare topological features of the electron density in molecules called non-nuclear attractors. This set me thinking about other molecules that might exhibit this and one of these is shown below.

The topology of the electron density is described by just four basic types, designed formally by the notation [3,-3], [3,-1], [3,1] and [3,3] and more colloquially by the terms nuclear attractor (NNA), line (or bond) critical point, a ring critical point and a cage critical point respectively. Mostly, the nuclear critical points coincide exactly with the actual nuclear positions, but more rarely these points are not found centered at a nucleus. It was such an NNA that was suggested as a comment on the post on Na2He. There I replied that another example of an NNA is to be found in H3+ and so its a short step to take a look at H42+ in a tetrahedral arrangement (DOI: 10.14469/hpc/2165). Since only two electrons are available for bonding, it is tempting to represent it as below, with dashed partial bonds connnecting the six edges of the tetrahedron and is associated with real normal vibrational modes; ν 416, 1490 and 1861 cm-1. A brief search on Scifinder, which appears to reference this species as hydrogen, ion (H42+), does not identify any publications associated with it (there are studies on H41+ however); if any reader here knows of any discussion please alert us!

Analysing the density however gives a different result. A NNA is located at the centre of the tetrahedron and a line (bond) critical point connects this to each of the four hydrogen nuclei. This result is similar to the obtained for H3+. It is rather odd that these non-nuclear attractors have not entered into the vocabulary used to describe the bonding in simple molecules, but this picture is certainly different from the more empirical dashed lines between the four nuclei that one is instinctively drawn to use (above).

The ELF analysis (performed using multiWFN) is more interesting. The nuclear basins associated with the hydrogens reveal each has 0.425e, but the central one (green below) has its own basin with 0.301e. 

The NICS value associated with the non-nuclear attractor is -27 ppm, which is indicative of strong spherical aromaticity.

All of which goes to show that even the simplest of molecular species may still have properties that are unexpected or certainly not well-known!

Tags:, , , , , , , , ,
Posted in Interesting chemistry | 11 Comments »

Na2He: a stable compound of helium and sodium at high pressure.

Saturday, February 11th, 2017

On February 6th I was alerted to this intriguing article[1] by a phone call, made 55 minutes before the article embargo was due to be released. Gizmodo wanted to know if I could provide an (almost)† instant‡ quote. After a few days, this report of a stable compound of helium and sodium still seems impressive to me and I now impart a few more thoughts here.

The discovery originates from 17 authors based in 17 different institutions, an impressive illustration of global science and cooperation. I illustrate with this diagram, to be found not in the main article body but in its supporting information and for which the caption reads:


Computed charge density (eÅ-3) of Na2He at 300 GPa, plotted in the [110] plane of the conventional cell. The color bar gives the scale.

The nuclei carry of course the greatest charge density, but the density labelled “2e” is not nuclear-centered. This is typical of species known as electrides, where positive cations are associated with just electrons acting as the counter-anion and about which there was an extensive debate earlier on this blog. There is much discussion in the article[1] about the essential role of the He atoms in bringing about the formation of such an electride, an effect that is summarised in a second diagram also found in the supporting information:

I found myself thinking that it would be great to have the first diagram represented as a movie, evolving as the pressure is increased from say ambient to 300 GPa, and presumably showing the “2e” feature (which means diamagnetic electrons) forming as the pressure increases. Would their evolution be abrupt (a step change) or gradual as the pressure increases and the interatomic distances all decrease? As I understand it, this chemical phenomenon is due not so much to the usual coulombic attraction between positive nuclei and negative charge density from the electronic wavefunction leading to e.g. covalent bonds, but to electron repulsions induced by decreasing nuclear separations resulting in electride-like ionisation and hence electron localisation into the “interstitial cavities” of the lattice. Without pressure, you would just have sodium and helium atoms!

The urge to obtain this intriguing electronic wavefunction for myself now appeared (wavefunctions are rarely if ever included in supporting information). To do this you must have atom coordinates available, But such data was not to be found in the supporting information. It was eventually tracked down (by a crystallographer; thanks Andrew!) to the caption in Figure 2.

However, you probably do need to be a crystallographer to convert this data into a set of coordinates. This was done and is here deposited as a CIF file for you to play with if you wish (DOI:10.14469/hpc/2154)[2]. I have reduced the packing of the unit cell obtained from this CIF file (198 atoms) to just 60 and you can enjoy them by clicking on the diagram below. I should point out that if one uses a program that can recognise the periodic lattice such as Crystal (used in the article discussed here), there is no need to make such reductions, but in this instance I wanted to use a program such as Gaussian in discrete (non-periodic) mode, for which the calculation (B3LYP/Def2-SVPD) has DOI: 10.14469/hpc/2156[3] and where you can also find a wavefunction file to play with if you wish.

Click for 3D model

An ELF analysis for this non-periodic wavefunction looks as below. The ELF basins labelled “2e” located in the centre of the cube show an integrated electron population of ~1.9e and correspond to the localised electron pairs noted in the article above.

Click for 3D

The basins on the boundaries of this non-periodic unit show reduced integrations (red arrows below, 0.08 – 1.7e) and are artefacts of the non-periodic approximation introduced.

The ionization into an electride is brought about by the close proximity of the atoms as induced by high pressure. Releasing the pressure would allow the ionized electrons to re-attach themselves to the valence shell of the sodium atoms, thus destroying the unique properties of the system. It is certainly true that this system challenges our normal concepts of what a molecule is. The presence of He is essential and yet its electrons are hardly involved in the re-organised wavefunction. I cannot wait for more examples to be discovered!

†To meet the 55 minute deadline, I was given about 15 minutes thinking time!

‡Instant responses on social media now seem a sine qua non of the political world, so why not the scientific one?!

References

  1. X. Dong, A.R. Oganov, A.F. Goncharov, E. Stavrou, S. Lobanov, G. Saleh, G. Qian, Q. Zhu, C. Gatti, V.L. Deringer, R. Dronskowski, X. Zhou, V.B. Prakapenka, Z. Konôpková, I.A. Popov, A.I. Boldyrev, and H. Wang, "A stable compound of helium and sodium at high pressure", Nature Chemistry, vol. 9, pp. 440-445, 2017. https://doi.org/10.1038/nchem.2716
  2. H. Rzepa, "Na2He: a stable compound of helium and sodium at high pressure.", 2017. https://doi.org/10.14469/hpc/2154
  3. H. Rzepa, "He20Na40", 2017. https://doi.org/10.14469/hpc/2156

Tags:, , , , , , , , , , , , , ,
Posted in Bond slam, crystal_structure_mining, Interesting chemistry | 11 Comments »

Quintuple bonds: resurfaced.

Sunday, January 31st, 2016

Six years ago, I posted on the nature of a then recently reported[1] Cr-Cr quintuple bond. The topic resurfaced as part of the discussion on a more recent post on NSF3, and a sub-topic on the nature of the higher order bonding in C2. The comment made a connection between that discussion and the Cr-Cr bond alluded to above. I responded briefly to that comment, but because I want to include 3D rotatable surfaces, I expand the discussion here and not in the comment.‡

Firstly, a quick update. Since the original post, quite a few Cr-Cr quintuple bonds have been reported. In searching the crystal structure database, I used the text "quintuple" as a text search term (since specifying a quintuple bond as such is not supported) along with a Boolean AND using the sub-structure Cr-Cr (with any type of bond allowed). The result is shown below. It is striking that in fact these "quintuple" bonds cluster into a set with a bond distance of ~1.74Å and another with 1.83Å. Are these valence bond isomers?

 

Now to the system shown at the top (one of the 1.74Å set). My original post discussed the results of a density functional evaluation of the properties of the electron density in the Cr-Cr region. Most striking was the value of the Laplacian ∇2ρ(r) of this density, the value of +1.45au being the largest ever reported for a pair of identical atoms. I should remind that ∇2ρ(r) is used as one measure of the character of a bond, being the balance between electronic kinetic energy density and potential energy density along a bond. But it is well recognised that the bonding between such transition metals has what is called multi-reference character; the wavefunction is not well described by just a single doubly occupied electronic configuration. More electronic configurations have to be included, and hence a MC-SCF (multi-configuration) self-consistent description of the wavefunction is needed. So as a response to the comment noted above, I decided to carry out CASSCF/6-311G(d) calculations, in which an active space of electrons and molecular orbitals is specified, and using the geometry previously obtained at the DFT level. Thus a CASSCF(8,8) calculation takes 8 electrons and evaluates all possible configurations arising from placing them into an active space of eight molecular orbitals. With metals unfortunately the active space is likely to be large, and so I decided to computed (10,10), (12,12) and (14,14) CASSCF as well to see if any convergence might occur. The last is close to the limit offered by the program. The values shown below are at the QTAIM line (bond) critical point along the Cr-Cr axis.

Active space ρ(r) ∇2ρ(r) Total energy, Hartree % of CS config Calculation DOI
8 .303 1.720 -2383.48049 63 [2]
10 .308 1.612 -2383.68830 61 [3]
12 .308 1.612 -2383.70398 60.6 [4]
14 .308 1.612 -2383.72161 59 [5]
DFT .313 1.45 – 100 [6]

From the trend above, we might safely conclude that the CASSCF active space IS convergent, at least for the density if not for the energy. Also convergent are the properties of the density such as ∇2ρ(r), and noteworthy is that the value of this property is even higher than was obtained using single-configuration DFT theory. So the claim that this system has a record such property does not change. Negative values of the Laplacian are normally taken to indicate a conventionally covalent bond, whereas +ve values show the bond has what is called charge-shift character.[7] So these Cr-Cr quintuple bonds must be amongst the most charge-shifted exemplars!

I show some surfaces (click on the image to get a rotatable model) computed from the CASSCF(14,14) density. Firstly the electron density ρ(r) itself, contoured at 0.25au, showing the high value between the chromium atoms.

Next, ∇2ρ(r) contoured at ±1.5, revealing its high value in the Cr-Cr region (blue = +ve, red = -ve) and then below at ± 0.25 which includes the covalent bonds of the ligands.

Finally, the ELF (electron localisation function) function which tries to gather the electron density into localised ELF basins (numbers are the integration of the electron density in this basin). This looks very similar to that shown previously and is striking because there is no basin in the Cr-Cr region. Instead, the localisation is along the Cr-N bonds. One might describe this as saying that the Cr-Cr region is very highly correlated.

‡It is a limitation of the WordPress system that such objects cannot be included in comments.

References

  1. C. Hsu, J. Yu, C. Yen, G. Lee, Y. Wang, and Y. Tsai, "Quintuply‐Bonded Dichromium(I) Complexes Featuring Metal–Metal Bond Lengths of 1.74 Å", Angewandte Chemie International Edition, vol. 47, pp. 9933-9936, 2008. https://doi.org/10.1002/anie.200803859
  2. H.S. Rzepa, "C 2 H 6 Cr 2 N 4", 2016. https://doi.org/10.14469/ch/191860
  3. H.S. Rzepa, "C 2 H 6 Cr 2 N 4", 2016. https://doi.org/10.14469/ch/191857
  4. H.S. Rzepa, "C 2 H 6 Cr 2 N 4", 2016. https://doi.org/10.14469/ch/191858
  5. H.S. Rzepa, "C2H6N2O2", 2016. https://doi.org/10.14469/ch/191855
  6. H.S. Rzepa, "C 2 H 6 Cr 2 N 4", 2010. https://doi.org/10.14469/ch/4156
  7. S. Shaik, D. Danovich, W. Wu, and P.C. Hiberty, "Charge-shift bonding and its manifestations in chemistry", Nature Chemistry, vol. 1, pp. 443-449, 2009. https://doi.org/10.1038/nchem.327

Tags:, , , , , , , , , , ,
Posted in General, Interesting chemistry | 6 Comments »

VSEPR Theory: A closer look at trifluorothionitrile, NSF3.

Saturday, January 16th, 2016

The post on applying VSEPR ("valence shell electron pair repulsion") theory to the geometry of ClF3 has proved perennially popular. So here is a follow-up on another little molecue, F3SN. As the name implies, it is often represented with an S≡N bond. Here I take a look at the conventional analysis.

trifluoorothionitrile

This is as follows:

  1. Six valence electrons on the central S atom.
  2. Three F atoms contribute one electron each.
  3. One electron from the N σ-bond.
  4. Donate two electrons from S to the two π-bonds.
  5. Eight electrons left around central S, ≡ four valence shell electron pairs.
  6. Hence a tetrahedral geometry.
  7. The bond-bond repulsions however are not all equal. The SN bond repels the three SF bonds more than the S-F bonds repel each-other.
  8. Hence the N-S-F angle is greater than the F-S-F angle, a distorted tetrahedron.

Now for a calculation[1];  ωB97XD/Def2-TZVP, where the wavefunction is analysed using ELF (electron localisation function), which is a useful way of locating the centroids of bonds and lone pairs (click on diagram below to see 3D model).

Trifluorosulfonitrile

  • At the outset one notes that there are six ELF disynaptic basins surrounding the central S, integrating to a total of 7.05e. The sulfur is NOT hypervalent; it does not exceed the octet rule.
  • These six "electron sub-pair" basins are arranged octahedrally around the sulfur. The coordination is NOT tetrahedral, as implied above.
  • The three S-N basins have slightly more electrons (1.25e) than the three S-F basins (1.10e), resulting in …
  • the angle subtended at the S for the SN basins being 96° (a bit larger than octahedral) whilst the angle subtended at the S for the SF basins being smaller (89.9°). This matches point 7 above, but is achieved in an entirely different manner.
  • As a result, the N-S-F angle (122.5°) is larger than the ideal tetrahedral angle and the F-S-F angle (93.9°) is smaller, an alternative way of expressing point 7 above.
  • The S≡N triple bond as shown above does have some reality;  it is a "banana bond" with three connectors rather than two. Each banana bond however has only 1.25e, so the bond order of this motif is ~four (not six) but nevertheless resulting in a short S-N distance (1.406Å) with multiple character.

So we have achieved the same result as classical VSEPR, but using partial rather than full electron pairs to do so. We got the same result with ClF3 before. So perhaps this variation could be called "valence shell partial electron pair repulsions" or VSPEPR.

References

  1. H.S. Rzepa, "F 3 N 1 S 1", 2016. https://doi.org/10.14469/ch/191808

Tags:, , , , , , , , , , ,
Posted in Hypervalency | 110 Comments »