Properties Of Molecular Crystals Research Articles

MolPackL: Quantification and Interpretation of Intermolecular Interactions Driven by Molecular Packing.

In organic optoelectronic devices, the properties of the aggregated organic materials depend not only on individual molecules or monomers but also significantly on their packing modes. Different from their inorganic counterparts linked by explicit covalent bonds, organic solids exhibit intricate and numerous intermolecular interactions (IMIs). Due to the intrinsic complexity and disorder of IMIs, identifying and understanding them is a formidable challenge in experimental, theoretical, and data-driven approaches. In this work, we constructed an innovative algorithm framework, Molecular Packing Learning (MolPackL), which can accurately quantify elusive IMIs using contact density histograms (CDHs) and efficiently extract intermolecular features for further property prediction of organic solids. It performs satisfactorily in training predictive models of IMI-related properties in molecular crystals. Particularly, the band gap predictive model based on MolPackL achieved the best-reported performance, with an MAE of 0.20 eV and an impressive R2 of 0.92. Class activation mapping (CAM) visually demonstrates MolPackL's accurate identification of effective interaction sites as the molecular packing changes. What is more, the elemental importance analysis verified that the superior score benefits from MolPackL's ability to comprehensively consider multiple influencing factors of IMIs. In summary, MolPackL provides a new framework for quantitative assessment and understanding of the effect of IMIs. The development of MolPackL marks a significant advancement in establishing predictive models of molecular aggregates, deepening the comprehension of IMIs on the material properties. Given the superior performance, we believe that MolPackL will also become a powerful tool in the design of high-performance organic optoelectronic materials.

Surface Engineering of the Mechanical Properties of Molecular Crystals via an Atomistic Model for Computing the Facet Stress Response of Solids.

Invited for the cover of this issue are Mubarak Almehairbi, Vikram C. Joshi, Changquan Calvin Sun and Sharmarke Mohamed. The image depicts the digital exploration of the mechanical properties of crystals on specific facets that may be of interest for materials applications by "dialing-in" their stress response. Read the full text of the article at 10.1002/chem.202400779.

Front Cover: Surface Engineering of the Mechanical Properties of Molecular Crystals via an Atomistic Model for Computing the Facet Stress Response of Solids (Chem. Eur. J. 38/2024)

Front Cover: Surface Engineering of the Mechanical Properties of Molecular Crystals via an Atomistic Model for Computing the Facet Stress Response of Solids (Chem. Eur. J. 38/2024)

Surface Engineering of the Mechanical Properties of Molecular Crystals via an Atomistic Model for Computing the Facet Stress Response of Solids.

Dynamic molecular crystals are an emerging class of crystalline materials that can respond to mechanical stress by dissipating internal strain in a number of ways. Given the serendipitous nature of the discovery of such crystals, progress in the field requires advances in computational methods for the accurate and high-throughput computation of the nanomechanical properties of crystals on specific facets which are exposed to mechanical stress. Here, we develop and apply a new atomistic model for computing the surface elastic moduli of crystals on any set of facets of interest using dispersion-corrected density functional theory (DFT-D) methods. The model was benchmarked against a total of 24 reported nanoindentation measurements from a diverse set of molecular crystals and was found to be generally reliable. Using only the experimental crystal structure of the dietary supplement, L-aspartic acid, the model was subsequently applied under blind test conditions, to correctly predict the growth morphology, facet and nanomechanical properties of L-aspartic acid to within the accuracy of the measured elastic stiffness of the crystal, 24.53±0.56 GPa. This work paves the way for the computational design and experimental realization of other functional molecular crystals with tailor-made mechanical properties.

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory.

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe-Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

Are Elastic Properties of Molecular Crystals within Reach of Density Functional Theory? Accuracy, Robustness, and Reproducibility of Current Approaches

Are Elastic Properties of Molecular Crystals within Reach of Density Functional Theory? Accuracy, Robustness, and Reproducibility of Current Approaches

Bifurcated Polymorphic Transition and Thermochromic Fluorescence of a Molecular Crystal Involving Three-Dimensional Supramolecular Gear Rotation.

The ability of molecules to move and rearrange in the solid state accounts for the polymorphic transition and stimuli-responsive properties of molecular crystals. However, how the crystal structure determines the molecular motion ability remains poorly understood. Here, we report that a three-dimensional (3D) supramolecular gear network in the green-emissive polymorph 1G of a dialkylamino-substituted anthracene-pentiptycene π-system (1) enables an unusual bifurcated polymorphic transition into a yellow-emissive polymorph (1Y) and a new green-emissive polymorph (1G*) via 3D correlated supramolecular rotation. The 90° forward correlated rotation causes the molecular conformation between the octyl and the anthracene units to change from syn to anti, the ladder-like supramolecular columns to constrict, and the gear network to disengage. This cooperative molecular motion is marked by the gradual formation of an intermediate state (1I) across the entire crystal from 170 to 230 °C, which then undergoes bifurcated (forward or backward rotation) and irreversible transitions to form polymorphs 1Y and 1G* at 230-235 °C. Notably, 1G* is similar to 1G but lacks gear engagement, preventing its transformation into 1Y. Nevertheless, 1G can be restored by grinding 1Y or 1G* or fuming with dichloromethane (DCM) vapor. This work illustrates the correlation between the crystal structure and solid-state molecular motion behavior and demonstrates how a 3D molecular gear system efficiently transmits thermal energy to drive the polymorphic transition and induce fluorochromism through significant conformational and packing changes.

Alkyl modification to optimize the ferroelastic properties of molecular crystals

Using alkyl modification, we obtained a ferroelastic material with three phase transitions, orange fluorescence and high-temperature ferroelasticity.

Predicting crystal form stability under real-world conditions

The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2–4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid–solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1–2 kJ mol−1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.

Low-Frequency Vibrations of Saccharides Using Terahertz Time-Domain Spectroscopy and Ab-Initio Simulations

The vibrational properties of molecular crystals in the terahertz range are controlled by the interplay of intermolecular and intramolecular interactions. The resulting delocalized normal modes span a substantial number of atoms within the unit cell of the molecular crystals and cannot be easily attributed to a specific stretching or bending mode. Several compounds were studied to understand the nature of normal modes. However, a systematic analysis of saccharides is still missing. This study investigates the terahertz vibrational properties of various saccharides, including glucose, galactose, lactose, cellobiose, and cellulose. The terahertz spectra were measured using terahertz time-domain spectroscopy. The samples were carefully characterized, and the residual Fabry–Perot oscillations were removed to obtain the absorption coefficient and refractive index of the saccharides. Density functional theory simulations were used to obtain theoretical terahertz spectra, considering hydrogen bonding interactions with an enhanced version of the van der Waals non-local density functional. The results revealed differences in low-energy vibrational modes frequencies, influenced by molecular structure properties, hydrogen bonding networks, and water content. Librations and internal vibrations were identified as dominant dynamics in the saccharides, with significant mixing between intermolecular and intramolecular vibrations. This comprehensive analysis sheds light on the vibrational behavior of saccharides in the terahertz range.

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals.

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew-Burke-Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (∼50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

Predicting the structures and vibrational spectra of molecular crystals containing large molecules with the generalized energy-based fragmentation approach.

The generalized energy-based fragmentation (GEBF) approach under periodic boundary conditions (PBCs) has been developed to facilitate calculations of molecular crystals containing large molecules. The PBC-GEBF approach can help predict structures and properties of molecular crystals at different theory levels by performing molecular quantum chemistry calculations on a series of non-periodic subsystems constructed from the studied systems. A more rigorous formula of the forces on translational vectors of molecular crystals was proposed and implemented, enabling more reliable predictions of crystal structures. Our benchmark results on several typical molecular crystals show that the PBC-GEBF approach could reproduce the forces on atoms and the translational vectors and the optimized crystal structures from the corresponding conventional periodic methods. The improved PBC-GEBF approach is then applied to predict the crystal structures and vibrational spectra of two molecular crystals containing large molecules. The PBC-GEBF approach can provide a satisfactory description on the crystal structure of a molecular crystal containing 312 atoms in a unit cell at density-fitting second-order Møller-Plesset perturbation theory and density functional theory (DFT) levels and the infrared vibrational spectra of another molecular crystal containing 864 atoms in a unit cell at the DFT level. The PBC-GEBF approach is expected to be a promising theoretical tool for electronic structure calculations on molecular crystals containing large molecules.

Pressure induced modification of the electronic properties of stilbene by two-photon spectroscopy.

Carbon nanothreads are the most exciting carbon based nanomaterials recently discovered. Obtained by compressing aromatics around 20GPa, they are characterized by potentially exceptional mechanical properties. The reaction mechanisms have been partly elucidated through computational studies and x-ray diffraction experiments. However, in all these studies, the electronic modifications to which the molecule is subjected with increasing pressure are neglected as also if, and to which extent, the electronic excited states are involved in the high-pressure reactivity. In fact, the pressure increase induces remarkable changes in the electronic properties of molecular crystals, which are often directly related to the reaction's onset and path. We report the pressure evolution of the two-photon induced emission spectrum of crystalline stilbene, the archetype of a class of molecules from which double-core nanothreads are obtained, with the twofold purpose of gaining insight into the reaction mechanism and monitoring if the structural changes observed in x-ray diffraction studies have a detectable counterpart in the electronic properties of the system. The freezing of the spectral diffusion observed on rising pressure is ascribed to a hampered conformational rearrangement because of the larger stiffness of the local environment. The transition to the high pressure phase where the nanothreads form is revealed by the slope change of the pressure shift of all spectral components, while the progressive intensification with pressure of the 0-0 transition suggests a strengthening of the ethylenic bond favoring the charge delocalization on the benzene moieties, which is likely the trigger of the chemical instability.

Mechanical properties and peculiarities of molecular crystals.

In the last century, molecular crystals functioned predominantly as a means for determining the molecular structures via X-ray diffraction, albeit as the century came to a close the response of molecular crystals to electric, magnetic, and light fields revealed that the physical properties of molecular crystals were as rich as the diversity of molecules themselves. In this century, the mechanical properties of molecular crystals have continued to enhance our understanding of the colligative responses of weakly bound molecules to internal frustration and applied forces. Here, the authors review the main themes of research that have developed in recent decades, prefaced by an overview of the particular considerations that distinguish molecular crystals from traditional materials such as metals and ceramics. Many molecular crystals will deform themselves as they grow under some conditions. Whether they respond to intrinsic stress or external forces or interactions among the fields of growing crystals remains an open question. Photoreactivity in single crystals has been a leading theme in organic solid-state chemistry; however, the focus of research has been traditionally on reaction stereo- and regio-specificity. However, as light-induced chemistry builds stress in crystals anisotropically, all types of motions can be actuated. The correlation between photochemistry and the responses of single crystals-jumping, twisting, fracturing, delaminating, rocking, and rolling-has become a well-defined field of research in its own right: photomechanics. The advancement of our understanding requires theoretical and high-performance computations. Computational crystallography not only supports interpretations of mechanical responses, but predicts the responses itself. This requires the engagement of classical force-field based molecular dynamics simulations, density functional theory-based approaches, and the use of machine learning to divine patterns to which algorithms can be better suited than people. The integration of mechanics with the transport of electrons and photons is considered for practical applications in flexible organic electronics and photonics. Dynamic crystals that respond rapidly and reversibly to heat and light can function as switches and actuators. Progress in identifying efficient shape-shifting crystals is also discussed. Finally, the importance of mechanical properties to milling and tableting of pharmaceuticals in an industry still dominated by active ingredients composed of small molecule crystals is reviewed. A dearth of data on the strength, hardness, Young's modulus, and fracture toughness of molecular crystals underscores the need for refinement of measurement techniques and conceptual tools. The need for benchmark data is emphasized throughout.

Can a Finite Chain of Hydrogen Cyanide Molecules Model a Crystal?

When calculating structural or spectroscopic properties of molecular crystals, the question arises whether it is sufficient to simulate only a single molecule or a small molecular cluster or whether the simulation of the entire crystal is indispensable. In this work we juxtapose calculations on the high-pressure structural properties of the (periodic) HCN crystal and chains of HCN molecules of finite length. We find that, in most cases, the behavior of the crystal can be reproduced by computational methods simulating only around 15 molecules. The pressure-induced lengthening of the C-H bond in HCN found in calculations on both the periodic and finite material are explained in terms of orbital interaction. Our results pave the way for a more thorough understanding of high-pressure structural properties of materials and give incentives for the design of materials that expand under pressure. In addition, they shed light on the complementarity between calculations on periodic materials and systems of finite size.

The Non-Classical Crystallization Mechanism of a Composite Biogenic Guanine Crystal.

Spectacular colors and visual phenomena in animals are produced by light interference from highly reflective guanine crystals. Little is known about how organisms regulate crystal morphology to tune the optics of these systems. By following guanine crystal formation in developing spiders, a crystallization mechanism is elucidated. Guanine crystallization is a "non-classical," multistep process involving a progressive ordering of states. Crystallization begins with nucleation of partially ordered nanogranules from a disordered precursor phase. Growth proceeds by orientated attachment of the nanogranules into platelets which coalesce into single crystals, via progressive relaxation of structural defects. Despite their prismatic morphology, the platelet texture is retained in the final crystals, which are composites of crystal lamellae and interlamellar sheets. Interactions between the macromolecular sheets and the planar face of guanine appear to direct nucleation, favoring platelet formation. These findings provide insights on how organisms control the morphology and optical properties of molecular crystals.

Fast Quantum Approach for Evaluating the Energy of Non-Covalent Interactions in Molecular Crystals: The Case Study of Intermolecular H-Bonds in Crystalline Peroxosolvates.

Energy/enthalpy of intermolecular hydrogen bonds (H-bonds) in crystals have been calculated in many papers. Most of the theoretical works used non-periodic models. Their applicability for describing intermolecular H-bonds in solids is not obvious since the crystal environment can strongly change H-bond geometry and energy in comparison with non-periodic models. Periodic DFT computations provide a reasonable description of a number of relevant properties of molecular crystals. However, these methods are quite cumbersome and time-consuming compared to non-periodic calculations. Here, we present a fast quantum approach for estimating the energy/enthalpy of intermolecular H-bonds in crystals. It has been tested on a family of crystalline peroxosolvates in which the H∙∙∙O bond set fills evenly (i.e., without significant gaps) the range of H∙∙∙O distances from ~1.5 to ~2.1 Å typical for strong, moderate, and weak H-bonds. Four of these two-component crystals (peroxosolvates of macrocyclic ethers and creatine) were obtained and structurally characterized for the first time. A critical comparison of the approaches for estimating the energy of intermolecular H-bonds in organic crystals is carried out, and various sources of errors are clarified.

Finding predictive models for singlet fission by machine learning

Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.

Correlation between two- and three-dimensional crystallographic lattices for epitaxial analysis. II. Experimental results.

While the crystal structure of the polymorph phase can be studied in three dimensions conveniently by X-ray methods like grazing-incidence X-ray diffraction (GIXD), the first monolayer is only accessible by surface-sensitive methods that allow the determination of a two-dimensional lattice. Here, GIXD measurements with sample rotation are compared with distortion-corrected low-energy electron diffraction (LEED) experiments on conjugated molecules: 3,4;9,10-perylenetetracarboxylic dianhydride (PTCDA), 6,13-pentacenequinone (P2O), 1,2;8,9-dibenzopentacene (trans-DBPen) and dicyanovinyl-quaterthiophene (DCV4T-Et2) grown by physical vapor deposition on Ag(111) and Cu(111) single crystals. For these molecular crystals, which exhibit different crystallographic lattices and crystal orientations as well as epitaxial properties, the geometric parameters of the three-dimensional lattice are compared with the corresponding geometry of the first monolayer. A comparison of the monolayer lattice from LEED investigations with the multilayer lattices determined by rotated GIXD experiments reveals a correlation between the first monolayer and the epitaxial growth of three-dimensional crystals together with lattice distortions and re-alignment of molecules. The selected examples show three possible scenarios of crystal growth on top of an ordered monolayer: (i) growth of a single polymorph, (ii) growth of three different polymorphs; in both cases the first monolayer serves as template. In the third case (iii) strong lattice distortion and distinct molecular re-alignments from the monolayer to epitaxially grown crystals are observed. This is the second part of our work concerning the correlation between two- and three-dimensional crystallographic lattices for epitaxial analysis. In the first part, the theoretical basis has been derived which provides a mathematical relationship between the six lattice parameters of the three-dimensional case and the three parameters obtained for the two-dimensional surface unit cell, together with their orientation to the single-crystalline substrate. In this work, a combined experimental approach of GIXD and LEED is introduced which can be used to investigate the effect of the epitaxial monolayer on the structural properties of molecular crystals grown on top.

Cleaving Method for Molecular Crystals and Its Application to Calculation of the Surface Free Energy of Crystalline β-d-Mannitol at Room Temperature.

Calculation of the surface free energy (SFE) is an important application of the thermodynamic integration (TI) methodology, which was mainly employed for atomic crystals (such as Lennard–Jones or metals). In this work, we present the calculation of the SFE of a molecular crystal using the cleaving technique which we implemented in the LAMMPS molecular dynamics package. We apply this methodology to a crystal of β-d-mannitol at room temperature and report the results for two types of force fields belonging to the GROMOS family: all atoms and united atoms. The results show strong dependence on the type of force field used, highlighting the need for the development of better force fields to model the surface properties of molecular crystals. In particular, we observe that the united-atoms force field, despite its higher degree of coarse graining compared to the all-atoms force field, produces SFE results in better agreement with the experimental results from inverse gas chromatography measurements.