We present an automated, open source toolkit for the first-principles screening and discovery of new inorganic molecules and intermolecular complexes. Challenges remain in the automatic generation of candidate inorganic molecule structures due to the high variability in coordination and bonding, which we overcome through a divide-and-conquer tactic that flexibly combines force-field preoptimization of organic fragments with alignment to first-principles-trained metal-ligand distances. Exploration of chemical space is enabled through random generation of ligands and intermolecular complexes from large chemical databases. We validate the generated structures with the root mean squared (RMS) gradients evaluated from density functional theory (DFT), which are around 0.02 Ha/au across a large 150 molecule test set. Comparison of molSimplify results to full optimization with the universal force field reveals that RMS DFT gradients are improved by 40%. Seamless generation of input files, preparation and execution of electronic structure calculations, and post-processing for each generated structure aids interpretation of underlying chemical and energetic trends. © 2016 Wiley Periodicals, Inc.

Articles you may be interested inWe estimate the prediction sensitivity with respect to Hartree-Fock exchange in approximate density functionals for representative Fe(II) and Fe(III) octahedral complexes. Based on the observation that the range of parameters spanned by the most widely employed functionals is relatively narrow, we compute electronic structure property and spin-state orderings across a relatively broad range of Hartree-Fock exchange (0%-50%) ratios. For the entire range considered, we consistently observe linear relationships between spin-state ordering that differ only based on the element of the direct ligand and thus may be broadly employed as measures of functional sensitivity in predictions of organometallic compounds. The role Hartree-Fock exchange in hybrid functionals is often assumed to play is to correct self-interaction error-driven electron delocalization (e.g., from transition metal centers to neighboring ligands). Surprisingly, we instead observe that increasing Hartree-Fock exchange reduces charge on iron centers, corresponding to effective delocalization of charge to ligands, thus challenging notions of the role of Hartree-Fock exchange in shifting predictions of spin-state ordering. C 2015 AIP Publishing LLC. [http://dx

Virtual high throughput screening, typically driven by first-principles, density functional theory calculations, has emerged as a powerful tool for the discovery of new materials. Although the computational materials science community has benefited from open source tools for the rapid structure generation, calculation, and analysis of crystalline inorganic materials, software and strategies to address the unique challenges of inorganic complex discovery have not been as widely available. We present a unified view of our recent developments in the open source molSimplify code for inorganic discovery. Building on our previous efforts in the automated generation of highly accurate inorganic molecular structures, first-principles simulation, and property analysis to accelerate high-throughput screening, we have recently incorporated a neural network that both improves structure generation and predicts electronic properties prior to first-principles calculation. We also provide an overview of how multi-million molecule organic libraries can be leveraged for inorganic discovery alongside cheminformatics concepts of molecular diversity in order to efficiently traverse chemical space. We demonstrate all of these tools on the discovery of design rules for octahedral Fe(II/III) redox couples with nitrogen ligands. Over a search of only approximately 40 new molecules, we obtain redox potentials relative to the Fc/Fc + couple ranging from-1 to 4.5 V in aqueous solution. Our new automated correlation analysis reveals heteroatom identity and the degree of structural branching to be key ligand descriptors in determining redox potential. This inorganic discovery toolkit provides a promising approach to advancing transition metal complex design.

Piecewise linearity of the energy with respect to fractional electron removal or addition is a requirement of an electronic structure method that necessitates the presence of a derivative discontinuity at integer electron occupation. Semi-local exchange-correlation (xc) approximations within density functional theory (DFT) fail to reproduce this behavior, giving rise to deviations from linearity with a convex global curvature that is evidence of many-electron, self-interaction error and electron delocalization. Popular functional tuning strategies focus on reproducing piecewise linearity, especially to improve predictions of optical properties. In a divergent approach, Hubbard U-augmented DFT (i.e., DFT+U) treats self-interaction errors by reducing the local curvature of the energy with respect to electron removal or addition from one localized subshell to the surrounding system. Although it has been suggested that DFT+U should simultaneously alleviate global and local curvature in the atomic limit, no detailed study on real systems has been carried out to probe the validity of this statement. In this work, we show when DFT+U should minimize deviations from linearity and demonstrate that a "+U" correction will never worsen the deviation from linearity of the underlying xc approximation. However, we explain varying degrees of efficiency of the approach over 27 octahedral transition metal complexes with respect to transition metal (Sc-Cu) and ligand strength (CO, NH3, and H2O) and investigate select pathological cases where the delocalization error is invisible to DFT+U within an atomic projection framework. Finally, we demonstrate that the global and local curvatures represent different quantities that show opposing behavior with increasing ligand field strength, and we identify where these two may still coincide.

Prediction of spin-state ordering in transition metal complexes is essential for understanding catalytic activity and designing functional materials. Semi-local approximations in density functional theory, such as the generalized-gradient approximation (GGA), suffer from several errors notably including delocalization error that give rise to systematic bias for more covalently bound low-spin electronic states. Incorporation of exact exchange is known to counteract this bias, instead favoring high-spin states, in a manner that has recently been identified to be ligand-field dependent. In this work, we introduce a tuning strategy to identify the effect of incorporating the Laplacian of the density (i.e., a meta-GGA) in exchange on spinstate ordering. We employ a diverse test set of M(II) and M(III) first-row transition metal ions from Ti to Cu as well as octahedral complexes of these ions with ligands of increasing field strength (i.e., H 2 O, NH 3 , and CO). We show that the sensitivity of spin-state ordering to meta-GGA exchange is highly ligand-field dependent, stabilizing high-spin states in strong-field (i.e., CO) cases and stabilizing low-spin states in weak-field (i.e., H 2 O, NH 3 , and isolated ions) cases. This diverging behavior leads to generally improved treatment of isolated ions and strong field complexes over a standard GGA but worsened treatment for the hexa-aqua or hexa-ammine complexes. These observations highlight the sensitivity of functional performance to subtle changes in chemical bonding.2