Chemically Separable Co(II) Spin-State Isomers

The phenomenon of spin crossover involves coordination complexes with switchable spin states. This spin state change is accompanied by significant geometric changes such that low and high spin forms of a complex are distinct isomers that exist in equilibrium with one another. Typically, spin-state isomers interconvert rapidly and are similar enough in polarity to prevent their independent separation and isolation. We report here the first example, to our knowledge, of cobalt(II) spin-state isomers that can be physically separated. The reaction of Mo2(dpa)4 (dpa = 2,2'-dipyridylamide) with CoBr2 produces a mixture of two heterometallic compounds with a linear, metal-metal-bonded Mo[Formula: see text]Mo-Co chain. The complexes, SC-[BrMo2(dpa)4Co]Br (SC-2) and HS-[BrMo2(dpa)4CoBr] (HS-2), have identical compositions (Mo2Co(dpa)4Br2) but different ground spin states and coordination geometries of the Co(II) ion. In the solid state, SC-2 undergoes incomplete spin crossover from an S = 1/2 state to an S = 3/2 state, and HS-2 has a high spin, S = 3/2, ground state, as confirmed by SQUID magnetometry and EPR spectroscopy. Crystallographic analyses of SC-2 and HS-2 show that SC-2 has an elongated Co-Br distance relative to HS-2 and is best described as the salt [BrMo2(dpa)4Co]Br. This limits SC-2's solubility in nonpolar solvents and allows for the physical separation of the two isomers. Solution studies of SC-2 and HS-2 indicate that SC-2 and HS-2 interconvert slowly relative to the NMR time scale. Additional solution-state EPR and UV-vis absorption measurements demonstrate that the choice of solvent polarity determines the predominant isomer present in solution.

Electrocatalytic Reduction of CO2 and H2O with Zn(II) Complexes Through Metal-Ligand Cooperation

Air and water-stable zinc (II) complexes of neutral pincer bis(diphenylphosphino)-2,6-di(amino)pyridine ("PN3P") ligands are reported. These compounds, [Zn(κ2-2,6-{Ph2PNR}2(NC5H3))Br2] (R=Me, 1; R=H, 2), were shown to be capable of electrocatalytic reduction of CO2 at -2.3 V vs. Fc+/0 to selectively yield CO in mixed water/acetonitrile solutions. These complexes also electrocatalytically generate H2 from water in acetonitrile solutions, at the same potential, with Faradaic efficiencies of up to 90 %. DFT computations support a proposed mechanism involving the first reduction of 1 or 2 occurring at the PN3P ligand. Furthermore, computational analysis suggested a mechanism involving metal-ligand cooperation of a Lewis acidic Zn(II) and a basic ligand.

Insights into the deviation from piecewise linearity in transition metal complexes from supervised machine learning models

Virtual high-throughput screening (VHTS) and machine learning (ML) with density functional theory (DFT) suffer from inaccuracies from the underlying density functional approximation (DFA). Many of these inaccuracies can be traced to the lack of derivative discontinuity that leads to a curvature in the energy with electron addition or removal. Over a dataset of nearly one thousand transition metal complexes typical of VHTS applications, we computed and analyzed the average curvature (i.e., deviation from piecewise linearity) for 23 density functional approximations spanning multiple rungs of "Jacob's ladder". While we observe the expected dependence of the curvatures on Hartree-Fock exchange, we note limited correlation of curvature values between different rungs of "Jacob's ladder". We train ML models (i.e., artificial neural networks or ANNs) to predict the curvature and the associated frontier orbital energies for each of these 23 functionals and then interpret differences in curvature among the different DFAs through analysis of the ML models. Notably, we observe spin to play a much more important role in determining the curvature of range-separated and double hybrids in comparison to semi-local functionals, explaining why curvature values are weakly correlated between these and other families of functionals. Over a space of 187.2k hypothetical compounds, we use our ANNs to pinpoint DFAs for which representative transition metal complexes have near-zero curvature with low uncertainty, demonstrating an approach to accelerate screening of complexes with targeted optical gaps.

Elucidating Two Distinct Pathways for Electrocatalytic Hydrogen Production Using CoII Pincer Complexes

Hydrogen gas is a sustainable energy source with water as the sole combustion product. As a result, efforts to catalyze H₂ production are pertinent and widespread. The electrocatalytic H₂ generating capabilities of two Co^II complexes, [Co(κ³ -2,6-{Ph₂ PNR}₂ (NC₅ H₃ ))Br₂ ] with R=H (I) or R=Me (II), were presented for a variety of proton sources including trifluoroacetic acid (TFA), acetic acid (AA), and trifluoroethanol (TFE). Cyclic voltammetry and controlled potential coulometry demonstrated that electrocatalysis from I and II occurred at two different potentials and are associated with different reduction processes. Density functional theory analysis provided insight into the identities of the catalyst and supported two distinct reaction pathways for electrocatalytic proton reduction. Specifically, stronger acids (e. g., AA, TFA) proceeded at -1.31 to -1.45 V through a M^I /M^III pathway while sources with higher pK_a values (e. g., TFE, H₂ O) generated hydrogen at -2.4 V via M⁰ /M^II ligand-assisted metal-centered reduction.

Slow-Relaxation Behavior of a Mononuclear Co(II) Complex Featuring Long Axial Co-O Bond

Co(II) mononuclear complex with different coordination geometry would display various of field-induced single-ion magnet (SIM) behaviors. Here, we identify a field-induced single-ion magnet in a mononuclear complex Co(H2DPA)2·H2O (H2DPA = 2,6-pyridine-dicarboxylic acid) by the hydrothermal method. The long axial Co-O coordination bond (Co1‧‧‧O3) can be formed by Co1 and O3. Therefore, Co(II) ion is six-coordinated in a distorted elongated octahedron. AC magnetization susceptibilities show that the effective energy barrier is up to 43.28 K. This is much larger than most mononuclear Co(II). The distorted elongated octahedron caused by the axial Co-O coordination bond is responsible for the high effective energy barrier. The distribution of electron density in Co1 and O3 atoms in the long axial bond would influence the magnetic relaxation process in turn. Our work deepens the relationship between the effective energy barrier and the weak change of ligand field by long axial bonds, which would facilitate constructing SIM with high energy temperature.

Catalytic and mechanistic studies of a highly active and E-selective Co(II) PNNH pincer catalyst system for transfer-semihydrogenation of internal alkynes

Herein we report the application of a Co(II) PNNH pincer catalyst system (PNNH = 2-(5-(t-butyl)-1H-pyrazol-3-yl)-6-(dialkylphosphinomethyl)pyridine) for the highly E-selective transfer semihydrogenation of internal diaryl alkynes using methanol and ammonia borane...

Deciphering Cryptic Behavior in Bimetallic Transition-Metal Complexes with Machine Learning

We demonstrate an alternative, data-driven approach to uncovering structure-property relationships for the rational design of heterobimetallic transition-metal complexes that exhibit metal-metal bonding. We tailor graph-based representations of the metal-local environment for these complexes for use in multiple linear regression and kernel ridge regression (KRR) models. We curate a set of 28 experimentally characterized complexes to develop a multiple linear regression model for oxidation potentials. We achieve good accuracy (mean absolute error of 0.25 V) and preserve transferability to unseen experimental data with a new ligand structure. We also train 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 these 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.

Machine learning to tame divergent density functional approximations: a new path to consensus materials design principles

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 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, “rungs” (e.g., semi-local to double hybrid) and basis sets on over 2000 TMCs. Although computed property values (e.g., spin state splitting and frontier orbital gap) 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, providing DFA-invariant, universal design rules. We devise a strategy to train artificial neural network (ANN) models informed by all 23 DFAs and use them to predict properties (e.g., spin-splitting energy) of over 187k TMCs. By requiring consensus of the ANN-predicted DFA properties, we improve correspondence of computational lead compounds with literature-mined, experimental compounds over the typically employed single-DFA approach.

Machine learning (ML)-based feature analysis reveals universal design rules regardless of density functional choices. Using the consensus among multiple functionals, we identify robust lead complexes in ML-accelerated chemical discovery.

Electrocatalytic H2 Generation from Water Relying on Cooperative Ligand Electron Transfer in "PN3 P" Pincer-Supported NiII Complexes

Water is the most sustainable source for H₂ production, and the efficient electrocatalytic production of H₂ from mixed water/acetonitrile solutions by using two new air-stable nickel(II) pincer complexes, [Ni(κ³ -2,6-{Ph₂ PNR}₂ (NC₅ H₃ )Br₂ ] (R=H I, Me II) is reported. Hydrogen generation from H₂ O/CH₃ CN solutions is initiated at -2 V against Fc^+/0 , and bulk electrocatalysis studies showed that the catalyst functions with an excellent Faradaic efficiency and a turnover frequency of 160 s^-1 . A DFT computational investigation of the reduction behavior of I and II revealed a correlation of H₂ formation with charge donation from electrons originating in a reduced ligand-localized orbital. As a result, these catalysts are proposed to proceed by a novel mechanism involving electron/proton transfer between a Ni^0I species bonded to an anionic PN³ P ligand ("L^- /Ni^0I ") and a Ni^I -hydride ("Ni-H"). Furthermore, these catalysts are able to reduce phenol and acetic acid, more active proton sources, at lower potentials that correlate with the substrate pK_a .

Rapid Detection of Strong Correlation with Machine Learning for Transition-Metal Complex High-Throughput Screening

Despite its widespread use in chemical discovery, approximate density functional theory (DFT) is poorly suited to many targets, such as those containing open-shell, 3d transition metals that can be expected to have strong multireference (MR) character. For discovery workflows to be predictive, we need automated, low-cost methods that can distinguish the regions of chemical space where DFT should be applied from those where it should not. We curate more than 4800 open-shell transition-metal complexes up to hundreds of atoms in size from prior high-throughput DFT studies and evaluate affordable, finite-temperature DFT fractional occupation number (FON)-based MR diagnostics. We show that intuitive measures of strong correlation (i.e., the HOMO-LUMO gap) are not predictive of MR character as judged by FON-based diagnostics. Analysis of independently trained machine learning (ML) models to predict HOMO-LUMO gaps and FON-based diagnostics reveals differences in the metal and ligand sensitivity of the two quantities. We use our trained ML models to rapidly evaluate MR character over a space of ∼187000 theoretical complexes, identifying large-scale trends in spin-state-dependent MR character and finding small HOMO-LUMO gap complexes while ensuring low MR character.

Single molecule magnet behaviour in a square planar S = 1/2 Co(ii) complex and spin-state assignment of multiple relaxation modes

We report the first example of field-induced single molecule magnet (SMM) behaviour in a square-planar S = 1/2 Co(ii) pincer complex [(PNNNP)CoBr]Br (2). The related five-coordinate complexes [(PCNCP)CoBr2] (1) and [(PONOP)CoBr2] (3) also exhibit SMM properties. Partial spin crossover displayed by 3 allows for assignment of distinct relaxation modes to each spin state.

Optimal diamagnetic dilution concentration for suppressing the dipole-dipole interaction in single-ion magnets

The effect of screening the Co^II moment of monomeric [Co^IIL₂(H₂O)] (L = 8-hydroxyquinaldine), having a trigonal bipyramid coordination, by diamagnetic Zn in Co_xZn_1-x solid solutions on its magnetic relaxation was explored using ac-susceptibility, high-field electron-spin-resonance measurements and CASPT2 calculations. The retention of the crystal structure for all the solid solutions was demonstrated using single crystal diffraction. The dc-magnetization and theoretical fittings of the susceptibility for Co₁ and Co_0.1Zn_0.9 gave a large zero-field-splitting (ZFS) D of 50 ± 6 cm^-1, and very weak dipole interaction between the nearest neighbors, while EPR and calculations confirmed the positive sign of the axial component (D). Consistent parameters were obtained from experiments and theory. Importantly, only field-induced relaxation was observed for the samples with less than 50% Co and a gradual change in the barrier energy to moment reversal and relaxation times was observed between 11% and 20% Co, while both were enhanced for higher dilutions. The results establish a clear barrier for extending the longevity of the magnetism for this type of single-ion species by lowering the intramolecular interactions. The results suggest that the magnetic interaction persists up to the second sphere, that is, for a dilution of 1 in 9 (11% Co). Importantly, this method is applicable to all single-ion magnet systems, that is, the optimum dilution concentration to restrain the dipole field can be given only by the single crystal structure.

Pogo-Stick Iron and Cobalt Complexes: Synthesis, Structures, and Magnetic Properties

The synthesis, structures, and magnetic properties of monomeric half-sandwich iron and cobalt imidazolin-2-iminato complexes have been comprehensively investigated. Salt metathesis reactions of [Cp'M(μ-I)]₂ (1-M, M = Fe, Co; Cp' = η⁵-1,2,4-tri-tert-butylcyclopentadienyl) with [Im^DippNLi]₂ (Im^DippN = 1,3-bis(2,6-diisopropylphenyl)imidazolin-2-iminato) furnishes the terminal half-sandwich compounds [Cp'M(NIm^Dipp)] (2-M, M = Fe, Co), which can be regarded as models for elusive half-sandwich iron and cobalt imido complexes. X-ray diffraction analysis confirmed the structure motif of a one-legged piano stool. Complex 2-Co can also be prepared by an acid-base reaction between [Cp'Co{N(SiMe₃)₂}] (3-Co) and Im^DippNH. The electronic and magnetic properties of 2-M and 3-Co were probed by ⁵⁷Fe Mössbauer spectroscopy (M = Fe), X-band EPR spectroscopy (M = Co), and solid-state magnetic susceptibility measurements. In particular, the central metal atom adopts a high-spin (S = 2) state in 2-Fe, while the cobalt complex 2-Co represents a rare example of a Co(II) species with a coordination number different from six displaying a low-spin to high-spin spin-crossover (SCO) behavior. The experimental observations are complemented by DFT calculations.

Designing in the Face of Uncertainty: Exploiting Electronic Structure and Machine Learning Models for Discovery in Inorganic Chemistry

Recent transformative advances in computing power and algorithms have made computational chemistry central to the discovery and design of new molecules and materials. First-principles simulations are increasingly accurate and applicable to large systems with the speed needed for high-throughput computational screening. Despite these strides, the combinatorial challenges associated with the vastness of chemical space mean that more than just fast and accurate computational tools are needed for accelerated chemical discovery. In transition-metal chemistry and catalysis, unique challenges arise. The variable spin, oxidation state, and coordination environments favored by elements with well-localized d or f electrons provide great opportunity for tailoring properties in catalytic or functional (e.g., magnetic) materials but also add layers of uncertainty to any design strategy. We outline five key mandates for realizing computationally driven accelerated discovery in inorganic chemistry: (i) fully automated simulation of new compounds, (ii) knowledge of prediction sensitivity or accuracy, (iii) faster-than-fast property prediction methods, (iv) maps for rapid chemical space traversal, and (v) a means to reveal design rules on the kilocompound scale. Through case studies in open-shell transition-metal chemistry, we describe how advances in methodology and software in each of these areas bring about new chemical insights. We conclude with our outlook on the next steps in this process toward realizing fully autonomous discovery in inorganic chemistry using computational chemistry.

A mononuclear five-coordinate Co(ii) single molecule magnet with a spin crossover between the S = 1/2 and 3/2 states

Although a great number of single-ion magnets (SIMs) and spin-crossover (SCO) compounds have been discovered, multifunctional materials with the combination of SCO and SIM properties are extremely scarce. Here magnetic studies have been carried out for a mononuclear, five-coordinate cobalt(ii) complex [Co(3,4-lut)4Br]Br (1) with square pyramidal geometry. Direct-current magnetic measurement confirms the spin transition between the S = 1/2 and 3/2 states in the range of 150-290 K with a small hysteresis loop. Frequency- and temperature-dependent alternating-current magnetic susceptibility reveals slow magnetization relaxation under an applied dc field of 3000 Oe. The work here presents the first instance of the five-coordinate mononuclear cobalt(ii)-based SIM exhibiting the thermally induced complete SCO.

Slow magnetic relaxation in five-coordinate spin-crossover cobalt(ii) complexes

We present the first examples of the coexistence of field induced slow magnetic relaxation and spin-crossover observed in five-coordinate cobalt(ii) complexes [Co(12-TMC)(CH₃CN)](X)₂ (12-TMC = 1,4,7,10-tetramethyl-1,4,7,10-tetraazacyclododecane, X = BF₄^-, 1; PF₆^-, 2). The direct-current (dc) magnetic data show a gradual and incomplete spin-crossover at high temperature. 1 and 2 display frequency- and temperature-dependent alternating-current (ac) magnetic susceptibility under an applied dc field of 2500 Oe, which originates from the S = 1/2 spin state of Co(ii) ions.

Ruthenium complexes featuring cooperative phosphine–pyridine–iminophosphorane (PNN) ligands: synthesis, reactivity and catalytic activity

Ruthenium complexes with an iminophosphorane based (PNN) ligand; the NP substituent influences the coordination and the reactivity of the formed complexes.

Electro- and Photocatalytic Generation of H2 Using a Distinctive CoII "PN3 P" Pincer Supported Complex with Water or Saturated Saline as a Hydrogen Source

Efficient electrocatalytic production of H₂ from mixed water/acetonitrile solutions was achieved using three new Co^II complexes supported by the neutral pincer ligand bis(diphenylphosphino)-2,6-di(methylamino)pyridine ("PN₃ P"). At -1.9 V vs. Fc/Fc⁺ , these catalysts showed 96 % Faradaic efficiency with added water or saturated aqueous saline at rates of up to 316 L(mol cat)^-1 (cm² )^-1  h^-1 using a glassy carbon working electrode. The complex [Co(κ³ -2,6-{Ph₂ PNMe}₂ (NC₅ H₃ )Br₂ ] (1) was also able to photocatalytically reduce water to hydrogen in the presence of a Ru(bpy)₃²⁺ photosensitizer and a reductant.

Crystal structure and spectroscopic characterization of a cobalt(II) tetraazamacrocycle: completing a series of first-row transition-metal complexes

The tetraazamacrocyclic ligand 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane (TMC) has been used to bind a variety of first-row transition metals but to date the crystal structure of the cobalt(II) complex has been missing from this series. The missing cobalt complex chlorido(1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane-κ⁴N)cobalt(II) chloride dihydrate, [CoCl(C₁₄H₃₂N₄)]Cl·2H₂O or [Co^IICl(TMC)]Cl·2H₂O, crystallizes as a purple crystal. This species adopts a distorted square-pyramidal geometry in which the TMC ligand assumes the trans-I configuration and the chloride ion binds in the syn-methyl pocket of the ligand. The Co^II ion adopts an S = 3/2 spin state, as measured by the Evans NMR method, and UV-visible spectroscopic studies indicate that the title hydrated salt is stable in solution. Density functional theory (DFT) studies reveal that the geometric parameters of [Co^IICl(TMC)]Cl·2H₂O are sensitive to the cobalt spin state and correctly predict a change in spin state upon a minor perturbation to the ligand environment.