Solution-processable polytriazoles from spirocyclic monomers for membrane-based hydrocarbon separations.

The thermal distillation of crude oil mixtures is an energy-intensive process, accounting for nearly 1% of global energy consumption. Membrane-based separations are an appealing alternative or tandem process to distillation due to intrinsic energy efficiency advantages. We developed a family of spirocyclic polytriazoles from structurally diverse monomers for membrane applications. The resulting polymers were prepared by a convenient step-growth method using copper-catalysed azide-alkyne cycloaddition, providing very fast reaction rates, high molecular weights and solubilities in common organic solvents and non-interconnected microporosity. Fractionation of whole Arabian light crude oil and atmospheric tower bottom feeds using these materials enriched the low-boiling-point components and removed trace heteroatom and metal impurities (comparable performance with the lighter feed as the commercial polyimide, Matrimid), demonstrating opportunities to reduce the energy cost of crude oil distillation with tandem membrane processes. Membrane-based molecular separation under these demanding conditions is made possible by high thermal stability and a moderate level of dynamic chain mobility, leading to transient interconnections between micropores, as revealed by the calculations of static and swollen pore structures.

A Combined Experimental and Computational Exploration of Heteroleptic cis-Pd2L2L'2 Coordination Cages through Geometric Complementarity.

Heteroleptic (mixed-ligand) coordination cages are of interest as host systems with more structurally and functionally complex cavities than homoleptic architectures. The design of heteroleptic cages, however, is far from trivial. In this work, we experimentally probed the self-assembly of Pd(II) ions with binary ligand combinations in a combinatorial fashion to search for new cis-Pd2L2L'2 heteroleptic cages. A hierarchy of computational analyses was then applied to these systems with the aim of elucidating key factors for rationalising self-assembly outcomes. Simple and inexpensive geometric analyses were shown to be effective in identifying complementary ligand pairs. Preliminary results demonstrated the viability of relatively rapid semi-empirical calculations for predicting the topology of thermodynamically favoured assemblies with rigid ligands, whilst more flexible systems proved challenging. Stemming from this, key challenges were identified for future work developing effective computational forecasting tools for self-assembled metallo-supramolecular systems.

Diastereoselective Self-Assembly of Low-Symmetry Pdn L2n Nanocages through Coordination-Sphere Engineering.

Metal-organic cages (MOCs) are popular host architectures assembled from ligands and metal ions/nodes. Assembling structurally complex, low-symmetry MOCs with anisotropic cavities can be limited by the formation of statistical isomer libraries. We set out to investigate the use of primary coordination-sphere engineering (CSE) to bias isomer selectivity within homo- and heteroleptic Pdn L2n cages. Unexpected differences in selectivities between alternative donor groups led us to recognise the significant impact of the second coordination sphere on isomer stabilities. From this, molecular-level insight into the origins of selectivity between cis and trans diastereoisomers was gained, highlighting the importance of both host-guest and host-solvent interactions, in addition to ligand design. This detailed understanding allows precision engineering of low-symmetry MOC assemblies without wholesale redesign of the ligand framework, and fundamentally provides a theoretical scaffold for the development of stimuli-responsive, shape-shifting MOCs.

Tetramine Aspect Ratio and Flexibility Determine Framework Symmetry for Zn8 L6 Self-Assembled Structures.

We derive design principles for the assembly of rectangular tetramines into Zn8 L6 pseudo-cubic coordination cages. Because of the rectangular, as opposed to square, geometry of the ligand panels, and the possibility of either Δ or Λ handedness of each metal center at the eight corners of the pseudo-cube, many different cage diastereomers are possible. Each of the six tetra-aniline subcomponents investigated in this work assembled with zinc(II) and 2-formylpyridine in acetonitrile into a single Zn8 L6 pseudo-cube diastereomer, however. Each product corresponded to one of four diastereomeric configurations, with T, Th , S6 or D3 symmetry. The preferred diastereomer for a given tetra-aniline subcomponent was shown to be dependent on its aspect ratio and conformational flexibility. Analysis of computationally modeled individual faces or whole pseudo-cubes provided insight as to why the observed diastereomers were favored.

Into the Unknown: How Computation Can Help Explore Uncharted Material Space.

Novel functional materials are urgently needed to help combat the major global challenges facing humanity, such as climate change and resource scarcity. Yet, the traditional experimental materials discovery process is slow and the material space at our disposal is too vast to effectively explore using intuition-guided experimentation alone. Most experimental materials discovery programs necessarily focus on exploring the local space of known materials, so we are not fully exploiting the enormous potential material space, where more novel materials with unique properties may exist. Computation, facilitated by improvements in open-source software and databases, as well as computer hardware has the potential to significantly accelerate the rational development of materials, but all too often is only used to postrationalize experimental observations. Thus, the true predictive power of computation, where theory leads experimentation, is not fully utilized. Here, we discuss the challenges to successful implementation of computation-driven materials discovery workflows, and then focus on the progress of the field, with a particular emphasis on the challenges to reaching novel materials.

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 (1C/C , 1Si/C and 1Si/Si ) are studied. Three distinct phases for 1C/C and four for 1Si/C and 1Si/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 1C/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.

Sterics and Hydrogen Bonding Control Stereochemistry and Self-Sorting in BINOL-Based Assemblies.

Here we demonstrate how the hydrogen-bonding ability of a BINOL-based dialdehyde subcomponent dictated the stereochemical outcome of its subsequent self-assembly into one diastereomeric helicate form when bearing free hydroxy groups, and another in the case of its methylated congener. The presence of methyl groups also altered the self-sorting behavior when mixed with another, short linear dialdehyde subcomponent, switching the outcome of the system from narcissistic to integrative self-sorting. In all cases, the axial chirality of the BINOL building block also dictated helicate metal center handedness during stereospecific self-assembly. A new family of stereochemically pure heteroleptic helicates were thus prepared using the new knowledge gained. We also found that switching from FeII to ZnII, or the incorporation of a longer linear ligand, favored heteroleptic structure formation.

stk: An extendable Python framework for automated molecular and supramolecular structure assembly and discovery.

Computational software workflows are emerging as all-in-one solutions to speed up the discovery of new materials. Many computational approaches require the generation of realistic structural models for property prediction and candidate screening. However, molecular and supramolecular materials represent classes of materials with many potential applications for which there is no go-to database of existing structures or general protocol for generating structures. Here, we report a new version of the supramolecular toolkit, stk, an open-source, extendable, and modular Python framework for general structure generation of (supra)molecular structures. Our construction approach works on arbitrary building blocks and topologies and minimizes the input required from the user, making stk user-friendly and applicable to many material classes. This version of stk includes metal-containing structures and rotaxanes as well as general implementation and interface improvements. Additionally, this version includes built-in tools for exploring chemical space with an evolutionary algorithm and tools for database generation and visualization. The latest version of stk is freely available at github.com/lukasturcani/stk.

High-Throughput Computational Evaluation of Low Symmetry Pd2 L4 Cages to Aid in System Design*.

Unsymmetrical ditopic ligands can self-assemble into reduced-symmetry Pd2 L4 metallo-cages with anisotropic cavities, with implications for high specificity and affinity guest-binding. Mixtures of cage isomers can form, however, resulting in undesirable system heterogeneity. It is paramount to be able to design components that preferentially form a single isomer. Previous data suggested that computational methods could predict with reasonable accuracy whether unsymmetrical ligands would preferentially self-assemble into single cage isomers under constraints of geometrical mismatch. We successfully apply a collaborative computational and experimental workflow to mitigate costly trial-and-error synthetic approaches. Our rapid computational workflow constructs unsymmetrical ligands and their Pd2 L4 cage isomers, ranking the likelihood for exclusively forming cis-Pd2 L4 assemblies. From this narrowed search space, we successfully synthesised four new, low-symmetry, cis-Pd2 L4 cages.

Enhanced Activity of Enzymes Encapsulated in Hydrophilic Metal-Organic Frameworks.

Encapsulation of biomacromolecules in metal-organic frameworks (MOFs) can preserve biological functionality in harsh environments. Despite the success of this approach, termed biomimietic mineralization, limited consideration has been given to the chemistry of the MOF coating. Here, we show that enzymes encapsulated within hydrophilic MAF-7 or ZIF-90 retain enzymatic activity upon encapsulation and when exposed to high temperatures, denaturing or proteolytic agents, and organic solvents, whereas hydrophobic ZIF-8 affords inactive catalase and negligible protection to urease.

Molecular Tectonics: A Node-and-Linker Building Block Approach to a Family of Hydrogen-Bonded Frameworks.

While numerous hydrogen-bonded organic frameworks (HOFs) have been reported, typically these cannot be prepared predictably or in a modular fashion. In this work, we report a family of nine diamondoid crystalline porous frameworks assembled via hydrogen bonding between poly-amidinium and poly-carboxylate tectons. The frameworks are prepared at room temperature in either water or water/alcohol mixtures. Importantly, both the cationic and anionic components can be varied and additional functionality can be incorporated into the frameworks, which show good stability including to prolonged heating in DMSO or water.

The effect of disorder in multi-component covalent organic frameworks.

We examined the effect of two different types of linker distribution-random or correlated distribution-on the pore size and shape within single-layers of three multi-component COFs. We reveal a relationship between linker distribution and the porosity of COF solid solutions. The methods presented in this paper are generalisable and could be used in further studies to examine the properties of disordered framework materials.

Systematic exploration of accessible topologies of cage molecules via minimalistic models.

Cages are macrocyclic structures with an intrinsic internal cavity that support applications in separations, sensing and catalysis. These materials can be synthesised via self-assembly of organic or metal-organic building blocks. Their bottom-up synthesis and the diversity in building block chemistry allows for fine-tuning of their shape and properties towards a target property. However, it is not straightforward to predict the outcome of self-assembly, and, thus, the structures that are practically accessible during synthesis. Indeed, such a prediction becomes more difficult as problems related to the flexibility of the building blocks or increased combinatorics lead to a higher level of complexity and increased computational costs. Molecular models, and their coarse-graining into simplified representations, may be very useful to this end. Here, we develop a minimalistic toy model of cage-like molecules to explore the stable space of different cage topologies based on a few fundamental geometric building block parameters. Our results capture, despite the simplifications of the model, known geometrical design rules in synthetic cage molecules and uncover the role of building block coordination number and flexibility on the stability of cage topologies. This leads to a large-scale and systematic exploration of design principles, generating data that we expect could be analysed through expandable approaches towards the rational design of self-assembled porous architectures.

Pore topology analysis in porous molecular systems.

Porous molecular materials are constructed from molecules that assemble in the solid-state such that there are cavities or an interconnected pore network. It is challenging to control the assembly of these systems, as the interactions between the molecules are generally weak, and subtle changes in the molecular structure can lead to vastly different intermolecular interactions and subsequently different crystal packing arrangements. Similarly, the use of different solvents for crystallization, or the introduction of solvent vapour, can result in different polymorphs and pore networks being formed. It is difficult to uniquely describe the pore networks formed, and thus we analyse 1033 crystal structures of porous molecular systems to determine the underlying topology of their void spaces and potential guest diffusion networks. Material-agnostic topology definitions are applied. We use the underlying topological nets to examine whether it is possible to apply isoreticular design principles to porous molecular materials. Overall, our automatic analysis of a large dataset gives a general insight into the relationships between molecular topologies and the topological nets of their pore network. We show that while porous molecular systems tend to pack similarly to non-porous molecules, the topologies of their pore distributions resemble those of more prominent porous materials, such as metal-organic frameworks and covalent organic frameworks.

Competitive aminal formation during the synthesis of a highly soluble, isopropyl-decorated imine porous organic cage.

The synthesis of a new porous organic cage decorated with isopropyl moieties (CC21) was achieved from the reaction of triformylbenzene and an isopropyl functionalised diamine. Unlike structurally analogous porous organic cages, its synthesis proved challenging due to competitive aminal formation, rationalised using control experiments and computational modelling. The use of an additional amine was found to increase conversion to the desired cage.

Unlocking the computational design of metal-organic cages.

Metal-organic cages are macrocyclic structures that can possess an intrinsic void that can hold molecules for encapsulation, adsorption, sensing, and catalysis applications. As metal-organic cages may be comprised from nearly any combination of organic and metal-containing components, cages can form with diverse shapes and sizes, allowing for tuning toward targeted properties. Therefore, their near-infinite design space is almost impossible to explore through experimentation alone and computational design can play a crucial role in exploring new systems. Although high-throughput computational design and screening workflows have long been known as powerful tools in drug and materials discovery, their application in exploring metal-organic cages is more recent. We show examples of structure prediction and host-guest/catalytic property evaluation of metal-organic cages. These examples are facilitated by advances in methods that handle metal-containing systems with improved accuracy and are the beginning of the development of automated cage design workflows. We finally outline a scope for how high-throughput computational methods can assist and drive experimental decisions as the field pushes toward functional and complex metal-organic cages. In particular, we highlight the importance of considering realistic, flexible systems.

Modeling the Effect of Defects and Disorder in Amorphous Metal-Organic Frameworks.

Amorphous metal-organic frameworks (aMOFs) are a class of disordered framework materials with a defined local order given by the connectivity between inorganic nodes and organic linkers, but absent long-range order. The rational development of function for aMOFs is hindered by our limited understanding of the underlying structure-property relationships in these systems, a consequence of the absence of long-range order, which makes experimental characterization particularly challenging. Here, we use a versatile modeling approach to generate in silico structural models for an aMOF based on Fe trimers and 1,3,5-benzenetricarboxylate (BTC) linkers, Fe-BTC. We build a phase space for this material that includes nine amorphous phases with different degrees of defects and local order. These models are analyzed through a combination of structural analysis, pore analysis, and pair distribution functions. Therefore, we are able to systematically explore the effects of the variation of each of these features, both in isolation and combined, for a disordered MOF system, something that would not be possible through experiment alone. We find that the degree of local order has a greater impact on structure and properties than the degree of defects. The approach presented here is versatile and allows for the study of different structural features and MOF chemistries, enabling the derivation of design rules for the rational development of aMOFs.

Orientational self-sorting in cuboctahedral Pd cages.

Cuboctahedral coordination cages of the general formula [Pd12L24]24+ (L = low-symmetry ligand) were analyzed theoretically and experimentally. With 350 696 potential isomers, the structural space of these assemblies is vast. Orientational self-sorting refers to the preferential formation of particular isomers within the pool of potential structures. Geometric and computational analyses predict the preferred formation of cages with a cis arrangement at the metal centers. This prediction was corroborated experimentally by synthesizing a [Pd12L24]24+ cage with a bridging 3-(4-(pyridin-4-yl)phenyl)pyridine ligand. A crystallographic analysis of this assembly showed exclusive cis coordination of the 3- and the 4-pyridyl donor groups at the Pd2+ ions.

N-Aryl-linked spirocyclic polymers for membrane separations of complex hydrocarbon mixtures.

The fractionation of crude-oil mixtures through distillation is a large-scale, energy-intensive process. Membrane materials can avoid phase changes in such mixtures and thereby reduce the energy intensity of these thermal separations. With this application in mind, we created spirocyclic polymers with N-aryl bonds that demonstrated noninterconnected microporosity in the absence of ladder linkages. The resulting glassy polymer membranes demonstrated nonthermal membrane fractionation of light crude oil through a combination of class- and size-based "sorting" of molecules. We observed an enrichment of molecules lighter than 170 daltons corresponding to a carbon number of 12 or a boiling point less than 200°C in the permeate. Such scalable, selective membranes offer potential for the hybridization of energy-efficient technology with conventional processes such as distillation.