Modelling the geometry of unbranched alkanes.

By about C17H36, the geometry of “cold-isolated” unbranched saturated alkenes is supposed not to contain any fully anti-periplanar conformations. [cite]10.1002/anie.201202894[/cite] Indeed, a (co-crystal) of C16H34 shows it to have two-gauche bends.[cite]10.1002/chem.200801428[/cite]. Surprisingly, the longest linear alkane I was able to find a crystal structure for, C28H58 appears to be fully extended[cite]10.1107/S0108768191011059[/cite],[cite]10.1107/S0567740876005025[/cite] (an early report of a low quality structure for C36H74[cite]10.1107/S0365110X5600111X[/cite] also appears to show it as linear). Here I explore how standard DFT theories cope with these structures.

I start with noting the use of a TZVP basis set. In a recent article[cite]10.1039/C3SC53416B[/cite] we noted that the basis-set-superposition-errors for this basis were about a quarter of that for the standard Pople-type 6-311G(d,p) basis that I tend to use for modelling in this blog. This matters, since the relative energy of a folded-conformation vs an extended linear one might depend on the quality of the basis set and its inherent BSSE. The DFT method is the classical B3LYP. I also modelled C58H118 as the hydrocarbon as being well beyond the region anticipated above for folding of the chain to have started (no, there is no crystal structure). The geometries of linear and bent forms are shown below.

003001

The relative free energy of the V-shaped bent form[cite]10.6084/m9.figshare.978501[/cite] emerges as 3.5 kcal/mol higher than the linear form[cite]10.6084/m9.figshare.978502[/cite]. Now, to add a Grimme-D3 dispersion correction to the energies. The V-shape of the bent form now adopts the hairpin mode,[cite]10.6084/m9.figshare.978832[/cite] and its energy is now 2.5 kcal/mol lower than the linear form.[cite]10.6084/m9.figshare.978833[/cite]

002

Note in the above the very slight strange oscillation (kink) that appears about 11 atoms away from the hairpin bend. I repeated this with the wB97XD DFT procedure (in which dispersion is implicit) and found the same result.

As triple-ζ basis quality modelling of molecules with >100 atoms becomes increasingly common, it is worth repeating yet again that the model should always contain dispersion (and solvent if appropriate) corrections as default. Indeed, it is probably also worth re-investigating much early modelling (by this I mean modelling done ten or more years ago) to see if such corrections significantly influence the conclusions.[cite]10.1039/C3SC53416B[/cite]


The searches cannot be carried out according to the formula CnH2n+2, but must be done individually for the value of n. I gave up at C50.

Tags: , , ,

Leave a Reply