Over 14 000 porous, three-dimensional metal–organic framework structures are compiled and analyzed as a part of an update to the Computation-Ready, Experimental Metal–Organic Framework Database (CoRE MOF Database). The updated database includes additional structures that were contributed by CoRE MOF users, obtained from updates of the Cambridge Structural Database and a Web of Science search, and derived through semiautomated reconstruction of disordered structures using a topology-based crystal generator. In addition, value is added to the CoRE MOF database through new analyses that can speed up future nanoporous materials discovery activities, including open metal site detection and duplicate searches. Crystal structures (only for the subset that underwent significant changes during curation), pore analytics, and physical property data are included with the publicly available CoRE MOF 2019 database.

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal–ligand and metal–metal cooperativity, as well as modeling complex catalytic systems such as metal–organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

The catalytic properties of the metal-organic framework Fe2(dobdc), containing open Fe(II) sites, include hydroxylation of phenol by pure Fe2(dobdc) and hydroxylation of ethane by its magnesium-diluted analogue, Fe0.1Mg1.9(dobdc). In earlier work, the latter reaction was proposed to occur through a redox mechanism involving the generation of an iron(IV)-oxo species, which is an intermediate that is also observed or postulated (depending on the case) in some heme and nonheme enzymes and their model complexes. In the present work, we present a detailed mechanism by which the catalytic material, Fe0.1Mg1.9(dobdc), activates the strong C-H bonds of ethane. Kohn-Sham density functional and multireference wave function calculations have been performed to characterize the electronic structure of key species. We show that the catalytic nonheme-Fe hydroxylation of the strong C-H bond of ethane proceeds by a quintet single-state σ-attack pathway after the formation of highly reactive iron-oxo intermediate. The mechanistic pathway involves three key transition states, with the highest activation barrier for the transfer of oxygen from N2O to the Fe(II) center. The uncatalyzed reaction, where nitrous oxide directly oxidizes ethane to ethanol is found to have an activation barrier of 280 kJ/mol, in contrast to 82 kJ/mol for the slowest step in the iron(IV)-oxo catalytic mechanism. The energetics of the C-H bond activation steps of ethane and methane are also compared. Dehydrogenation and dissociation pathways that can compete with the formation of ethanol were shown to involve higher barriers than the hydroxylation pathway.

In the garden of dispersion: High-accuracy ab initio calculations are performed to determine the nature of the interactions and the most favorable geometries between CO(2) and heteroaromatic molecules containing nitrogen (see figure). Dispersion forces play a key role in the stabilization of the dimer, because correlation effects represent about 50 % of the total interaction energy. The interactions between carbon dioxide and organic heterocyclic molecules containing nitrogen are studied by using high-accuracy ab initio methods. Various adsorption positions are examined for pyridine. The preferred configuration is an in-plane configuration. An electron donor-electron acceptor (EDA) mechanism between the carbon of CO(2) and the nitrogen of the heterocycle and weak hydrogen bonds stabilize the complex, with important contributions from dispersion and induction forces. Quantitative results of the binding energy of CO(2) to pyridine (C(5)H(5)N), pyrimidine, pyridazine, and pyrazine (C(4)H(4)N(2)), triazine (C(3)H(3)N(3)), imidazole (C(3)H(4)N(2)), tetrazole (CH(2)N(4)), purine (C(5)H(4)N(4)), imidazopyridine (C(6)H(5)N(3)), adenine (C(5)H(5)N(5)), and imidazopyridamine (C(6)H(6)N(4)) for the in-plane configuration are presented. For purine, three different binding sites are examined. An approximate coupled-cluster model including single and double excitations with a perturbative estimation of triple excitations (CCSD(T)) is used for benchmark calculations. The CCSD(T) basis-set limit is approximated from explicitly correlated second-order Møller-Plesset (MP2-F12) calculations in the aug-cc-pVTZ basis in conjunction with contributions from single, double, and triple excitations calculated at the CCSD(T)/6-311++G** level of theory. Extrapolations to the MP2 basis-set limit coincide with the MP2-F12 calculations. The results are interpreted in terms of electrostatic potential maps and electron density redistribution plots. The effectiveness of density functional theory with the empirical dispersion correction of Grimme (DFT-D) is also examined.

A new large-scale parallel multiconfigurational self-consistent field (MCSCF) implementation in the open-source NWChem computational chemistry code is presented. The generalized active space approach is used to partition large configuration interaction (CI) vectors and generate a sufficient number of batches that can be distributed to the available cores. Massively parallel CI calculations with large active spaces can be performed. The new parallel MCSCF implementation is tested for the chromium trimer and for an active space of 20 electrons in 20 orbitals, which can now routinely be performed. Unprecedented CI calculations with an active space of 22 electrons in 22 orbitals for the pentacene systems were performed and a single CI iteration calculation with an active space of 24 electrons in 24 orbitals for the chromium tetramer was possible. The chromium tetramer corresponds to a CI expansion of one trillion Slater determinants (914 058 513 424) and is the largest conventional CI calculation attempted up to date.

Utilization of self‐healing chemistry to develop synthetic polymer materials that can heal themselves with restored mechanical performance and functionality is of great interest. Self‐healable polymer elastomers with tunable mechanical properties are especially attractive for a variety of applications. Herein, a series of urea functionalized poly(dimethyl siloxane)‐based elastomers (U‐PDMS‐Es) are reported with extremely high stretchability, self‐healing mechanical properties, and recoverable gas‐separation performance. Tailoring the molecular weights of poly(dimethyl siloxane) or weight ratio of elastic cross‐linker offers tunable mechanical properties of the obtained U‐PDMS‐Es, such as ultimate elongation (from 984% to 5600%), Young's modulus, ultimate tensile strength, toughness, and elastic recovery. The U‐PDMS‐Es can serve as excellent acoustic and vibration damping materials over a broad range of temperature (over 100 °C). The strain‐dependent elastic recovery behavior of U‐PDMS‐Es is also studied. After mechanical damage, the U‐PDMS‐Es can be healed in 120 min at ambient temperature or in 20 min at 40 °C with completely restored mechanical performance. The U‐PDMS‐Es are also demonstrated to exhibit recoverable gas‐separation functionality with retained permeability/selectivity after being damaged.

A dicobalt complex catalyzes N2 silylation with Me3SiCl and KC8 under 1 atm N2 at ambient temperature. Tris(trimethylsilyl)amine is formed with an initial turnover rate of 1 N(TMS)3/min, ultimately reaching a turnover number of ∼200. The dicobalt species features a metal-metal interaction, which we postulate is important to its function. Although N2 functionalization occurs at a single cobalt site, the second cobalt center modifies the electronics at the active site. Density functional calculations reveal that the Co-Co interaction evolves during the catalytic cycle: weakening upon N2 binding, breaking with silylation of the metal-bound N2 and reforming with expulsion of [N2(SiMe3)3](-).