Tandem High-Pressure Crystallography-Optical Spectroscopy Unpacks Noncovalent Interactions of Piezochromic Fluorescent Molecular Rotors.

To develop luminescent molecular materials with predictable and stimuli-responsive emission, it is necessary to correlate changes in their geometries, packing structures, and noncovalent interactions with the associated changes in their optical properties. Here, we demonstrate that high-pressure single-crystal X-ray diffraction can be combined with high-pressure UV-visible absorption and fluorescence emission spectroscopies to elucidate how subtle changes in structure influence optical outputs. A piezochromic aggregation-induced emitter, sym-heptaphenylcycloheptatriene (Ph7C7H), displays bathochromic shifts in its absorption and emission spectra at high pressure. Parallel X-ray measurements identify the pressure-induced changes in specific phenyl-phenyl interactions responsible for the piezochromism. Pairs of phenyl rings from neighboring molecules approach the geometry of a stable benzene dimer, while conformational changes alter intramolecular phenyl-phenyl interactions correlated with a relaxed excited state. This tandem crystallographic and spectroscopic analysis provides insights into how subtle structural changes relate to the photophysical properties of Ph7C7H and could be applied to a library of similar compounds to provide general structure-property relationships in fluorescent organic molecules with rotor-like geometries.

Quantifying Mechanical Properties of Molecular Crystals: A Critical Overview of Experimental Elastic Tensors.

This review presents a critical and comprehensive overview of current experimental measurements of complete elastic constant tensors for molecular crystals. For a large fraction of these molecular crystals, detailed comparisons are made with elastic tensors obtained using the corrected small basis set Hartree-Fock method S-HF-3c, and these are shown to be competitive with many of those obtained from more sophisticated density functional theory plus dispersion (DFT-D) approaches. These detailed comparisons between S-HF-3c, experimental and DFT-D computed tensors make use of a novel rotation-invariant spherical harmonic description of the Young's modulus, and identify outliers among sets of independent experimental results. The result is a curated database of experimental elastic tensors for molecular crystals, which we hope will stimulate more extensive use of elastic tensor information-experimental and computational-in studies aimed at correlating mechanical properties of molecular crystals with their underlying crystal structure.

Global Analysis of the Mechanical Properties of Organic Crystals.

Organic crystals, although widely studied, have not been considered nascent candidate materials in engineering design. Here we summarize the mechanical properties of organic crystals that have been reported over the past three decades, and we establish a global mechanical property profile that can be used to predict and identify mechanically robust organic crystals. Being composed of light elements, organic crystals populate a narrow region in the mechanical property-density space between soft, disordered organic materials and stiff, ordered materials. Two subsets of extraordinarily stiff and hard organic crystalline materials were identified and rationalized by the normalized number density, strength, and directionality of their intermolecular interactions. We conclude that future lightweight, soft, all-organic components in devices should capitalize on the greatest asset of organic single crystals-namely, the combination of long-range structural order and softness.

Measurement of Electric Fields Experienced by Urea Guest Molecules in the 18-Crown-6/Urea (1:5) Host-Guest Complex: An Experimental Reference Point for Electric-Field-Assisted Catalysis.

High-resolution synchrotron and neutron single-crystal diffraction data of 18-crown-6/(pentakis)urea measured at 30 K are combined, with the aim of better appreciating the electrostatics associated with intermolecular interactions in condensed matter. With two 18-crown-6 molecules and five different urea molecules in the crystal, this represents the most ambitious combined X-ray/synchrotron and neutron experimental charge density analysis to date on a cocrystal or host-guest system incorporating such a large number of unique molecules. The dipole moments of the five urea guest molecules in the crystal are enhanced considerably compared to values determined for isolated molecules, and 2D maps of the electrostatic potential and electric field show clearly how the urea molecules are oriented with dipole moments aligned along the electric field exerted by their molecular neighbors. Experimental electric fields in the range of 10-19 GV m-1, obtained for the five different urea environments, corroborate independent measurements of electric fields in the active sites of enzymes and provide an important experimental reference point for recent discussions focused on electric-field-assisted catalysis.

Revisiting a Historical Concept by Using Quantum Crystallography: Are Phosphate, Sulfate and Perchlorate Anions Hypervalent?

There are many examples of atoms in molecules that violate Lewis' octet rule, because they have more than four electron pairs assigned to their valence. These atoms are referred to as hypervalent. However, hypervalency may be regarded as an artifact arising from Lewis' description of molecules, which is based on the assumption that electrons are localized in two-center two-electron bonds and lone pairs. In the present paper, the isoelectronic phosphate (PO4 3- ), sulfate (SO4 2- ) and perchlorate (ClO4 - ) anions were examined with respect to the concept of hypervalency. Lewis formulas containing a hypervalent central atom exist for all three anions. Based on X-ray wavefunction refinements of high-resolution X-ray diffraction data of representative crystal structures (MgNH4 PO4 ⋅6 H2 O, Li2 SO4 ⋅H2 O, and KClO4 ), complementary bonding analyses were performed. In this way, experimental information from the new field of quantum crystallography validate long-known facts, or refute long-standing misunderstandings. It is shown that the P-O and S-O bonds are highly polarized covalent bonds and, thus, the increase in the valence population following three-center four-electron bonding is not sufficient to yield hypervalent phosphorus or sulfur atoms, respectively. However, for the highly covalent Cl-O bond, most bonding indicators imply a hypervalent chlorine atom.

Investigation of an Unusual Crystal Habit of Hydrochlorothiazide Reveals Large Polar Enantiopure Domains and a Possible Crystal Nucleation Mechanism.

The observation of an unusual crystal habit in the common diuretic drug hydrochlorothiazide (HCT), and identification of its subtle conformational chirality, has stimulated a detailed investigation of its crystalline forms. Enantiomeric conformers of HCT resolve into an unusual structure of conjoined enantiomorphic twin crystals comprising enantiopure domains of opposite chirality. The purity of the domains and the chiral molecular conformation are confirmed by spatially revolved synchrotron micro-XRD experiments and neutron diffraction, respectively. Macroscopic inversion twin symmetry observed between the crystal wings suggests a pseudoracemic structure that is not a solid solution or a layered crystal structure, but an unusual structural variant of conglomerates and racemic twins. Computed interaction energies for molecular pairs in the racemic and enantiopure polymorphs of HCT, and the observation of large opposing unit-cell dipole moments for the enantiopure domains in these twin crystals, suggest a plausible crystal nucleation mechanism for this unusual crystal habit.

Bridging Crystal Engineering and Drug Discovery by Utilizing Intermolecular Interactions and Molecular Shapes in Crystals.

Most structure-based drug discovery methods utilize crystal structures of receptor proteins. Crystal engineering, on the other hand, utilizes the wealth of chemical information inherent in small-molecule crystal structures in the Cambridge Structural Database (CSD). We show that the interaction surfaces and shapes of molecules in experimentally determined small-molecule crystal structures can serve as effective tools in drug discovery. Our description of the shape and interaction propensities of molecules in their crystal structures can be used to screen them for specific binding compatibility with protein targets, as demonstrated through the high-throughput profiling of around 138 000 small-molecule structures in the CSD and a series of drug-protein crystal structures. Electron-density-based intermolecular boundary surfaces in small-molecule crystal structures and in target-protein pockets are utilized to identify potential ligand molecules from the CSD based on 3D shape and intermolecular interaction matching.

Intermolecular interactions in molecular crystals: what's in a name?

Structure-property relationships are the key to modern crystal engineering, and for molecular crystals this requires both a thorough understanding of intermolecular interactions, and the subsequent use of this to create solids with desired properties. There has been a rapid increase in publications aimed at furthering this understanding, especially the importance of non-canonical interactions such as halogen, chalcogen, pnicogen, and tetrel bonds. Here we show how all of these interactions - and hydrogen bonds - can be readily understood through their common origin in the redistribution of electron density that results from chemical bonding. This redistribution is directly linked to the molecular electrostatic potential, to qualitative concepts such as electrostatic complementarity, and to the calculation of quantitative intermolecular interaction energies. Visualization of these energies, along with their electrostatic and dispersion components, sheds light on the architecture of molecular crystals, in turn providing a link to actual crystal properties.

The Elusive Structural Origin of Plastic Bending in Dimethyl Sulfone Crystals with Quasi-isotropic Crystal Packing.

Bending in molecular crystals is typically associated with the anisotropy of intermolecular interactions. The intriguing observation is reported of plastic bending in dimethyl sulfone, which exhibits nearly isotropic crystal packing and interaction topology, defying the known structural models of bending crystals. The origin of the bending phenomenon has been explored in terms of intermolecular interaction energies, experimental X-ray charge density analysis, and variable temperature neutron diffraction studies. H⋅⋅⋅H dihydrogen interactions and differences in electrostatic complementarity between molecular layers are found to facilitate the bending behavior.

Structural Collapse of the Hydroquinone-Formic Acid Clathrate: A Pressure-Medium-Dependent Phase Transition.

The energy landscape governing a new pressure-induced phase transition in the hydroquinone-formic acid clathrate is reported in which the host structure collapses, opening up the cavity channels within which the guest molecules migrate and order. The reversible isosymmetric phase transition causes significant changes in the morphology and the birefringence of the crystal. The subtle intermolecular interaction energies in the clathrate are quantified at varying pressures using novel model energies and energy frameworks. These calculations show that the high-pressure phase forms a more stable host network at the expense of less-stable host-guest interactions. The phase transition can be kinetically hindered using a nonhydrostatic pressure-transmitting medium, enabling the comparison of intermolecular energies in two polymorphic structures in the same pressure range. Overall this study illustrates a need for accurate intermolecular energies when analyzing self-assembly structures and supramolecular aggregates.

Energy frameworks and a topological analysis of the supramolecular features in in situ cryocrystallized liquids: tuning the weak interaction landscape via fluorination.

Weak intermolecular interactions observed in crystalline materials are often influenced or forced by stronger interactions such as classical hydrogen bonds. Room temperature liquids offer a scenario where such strong interactions are absent so that the role and nature of the weak interactions can be studied more reliably. In this context, we have analyzed the common organic reagent benzoyl chloride (BC) and a series of its fluorinated derivatives using in situ cryocrystallography. The intermolecular interaction energies have been estimated and their topologies explored using energy framework analysis in a series of ten benzoyl chloride analogues, which reveal that the ππ stacking interactions serve as the primary building blocks in these crystal structures. The crystal packing is also stabilized by a variety of interaction motifs involving weak C-HO/F/Cl hydrogen bonds and FF, FCl, and ClCl interactions. It is found that fluorination alters the electrostatic nature of the benzoyl chlorides, with subsequent changes in the formation of different weak interaction motifs. The effects of fluorination on these weak intermolecular interactions have been systematically analyzed further via detailed inputs from a topological analysis of the electron density and Hirshfeld surface analysis.

Host perturbation in a ß-hydroquinone clathrate studied by combined X-ray/neutron charge-density analysis: implications for molecular inclusion in supramolecular entities.

X-ray/neutron (X/N) diffraction data measured at very low temperature (15 K) in conjunction with ab initio theoretical calculations were used to model the crystal charge density (CD) of the host-guest complex of hydroquinone (HQ) and acetonitrile. Due to pseudosymmetry, information about the ordering of the acetonitrile molecules within the HQ cavities is present only in almost extinct, very weak diffraction data, which cannot be measured with sufficient accuracy even by using the brightest X-ray and neutron sources available, and the CD model of the guest molecule was ultimately based on theoretical calculations. On the other hand, the CD of the HQ host structure is well determined by the experimental data. The neutron diffraction data provide hydrogen anisotropic thermal parameters and positions, which are important to obtain a reliable CD for this light-atom-only crystal. Atomic displacement parameters obtained independently from the X-ray and neutron diffraction data show excellent agreement with a |ΔU| value of 0.00058 Å(2) indicating outstanding data quality. The CD and especially the derived electrostatic properties clearly reveal increased polarization of the HQ molecules in the host-guest complex compared with the HQ molecules in the empty HQ apohost crystal structure. It was found that the origin of the increased polarization is inclusion of the acetonitrile molecule, whereas the change in geometry of the HQ host structure following inclusion of the guest has very little effect on the electrostatic potential. The fact that guest inclusion has a profound effect on the electrostatic potential suggests that nonpolarizable force fields may be unsuitable for molecular dynamics simulations of host-guest interaction (e.g., in protein-drug complexes), at least for polar molecules.

Simulations of guest transport in clathrates of Dianin's compound and hydroquinone.

Clathrates have been proposed for use in a variety of applications including gas storage, mixture separation and catalysis due to the potential for controlled guest diffusion through their porous lattices. Here molecular dynamics simulations are employed to study guest transport in clathrates of hydroquinone (HQ) and Dianin's compound (DC). Systems investigated were HQ with methanol and acetonitrile, and DC with methanol and ethanol. Simulations were set up with one guest in the pore, two guests in the pore and one vacancy in the pore and a filled pore, and free-energy barriers for movement between cavities of the pore were estimated for all cases. Comparison between these simulations indicates that guest transport most likely proceeds by molecules moving from full to empty cavities consecutively, one by one, rather than in a concerted manner. Thus, the presence of empty cavities is very important for guest transport, which becomes more energetically demanding in fully loaded systems. Flexibility of the host can assist guest transport. In the studied DC clathrates transport occurs via an intermediate conformation in which the hydroxyl group of the alcohol guest molecule participates in the hydrogen-bonded ring of the host. We also address the issue of the number of methanol guest molecules that DC accommodates, for which conflicting information exists. We found that this is likely to be temperature dependent and suggest that under some conditions the system is most likely non-stoichiometric.

Visualizing lithium-ion migration pathways in battery materials.

The understanding of lithium-ion migration through the bulk crystal structure is crucial in the search for novel battery materials with improved properties for lithium-ion conduction. In this paper, procrystal calculations are introduced as a fast, intuitive way of mapping possible migration pathways, and the method is applied to a broad range of lithium-containing materials, including the well-known battery cathode materials LiCoO2 , LiMn2 O4 , and LiFePO4 . The outcome is compared with both experimental and theoretical studies, as well as the bond valence site energy approach, and the results show that the method is not only a strong, qualitative visualization tool, but also provides a quantitative measure of electron-density thresholds for migration, which are correlated with theoretically obtained activation energies. In the future, the method may be used to guide experimental and theoretical research towards materials with potentially high ionic conductivity, reducing the time spent investigating nonpromising materials with advanced theoretical methods.

Bonded radii and the contraction of the electron density of the oxygen atom by bonded interactions.

The bonded radii for more than 700 bonded pairs of atoms, comprising more than 50 oxide crystals, extracted from experimental and theoretical electron density distributions, are averaged and compared with the ionic radii for first, second, and third row atoms. At odds with the assumption of a "fixed" ionic radius of 1.40 Å for the oxide anion, the bonded radius for the anion, r(b)(O), decreases systematically from 1.40 to 0.65 Å as the electron density distribution of the atom is progressively polarized and contracted by its bonded interactions. The radii for the more electropositive metal atoms agree with the ionic radii when the electron density distribution of the anion is largely unpolarized by its bonded interactions. However, those for the more electronegative metal atoms are progressively larger than the ionic radii as the electron density distribution of the anion is progressively polarized and contracted along the bond vectors with decreasing bond length. The progressive decrease of r(b)(O) indicates that the compilation of sets of ionic radii, based on a fixed radius for the oxide anion, is problematic and impacts the accuracy of the ionic radii for the metal atoms. The assumption of a "fixed" radius for the anion, made in the derivation of sets of radii, not only results in unrealistic negative ionic radii for the more electronegative atoms but also in ionic radii that are as much as 0.5 Å smaller than the bonded radii, particularly for the more electronegative M atoms. The lack of agreement between the ionic and the bonded radii for the more shared bonded interactions is ascribed to the progressive increase in the polarization and contraction of the electron density of the oxide anion by the bonded interactions with a concomitant decrease in the radius of the anion, a factor that was largely neglected in the compilation of the ionic radii for fluoride, oxide, sulfide, and nitride crystals. The close agreement of the bonded radii and procrystal bonded radii is consistent with the argument that the chemical forces that govern the electron density distributions and bonded radii are largely atomic in nature, resulting in comparable electron density distributions.

A transferable quantum mechanical energy model for intermolecular interactions using a single empirical parameter.

The calculation of intermolecular interactions in molecular crystals using model energies provides a unified route to understanding the complex interplay of driving forces in crystallization, elastic properties and more. Presented here is a new single-parameter interaction energy model (CE-1p), extending the previous CrystalExplorer energy model and calibrated using density functional theory (DFT) calculations at the ωB97M-V/def2-QZVP level over 1157 intermolecular interactions from 147 crystal structures. The new model incorporates an improved treatment of dispersion interactions and polarizabilities using the exchange-hole dipole model (XDM), along with the use of effective core potentials (ECPs), facilitating application to molecules containing elements across the periodic table (from H to Rn). This new model is validated against high-level reference data with outstanding performance, comparable to state-of-the-art DFT methods for molecular crystal lattice energies over the X23 set (mean absolute deviation 3.6 kJ mol-1) and for intermolecular interactions in the S66x8 benchmark set (root mean-square deviation 3.3 kJ mol-1). The performance of this model is further examined compared to the GFN2-xTB tight-binding model, providing recommendations for the evaluation of intermolecular interactions in molecular crystal systems.

Application of atomic Hirshfeld surface analysis to intermetallic systems: is Mn in cubic CeMnNi4 a thermoelectric rattler atom?

The Mn atom in the cubic polymorph of CeMnNi(4) appears to be located in an oversized cage-like structure, and anomalously large atomic displacement parameters (ADPs) for the Mn atom indicate that it is a potential "rattler" atom. Here, multitemperature synchrotron powder X-ray diffraction data measured between 110 and 900 K are used to estimate ADPs for the Mn "guest" atom and the "host" structure atoms in cubic CeMnNi(4). The ADPs are subsequently fitted with Debye and Einstein models, giving Θ(D) = 301(2) K for the "host" structure and Θ(E) = 165(2) K for the Mn atom. This is higher than typical Einstein temperatures for rattlers in thermoelectric skutterudites and clathrates (Θ(E) = 50-80 K), indicating that the Mn atom in cubic CeMnNi(4) is more strongly bonded. In order to probe the chemical interactions of the potential Mn rattler atom, atomic Hirshfeld surface (AHS) analysis is carried out and compared with AHS analysis of well-established guest atom rattlers in archetypical skutterudites, MCoSb(3). Surprisingly, the skutterudite rattlers have more deformed AHSs than the Mn atom in cubic CeMnNi(4). This is related to the highly ionic nature of the skutterudite rattlers, which is not taken into account in the neutral spherical atom approach of the AHS. Additionally, visualization of void spaces in the two materials using the procrystal electron density shows that while the Mn atom is tightly fitting in the CeMnNi(4) structure then the La atom in the skutterudite is truly situated in an oversized cage of the host structure. Overall, we conclude that the Mn atom in cubic CeMnNi(4) cannot be coined a rattler.

Computational study of methyl group dynamics in the hydroquinone clathrate of acetonitrile.

We report molecular dynamics simulations of the acetonitrile clathrate of hydroquinone, with a focus on the dynamics of acetonitrile methyl groups. There are three inequivalent acetonitrile molecules in the unit cell, one with its dipole parallel to the c-axis, and the other two antiparallel. Although these three guest molecules have previously been found to exhibit two slightly different frequencies of rotation over a wide range of temperatures, the frequencies could not be assigned to specific methyl groups. Perhaps counterintuitively, our simulations suggest that the molecule with the lower frequency is one of the two molecules oriented the same way, the different dynamical behaviour being due to subtle differences in the environments of the molecules.

Solution-phase decomposition of ferrocene into wüstite-iron oxide core-shell nanoparticles.

We report an improved method for the controlled solvent-phase decomposition of ferrocene into highly crystalline monodisperse iron oxide nanoparticles at relatively low temperatures. Solution-phase decomposition of ferrocene into nanoparticles has received little attention in the literature, due to the percieved stability of ferrocene. However, we synthesised wüstite FeO-iron oxide core-shell nanoparticles by thermally decomposing ferrocene in 1-octadecene solvent and in the presence of oleic acid and oleylamine, as surfactants. We report procedures that provide cubic and spherical core-shell iron oxide nanoparticles whose size (29.3 ± 2.3 nm for spheres, 38.6 ± 6.9 nm for distorted cubes and 23.5 ± 2.4 nm for distorted cubes with concave faces) and shape can be controlled through simple adjustments to reaction parameters. Transmission electron microscopy, scanning transmission electron microscopy, energy dispersive X-ray spectroscopy, electron energy-loss spectroscopy and powder X-ray diffraction analysis methods were used to characterise the nanoparticles.

Mechanism of concerted hydrogen bond reorientation in clathrates of Dianin's compound and hydroquinone.

Molecular dynamics provides a means to examine the mechanism of reorientation of hydrogen bond networks that are present in a range of biological and crystalline materials. Simulations of hydroxyl reorientation in the six-membered hydrogen bonded rings in crystalline clathrates of Dianin's compound (DC) and hydroquinone (HQ) reveal that in the clathrate of Dianin's compound with ethanol (DC:ethanol), hydroxyl groups perform single independent flips, and occasionally all six hydroxyls in a ring reorient following a sequential mechanism with participation of the guest ethanol molecule. The free energy estimated for this process agrees well with experimental results. The simulations suggest that hydroxyl reorientation occurs in the empty DC lattice as well, but at a higher energy cost, from which we conclude that it is the participation of ethanol that lowers the barrier of reorientation. Single independent flips of hydroxyl groups are observed to be more frequent in the hydroquinone clathrate with methanol (HQ:methanol) than in DC:ethanol, but reorientation of all six hydroxyls does not occur. This is attributed to the larger difference in energy between the original and reoriented positions of hydroxyl hydrogen atoms in HQ:methanol compared to DC:ethanol.