Posts Tagged ‘Nitrogen’

Hypervalent or not? A fluxional triselenide.

Saturday, February 24th, 2018

Another post inspired by a comment on an earlier one; I had been discussing compounds of the type I.In (n=4,6) as possible candidates for hypervalency. The comment suggests the below as a similar analogue, deriving from observations made in 1989.[1]

This compound was investigated using 77Se NMR, with the following conclusions:

  1. The compound is fluxional, with the lines at room temperature broadened compared to those at -50°C.
  2. At -50°C the peaks are sharp enough to discern 1JSe-Se couplings, with multiplicities and integrations that suggest a central Se is surrounded by four equivalent further Se atoms, with shifts of 655.1 and 251.2 ppm.
  3. The magnitude of this 1JSe-Se coupling (391 Hz) leads to the suggestion of a considerable contribution of a resonance form with Se=Se bonds (structure 2 above).
  4. This was supported by 2J13C-77Se couplings which also imply a symmetrically coordinated central  Se.
  5. Thus the two resonance forms 1 or 2 above were suggested as the predominant form at -50°C, with an increasing incursion of the open chain isomer 3 at higher temperatures giving rise to the observed fluxional dynamic behaviour.
  6. One may surmise from these results that the central Se is certainly hypercoordinated and by the classical interpretations hypervalent.

Here are some calculations (R=H), at the ωB97XD/Def2-TZVPP/SCRF=chloroform level.‡ In red are the calculated Wiberg Se-Se bond orders, which give little indication of any Se=Se double bond character. 

The calculated 77Se shifts are shown in magenta, with the observed values being 655 and 255 ppm. The match is not good, the errors were 120 and 20.5 ppm.†  However calculated shifts for elements adjacent to e.g. Se or Br etc suffer from relativistic effects such as spin orbit coupling.[2] Thus the shift for the central Se, surrounded by four other Se atoms is likely to have a significant error, but the error for the four other Se atoms should be less. The reverse is true.

However, all the calculations of this species (up to Def2-TZVPPD basis set) showed this symmetric form of D2h symmetry to actually be a transition state, as per below.

There is a minimum with the structure below in which one pair of Se-Se lengths are longer than the other pair and for which the free energy is 6.5 kcal/mol lower. The Wiberg bond orders for the two sets of Se-Se bonds are now 0.16 and 0.86, which very much corresponds to structure 3 above.

Assuming that this compound is fluxional even at -50°C, the average of the pairs of Se atoms gives calculated shifts of 667 ppm (655 obs) whilst the central Se is 204.6 ppm (251 obs). The latter, influenced by two especially short Se-Se distances, is likely to have a very large spin-orbit coupling error, whilst for the former the error will be smaller (13C shifts adjacent to one Br typically have induced calculated errors of about 14 ppm[2]).

At this point I searched the Cambridge structure database for Se coordinated by four other Se atoms. A close analogue[3] has the structure shown below, in which pairs of Se-Se interactions have unequal bond lengths, the shorter being ~2.45Å. This matches the calculation above reasonably well.

Reconciling these various observations, we might assume that even at -50°C the fluxional behaviour has not been frozen out. Given that the fluxional barrier is only 6.5 kcal/mol, it is unlikely that the spectrum could be measured at a sufficiently low temperature to reveal not two sets of Se signals in the ratio 4:1 but three in the ratio 2:2:1. The spin-spin couplings reported presumably are a result of averaging a genuine 1JSe-Se coupling with a through space coupling.

So it appears that the analysis of the 77Se NMR reported in this article [1] may not be quite what it seems. A better interpretation is that structure 3 is the most realistic. This means no hypercoordination for the Se, never mind hypervalence!

‡FAIR data at DOI: 10.14469/hpc/3724. †The original reference, Me2Se was incorrectly calculated without solvation by chloroform. The values shown here are now corrected from those shown in the original post.

References

  1. Y. Mazaki, and K. Kobayashi, "Structure and intramolecular dynamics of bis(diisobutylselenocarbamoyl) triselenide as identified in solution by the 77Se-NMR spectroscopy", Tetrahedron Letters, vol. 30, pp. 2813-2816, 1989. https://doi.org/10.1016/s0040-4039(00)99132-9
  2. D.C. Braddock, and H.S. Rzepa, "Structural Reassignment of Obtusallenes V, VI, and VII by GIAO-Based Density Functional Prediction", Journal of Natural Products, vol. 71, pp. 728-730, 2008. https://doi.org/10.1021/np0705918
  3. R.O. Gould, C.L. Jones, W.J. Savage, and T.A. Stephenson, "Crystal and molecular structure of bis(NN-diethyldiselenocarbamato)-selenium(II)", Journal of the Chemical Society, Dalton Transactions, pp. 908, 1976. https://doi.org/10.1039/dt9760000908

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

Ammonide: an alkalide formed from ammonia and resembling an electride.

Sunday, December 17th, 2017

Alkalides are anionic alkali compounds containing e.g. sodide (Na–), kalide (K–), rubidide (Rb–) or caeside (Cs–). Around 90 examples can be found in the Cambridge structure database (see DOI: 10.14469/hpc/3453  for the search query and results). So what about the ammonium analogue, ammonide (NH4–)? A quick search of Scifinder drew a blank! So here I take a look at this intriguingly simple little molecule.‡

It can be formed by adding two electrons to the ammonium cation; NH4+ + 2e ↠ NH4–. One might be encouraged to do this since the LUMO (lowest unoccupied molecular orbital, below) of the ammonium cation has A1 symmetry and so can accept two electrons without the penalty of Jahn-Teller distortions. These electrons will apparently expand the valence electron “octet” around the nitrogen from 8 to 10; a hypervalent species then?

So what are the (calculated) properties of NH4–? The energy of the now HOMO (highest occupied molecular orbital) at the ωB97XD/Def2-TZVPPD/solvent=water level is -3.6eV, a respectable electron affinity (the additional electrons are said to be bound). More insight can be obtained from the NBO analysis, which produces localized versions of the molecular orbitals. There are four equivalent NBOs, one of which is shown below.

Each is bonding along one H-N bond, mildly anti-bonding along the other three N-H bonds, but again bonding in the H-H regions! This matches the observations made earlier that when more electrons are pumped into normally valent main group molecules, they will occupy the antibonding levels. This is accompanied by a reduction in the bond orders associated with the central atom. In this case, the N-H bond orders are reduced from 0.79 to 0.602 and the total bond index at the nitrogen is reduced from 3.16 to 2.408. The bond index at hydrogen is at first sight increased from 0.79 to a surprising 1.0003, but this is explained because the H-H bond orders are 0.1328 (three for each H), which brings the H index up to 1.0. The N-H vibration (A1 symmetric) is 3466 cm-1 for NH4+  which is reduced to 2659 for NH4–.†

So it appears that adding two electrons to the ammonium cation induces H-H bonding! More insight can be obtained from an ELF analysis of the electron density basins.

The above shows four attractors (as they are called) centered at the hydrogen nuclei, with 2.053e each (4*2.053 = 8.212e). The remaining ~2e are located in basins (green) centered at two different types of attractors. One is along the axis of each N-H bond and exo to it (0.316e) and the other sits on top of any set of three hydrogens (0.103e), 1.68e in total. The value of the ELF function at the attractor is shown above. You should realize that ELF=1.0 corresponds to perfectly localized electrons (for which the kinetic energy density is zero) and ELF=0.5 would correspond to a free-electron gas. The ELF value of e.g. 0.77 is getting close to an electron gas, and in fact corresponds to what we call an electride.

So, the nitrogen valence shell electron octet is not actually exceeded! The additional two electrons in ammonide sit beyond the nitrogen, in H-H regions. Whether or not it is a viable species for detection remains to be established, but even its computed bonding properties have proved interesting and it deserves to join the alkalide family. 

Postscript

Ammonide exists in a shallow well in the potential energy surface, shown below, with the dissociation to ammonia and hydride anion being exothermic.

The intrinsic reaction coordinate shows one interesting feature at  IRC ~-1.1 which corresponds to repulsion between the hydride and the lone pair of the nitrogen atom resulting in inversion of configuration during the latter stages of the IRC.

‡FAIR data collection; 10.14469/hpc/3455. †Perhaps unsurprisingly, these values are somewhat basis set dependent. Thus a ωB97XD/Def2-QZVPPD/Water calculation gives the following values: bond index at N, 1.998, N-H bond index, 0.4995, H-H bond index 0.1689, H bond index 1.0062, total Rydberg population, 0.2738, ν(A1) 2686 cm-1. The ELF basins are H, 2.039, exo-basins 0.282 and 0.141 (total 1.692). The improved basis set better describes the diffuse regions beyond the N-H bonds.

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

First, hexacoordinate carbon – now pentacoordinate oxygen?

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.

  1. 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.
  2. as does the system  from cyclobutadienyl di-anion and O=P4+(DOI: 10.14469/hpc/2358)
  3. and the system from cyclobutadienyl di-anion and HS5+(DOI: 10.14469/hpc/2360).
  4. 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:, , , , ,
Posted in Bond slam, Hypervalency | 4 Comments »

Anomeric effects at carbon involving lone pairs originating from one or two nitrogens.

Friday, July 8th, 2016

The previous post looked at anomeric effects set up on centres such as B, Si or P, and involving two oxygen groups attached to these atoms. Here I vary the attached groups to include either one or two nitrogen atoms.[1]

.aminol-sq

The plot below shows aminols, C(NHR)(OR”). A torsion along either the C-O or C-N bond of ~60° implies that (at two coordinate oxygen or three coordinated nitrogen) there may be a lone pair with a torsion of 180°, which would set up an antiperiplanar alignment between that lone pair and the adjacent C-O or C-N bond (the anomeric effect). The clear hotspot is at angles of ~80°, which does raise the issue of why it deviates from 60°. Only a location of the lone pair centroid (using eg the ELF quantum mechanical technique) would cast light on that. There is a less distinct region for which the C-N torsion is 60° and the C-O torsion 180°, and an even less distinct region for the reverse (C-O torsion is 60° and the C-N torsion 180°). This tends to imply that a nitrogen lone pair is a better donor into a C-O bond than the reverse. Electronegativity suggests this should indeed be so, with the N lone pair less bound by the N nucleus and hence easier to release into a C-Oσ* orbital which is a better acceptor than then equivalent C-Nσ* orbital. aminol This plot is where both heteroatoms are nitrogen (geminal diamines). There are about twice as many examples, resulting in more distinct clustering. The anomeric hotspot is now around 70°  and there are equally populated clusters where only one torsion is ~70°. There is another cluster for which both torsions are 180° (no stereoelectronic alignment of lone pairs) and three small clusters where the torsions are either 180° or 0°. There is finally an intriguing cluster for which both torsions at ~120° (again no stereoelectronics). diamine

Searches like this seem to be good at creating more questions than they answer. Clearly, the origins of the various hotspots need to be investigated, ideally using quantum mechanics to quantify the stereoelectronic interactions involved. So this sort of (ten minute) exercise is very good at raising research project investigations.

References

  1. H. Rzepa, "Anomeric effects at carbon, involving lone pairs originating from one or two nitrogens", 2016. https://doi.org/10.14469/hpc/936

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