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.

The air-free reaction of CoCl2 with 1,3,5-tri(1H-1,2,3-triazol-5-yl)benzene (H3BTTri) in N,N-dimethylformamide (DMF) and methanol leads to the formation of Co-BTTri (Co3[(Co4Cl)3(BTTri)8]2·DMF), a sodalite-type metal-organic framework. Desolvation of this material generates coordinatively unsaturated low-spin cobalt(II) centers that exhibit a strong preference for binding O2 over N2, with isosteric heats of adsorption (Qst) of -34(1) and -12(1) kJ/mol, respectively. The low-spin (S = 1/2) electronic configuration of the metal centers in the desolvated framework is supported by structural, magnetic susceptibility, and computational studies. A single-crystal X-ray structure determination reveals that O2 binds end-on to each framework cobalt center in a 1:1 ratio with a Co-O2 bond distance of 1.973(6) Å. Replacement of one of the triazolate linkers with a more electron-donating pyrazolate group leads to the isostructural framework Co-BDTriP (Co3[(Co4Cl)3(BDTriP)8]2·DMF; H3BDTriP = 5,5'-(5-(1H-pyrazol-4-yl)-1,3-phenylene)bis(1H-1,2,3-triazole)), which demonstrates markedly higher yet still fully reversible O2 affinities (Qst = -47(1) kJ/mol at low loadings). Electronic structure calculations suggest that the O2 adducts in Co-BTTri are best described as cobalt(II)-dioxygen species with partial electron transfer, while the stronger binding sites in Co-BDTriP form cobalt(III)-superoxo moieties. The stability, selectivity, and high O2 adsorption capacity of these materials render them promising new adsorbents for air separation processes.

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.