Posts Tagged ‘calculated free energy barrier’

Pyrophoric metals + the mechanism of thermal decomposition of magnesium oxalate.

Sunday, March 19th, 2017

A pyrophoric metal is one that burns spontaneously in oxygen; I came across this phenomenon as a teenager doing experiments at home. Pyrophoric iron for example is prepared by heating anhydrous iron (II) oxalate in a sealed test tube (i.e. to 600° or higher). When the tube is broken open and the contents released, a shower of sparks forms. Not all metals do this; early group metals such as calcium undergo a different reaction releasing carbon monoxide and forming calcium carbonate and not the metal itself. Here as a prelude to the pyrophoric reaction proper, I take a look at this alternative mechanism using calculations.


There are ~60 crystal structures of metal oxalates, of which several are naturally occurring minerals (Fe, humboldtine[1], Ca, Weddellite[2], Li[3], Na[4], K[5], Cs[6]. The natural geometry of the oxalate di-anion is planar (torsion 0 or 180°) but a small number are twisted such as the caesium oxalate.

The kinetics of pyrolysis of a number of metal  oxalates were studied some years ago (Ca[7], Li[8]) indicating barriers ranging from 53-68 kcal/mol. One proposed mechanism is as shown in this article.[7]

It was surmised from the kinetic analysis that the k1 activation step (rotation about the C-C bond from planar to twisted) was ~12 ± 20 kcal/mol, whilst steps k2 or k3 had the much higher activation energy noted above. A search (of Scifinder) for quantum mechanical “reality checks” of this mechanism revealed a blank and so I apply such a check here using Mg as the metal.

The carbonyl extrusion step (ωB97XD/Def2-TZVPPD/SCRF=water, DOI: 10.14469/hpc/2320) was studied with a water solvent field applied in an effort to mimic the solid state crystal structure of the species as a better representation of the ionic lattice than a pure vacuum calculation.An IRC (intrinsic reaction coordinate, DOI: 10.14469/hpc/2324) reveals the start-point geometry still has a very small negative force constant (-38 cm-1, DOI: 10.14469/hpc/2321) which now corresponds to a small rotation about the C-C bond to give a C2-symmetric conformation:

But the barrier for this process is tiny and nothing like the ~12 ± 20 kcal/mol inferred from the kinetic analysis. Perhaps most of the incentive to pack into a totally planar geometry comes from the interactions in the ionic lattice. The calculated free energy barrier (ΔG298‡ 54.7 kcal/mol, ΔG755‡ 55.1 kcal/mol) is within the reported measured range.

The mechanism for production of pyrophoric metal itself is likely to be far more complex, involving (inter alia) electron transfer from oxygen to metal. If I find anything I will report back here.

References

  1. T. Echigo, and M. Kimata, "Single-crystal X-ray diffraction and spectroscopic studies on humboldtine and lindbergite: weak Jahn–Teller effect of Fe2+ ion", Physics and Chemistry of Minerals, vol. 35, pp. 467-475, 2008. https://doi.org/10.1007/s00269-008-0241-7
  2. C. Sterling, "Crystal structure analysis of weddellite, CaC2O4.(2+x)H2O", Acta Crystallographica, vol. 18, pp. 917-921, 1965. https://doi.org/10.1107/s0365110x65002219
  3. https://doi.org/
  4. G.A. Jeffrey, and G.S. Parry, "The Crystal Structure of Sodium Oxalate", Journal of the American Chemical Society, vol. 76, pp. 5283-5286, 1954. https://doi.org/10.1021/ja01650a007
  5. Dinnebier, R.E.., Vensky, S.., Panthofer, M.., and Jansen, M.., "CCDC 192180: Experimental Crystal Structure Determination", 2003. https://doi.org/10.5517/cc6fzcy
  6. Dinnebier, R.E.., Vensky, S.., Panthofer, M.., and Jansen, M.., "CCDC 192182: Experimental Crystal Structure Determination", 2003. https://doi.org/10.5517/cc6fzf0
  7. F.E. Freeberg, K.O. Hartman, I.C. Hisatsune, and J.M. Schempf, "Kinetics of calcium oxalate pyrolysis", The Journal of Physical Chemistry, vol. 71, pp. 397-402, 1967. https://doi.org/10.1021/j100861a029
  8. D. Dollimore, and D. Tinsley, "The thermal decomposition of oxalates. Part XII. The thermal decomposition of lithium oxalate", Journal of the Chemical Society A: Inorganic, Physical, Theoretical, pp. 3043, 1971. https://doi.org/10.1039/j19710003043

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

Ménage à deux: Non-classical SC bonds.

Wednesday, December 30th, 2009

A previous post posed the question; during the transformation of one molecule to another, what is the maximum number of electron pairs that can simultaneously move either to or from any one atom-pair bond as part of the reaction? A rather artificial example (atom-swapping between three nitrosonium cations) was used to illustrate the concept, in which three electron pairs would all move from a triple bond to a region not previously containing any electrons to form new triple bonds and destroy the old. Here is a slightly more realistic example of the phenomenon, illustrated by the (narcisistic) reaction below of a bis(sulfur trifluoride) carbene. Close relatives of this molecule are actually known, with either one SF3 of the units replaced by a CF3 group or a SF5 replacing the SF3 (DOI: 10.1021/ja00290a038 ).

F3SCSF3 and the nature of its S-C bonds

The two C-S bonds in this molecule are not the same (and similarly for the CF3 analogue), one being long (single), the other short (assumed triple), and the angle subtended at the central carbon is around 150° (B3LYP/cc-pVTZ calculation, DOI: 10042/to-3643). The transition state for interconverting one form to the other would presumably correspond to the concerted movement of two pairs of electrons from one CS region to the other as shown above, not so much a Ménage à trois, as a Ménage à deux! The transition state itself (DOI: 10042/to-3644) has C2 symmetry, with a calculated free energy barrier of 31 kcal/mol and ν 284i cm-1 for the bond shifting process.

Transition state for bond equalisation

Transition state for bond equalisation. Click for animation

The molecule above does have a further point of interest; one of the sulfur atoms (the triply bonded one) is approximately tetrahedral in coordination, whilst the other has a “T-shape”. An inorganic chemist would describe one sulfur as tetravalent (oxidation state IV), the other as hexavalent (oxidation state VI) and the equilibrium between them a dismutation of the two oxidation states. Does this have any reality? The ELF method has been mentioned a number of times in these posts, and it is applied here to seek an answer. The ELF basin centroids are shown below.

The ELF function, as isosurfaces contoured at various thresholds

ELF basins for F3SCSF3. Click for 3D

The integrations are as follows: 14 = 2.24 (a single C-S bond), 30=1.66 (an incipient carbene forming, as implied above), 13+15+16 = 4.34 (a reasonably persuasive triple bond, comprising, unusually, three separated basins). The fluorines 2, 3 and 6 all exhibit bonding basins to the S (respectively 2.17, 2.17 and 2.09), but fluorines 1,5 and 4 do not! Sulfur 8 additionally has a lone pair, 29=2.31, but sulfur 9 does not. One aspect of this analysis is the nature of the triple bond between S9-C7. Because the three basins are separate, does that mean that the bond cannot rotate about its axis?

AIM Analysis of F3SCSF3

An alternative AIM analysis is shown above. Now, the CS triple bond is reduced to a single bond critical point (BCP), labelled 10. AIM allows a property known as bond ellipticity to be computed at that BCP. Typically, single and triple bonds have ellipticities close to zero, whilst double bonds have a value of around 0.4 to 0.5. That for point 10 is 0.18, which seems to support the ELF analysis above. Pretty unsual bonding it would have to be agreed!

ELF centroids for transition state for dismutation.

But what of the original question posed at the start in the diagram; do two pairs of electrons move away together from one triple bond to form another. A further ELF analysis at the transition state for this process reveals that in effect the two pairs do different things. One localizes onto the carbon, to form a proper carbene, the other becomes a sulfur lone pair. So the valence dismutation involves three pairs of electrons, not two as shown at the start, with each pair doing its own thing.

Six-electron model for valence isomerism in F3SCSF3

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