Posts Tagged ‘Organometallic chemistry’

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

Tags:, , , , , , , , , , ,
Posted in crystal_structure_mining | 1 Comment »

What’s in a name? Carbenes: a reality check.

Sunday, September 11th, 2016

To quote from Wikipedia: in chemistry, a carbene is a molecule containing a neutral carbon atom with a valence of two and two unshared valence electrons. The most ubiquitous type of carbene of recent times is the one shown below as 1, often referred to as a resonance stabilised or persistent carbene. This type is of interest because of its ability to act as a ligand to an astonishingly wide variety of metals, with many of the resulting complexes being important catalysts. The Wiki page on persistent carbenes shows them throughout in form 1 below, thus reinforcing the belief that they have a valence of two and by implication six (2×2 shared + 2 unshared) electrons in the valence shell of carbon. Here I consider whether this name is really appropriate.

carbenes

Let us start by counting the electrons in the 2p atomic orbitals on the ring atoms of 1, forming what we call a π-system. There are six; two from the carbons shown connected by a double bond, C=C and a further four from the two nitrogen lone pairs. Now in benzene, we also have six π-electrons in a ring and this molecule is of course famously aromatic due to the diatropic ring current created by the circulation of these six electrons. Moreover, all the C-C bonds are equal in length, ~1.4Å long (although the reasons for this equality are subtle).

So does 1 behave similarly? A ωB97XD/Def2-TZVPP calculation[1] shows the following calculated structure, in which all the bonds are clearly intermediate between single and double. The N-C(“carbene”) length of 1.357Å in particular is much shorter than a C-N single bond (~1.44AÅ), which tends to suggest that the resonance form 2 is a better representation than 1. This form is also pretty similar to pyrrole, itself a well-known hetero-aromatic species.nhc1

An alternative reality check is crystal structures. There are 42 examples (no errors, no disorder, R < 0.05) in the Cambridge structure database (CSD) and the distribution of C-N bond lengths below is indeed quite similar to the calculation shown above for the unsubstituted parent, with the lhs “hot-spot” almost exactly coincident. The C-C length similarly corresponds.

nhc2

nhc3

Let us try a technique for explicitly counting electrons, the ELF (electron localisation method), which works directly on a function of the electron density to identify the centroids of localized “basins” containing the integrated density. The three surrounding the “carbene” atom sum to 7.54e (with small seepage also into the carbon 1s core; 2.08e). A “normal” carbon on the C=C bond is 7.65e. The localization below turns out to closely resemble resonance structure 2 above.

nhc4

Further in-silico experiments can be carried out with species 3 and 4, in which a carbon atom replaces each of the nitrogens. This reduces the total electron count by two and now this poor molecule has a difficult choice to make. Should it be the π-system that sacrifices these two electrons, or could it be the σ-lone-pair found on the two-coordinate carbon? We will let the quantum mechanical solution decide[2] (with a constraint that the molecule be planar). The electrons arrange themselves to resemble the resonance form 4, choosing to retain the six π-electrons and sacrifice the carbene “unshared pair”. The 2-coordinate carbon as a vinyl cation now does have ~6 valence electrons (ELF indicates 5.23e). nhc5

What about the other choice? By promoting two electrons from HOMO to LUMO one can also calculate 3 (again constrained to planarity)[3] which finally does correspond to the classical description of a carbene.

nhc6

The arrow connecting 3 and 4 in the scheme at the top is NOT in this case an electronic resonance, but a a real equilibrium between two different species separated by an energy barrier. With only four π-electrons in a cycle it is also antiaromatic, and so the two localised alkene bonds avoid any conjugation with each other. This form has a free energy some 5.7 kcal/ml higher than the aromatic form. In fact, the molecule is very keen to avoid all antiaromaticity and hence if the planar constraint is lifted, it will distort with no activation to a non-planar diene (just as cyclo-octatetraene does to a non-planar tetra-ene). And to complete the tale, even though 4 is aromatic, it too distorts without activation to an odd-looking non-planar form with no symmetry[4],[5],[6] (but that is another story).

The final word should be that the naming of these types of persistent carbene does need a reality check; they should not be called this at all! They are really dipolar species or carbon-ylides as shown in 2. As it happens, a very closely related species in which one sulfur replaces one nitrogen is a very familiar compound, vitamin B1 or thiamine. The only example of a stable deprotonated thiamine derivative is referred to as a carbene[7], perhaps because with an acid catalyst it can dimerise in the manner expected of a real carbene. Significantly however, without acid catalyst this does not happen; a true carbene would not require such a catalyst.

References

  1. H. Rzepa, "NHC wfn", 2016. https://doi.org/10.14469/hpc/1473
  2. H. Rzepa, "butadiene carbene aromatic -192.700746", 2016. https://doi.org/10.14469/hpc/1581
  3. H. Rzepa, "butadiene carbene antiaromatic guess=alter -192.691607", 2016. https://doi.org/10.14469/hpc/1582
  4. H. Rzepa, "C5H4 non-planar, Cs symmetry", 2016. https://doi.org/10.14469/hpc/1583
  5. H. Rzepa, "C5H4 non-planar, C2 symmetry", 2016. https://doi.org/10.14469/hpc/1584
  6. H. Rzepa, "C5H4 non-planar, no symmetry", 2016. https://doi.org/10.14469/hpc/1585
  7. A.J. Arduengo, J.R. Goerlich, and W.J. Marshall, "A Stable Thiazol‐2‐ylidene and Its Dimer", Liebigs Annalen, vol. 1997, pp. 365-374, 1997. https://doi.org/10.1002/jlac.199719970213

Tags:, , , , , , , , , , , , ,
Posted in crystal_structure_mining, General | No Comments »

What's in a name? Carbenes: a reality check.

Sunday, September 11th, 2016

To quote from Wikipedia: in chemistry, a carbene is a molecule containing a neutral carbon atom with a valence of two and two unshared valence electrons. The most ubiquitous type of carbene of recent times is the one shown below as 1, often referred to as a resonance stabilised or persistent carbene. This type is of interest because of its ability to act as a ligand to an astonishingly wide variety of metals, with many of the resulting complexes being important catalysts. The Wiki page on persistent carbenes shows them throughout in form 1 below, thus reinforcing the belief that they have a valence of two and by implication six (2×2 shared + 2 unshared) electrons in the valence shell of carbon. Here I consider whether this name is really appropriate.

carbenes

Let us start by counting the electrons in the 2p atomic orbitals on the ring atoms of 1, forming what we call a π-system. There are six; two from the carbons shown connected by a double bond, C=C and a further four from the two nitrogen lone pairs. Now in benzene, we also have six π-electrons in a ring and this molecule is of course famously aromatic due to the diatropic ring current created by the circulation of these six electrons. Moreover, all the C-C bonds are equal in length, ~1.4Å long (although the reasons for this equality are subtle).

So does 1 behave similarly? A ωB97XD/Def2-TZVPP calculation[1] shows the following calculated structure, in which all the bonds are clearly intermediate between single and double. The N-C(“carbene”) length of 1.357Å in particular is much shorter than a C-N single bond (~1.44AÅ), which tends to suggest that the resonance form 2 is a better representation than 1. This form is also pretty similar to pyrrole, itself a well-known hetero-aromatic species.nhc1

An alternative reality check is crystal structures. There are 42 examples (no errors, no disorder, R < 0.05) in the Cambridge structure database (CSD) and the distribution of C-N bond lengths below is indeed quite similar to the calculation shown above for the unsubstituted parent, with the lhs “hot-spot” almost exactly coincident. The C-C length similarly corresponds.

nhc2

nhc3

Let us try a technique for explicitly counting electrons, the ELF (electron localisation method), which works directly on a function of the electron density to identify the centroids of localized “basins” containing the integrated density. The three surrounding the “carbene” atom sum to 7.54e (with small seepage also into the carbon 1s core; 2.08e). A “normal” carbon on the C=C bond is 7.65e. The localization below turns out to closely resemble resonance structure 2 above.

nhc4

Further in-silico experiments can be carried out with species 3 and 4, in which a carbon atom replaces each of the nitrogens. This reduces the total electron count by two and now this poor molecule has a difficult choice to make. Should it be the π-system that sacrifices these two electrons, or could it be the σ-lone-pair found on the two-coordinate carbon? We will let the quantum mechanical solution decide[2] (with a constraint that the molecule be planar). The electrons arrange themselves to resemble the resonance form 4, choosing to retain the six π-electrons and sacrifice the carbene “unshared pair”. The 2-coordinate carbon as a vinyl cation now does have ~6 valence electrons (ELF indicates 5.23e). nhc5

What about the other choice? By promoting two electrons from HOMO to LUMO one can also calculate 3 (again constrained to planarity)[3] which finally does correspond to the classical description of a carbene.

nhc6

The arrow connecting 3 and 4 in the scheme at the top is NOT in this case an electronic resonance, but a a real equilibrium between two different species separated by an energy barrier. With only four π-electrons in a cycle it is also antiaromatic, and so the two localised alkene bonds avoid any conjugation with each other. This form has a free energy some 5.7 kcal/ml higher than the aromatic form. In fact, the molecule is very keen to avoid all antiaromaticity and hence if the planar constraint is lifted, it will distort with no activation to a non-planar diene (just as cyclo-octatetraene does to a non-planar tetra-ene). And to complete the tale, even though 4 is aromatic, it too distorts without activation to an odd-looking non-planar form with no symmetry[4],[5],[6] (but that is another story).

The final word should be that the naming of these types of persistent carbene does need a reality check; they should not be called this at all! They are really dipolar species or carbon-ylides as shown in 2. As it happens, a very closely related species in which one sulfur replaces one nitrogen is a very familiar compound, vitamin B1 or thiamine. The only example of a stable deprotonated thiamine derivative is referred to as a carbene[7], perhaps because with an acid catalyst it can dimerise in the manner expected of a real carbene. Significantly however, without acid catalyst this does not happen; a true carbene would not require such a catalyst.

References

  1. H. Rzepa, "NHC wfn", 2016. https://doi.org/10.14469/hpc/1473
  2. H. Rzepa, "butadiene carbene aromatic -192.700746", 2016. https://doi.org/10.14469/hpc/1581
  3. H. Rzepa, "butadiene carbene antiaromatic guess=alter -192.691607", 2016. https://doi.org/10.14469/hpc/1582
  4. H. Rzepa, "C5H4 non-planar, Cs symmetry", 2016. https://doi.org/10.14469/hpc/1583
  5. H. Rzepa, "C5H4 non-planar, C2 symmetry", 2016. https://doi.org/10.14469/hpc/1584
  6. H. Rzepa, "C5H4 non-planar, no symmetry", 2016. https://doi.org/10.14469/hpc/1585
  7. A.J. Arduengo, J.R. Goerlich, and W.J. Marshall, "A Stable Thiazol‐2‐ylidene and Its Dimer", Liebigs Annalen, vol. 1997, pp. 365-374, 1997. https://doi.org/10.1002/jlac.199719970213

Tags:, , , , , , , , , , , , ,
Posted in crystal_structure_mining, General | No Comments »