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

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.

Supramolecular polynuclear clusters sustained cubic hydrogen bonded frameworks with octahedral cages for reversible photochromism

Developing supramolecular porous crystalline frameworks with tailor-made architectures from advanced secondary building units (SBUs) remains a pivotal challenge in reticular chemistry. Particularly for hydrogen-bonded organic frameworks (HOFs), construction of geometrical cavities through secondary units has been rarely achieved. Herein, a body-centered cubic HOF (TCA_NH4) with octahedral cages was constructed by a C3-symmetric building block and NH4+ node-assembled cluster (NH4)4(COOH)8(H2O)2 that served as supramolecular secondary building units (SSBUs), akin to the polynuclear SBUs in reticular chemistry. Specifically, the octahedral cages could encapsulate four homogenous haloforms including CHCl3, CHBr3, and CHI3 with truncated octahedron configuration. Crystallographic evidence revealed the cages served as spatially-confined nanoreactors, enabling fast, broadband photochromic effect associated with the reversible photo/thermal transformation between encapsulated CHI3 and I2. Overall, this work provides a strategy by shaping SSBUs to expand the framework topology of HOFs and a prototype of hydrogen-bonded nanoreactors to accommodate reversible photochromic reactions.

Machine-Learned Potentials by Active Learning from Organic Crystal Structure Prediction Landscapes

A primary challenge in organic molecular crystal structure prediction (CSP) is accurately ranking the energies of potential structures. While high-level solid-state density functional theory (DFT) methods allow for mostly reliable discrimination of the low-energy structures, their high computational cost is problematic because of the need to evaluate tens to hundreds of thousands of trial crystal structures to fully explore typical crystal energy landscapes. Consequently, lower-cost but less accurate empirical force fields are often used, sometimes as the first stage of a hierarchical scheme involving multiple stages of increasingly accurate energy calculations. Machine-learned interatomic potentials (MLIPs), trained to reproduce the results of ab initio methods with computational costs close to those of force fields, can improve the efficiency of the CSP by reducing or eliminating the need for costly DFT calculations. Here, we investigate active learning methods for training MLIPs with CSP datasets. The combination of active learning with the well-developed sampling methods from CSP yields potentials in a highly automated workflow that are relevant over a wide range of the crystal packing space. To demonstrate these potentials, we illustrate efficiently reranking large, diverse crystal structure landscapes to near-DFT accuracy from force field-based CSP, improving the reliability of the final energy ranking. Furthermore, we demonstrate how these potentials can be extended to more accurately model structures far from lattice energy minima through additional on-the-fly training within Monte Carlo simulations.

Influence of donor point modifications on the assembly of chalcogen-bonded organic frameworks

Incremental, single-atom substitutions of Se-based chalcogen bond (Ch-bond) donors with stronger donating Te centers were implemented in two new triptycene tris(1,2,5-chalcogenadiazole) tectons. The appreciably more favorable Ch-bonding ability of the Te-based donors promotes assembly of low-density networks and more stable Ch-bonded organic frameworks (ChOFs).

Continuous chiral distances for two-dimensional lattices

Chirality was traditionally considered a binary property of periodic lattices and crystals. However, the classes of two-dimensional lattices modulo rigid motion form a continuous space, which was recently parametrized by three geographic-style coordinates. The four non-oblique Bravais classes of two-dimensional lattices form low-dimensional singular subspaces in the full continuous space. Now, the deviations of a lattice from its higher symmetry neighbors can be continuously quantified by real-valued distances satisfying metric axioms. This article analyzes these and newer G-chiral distances for millions of two-dimensional lattices that are extracted from thousands of available two-dimensional materials and real crystal structures in the Cambridge Structural Database.

Soft Hydrogen-Bonded Organic Frameworks Constructed Using a Flexible Organic Cage Hinge

Soft porous crystals combine flexibility and porosity, allowing them to respond structurally to external physical and chemical environments. However, striking the right balance between flexibility and sufficient rigidity for porosity is challenging, particularly for molecular crystals formed by using weak intermolecular interactions. Here, we report a flexible oxygen-bridged prismatic organic cage molecule, Cage-6-COOH, which has three pillars that exhibit "hinge-like" rotational motion in the solid state. Cage-6-COOH can form a range of hydrogen-bonded organic frameworks (HOFs) where the "hinge" can accommodate a remarkable 67° dihedral angle range between neighboring units. This stems both from flexibility in the noncovalent hydrogen-bonding motifs in the HOFs and the molecular flexibility in the oxygen-linked cage hinge itself. The range of structures for Cage-6-COOH includes two topologically complex interpenetrated HOFs, CageHOF-2α and CageHOF-2β. CageHOF-2α is nonporous, while CageHOF-2β has permanent porosity and a surface area of 458 m2 g-1. The flexibility of Cage-6-COOH allows this molecule to rapidly transform from a low-crystallinity solid into the two crystalline interpenetrated HOFs, CageHOF-2α and CageHOF-2β, under mild conditions simply by using acetonitrile or ethanol vapor, respectively. This self-healing behavior was selective, with the CageHOF-2β structure exhibiting structural memory behavior.

Pore Modulation of Hydrogen-Bonded Organic Frameworks for Efficient Separation of Propylene

Developing hydrogen-bonded organic frameworks (HOFs) that combine functional sites, size control, and storage capability for targeting gas molecule capture is a novel and challenging venture. However, there is a lack of effective strategies to tune the hydrogen-bonded network to achieve high-performance HOFs. Here, a series of HOFs termed as HOF-ZSTU-M (M=1, 2, and 3) with different pore structures are obtained by introducing structure-directing agents (SDAs) into the hydrogen-bonding network of tetrakis (4-carboxyphenyl) porphyrin (TCPP). These HOFs have distinct space configurations with pore channels ranging from discrete to continuous multi-dimensional. Single-crystal X-ray diffraction (SCXRD) analysis reveals a rare diversity of hydrogen-bonding models dominated by SDAs. HOF-ZSTU-2, which forms a strong layered hydrogen-bonding network with ammonium (NH4 + ) through multiple carboxyl groups, has a suitable 1D "pearl-chain" channel for the selective capture of propylene (C3 H6 ). At 298 K and 1 bar, the C3 H6 storage density of HOF-ZSTU-2 reaches 0.6 kg L-1 , representing one of the best C3 H6 storage materials, while offering a propylene/propane (C3 H6 /C3 H8 ) selectivity of 12.2. Theoretical calculations and in situ SCXRD provide a detailed analysis of the binding strength of C3 H6 at different locations in the pearl-chain channel. Dynamic breakthrough tests confirm that HOF-ZSTU-2 can effectively separate C3 H6 from multi-mixtures.

From Supramolecular Organic Cages to Porous Covalent Organic Frameworks for Enhancing Iodine Adsorption Capability by Fully Exposed Nitrogen-Rich Sites

In order to overcome the limitations of supramolecular organic cages for their incomplete accessibility of active sites in the solid state and uneasy recyclability in liquid solution, herein a nitrogen-rich organic cage is rationally linked into framework systems and four isoreticular covalent organic frameworks (COFs), that is, Cage-TFB-COF, Cage-NTBA-COF, Cage-TFPB-COF, and Cage-TFPT-COF, are successfully synthesized. Structure determination reveals that they are all high-quality crystalline materials derived from the eclipsed packing of related isoreticular two-dimensional frameworks. Since the nitrogen-rich sites usually have a high affinity toward iodine species, iodine adsorption investigations are carried out and the results show that all of them display an enhancement in iodine adsorption capacities. Especially, Cage-NTBA-COF exhibits an iodine adsorption capacity of 304 wt%, 14-fold higher than the solid sample packed from the cage itself. The strong interactions between the nitrogen-rich sites and the adsorbed iodine species are revealed by spectral analyses. This work demonstrates that, utilizing the reticular chemistry strategy to extend the close-packed supramolecular organic cages into crystalline porous framework solids, their inherent properties can be greatly exploited for targeted applications.

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.

Thermally Crosslinked Hydrogen-Bonded Organic Framework Membranes for Highly Selective Ion Separation

The weak bonding energy and flexibility of hydrogen bonds can hinder the long-term use of hydrogen-bonded organic framework (HOF) materials under harsh conditions. Here we invented a thermal-crosslinking method to form polymer materials based on a diamino triazine (DAT) HOF (FDU-HOF-1), containing high-density hydrogen bonding of N-H⋯N. With the increase of temperature to 648 K, the formation of -NH- bonds between neighboring HOF tectons by releasing NH₃ was observed based on the disappearance of the characteristic peaks of amino groups on FDU-HOF-1 in the Fourier transform infrared (FTIR) and solid-state nuclear magnetic resonance (ss-NMR). The variable temperature PXRD indicated the formation of a new peak at 13.2° in addition to the preservation of the original diffraction peaks of FDU-HOF-1. The water adsorption, acid-base stability (12 M HCl to 20 M NaOH) and solubility experiments concluded that the thermally crosslinked HOFs (TC-HOFs) are highly stable. The membranes fabricated by TC-HOF demonstrate the permeation rate of K⁺ ions as high as 270 mmol m^-2 h^-1 as well as high selectivity of K⁺/Mg²⁺ (50) and Na⁺/Mg²⁺ (40), which was comparable to Nafion membranes. This study provides guidance for the future design of highly stable crystalline polymer materials based on HOFs.

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.

Geographic style maps for two-dimensional lattices

This paper develops geographic style maps containing two-dimensional lattices in all known periodic crystals parameterized by recent complete invariants. Motivated by rigid crystal structures, lattices are considered up to rigid motion and uniform scaling. The resulting space of two-dimensional lattices is a square with identified edges or a punctured sphere. The new continuous maps show all Bravais classes as low-dimensional subspaces, visualize hundreds of thousands of lattices of real crystal structures from the Cambridge Structural Database, and motivate the development of continuous and invariant-based crystallography.

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.

Hydrogen-Bonded Organic Frameworks: Structural Design and Emerging Applications

Constructing well-organized organic frameworks with tailor-made functionalities potentially boost multi-domain applications. Hydrogen bonding (H-bonding) is a category of general and weak intermolecular interactions when compared with covalent bonding or metal-ligand coordination. Porous frameworks mainly assembled by H-bonding (named hydrogen-bonded organic frameworks, HOFs) are intrinsically capable of decomposing and regenerating, a distinctive advantage to improve their processability while expanding the applicability. This paper summarizes the basic building concepts of HOFs, including feasible hydrogen bonded motifs, effective molecular structures, and their emerging applications.

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.

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.