Archive for the ‘crystal_structure_mining’ Category

Molecules of the year 2019: Hexagonal planar crystal structures.

Thursday, January 23rd, 2020

Here is another selection from the Molecules-of-the-Year shortlist published by C&E News, in which hexagonal planar transition metal coordination is identified. This was a mode of metal coordination first mooted more than 100 years ago,[1] but with the first examples only being discovered recently. The C&E News example comprises a central palladium atom surrounded by three hydride and three magnesium atoms, all seven atoms being in the same plane.

As the original article makes clear,[1] the relative orbital simplicity of these early main group based ligands allows the bonding to be better understood, hence itself allowing “additional design principles” to be introduced for transition metal complexes. Here I thought I might extend the scope of this motif by a generalised crystal search for any other hexagonal planar structures to be found in the Cambridge crystal structure database.

A search query can be constructed by defining a plane using the six ligand atoms and then constraining the perpendicular distance between this plane and the transition metal atom at the centre to < 0.1Å. Six angles between adjacent ligands are then themselves constrained to the range 0-70° and the coordination of the central atom can also be constrained to either 6, 7 or 8. All searches are also defined by no disorder and no errors. The search queries can be found at DOI: 10.14469/hpc/6731

I carried out a number of separate searches. The least constrained (any coordination number at the central atom, and any type of attached ligand atom) produced 62 hits, exhibiting a variety of sometimes complex coordination modes. To simplify the search, I separated the searches into specific types. You can view 3D models of any of the molecules below by clicking on the static image.

  1. The first restricted the transition metal atom to 6-coordinate and for which the ligands all derive from the early main group periods (1A, 2A and H). The sole example (NORLOY[[1], DOI: 10.5517/ccdc.csd.cc2235v7) is the one noted above.
  2. The next search restricted all the ligands to a transition metal connected to the central atom, along with 6-coordination. Again just one hit (RISQEP[2] dating from 1997 and comprising a central Au atom surrounded by four Au-based ligands and two Fe-based ligands. This type of molecule is a member of a class known as a hexagonal planar metal cluster.

  3. This search now constrains all six ligands as comprising late main group atoms bonding to the central metal. Two hits, VOVZOV, dating from 1992[3] with a P-based ligand Ni central atom and ZUDWUQ from 1996[4] using As and Ni.

  4. The next category combines the previous two, with ligands either from the transition series or the late main group series, resulting in four more hits (HACBOF and HACBUL[5] (DOI: 10.5517/ccdc.csd.cc1n26pm and 10.5517/ccdc.csd.cc1n26qn) using a combination of three carbons (as acetylide)  and three Ag ligands, with Cu as the central ligand.

  5. Next, any ligand is allowed, together with a 7-coordinate central atom. Three examples, including  VAPZEU (2016, [6], DOI: 10.5517/ccdc.csd.cc1md045) with three Pd and three Si hexagonal ligands, with an additional 7th Cu surrounding a central Pd.

     

  6. One example (of eight found) with 8-coordination at central atom (1996, NANPOH[7]), a central Cd, six oxygen ligands and two further cyanide axial ligands.

  7. To finish, a rather wacky polymeric example with Ti at the centre, six hydrogens deriving from a terminal borane as the hexacoordinate planar motif and two axial P ligands (PEDJOY[8], DOI: 10.5517/ccb0d6x).

I hope this short journey through hexacoordinate planar transition metal complexes has revealed at least a flavour of the diversity in this category. I am also going to take some gentle issue with the C&E News reporting of this molecule, “Scientists proposed a hexagonal planar geometry more than 100 years ago, but it has never been captured in crystal form until now” (referring to the 2019 article which inspired this blog[1]). As I hope I have shown, a number of the examples above in crystal form actually emerged rather earlier than 2019!

References

  1. M. Garçon, C. Bakewell, G.A. Sackman, A.J.P. White, R.I. Cooper, A.J. Edwards, and M.R. Crimmin, "A hexagonal planar transition-metal complex", Nature, vol. 574, pp. 390-393, 2019. https://doi.org/10.1038/s41586-019-1616-2
  2. V.G. Albano, M.C. Iapalucci, G. Longoni, L. Manzi, and M. Monari, "Synthesis of [Au<sub>3</sub>Fe<sub>2</sub>(CO)<sub>8</sub>(dppm)]<sup>-</sup> and [Au<sub>5</sub>Fe<sub>2</sub>(CO)<sub>8</sub>(dppm)<sub>2</sub>]<sup>+</sup>:  X-ray Structures of [NEt<sub>4</sub>][Au<sub>3</sub>Fe<sub>2</sub>(CO)<sub>8</sub>(dppm)] and [Au<sub>5</sub>Fe<sub>2</sub>(CO)<sub>8</sub>(dppm)<sub>2</sub>][BF<sub>4</sub>]", Organometallics, vol. 16, pp. 497-499, 1997. https://doi.org/10.1021/om960850g
  3. R. Ahlrichs, D. Fenske, H. Oesen, and U. Schneider, "Synthesis and Structure of [Ni(P<i>t</i>Bu<sub>6</sub>)] and [Ni<sub>5</sub>(P<i>t</i>Bu)<sub>6</sub>(CO)<sub>5</sub>] and Calculations on the Electronic Structure of [Ni(P<i>t</i>Bu)<sub>6</sub>] and (PR)<sub>6</sub>, R = <i>t</i>Bu,Me", Angewandte Chemie International Edition in English, vol. 31, pp. 323-326, 1992. https://doi.org/10.1002/anie.199203231
  4. E. Hey‐Hawkins, M. Pink, H. Oesen, and D. Fenske, "Synthesen und Charakterisierung von [Ni(<i>t</i>BuAs)<sub>6</sub>] und [Pd(<i>t</i>BuAs)<sub>6</sub>]", Zeitschrift für anorganische und allgemeine Chemie, vol. 622, pp. 689-691, 1996. https://doi.org/10.1002/zaac.19966220420
  5. S.C.K. Hau, M.C. Yeung, V.W. Yam, and T.C.W. Mak, "Assembly of Heterometallic Silver(I)–Copper(I) Alkyl-1,3-diynyl Clusters via Inner-Core Expansion", Journal of the American Chemical Society, vol. 138, pp. 13732-13739, 2016. https://doi.org/10.1021/jacs.6b08674
  6. M. Tanabe, R. Yumoto, T. Yamada, T. Fukuta, T. Hoshino, K. Osakada, and T. Tanase, "Planar PtPd<sub>3</sub> Complexes Stabilized by Three Bridging Silylene Ligands", Chemistry – A European Journal, vol. 23, pp. 1386-1392, 2016. https://doi.org/10.1002/chem.201604502
  7. J. Kim, and K. Kim, "NEW THREE DIMENSIONAL [Cd(CN)<sub>2</sub>]<sub>n</sub> FRAMEWORK FORMED WITH CADMIUM CYANIDE AND Cd(CN)<sub>2</sub>·(18-CROWN-6): CRYSTAL STRUCTURE OF [Cd(CN)<sub>2</sub>]·1/2[Cd(CN)<sub>2</sub> (18-CROWN-6)]·3/2EtOH<sup>+</sup>", Journal of Coordination Chemistry, vol. 37, pp. 7-15, 1996. https://doi.org/10.1080/00958979608023536
  8. D.M. Goedde, and G.S. Girolami, "Titanium(II) and Titanium(III) Tetrahydroborates. Crystal Structures of [Li(Et<sub>2</sub>O)<sub>2</sub>][Ti<sub>2</sub>(BH<sub>4</sub>)<sub>5</sub>(PMe<sub>2</sub>Ph)<sub>4</sub>], Ti(BH<sub>4</sub>)<sub>3</sub>(PMe<sub>2</sub>Ph)<sub>2</sub>, and Ti(BH<sub>4</sub>)<sub>3</sub>(PEt<sub>3</sub>)<sub>2</sub>", Inorganic Chemistry, vol. 45, pp. 1380-1388, 2006. https://doi.org/10.1021/ic051556w

Can a carbon radical act as a hydrogen bond acceptor?

Saturday, December 28th, 2019

Having shown that carbon as a carbene centre, C: can act as a hydrogen bond acceptor, as seen from a search of crystal structures, I began to wonder if there is any chance that carbon as a radical centre, C could do so as well. Definitely a subversive thought, since radical centres are supposed to abstract hydrogens rather than to hydrogen bond to them.

One molecule that emerges from such a search (Query: 10.14469/hpc/6572) was reported recently as a resonance stabilized radical,[1] with the intermolecular hydrogen bond that emerges being from an aryl C-H directed at the carbon radical centre. The length (after correction by -0.1Å) is typical (this interaction is not noted in the article itself). Most of the 31 hits are in fact intra-molecular.

Click image to view 3D model

The spin density arising from the unpaired electron of the radical is indeed delocalised, although the largest part is in a pπ orbital on the carbon radical centre.

The molecular electrostatic potential (MEP) shows a negative π-potential (attractive to a proton) covering the radical carbon, but also the adjacent boron.

These types of hydrogen bond to a carbon radical acceptor are clearly weak (if indeed they are real), but perhaps a balance has to be achieved between two effects: less delocalised carbon radicals might form stronger hydrogen bonds but they will also abstract hydrogen atoms from potential hydrogen bond donors. More highly delocalised radicals are less likely to abstract, but probably also less likely to form strong hydrogen bond acceptors. Nonetheless, one can ask whether a stronger carbon radical hydrogen bond acceptor might be found that exists in that region where abstraction does not occur. As I noted at the start, I am trying to be provocative!

References

  1. T. Kushida, S. Shirai, N. Ando, T. Okamoto, H. Ishii, H. Matsui, M. Yamagishi, T. Uemura, J. Tsurumi, S. Watanabe, J. Takeya, and S. Yamaguchi, "Boron-Stabilized Planar Neutral π-Radicals with Well-Balanced Ambipolar Charge-Transport Properties", Journal of the American Chemical Society, vol. 139, pp. 14336-14339, 2017. https://doi.org/10.1021/jacs.7b05471

CH…O hydrogen bonding competing with layered dispersion attractions.

Friday, July 19th, 2019

I have previously looked at the topic of hydrogen bonding interactions from the hydrogen of chloroform Here I generalize C-H…O interactions by conducting searches of the CSD (Cambridge structure database) as a function of the carbon hybridisation. I am going to jump straight to a specific molecule XEVJIR (DOI: 10.5517/cc5fgpq) identified from the searches appended to this post as interesting for further inspection.[1]

The distances from the carbonyl oxygen to CH groups of an adjacent intermolecular molecule are shown, revealing a bifurcated strong + weaker CH…O interaction. I would note that the CH…O distances are un-normalized, in the sense that a C-H distance obtained from X-ray diffraction data is normally about 0.1Å too short. A corrected value for the H…O distance is probably closer to 1.994Å. Next, a B3LYP+G3BJ/Def2-TZVPP calculation of just this dimeric interaction, which shows a somewhat different pattern, particularly from the carbonyl to the sp3-C-H (FAIR data DOI: 10.14469/hpc/5943) with one distance being shorter and one longer.

Click to load 3D model

A QTAIM analysis reveals the electron density ρ(r) of 0.021au, a relatively high value indicating a relatively strong interaction.

Side-views reveals a possible reason for why the calculation does not match the crystal structure. In the crystal structure, the sp3-CH2 group adopts a different conformation from that computed for just two interacting molecules, since this shape allows more efficient stacking of layers and hence allowing stabilizing dispersion energy between the layers to overcome some loss of hydrogen bonding energies in the plane of the layer. If this packing constraint is removed in the pure dimer, one sp3-CH moves into the plane allowing a shorter interaction to the carbonyl oxygen and the other sp3-CH adopts a pure axial position, unconstrained by any packed layer above it.The absence of layered dispersion attractions is hence compensated by forming strong CH…O interactions.

A calculation using six molecules arranged in three layers of two is an attempt to add back at least some of the layering dispersion terms (a full periodic boundary lattice calculation is the proper way of doing this calculation, but at the level chosen here would take far too much computer time!). The new CH…O distances are now 2.018 and 2.384Å (compared to 2.036 and 2.199Å for a model with just two molecules). Probably, more layers would be needed to replicate the crystal structure more accurately.


And now for the searches. The first is for sp-hybridised carbon, as an intermolecular interaction (R < 0.05, no errors, no disorder, T=<150K, H-position normalised for distances shorter than the sum of the vdW radii -0.4), for which a clear hot spot occurs at a H…O distance of ~2.1Å

Intermolecular to sp carbon

Next, sp2-C as an intermolecular interaction (T=<90K), where the hot spot is less distinct, being at the distance cut-off specified for the search. The shortest distance is ~2.0Å. I will return to this example shortly.

Intermolecular to sp2 carbon

An intramolecular version of this search shows a clearer hotspot, again at ~2.15Å

Intramolecular to sp2 carbon

Next, intramolecular sp3 hybridisation, for which there few examples with no clear hotspot.

Intramolecular to sp3 carbon

Finally, intermolecular sp3 hybridisation. The H…O distance hotspot is very slightly longer, as might be expected for a less acidic hydrogen. Nonetheless, the variation in the H…O distances with hybridisation is perhaps unexpectedly small.


To summarise, by performing a general search of the crystal structure database, one can identify general trends and then go to inspect outliers. In this case, this brought the focus onto an (dare I say otherwise umremarkable) molecule in which layers of aromatic molecules set up a competition between intra-layer CH…O hydrogen bonding and inter-layer dispersion stabilizations. I suspect this competition between these two type of weak interactions is far more common than is generally recognised.

 

References

  1. K.S. Huang, M.J. Haddadin, M.M. Olmstead, and M.J. Kurth, "Synthesis and Reactions of Some Heterocyclic Azacyanines<sup>1</sup>", The Journal of Organic Chemistry, vol. 66, pp. 1310-1315, 2001. https://doi.org/10.1021/jo001484k

The shortest known CF…HO hydrogen 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

  1. 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
  2. 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

How FAIR are the data associated with the 2017 Molecules-of-the-Year?

Wednesday, March 7th, 2018

C&EN has again run a vote for the 2017 Molecules of the year. Here I take a look not just at these molecules, but at how FAIR (Findable, Accessible, Interoperable and Reusable) the data associated with these molecules actually is.

I went about finding out as follows:

  1. The article DOI for all seven candidates was linked to the C&EN site.
  2. From there I manually tracked down the Supporting information
  3. Some of this SI gave a CCDC deposition number for crystal structure data for the molecule in question. The easiest way of going directly to the data was to use the search.datacite.org search engine and to enter the keywords CCDC + deposition number. This gives a DOI for the data, examples of which are included in the table below.
  4. In other examples, I used the CSD Conquest search program and entered the names of 2-3 of the authors of the articles. This also worked well.
  5. Most of the SI files, downloaded as PDF files also had static images of NMR spectra included. This is not active data, and hence does not fulfil the F and I of FAIR, and probably the A as well. None of it is FAIR as defined by my post here although it is actually really easy to make it so. One of the examples had ~116 spectra so unFAIRed.
  6. In another example there was also computational data, included simply as a set of XYZ coordinates and again contained in the PDF file. This too is not really FAIR, since one has to know how to extract it from this container and repurpose it. It also represents a tiny subset of the data potentially available.
How FAIR are the data associated with the 2017 Molecules-of-the-Year?
# Title Article DOI Data DOI
1 Persulfurated Coronene: A New Generation of “Sunflower” 10.1021/jacs.6b12630 Data available only as PDF
Hosted by Figshare
The SI also has its own DOI:
10.1021/jacs.6b12630.s001
2 A Truncated Molecular Star 10.1021/jacs.6b12630 Crystal structure data:
10.5517/ccdc.csd.cc1nb303
3 Synthesis of trinorbornane 10.1039/c7cc06273g Crystal structure data:
10.5517/ccdc.csd.cc1p7806
4 Braiding a molecular knot with eight crossings 10.1126/science.aal1619 Crystal structure data:
10.5517/ccdc.csd.cc1m85y0
5 Unique physicochemical and catalytic properties dictated by the B3NO2 ring system 10.1038/nchem.2708 Crystal structure data:
10.5517/ccdc.csd.cc1lkff0
6 Total synthesis of mycobacterial arabinogalactan containing 92 monosaccharide units 10.1038/ncomms148510 116 NMR spectra available only as PDF. No crystal structure
7 Nitrogen Lewis Acids 10.1021/jacs.6b12360 NMR spectra available only as PDF.
Computed coordinates available only as PDF
Crystal structures data:
CCDC 1457983-1457987,1458000-1458001
e.g. 10.5517/ccdc.csd.cc1ky4qc
10.5517/ccdc.csd.cc1ky4rd

The FAIRness of the data for these molecules of the year is largely rescued by the crystal structure data deposited with the CCDC in their CSD database and rendered F of FAIR by the persistent identifiers such as the (parochial) deposition numbers or the more general DOI. Now if the NMR and computational data were also covered in this way, we would be making great progress. There are of course many other types of data included with these examples, and procedures for making such data also FAIR have to be worked out by the community.

In order to construct the table above, I had to put about two hours of effort into tracking down the items (and this only because I have done this sort of search before). Perhaps next year I might persuade C&EN to include such a table in their own article!

FAIR data ⇌ Raw data.

Thursday, December 7th, 2017

FAIR data is increasingly accepted as a description of what research data should aspire to; Findable, Accessible, Inter-operable and Re-usable, with Context added by rich metadata (and also that it should be Open). But there are two sides to data, one of which is the raw data emerging from say an instrument or software simulations and the other in which some kind of model is applied to produce semi- or even fully processed/interpreted data. Here I illustrate a new example of how both kinds of data can be made to co-exist.

I will start with a recent publication[1] with the title Crystallographic Snapshot of an Arrested Intermediate in the Biomimetic Activation of CO2The nature of this intermediate caught the eye of another research group, who responded with their own critique[2] along with the comment “However, since we have no access to the original crystallographic data …” They might have been referring to the semi-processed data (containing the so-called hkl structure factors) but they may also have been alluding to the raw image data captured directly from the diffractometer cameras. That traditionally has not been available via the CSD (Cambridge structural database), but would be required for a complete re-analysis of the crystal structure. Now the first example of how both FAIR (processed) data and raw data can co-exist has appeared.

The latest version of the CSD database shows an entry resulting from the following publication[3] and the deposited data has its own DOI there (10.5517/ccdc.csd.cc1n9ppb). That entry in turn has a DOI pointer to the Raw data (10.14469/hpc/2300) held in a different location and the pointer is reciprocated (⇌) with the latter pointing back to the former. Both datasets point to the original article, thus completing a holy triangle.

There is more. The Raw dataset (10.14469/hpc/2300) declares it is a member of a superset, called Crystal structure data for Synthesis and Reactions of Benzannulated Spiroaminals; Tetrahydrospirobiquinolines (10.14469/hpc/2297where you can find information about six other related structures. That collection is in turn a member of a superset called Synthesis and Reactions of Benzannulated Spiroaminals; Tetrahydrospirobiquinolines (10.14469/hpc/2099where DOIs to other types of data associated with this project can be found, such as Computational data (10.14469/hpc/2098) and NMR data (10.14469/hpc/2294). Although a human can with some determination follow these associations up, down and across, the system is designed to also be followed by automated algorithms that could traverse this web quickly and efficiently.

So you can now see that a crystal structure held in the CSD could be the starting point for a journey of FAIR data discovery, in manner that has not hitherto been possible. How quickly the CSD will become populated by links to Raw (and other) data remains to be seen. I have not yet discovered any mechanism for specifying a CSD query which stipulates that Raw data must be available, but no doubt this will come.

To end, back to the Biomimetic Activation of CO2 referred to at the start. With no access to the original data, recourse was made to computational modelling.[2] Which where  I came in, since I wanted access to the original (computational) data. Sadly it did not appear to be available with the article,[2] in much the same manner as the original complaint. Perhaps, when FAIR data becomes fully accepted as part of how science is done nowadays, such complaints will become ever rarer!


In fact the original authors did respond[4] with an acknowledgement that their original conclusions were not correct.

Almost. The article [3] cites DOI: 10.14469/hpc/2099 (Ref 28), but it does not cite DOI: 10.5517/ccdc.csd.cc1n9ppb because the latter had not been minted yet at the time the final proofs were corrected, and there is no mechanism to add it at a later stage.

References

  1. S.L. Ackermann, D.J. Wolstenholme, C. Frazee, G. Deslongchamps, S.H.M. Riley, A. Decken, and G.S. McGrady, "Crystallographic Snapshot of an Arrested Intermediate in the Biomimetic Activation of CO<sub>2</sub>", Angewandte Chemie International Edition, vol. 54, pp. 164-168, 2014. https://doi.org/10.1002/anie.201407165
  2. J. Hurmalainen, M.A. Land, K.N. Robertson, C.J. Roberts, I.S. Morgan, H.M. Tuononen, and J.A.C. Clyburne, "Comment on “Crystallographic Snapshot of an Arrested Intermediate in the Biomimetic Activation of CO<sub>2</sub>”", Angewandte Chemie International Edition, vol. 54, pp. 7484-7487, 2015. https://doi.org/10.1002/anie.201411654
  3. J. Almond-Thynne, A.J.P. White, A. Polyzos, H.S. Rzepa, P.J. Parsons, and A.G.M. Barrett, "Synthesis and Reactions of Benzannulated Spiroaminals: Tetrahydrospirobiquinolines", ACS Omega, vol. 2, pp. 3241-3249, 2017. https://doi.org/10.1021/acsomega.7b00482
  4. S.L. Ackermann, D.J. Wolstenholme, C. Frazee, G. Deslongchamps, S.H.M. Riley, A. Decken, and G.S. McGrady, "Corrigendum: Crystallographic Snapshot of an Arrested Intermediate in the Biomimetic Activation of CO<sub>2</sub>", Angewandte Chemie International Edition, vol. 54, pp. 7470-7470, 2015. https://doi.org/10.1002/anie.201504197

Elongating an N-B single bond is much easier than stretching a C-C single bond.

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.

 
 

Dispersion "bonds" not involving hydrogen. A Cl…Cl candidate?

Thursday, June 29th, 2017

In the previous post, I noted the crystallographic detection of an unusually short non-bonded H…H contact of ~1.5Å, some 0.9Å shorter than twice the van der Waals radius of hydrogen (1.2Å, although some sources quote 1.1Å which would make the contraction ~0.7Å). This was attributed to dispersion attractions accumulating in the rest of the molecule. I asked myself what the potential might be for other elements to reveal significantly contracted non-bonded distances as a result of dispersive attractions.

Here is a simple search of the CSD (Cambridge structure database, query DOI: 10.14469/hpc/2700) for the monovalent halogens as in C-X. The other constraints are R<0.05, no errors, no disorder, normalised hydrogen positions and a non-bonded intermolecular contacts ≥0.4Å shorter than the sum of the van der Waals radii (a recent compendium of atomic values is available here[1]). I should state at the outset that this sort of search for non-bonded contacts is quite effective at revealing errors in the crystallographic determination, despite the search request that there be none. None, that is, that have been identified; there are many that have not been! So this survey only aims to tease out any broad trends, rather than a specific focus on individual systems.

Hits

Comments
15,500 H…H, vdW sum = 2.4Å; The accepted shortest H…H contacts are around 1.5Å; it is likely the majority of entries (all?) with shorter values are crystallographic errors. There is a curious maximum at ~1.78Å which probably relates to the H…H non-bonded distance in CH2 groups, despite the attempted constraint that the interaction be intermolecular.
31 F…F, vdW sum = 2.94Å; there are few examples in total and even fewer with the non-bonded distance contracted by ≥0.5Å. Those with ≥0.8Å are probably crystal errors.
70 Cl…Cl vdW sum =3.5Å. Examples at distances of <2.0Å (a contraction of 1.5Å) are almost certainly errors, even the contractions ~1.0Å are suspect. Does this however indicate Cl is more polarisable than F?
26 Br…Br vdW sum = 3.7Å. Examples contracted by ≥0.4Å are isolated (errors?).
10 I…I vdW sum = 3.96Å. There are no examples contracted by ≥0.5Å; shorter contractions are again errors? 
1164 H…F vdW sum = 2.67Å. The separator between real values and probable errors is ~2.1Å, which indicates contractions of  0.6Å are probably real.
222 H…Cl vdW sum = 2.95Å. The separator between real values and probable errors is ~2.4Å; contractions of  ~0.5Å are probably real.
37 H…Br vdW sum = 3.05Å. The separator between real values and probable errors is ~2.8Å; contractions of  ~0.3Å are probably real.
4 H…I vdW sum = 3.18Å. More data is needed!

We see from these results that the number of short H…H contacts far exceeds those found for the halogen series. Clearly whilst the electron density surrounding H is low, that for the halogens is far higher and hence that electron repulsions are going to be far greater. The effect is attenuated if one partner is H itself, with van der Waals contractions in-between those for X…X  and H…H contacts.

This brief survey suggests that  H…H distances are by far the best for probing close contacts brought about by dispersion attractions. Perhaps the next element to focus on for such effects might be chlorine rather than fluorine. For example, see DOI: 10.5517/cczjyhx, being Ph3C-Cl…Cl-CPh3 where the Cl…Cl distance is contracted by ~0.3Å. What would happen if the Ph groups were adorned with t-butyl groups to increase the dispersion attractions?


It might of course also mean that the van der Waals radius for H is set too high.

References

  1. S. Alvarez, "A cartography of the van der Waals territories", Dalton Transactions, vol. 42, pp. 8617, 2013. https://doi.org/10.1039/c3dt50599e

Dispersion “bonds” not involving hydrogen. A Cl…Cl candidate?

Thursday, June 29th, 2017

In the previous post, I noted the crystallographic detection of an unusually short non-bonded H…H contact of ~1.5Å, some 0.9Å shorter than twice the van der Waals radius of hydrogen (1.2Å, although some sources quote 1.1Å which would make the contraction ~0.7Å). This was attributed to dispersion attractions accumulating in the rest of the molecule. I asked myself what the potential might be for other elements to reveal significantly contracted non-bonded distances as a result of dispersive attractions.

Here is a simple search of the CSD (Cambridge structure database, query DOI: 10.14469/hpc/2700) for the monovalent halogens as in C-X. The other constraints are R<0.05, no errors, no disorder, normalised hydrogen positions and a non-bonded intermolecular contacts ≥0.4Å shorter than the sum of the van der Waals radii (a recent compendium of atomic values is available here[1]). I should state at the outset that this sort of search for non-bonded contacts is quite effective at revealing errors in the crystallographic determination, despite the search request that there be none. None, that is, that have been identified; there are many that have not been! So this survey only aims to tease out any broad trends, rather than a specific focus on individual systems.

Hits

Comments
15,500 H…H, vdW sum = 2.4Å; The accepted shortest H…H contacts are around 1.5Å; it is likely the majority of entries (all?) with shorter values are crystallographic errors. There is a curious maximum at ~1.78Å which probably relates to the H…H non-bonded distance in CH2 groups, despite the attempted constraint that the interaction be intermolecular.
31 F…F, vdW sum = 2.94Å; there are few examples in total and even fewer with the non-bonded distance contracted by ≥0.5Å. Those with ≥0.8Å are probably crystal errors.
70 Cl…Cl vdW sum =3.5Å. Examples at distances of <2.0Å (a contraction of 1.5Å) are almost certainly errors, even the contractions ~1.0Å are suspect. Does this however indicate Cl is more polarisable than F?
26 Br…Br vdW sum = 3.7Å. Examples contracted by ≥0.4Å are isolated (errors?).
10 I…I vdW sum = 3.96Å. There are no examples contracted by ≥0.5Å; shorter contractions are again errors? 
1164 H…F vdW sum = 2.67Å. The separator between real values and probable errors is ~2.1Å, which indicates contractions of  0.6Å are probably real.
222 H…Cl vdW sum = 2.95Å. The separator between real values and probable errors is ~2.4Å; contractions of  ~0.5Å are probably real.
37 H…Br vdW sum = 3.05Å. The separator between real values and probable errors is ~2.8Å; contractions of  ~0.3Å are probably real.
4 H…I vdW sum = 3.18Å. More data is needed!

We see from these results that the number of short H…H contacts far exceeds those found for the halogen series. Clearly whilst the electron density surrounding H is low, that for the halogens is far higher and hence that electron repulsions are going to be far greater. The effect is attenuated if one partner is H itself, with van der Waals contractions in-between those for X…X  and H…H contacts.

This brief survey suggests that  H…H distances are by far the best for probing close contacts brought about by dispersion attractions. Perhaps the next element to focus on for such effects might be chlorine rather than fluorine. For example, see DOI: 10.5517/cczjyhx, being Ph3C-Cl…Cl-CPh3 where the Cl…Cl distance is contracted by ~0.3Å. What would happen if the Ph groups were adorned with t-butyl groups to increase the dispersion attractions?


It might of course also mean that the van der Waals radius for H is set too high.

References

  1. S. Alvarez, "A cartography of the van der Waals territories", Dalton Transactions, vol. 42, pp. 8617, 2013. https://doi.org/10.1039/c3dt50599e

Chemistry rich diagrams: do crystal structures carry spin information? Iron-di-imine complexes.

Sunday, June 18th, 2017

The iron complex shown below forms the basis for many catalysts.[1] With iron, the catalytic behaviour very much depends on the spin-state of the molecule, which for the below can be either high (hextet) or medium (quartet) spin, with a possibility also of a low spin (doublet) state. Here I explore whether structural information in crystal structures can reflect such spin states.

We studied this a few years back and the talk I gave on the topic included some of our first statistical explorations of the CSD (Cambridge structure database). Here I update those searches, using the search query (DOI: 10.14469/hpc/2675) shown below. The di-imine ligand contains only 3-coordinate atoms, whilst the iron is 5-coordinate. The angles subtended at the Fe and group X=NM (any non-metal atom) are as defined below.

The resulting scatterplot is shown below and contains a rich variety of phenomena.

  1. In the bond length region of 1.85-1.95Å one sees three clusters, one arranged on the diagonal indicating both N-Fe lengths are the same and two off the diagonal which indicates one length is ~0.1Å longer than the other.
    • To explain this, one needs to know that 5-coordinate Fe has a trigonal bipyramidal shape in which one X=NM group subtends an (anti-periplanar)  angle of ~180° at Fe with one of the ring nitrogens and the other two X=NM groups each subtend an angle of <120° with the other ring nitrogen. The result is that if the group X has the appropriate (electron withdrawing) properties, the two N-Fe bond lengths are no longer equal. If group X is more passive, the two N-Fe bond lengths may remain more equal.
  2. A second cluster occurs at ~2.00-2.1Å, mostly along the diagonal but with hints of smaller off-diagonal clusters.
  3. A third feature occurs at ~2.1-2.3Å, where now the off-diagonal clusters contain more examples than are on the diagonal itself.

Clearly, there is more going on here than can be explained simply by the orientation of X=NM with respect to the Fe-N bond axis. That something is the spin-multiplicity of the molecule. With the Fe complex shown above, this can be one of doublet (one unpaired electron), quartet (three unpaired electrons) or hextet (five unpaired electrons). To gain insight into how this affects the bond lengths, some calculations are needed, using X=Cl, R=H. Here they are done at the TPSSH/Def2-TZVPP level. In fact it is well-known[2] that the energy separations of low/medium/high spin Fe complexes are highly sensitive to the functional, but TPSSH seems to be amongst the best. This shows that the energy ordering of the three states using this particular method is hextet (0.0, DOI: 10.14469/hpc/2676) < quartet (10.5, DOI: 10.14469/hpc/2677) < doublet (13.2 kcal/mol, DOI: 10.14469/hpc/2678), with the bond lengths shown below (for X=Cl).

We might make tentative hypotheses based on these values:

  1. The off-diagonal bottom left clusters (1 in list above) might arise from doublet states.
  2. The off-diagonal top right clusters (3 in list above) might arise from sextet states.
  3. The cluster (2 in list above) might be quartet states for which X is not sufficiently electronegative to induce bond length discriminations.
  4. It is worth noting that the energy span between the three states for the above molecule is only ~13 kcal/mol, which is small enough to be altered by substituents.

Testing these hypotheses requires knowledge of the spin state of all the entries in any cluster. This information is unfortunately not carried by the CSD, which has relatively little information over and above structural data. Each entry would have to be individually inspected. Indeed the spin state of many of these complexes may not even be known. Nevertheless, it would be great to repeat the graphs shown above as a function of known spin state so that the (again I repeat tentative) hypotheses might be confirmed or refuted.


This article evaluates a whole host of functionals against e.g. the spin-state energy separations of the Fe2+ ion. As it happens, TPSSH was not one that was evaluated, but in fact it gives more or less the best match to experiment. Thus Esinglet-Equintet obs = 85.6 kcal/mol, calc 92.4; Etriplet-Equintet obs 56.1, calc 59.5 kcal/mol. A hypothesis therefore is that the TPSSH functional is a reasonable one to go exploring such high-spin species.

References

  1. M.P. Shaver, L.E.N. Allan, H.S. Rzepa, and V.C. Gibson, "Correlation of Metal Spin State with Catalytic Reactivity: Polymerizations Mediated by α‐Diimine–Iron Complexes", Angewandte Chemie International Edition, vol. 45, pp. 1241-1244, 2006. https://doi.org/10.1002/anie.200502985
  2. P. Verma, Z. Varga, J.E.M.N. Klein, C.J. Cramer, L. Que, and D.G. Truhlar, "Assessment of electronic structure methods for the determination of the ground spin states of Fe(<scp>ii</scp>), Fe(<scp>iii</scp>) and Fe(<scp>iv</scp>) complexes", Physical Chemistry Chemical Physics, vol. 19, pp. 13049-13069, 2017. https://doi.org/10.1039/c7cp01263b