Posts Tagged ‘Simple aromatic rings’

Cyclopropenium cyclopentadienide: a strangely neutral ion-pair?

Sunday, April 9th, 2017

Both the cyclopropenium cation and the cyclopentadienide anion are well-known 4n+2-type aromatic ions, but could the two together form an ion-pair?

A search of the Cambridge structure database reveals 52 instances of the cyclopropenium cation with a variety of counter-anions, 77 cyclopentadienide anions with a variety of counter-cations and one (SOWMOG, private communication to CSD) where the two sub-structures are common. The pyridinium-cyclopropenium fragment is actually a di-cation stabilized with dimethylamino substituents, with these charges balanced by two cyclopentadienide anions stabilized with ester substituents. The stacking distance between the ion-pairs is ~3.5-3.6Å, a bit larger than normal π-π stacking distances of 3.2-3.3Å

So could a “pure” cyclopropenium cyclopentadienide ion-pair exist, and if so what would its π-π stacking distance be? A ωB97XD/Def2-TZVPPD/SCRF=water calculation (DOI: 10.14469/hpc/2442) provides one answer to this question; 2.57Å!‡ It is a true minimum in the potential energy surface (all +ve force constants) with a calculated dipole moment of only 7.57D. This species is “only” 27.1 kcal/mol higher in ΔG than the neutral hydrocarbon (DOI: 10.14469/hpc/2443), a difference which is as low as it is because of the gain in aromatic stabilization of two rings upon ion-pair formation.

A few posts back, I was considering candidates for the most polar neutral compound synthesized and I suggested a candidate with a dipole moment of ~22D, based as it happens on cyclopropenium and cyclopentadienide rings directly connected by a bond. So when this bond is removed and the two rings are allowed to stack one above the other, we now have an interesting inversion of the original challenge: what is the least-polar ionic organic compound (ionic in the sense of being an unconnected ion-pair)?

Here are some more properties of this intriguing “neutral” ion-pair.

  1. It has a number of low-frequency modes with correspond to the two rings moving with respect to each other (ν 216 cm-1)
  2. The molecular electrostatic potential illustrates the sense of polarization, with negative region (orange) residing on the 5-membered ring:
  3. The most stable π-type molecular orbital (below) reminds of the π-complex formed in the benzidine rearrangement and that in fact modelling this ion-pair may require a multi-reference (CASSCF) wavefunction, with the single-determinantal one used here only being a first approximation.
  4. A QTAIM analysis of the electron density topology shows only weak “bond” connectors between the two rings, with ρ(r) being typical of weak interactions such as hydrogen bonds.
  5. An ELF (electron localisation function) analysis also holds no surprises, with all the electron density basins (purple) confined to the two rings, just as expected of an ion-pair.
  6. I will leave one further question to a future discussion; what happens to the aromaticity and ring currents of the two individual rings as they combine to form this ion-pair? Might this property be connected to the very close separation between the two rings?

So we have a remarkably “neutral” ionic hydrocarbon to match the “ionic” neutral organic molecules previously discussed. This ion-pair may yet prove to have interesting properties, even if is unlikely to be synthesized without the addition of stabilising substituents.

‡ For example, the stacking distance in graphite is 3.35Å.

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

How does an OH or NH group approach an aromatic ring to hydrogen bond with its π-face?

Wednesday, June 22nd, 2016

I previously used data mining of crystal structures to explore the directing influence of substituents on aromatic and heteroatomatic rings. Here I explore, quite literally, a different angle to the hydrogen bonding interactions between a benzene ring and OH or NH groups.

aromatic-pi-query

I start by defining a benzene ring with a centroid. The distance is from that centroid to the H atom of an OH or NH group and the angle is C-centroid-H. To limit the search to approach of the OH or NH group more or less orthogonal to the ring, the absolute value of the torsion between the centroid-H vector and the ring C-C vector is constrained to lie between 70-100° (the other constraints being no disorder, no errors, T < 140K and R < 0.05).[1]

aromatic-pi-HN-140

The above shows the results for NH groups interacting with the aromatic ring. The maximum distance 2.8Å is more or less the van der Waals contact distance between a hydrogen and a carbon and as you can see the contacts "funnel down" to the centroid at < 2.1Å. The shortest distance[2] is for ammonium tetraphenylborate, which you can view in e.g. spacefill mode here[3]

390

The other interesting close contact derives from a protonated pyridine[4], which can in turn be viewed here.[5] The main message from the distribution shown above is that as the distances between the HN and the centroid get shorter, the "trajectory" of approach remains orthogonal to the ring (the angle defined above remains ~90°) and heads towards the centroid of the π-cloud. The hotspot itself (red, ~2.6Å) also lies along this trajectory.

Recollect that when I used such hydrogen bonding to see if crystal structures discriminate between the ortho or meta positions of a ring carrying an electron donating substituent, it was the distance from a HO to the carbon that was measured as the discriminator. So it's a faint surprise to find that with HN, and without the necessary perturbation of an electron donating substituent, the intrinsic preference seems to be for the ring centroid and not any specific carbon atom of the ring.

So how about the OH group? There are in fact rather fewer examples, and so the statistics are a bit less clear-cut. But there is a tantalising suggestion that this time, the trajectory is not ~90° but rather less, implying that the destination is no longer the centroid of the π-cloud but one of the carbon atoms of the ring itself. For those who like to "read between the lines" and spot things that are absent rather than present, you may have asked yourself why I did not use NH probes in my earlier post. Well, it appears that the NH group is less effective at e.g. o/p discrimination than is an OH group.

aromatic-pi-OH-140

I can only speculate as to the origins (real or not) of the difference in behaviour between OH and NH groups towards a phenyl π-face. Perhaps it is simply bias in the CSD database? Or might there be electronic origins? Time to end with that phrase "watch this space".

 

References

  1. H. Rzepa, "How does an OH or NH group approach an aromatic ring to hydrogen bond with its π-face?", 2016. https://doi.org/10.14469/hpc/673
  2. T. Steiner, and S.A. Mason, "Short N<sup>+</sup>—H...Ph hydrogen bonds in ammonium tetraphenylborate characterized by neutron diffraction", Acta Crystallographica Section B Structural Science, vol. 56, pp. 254-260, 2000. https://doi.org/10.1107/s0108768199012318
  3. Steiner, T.., and Mason, S.A.., "CCDC 144361: Experimental Crystal Structure Determination", 2000. https://doi.org/10.5517/cc4v6tz
  4. O. Danylyuk, B. Leśniewska, K. Suwinska, N. Matoussi, and A.W. Coleman, "Structural Diversity in the Crystalline Complexes of <i>para</i>-Sulfonato-calix[4]arene with Bipyridinium Derivatives", Crystal Growth & Design, vol. 10, pp. 4542-4549, 2010. https://doi.org/10.1021/cg100831c
  5. Danylyuk, O.., Lesniewska, B.., Suwinska, K.., Matoussi, N.., and Coleman, A.W.., "CCDC 819118: Experimental Crystal Structure Determination", 2011. https://doi.org/10.5517/ccwhc5w

Tags:, , , , , , , , , , , , ,
Posted in Chemical IT, crystal_structure_mining | 3 Comments »

Exploring the electrophilic directing influence of heteroaromatic rings using crystal structure data mining.

Tuesday, June 21st, 2016

This is a follow-up to the post on exploring the directing influence of (electron donating) substituents on benzene[1] with the focus on heteroaromatic rings such indoles, pyrroles and group 16 analogues (furans, thiophenes etc).

s-cis-ester1

The search query is shown above (and is available here[2]). As before, the distance is compared from an electrophile, modelled as the hydrogen atom of an OH group, to both the carbon next to the heteroatom (C2) and the C3 carbon. The torsion is defined so as to ensure that the OH group is approaching the π-face of the ring. The other constraints are R < 0.1, no disorder and no errors and normalised H positions.

Firstly, indoles (as above). There are only a few hits, but even so one can see that they all cluster in the top left triangle, where the distance to C2 is always longer than to C3. Indeed, this is the known position for electrophilic substitution of indoles.

s-cis-ester1

The search can be extended by removing the benzo group so as to also include pyrroles. More hits are obtained, and again most of them collect in the top left triangle. The hot spot indicates that the difference in lengths is ~0.3Å in favour of the 3-position, a very similar discrimination to that previously found for benzene groups with an electron donating substituent.

s-cis-ester1

Next, the N atom is replaced by any atom from group 16 of the periodic table (i.e. O, S, etc). The scatter is now in both top left and bottom right triangles, which suggest much weaker discrimination between C2 and C3;  if anything in favour of C2 (often the observed regiospecificities for such compounds).

s-cis-ester1

Finally, pyridines. Only a slight bias towards the C2 position. With pyridines of course, the electrophile in fact first interacts with the nitrogen lone pair in the plane of the molecule, which perturbs the eventual outcome. So this crystallographic method is perhaps a better intrinsic probe than kinetic reactivity.

s-cis-ester1

References

  1. H.S. Rzepa, "Discovering More Chemical Concepts from 3D Chemical Information Searches of Crystal Structure Databases", Journal of Chemical Education, vol. 93, pp. 550-554, 2015. https://doi.org/10.1021/acs.jchemed.5b00346
  2. H. Rzepa, "Search for HO interactions to indoles, pyrroles, furans, and thiophenes", 2016. https://doi.org/10.14469/hpc/665

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

I’ve started so I’ll finish. The ionisation mechanism and kinetic isotope effects for 1,3-dimethylindolin-2 one

Thursday, January 7th, 2016

This is the third and final study deriving from my Ph.D.[1]. The first two topics dealt with the mechanism of heteroaromatic electrophilic attack using either a diazonium cation or a proton as electrophile, followed by either proton abstraction or carbon dioxide loss from the resulting Wheland intermediate. This final study inverts this sequence by starting with the proton abstraction from an indolinone by a base to create/aromatize to a indole-2-enolate intermediate, which only then is followed by electrophilic attack (by iodine).  Here I explore what light quantum chemical modelling might cast on the mechanism.

Indole diazocoupling

The concentration of I3– is used to follow the reaction, given by the expression:  [I3–] = k1[B][indolinone]t – k-1/k2*ln[I3–] + const, where  k2* = k2/715[I–] + k2' , the latter being the rate coefficient for the reaction between the enolate intermediate and I3–. With appropriate least squares analysis of this rate equation,‡ a value for k1 using either 1H or 2H (≡ D) isotopes can be extracted and this gives an isotope effect k1H/k1D of 6.3 ± 0.6. Note that this value does NOT depend on [B]. Here, I am going to try to see if I can construct a quantum mechanical model which reproduces this value.

Indole diazocoupling

  1. Model 1 uses just three water molecules as a proton relay (B3LYP+D3/Def2-TZVP/SCRF=water).
  2. Model 2 uses 2H2O.NaOH solvated by two extra passive water molecules. Since under these conditions, the NaOH is largely ionic, [B] ≡ [OH–]
Model ΔG‡298 (ΔH‡298) kH/kD (298K) DataDOIs
1 28.0 (22.9) 10.3 [2],[3],[4]
2 2.5 (2.8) 4.4 [5],[6],[7]

The plot of rate vs [B] shows[1] that the uncatalysed (water) rate is very slow (intercept passes more or less through zero) and the calculated free energy barrier (28.0 kcal/mol) confirms a slow rate at ambient temperatures. Note in the final (aromatized) product, there is a noticeable hydrogen bond between the 3-carbon and a water molecule (2.14Å). The calculated kinetic isotope effect[8] is substantially larger than observed experimentally for the base catalysed contribution.

Indolineone ionization using 3 water molecules

In the presence of NaOH (standard state = 1 atm = 0.044M), the enthalpy barrier drops very substantially to 2.8 kcal/mol and the free energy to 2.5 kcal/mol. Similar behaviour was noted previously on this blog for the hydrolysis of thalidomide. Although the magnitude of the reduction in barrier in fact implies an extremely fast reaction, recollect that [B]=[OH–] appears in the rate equation  and since its value is very much less than 0.044M, the observed rate is relatively slow.

Indolineone ionization using 3 water molecules + NaOH

The calculated KIE for the hydroxide catalysed mechanism is much smaller that for the water route, but also smaller than is observed. This is a value uncorrected for tunnelling, which given the small barrier might be significant. 

These calculations show how a model for ionization of indolinone can be constructed, and used to e.g. probe the sensitivity of KIE to perturbations induced by ring substituents, which may form the basis of a future post.

‡This is a non-linear equation with kinetics that straddle zero and first order behaviour. In 1972, it was not easily possible to graph such functions in a manner where the slope of a linear plot would yield the rate constant. It was only computers and languages such as Fortran which allowed such non-linear least squares analysis of the rate. In the event, it turned out that the presence of 50% methanol in the mixed aqueous solutions was the cause; in other solvents the kinetics approximated zero order behavour very well.

References

  1. B.C. Challis, and H.S. Rzepa, "Heteroaromatic hydrogen exchange reactions. Part VIII. The ionisation of 1,3-dimethylindolin-2-one", Journal of the Chemical Society, Perkin Transactions 2, pp. 1822, 1975. https://doi.org/10.1039/p29750001822
  2. H.S. Rzepa, "C 10 H 17 N 1 O 4", 2016. https://doi.org/10.14469/ch/191786
  3. H.S. Rzepa, "C 10 H 17 N 1 O 4", 2016. https://doi.org/10.14469/ch/191765
  4. H.S. Rzepa, "C10H17NO4", 2016. https://doi.org/10.14469/ch/191784
  5. H.S. Rzepa, "C 10 H 20 N 1 Na 1 O 6", 2016. https://doi.org/10.14469/ch/191787
  6. H.S. Rzepa, "C 10 H 20 N 1 Na 1 O 6", 2016. https://doi.org/10.14469/ch/191782
  7. H.S. Rzepa, "C10H20NNaO6", 2016. https://doi.org/10.14469/ch/191785
  8. H. Rzepa, "Mechanisms and kinetic isotope effects for the base catalysed ionisation of 1,3-dimethyl indolinone.", 2016. https://doi.org/10.14469/hpc/202

Tags:, , , , , , , , , , , , ,
Posted in Historical, reaction mechanism | No Comments »