- Research article
- Open Access
A mechanophysical phase transition provides a dramatic example of colour polymorphism: the tribochromism of a substituted tri(methylene)tetrahydrofuran-2-one
© Asiri et al.; licensee springer. 2014
- Received: 28 July 2014
- Accepted: 30 October 2014
- Published: 16 December 2014
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.
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.
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.
- Crystal structure analysis
- Molecular and lattice energy calculations
- Colour polymorphism
One compound, identified as I, was obtained as yellow crystals on crystallization from a deep red solution of 3:7 ethyl acetate and petroleum. Since the yellow form gave deep red solutions when re-dissolved, a photochromic change in solution on exposure to ultraviolet light could not be established. With a view to determine whether the yellow form was photochromic in the solid state, a crystal was crushed on a filter paper with the intention then of irradiating the expected yellow powder with ultraviolet light. By simply crushing the crystal, a dark red powder was produced which did not revert to the yellow form unless the red powder was re-crystallised from 3:7 ethyl acetate and petroleum. This seemingly irreversible solid state process of colour change by rubbing, grinding or fracturing the crystals, was termed by us as “tribochromism” from the Greek ‘tribos’ - rubbing and ‘chroma’ - colour. In order to explore this phenomenon further, additional compounds were synthesised with different substituents at the isopropylidene site, (see Scheme 2), but these exhibited variable properties. Thus, when the isopropylidene was replaced by a cyclopentylidene, the product II was obtained directly as deep red crystals. No yellow crystals could be obtained for this compound. A product which could be obtained as both yellow and red crystals was eventually obtained, when the methyl groups of compound I were replaced by cyclopropyl groups, giving compound III. In order to explore the relationship between structure and colour, we set out to determine the crystal structures of key products. These studies showed that the yellow (IIIa) and red (IIIb) forms from the solution of compound III are formally polymorphs of III, in which the molecules have different conformations. The structure of a fourth compound, IV, a very close analogue of III, but which gave only yellow crystals, was also determined, to check its molecular conformation, and also with a view to see how the slight difference in molecular structure would affect the crystal packing. The syntheses and structure determinations of IIIa and IIIb have been briefly described in a previous, short communication . For this current paper, we have now had the opportunity to examine the relationship between the structures using computational chemistry techniques to ascertain the electronic structures of the molecules in the two forms, a and b, of compound III, and calculate lattice and related energies of these and the other, related crystal structures.
a). Syntheses and crystallizations
All solvents were carefully dried before use. Petroleum refers to the fraction b.p. 60-80°C. The preparations of the first stage succinic anhydrides, and all characterisation data are described in the supplementary information, Additional file 1.
Compound I. 5-dicyanomethylene-4-diphenylmethylene-3-isopropylidenetetrahydrofuran-2-one
Diethylamine (0.48 g, 6.6 mmoles) was added to a mixture of diphenylmethylene(isopropylidene) succinic anhydride (1 g, 3.3 mmoles) and malononitrile (0.22 g, 3.3 mmoles) in tetrahydrofuran (15 ml) at 0°C and then stirred at room temperature (20 h). Solvent was removed and the residue treated with ether (5 ml). A colourless salt separated which was cyclised by stirring with acetyl chloride (1 ml) in dichloromethane (1 ml) for 6 hours. The residue was purified by chromatography on a column of silica gel (100 g) using a 3:7 mixture of ethyl acetate and petroleum as eluant. The dicyanomethylene compound was obtained as yellow crystals on evaporation of the deep red solution, m.p. 182-183°C (0.64 g, 55%).
Compound II. 5-dicyanomethylene-4-diphenylmethylene-3-cyclopentanylidentetrahydrofuran-2-one
Diethylamine (0.67 g, 9.2 mmoles) was added to a mixture of cyclopentanylidene (diphenylmethylene)succinic anhydride (1.5 g, 4.5 mmoles) and malononitrile (0.3 g, 4.5 mmoles) in tetrahydrofuran (15 ml) at 0°C and then stirred at room temperature (20 h). Solvent was removed and the residue treated with ether (5 ml). A colourless salt separated which was cyclised by stirring with acetyl chloride (1 ml) in dichloromethane (1 ml) for 6 hours. The residue was purified by chromatography on a column of silica gel (100 g) using a 3:7 mixture of ethyl acetate and petroleum as eluant. The dicyanomethylene compound was obtained as dark red crystals on partial evaporation of the solvent, m.p. 162-164°C (0.64 g, 55%) and was recrystallised from dichloromethane and petroleum.
Compound IIIa,b. 5-dicyanomethylene-4-diphenylmethylene-3-dicyclopropylmethylenetetrahydrofuran-2-one
Diethylamine (0.53 g, 7 mmoles) was added to a mixture of diphenylmethylene (isopropylidene) succinic anhydride (1.26 g, 3.5 moles) and malononitrile (0.24 g, 3.6 mmoles) in tetrahydrofuran (10 ml) at 0°C and then stirred at room temperature (12 h). Solvent was removed. On treatment of the residue with ether (5 ml), a colourless salt separated, which was cyclised as described above. The residue was purified by chromatography on a column of silica gel (100 g) using a 3:7 mixture of ethyl acetate and petroleum as eluant. On slow evaporation of the deep red eluant, yellow crystals of form IIIa of the dicyanomethylene compound (0.3 g, 21% yield) m.p. 183-184°C were initially obtained.
Further evaporation of the concentrated deep red solution gave dark red crystals of form IIIb of the same dicyanomethylene compound (0.1 g 7% yield) m.p. 188-189°C.
Compound IV(c,d). E- and Z- 5-dicyanomethylene-3-dicyclopropylmethylene-4-[4-methylphenyl)(phenyl)methylene]tetrahydrofuran-2-one
Diethylamine (0.59 g, 8.1 moles) was added to a mixture of 1:1 E: Z mixture of 4-methylphenyl (phenyl)methylene(dicyclopropylidene)succinic anhydride (1.50 g, 4.05 mmoles) and malononitrile (0.26 g, 4.05 mmoles) in tetrahydrofuran (10 ml) at 0°C and then stirred at room temperature (20 h). The solvent was removed, the residue dissolved in DCM (10.ml) and treated with acetyl chloride (10 ml) to give the dicyanomethylene compound as yellow crystals only. (M.p. 183–185°C). The single crystal selected for X-ray crystallography turned out to be the Z-isomer, IV(d).
b). X-ray crystallography
Selected crystal structure data
CCDC Dep/Ref code
c). Computational studies
i) Lattice energy calculations
Solid state calculations of the lattice energies and the crystal structures of I, II, IIIa, IIIb and IV were carried out using the GRACE package . GRACE provides an efficient algorithm for optimising the co-ordinates and lattice parameters of molecular crystals. Lattices energies and gradients are provided by a solid state density functional theory method with corrections for dispersive interactions (the DFT-D method) . GRACE uses the VASP  program with the PW91 density functional  to calculate the lattice energy and its gradients. The dispersive correction is provided by a damped molecular mechanical potential. Starting with the experimental crystal structures, the unit cells and molecular geometries were fully optimised within the constraints of the experimental space group symmetries.
The DFT-D lattice energy is defined by the energy of the process of going from molecules in the ‘gas phase’ to the crystal, with the lowest energy conformer of each molecule taken as the reference energy. It is possible to break down the contributions to the lattice energy into terms arising from conformational, conformational deformation, van der Waals and Coulombic components. In the following description of the solid state calculations a ‘molecular’ or ‘gas phase’ calculation refers to a DFT-D calculation of a single molecule in a unit cell which is large enough that a molecule does not interact with images of itself. Typically unit cells with dimensions of over 23 Å were needed to ensure convergence of the molecular calculations. In these calculations the molecular structure was optimised but the unit cell and the orientation of the molecule in the unit cell was held fixed. A ‘crystal calculation’ refers to a DFT-D optimisation of the crystal structure in which the crystal structure, comprising the molecular geometry, the molecular orientation and the unit cell dimensions were fully optimised. In order to estimate the Coulombic contribution to the DFT-D lattice energy, atomic charges which reproduce the molecular electrostatic potential associated with the optimised molecule were calculated using GRACE. These point charges were then used to calculate the change in Coulomb energy on going from the molecule in the ‘gas-phase’ to the solid state. The van der Waals' contribution was estimated from the difference in the molecular and solid state dispersive energy corrections used by the DFT-D method. The conformational contribution arises from the difference in energies between molecularly optimised conformations of the same molecule, whilst the conformational deformation energy arises from the change in energy of a conformation on moving from the gas phase to the solid state. The conformational deformation energy was calculated from molecular calculations of the energy difference between a molecule optimised in the solid state and in the gas phase. The sum of the conformational deformation and conformational energy terms represents the increase in internal energy of the molecule due to its presence in the crystal. It follows that the ‘deformation’ term and the conformational energy term are always positive energies. The remaining contributions to lattice energy, including terms arising from polarisation and induction, are calculated from the differences between the total lattice energies and the sum of all the other contributions; conformational, conformational deformation, van der Waals and Coulomb.
(ii) Molecular calculations
Molecular calculations were carried out with the ORCA package . Geometry optimisations were performed starting with molecular structures extracted from the experimental crystal structures. The structures were optimised using the TZVPP basis  with the BLYP density functional, including a dispersive correction . The basis set consists of 5 s, 3p, 2d and 1f functions on C, N and O, and 3 s, 2p and 1d functions on H. Calculations of the electronic excitation energies for the lowest excited singlet states for each molecule were performed at the resulting optimised geometries. The time dependent density functional method using the B2PLYP functional  was used to calculate the excited state energies and wavefunctions. The B2PLYP functional is a double hybrid functional which is a combination of the generalized gradient approximations for the correlation and exchange, the Hartree-Fock exchange and a second order correlation term which is included by perturbation theory . It has been shown to provide reasonable predictions of the excited state energies of organic molecules ,. A TZVP basis was used  for the calculation of the excited state energies. The basis set consists of 5 s, 3p and 1d functions on C, N and O, and 3 s and 1p functions on H and is the same as the basis set used in the study of Goerigk et al. . To estimate the effect of the crystal environment on the absorption spectrum, molecular geometry optimisations and calculations of the first excite states were repeated using the COSMO solvation model  using an effective medium with a refractive index of 75 and a dielectric constant of 3.0 which is representative of an organic molecular crystalline environment .
Excited state calculations at the geometry of the optimised ground state give vertical excitation energies. In order to estimate the adiabatic excitation energies, which are more appropriate for comparison with experiment, calculations of the optimised ground and excited states were performed using the TZVP basis with the BLYP functional. For these calculations no dispersive correction was employed. For each molecule an adiabatic correction was determined by calculating the energy change on optimising the first excited state starting at the optimised ground state geometry.
d) Molecular packing studies
Molecular packing analysis
(i) Crystal structure calculations
Comparison of DFT-D optimised and experimental unit cells
Unit cell parameters
Breakdown of contributions to the total DFT-D lattice energies
Conformational energy (kcal mol −1)
Conformational deformation energy (kcal mol −1)
van der Waals energy (kcal mol −1)
Coulomb energy (kcal mol −1)
Remaining energy (kcal mol −1)
Lattice energy (kcal mol −1)
In the case of the two polymorphs of III, form IIIb is calculated to be 4.5 kcal mol−1 more stable than IIIa. The origin of the stability difference is complex. Although the molecular conformation found in IIIa is 1.6 kcal mol−1 higher in energy than that in IIIb, this is compensated for, in part, by the van der Waals' contribution which favours IIIa by 1.2 kcal mol−1. The stability of IIIb relative to IIIa seems to come from an increased Coulombic contribution (2.1 kcal mol−1) along with a further 2.1 kcal mol−1 related to other contributions such as induction effects.
For all crystals the conformational deformation contribution to the lattice energy is quite small (less than 1.4 kcal mol−1). The van der Waals' contribution is the largest contribution (between 77 and 89% of the total lattice energy) and generally gets larger as the molecules get larger. The Coulombic contribution varies between 20 and 25% of the total and the contribution whose origin has not been attributed varies between 4 and 7% of the total.
The conformational energy of the molecule in structure IIIa is marginally higher, and the conformation thus less stable, than that of the molecule in structure IIIb
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, however, by the contribution of a more stabilising Coulomb energy to the overall more favourable lattice energy of structure IIIb.
The break-down final lattice energies proposed here is not unique. A method based on gas phase optimisation of a conformer starting from its experimental structure as found in the crystal has been used by Bernstein and Cruz-Cabeza . The Pixel method , which uses the gas phase electron density of the molecule to calculate the interaction energies in the experimental or calculated crystal structure has also been used.
(ii) Molecular calculations
Molecular energies, HUMO-LUMO gap and the predicted wavelengths associated with the first excited states
Absolute energies a(Hartree)
Relative energies akcal/mol
HOMO LUMO aGap (eV)
Gas phase 1 stExcited state b(eV)
COSMO 1 stExcited state b(eV)
Adiabatic Correction c(eV)
For all molecules the wavefunction of the first excited state is dominated by a single electronic configuration, which arises from the transfer of an electron from the HOMO to the LUMO. It is not surprising therefore that there is a good correlation between the calculated HOMO-LUMO energy gap and the energy of the first excited state. The partial inclusion of the influence of the crystal environment through the COSMO model introduces a reduction in the calculated energy of the first excited stated by between 0.20 and 0.47 eV. Attempts were made to calculate the optimised geometries of the excited states in order to estimate the adiabatic excitation energies. In the cases of molecules I, IIIa and IV, optimised geometries were determined. For molecules II and IIIb, both of which have twisted conformations, geometry optimisations failed because large changes in the molecular geometries were causing excited states to cross during the optimisation processes.
(iii). Molecular geometry comparisons
Bond Lengths. In the 5-membered ring, the O(1)-C(2) bond is consistently shorter than the C(1)-O(1) bond, indicating a strengthened interaction between O(1) and the dicyanomethylene carbon C(2). This is likely to be a feature of the fundamental electronic structure of the trimethylene substituted tetrahydrofuran-2-one, and particularly the dicyanomethylene group, since this bonding distribution is also found in the structurally related molecules LORMAH , MUZBOY  and YETBEE , each of which contains this component. The calculated, minimised structures IIIa(C) and IIIb(C) both capture this result also.
A second feature, again captured by the calculations, is the small, but consistent lengthening – 0.03-0.04 Å, of the diphenylmethylene double bond, C(3) = C(8), in the twisted, red forms, compared with the folded forms.
Bond Angles. Here we examine the angles involving the tetrahydrofuran-2-one ring and its organic substituents. A first, simple calculation of the sum of the internal angles of the ring shows that the twisted structures of II and IIIb retain quite planar rings (sum of internal angles are 539.5° and 539.6°, whilst the folded structures have rings with internal angles summing to 533.1° and 535.2°, reflecting the formal reduction of ring angles at a fold. An additional point of interest is the comparison of the external angles at each of the methylene substituents, where distortions may occur due to differences in intramolecular interactions with neighbouring groups. Thus, we find that for the dicyanomethylene function, the O(1)-C(2) = C(5) and C(3)-C(2) = C(5) angles show a difference of 11-20°, with the larger angles on the side which clashes with the diphenylmethylene group. For this latter group itself, the angles differ only by 1-4°, whilst the methylene group at C(4) shows an intermediate distortion of 8-12°, again with the larger angle due to clashes with the diphenyl functions.
Dihedral and Torsion Angles. Clearly, the critical parameters for the major distortions in the folded, yellow molecules are the angles of fold across C(2)…C(4), which have values of 23.7(3)o, 20.9(10)o, 23.7(2)o for I, IIIa and IV, respectively. Whilst it logical to quantify the twist about the C(3) = C(8) bonds in II and IIIb by the dihedral angles between the two planes C(2),C(3),C(4) and C(9),C(8),C(15), which have values of 37.4(2)o and 35.7(3)o for structures II and IIIb, these do not give the full story, and it is interesting to compare the important torsion angles also. Values of these angles will indicate whether the distortions in the bond angles mentioned above are essentially in-plane or whether any unbalanced intramolecular repulsions again occur, which may “distort the distortions”. The torsion angles we are interested are the cis bonds at the ends of each of the carbon-carbon double bonds, typified by the pair O(1)-C(2) = C(5)-C(6) and C(3)-C(2) = C(5)-C(7). The relevant data is presented in the last section of Table 2, and shows that the differences vary from 0.5° to 17°. There is no marked difference related to the folded or twisted conformations, but, in general, the larger values are obtained for the structures containing the bis-cyclopropyl substituent, indicating the effect of the accommodation of the two bulky ligands in these molecules.
Selected, comparative molecular geometry parameters for the five structures; X denotes X-ray results, C the results from the computational study
C(2) = C(5)
C(3) = C(8)
C(4) = C(21)
O(1)-C(2) = C(5)
C(3)-C(2) = C(5)
C(2)-C(3) = C(8)
C(4)-C(3) = C(8)
C(3)-C(4) = C(21)
C(1)-C(4) = C(21)
O(1)-C(2) = C(5)-C(6)
C(3)-C(2) = C(5)-C(7)
C(2)-C(3) = C(8)-C(9)
C(4)-C(3) = C(8)-C(15)
C(3)-C(4) = C(21)-C(23)
C(1)-C(4) = C(21)-C(22)
The theoretical calculations of the parameters defining the crystal and molecular structures, map very well with those found experimentally. The molecular and lattice energy calculations indicate that the folded yellow conformation of compound III has a higher molecular energy than the twisted red conformation, but that the crystal packing energy is lower, and supports the less stable conformation in structure IIIb and presumably in structure II, for which no red crystals could be obtained. However, we presume that the twisted form is predominant in all solutions, which are red. These results will, no doubt, be of particular interest to colleagues engaged in structure prediction calculations.
We thank Avant-garde Materials Simulation for the provision of a courtesy licence to the GRACE software package, and EPSRC for support of the Crystallography facilities.
- Fan M-G, Yu L, Wei Z: Fulgide Family Compounds: Synthesis, Photochromism and Applications in Photochromic and Thermochromic Compounds. 1999, Plenum Press, New YorkGoogle Scholar
- Arunasalam VC, Baxter I, Drake SR, Hursthouse MB, Malik KMA, Otway DJ: Synthesis and characterization of the first group 2 mixed β-diketonate alkoxide complexes - X-ray crystal-structures of [Ca4(tmhd)4(OEt)4(EtOH)4], [Ca4(tmhd)6(OCH2CH2NMe2)2], and [H2Ba4(tmhd)6(OCH2CH2OPri)4]. Inorg Chem. 1995, 34: 5295-5306. 10.1021/ic00125a032.View ArticleGoogle Scholar
- Avante-garde Materials Simulation: GRACE (Generation, Ranking, and Characterization Engine). Deutschland GmbH, Merzhauser Str. 177, 79100 Freiburg, Germany.Google Scholar
- Neumann MA, Perrin MA: Energy ranking of molecular crystals using density functional theory calculations and an empirical van der Waals correction. J Phys Chem B. 2005, 109: 15531-15541. 10.1021/jp050121r.View ArticleGoogle Scholar
- Kresse G, Hafner J: Ab initio molecular-dynamics for liquid-metals. Phys Rev B. 1993, 47: 558-561. 10.1103/PhysRevB.47.558. b. Kresse G, Furthmuller J: Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Physical Review B 1996, 54:11169–11186.,View ArticleGoogle Scholar
- Wang Y, Perdew JP: Spin scaling of the electron-gas correlation-energy in the high-density limit. Phys Rev B. 1991, 43: 8911-8916. 10.1103/PhysRevB.43.8911.View ArticleGoogle Scholar
- Neese F: The ORCA program system. Wiley Interdiscip Rev. 2012, 2: 73-78.Google Scholar
- Schafer A, Horn H, Ahlrichs R: Fully optimized contracted Gaussian-basis sets for atoms Li to Kr. J Chem Phys. 1992, 97: 2571-2577. 10.1063/1.463096. b. Weigend F, Ahlrichs R: Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Physical Chemistry Chemical Physics 2005, 7:3297–3305.,View ArticleGoogle Scholar
- Grimme S, Antony J, Ehrlich S, Krieg H: A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J Chem Phys. 2010, 132: 19-10.1063/1.3382344.View ArticleGoogle Scholar
- Schwabe T, Grimme S: Double-hybrid density functionals with long-range dispersion corrections: higher accuracy and extended applicability. Phys Chem Chem Phys. 2007, 9: 3397-3406. 10.1039/b704725h.View ArticleGoogle Scholar
- Grimme S, Neese F: Double-hybrid density functional theory for excited electronic states of molecules. J Chem Phys. 2007, 127: 18-10.1063/1.2772854.View ArticleGoogle Scholar
- Goerigk L, Moellmann J, Grimme S: Computation of accurate excitation energies for large organic molecules with double-hybrid density functionals. Phys Chem Chem Phys. 2009, 11: 4611-4620. 10.1039/b902315a.View ArticleGoogle Scholar
- Sinnecker S, Rajendran A, Klamt A, Diedenhofen M, Neese F: Calculation of solvent shifts on electronic g-tensors with the conductor-like screening model (COSMO) and its self-consistent generalization to real solvents (Direct COSMO-RS). J Phys Chem A. 2006, 110: 2235-2245. 10.1021/jp056016z.View ArticleGoogle Scholar
- Cooper TG, Hejczyk KE, Jones W, Day GM: Molecular polarization effects on the relative energies of the real and putative crystal structures of valine. J Chem Theory Comput. 2008, 4: 1795-1805. 10.1021/ct800195g.View ArticleGoogle Scholar
- Gelbrich, T, Hursthouse, MB: A versatile procedure for the identification, description and quantification of structural similarity in molecular crystals. Crystengcomm 2005, 7; 324–336.,Google Scholar
- Asiri AMA, Heller HG, Hursthouse MB, Karalulov A: Tribochromic compounds, exemplified by 3-dicyclopropylmethylene-5-dicyanomethylene-4-diphenylmethylenetetrahydrofuran-2-one. Chemical Communications 2000:799–800,Google Scholar
- Macrae CF, Bruno IJ, Chisholm JA, Edgington PR, McCabe P, Pidcock E, Rodriguez-Monge L, Taylor R, van de Streek J, Wood PA: Mercury CSD 2.0 - New features for the visualization and investigation of crystal structures . J. Appl. Cryst. 2008, 41: 466-470. 10.1107/S0021889807067908.View ArticleGoogle Scholar
- Csikos E, Ferenczy GG, Angyan JG, Bocskei Z, Simon K, Gonczi C, Hermecz I: Colour polymorphism of a bis(quinoxaline) compound. Eur J Organ Chem 1999:2119–2125,Google Scholar
- Imai Y, Kinuta T, Nagasaki K, Harada T, Sato T, Tajima N, Sasaki Y, Kuroda R, Matsubara Y: Conformational and color polymorphism of achiral 2-methyl-3-(2-naphthalenylthio)-1,4-naphthalenedione. Crystengcomm. 2009, 11: 1223-1226. 10.1039/b900806c.View ArticleGoogle Scholar
- Braun DE, Gelbrich T, Jetti RKR, Kahlenberg V, Price SL, Griesser UJ: Colored polymorphs: thermochemical and structural features of N-picryl-p-toluidine polymorphs and solvates. Crystal Growth Design. 2008, 8: 1977-1989. 10.1021/cg8000224.View ArticleGoogle Scholar
- Cruz-Cabeza AJ, Bernstein J: Conformational polymorphism. Chem. Rev. 2014, 114: 2170-91. 10.1021/cr400249d.View ArticleGoogle Scholar
- Chan HCS, Kendrick J, Leusen FJJ: Molecule VI, a benchmark crystal-structure-prediction sulfonimide: are its polymorphs predictable?. Angew Chem Int Ed. 2011, 50: 2979-2981. 10.1002/anie.201007488.View ArticleGoogle Scholar
- Kendrick J, Leusen FJJ, Neumann MA, van de Streek J: Progress in crystal structure prediction. Chem A Eur J. 2011, 17: 10736-10744. 10.1002/chem.201100689.View ArticleGoogle Scholar
- Montis R, Hursthouse MB, Chan HCS, Kendrick J, Leusen FJJ: Experimental and theoretical investigations of the polymorphism of 5-chloroacetoxybenzoic acid (5-chloroaspirin). Crystengcomm. 2012, 14: 1672-1680. 10.1039/c2ce06313a.View ArticleGoogle Scholar
- Kendrick J, Stephenson GA, Neumann MA, Leusen FJJ: Crystal structure prediction of a flexible molecule of pharmaceutical interest with unusual polymorphic behavior. Crystal Growth Design. 2013, 13: 581-589. 10.1021/cg301222m.View ArticleGoogle Scholar
- Chan HCS, Kendrick J, Neumann MA, Leusen FJJ: Towards ab initio screening of co-crystal formation through lattice energy calculations and crystal structure prediction of nicotinamide, isonicotinamide, picolinamide and paracetamol multi-component crystals. Crystengcomm. 2013, 15: 3799-3807. 10.1039/c3ce40107c.View ArticleGoogle Scholar
- Kendrick J, Montis R, Hursthouse MB, Leusen FJJ: “In-silico seeding”: Isostructurality and pseudoisostructurality in a family of aspirin derivatives. Crystal Growth & Design. 2013, 13: 2906-2915. 10.1021/cg4003324.View ArticleGoogle Scholar
- Gavezzotti A: Calculation of lattice energies of organic crystals: the PIXEL integration method in comparison with more traditional methods. Zeitschrift Fur Krist. 2005, 220: 499-510. 10.1524/zkri.220.5.499.65063.Google Scholar
- Asiri AM, Ng SW: (E, E)-2-3,4-bis(4-methylbenzylidene)-5-oxotetrahydrofuran-2-ylidene propanedinitrile . Acta Crystallogr Sect E Struct Rep Online. 2009, 65: 0761-U2055. 10.1107/S1600536809008526.View ArticleGoogle Scholar
- Liang YC, Dvornikov AS, Rentzepis PM: New near infrared-sensitive photochromic fluorescing molecules. J Mater Chem. 2003, 13: 286-290. 10.1039/b206577k.View ArticleGoogle Scholar
- Heller HG, Hughes DS, Hursthouse MB, Koh KVS: A new class of photochromic compounds exemplified by E-5-dicyanomethylene-4-(dialkyl and dicycloalkyl)methylene 1-(2,5-methyl-3-furyl) and (2-methyl-5-phenyl-3-thienyl)ethylidene tetrahydrofuran-2-ones. Journal of the Chemical Society-Chemical Communications 1994:2713–2714,Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.