Search

To accelerate exploration of chemical space, it is necessary to identify the compounds that will provide the most additional information or value. A large-scale analysis of mononuclear octahedral transition metal complexes deposited in an experimental database confirms an under-representation of lower-symmetry complexes. From a set of around 1000 previously studied Fe(II) complexes, we show that the theoretical space of synthetically accessible complexes formed from the relatively small number of unique ligands is significantly (ca. 816k) larger. For the properties of these complexes, we validate the concept of ligand additivity by inferring heteroleptic properties from a stoichiometric combination of homoleptic complexes. An improved interpolation scheme that incorporates information about cis and trans isomer effects predicts the adiabatic spin-splitting energy to around 2 kcal/mol and the HOMO level to less than 0.2 eV. We demonstrate a multi-stage strategy to discover leads from the 816k Fe(II) complexes within a targeted property region. We carry out a coarse interpolation from homoleptic complexes that we refine over a subspace of ligands based on the likelihood of generating complexes with targeted properties. We validate our approach on 9 new binary and ternary complexes predicted to be in a targeted zone of discovery, suggesting opportunities for efficient transition metal complex discovery.

Accurate virtual high-throughput screening (VHTS) of transition metal complexes (TMCs) remains challenging due to the possibility of high multi-reference (MR) character that complicates property evaluation. We compute MR diagnostics for over 5,000 ligands present in previously synthesized transition metal complexes in the Cambridge Structural Database (CSD). To accomplish this task, we introduce an iterative approach for consistent ligand charge assignment for ligands in the CSD. Across this set, we observe that MR character correlates linearly with the inverse value of the averaged bond order over all bonds in the molecule. We then demonstrate that ligand additivity of MR character holds in TMCs, which suggests that the TMC MR character can be inferred from the sum of the MR character of the ligands. Encouraged by this observation, we leverage ligand additivity and develop a ligand-derived machine learning representation to train neural networks to predict the MR character of TMCs from properties of the constituent ligands. This approach yields models with excellent performance and superior transferability to unseen ligand chemistry and compositions.

Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition metal chemistry. Here, we demonstrate the judiciously-modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition metal complexes. We curate a set of nine representative Ti(III) and V(IV) $d^1$ transition metal complexes and evaluate their flat plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals for both the case of fractional spins and fractional charges. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of $d^1$ MOPs to be suitable for nearly eliminating all energetic delocalization error and static correlation error. In all cases, the MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat plane errors at non-empirical values. Unlike widely employed DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization error and static correlation error within a non-empirical framework.

The rational tailoring of transition metal complexes is necessary to address outstanding challenges in energy utilization and storage. Heterobimetallic transition metal complexes that exhibit metal-metal bonding in stacked "double decker" ligand structures are an emerging, attractive platform for catalysis, but their properties are challenging to predict prior to laborious synthetic efforts. We demonstrate an alternative, data-driven approach to uncovering structure-property relationships for rational bimetallic complex design. We tailor graph-based representations of the metal-local environment for these heterobimetallic complexes for use in training of multiple linear regression and kernel ridge regression (KRR) models. Focusing on oxidation potentials, we obtain a set of 28 experimentally characterized complexes to develop a multiple linear regression model. On this training set, we achieve good accuracy (mean absolute error, MAE, of 0.25 V) and preserve transferability to unseen experimental data with a new ligand structure. We trained a KRR model on a subset of 330 structurally characterized heterobimetallics to predict the degree of metal-metal bonding. This KRR model predicts relative metal-metal bond lengths in the test set to within 5%, and analysis of key features reveals the fundamental atomic contributions (e.g., the valence electron configuration) that most strongly influence the behavior of complexes. Our work provides guidance for rational bimetallic design, suggesting that properties including the formal shortness ratio should be transferable from one period to another.

Modifying the optoelectronic properties of nanostructured materials through introduction of dopant atoms has attracted intense interest. Nevertheless, the approaches employed are often trial and error, preventing rational design. We demonstrate the power of large-scale electronic structure calculations with density functional theory (DFT) to build an atlas of preferential dopant sites for a range of M(II) and M(III) dopants in the representative III-V InP magic sized cluster (MSC). We quantify the thermodynamic favorability of dopants, which we identify to be both specific to the sites within the MSC (i.e., interior vs surface) and to the nature of the dopant atom (i.e., smaller Ga(III) vs larger Y(III) or Sc(III)). These observations motivate development of maps of the most and least favorable doping sites, which are consistent with some known experimental expectations but also yield unexpected observations. For isovalent doping (i.e., Y(III)/Sc(III) or Ga(III), we observed stronger sensitivity of the predicted energetics to the type of ligand orientation on the surface than to the dopant type, but divergent behavior is observed for whether interior doping is favorable. For charge balancing with M(II) (i.e., Zn or Cd) dopants, we show that the type of ligand removed during the doping reaction is critical. We show that limited cooperativity with dopants up to moderate concentrations occurs, indicating rapid single-dopant estimations of favorability from DFT can efficiently guide rational design. Our work emphasizes the strong importance of ligand chemistry and surface heterogeneity in determining paths to favorable doping in quantum dots, an observation that will be general to other III-V and II-VI quantum dot systems generally synthesized with carboxylate ligands.

Computational virtual high-throughput screening (VHTS) with density functional theory (DFT) and machine-learning (ML)-acceleration is essential in rapid materials discovery. By necessity, efficient DFT-based workflows are carried out with a single density functional approximation (DFA). Nevertheless, properties evaluated with different DFAs can be expected to disagree for the cases with challenging electronic structure (e.g., open shell transition metal complexes, TMCs) for which rapid screening is most needed and accurate benchmarks are often unavailable. To quantify the effect of DFA bias, we introduce an approach to rapidly obtain property predictions from 23 representative DFAs spanning multiple families and "rungs" (e.g., semi-local to double hybrid) and basis sets on over 2,000 TMCs. Although computed properties (e.g., spin-state ordering and frontier orbital gap) naturally differ by DFA, high linear correlations persist across all DFAs. We train independent ML models for each DFA and observe convergent trends in feature importance; these features thus provide DFA-invariant, universal design rules. We devise a strategy to train ML models informed by all 23 DFAs and use them to predict properties (e.g., spin-splitting energy) of over 182k TMCs. By requiring consensus of the ANN-predicted DFA properties, we improve correspondence of these computational lead compounds with literature-mined, experimental compounds over the single-DFA approach typically employed. Both feature analysis and consensus-based ML provide efficient, alternative paths to overcome accuracy limitations of practical DFT.