Laura Gagliardi

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas–Kroll–Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC‐PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large‐scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Cooperative insertion of CO2 in diamine-appended metal-organic frameworks

The process of carbon capture and sequestration has been proposed as a method of mitigating the build-up of greenhouse gases in the atmosphere. If implemented, the cost of electricity generated by a fossil fuel-burning power plant would rise substantially, owing to the expense of removing CO2 from the effluent stream. There is therefore an urgent need for more efficient gas separation technologies, such as those potentially offered by advanced solid adsorbents. Here we show that diamine-appended metal-organic frameworks can behave as ‘phase-change’ adsorbents, with unusual step-shaped CO2 adsorption isotherms that shift markedly with temperature. Results from spectroscopic, diffraction and computational studies show that the origin of the sharp adsorption step is an unprecedented cooperative process in which, above a metal-dependent threshold pressure, CO2 molecules insert into metal-amine bonds, inducing a reorganization of the amines into well-ordered chains of ammonium carbamate. As a consequence, large CO2 separation capacities can be achieved with small temperature swings, and regeneration energies appreciably lower than achievable with state-of-the-art aqueous amine solutions become feasible. The results provide a mechanistic framework for designing highly efficient adsorbents for removing CO2 from various gas mixtures, and yield insights into the conservation of Mg2+ within the ribulose-1,5-bisphosphate carboxylase/oxygenase family of enzymes.

The restricted active space followed by second-order perturbation theory method: Theory and application to the study of CuO2 and Cu2O2 systems

A multireference second-order perturbation theory using a restricted active space self-consistent field wave function as reference (RASPT2/RASSCF) is described. This model is particularly effective for cases where a chemical system requires a balanced orbital active space that is too large to be addressed by the complete active space self-consistent field model with or without second-order perturbation theory (CASPT2 or CASSCF, respectively). Rather than permitting all possible electronic configurations of the electrons in the active space to appear in the reference wave function, certain orbitals are sequestered into two subspaces that permit a maximum number of occupations or holes, respectively, in any given configuration, thereby reducing the total number of possible configurations. Subsequent second-order perturbation theory captures additional dynamical correlation effects. Applications of the theory to the electronic structure of complexes involved in the activation of molecular oxygen by mono- and binuclear copper complexes are presented. In the mononuclear case, RASPT2 and CASPT2 provide very similar results. In the binuclear cases, however, only RASPT2 proves quantitatively useful, owing to the very large size of the necessary active space.

Quantum chemical calculations show that the uranium molecule U2 has a quintuple bond

Covalent bonding is commonly described by Lewiss theory, with an electron pair shared between two atoms constituting one full bond. Beginning with the valence bond description for the hydrogen molecule, quantum chemists have further explored the fundamental nature of the chemical bond for atoms throughout the periodic table, confirming that most molecules are indeed held together by one electron pair for each bond. But more complex binding may occur when large numbers of atomic orbitals can participate in bond formation. Such behaviour is common with transition metals. When involving heavy actinide elements, metal–metal bonds might prove particularly complicated. To date, evidence for actinide–actinide bonds is restricted to the matrix-isolation of uranium hydrides, including H2U–UH2, and the gas-phase detection and preliminary theoretical study of the uranium molecule, U2. Here we report quantum chemical calculations on U2, showing that, although the strength of the U2 bond is comparable to that of other multiple bonds between transition metals, the bonding pattern is unique. We find that the molecule contains three electron-pair bonds and four one-electron bonds (that is, 10 bonding electrons, corresponding to a quintuple bond), and two ferromagnetically coupled electrons localized on one U atom each—so all known covalent bonding types are contributing.

Local properties of quantum chemical systems: The LoProp approach

A new method is presented, which makes it possible to partition molecular properties like multipole moments and polarizabilities, into atomic and interatomic contributions. The method requires a subdivision of the atomic basis set into occupied and virtual basis functions for each atom in the molecular system. The localization procedure is organized into a series of orthogonalizations of the original basis set, which will have as a final result a localized orthonormal basis set. The new localization procedure is demonstrated to be stable with various basis sets, and to provide physically meaningful localized properties. Transferability of the methyl properties for the alkane series and of the carbon and hydrogen properties for the benzene, naphtalene, and anthracene series is demonstrated.

Oxidation of ethane to ethanol by N2O in a metal–organic framework with coordinatively unsaturated iron(II) sites

Enzymatic haem and non-haem high-valent iron-oxo species are known to activate strong C-H bonds, yet duplicating this reactivity in a synthetic system remains a formidable challenge. Although instability of the terminal iron-oxo moiety is perhaps the foremost obstacle, steric and electronic factors also limit the activity of previously reported mononuclear iron(IV)-oxo compounds. In particular, although natures non-haem iron(IV)-oxo compounds possess high-spin S = 2 ground states, this electronic configuration has proved difficult to achieve in a molecular species. These challenges may be mitigated within metal-organic frameworks that feature site-isolated iron centres in a constrained, weak-field ligand environment. Here, we show that the metal-organic framework Fe2(dobdc) (dobdc(4-) = 2,5-dioxido-1,4-benzenedicarboxylate) and its magnesium-diluted analogue, Fe0.1Mg1.9(dobdc), are able to activate the C-H bonds of ethane and convert it into ethanol and acetaldehyde using nitrous oxide as the terminal oxidant. Electronic structure calculations indicate that the active oxidant is likely to be a high-spin S = 2 iron(IV)-oxo species.

Atomic Cholesky decompositions: A route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency

Cholesky decomposition of the atomic two-electron integral matrix has recently been proposed as a procedure for automated generation of auxiliary basis sets for the density fitting approximation [F. Aquilante et al., J. Chem. Phys. 127, 114107 (2007)]. In order to increase computational performance while maintaining accuracy, we propose here to reduce the number of primitive Gaussian functions of the contracted auxiliary basis functions by means of a second Cholesky decomposition. Test calculations show that this procedure is most beneficial in conjunction with highly contracted atomic orbital basis sets such as atomic natural orbitals, and that the error resulting from the second decomposition is negligible. We also demonstrate theoretically as well as computationally that the locality of the fitting coefficients can be controlled by means of the decomposition threshold even with the long-ranged Coulomb metric. Cholesky decomposition-based auxiliary basis sets are thus ideally suited for local density fitting approximations.

Metal-Organic Framework Nodes as Nearly Ideal Supports for Molecular Catalysts: NU-1000- and UiO-66-Supported Iridium Complexes

Metal-organic frameworks with Zr6 nodes, UiO-66 and NU-1000, were investigated as supports for Ir(CO)2 and Ir(C2H4)2 complexes. A single bonding site for the iridium is identified on the nodes of NU-1000, whereas two sites are identified on UiO-66, although at low iridium loadings only one site is occupied. Density functional theory calculations provide structural results that are in good agreement with infrared and X-ray absorption fine-structure spectra. The reactivity of node-supported Ir(CO)2 with C2H4 and the catalytic activity and selectivity of the species initially present as Ir(C2H4)2 for ethylene hydrogenation and dimerization were investigated both experimentally and computationally and shown to be strongly influenced by the node.

The generalized active space concept in multiconfigurational self-consistent field methods.

A multiconfigurational self-consistent field method based on the concept of generalized active space (GAS) is presented. GAS wave functions are obtained by defining an arbitrary number of active spaces with arbitrary occupation constraints. By a suitable choice of the GAS spaces, numerous ineffective configurations present in a large complete active space (CAS) can be removed, while keeping the important ones in the CI space. As a consequence, the GAS self-consistent field approach retains the accuracy of the CAS self-consistent field (CASSCF) ansatz and, at the same time, can deal with larger active spaces, which would be unaffordable at the CASSCF level. Test calculations on the Gd atom, Gd(2) molecule, and oxoMn(salen) complex are presented. They show that GAS wave functions achieve the same accuracy as CAS wave functions on systems that would be prohibitive at the CAS level.

Sintering-resistant Single-Site Nickel Catalyst Supported by Metal-Organic Framework

Developing supported single-site catalysts is an important goal in heterogeneous catalysis since the well-defined active sites afford opportunities for detailed mechanistic studies, thereby facilitating the design of improved catalysts. We present herein a method for installing Ni ions uniformly and precisely on the node of a Zr-based metal-organic framework (MOF), NU-1000, in high density and large quantity (denoted as Ni-AIM) using atomic layer deposition (ALD) in a MOF (AIM). Ni-AIM is demonstrated to be an efficient gas-phase hydrogenation catalyst upon activation. The structure of the active sites in Ni-AIM is proposed, revealing its single-site nature. More importantly, due to the organic linker used to construct the MOF support, the Ni ions stay isolated throughout the hydrogenation catalysis, in accord with its long-term stability. A quantum chemical characterization of the catalyst and the catalytic process complements the experimental results. With validation of computational modeling protocols, we further targeted ethylene oligomerization catalysis by Ni-AIM guided by theoretical prediction. Given the generality of the AIM methodology, this emerging class of materials should prove ripe for the discovery of new catalysts for the transformation of volatile substrates.