Measuring Shared Electrons in Extended Molecular Systems: Covalent Bonds from Plane-Wave Representation of Wave Function.

In the study of materials and macromolecules by first-principle methods, the bond order is a useful tool to represent molecules, bulk materials and interfaces in terms of simple chemical concepts. Despite the availability of several methods to compute the bond order, most applications have been limited to small systems because a high spatial resolution of the wave function and an all-electron representation of the electron density are typically required. Both limitations are critical for large-scale atomistic calculations, even within approximate density-functional theory (DFT) approaches. In this work, we describe our methodology to quickly compute delocalization indices for all atomic pairs, while keeping the same representation of the wave function used in most compute-intensive DFT calculations on high-performance computing equipment. We describe our implementation into a post-processing tool, designed to work with Quantum ESPRESSO, a popular open-source DFT package. In this way, we recover a description in terms of covalent bonds from a representation of wave function containing no explicit information about atomic types and positions.

Metal-Organic Frameworks Invert Molecular Reactivity: Lewis Acidic Phosphonium Zwitterions Catalyze the Aldol-Tishchenko Reaction.

The influence of metal-organic frameworks (MOFs) as additives is herein described for the reaction of n-alkyl aldehydes in the presence of methylvinylketone and triphenylphosphine. In the absence of a MOF, the expected Morita-Baylis-Hillman product, a ß-hydroxy enone, is observed. In the presence of MOFs with UMCM-1 and MOF-5 topologies, the reaction is selective to Aldol-Tishchenko products, the 1 and 3 n-alkylesters of 2-alkyl-1,3-diols, which is unprecedented in organocatalysis. The (3-oxo-2-butenyl)triphenylphosphonium zwitterion, a commonly known nucleophile, is identified as the catalytic active species. This zwitterion favors nucleophilic character in solution, whereas once confined within the framework, it becomes an electrophile yielding Aldol-Tishchenko selectivity. Computational investigations reveal a structural change in the phosphonium moiety induced by the steric confinement of the framework that makes it accessible and an electrophile.

Mixed-linker UiO-66: structure-property relationships revealed by a combination of high-resolution powder X-ray diffraction and density functional theory calculations.

The use of mixed-linker metal-organic frameworks (MIXMOFs) is one of the most effective strategies to modulate the physical-chemical properties of MOFs without affecting the overall crystal structure. In many instances, MIXMOFs have been recognized as solid solutions, with random distribution of ligands, in agreement with the empirical rule known as Vegard's law. In this work, we have undertaken a study combining high-resolution powder X-ray diffraction (HR-PXRD) and density functional theory (DFT) calculations with the aim of understanding the reasons why UiO-66-based amino- and bromo-functionalized MIXMOFs (MIXUiO-66) undergo cell expansion obeying Vegard's law and how this behaviour is related to their physical-chemical properties. DFT calculations predict that the unit cell in amino-functionalized UiO-66 experiences only minor expansion as a result of steric effects, whereas major modification to the electronic features of the framework leads to weaker metal-linker interaction and consequently to the loss of stability at higher degrees of functionalization. For bromo-functionalized UiO-66, steric repulsion due to the size of bromine yields a large cell expansion, but the electronic features remain very similar to pristine UiO-66, preserving the stability of the framework upon functionalization. MIXUiO-66 obtained by either direct synthesis or by post-synthetic exchange shows Vegard-like behaviour, suggesting that both preparation methods yield solid solutions, but the thermal stability and the textural properties of the post-synthetic exchanged materials do not display a clear dependence on the chemical composition, as observed for the MOFs obtained by direct synthesis.

An interacting quantum atoms analysis of the metal-metal bond in [M2(CO)8]n systems.

The metal-metal interaction in policarbonyl metal clusters remains one of the most challenging and controversial issues in metal-organic chemistry, being at heart of a generalized understanding of chemical bonding and of specific applications of these molecules. In this work, the interacting quantum atoms (IQA) approach is used to study the metal-metal interaction in dimetal polycarbonyl dimers, analyzing bridged (Co2(CO)8)), semibridged ([FeCo(CO)8](-)) and unbridged (Co2(CO)8, [Fe2(CO)8](2-)) clusters. In all systems, a delocalized covalent bond is found to occur, involving the metals and the carbonyls, but the global stability of the dimers mainly originates from the Coulombic attraction between the metals and the oxygens.

Influence of the exchange-correlation functional on the quasi-harmonic lattice dynamics of II-VI semiconductors.

We perform a systematic comparison of the finite-temperature structure and properties of four bulk semiconductors (PbS, PbTe, ZnS, and ZnTe) predicted by eight popular exchange-correlation functionals from quasi-harmonic lattice-dynamics calculations. The performance of the functionals in reproducing the temperature dependence of a number of material properties, including lattice parameters, thermal-expansion coefficients, bulk moduli, heat capacities, and phonon frequencies, is evaluated quantitatively against available experimental data. We find that the phenomenological over- and under-binding characteristics of the local-density approximation and the PW91 and Perdew-Burke-Enzerhof (PBE) generalised-gradient approximation (GGA) functionals, respectively, are exaggerated at finite temperature, whereas the PBEsol GGA shows good general performance across all four systems. The Tao-Perdew-Staroverov-Scuseria (TPSS) and revTPSS meta-GGAs provide relatively small improvements over PBE, with the latter being better suited to calculating structural and dynamical properties, but both are considerably more computationally demanding than the simpler GGAs. The dispersion-corrected PBE-D2 and PBE-D3 functionals perform well in describing the lattice dynamics of the zinc chalcogenides, whereas the lead chalcogenides appear to be challenging for these functionals. These findings show that quasi-harmonic calculations with a suitable functional can predict finite-temperature structure and properties with useful accuracy, and that this technique can serve as a means of evaluating the performance of new functionals in the future.

Tunable trimers: using temperature and pressure to control luminescent emission in gold(I) pyrazolate-based trimers.

A systematic investigation into the relationship between the solid-state luminescence and the intermolecular Au⋅⋅⋅Au interactions in a series of pyrazolate-based gold(I) trimers; tris(µ2 -pyrazolato-N,N')-tri-gold(I) (1), tris(µ2 -3,4,5- trimethylpyrazolato-N,N')-tri-gold(I) (2), tris(µ2 -3-methyl-5-phenylpyrazolato-N,N')-tri-gold(I) (3) and tris(µ2 -3,5-diphenylpyrazolato-N,N')-tri-gold(I) (4) has been carried out using variable temperature and high pressure X-ray crystallography, solid-state emission spectroscopy, Raman spectroscopy and computational techniques. Single-crystal X-ray studies show that there is a significant reduction in the intertrimer Au⋅⋅⋅Au distances both with decreasing temperature and increasing pressure. In the four complexes, the reduction in temperature from 293 to 100 K is accompanied by a reduction in the shortest intermolecular Au⋅⋅⋅Au contacts of between 0.04 and 0.08 Å. The solid-state luminescent emission spectra of 1 and 2 display a red shift with decreasing temperature or increasing pressure. Compound 3 does not emit under ambient conditions but displays increasingly red-shifted luminescence upon cooling or compression. Compound 4 remains emissionless, consistent with the absence of intermolecular Au⋅⋅⋅Au interactions. The largest pressure induced shift in emission is observed in 2 with a red shift of approximately 630 cm(-1) per GPa between ambient and 3.80 GPa. The shifts in all the complexes can be correlated with changes in Au⋅⋅⋅Au distance observed by diffraction.

Computational screening of structural and compositional factors for electrically conductive coordination polymers.

The combination of organic and inorganic chemical building blocks to form metal-organic frameworks (MOFs) offers opportunities for producing functional materials suitable for energy generation, storage and conversion. However, such applications rely on robust electron transport and the design of conductive hybrid materials is still in its infancy. Here we apply density functional theory to assess the important structural and compositional factors for forming conducting MOFs. We focus on 1D metal-organic polymers as a model system and assess the choice of organic, inorganic and linking units. The results demonstrate that electronic communication is sensitive to the energy and symmetry of the frontier orbitals associated with the organic and inorganic building blocks and offers guidance on how to optimise electrical conduction in hybrid materials.

Synthesis and characterisation of antimicrobial metal-organic frameworks as multi-drug carriers.

Antibiotic resistance is a significant global concern, necessitating the development of either new antibiotics or advanced delivery methods. With this in mind, we report on the synthesis and characterisation of a new family of Metal-Organic Frameworks (MOFs), OnG6 MOFs, designed to act as multi-drug carriers for bacterial infection treatment. OnG6 is based on the pro-drug 4,4'-azodisalicylic acid (AZDH4), which in vivo produces two equivalents of para-aminosalicylic acid (ASA), a crucial drug for M. tuberculosis treatment. X-ray and computational studies revealed that OnG6 MOFs are mesoporous MOFs with etb topology and an [M2(AZD)] formula (M = Zn, OnG6-Zn; Mg, OnG6-Mg; Cu, OnG6-Cu; and Co, OnG6-Co), featuring 1-dimensional channel type pores of 25 Å diameter. OnG6 MOFs are the first reported MOFs bearing the ligand AZDH4, joining the family of mesoporous MOFs arranged in a honeycomb pattern. They absorb isoniazid (INH) and ciprofloxacin (CIPRO) with the former being a specific antibiotic for M. tuberculosis, and the latter being a broader-spectrum antibiotic. The stability of the MOFs and their capacity for antibiotic uptake depend on the nature of the metal ion, with OnG6-Mg demonstrating the highest drug absorption. The antimicrobial activity of these species was assessed against S. aureus and E. coli, revealing that the carriers containing CIPRO displayed optimal efficacy.

Engineering the optical response of the titanium-MIL-125 metal-organic framework through ligand functionalization.

Herein we discuss band gap modification of MIL-125, a TiO2/1,4-benzenedicarboxylate (bdc) metal-organic framework (MOF). Through a combination of synthesis and computation, we elucidated the electronic structure of MIL-125 with aminated linkers. The band gap decrease observed when the monoaminated bdc-NH2 linker was used arises from donation of the N 2p electrons to the aromatic linking unit, resulting in a red-shifted band above the valence-band edge of MIL-125. We further explored in silico MIL-125 with the diaminated linker bdc-(NH2)2 and other functional groups (-OH, -CH3, -Cl) as alternative substitutions to control the optical response. The bdc-(NH2)2 linking unit was predicted to lower the band gap of MIL-125 to 1.28 eV, and this was confirmed through the targeted synthesis of the bdc-(NH2)2-based MIL-125. This study illustrates the possibility of tuning the optical response of MOFs through rational functionalization of the linking unit, and the strength of combined synthetic/computational approaches for targeting functionalized hybrid materials.

Exploring the Crystal Structure Landscape of Sulfasalazine through Various Multicomponent Crystals.

Sulfasalazine is used as an anti-inflammatory drug to treat large intestine diseases and atrophic arthritis. In the solid state, two tautomers are known: an amide tautomer (triclinic polymorph) and an imide tautomer (monoclinic polymorph). Crystallization of six new multicomponent solids of sulfasalazine with three cocrystal formers and three salt formers has been achieved by slurry, liquid-assisted grinding and slow evaporation methods. All of the solid forms are characterized by X-ray diffraction techniques, thermal analysis, and Fourier transform infrared spectroscopy. The crystal structural analysis reveals that two sulfasalazine molecules or anions arrange in a head-to-head fashion involving their pyridyl, amide, and sulfonyl groups in an R22(7):R22(8):R22(7) motif. This is the key structural unit appearing in both sulfasalazine imide polymorph and all six multicomponent crystals. In addition, sulfasalazine exists in the amide form in all unsolvated multicomponent crystals obtained in this work and adopts the imide tautomer in the solvated cocrystals and salt. Hirshfeld surface analysis and the associated two-dimensional (2D) fingerprint plots demonstrate that sulfasalazine has significant hydrogen bond donor capability when cocrystallized and is a significant hydrogen bond acceptor in the salts. The frontier molecular orbital analysis indicates that sulfasalazine cocrystals are chemically more stable than the salts.

Experimental and Theoretical Investigation of Hydrogen-Bonding Interactions in Cocrystals of Sulfaguanidine.

Pharmaceutical cocrystals, a type of multicomponent crystalline material incorporating two or more molecular and/or ionic compounds connected by noncovalent interactions (such as hydrogen bonds, π-π interactions, and halogen bonds), are attracting increasing attention in crystal engineering. Sulfaguanidine (SGD), one of the most frequently used sulfonamide compounds, was chosen as a model compound in this work to further investigate the hydrogen bond interactions in cocrystals, since it possesses various hydrogen bond donor and acceptor sites. Five cocrystals of SGD, synthesized successfully by slurry and slow evaporation methods, were fully characterized by thermal analysis, X-ray techniques, and Fourier transform infrared spectroscopy. To gain insight into the nature of hydrogen-bonding interactions, theoretical calculations including the analysis of Hirshfeld surface, MEPS (molecular electrostatic potential surface), and QTAIM (quantum theory of atoms in molecules) were conducted. The results are a part of a systematic study of cocrystals of sulfonamides that aims to establish synthon hierarchies in cocrystals containing molecules with multiple hydrogen-bonding functional groups.

Conductive metal-organic frameworks and networks: fact or fantasy?

Electrical conduction is well understood in materials formed from inorganic or organic building blocks, but their combination to produce conductive hybrid frameworks and networks is an emerging and rapidly developing field of research. Self-assembling organic-inorganic compounds offer immense potential for functionalising material properties for a wide scope of applications including solar cells, light emitters, gas sensors and bipolar transparent conductors. The flexibility of combining two distinct material classes into a single solid-state system provides an almost infinite number of chemical and structural possibilities; however, there is currently no systematic approach established for designing new compositions and configurations with targeted electronic or optical properties. We review the current status in the field, in particular, the range of hybrid systems reported to date and the important role of materials modelling in the field. From theoretical arguments, the Mott insulator-to-metal transition should be possible in semiconducting metal-organic frameworks, but has yet to be observed. The question remains as to whether electro-active hybrid materials will evolve from chemical curiosities towards practical applications in the near term.

Restoring orbital thinking from real space descriptions: bonding in classical and non-classical transition metal carbonyls.

A combined strategy that unifies our interacting quantum atoms approach (IQA), a chemically intuitive energetic perspective within the quantum theory of atoms in molecules (QTAIM), the domain natural orbitals obtained by the diagonalization of the charge-weighted domain-averaged Fermi hole (DAFH), and the statistical analyses of chemical bonding provided by the electron number distribution functions (EDF) is presented. As shown, it allows for recovering traditional orbital images from the orbital invariant descriptions of QTAIM. It does also provide bonding indices (like bond orders) and bond energetics, all in a per orbital basis, still invariant manner, using a single unified framework. The procedure is applied to show how the Dewar, Chatt, and Ducanson model of bonding in simple transition metal carbonyls may be recovered in the real space. The balance between the number of σ-donated and π-backdonated electrons is negative in classical compounds and positive in non-classical ones. The energetic strength of backdonation is, however, smaller than that of donation. Our technique surpasses conventional orbital models by providing physically sound, quantitative energetics of chemical bonds (or interactions) together with effective one-electron pictures, all for arbitrary wavefunctions.

Metal-organic frameworks as kinetic modulators for branched selectivity in hydroformylation.

Finding heterogeneous catalysts that are superior to homogeneous ones for selective catalytic transformations is a major challenge in catalysis. Here, we show how micropores in metal-organic frameworks (MOFs) push homogeneous catalytic reactions into kinetic regimes inaccessible under standard conditions. Such property allows branched selectivity up to 90% in the Co-catalysed hydroformylation of olefins without directing groups, not achievable with existing catalysts. This finding has a big potential in the production of aldehydes for the fine chemical industry. Monte Carlo and density functional theory simulations combined with kinetic models show that the micropores of MOFs with UMCM-1 and MOF-74 topologies increase the olefins density beyond neat conditions while partially preventing the adsorption of syngas leading to high branched selectivity. The easy experimental protocol and the chemical and structural flexibility of MOFs will attract the interest of the fine chemical industries towards the design of heterogeneous processes with exceptional selectivity.

Taking lanthanides out of isolation: tuning the optical properties of metal-organic frameworks.

Metal organic frameworks (MOFs) are increasingly used in applications that rely on the optical and electronic properties of these materials. These applications require a fundamental understanding on how the structure of these materials, and in particular the electronic interactions of the metal node and organic linker, determines these properties. Herein, we report a combined experimental and computational study on two families of lanthanide-based MOFs: Ln-SION-1 and Ln-SION-2. Both comprise the same metal and ligand but with differing structural topologies. In the Ln-SION-2 series the optical absorption is dominated by the ligand and using different lanthanides has no impact on the absorption spectrum. The Ln-SION-1 series shows a completely different behavior in which the ligand and the metal node do interact electronically. By changing the lanthanide in Ln-SION-1, we were able to tune the optical absorption from the UV region to absorption that includes a large part of the visible region. For the early lanthanides we observe intraligand (electronic) transitions in the UV region, while for the late lanthanides a new band appears in the visible. DFT calculations showed that the new band in the visible originates in the spatial orbital overlap between the ligand and metal node. Our quantum calculations indicated that Ln-SION-1 with late lanthanides might be (photo)conductive. Experimentally, we confirm that these materials are weakly conductive and that with an appropriate co-catalysts they can generate hydrogen from a water solution using visible light. Our experimental and theoretical analysis provides fundamental insights for the rational design of Ln-MOFs with the desired optical and electronic properties.

Using pseudopotentials within the interacting quantum atoms approach.

A general strategy to extend the interacting quantum atoms (IQA) approach to pseudopotential or effective core potential electronic structure calculations is presented. With the protocol proposed here, the scope of IQA thinking opens to chemical bonding problems in heavy-atom systems, as well as to larger molecules than those presently allowed by computational limitations. We show that, provided that interatomic surfaces are computed from core-reconstructed densities, reasonable results are obtained by integrating reduced density matrices built from the pseudowave functions. Comparison with all-electron results in a few test systems shows that exchange-correlation energies are better reproduced than Coulombic contributions, an effect which is traced to inadequate atomic populations and leakage of the core population into the surrounding quantum atoms.

Rapid, Selective Heavy Metal Removal from Water by a Metal-Organic Framework/Polydopamine Composite.

Drinking water contamination with heavy metals, particularly lead, is a persistent problem worldwide with grave public health consequences. Existing purification methods often cannot address this problem quickly and economically. Here we report a cheap, water stable metal-organic framework/polymer composite, Fe-BTC/PDA, that exhibits rapid, selective removal of large quantities of heavy metals, such as Pb2+ and Hg2+, from real world water samples. In this work, Fe-BTC is treated with dopamine, which undergoes a spontaneous polymerization to polydopamine (PDA) within its pores via the Fe3+ open metal sites. The PDA, pinned on the internal MOF surface, gains extrinsic porosity, resulting in a composite that binds up to 1634 mg of Hg2+ and 394 mg of Pb2+ per gram of composite and removes more than 99.8% of these ions from a 1 ppm solution, yielding drinkable levels in seconds. Further, the composite properties are well-maintained in river and seawater samples spiked with only trace amounts of lead, illustrating unprecedented selectivity. Remarkably, no significant uptake of competing metal ions is observed even when interferents, such as Na+, are present at concentrations up to 14 000 times that of Pb2+. The material is further shown to be resistant to fouling when tested in high concentrations of common organic interferents, like humic acid, and is fully regenerable over many cycles.

Origin of the Strong Interaction between Polar Molecules and Copper(II) Paddle-Wheels in Metal Organic Frameworks.

The copper paddle-wheel is the building unit of many metal organic frameworks. Because of the ability of the copper cations to attract polar molecules, copper paddle-wheels are promising for carbon dioxide adsorption and separation. They have therefore been studied extensively, both experimentally and computationally. In this work we investigate the copper-CO2 interaction in HKUST-1 and in two different cluster models of HKUST-1: monocopper Cu(formate)2 and dicopper Cu2(formate)4. We show that density functional theory methods severely underestimate the interaction energy between copper paddle-wheels and CO2, even including corrections for the dispersion forces. In contrast, a multireference wave function followed by perturbation theory to second order using the CASPT2 method correctly describes this interaction. The restricted open-shell Møller-Plesset 2 method (ROS-MP2, equivalent to (2,2) CASPT2) was also found to be adequate in describing the system and used to develop a novel force field. Our parametrization is able to predict the experimental CO2 adsorption isotherms in HKUST-1, and it is shown to be transferable to other copper paddle-wheel systems.

Adsorbate-induced lattice deformation in IRMOF-74 series.

IRMOF-74 analogues are among the most widely studied metal-organic frameworks (MOFs) for adsorption applications because of their one-dimensional channels and high metal density. Most studies involving the IRMOF-74 series assume that the crystal lattice is rigid. This assumption guides the interpretation of experimental data, as changes in the crystal symmetry have so far been ignored as a possibility in the literature. Here, we report a deformation pattern, induced by the adsorption of argon, for IRMOF-74-V. This work has two main implications. First, we use molecular simulations to demonstrate that the IRMOF-74 series undergoes a deformation that is similar to the mechanism behind breathing MOFs, but is unique because the deformation pattern extends beyond a single unit cell of the original structure. Second, we provide an alternative interpretation of experimental small-angle X-ray scattering profiles of these systems, which changes how we view the fundamentals of adsorption in this MOF series.

Correction: A universal chemical potential for sulfur vapours.

Correction for 'A universal chemical potential for sulfur vapours' by Adam J. Jackson et al., Chem. Sci., 2016, 7, 1082-1092.