Posts Tagged ‘Chemical bond’
Sunday, March 24th, 2019
There is a predilection amongst chemists for collecting records; one common theme is the length of particular bonds, either the shortest or the longest. A particularly baffling type of bond is that between the very electronegative F atom and an acid hydrogen atom such as that in OH. Thus short C-N…HO hydrogen bonds are extremely common, as are C-O…HO.‡ But F atoms in C-F bonds are largely thought to be inert to hydrogen bonding, as indicated by the use of fluorine in many pharmaceuticals as inert isosteres.[1] Here I do an up-to-date search of the CSD crystal structure database, which is now on the verge of accumulating 1 million entries, to see if any strong C-F…HO hydrogen bonding may have been recently discovered.
The search query uses the CF…HO distance as one variable, and the C-F-H angle as the second. The first diagram shows just intermolecular interactions, up to a distance of 2.7Å which is the sum of the van der Waals radii of the two elements. The hot spot occurs at this value, and an angle of ~95°.
The intra-molecular plot shows a similar value for the most common F…H distance, with the interesting variation that the angle subtended at F is about 80°.
The outlier at the short end of the spectrum (arrow) was observed in 2014[2] with the structure shown below. It is indeed the current record holder by some margin! This length by the way is however a great deal longer than the shortest O…HO hydrogen bonds, which can be in the region of 1.2Å (with the proton sometimes symmetrically disposed between the two oxygen atoms). The value is also very similar to the record holder for the shortest C-H…H-C interaction.
It is always useful to check up on crystallographic hydrogen atom positions using a quantum calculation, so here is one at the ωB97XD/Def2-TZVPP level (Data DOI: 10.14469/hpc/5131) which replicates the values nicely.

ωB97XD/Def2-TZVPP Calculation
A QTAIM analysis of the critical points shows that the F…H BCP has a high value of ρ(r) (most hydrogen bonds only reach about 0.03 au).

NBO analysis indicates the E(2) perturbation energy for donation from an F lone pair into the H-O σ* orbital is 21.2 kcal/mol, which indicates a strong H-bond (typical C-O…HO values are 18-22 kcal/mol). The F…H bond order is 0.05.
This molecule has another interesting property, also noted in the original article;[2] the shift in wavenumber of the O-H stretching vibration. Most hydrogen bonds are characterised by the shift (mostly red and recently discovered blue shifts) that occurs in the OH group when it hydrogen bonds. These shifts are typically 100-200 cm-1 but in this molecule there is no shift, which is described as “exceptional”.
The 1H NMR shift of the OH proton is observed at δ 4.8 ppm, with the value calculated here (ωB97XD/Def2-TZVPP) being 4.75 ppm. A very large H-F coupling was observed of 68 Hz, again a very high value for a “through space” hydrogen bond.
So another record for the molecule makers to try to break!
‡Respectively 7142 and 31428 intermolecular (3859 and 10602 intra) examples using the same search parameters as above, with the shortest values being ~1.28 and ~1.2Å.

References
- S. Purser, P.R. Moore, S. Swallow, and V. Gouverneur, "Fluorine in medicinal chemistry", Chem. Soc. Rev., vol. 37, pp. 320-330, 2008. https://doi.org/10.1039/b610213c
- M.D. Struble, C. Kelly, M.A. Siegler, and T. Lectka, "Search for a Strong, Virtually “No‐Shift” Hydrogen Bond: A Cage Molecule with an Exceptional OH⋅⋅⋅F Interaction", Angewandte Chemie International Edition, vol. 53, pp. 8924-8928, 2014. https://doi.org/10.1002/anie.201403599
Tags:Chemical bond, chemical bonding, Chemical elements, Chemistry, Fluorine, Hydrogen, Hydrogen bond, Intermolecular forces, Natural sciences, perturbation energy, pharmaceuticals, Physical sciences, Refrigerants, search parameters, search query, Supramolecular chemistry
Posted in crystal_structure_mining | No Comments »
Sunday, March 4th, 2018
A bond index (BI) approximately measures the totals of the bond orders at any given atom in a molecule. Here I ponder what the maximum values might be for elements with filled valence shells.
Following Lewis in 1916[1] who proposed that the full valence shell for main group elements should be 2 (for the first two elements) and 8 (the “octet“), Bohr (1922[2]), Langmuir (1919-1921[3]) and Bury (1921[4]) extended this rule to include 18 (the transition series) and 32 (the lanthanides and actinides). If we assume no contributions from higher Rydberg shells (thus 3s, 3p, 3d for carbon etc) and an electron pair model for orbital population (which amounts to the single-determinantal model), then the maximum bond index for hydrogen (and helium) would be 1, it would be 4 for main group elements, and then what?
For the special case of hydrogen, I have previously identified (for a hypothetical species) a bond index of 1.33, due mostly to a high Rydberg occupancy of 1.19e. The more normal BI is <1.0, as noted for this hexacoordinated hydride system. My current estimate for the maximum bond index for main group elements is <4.5. Thus for SF6, it has the value of ~4.33 and that includes a modest occupancy of Rydberg shells of 0.36e = 0.18 BI. Exclude these and it is close to 4.
Move on from group 16 to group 6 and you get compounds such as Me4CrCrMe44- or ReMe82- where the metal bond indices are ~6.5.‡ Compounds such as Cr(Me)6 (BI = 5.6) and W(Me)6 (BI = 6.1) are rather lower. This is a long way from 18/2 = 9. The lanthanides and actinides[5] are unlikely to reveal many large BIs (32/2= 16 maximum value) since they are often ionic and the wavefunctions may be too complex to allow a simple index such as a BI to be safely computed.
So if we are hunting for record BIs, the transition elements are the place to hunt. Can a BI of 6.5 be beaten? Can it even approach 9, its maximum value? Does anyone know of candidate molecules?
‡FAIR Data doi: 10.14469/hpc/3352.
References
- G.N. Lewis, "THE ATOM AND THE MOLECULE.", Journal of the American Chemical Society, vol. 38, pp. 762-785, 1916. https://doi.org/10.1021/ja02261a002
- N. Bohr, "Der Bau der Atome und die physikalischen und chemischen Eigenschaften der Elemente", Zeitschrift f�r Physik, vol. 9, pp. 1-67, 1922. https://doi.org/10.1007/bf01326955
- I. Langmuir, "Types of Valence", Science, vol. 54, pp. 59-67, 1921. https://doi.org/10.1126/science.54.1386.59
- C.R. Bury, "LANGMUIR'S THEORY OF THE ARRANGEMENT OF ELECTRONS IN ATOMS AND MOLECULES.", Journal of the American Chemical Society, vol. 43, pp. 1602-1609, 1921. https://doi.org/10.1021/ja01440a023
- P. Pyykkö, C. Clavaguéra, and J. Dognon, "The 32‐Electron Principle", Computational Methods in Lanthanide and Actinide Chemistry, pp. 401-424, 2015. https://doi.org/10.1002/9781118688304.ch15
Tags:Atom, Chemical bond, chemical bonding, chemical properties, Chemistry, metal bond indices, Molecule, Nature, Quantum chemistry, Residential REITs, Resonance, Tennessine, Valence, Valence electron
Posted in Interesting chemistry | No Comments »
Monday, November 27th, 2017
Previously: “Non-polar” species such as SeMe6, SMe6, ClMe3, ClMe5 all revealed interesting properties for the Se-C, S-C or Cl-C “single” bonds. The latter two examples in particular hinted at internal structures for these single bonds, as manifested by two ELF basins for some of the bonds. Here I take a look at the related molecule where a formal double bond between carbon and the central sulfur atom replacing the single-bond might also hint at octet expansions and hypervalence.

Starting with X=Y=Z=CH2,‡ the calculated (ωB97Xd/Def2-TZVPP) geometry has an interesting chiral D3-symmetric form. The density based ELF-basin centroids are shown below, with each formal C=S π-double bond represented by two ELF basins above and below the C-S axis and with each pair of ELF basins being twisted by 48° with respect to the other two pairs. The total valence shell count around the S is 10.98e and the octet is “expanded” (by ~3e).

The orbital-based NBO approach indicates little utilisation of higher (Rydberg) atomic orbital shells (S: [core]3S(1.13)3p(3.35)3d(0.11)4p(0.02); C: [core]2S(1.15)2p(3.77)3p(0.01)3d(0.01) ). Each S-C bond has a Wiberg bond order of 1.36 (significantly less than a double bond), and the central S has an overall bond index of 4.102. There is a real mis-match between the orbital partitioning (2*1.36 = 2.72e) and the ELF partitioning (2*1.83 = 3.66e) into the S-C bonds. The former indicates that ~two of the twelve valence electrons are entering into anti-bonding orbitals to reduce the total bond index from a possible six to just four, but that they still contribute to the electron-density based ELF disynaptic C-S basins. To cast light on this behaviour, successively one to three of the CH2 groups can be replaced by O.

For each “S=O” bond, we find the ELF basin population more or less halves and electrons instead populate the non-bonding O “lone pairs”. The S-C ELF populations in contrast remain approximately constant. These species therefore have “double” S=C bonds but just “single” S-O bonds. The Rydberg population increases slightly; S: [core]3S(1.06)3p(2.95)3d(0.16)4p(0.02)) and the S bond index is 4.18 for one oxygen and S: [core]3S(0.99)3p(2.67)3d(0.19)4p(0.02) and S bond index 4.16 for two oxygens.

Sulfur trioxide (below) seems best represented by S-O rather than S=O bonds. The Rydberg population is S: [core]3S(0.91)3p(2.41)3d(0.21)4p(0.03) and the S bond index is 4.32.

Just for good measure sulfur trisulfide S(S)3 shows rather lower lone pair population because of course it is less electronegative than oxygen, and hence has a slightly greater S-S ELF basin population. Rydberg, S: [core]3S(1.43)3p(3.73)3d(0.21)4p(0.03) and central S bond index 4.04.

It seems molecules where the electrons in a valence shell exceed the “octet” are only too happy to let the excess electrons leak out into adjacent electronegative atoms as lone pairs, where they are no longer classified as “shared”. Trimethylene-λ6-sulfane does not have this option and the excess electrons remain in the region of the valence shell, but here they do not contribute to augmenting the bond index at the central atom. In this specific interpretation, the octet is exceeded, but hypervalence is not induced. It is a slippery concept; one where general agreement about its properties may indeed be difficult to achieve!
‡The FAIR data DOI collection for this post is 10.14469/hpc/3316.
Tags:Chemical bond, chemical bonding, Chemical polarity, Chemistry, double bond, Hypervalent molecule, Nature, single bond, Tetravalence, Valence
Posted in Hypervalency | No Comments »
Sunday, November 12th, 2017
A few years back, I took a look at the valence-shell electron pair repulsion approach to the geometry of chlorine trifluoride, ClF3 using so-called ELF basins to locate centroids for both the covalent F-Cl bond electrons and the chlorine lone-pair electrons. Whereas the original VSEPR theory talks about five “electron pairs” totalling an octet-busting ten electrons surrounding chlorine, the electron density-based ELF approach located only ~6.8e surrounding the central chlorine and no “octet-busting”. The remaining electrons occupied fluorine lone pairs rather than the shared Cl-F regions. Here I take a look at ClMe3, as induced by the analysis of SeMe6.
The difference between ClF3 and ClMe3 is that octet-excess electrons (two in this case) in the former can relocate into fluorine lone pairs by occupying in effect anti-bonding orbitals and hence end up not contributing to the central atom valence shell.‡ With ClMe3 the methyl groups cannot apparently sustain such lone pairs, at least not distinct from the Cl-C bond region. So might we get an octet-busting example with this molecule? A ClMe3 calculation (ωb97xd/6-311++g(d,p)) reveals a molecule with all real vibrational modes (i.e. a minimum, FAIR data DOI: 10.14469/hpc/3241) and ELF (FAIR data DOI 10.14469/hpc/3242)† basins as shown below:

Density-derived approach: Two of the C-Cl bonds each exhibit two ELF basins; one disynaptic basin (0.94e) and one monosynaptic basin (0.20e) closer to the chlorine. The former pair integrate to 1.88e, density which largely arises from carbon (natural charge -0.84) and which contribute to a total integration for these carbons of 7.17e. The latter pair contributes to a total chlorine integration of 7.19e. The angle subtended at chlorine for the two 2.68e “lone pair” basins is 141°. Thus an inner, octet-compliant, valence-shell for chlorine is revealed, plus an expanded-octet outer one into which the two additional electrons go. The latter contribute to forming an octet-compliant carbon valence shell, but may be considered as not contributing to the valence shell of the other atom of the pair, the chlorine. An endo lone-pair rather than the more usual exo lone-pair if you will. These results reveal that the molecular feature we know as a (single) “bond” may in fact have more complex inner structures or zones, something we do not normally consider bonds as having. In this model, these zones are not invariably considered as shared between both the atoms comprising the bond.
Orbital-derived approach: NBO analysis (FAIR data DOI: 10.14469/hpc/3241) reveals the chlorine electronic configuration as [core]3S(1.83)3p(4.67)4S(0.01)3d(0.03)5p(0.02,) showing very little population of the Rydberg shells (4s, 3d, 5p) occurs (0.13e in total). This method of partitioning the electrons allocates a chlorine Wiberg bond index of 2.00 and the methyl carbon bond index of 3.83. If the regular valence of Cl is taken as 1, then the central chlorine can be regarded as non-Rydberg hypervalent (the electrons in the 0.94e basins are taken as contributing to the chlorine bond index).
The carbon-halogen bond internal structures simplify for Br (DOI: 10.14469/hpc/3248, 10.14469/hpc/3250) and I (DOI: 10.14469/hpc/3249, 10.14469/hpc/3247); for each only a single ELF basin is located and the NBO Br and I bond indices are respectively 2.10 and 2.1. This is not due to incursion of Rydberg hypervalence (Br: [core]4S(1.83)4p(4.46)5S(0.02)4d(0.03)6p( 0.01); I: [core]5S(1.82)5p(4.29)6S(0.02)5d(0.02)6p(0.01) ) but of a merging of the carbon and halogen valence basin such that the ELF contributions to each cannot be deconvoluted. In each case the NBO bond indices of ~2 suggest hypervalency for the halogen.


What have we learnt? That the shared electron (covalent) bond can have complex internal features, such as two discrete basins for the apparently shared electrons. How one partitions these electrons can influence the value one obtains for the total shared electron count and hence whether the octet is retained or expanded for main group elements such as the halogens. And finally, that hypervalence and hyper-coordination are related in the orbital model at least. Thus along the series MenI (n= coordination number 1,3,5,7), the values of the Wiberg bond index at the halogen progress as 1.0, 2.1, 3.1 (DOI: 10.14469/hpc/3236) and 4.01 (DOI: 10.14469/hpc/3238), or one extra atom bond index per electron pair. Given this, it seems useful to retain the distinction between the terms hypervalence and hyper-coordination, but also recognize that we still may have much to learn about the former.
‡See the previous post on SeMe6 for a more detailed discussion.
† The FAIR Data accompanying this blog post is organised in a new way here. All the calculations are collected together with an over-arching DOI: 10.14469/hpc/3252 associated with this post, with individual entries accessible directly using the DOIs given above. The post itself has a DOI: 10.14469/hpc/3255 and the two identifiers are associated with each-other via their respective metadata. A set of standards (https://jats.nlm.nih.gov) with implementation guidelines for e.g. repositories, authors and publishers-editors are expected in the future to establish infra-structures for cross-linking narratives/stories with the data on which they are based.
Tags:Chemical bond, chemical bonding, Chemistry, Chlorine, Covalent bond, Lone pair, Oxidizing agents, Quantum chemistry, Stereochemistry, Valence, VSEPR theory
Posted in Chemical IT, Hypervalency | 5 Comments »
Tuesday, October 24th, 2017
An N-B single bond is iso-electronic to a C-C single bond, as per below. So here is a simple question: what form does the distribution of the lengths of these two bonds take, as obtained from crystal structures?
The Conquest search query is very simple (no disorder, no errors).
When applied to the Cambridge structure database (CSD) the following two distributions are obtained. That for carbon is pretty symmetric with the peak at ~1.53Å but with rather faster decay in the region >1.6Å compared with the region <1.46Å (the latter may be caused by hyperconjugation shortening the C-C bond).

In contrast, the iso-electronic N-B distribution is more asymmetric about the peak of 1.56Å, exhibiting a long tail beyond 1.63Å, up to a value of 1.825Å.

The molecule with that longest N-B bond (1.825Å) is shown below; UWOHUK, Data DOI: 10.5517/ccwcwlp. This by the way is no crystal artefact; a calculation (ωB97XD/6-311G(d,p), Data DOI: 10.14469/hpc/3202) gives a calculated length of 1.81Å, with a N-B bond order of 0.48.
Stretching a C-C bond heterolytically requires charge separation (a relatively unfavourable process) and likewise homolytic stretching would tend to form a biradical, in effect an excited state and again not favourable. In contrast, elongating the N-B bond reduces (at least formally) any charge separation and allows this heteronuclear pair to sustain (single) bond lengths over the much wider range of ~0.4Å without requiring biradical formation.
One might wonder what other single-bonded atoms pairs give such unusually large spans in their bond length distributions.
Tags:bond, Bond valence method, Chemical bond, chemical bonding, Chemistry, Covalent bond, crystal structure, Nature, Quantum chemistry, search query
Posted in crystal_structure_mining | 3 Comments »
Saturday, September 16th, 2017
Early in 2011, I wrote about how the diatomic molecule Be2 might be persuaded to improve upon its normal unbound state (bond order ~zero) by a double electronic excitation to a strongly bound species. I yesterday updated this post with further suggestions and one of these inspired this follow-up.
The standard molecular orbital diagram for Be2 below shows two electrons in both the 2s Σg and Σu levels, the first being considered bonding and the second antibonding. By exciting the two electrons from the Σu into the Πu MO to form a triplet, one converts one antibonding occupancy into two bonding occupancies, in the process changing the total formal bond order from zero to two.

The triplet excited state of diberyllium
You can see the results of my playing with these ideas both in my appended comments to the original post and the table below. This shows that the calculated bond order for the excited triplet state of Be2 is actually closer to 1.50 rather than to two, but definitely not zero!
The games above represent isoelectronic substitutions and here I try one more, namely that Li22- is isoelectronic with Be2. Unlike the latter, there is no need to force an electronic excitation (ωB97XD/Def2-QZVPPD/SCRF=water) to achieve the required occupancies with Li22-.
| System |
Wiberg bond order |
Bond length |
FAIR Data |
| Li22- triplet |
1.501 |
2.381 |
10.14469/hpc/3087 |
I also checked what crystal structures could tell us about Li-Li bonds and it seems 2.38Å is about as short as they get. 
At this point, the NBO analysis of the Li22- localised orbitals alerted me to another feature, which is that the Rydberg occupancy amounted to 2.18e. This in turn reminded me of the previous post which dealt with such occupancy in another small molecule, CH3F2-, but here the Rydberg occupancy involved the 3s/3p AOs of the carbon and the fluorine. With Li22- triplet, it is of the lithium 2p AO (2.18e) and only a tiny occupancy of 3d (0.03). By definition, for alkali metals such as Li the normal valence shell is just 2s, whereas 2p occupancy is considered a Rydberg state; a hypervalent state if you will. So Li22- triplet has a Li-Li hyper-bond!‡ Of course, by this definition most Li compounds are then hypervalent, since many have populated 2p shells.
Even if use of the term hyper-bond to describe Li22- triplet is rather artificial, this example does reveal the games one can play with the first row elements Li-B (see table above). Given that most introductory text books on bonding normally only explain the diatomics formed from N-Ne (occasionally including C), I might suggest that these earlier elements are equally instructive and fun to play with.
‡ This species is 36.0 kcal/mol higher in free energy than two separated Li– anions.
Tags:Be-Be double bond, Be-Be triple bond, Chemical bond, Chemistry, Cs-Cs double bond, Diatomic molecule, free energy, General chemistry, K-K double bond, Li-Li double bond, Molecular geometry, Oxygen, Province/State: Be2, Quantum chemistry, Rb-Rb double bond, Stereochemistry
Posted in Interesting chemistry | 3 Comments »
Saturday, September 16th, 2017
Early in 2011, I wrote about how the diatomic molecule Be2 might be persuaded to improve upon its normal unbound state (bond order ~zero) by a double electronic excitation to a strongly bound species. I yesterday updated this post with further suggestions and one of these inspired this follow-up.
The standard molecular orbital diagram for Be2 below shows two electrons in both the 2s Σg and Σu levels, the first being considered bonding and the second antibonding. By exciting the two electrons from the Σu into the Πu MO to form a triplet, one converts one antibonding occupancy into two bonding occupancies, in the process changing the total formal bond order from zero to two.

The triplet excited state of diberyllium
You can see the results of my playing with these ideas both in my appended comments to the original post and the table below. This shows that the calculated bond order for the excited triplet state of Be2 is actually closer to 1.50 rather than to two, but definitely not zero!
The games above represent isoelectronic substitutions and here I try one more, namely that Li22- is isoelectronic with Be2. Unlike the latter, there is no need to force an electronic excitation (ωB97XD/Def2-QZVPPD/SCRF=water) to achieve the required occupancies with Li22-.
| System |
Wiberg bond order |
Bond length |
FAIR Data |
| Li22- triplet |
1.501 |
2.381 |
10.14469/hpc/3087 |
I also checked what crystal structures could tell us about Li-Li bonds and it seems 2.38Å is about as short as they get. 
At this point, the NBO analysis of the Li22- localised orbitals alerted me to another feature, which is that the Rydberg occupancy amounted to 2.18e. This in turn reminded me of the previous post which dealt with such occupancy in another small molecule, CH3F2-, but here the Rydberg occupancy involved the 3s/3p AOs of the carbon and the fluorine. With Li22- triplet, it is of the lithium 2p AO (2.18e) and only a tiny occupancy of 3d (0.03). By definition, for alkali metals such as Li the normal valence shell is just 2s, whereas 2p occupancy is considered a Rydberg state; a hypervalent state if you will. So Li22- triplet has a Li-Li hyper-bond!‡ Of course, by this definition most Li compounds are then hypervalent, since many have populated 2p shells.
Even if use of the term hyper-bond to describe Li22- triplet is rather artificial, this example does reveal the games one can play with the first row elements Li-B (see table above). Given that most introductory text books on bonding normally only explain the diatomics formed from N-Ne (occasionally including C), I might suggest that these earlier elements are equally instructive and fun to play with.
‡ This species is 36.0 kcal/mol higher in free energy than two separated Li– anions.
Tags:Be-Be double bond, Be-Be triple bond, Chemical bond, Chemistry, Cs-Cs double bond, Diatomic molecule, free energy, General chemistry, K-K double bond, Li-Li double bond, Molecular geometry, Oxygen, Province/State: Be2, Quantum chemistry, Rb-Rb double bond, Stereochemistry
Posted in Interesting chemistry | 3 Comments »
Wednesday, September 6th, 2017
The chemical bond zoo is relatively small (the bond being a somewhat fuzzy concept, I am not sure there is an actual count of occupants). So when two new candidates come along, it is worth taking notice. I have previously noted the Chemical Bonds at the 21st Century-2017: CB2017 Aachen conference, where both were discussed.
- The first now has a name, the exo-bond, one example of which is the C2 diatomic. The hint that a quadruple bond could be formulated between the two carbon atoms goes back a little while[1] (see Table 1), but revived interest really took off after ~2010, around the time my blog on the topic also appeared. You can see the abundance of post 2010 articles in the bibliography at the Aachen bond-slam. At the conference, four speakers all agreed using rather different methods that there was indeed “something” additional to a C≡C triple bond and that this “something” might be worth ~15-30 kcal/mol of stabilization.† The debate centered around whether this term deserved to be called a bond, or whether it should be downgraded to merely that of biradicaloid stabilizations.‡ The more conventional population of a σ*-antibond it was argued would not result in such stabilizations. Since many kinds of bonds have stabilization energies of similar magnitude, not least the weaker hydrogen bonds, agostic bonds, halogen bonds etc, let us for the sake of argument call it a bond here. Because four electrons might occupy the same space along the σ-symmetric C-C axis, they experience significant so-called static correlation∞ which results in partition into one electron pair occupying the central region (the endo-bond) and the other pair in the outer region (the exo-bond). This separation decreases the Pauli electron repulsions along the entire C-C axis region. An example of an exo-bond is found in [1.1.1] propellane, where the notional central C-C bond is thought to actually occupy the region outside the central C-C bond axis, but largely in this example because of angular strains. In this case however, the propellane bond is not competing with an endo-bond along the same axis. We might conclude therefore that the convention of characterising a bond using the separation between the two nuclei (the bond-length) is rather stressed when one has two different bonds along the axis of the nuclei, one of which is obviously “longer” than the other.
Which brings us to representations; e.g. Chemdraw now allows drawing of quadruple bonds and so it can be drawn thus quite simply.

The second form breaks the century-old convention that all bonds along a diatomic axis are drawn in the same manner, by isolating the exo-bond to make the point clear. Perhaps we should stick to the first, but be prepared to explain the underlying complexity of the quantum mechanical symmetries as we do to students with σ/π/δ/φ bonds, which are another mechanism for avoiding having bonds in exactly the same regions. I know the story has not yet ended; but is it time to at least speculate when the text-books will start to reflect/discuss the exo-bond?
- The second I dub the hyper-bond. This goes back to G. N. Lewis and his famous octet rule for main group elements, the expansion of which was subsequently described by the term hypervalent. That term has become rather confused with hypercoordinate, since hypervalent is often used to describe hypercoordinate species such as PCl5, SF6 or I.I7. But this does rather break the original definition, since few if indeed any of these hypercoordinate molecules have a significantly expanded octet shell. At the Aachen meeting, a molecule fitting the original definition was presented, appearing first in early form on this blog. Put simply, a wavefunction for CH3F2- can be calculated (ωB97XD/Def2-QZVPPD/SCRF=water, DOI: cb3n ) for which the two additional electrons populate a molecular orbital with significant contributions from the 3s/3p valence shell AOs (atomic orbitals) for both carbon and fluorine. The alternative would have been to populate the anti-bonding C-H or C-F orbitals composed of 2s/2p valence shell AOs. The former results in a total population of these higher valence shells of 1.55e and makes the C-F (Wiberg) bond order >1 (1.14) and the total Wiberg bond indices >4 for carbon (4.162) and >1 for F (1.275). The resulting HOMO (highest occupied molecular orbital) or NBO (they are very similar) looks as below. It takes the approximate form of a torus or cylinder wrapping the inner C-F bond, a second layer to the C-F bond if you wish.
 Normal valence shell F-C σ-orbital defining the regular C-F bond.
|
 Higher valence shell F-C σ-orbital defining the C-F hyper-bond.
|
Rather than the entire molecule being defined as hypervalent, only one (in this case localized) orbital is given the term and the other orbitals are conventional.
In both cases the molecules are either very reactive (C2) or with such a low barrier to fragmentation (into CH3– and F– for CH3F2-) that detection of the latter is unlikely. But these are interesting Gedanken experiments in quantum mechanics, which in turn catalyse the development of new techniques and in some cases might even lead to the design and isolation of new types of molecules.
†The known thermochemistry of the two reactions; HC≡CH → HC≡C + H•; HC≡C• → CC + H• is ~17 kcal/mol less endothermic for the second step, suggesting some factor is needed to account for the additional stabilization when CC is formed.
‡The singlet to triplet excitation energy for C2 is ~+30 kcal/mol, so the biradicaloid electrons are certainly spin-coupled.
∞Other “difficult” correlated molecules include Be2 and B2.
References
- R.S. Mulliken, "Note on Electronic States of Diatomic Carbon, and the Carbon-Carbon Bond", Physical Review, vol. 56, pp. 778-781, 1939. https://doi.org/10.1103/physrev.56.778
Tags:bond, Chemical bond, chemical bond zoo, City: Aachen, Concepts, Dialectic, Fuzzy concept, Non-classical logic, Psychometrics, triplet excitation energy
Posted in Bond slam, Interesting chemistry | No Comments »
Thursday, April 6th, 2017
Enols are simple compounds with an OH group as a substituent on a C=C double bond and with a very distinct conformational preference for the OH group. Here I take a look at this preference as revealed by crystal structures, with the theoretical explanation.

First, a search of the Cambridge structure database (CDS), using the search query shown below (DOI: 10.14469/hpc/2429)


The first search (no errors, no disorder, R < 0.05) is unconstrained in the sense that the HO group is free to hydrogen bond itself. The syn conformer has the torsion of 0° and it has a distinct preponderance over the anti isomer (180°). There is the first hint that the most probable C=C distance for the syn isomer may be longer than that for the anti, but this is not yet entirely convincing.
To try to make it so, a constrained search is now performed, in which only structures where the HO group has no contact (hydrogen bonding) interaction are included. This is achieved using a “Boolean” search;

The number of hits approximately halves, but the proportion of syn examples increases considerably. There is an interesting double “hot-spot” distribution, which amplifies the lengthening of the C=C bond compared to the anti orientation.

The next constraint added is that the data collection must be <100K (to reduce thermal noise) which reduces the hits considerably but now shows the lengthening of the C=C bond for the syn isomer very clearly.
A final plot is of the C=C length vs the C-O length (no temperature, but HO interaction constraint). If there were no correlation, the distribution would be ~circular. In fact it clearly shows that as the C=C bond lengthens, the C-O bond contracts.

Now for some calculations (ωB97XD/Def2-TZVPP, DOI: 10.14469/hpc/2429) which reveal the following:
- The free energy of the syn isomer is 1.2 kcal/mol lower than that of the syn. The effect is small, and hence easily masked by other interactions such as hydrogen bonding to the OH group. Hence the reason why removing such interactions from the search above increased the syn population compared to anti.
- The syn C=C bond length (1.325Å) is longer than the anti (1.322Å).
- The syn isomer has one unique σO-Lp/σ*C-C NBO orbital interaction (below) with a value of E(2) 7.7 kcal/mol, which is absent in the anti form. As it happens, a πO/π*C=C interaction is present in both forms but is also stronger in the syn isomer (E(2)= 46.8 vs 44.2 kcal/mol).
| unoccupied NBO, σ*C-C |
|
| Occupied NBO, σO-Lp |
 |
- The overlap of the filled σO-Lp with the empty σ*C-C orbital is shown below (blue overlaps with purple, red overlaps with orange).

To view the overlap in rotatable 3D, click on any of the colour diagrams above.‡
It is nice to see how experiment (crystal structures) and theory (the calculation of geometries and orbital interactions) can quickly and simply be reconciled. Both these searches and the calculations can be done in just one day of “laboratory time” and I think it would make for an interesting undergraduate chemistry lab experiment.
‡ This visualisation uses Java. Increasingly this browser plugin is becoming more onerous to activate (because of increased security) and some browsers do not support it at all. The macOS Safari browser is one that still does, but you do have to allow it via the security permissions.
Tags:Chemical bond, chemical bonding, Chemistry, Conformational isomerism, constrained search, Enol, free energy, Gauche effect, Hydrogen bond, Isomerism, Java, Physical organic chemistry, search query, Stereochemistry, Supramolecular chemistry
Posted in crystal_structure_mining, reaction mechanism | 2 Comments »
Saturday, March 25th, 2017
The previous post demonstrated the simple iso-electronic progression from six-coordinate carbon to five coordinate nitrogen. Here, a further progression to oxygen is investigated computationally.
The systems are formally constructed from a cyclobutadienyl di-anion and firstly the HO5+ cation, giving a tri-cationic complex. There are no examples of the resulting motif in the Cambridge structure database. A ωB97XD/Def2-TZVPP calculation (DOI: 10.14469/hpc/2350) shows it is again a stable minimum, with a Kekule mode of 1203 cm-1.
A QTAIM topological analysis of the electron density shows it differs from the nitrogen analogue in now having the ring topological feature for the basal four carbons, which in turn gives rise to a cage critical point (blue dot). The values of the electron density are lower than for N.
The ELF basin analysis shows the C-C bonds are regular single ones (2.01e), whereas the C-O bonds have a slightly greater electron population than the C-N bonds discussed in the previous post.

I suspect the prospects of making a stable tri-cation in such a small molecule are lower than the crystal di-cation achieved with carbon as the apical atom. But the charge can be reduced to a di-cation by replacing the HO5+ above with S–-O5+; the animation below showing the Kekule mode (1140 cm-1, DOI: 10.14469/hpc/2356).


And for some (negative) loose ends.
- The P equivalent constructed from cyclobutadienyl di-anion and HP4+ is now unremarkably 5-coordinate. But in fact it is not a stable minimum (DOI: 10.14469/hpc/2357), having two negative force constants.
- as does the system from cyclobutadienyl di-anion and O=P4+(DOI: 10.14469/hpc/2358)
- and the system from cyclobutadienyl di-anion and HS5+(DOI: 10.14469/hpc/2360).
- Transposition of S/O to give O–-S5+ likewise (DOI: 10.14469/hpc/2359).
So the family of hyper-coordinate 2nd row main group elements now comprises the experimentally verified C, with N and O now open to such verification.
Tags:animation, Chemical bond, Chemistry, Matter, Nitrogen, Quantum chemistry
Posted in Bond slam, Hypervalency | 4 Comments »