Mining predicted crystal structure landscapes with high throughput crystallisation: old molecules, new insights

Structural Design and Energy and Environmental Applications of Hydrogen-Bonded Organic Frameworks: A Systematic Review

Hydrogen-bonded organic frameworks (HOFs) are emerging porous materials that show high structural flexibility, mild synthetic conditions, good solution processability, easy healing and regeneration, and good recyclability. Although these properties give them many potential multifunctional applications, their frameworks are unstable due to the presence of only weak and reversible hydrogen bonds. In this work, the development history and synthesis methods of HOFs are reviewed, and categorize their structural design concepts and strategies to improve their stability. More importantly, due to the significant potential of the latest HOF-related research for addressing energy and environmental issues, this work discusses the latest advances in the methods of energy storage and conversion, energy substance generation and isolation, environmental detection and isolation, degradation and transformation, and biological applications. Furthermore, a discussion of the coupling orientation of HOF in the cross-cutting fields of energy and environment is presented for the first time. Finally, current challenges, opportunities, and strategies for the development of HOFs to advance their energy and environmental applications are discussed.

Hydrogen-bonded organic frameworks for membrane separation

Hydrogen-bonded organic frameworks (HOFs) are a new class of crystalline porous materials that are formed through the interconnection of organic or metal-organic building units via intermolecular hydrogen bonds. The remarkable flexibility and reversibility of hydrogen bonds, coupled with the customizable nature of organic units, endow HOFs with mild synthesis conditions, high crystallinity, solvent processability, and facile self-healing and regeneration properties. Consequently, these features have garnered significant attention across various fields, particularly in the realm of membrane separation. Herein, we present an overview of the recent advances in HOF-based membranes, including their advanced fabrication strategies and fascinating applications in membrane separation. To attain the desired HOF-based membranes, careful consideration is dedicated to crucial factors such as pore size, stability, hydrophilicity/hydrophobicity, and surface charge of the HOFs. Additionally, diverse preparation methods for HOF-based membranes, including blending, in situ growth, solution-processing, and electrophoretic deposition, have been analyzed. Furthermore, applications of HOF-based membranes in gas separation, water treatment, fuel cells, and other emerging application areas are presented. Finally, the challenges and prospects of HOF-based membranes are critically pointed out.

Hierarchical Self-Assembly of Capsule-Shaped Zirconium Coordination Cages with Quaternary Structure

Biological macromolecules exhibit emergent functions through hierarchical self-assembly, a concept that is extended to design artificial supramolecular assemblies. Here, the first example of breaking the common parallel arrangement of capsule-shaped zirconium coordination cages is reported by constructing the hierarchical porous framework ZrR-1. ZrR-1 adopts a quaternary structure resembling protein and contains 12-connected chloride clusters, representing the highest connectivity for zirconium-based cages reported thus far. Compared to the parallel framework ZrR-2, ZrR-1 demonstrated enhanced stability in acidic aqueous solutions and a tenfold increase in BET surface area (879 m2  g-1 ). ZrR-1 also exhibits excellent proton conductivity, reaching 1.31 × 10-2 S·cm-1 at 353 K and 98% relative humidity, with a low activation energy of 0.143 eV. This finding provides insights into controlling the hierarchical self-assembly of metal-organic cages to discover superstructures with emergent properties.

Modular, multi-robot integration of laboratories: an autonomous workflow for solid-state chemistry

Automation can transform productivity in research activities that use liquid handling, such as organic synthesis, but it has made less impact in materials laboratories, which require sample preparation steps and a range of solid-state characterization techniques. For example, powder X-ray diffraction (PXRD) is a key method in materials and pharmaceutical chemistry, but its end-to-end automation is challenging because it involves solid powder handling and sample processing. Here we present a fully autonomous solid-state workflow for PXRD experiments that can match or even surpass manual data quality, encompassing crystal growth, sample preparation, and automated data capture. The workflow involves 12 steps performed by a team of three multipurpose robots, illustrating the power of flexible, modular automation to integrate complex, multitask laboratories.

Controlling the Crystallisation and Hydration State of Crystalline Porous Organic Salts

Crystalline porous organic salts (CPOS) are a subclass of molecular crystals. The low solubility of CPOS and their building blocks limits the choice of crystallisation solvents to water or polar alcohols, hindering the isolation, scale-up, and scope of the porous material. In this work, high throughput screening was used to expand the solvent scope, resulting in the identification of a new porous salt, CPOS-7, formed from tetrakis(4-sulfophenyl)methane (TSPM) and tetrakis(4-aminophenyl)methane (TAPM). CPOS-7 does not form with standard solvents for CPOS, rather a hydrated phase (Hydrate2920) previously reported is isolated. Initial attempts to translate the crystallisation to batch led to challenges with loss of crystallinity and Hydrate2920 forming favorably in the presence of excess water. Using acetic acid as a dehydrating agent hindered formation of Hydrate2920 and furthermore allowed for direct conversion to CPOS-7. To allow for direct formation of CPOS-7 in high crystallinity flow chemistry was used for the first time to circumvent the issues found in batch. CPOS-7 and Hydrate2920 were shown to have promise for water and CO2 capture, with CPOS-7 having a CO2 uptake of 4.3 mmol/g at 195 K, making it one of the most porous CPOS reported to date.

Crystal structure of 2-methyl-1H-imidazol-3-ium 3,5-di-carb-oxy-benzoate

The structure of the title salt, C4H7N2 +·C9H5O6 - (1), is reported. The compound is built from a protonated 2-methyl-imidazole and a singly deprotonated trimesic acid. Detailed analysis of bond distances and angles for both ions reveals subtle differences compared with their neutral mol-ecule counterpart. Analysis of the crystal packing in compound 1 reveals the formation of undulating chains by the ions through hydrogen bonding. The chains stack along the b axis through π-π inter-actions and inter-connect with other chains in an out-of-phase arrangement along the ac plane through further hydrogen-bonding inter-actions.

Selective adsorption of polycyclic aromatic hydrocarbons by isostructural hydrogen-bonded organic frameworks

Two isostructural hydrogen-bonded organic frameworks (HOFs) with 1-D hexagonal-shaped pores were crystallised using the molecules biphenyl-3,3',5,5'-tetracarboxylic acid (BPTCA) and [1,1':4',1'']terphenyl- 3,3'',5,5''-tetracarboxylic acid (TPTCA). The desolvated HOFs, named BPTCA-2 and TPTCA-2, exhibited selective adsorption towards naphthalene and anthracene, respectively, during competitive adsorption experiments.

From liquid to crystal via mechanochemical grinding: unique host-guest (HOF) cocrystal

Mechanochemical synthesis via grinding of trimesic acid (TA, C9H6O6) and 4-chlorophenyl diphenyl phosphate (4CDP, C18H14ClO4P) (liquid at room temperature) in a 1:1 ratio resulted in the formation of an inclusion type of cocrystal. The crystallization of this phase via slow evaporation at low temperature (276-277 K) from methanol resulted in a rare `stairstep morphology' during the process of crystal growth. This morphology was not observed after crystallization of the compound from other solvents like toluene, dichloromethane, acetone, hexane and isooctane, and hence this was characteristically observed in methanol only. The characterization from single-crystal X-ray diffraction revealed the formation of a cocrystal with five molecules of TA and two molecules of 4CDP in the asymmetric unit. The trimesic acid molecules form hydrogen-bonded dimers resulting in hexagonal rings, and these rings are stacked through π-π intermolecular interactions to make a hexagonal honeycomb-like structure. The phosphate molecules, 4CDP, were found to be trapped as guests in these hexagonal channels. The similarity in the packing of trimesic acid is compared in the cocrystal and the free acid quantitatively via Xpac analysis, which establishes the relationship of a `2D supramolecular construct' between them. This signifies a unique type of arrangement in which the voids created by the trimesic acid moiety do not undergo distortion by the inclusion of the guest molecules. The quantitative analysis of the intermolecular interactions using Hirshfeld surfaces and fingerprint plots deciphers the role of both strong O-H...O hydrogen bonds and weak intermolecular interactions in the crystal packing.

Structure prediction in low dimensions: concepts, issues and examples

Structure prediction of stable and metastable polymorphs of chemical systems in low dimensions has become an important field, since materials that are patterned on the nano-scale are of increasing importance in modern technological applications. While many techniques for the prediction of crystalline structures in three dimensions or of small clusters of atoms have been developed over the past three decades, dealing with low-dimensional systems-ideal one-dimensional and two-dimensional systems, quasi-one-dimensional and quasi-two-dimensional systems, as well as low-dimensional composite systems-poses its own challenges that need to be addressed when developing a systematic methodology for the determination of low-dimensional polymorphs that are suitable for practical applications. Quite generally, the search algorithms that had been developed for three-dimensional systems need to be adjusted when being applied to low-dimensional systems with their own specific constraints; in particular, the embedding of the (quasi-)one-dimensional/two-dimensional system in three dimensions and the influence of stabilizing substrates need to be taken into account, both on a technical and a conceptual level. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Reducing overprediction of molecular crystal structures via threshold clustering

Crystal structure prediction is becoming an increasingly valuable tool for assessing polymorphism of crystalline molecular compounds, yet invariably, it overpredicts the number of polymorphs. One of the causes for this overprediction is in neglecting the coalescence of potential energy minima, separated by relatively small energy barriers, into a single basin at finite temperature. Considering this, we demonstrate a method underpinned by the threshold algorithm for clustering potential energy minima into basins, thereby identifying kinetically stable polymorphs and reducing overprediction.

Charge-Assisted Hydrogen-Bonded Organic Frameworks with Inorganic Ammonium Regulated Switchable Open Polar Sites

Surface open polar sites within the voids of porous molecular crystals define the localized physicochemical environment for critical functions such as gas separation and molecular recognition. This study presents a new charge-assisted hydrogen bonding (H-bonding) motif, by exploiting inorganic ammonium (NH₄ ⁺ ) cations as H-bond donors, to regulate the assembly of C₂ -symmetric carboxylic tectons for building robust H-bonded frameworks with permanent ultra-micropores and open oxygen sites. Diverse building blocks are bridged by tetrahedral NH₄ ⁺ to expand distinctive H-bonded networks with varied pore architectures. Particularly, the open polar oxygen sites can be switched by altering NH₄ ⁺ sources to tune the deprotonation of carboxyl-containing tectons. The activated porous PTBA·NH₄ ·DMF preserves the pore architecture and open polar oxygen sites, exhibiting remarkably selective sorption of CO₂ (107.8 cm³ g^-1 ,195 K) over N₂ (11.2 cm³ g^-1 , 77 K) and H₂ (1.4 cm³ g^-1 , 77 K).

A theoretical framework for the design of molecular crystal engines

Photomechanical molecular crystals have garnered attention for their ability to transform light into mechanical work, but difficulties in characterizing the structural changes and mechanical responses experimentally have hindered the development of practical organic crystal engines. This study proposes a new computational framework for predicting the solid-state crystal-to-crystal photochemical transformations entirely from first principles, and it establishes a photomechanical engine cycle that quantifies the anisotropic mechanical performance resulting from the transformation. The approach relies on crystal structure prediction, solid-state topochemical principles, and high-quality electronic structure methods. After validating the framework on the well-studied [4 + 4] cycloadditions in 9-methyl anthracene and 9-tert-butyl anthracene ester, the experimentally-unknown solid-state transformation of 9-carboxylic acid anthracene is predicted for the first time. The results illustrate how the mechanical work is done by relaxation of the crystal lattice to accommodate the photoproduct, rather than by the photochemistry itself. The large ∼10⁷ J m^-3 work densities computed for all three systems highlight the promise of photomechanical crystal engines. This study demonstrates the importance of crystal packing in determining molecular crystal engine performance and provides tools and insights to design improved materials in silico.

Synthesis, Crystal Packing Aspects and Pseudosymmetry in Coordination Compounds with a Phosphorylamide Ligand

This work reports the synthesis and crystal structure of new closely related coordination compounds, [ML2]·nTHF, where M is Zn or Mn; L is a phosphorylmethylamide derivative of benzothiadiazole; n = 1 (M = Zn) and 1, 2 (M = Mn); and THF is tetrahydrofuran. The zinc compound, 1·THF, crystallizes in a high-symmetry space group, I41/a, that is relatively rare for compounds with organic ligands. The corresponding manganese congener, 2·THF, with a similar crystal packing, features a pseudosymmetrical structure P21/c of the doubled volume of the unit cell as compared to 1·THF. The main difference between the structures lies in a different arrangement of solvate THF molecules, which likely modulates the crystal packing of the complexes. Another manganese solvatomorph, 2·2THF, reveals a fundamentally different crystal packing while exhibiting a similar geometry of the complex. We consider the problem of localization of solvate THF molecules and the types of their disorder by the example of compounds 1–2.

Multifunctional Porous Hydrogen-Bonded Organic Frameworks: Current Status and Future Perspectives

Hydrogen-bonded organic frameworks (HOFs), self-assembled from organic or metalated organic building blocks (also termed as tectons) by hydrogen bonding, π-π stacking, and other intermolecular interactions, have become an emerging class of multifunctional porous materials. So far, a library of HOFs with high porosity has been synthesized based on versatile tectons and supramolecular synthons. Benefiting from the flexibility and reversibility of H-bonds, HOFs feature high structural flexibility, mild synthetic reaction, excellent solution processability, facile healing, easy regeneration, and good recyclability. However, the flexible and reversible nature of H-bonds makes most HOFs suffer from poor structural designability and low framework stability. In this Outlook, we first describe the development and structural features of HOFs and summarize the design principles of HOFs and strategies to enhance their stability. Second, we highlight the state-of-the-art development of HOFs for diverse applications, including gas storage and separation, heterogeneous catalysis, biological applications, sensing, proton conduction, and other applications. Finally, current challenges and future perspectives are discussed.

Water Sorption Controls Extreme Single-Crystal-to-Single-Crystal Molecular Reorganization in Hydrogen Bonded Organic Frameworks

As hydrogen bonded frameworks are held together by relatively weak interactions, they often form several different frameworks under slightly different synthesis conditions and respond dynamically to stimuli such as heat and vacuum. However, these dynamic restructuring processes are often poorly understood. In this work, three isoreticular hydrogen bonded organic frameworks assembled through charge-assisted amidinium⋅⋅⋅carboxylate hydrogen bonds (1^C/C , 1^Si/C and 1^Si/Si ) are studied. Three distinct phases for 1^C/C and four for 1^Si/C and 1^Si/Si are fully structurally characterized. The transitions between these phases involve extreme yet recoverable molecular-level framework reorganization. It is demonstrated that these transformations are related to water content and can be controlled by humidity, and that the non-porous anhydrous phase of 1^C/C shows reversible water sorption through single crystal to crystal restructuring. This mechanistic insight opens the way for the future use of the inherent dynamism present in hydrogen bonded frameworks.

Global analysis of the energy landscapes of molecular crystal structures by applying the threshold algorithm

Polymorphism in molecular crystals has important consequences for the control of materials properties and our understanding of crystallization. Computational methods, including crystal structure prediction, have provided important insight into polymorphism, but have usually been limited to assessing the relative energies of structures. We describe the implementation of the Monte Carlo threshold algorithm as a method to provide an estimate of the energy barriers separating crystal structures. By sampling the local energy minima accessible from multiple starting structures, the simulations yield a global picture of the crystal energy landscapes and provide valuable information on the depth of the energy minima associated with crystal structures. We present results from applying the threshold algorithm to four polymorphic organic molecular crystals, examine the influence of applying space group symmetry constraints during the simulations, and discuss the relationship between the structure of the energy landscape and the intermolecular interactions present in the crystals.

Machine-Learning-Guided Identification of Coordination Polymer Ligands for Crystallizing Separation of Cs/Sr

Separation of Cs/Sr is one of many coordination-chemistry-centered processes in the grand scheme of spent nuclear fuel reprocessing, a critical link for a sustainable nuclear energy industry. To deploy a crystallizing Cs/Sr separation technology, we planned to systematically screen and identify candidate ligands that can efficiently and selectively bind to Sr2+ and form coordination polymers. Therefore, we mined the Cambridge Structural Database for characteristic structural information and developed a machine-learning-guided methodology for ligand evaluation. The optimized machine-learning model, correlating the molecular structures of the ligands with the predicted coordinative properties, generated a ranking list of potential compounds for Cs/Sr selective crystallization. The Sr2+ sequestration capability and selectivity over Cs+ of the promising ligands identified (squaric acid and chloranilic acid) were subsequently confirmed experimentally, with commendable performances, corroborating the artificial-intelligence-guided strategy.

Analogy Powered by Prediction and Structural Invariants: Computationally Led Discovery of a Mesoporous Hydrogen-Bonded Organic Cage Crystal

Mesoporous molecular crystals have potential applications in separation and catalysis, but they are rare and hard to design because many weak interactions compete during crystallization, and most molecules have an energetic preference for close packing. Here, we combine crystal structure prediction (CSP) with structural invariants to continuously qualify the similarity between predicted crystal structures for related molecules. This allows isomorphous substitution strategies, which can be unreliable for molecular crystals, to be augmented by a priori prediction, thus leveraging the power of both approaches. We used this combined approach to discover a rare example of a low-density (0.54 g cm^–3) mesoporous hydrogen-bonded framework (HOF), 3D-CageHOF-1. This structure comprises an organic cage (Cage-3-NH₂) that was predicted to form kinetically trapped, low-density polymorphs via CSP. Pointwise distance distribution structural invariants revealed five predicted forms of Cage-3-NH₂ that are analogous to experimentally realized porous crystals of a chemically different but geometrically similar molecule, T2. More broadly, this approach overcomes the difficulties in comparing predicted molecular crystals with varying lattice parameters, thus allowing for the systematic comparison of energy–structure landscapes for chemically dissimilar molecules.

HOFs Built from Hexatopic Carboxylic Acids: Structure, Porosity, Stability, and Photophysics

Hydrogen-bonded organic frameworks (HOFs) have attracted renewed attention as another type of promising candidates for functional porous materials. In most cases of HOF preparation, the applied molecular design principle is based on molecules with rigid π-conjugated skeleton together with more than three H-bonding groups to achieve 2D- or 3D-networked structures. However, the design principle does not always work, but results in formation of unexpected structures, where subtle structural factors of which we are not aware dictate the entire structure of HOFs. In this contribution, we assess recent advances in HOFs, focusing on those composed of hexatopic building block molecules, which can provide robust frameworks with a wide range of topologies and properties. The HOFs described in this work are classified into three types, depending on their H-bonded structural motifs. Here in, we focus on: (1) the chemical aspects that govern their unique fundamental chemistry and structures; and (2) their photophysics at the ensemble and single-crystal levels. The work addresses and discusses how these aspects affect and orient their photonic applicability. We trust that this contribution will provide a deep awareness and will help scientists to build up a systematic series of porous materials with the aim to control both their structural and photodynamical assets.

The extensive solid-form landscape of sulfathiazole: hydrogen-bond topology and node shape

Patterns of hydrogen bonds are described in a set of 101 crystal structures containing sulfathiazole. Topological analysis of the hydrogen-bond nets is augmented by comparison of the shapes of the nodes extracted from each net.

Luminescent organic porous crystals from non-cyclic molecules and their applications

Organic porous crystals from small and non-cyclic organic molecules can be constructed by various intermolecular weak interactions. Owing to their precise stacking types, intermolecular interaction and pore microstructure, the relationship...

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.

Enabling Technology for Supramolecular Chemistry

Supramolecular materials-materials that exploit non-covalent interactions-are increasing in structural complexity, selectivity, function, stability, and scalability, but their use in applications has been comparatively limited. In this Minireview, we summarize the opportunities presented by enabling technology-flow chemistry, high-throughput screening, and automation-to wield greater control over the processes in supramolecular chemistry and accelerate the discovery and use of self-assembled systems. Finally, we give an outlook for how these tools could transform the future of the field.

Recent progress and future challenges in the supramolecular polymerization of metal-containing monomers

The self-assembly of discrete molecular entities into functional nanomaterials has become a major research area in the past decades. The library of investigated compounds has diversified significantly, while the field as a whole has matured. The incorporation of metal ions in the molecular design of the (supra-)molecular building blocks greatly expands the potential applications, while also offering a promising approach to control molecular recognition and attractive and/or repulsive intermolecular binding events. Hence, supramolecular polymerization of metal-containing monomers has emerged as a major research focus in the field. In this perspective article, we highlight recent significant advances in supramolecular polymerization of metal-containing monomers and discuss their implications for future research. Additionally, we also outline some major challenges that metallosupramolecular chemists (will) have to face to produce metallosupramolecular polymers (MSPs) with advanced applications and functionalities.

Accelerating computational discovery of porous solids through improved navigation of energy-structure-function maps

While energy-structure-function (ESF) maps are a powerful new tool for in silico materials design, the cost of acquiring an ESF map for many properties is too high for routine integration into high-throughput virtual screening workflows. Here, we propose the next evolution of the ESF map. This uses parallel Bayesian optimization to selectively acquire energy and property data, generating the same levels of insight at a fraction of the computational cost. We use this approach to obtain a two orders of magnitude speedup on an ESF study that focused on the discovery of molecular crystals for methane capture, saving more than 500,000 central processing unit hours from the original protocol. By accelerating the acquisition of insight from ESF maps, we pave the way for the use of these maps in automated ultrahigh-throughput screening pipelines by greatly reducing the opportunity risk associated with the choice of system to calculate.

Crystallography companion agent for high-throughput materials discovery

The discovery of new structural and functional materials is driven by phase identification, often using X-ray diffraction (XRD). Automation has accelerated the rate of XRD measurements, greatly outpacing XRD analysis techniques that remain manual, time-consuming, error-prone and impossible to scale. With the advent of autonomous robotic scientists or self-driving laboratories, contemporary techniques prohibit the integration of XRD. Here, we describe a computer program for the autonomous characterization of XRD data, driven by artificial intelligence (AI), for the discovery of new materials. Starting from structural databases, we train an ensemble model using a physically accurate synthetic dataset, which outputs probabilistic classifications-rather than absolutes-to overcome the overconfidence in traditional neural networks. This AI agent behaves as a companion to the researcher, improving accuracy and offering substantial time savings. It is demonstrated on a diverse set of organic and inorganic materials characterization challenges. This method is directly applicable to inverse design approaches and robotic discovery systems, and can be immediately considered for other forms of characterization such as spectroscopy and the pair distribution function.

Integrating Computational and Experimental Workflows for Accelerated Organic Materials Discovery

Organic materials find application in a range of areas, including optoelectronics, sensing, encapsulation, molecular separations, and photocatalysis. The discovery of materials is frustratingly slow however, particularly when contrasted to the vast chemical space of possibilities based on the near limitless options for organic molecular precursors. The difficulty in predicting the material assembly, and consequent properties, of any molecule is another significant roadblock to targeted materials design. There has been significant progress in the development of computational approaches to screen large numbers of materials, for both their structure and properties, helping guide synthetic researchers toward promising materials. In particular, artificial intelligence techniques have the potential to make significant impact in many elements of the discovery process. Alongside this, automation and robotics are increasing the scale and speed with which materials synthesis can be realized. Herein, the focus is on demonstrating the power of integrating computational and experimental materials discovery programmes, including both a summary of key situations where approaches can be combined and a series of case studies that demonstrate recent successes.

Digital navigation of energy-structure-function maps for hydrogen-bonded porous molecular crystals

Energy-structure-function (ESF) maps can aid the targeted discovery of porous molecular crystals by predicting the stable crystalline arrangements along with their functions of interest. Here, we compute ESF maps for a series of rigid molecules that comprise either a triptycene or a spiro-biphenyl core, functionalized with six different hydrogen-bonding moieties. We show that the positioning of the hydrogen-bonding sites, as well as their number, has a profound influence on the shape of the resulting ESF maps, revealing promising structure-function spaces for future experiments. We also demonstrate a simple and general approach to representing and inspecting the high-dimensional data of an ESF map, enabling an efficient navigation of the ESF data to identify 'landmark' structures that are energetically favourable or functionally interesting. This is a step toward the automated analysis of ESF maps, an important goal for closed-loop, autonomous searches for molecular crystals with useful functions.

What's in an Atom? A Comparison of Carbon and Silicon-Centred Amidinium⋅⋅⋅Carboxylate Frameworks*

Despite their apparent similarity, framework materials based on tetraphenylmethane and tetraphenylsilane building blocks often have quite different structures and topologies. Herein, we describe a new silicon tetraamidinium compound and use it to prepare crystalline hydrogen bonded frameworks with carboxylate anions in water. The silicon-containing frameworks are compared with those prepared from the analogous carbon tetraamidinium: when biphenyldicarboxylate or tetrakis(4-carboxyphenyl)methane anions were used similar channel-containing networks are observed for both the silicon and carbon tetraamidinium. When terephthalate or bicarbonate anions were used, different products form. Insights into possible reasons for the different products are provided by a survey of the Cambridge Structural Database and quantum chemical calculations, both of which indicate that, contrary to expectations, tetraphenylsilane derivatives have less geometrical flexibility than tetraphenylmethane derivatives, that is, they are less able to distort away from ideal tetrahedral bond angles.

Dynamic structural reconstruction of (guanidinium+)2(benzene-1,4-disulfonate2−) host crystal by guest adsorption

Guanidinium (G⁺) and benzene-1,4-disulfonate (BDS²⁻) form a rigid electrostatic cation–anion crystal lattice, which undergoes an interesting dynamic structural reconstruction through guest adsorption–desorption processes.

Reactive crystallization: a review

Reactive crystallization is not new, but there has been recent growth in its use as a means of improving performance and sustainability of industrial processes.

Structure prediction of crystals, surfaces and nanoparticles

We review the current techniques used in the prediction of crystal structures and their surfaces and of the structures of nanoparticles. The main classes of search algorithm and energy function are summarized, and we discuss the growing role of methods based on machine learning. We illustrate the current status of the field with examples taken from metallic, inorganic and organic systems. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.

Hidden Solvates and Transient Forms of Trimesic Acid

This article discusses the formation of trimesic acid (TMA) solvates with ethanol, isopropyl alcohol and dimethylformamide via liquid-assisted grinding and slurry experiments. Through the use of X-ray diffraction methods, we highlight the formation of a new ethanol solvate of TMA that completes the series of alcohol solvates observed, a temperature-induced phase transition in the isopropyl alcohol solvate between 233 K and 243 K, and a transient 1:3 solvate with dimethylformamide that mimics a previously identified dimethylsulfoxide solvate. The alcohol structures possess a TMA framework that is geometrically similar where the intermolecular energies between TMA molecules are equivalent. We have observed that increasing the length of the alcohol induces an increase in the distortion of the TMA framework to accommodate the longer alkyl tails.

Minimizing Polymorphic Risk through Cooperative Computational and Experimental Exploration

We combine state-of-the-art computational crystal structure prediction (CSP) techniques with a wide range of experimental crystallization methods to understand and explore crystal structure in pharmaceuticals and minimize the risk of unanticipated late-appearing polymorphs. Initially, we demonstrate the power of CSP to rationalize the difficulty in obtaining polymorphs of the well-known pharmaceutical isoniazid and show that CSP provides the structure of the recently obtained, but unsolved, Form III of this drug despite there being only a single resolved form for almost 70 years. More dramatically, our blind CSP study predicts a significant risk of polymorphism for the related iproniazid. Employing a wide variety of experimental techniques, including high-pressure experiments, we experimentally obtained the first three known nonsolvated crystal forms of iproniazid, all of which were successfully predicted in the CSP procedure. We demonstrate the power of CSP methods and free energy calculations to rationalize the observed elusiveness of the third form of iproniazid, the success of high-pressure experiments in obtaining it, and the ability of our synergistic computational-experimental approach to “de-risk” solid form landscapes.

High-Throughput Approaches for the Discovery of Supramolecular Organic Cages

The assembly of complex molecules, such as organic cages, can be achieved through supramolecular and dynamic covalent strategies. Their use in a range of applications has been demonstrated, including gas uptake, molecular separations, and in catalysis. However, the targeted design and synthesis of new species for particular applications is challenging, particularly as the systems become more complex. High-throughput computation-only and experiment-only approaches have been developed to streamline the discovery process, although are still not widely implemented. Additionally, combined hybrid workflows can dramatically accelerate the discovery process and lead to the serendipitous discovery and rationalisation of new supramolecular assemblies that would not have been designed based on intuition alone. This Minireview focuses on the advances in high-throughput approaches that have been developed and applied in the discovery of supramolecular organic cages.

Hydrogen-Bonded Organic Frameworks as a Tunable Platform for Functional Materials

As a novel class of porous crystalline materials, hydrogen-bonded organic frameworks (HOFs), self-assembled from organic or metal-organic building blocks through intermolecular hydrogen-bonding interactions, have attracted more and more attention. Over the past decade, a number of porous HOFs have been constructed through judicious selection of H-bonding motifs, which are further enforced by other weak intermolecular interactions such as π-π stacking and van der Waals forces and framework interpenetration. Since the H-bonds are weaker than coordinate and covalent bonds used for the construction of metal-organic frameworks (MOFs) and covalent organic frameworks (COFs), HOFs have some unique features such as mild synthesis condition, solution processability, easy healing, and regeneration. These features enable HOFs to be a tunable platform for the construction of functional materials. Here, we review the H-bonding motifs used for constructing porous HOFs and highlight some of their applications, including gas separation and storage, chiral separation and structure determination, fluorescent sensing, heterogeneous catalysis, biological applications, proton conduction, photoluminescent materials, and membrane-based applications.

An Expandable Hydrogen-Bonded Organic Framework Characterized by Three-Dimensional Electron Diffraction

A molecular crystal of a 2-D hydrogen-bonded organic framework (HOF) undergoes an unusual structural transformation after solvent removal from the crystal pores during activation. The conformationally flexible host molecule, ABTPA, adapts its molecular conformation during activation to initiate a framework expansion. The microcrystalline activated phase was characterized by three-dimensional electron diffraction (3D ED), which revealed that ABTPA uses out-of-plane anthracene units as adaptive structural anchors. These units change orientation to generate an expanded, lower density framework material in the activated structure. The porous HOF, ABTPA-2, has robust dynamic porosity (SA_BET = 1183 m² g^-1) and exhibits negative area thermal expansion. We use crystal structure prediction (CSP) to understand the underlying energetics behind the structural transformation and discuss the challenges facing CSP for such flexible molecules.

Efficient Screening for Ternary Molecular Ionic Cocrystals Using a Complementary Mechanosynthesis and Computational Structure Prediction Approach

The discovery of molecular ionic cocrystals (ICCs) of active pharmaceutical ingredients (APIs) widens the opportunities for optimizing the physicochemical properties of APIs whilst facilitating the delivery of multiple therapeutic agents. However, ICCs are often observed serendipitously in crystallization screens and the factors dictating their crystallization are poorly understood. We demonstrate here that mechanochemical ball milling is a versatile technique for the reproducible synthesis of ternary molecular ICCs in less than 30 min of grinding with or without solvent. Computational crystal structure prediction (CSP) calculations have been performed on ternary molecular ICCs for the first time and the observed crystal structures of all the ICCs were correctly predicted. Periodic dispersion-corrected DFT calculations revealed that all the ICCs are thermodynamically stable (mean stabilization energy=-2 kJ mol-1 ) relative to the crystallization of a physical mixture of the binary salt and acid. The results suggest that a combined mechanosynthesis and CSP approach could be used to target the synthesis of higher-order molecular ICCs with functional properties.

Morphological Evolution of Two-Dimensional Porous Hexagonal Trimesic Acid Framework

Hexagonal single crystal structure (Form II) of trimesic acid (TMA) has been isolated by dissolving the interpenetrated Form I of TMA in tetrahydrofuran. Form II (hexagonal) was converted to Form I (interpenetrated) at room temperature through some intermediate structures. A detailed time-dependent FESEM study shows that the external morphology of Form II (hexagonal) is a hollow hexagonal tube that mimics its crystal structure. The block-shaped (morphology) of Form I (interpenetrated) was converted to the hollow hexagonal tube through some intermediate morphologies which are corresponding to particular crystal structures. Here, we have established a strong correlation between crystal structures with the morphology. These hollow tubes have been employed for Rhodamine B dye adsorption studies and showed an uptake of 82%, much more significant than Form I (interpenetrated) (39%) due to the presence of a pore channel in the crystal structure.

In vivo real-time tracking of tumor-specific biocatalysis in cascade nanotheranostics enables synergistic cancer treatment

Glucose oxidase (GOD)-based synergistic cancer therapy has aroused great research interest in the context of cancer treatment due to the inherent biocompatibility and biodegradability. However, this emerging therapeutic system still lacks a strategy to predict and regulate the in vivo biocatalytic behavior of GOD in real time to minimize the side effects on normal tissues. Herein, we developed a tumor-specific cascade nanotheranostic system (BNG) that combines GOD-catalyzed oxidative stress and dual-channel fluorescent sensing, significantly improving the synergistic therapeutic efficacy with real-time feedback information. The nanotheranostic system remains completely silent in the blood circulatory system and selectively releases GOD enzymes in the tumor site, with enhanced near-infrared (NIR) fluorescence at 825 nm. Subsequently, GOD catalyzes H₂O₂ production, triggering cascade reactions with NIR fluorescence at 650 nm as an optical output, along with GSH depletion, enabling synergistic cancer treatment. The designed nanotheranostic system, integrated with tumor-activated cascade reactions and triggering a dual-channel output at each step, represents an insightful paradigm for precise cooperative cancer therapy.

Herein, we described the strategy of cascade nanotheranostics for in real-time tracking of GOD release and thus provide the feedback information of biocatalysis cascade process with significantly enhanced in vivo therapeutic efficiency.

Hydrogen-bonded porous frameworks constructed by rigid π-conjugated molecules with carboxy groups

This review covers construction and properties of porous molecular crystals (PMCs) constructed through hydrogen-bonding of C3-symmetric, rigid, π-conjugated molecular building blocks possessing carboxyaryl groups, which was reported in the last 5 years by the author’s group. PMCs with well-defined, self-standing pores have been attracted attention due to various functionalities provided by selective and reversible inclusion of certain chemical species into the pores. However, it has been recognized for long time that construction of PMCs with permanent porosity is not easy due to weakness of noncovalent intermolecular interactions. Systematic construction of PMCs have been limited so far. To overcome this problem, the author has proposed a unique molecular design concept based on C3-symmetric π-conjugated molecules (C3PIs) possessing o-bis(4-carboxyphenyl)benzene moieties in their periphery and demonstrated that C3PIs systematically yielded hydrogen-bonded organic frameworks (HOFs) composed of H-bonded 2D hexagonal networks (H-HexNets) or interpenetrated 3D pcu-networks, which exhibit permanent porosity, significant thermal stability, polar solvent durability, robustness/flexibility, and/or multifunctionality.