Hydrogen-bond-modulated negative linear compressibility in a V-shaped molecular crystal

A material with the "hidden" negative linear compressibility (NLC) will expand along a specific crystal direction upon uniformly compression to a critical pressure; such materials are thought to be promising candidates for non-linear actuators, switches and sensors. Herein, we use density functional theory (DFT) calculations to uncover the hidden NLC in a V-shaped molecular crystal, bis(5-amino-1,2,4-triazol-3-yl)methane (BATZM). The calculations indicate that the crystal is normally compressed over the pressure range of 0-3 GPa while it expands along the b-axis when the external hydrostatic pressure exceeds 3 GPa. The compressive behavior of the BATZM crystal is modulated by inter-molecular hydrogen bonds, which act as highly compressible springs at low pressures but robust struts at high pressures. Hence, the crystal prefers to compress the hydrogen bonds coupled with PLC at first and flatten the molecules, coupled with later NLC to resist the increasing external pressure. The compressive behavior of BATZM provides a strategy to design more hidden NLC materials via the rational use of the hydrogen bonds.

ROY Crystallization on Poly(ethylene) Fibers, a Model for Bed Net Crystallography

Many long-lasting insecticidal bed nets for protection against disease vectors consist of poly(ethylene) fibers in which insecticide is incorporated during manufacture. Insecticide molecules diffuse from within the supersaturated polymers to surfaces where they become bioavailable to insects and often crystallize, a process known as blooming. Recent studies revealed that contact insecticides can be highly polymorphic. Moreover, insecticidal activity is polymorph-dependent, with forms having a higher crystal free energy yielding faster insect knockdown and mortality. Consequently, the crystallographic characterization of insecticide crystals that form on fibers is critical to understanding net function and improving net performance. Structural characterization of insecticide crystals on bed net fiber surfaces, let alone their polymorphs, has been elusive owing to the minute size of the crystals, however. Using the highly polymorphous compound ROY (5-methyl-2-[(2-nitrophenyl)-amino]thiophene-3-carbonitrile) as a proxy for insecticide crystallization, we investigated blooming and crystal formation on the surface of extruded poly(ethylene) fibers containing ROY. The blooming rates, tracked from the time of extrusion, were determined by UV-vis spectroscopy after successive washes. Six crystalline polymorphs (of the 13 known) were observed on poly(ethylene) fiber surfaces, and they were identified and characterized by Raman microscopy, scanning electron microscopy, and 3D electron diffraction. These observations reveal that the crystallization and phase behavior of polymorphs forming on poly(ethylene) fibers is complex and dynamic. The characterization of blooming and microcrystals underscores the importance of bed net crystallography for the optimization of bed net performance.

Frontiers of molecular crystal structure prediction for pharmaceuticals and functional organic materials

The reliability of organic molecular crystal structure prediction has improved tremendously in recent years. Crystal structure predictions for small, mostly rigid molecules are quickly becoming routine. Structure predictions for larger, highly flexible molecules are more challenging, but their crystal structures can also now be predicted with increasing rates of success. These advances are ushering in a new era where crystal structure prediction drives the experimental discovery of new solid forms. After briefly discussing the computational methods that enable successful crystal structure prediction, this perspective presents case studies from the literature that demonstrate how state-of-the-art crystal structure prediction can transform how scientists approach problems involving the organic solid state. Applications to pharmaceuticals, porous organic materials, photomechanical crystals, organic semi-conductors, and nuclear magnetic resonance crystallography are included. Finally, efforts to improve our understanding of which predicted crystal structures can actually be produced experimentally and other outstanding challenges are discussed.

Polymorphs of Substituted p-Toluenesulfonanilide: Synthesis, Single-Crystal Analysis, Hirshfeld Surface Exploration, and Theoretical Investigation

Polymorphism is an exciting feature of chemical systems where a compound can exist in different crystal forms. The present investigation is focused on the two polymorphic forms, triclinic (MSBT) and monoclinic (MSBM), of ethyl 3-iodo-4-((4-methylphenyl)sulfonamido)benzoate prepared from ethyl 4-amino-3-iodobenzoate. The prepared polymorphs were unambiguously confirmed by single-crystal X-ray diffraction (SC-XRD) analysis. According to the SC-XRD results, the molecular configurations of both structures are stabilized by intramolecular N-H···I and C-H···O bonding. The crystal packing of MSBT is different as compared to the crystal packing of MSBM because MSBT is crystallized in the triclinic crystal system with the space group P1̅, whereas MSBM is crystallized in the monoclinic crystal system with the space group P21/c. The molecules of MSBT are interlinked in the form of dimers through N-H···O bonding to form R22(8) loops, while the MSBM molecules are connected with each other in the form of an infinite chain through C-H···O bonding. The crystal packing of both compounds is further stabilized by off-set π···π stacking interactions between phenyl rings, which is found stronger in MSBM as compared to in MSBT. Moreover, Hirshfeld surface exploration of the polymorphs was carried out, and the results were compared with the closely related literature structure. Accordingly, the supramolecular assembly of these polymorphs is mainly stabilized by noncovalent interactions or intermolecular interactions. Furthermore, a density functional theory (DFT) study was also carried out, which provided good support for the SC-XRD and Hirshfeld studies, suggesting the formation of both intramolecular and intermolecular interactions for both compounds.

Behavior of Occupied and Void Space in Molecular Crystal Structures at High Pressure

We report a Monte Carlo algorithm for calculation of occupied ("network") and unoccupied ("void") space in crystal structures. The variation of the volumes of the voids and the network of intermolecular contacts with pressure sensitively reveals discontinuities associated with first- and second-order phase transitions, providing insights into the effect of compression (and, in principle, other external stimuli) at a level between those observed in individual contact distances and the overall unit cell dimensions. The method is shown to be especially useful for the correlation of high-pressure crystallographic and spectroscopic data, illustrated for naphthalene, where a phase transition previously detected by vibrational spectroscopy, and debated in the literature for over 80 years, has been revealed unambiguously in crystallographic data for the first time. Premonitory behavior before a phase transition and crystal collapse at the end of a compression series has also been detected. The network and void volumes for 129 high-pressure studies taken from the Cambridge Structural Database (CSD) were fitted to equation of state to show that networks typically have bulk moduli between 40 and 150 GPa, while those of voids fall into a much smaller range, 2-5 GPa. These figures are shown to reproduce the narrow range of overall bulk moduli of molecular solids (ca. 5-20 GPa). The program, called CellVol, has been written in Python using the CSD Python API and can be run through the command line or through the Cambridge Crystallographic Data Centre's Mercury interface.

Use of a miniature diamond-anvil cell in a joint X-ray and neutron high-pressure study on copper sulfate pentahydrate

Single-crystal X-ray and neutron diffraction data are usually collected using separate samples. This is a disadvantage when the sample is studied at high pressure because it is very difficult to achieve exactly the same pressure in two separate experiments, especially if the neutron data are collected using Laue methods where precise absolute values of the unit-cell dimensions cannot be measured to check how close the pressures are. In this study, diffraction data have been collected under the same conditions on the same sample of copper(II) sulfate pentahydrate, using a conventional laboratory diffractometer and source for the X-ray measurements and the Koala single-crystal Laue diffractometer at the ANSTO facility for the neutron measurements. The sample, of dimensions 0.40 × 0.22 × 0.20 mm3 and held at a pressure of 0.71 GPa, was contained in a miniature Merrill-Bassett diamond-anvil cell. The highly penetrating diffracted neutron beams passing through the metal body of the miniature cell as well as through the diamonds yielded data suitable for structure refinement, and compensated for the low completeness of the X-ray measurements, which was only 24% on account of the triclinic symmetry of the sample and the shading of reciprocal space by the cell. The two data-sets were combined in a single 'XN' structure refinement in which all atoms, including H atoms, were refined with anisotropic displacement parameters. The precision of the structural parameters was improved by a factor of up to 50% in the XN refinement compared with refinements using the X-ray or neutron data separately.

How many more polymorphs of ROY remain undiscovered

With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z′ = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs.

Crystal structure prediction suggests that the low-energy polymorphs of ROY have already been found, but a new high-pressure form is predicted.

Direct evidence for distinct colour origins in ROY polymorphs

ROY is one of the most well-studied families of crystal structures owing to it being the most polymorphic organic material on record. The various red, orange, and yellow colours of its crystal structures are widely-believed to originate from molecular conformation, though the orange needle (ON) polymorph is thought to be an exception. We report high-pressure, single-crystal X-ray measurements which provide direct experimental evidence that the colour origin in ON is intermolecular, revealing that the molecule undergoes minimal deformation but still exhibits a pronounced, reversible, pale orange → dark red colour change between ambient pressure and 4.18 GPa. Our experimental data are rationalised with band structures, calculated using an accurate hybrid DFT approach, where we are able to account for the variation in colour for five polymorphs of ROY. We highlight the outlier behaviour of ON which shows marked π⋯π stacking interactions that are directly modified through application of pressure. Band structure calculations confirm these intermolecular interactions as the origin of the colour change.

Prolific Polymorph Generator ROY in Its Liquid and Glass: Two Conformational Populations Mirroring the Crystalline-State Distribution

5-Methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, dubbed ROY for its numerous crystal polymorphs of red, orange, and yellow colors, has been studied in its liquid and glassy state by infrared spectroscopy. Two populations of conformers are observed, whose equilibrium is characterized by ΔH = 2.4 kJ/mol and ΔS = 8.0 J/K/mol. The two populations correspond to the global and local minima of the torsional energy surface and to the conformational preference of the 13 crystal polymorphs. The local minimum features a more coplanar arrangement of the two aromatic rings, greater π conjugation, and lower CN stretch frequency. In the gas phase, the lowest-energy path between the two minima has an energy barrier 3.9 kJ/mol above the global minimum, consistent with the rapid equilibration between the two populations. The relevance of our result for understanding the prolific polymorphism of ROY is discussed.

Accurate and Efficient Computation of Optical Absorption Spectra of Molecular Crystals: The Case of the Polymorphs of ROY

When calculating the optical absorption spectra of molecular crystals from first principles, the influence of the crystalline environment on the excitations is of significant importance. For such systems, however, methods to describe the excitations accurately can be computationally prohibitive due to the relatively large system sizes involved. In this work, we demonstrate a method that allows optical absorption spectra to be computed both efficiently and at high accuracy. Our approach is based on the spectral warping method successfully applied to molecules in solvent. It involves calculating the absorption spectrum of a supercell of the full molecular crystal using semi-local time-dependent density functional theory (TDDFT), before warping the spectrum using a transformation derived from smaller-scale semi-local and hybrid TDDFT calculations on isolated dimers. We demonstrate the power of this method on three polymorphs of the well-known color polymorphic compound ROY and find that it outperforms both small-scale hybrid TDDFT dimer calculations and large-scale semi-local TDDFT supercell calculations, when compared to the experiment.

Linear negative thermal expansion in Pd(acac)2

We propose that chemical interactions in the coordination compound Pd(acac)2 can be mapped on a wine-rack network, explaining linear negative thermal expansion in a classic coordination compound based on concepts known from network materials.

Theoretical investigation of polymorph- and coformer-dependent photoluminescence in molecular crystals

A novel density-functional approach provides accurate predictions for the colour zoning of ROY polymorphs and the fluorescence energies of a family of 9-acetylanthracene cocrystals.

High-pressure polymorphism in pyridine

Single crystals of the high-pressure phases II and III of pyridine have been obtained by in situ crystallization at 1.09 and 1.69 GPa, revealing the crystal structure of phase III for the first time using X-ray diffraction. Phase II crystallizes in P212121 with Z' = 1 and phase III in P41212 with Z' = ½. Neutron powder diffraction experiments using pyridine-d5 establish approximate equations of state of both phases. The space group and unit-cell dimensions of phase III are similar to the structures of other simple compounds with C 2v molecular symmetry, and the phase becomes stable at high pressure because it is topologically close-packed, resulting in a lower molar volume than the topologically body-centred cubic phase II. Phases II and III have been observed previously by Raman spectroscopy, but have been mis-identified or inconsistently named. Raman spectra collected on the same samples as used in the X-ray experiments establish the vibrational characteristics of both phases unambiguously. The pyridine molecules interact in both phases through CH⋯π and CH⋯N interactions. The nature of individual contacts is preserved through the phase transition between phases III and II, which occurs on decompression. A combination of rigid-body symmetry mode analysis and density functional theory calculations enables the soft vibrational lattice mode which governs the transformation to be identified.

Pressure-induced inclusion of neon in the crystal structure of a molecular Cu2(pacman) complex at 4.67 GPa

Crystals of Cu₂(pacman) inflate on taking up neon at 46 000 atm through a switch in the ligand conformation.

The Effect of Pressure on Halogen Bonding in 4-Iodobenzonitrile

The crystal structure of 4-iodobenzonitrile, which is monoclinic (space group I2/a) under ambient conditions, contains chains of molecules linked through C≡N···I halogen-bonds. The chains interact through CH···I, CH···N and π-stacking contacts. The crystal structure remains in the same phase up to 5.0 GPa, the b axis compressing by 3.3%, and the a and c axes by 12.3 and 10.9 %. Since the chains are exactly aligned with the crystallographic b axis these data characterise the compressibility of the I···N interaction relative to the inter-chain interactions, and indicate that the halogen bond is the most robust intermolecular interaction in the structure, shortening from 3.168(4) at ambient pressure to 2.840(1) Å at 5.0 GPa. The π∙∙∙π contacts are most sensitive to pressure, and in one case the perpendicular stacking distance shortens from 3.6420(8) to 3.139(4) Å. Packing energy calculations (PIXEL) indicate that the π∙∙∙π interactions have been distorted into a destabilising region of their potentials at 5.0 GPa. The structure undergoes a transition to a triclinic ( P 1 ¯ ) phase at 5.5 GPa. Over the course of the transition, the initially colourless and transparent crystal darkens on account of formation of microscopic cracks. The resistance drops by 10% and the optical transmittance drops by almost two orders of magnitude. The I···N bond increases in length to 2.928(10) Å and become less linear [<C-I∙∙∙N = 166.2(5)°]; the energy stabilises by 2.5 kJ mol^-1 and the mixed C-I/I..N stretching frequency observed by Raman spectroscopy increases from 249 to 252 cm^-1. The driving force of the transition is shown to be relief of strain built-up in the π∙∙∙π interactions rather than minimisation of the molar volume. The triclinic phase persists up to 8.1 GPa.

Pressure-induced isosymmetric phase transition in biurea

An isosymmetric phase transition of biurea has been found in which the molecules appear to “snap” from one conformer to another upon compression.

Living in the salt-cocrystal continuum: indecisive organic complexes with thermochromic behaviour

The thermochromic behaviour of the haloaniline – 3,5-dinotrobenzoic acid cocrystals over the temperature range from 30–108 °C.

Structural behaviour of OP-ROY at extreme conditions

A highly-flexible organic molecule, with numerous polymorphs at ambient conditions, shows remarkable metastability, resisting the effects of extreme pressure and temperature.

Role of halogen substituents in a series of polymorphic 2,5-diamino-3,6-dicyanopyrazine derivatives with highly flexible groups

The crystal structures of the ortho-X-benzyl derivatives, where X=F, Cl, Br, I, and Me, of 2,5-bis(N,N-dibenzylamino)-3,6-dicyanopyrazine dyes (C34H24N6X4) were analysed to evaluate the effect of a systematic series of structures on the occurrence of polymorphism. Detailed crystal structure analysis indicated that the thermally stable forms of the polymorphic derivatives (Cl and Br derivatives) were close-packed, whereas those of the non-polymorphic derivatives (F and I derivatives) were stabilised by an intermolecular interaction involving the ortho-substituents. In the thermally metastable forms of the polymorphic derivative, halogen-halogen and halogen-nitrogen interactions contributed to the stabilisation of these crystals in the same way as the thermally stable form of the non-polymorphic derivatives. This indicated that the ease of polymorph occurrence would require an appropriate balance between the crystal energy of the close-packed structure and that of the crystal structure generated mainly by the electrostatic interactions involving the halogens in these halogenated pyrazine derivatives. In addition, the similar tendency of the occurrence of polymorphs in these halogenated pyrazine derivatives was found in 19 sets of halogenated compounds having known crystal structures of F, Cl, Br and I derivatives including at least one polymorphic derivative in the crystal structure database.

Negative Linear Compressibility Due to Layer Sliding in a Layered Metal-Organic Framework

Negative linear compressibility (NLC) is a rare and counterintuitive phenomenon because materials with this property would expand along one specific direction when uniformly compressed. NLC materials have a broad range of potential applications in designing pressure sensors, artificial muscles, and so on. Designing and searching for systems with NLC is desired and crucial for material and compression science. Herein, with the help of high-pressure X-ray diffraction measurements and density functional theory calculations, we find that the 2D layered Co(SCN)₂(pyrazine)₂ exhibits NLC with a new mechanism: layer sliding. When compressed, the ab planes slide along the a axis, leading to the decrease of lattice parameter β, which results in the NLC effect along principal axis X₃ (≈ -0.84a - 0.55c). The layer sliding mechanism opens exciting opportunities for seeking, designing, and synthesizing new classes of materials with anomalous mechanical properties in monoclinic layered or other related systems.

Piezochromism and hydrochromism through electron transfer: new stories for viologen materials

A pyridinium-carboxylate compound undergoes reversible color change under pressure owing to the formation of radicals via electron transfer; dehydration and hydration can also trigger electron transfer.

While viologen derivatives have long been known for electrochromism and photochromism, here we demonstrated that a viologen-carboxylate zwitterionic molecule in the crystalline state exhibits piezochromic and hydrochromic behaviors. The yellow crystal undergoes a reversible color change to red under high pressure, to green after decompression, and finally back to yellow upon standing at ambient pressure. Ultraviolet-visible spectroscopy, X-ray photoelectron spectroscopy, electron paramagnetic resonance X-ray diffraction and DFT calculations suggested that the piezochromism is due to the formation of radicals via pressure-induced electron transfer from carboxylate to pyridinium, without a crystallographic phase transition. It was proposed that electron transfer is induced by pressure-forced reduction of intermolecular donor–acceptor contacts. The electron transfer can also be induced by dehydration, which gives a stable green anhydrous radical phase. The color change is reversible upon reabsorption of water, which triggers reverse electron transfer. The compound not only demonstrates new chromic phenomena for viologen compounds, but also represents the first example of organic mechanochromism and hydrochromism associated with radical formation via electron transfer.

Pharmaceutical polymorph control in a drug-mimetic supramolecular gel

A supramolecular gel designed to chemically mimic the structure of a pharmaceutical compound controls the polymorphic outcome of the crystallization of the substrate.

We report the synthesis of a bis(urea) gelator designed to specifically mimic the chemical structure of the highly polymorphic drug substance ROY. Crystallization of ROY from toluene gels of this gelator results in the formation of the metastable red form instead of the thermodynamic yellow polymorph. In contrast, all other gels and solution control experiments give the yellow form. Conformational and crystal structure prediction methods have been used to propose the structure of the gel and show that the templation of the red form by the targeted gel results from conformational matching of the gelator to the ROY substrate coupled with overgrowth of ROY onto the local periodic structure of the gel fibres.