Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes

Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.

Spectroscopic characterization of carbon monoxide activation by neutral chromium carbides

Spectroscopic characterization of carbon monoxide activation by neutral metal carbides is of essential importance for understanding the structure-reactivity relationships of catalytic sites, but has been proven to be very challenging owing to the difficulty in size selection. Here, we report a size-specific infrared-vacuum ultraviolet spectroscopic study of the reactions between carbon monoxide with neutral chromium carbides. Quantum chemical calculations were carried out to identify the low-lying structures and to interpret the experimental features. The results reveal that the most stable structure of CrC3(CO)2 consists of a CCO ketenylidene unit and that of CrC4(CO)2 has a semi-bridging CO with a very low CO stretching vibrational frequency at 1821 cm-1. The electron structure analyses show that this semi-bridging CO is highly activated through the delocalized Cr-C-C three-center two-electron (3c-2e) interaction between the antibonding orbitals of CO and the metal carbide skeleton. The formation of these metal carbide carbonyls is found to be both thermodynamically exothermic and kinetically facile in the gas phase. The present findings have important implications for the mechanical understanding of the catalytic processes with isolated metal atoms/clusters dispersed on supports.

Local Oxidation States in {FeNO}6-8 Porphyrins: Insights from DMRG/CASSCF-CASPT2 Calculations

A first DMRG/CASSCF-CASPT2 study of a series of paradigmatic {FeNO}6, {FeNO}7, and {FeNO}8 heme-nitrosyl complexes has led to substantial new insight as well as uncovered key shortcomings of the DFT approach. By virtue of its balanced treatment of static and dynamic correlation, the calculations have provided some of the most authoritative information available to date on the energetics of low- versus high-spin states of different classes of heme-nitrosyl complexes. Thus, the calculations indicate low doublet-quartet gaps of 1-4 kcal/mol for {FeNO}7 complexes and high singlet-triplet gaps of ≳20 kcal/mol for both {FeNO}6 and {FeNO}8 complexes. In contrast, DFT calculations yield widely divergent spin state gaps as a function of the exchange-correlation functional. DMRG-CASSCF calculations also help calibrate DFT spin densities for {FeNO}7 complexes, pointing to those obtained from classic pure functionals as the most accurate. The general picture appears to be that nearly all the spin density of Fe[P](NO) is localized on the Fe, while the axial ligand imidazole (ImH) in Fe[P](NO)(ImH) pushes a part of the spin density onto the NO moiety. An analysis of the DMRG-CASSCF wave function in terms of localized orbitals and of the resulting configuration state functions in terms of resonance forms with varying NO(π*) occupancies has allowed us to address the longstanding question of local oxidation states in heme-nitrosyl complexes. The analysis indicates NO(neutral) resonance forms [i.e., Fe(II)-NO0 and Fe(III)-NO0] as the major contributors to both {FeNO}6 and {FeNO}7 complexes. This finding is at variance with the common formulation of {FeNO}6 hemes as Fe(II)-NO+ species but is consonant with an Fe L-edge XAS analysis by Solomon and co-workers. For the {FeNO}8 complex {Fe[P](NO)}-, our analysis suggests a resonance hybrid description: Fe(I)-NO0 ↔ Fe(II)-NO-, in agreement with earlier DFT studies. Vibrational analyses of the compounds studied indicate an imperfect but fair correlation between the NO stretching frequency and NO(π*) occupancy, highlighting the usefulness of vibrational data as a preliminary indicator of the NO oxidation state.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Perspective from a Hubbard U-density corrected scheme towards a spin crossover-mediated change in gas affinity

By employing a recently proposed Hubbard U density-corrected scheme within density functional theory, we provide design principles towards the design of materials exhibiting a spin crossover-assisted gas release. Small molecular fragments are used as case study to identify two main mechanisms behind the change in binding energy upon spin transitions. The feasibility of the proposed mechanism in porous crystals is assessed by correlating the change in binding energy of CO2, CO, N2, and H2, upon spin crossover, with the adiabatic energy difference associated with the spin state change of the square-planar metal in Hofmann-type clathrates (M = Fe, Mn, Ni). A few promising cases are identified for the adsorption of intermediate ligand field strength molecules such as N2 and H2. The latter stands out as the most original result as the strong interaction in low spin, as expected from a Kubas mechanism, results in a large change in binding energy. This work provides a general perspective towards the engineering of open-metal site frameworks exhibiting local environments designed to have a spin crossover upon adsorption of specific gas molecules.

Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes

Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu²⁺-containing complexes. For metal-nucleic acid systems that do not contain Cu²⁺, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.

Transition Metal-Decorated B12N12-X (X = Au, Cu, Ni, Os, Pt, and Zn) Nanoclusters as Biosensors for Carboplatin

Theoretical studies on the adsorption, sensibility, and reactivity of a boron nitride nanocage decorated with Au, Cu, Ni, Os, Pt, and Zn metals as a biosensor material were carried out for the adsorption of carboplatin by applying the density functional theory computation at the B3LYP-GD3BJ/def2svp level of theory. All the optimized structures, as well as the calculations as regards the studied objective including electronic properties, geometry optimization parameters, adsorption energy studies, natural bond orbital analysis, topology studies, sensor mechanistic parameters, and thermodynamic properties (ΔG and ΔH), were investigated herein. As a result, the noticeable change in the energy gap of the studied surfaces when interacting with carboplatin accounted for the surfaces' reactivity, stability, conductivity, work function, and overall adsorption ability, implying that the studied decorated surfaces are good sensor materials for sensing carboplatin. Furthermore, the negative adsorption energies obtained for interacting surfaces decorated with Cu, Ni, Os, and Zn suggest that the surface has a superior ability to sense carboplatin as chemisorption was seen. Substantially, the geometric short adsorption bond length after adsorption, thermodynamically spontaneous reactions, and acceptable sensor mechanism results demonstrate that the investigated surfaces have strong sensing characteristics for sensing carboplatin.

Photodissociation Spectroscopy and Photofragment Imaging of the Fe+(Acetylene) Complex

Tunable laser photodissociation spectroscopy in the 700-400 nm region and photofragment imaging experiments are employed to investigate the Fe⁺(acetylene) ion-molecule complex. At energies above a threshold at 679 nm, continuous dissociation is detected throughout the visible wavelength region, with regions of broad structure. Comparison to the spectrum predicted by time-dependent density functional theory (TD-DFT) indicates that the complex has a quartet ground state. The dissociation threshold for Fe⁺(acetylene) at 679 nm provides the dissociation energy on the quartet potential energy surface. Correction for the atomic quartet-sextet spin state energy difference provides an adiabatic dissociation energy of 36.8 ± 0.2 kcal/mol. Photofragment imaging of the Fe⁺ photoproduct produced at 603.5 nm produces significant kinetic energy release (KER). The photon energy and the maximum value of the KER provide an upper limit on the dissociation energy of D₀ ≤ 34.6 ± 3.2 kcal/mol. The dissociation energies determined from the spectroscopy and photofragment imaging experiments agree nicely with the value determined previously by collision-induced dissociation (38.0 ± 2.6 kcal/mol). However, both values are significantly lower than those produced by computational chemistry at the DFT level using different functionals recommended for transition-metal chemistry.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Toward a Reliable Theory Benchmark

A new composite method for the calculation of spin-crossover energies in 3d transition-metal complexes based on multireference methods is presented. The method reduces to MRCISD+Q at the complete-basis-set (CBS) level for atomic ions, for which it gives excitation energies with a mean absolute error of only ca. 0.01 eV. For molecular complexes, the CASPT2+δMRCI composite approach corresponds to a CASPT2/CBS calculation augmented by a high-level MRCISD+Q-CASPT2 correction with a smaller ligand basis set. For a set of [Fe(He)6]n+ test complexes, the approach reproduces full MRCISD+Q/CBS results to within better than 0.04 eV, without depending on any arbitrary IPEA shifts. The high-quality CASPT2+δMRCI method has then been applied to a series of 3d transition-metal hexaqua complexes in aqueous solution, augmented by an elaborate 3D-RISM-SCF solvent treatment of the underlying structures. It provides unprecedented agreement with experiment for the lowest-lying vertical spin-flip excitation energies, except for the Fe³⁺ system. Closer examination of the latter case provides strong evidence that the observed lowest-energy excitation at 1.56 eV, which has been used frequently for evaluating quantum-chemical methods, does not arise from the iron(III) hexaqua complex in solution, but from its singly deprotonated counterpart, [Fe(H2O)5OH]2+.

A transferable recommender approach for selecting the best density functional approximations in chemical discovery

Approximate density functional theory has become indispensable owing to its balanced cost-accuracy trade-off, including in large-scale screening. To date, however, no density functional approximation (DFA) with universal accuracy has been identified, leading to uncertainty in the quality of data generated from density functional theory. With electron density fitting and Δ-learning, we build a DFA recommender that selects the DFA with the lowest expected error with respect to the gold standard (but cost-prohibitive) coupled cluster theory in a system-specific manner. We demonstrate this recommender approach on the evaluation of vertical spin splitting energies of transition metal complexes. Our recommender predicts top-performing DFAs and yields excellent accuracy (about 2 kcal mol-1) for chemical discovery, outperforming both individual Δ-learning models and the best conventional single-functional approach from a set of 48 DFAs. By demonstrating transferability to diverse synthesized compounds, our recommender potentially addresses the accuracy versus scope dilemma broadly encountered in computational chemistry.

Coordination and Spin States in Fe+(C2H2)n Complexes Studied with Selected-Ion Infrared Spectroscopy

Fe⁺(acetylene)_n ion-molecule complexes are produced in a supersonic molecular beam with pulsed laser vaporization. These ions are mass selected and studied with infrared photodissociation spectroscopy in the C-H stretching region, complemented by computational chemistry calculations. All C-H stretch vibrations are shifted to frequencies lower than the vibrations of isolated acetylene because of the charge transfer that occurs between the metal ion and the molecules. Complexes in the size range of n = 1-4 are found to have structures with individual acetylene molecules bound to the core metal ion via cation-π interactions. The coordination is completed with four ligands in a structure close to a distorted tetrahedron. Larger complexes in the range of n = 5-8 have external acetylene molecules solvating this n = 4 core ion via CH-π bonding to inner-shell ligands. DFT computations predict that quartet spin states are more stable for all complex sizes, but infrared spectra for quartet and doublet spin states are quite similar, precluding definitive determination of the spin states. There is no evidence for any of these complexes having acetylenes coupled into reacted structures. This is consistent with computed thermochemistry, which finds significant activation barriers to such reactions.

Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes

The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = E_quintet - E_singlet) of [Fe^IIL₆]²⁺ complexes (L = CNH, CO, NCH, NH₃, H₂O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH₃, H₂O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.

Mechanistic views and computational studies on transition-metal-catalyzed reductive coupling reactions

Transition-metal-catalyzed reductive coupling reactions have been considered as a powerful tool to convert two electrophiles into value-added products. Numerous related reports have shown the fascinating potential. Mechanistic studies, especially theoretical studies, can provide important implications for the design of novel reductive coupling reactions. In this review, we summarize the representative advancements in theoretical studies on transition-metal-catalyzed reductive coupling reactions and systematically elaborate the mechanisms for the key steps of reductive coupling reactions. The activation modes of electrophiles and the deep insights of selectivity generation are mechanistically discussed. In addition, the mechanism of the reduction of high-oxidation-state catalysts and further construction of new chemical bonds are also described in detail.

Accurate computed singlet-triplet energy differences for cobalt systems: implication for two-state reactivity

Accurate singlet-triplet energy differences for cobalt and rhodium complexes were calculated by using several wave function methods, such as MRCISD, CASPT2, CCSD(T) and BCCD(T). Relaxed energy differences were obtained by considering the singlet and triplet complexes, each at the minimum of their potential energy surfaces. Active spaces for multireference calculations were carefully checked to provide accurate results. The considered systems are built by increasing progressively the first coordination sphere around the metal. We included in our set two CpCoX complexes (Cp = cyclopentadienyl, X = alkenyl ligand), which have been suggested as intermediates in cycloaddition reactions. Indeed, cobalt systems have been used for more than a decade as active species in this kind of transformations, for which a two-state reactivity has been proposed. Most of the considered systems display a triplet ground state. However, in the case of a reaction intermediate, while a triplet ground state was predicted on the basis of Density Functional Theory results, our calculations suggest a singlet ground state. This stems from the competition between the exchange term (stabilising the triplet) and the accessibility of an intramolecular coordination (stabilising the singlet). This finding has an impact on the general mechanism of the cycloaddition reaction. Analogous rhodium systems were also studied and, as expected, they have a larger tendency to electron pairing than cobalt species.

The role of the intermediate triplet state in iron-catalyzed multi-state C-H activation

Efficient activation and functionalization of the C-H bond under mild conditions are of a great interest in chemical synthesis. We investigate the previously proposed spin-accelerated activation of the C(sp²)-H bond by a Fe(II)-based catalyst to clarify the role of the intermediate triplet state in the reaction mechanism. High-level electronic structure calculations on a small model of a catalytic system utilizing the coupled cluster with the single, double, and perturbative triple excitations [CCSD(T)] are used to select the density functional for the full-size model. Our analysis indicates that the previously proposed two-state quintet-singlet reaction pathway is unlikely to be efficient due to a very weak spin-orbit coupling between these two spin states. We propose a more favorable multi-state quintet-triplet-singlet reaction pathway and discuss the importance of the intermediate triplet state. This triplet state facilitates a spin-accelerated reaction mechanism by strongly coupling to both quintet and singlet states. Our calculations show that the C-H bond activation through the proposed quintet-triplet-singlet reaction pathway is more thermodynamically favorable than the single-state quintet and two-state singlet-quintet mechanisms.

Benchmarking CASPT3 Vertical Excitation Energies

Based on 280 reference vertical transition energies of various natures (singlet, triplet, valence, Rydberg, n → π∗, π → π∗, and double excitations) extracted from the QUEST database, we assess the accuracy of third-order multireference perturbation theory, CASPT3, in the context of molecular excited states. When one applies the disputable ionization- potential-electron-affinity (IPEA) shift, we show that CASPT3 provides a similar accuracy as its second-order counterpart, CASPT2, with the same mean absolute error of 0.11 eV. However, as already reported, we also observe that the accuracy of CASPT3 is almost insensitive to the IPEA shift, irrespective of the transition type and system size, with a small reduction of the mean absolute error to 0.09 eV when the IPEA shift is switched off.

Minimal Active Space for Diradicals Using Multistate Density Functional Theory

This work explores the electronic structure as well as the reactivity of singlet diradicals, making use of multistate density functional theory (MSDFT). In particular, we show that a minimal active space of two electrons in two orbitals is adequate to treat the relative energies of the singlet and triplet adiabatic ground state as well as the first singlet excited state in many cases. This is plausible because dynamic correlation is included in the first place in the optimization of orbitals in each determinant state via block-localized Kohn–Sham density functional theory. In addition, molecular fragment, i.e., block-localized Kohn–Sham orbitals, are optimized separately for each determinant, providing a variational diabatic representation of valence bond-like states, which are subsequently used in nonorthogonal state interactions (NOSIs). The computational procedure and its performance are illustrated on some prototypical diradical species. It is shown that NOSI calculations in MSDFT can be used to model bond dissociation and hydrogen-atom transfer reactions, employing a minimal number of configuration state functions as the basis states. For p- and s-types of diradicals, the closed-shell diradicals are found to be more reactive than the open-shell ones due to a larger diabatic coupling with the final product state. Such a diabatic representation may be useful to define reaction coordinates for electron transfer, proton transfer and coupled electron and proton transfer reactions in condensed-phase simulations.

FCIQMC-Tailored Distinguishable Cluster Approach: Open-Shell Systems

A recently proposed tailored approach based on the distinguishable cluster method and the stochastic FCI solver, FCIQMC [J. Chem. Theory Comput. 2020, 16, 5621], is extended to open-shell molecular systems. The method is employed to calculate spin gaps of various Fe(II) complexes, including a Fe(II) porphyrin model system. Both distinguishable cluster and fully relaxed CASSCF natural orbitals were used in this work as reference for the subsequent tailored distinguishable cluster calculations. The distinguishable cluster natural orbitals occupation numbers were also used as an aid to the selection of the active space. The effect of the active space sizes and of the explicit correlation correction (F12) onto the predicted spin gaps is investigated. The tailored distinguishable cluster with singles and doubles yields consistently more accurate results compared to the tailored coupled cluster with singles and doubles.

Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches

Hydrated transition metal ions are prototypical systems that can be used to model properties of transition metals in complex chemical environments. These seemingly simple systems present challenges for computational chemistry and are thus crucial in evaluations of quantum chemical methods for spin-state and redox energetics. In this work, we explore the applicability of the domain-based pair natural orbital implementation of coupled cluster (DLPNO-CC) theory to the calculation of ionization energies and redox potentials for hydrated ions of all first transition row (3d) metals in the 2+/3+ oxidation states, in connection with various solvation approaches. In terms of model definition, we investigate the construction of a minimally explicitly hydrated quantum cluster with a first and second hydration layer. We report on the convergence with respect to the coupled cluster expansion and the PNO space, as well as on the role of perturbative triple excitations. A recent implementation of the conductor-like polarizable continuum model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent redox potentials at the coupled cluster level. Our results establish conditions for the convergence of DLPNO-CCSD(T) energetics and stress the absolute necessity to explicitly consider the second solvation sphere even when CPCM is used. The achievable accuracy for redox potentials of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore, multilayer approaches that combine a higher-level DLPNO-CCSD(T) description of the first solvation sphere with a lower-level description of the second solvation layer are investigated. The present work establishes optimal and transferable methodological choices for employing DLPNO-based coupled cluster theory, the associated CPCM implementation, and cost-efficient multilayer derivatives of the approach for open-shell transition metal systems in complex environments.

Role of Magnetization on Catalytic Pathways of Non-Oxidative Methane Activation on Neutral Iron Carbide Clusters

Methane has emerged as a promising fuel due to its abundance and clean combustion properties. It is also a raw material for various value added chemicals. However, the conversion of...

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.

Improved Spin-State Energy Differences of Fe(II) Molecular and Crystalline Complexes via the Hubbard U-Corrected Density

We recently showed that the DFT+U approach with a linear-response U yields adiabatic energy differences biased toward high spin [Mariano et al. J. Chem. Theory Comput. 2020, 16, 6755-6762]. Such bias is removed here by employing a density-corrected DFT approach where the PBE functional is evaluated on the Hubbard U-corrected density. The adiabatic energy differences of six Fe(II) molecular complexes computed using this approach, named PBE[U] here, are in excellent agreement with coupled cluster-corrected CASPT2 values for both weak- and strong-field ligands resulting in a mean absolute error (MAE) of 0.44 eV, smaller than that of the recently proposed Hartree-Fock density-corrected DFT (1.22 eV) and any other tested functional, including the best performer TPSSh (0.49 eV). We take advantage of the computational efficiency of this approach and compute the adiabatic energy differences of five molecular crystals using PBE[U] with periodic boundary conditions. The results show, again, an excellent agreement (MAE = 0.07 eV) with experimentally extracted values and a superior performance compared with the best performers M06-L (MAE = 0.08 eV) and TPSSh (MAE = 0.31 eV) computed on molecular fragments.

Zeolites at the Molecular Level: What Can Be Learned from Molecular Modeling

This review puts the development of molecular modeling methods in the context of their applications to zeolitic active sites. We attempt to highlight the utmost necessity of close cooperation between theory and experiment, resulting both in advances in computational methods and in progress in experimental techniques.

Spin-Orbit Coupling Changes the Identity of the Hyper-Open-Shell Ground State of Ce+, and the Bond Dissociation Energy of CeH+ Proves to Be Challenging for Theory

Cerium (Ce) plays important roles in catalysis. Its position in the sixth period of the periodic table leads to spin-orbit coupling (SOC) and other open-shell effects that make the quantum mechanical calculation of cerium compounds challenging. In this work, we investigated the low-lying spin states of Ce⁺ and the bond energy of CeH⁺, both by multiconfigurational methods, in particular, SA-CASSCF, MC-PDFT, CASPT2, XMS-PDFT, and XMS-CASPT2, and by single-configurational methods, namely, Hartree-Fock theory and unrestricted Kohn-Sham density functional theory with 34 choices of the exchange-correlation functional. We found that only CASPT2, XMS-CASPT2, and SA-CASSCF (among the five multiconfigurational methods) and GAM, HCTH, SOGGA11, and OreLYP (among the 35 single-configuration methods) successfully predict that the SOC-free ground spin state of Ce⁺ is a doublet state, and CASPT2 and GAM give the most accurate multireference and single-reference calculations, respectively, of the excitation energy of the first SOC-free excited state for Ce⁺. We calculated that the ground doublet state of Ce⁺ is an intra-atomic hyper-open-shell state. We calculated the spin-orbit energy (E^SO) of Ce⁺ by the five multiconfigurational methods and found that E^SO calculated by CASPT2 is the closest to the experimental value. Taking advantage of the availability of an experimental D₀ for CeH⁺ as a way to provide a unique test of theory, we showed that all the multiconfigurational methods overestimate D₀ by at least 246 meV (5.7 kcal/mol), and only three functionals, namely, SOGGA, MN15, and GAM, have an error of D₀ that is less than 200 meV (5 kcal/mol).

Theoretical modeling of the singlet-triplet spin transition in different Ni(II)-diketo-pyrphyrin-based metal-ligand octahedral complexes

The structural stability, charge transfer effects and strength of the spin-orbit couplings in different Ni(ii)-ligand complexes have been studied at the DFT (B3LYP and CAM-B3LYP) and coupled cluster (DLPNO-CCSD(T)) levels of theory. Accordingly, two different, porphyrin- and diketo-pyrphyrin-based four-coordination macrocycles as planar ligands as well as pyridine (or pyrrole) and mesylate anion molecular groups as vertical ligands were considered in order to build metal-organic complexes with octahedral coordination configurations. For each molecular system, the identification of equilibrium geometries and the intersystem crossing (the minimum energy crossing) points between the potential energy surfaces of the singlet and triplet spin states is followed by computing the spin-orbit couplings between the two spin states. Structures, based on the diketo-pyrphyrin macrocycle as the planar ligand, show stronger six-coordination metal-organic complexes due to the extra electrostatic interaction between the positively charged central metal cation and the negatively charged vertical ligands. The results also show that the magnitude of the spin-orbit coupling is influenced by the atomic positions of deprotonations of the ligands, and implicitly the direction of the charge transfer between the ligand and the central metal ion.

Spin-state energetics of metallocenes: How do best wave function and density functional theory results compare with the experimental data?

We benchmark the accuracy of quantum-chemical methods, including wave function theory methods [coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation-theory (CASPT2, NEVPT2) and internally contracted multireference configuration interaction (MRCI)] and 30 density functional theory (DFT) approximations, in reproducing the spin-state splittings of metallocenes. The reference values of the electronic energy differences are derived from the experimental spin-crossover enthalpy for manganocene and the spectral data of singlet-triplet transitions for ruthenocene, ferrocene, and cobaltocenium. For ferrocene and cobaltocenium we revise the previous experimental interpretations regarding the lowest triplet energy; our argument is based on the comparison with the lowest singlet excitation energy and herein reported, carefully determined absorption spectrum of ferrocene. When deriving vertical energies from the experimental band maxima, we go beyond the routine vertical energy approximation by introducing vibronic corrections based on simulated vibrational envelopes. The benchmarking result confirms the high accuracy of the CCSD(T) method (in particular, for UCCSD(T) based on Hartree-Fock orbitals we find for our dataset: maximum error 0.12 eV, weighted mean absolute error 0.07 eV, weighted mean signed error 0.01 eV). The high accuracy of the single-reference method is corroborated by the analysis of a multiconfigurational character of the complete active space wave function for the triplet state of ferrocene. On the DFT side, our results confirm the non-universality problem with approximate functionals. The present study is an important step toward establishing an extensive and representative benchmark set of experiment-derived spin-state energetics for transition metal complexes.

Improvement of d-d interactions in density functional tight binding for transition metal ions with a ligand field model: assessment of a DFTB3+U model on nickel coordination compounds

To improve the description of interactions among the localized d, f electrons in transition metals, we have introduced a ligand-field motivated contribution into the Density Functional Tight Binding (DFTB) model. Referred to as DFTB3+U, the approach treats the d, f electron repulsions with rotationally invariant orbital-orbital interactions and a Hartree-Fock model; this represents a major conceptual improvement over the original DFTB3 approach, which treats the d, f-shell interactions in a highly averaged fashion without orbital level of description. The DFTB3+U approach is tested using a series of nickel compounds that feature Ni(ii) and Ni(iii) oxidation states. By using parameters developed with the original DFTB3 Hamiltonian and empirical +U parameters (F0/2/4 Slater integrals), we observe that the DFTB3+U model indeed provides substantial improvements over the original DFTB3 model for a number of properties of the nickel compounds, including the population and spin polarization of the d-shell, nature of the frontier orbitals, ligand field splitting and the energy different between low and high spin states at OPBE optimized structures. This proof-of-concept study suggests that with self-consistent parameterization of the electronic and +U parameters, the DFTB3+U model can develop into a promising model that can be used to efficiently study reactive events involving transition metals ion condensed phase systems. The methodology can be integrated with other approximate QM methods as well, such as the extended tight binding (xTB) approach.

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.