Towards Computational Screening for New Energetic Molecules: Calculation of Heat of Formation and Determination of Bond Strengths by Local Mode Analysis.

The reliable determination of gas-phase and solid-state heats of formation are important considerations in energetic materials research. Herein, the ability of PM7 to calculate the gas-phase heats of formation for CNHO-only and inorganic compounds has been critically evaluated, and for the former, comparisons drawn with isodesmic equations and atom equivalence methods. Routes to obtain solid-state heats of formation for a range of single-component molecular solids, salts, and co-crystals were also evaluated. Finally, local vibrational mode analysis has been used to calculate bond length/force constant curves for seven different chemical bonds occurring in CHNO-containing molecules, which allow for rapid identification of the weakest bond, opening up great potential to rationalise decomposition pathways. Both metrics are important tools in rationalising the design of new energetic materials through computational screening processes.

Inversion twinning in a second polymorph of the hydrochloride salt of the recreational drug ethylone.

A second polymorph of the hydrochloride salt of the recreational drug ethylone, C12H16NO3(+)·Cl(-), is reported [systematic name: (±)-2-ethylammonio-1-(3,4-methylenedioxyphenyl)propane-1-one chloride]. This polymorph, denoted form (A), appears in crystallizations performed above 308 K. The originally reported form (B) [Wood et al. (2015). Acta Cryst. C71, 32-38] crystallizes preferentially at room temperature. The conformations of the cations in the two forms differ by a 180° rotation about the C-C bond linking the side chain to the aromatic ring. Hydrogen bonding links the cations and anions in both forms into similar extended chains in which any one chain contains only a single enantiomer of the chiral cation, but the packing of the ions is different. In form (A), the aromatic rings of adjacent chains interleave, but pack equally well if neighbouring chains contain the same or opposite enantiomorph of the cation. The consequence of this is then near perfect inversion twinning in the structure. In form (B), neighbouring chains are always inverted, leading to a centrosymmetric space group. The question as to why the polymorphs crystallize at slightly different temperatures has been examined by density functional theory (DFT) and lattice energy calculations and a consideration of packing compactness. The free energy (ΔG) of the crystal lattice for polymorph (A) lies some 52 kJ mol(-1) above that of polymorph (B).

A mechanophysical phase transition provides a dramatic example of colour polymorphism: the tribochromism of a substituted tri(methylene)tetrahydrofuran-2-one.

BACKGROUND: Derivatives of fulgides have been shown to have interesting photochromic properties. We have synthesised a number of such derivatives and have found, in some cases, that crystals can be made to change colour on crushing, a phenomenon we have termed "tribochromism". We have studied a number of derivatives by X-ray crystallography, to see if the colour is linked to molecular structure or crystal packing, or both, and our structural results have been supported by calculation of molecular and lattice energies. RESULTS: A number of 5-dicyanomethylene-4-diphenylmethylene-3-disubstitutedmethylene-tetrahydrofuran-2-one compounds have been prepared and structurally characterised. The compounds are obtained as yellow or dark red crystals, or, in one case, both. In two cases where yellow crystals were obtained, we found that crushing the crystals gave a deep red powder. Structure determinations, including those of the one compound which gave both coloured forms, depending on crystallisation conditions, showed that the yellow crystals contained molecules in which the structure comprised a folded conformation at the diphenylmethylene site, whilst the red crystals contained molecules in a twisted conformation at this site. Lattice energy and molecular conformation energies were calculated for all molecules, and showed that the conformational energy of the molecule in structure IIIa (yellow) is marginally higher, and the conformation thus less stable, than that of the molecule in structure IIIb (red). However, the van der Waals energy for crystal structure IIIa, is slightly stronger than that of structure IIIb - which may be viewed as a hint of a metastable packing preference for IIIa, overcome by the contribution of a more stabilising Coulomb energy to the overall more favourable lattice energy of structure IIIb. CONCLUSIONS: Our studies have shown that the crystal colour is correlated with one of two molecular conformations which are different in energy, but that the less stable conformation can be stabilised by its host crystal lattice. Graphical abstractGraphical representation of the structural and colour change in the tribochromic compound (III).