Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Prediction and Measurement of Resonant and Nonresonant Shake Effects in the Core-Level X-ray Emission Spectra of 3d0 Transition Metal Compounds

Shake effects, resulting from sudden core potential changes during photoexcitation, are well-known in X-ray photoelectron spectroscopy (XPS) and often produce satellite peaks due to many-body excitations. It has been thought, however, that they are negligible in core-to-core X-ray emission spectroscopy (CTC-XES), where the difference in core-hole potentials upon radiative decay is rather small. We demonstrate that shake effects are significant in Kα XES from 3d transition metal systems with nominally zero valence electrons. We show that valence level shake satellites are amplified via interference due to a resonance between the 2p3/2-hole (Kα1) plus valence level shake state and the 2p1/2-hole (Kα2) state. Additionally, while the Kα2 shake satellite is indeed predicted to be weak, we observe it experimentally, providing further independent verification of our model. This prediction includes a detailed analysis of 2p to 1s Kα XES using density functional theory (DFT)-augmented multiplet ligand field theory (MLFT).

Ichor: A Python library for computational chemistry data management and machine learning force field development

We present ichor, an open-source Python library that simplifies data management in computational chemistry and streamlines machine learning force field development. Ichor implements many easily extensible file management tools, in addition to a lazy file reading system, allowing efficient management of hundreds of thousands of computational chemistry files. Data from calculations can be readily stored into databases for easy sharing and post-processing. Raw data can be directly processed by ichor to create machine learning-ready datasets. In addition to powerful data-related capabilities, ichor provides interfaces to popular workload management software employed by High Performance Computing clusters, making for effortless submission of thousands of separate calculations with only a single line of Python code. Furthermore, a simple-to-use command line interface has been implemented through a series of menu systems to further increase accessibility and efficiency of common important ichor tasks. Finally, ichor implements general tools for visualization and analysis of datasets and tools for measuring machine-learning model quality both on test set data and in simulations. With the current functionalities, ichor can serve as an end-to-end data procurement, data management, and analysis solution for machine-learning force-field development.

Performance of quantum chemistry methods for a benchmark set of spin-state energetics derived from experimental data of 17 transition metal complexes (SSE17)

Accurate prediction of spin-state energetics for transition metal (TM) complexes is a compelling problem in applied quantum chemistry, with enormous implications for modeling catalytic reaction mechanisms and computational discovery of materials. Computed spin-state energetics are strongly method-dependent and credible reference data are scarce, making it difficult to conduct conclusive computational studies of open-shell TM systems. Here, we present a novel benchmark set of first-row TM spin-state energetics, which is derived from experimental data of 17 complexes containing FeII, FeIII, CoII, CoIII, MnII, and NiII with chemically diverse ligands. The estimates of adiabatic or vertical spin-state splittings, which are obtained from spin crossover enthalpies or energies of spin-forbidden absorption bands, suitably back-corrected for the vibrational and environmental effects, are employed as reference values for benchmarking density functional theory (DFT) and wave function methods. The results demonstrate a high accuracy of the coupled-cluster CCSD(T) method, which features the mean absolute error (MAE) of 1.5 kcal mol-1 and maximum error of -3.5 kcal mol-1, and outperforms all the tested multireference methods: CASPT2, MRCI+Q, CASPT2/CC and CASPT2+δMRCI. Switching from Hartree-Fock to Kohn-Sham orbitals is not found to consistently improve the CCSD(T) accuracy. The best performing DFT methods are double-hybrids (PWPB95-D3(BJ), B2PLYP-D3(BJ)) with the MAEs below 3 kcal mol-1 and maximum errors within 6 kcal mol-1, whereas the DFT methods so far recommended for spin states (e.g., B3LYP*-D3(BJ) and TPSSh-D3(BJ)) are found to perform much worse with the MAEs of 5-7 kcal mol-1 and maximum errors beyond 10 kcal mol-1. This work is the first such extensive benchmark study of quantum chemistry methods for TM spin-state energetics making use of experimental reference data. The results are relevant for the proper choice of methods to characterize TM systems in computational catalysis and (bio)inorganic chemistry, and may also stimulate new developments in quantum-chemical or machine learning approaches.

Similarity between oxygen evolution in photosystem II and oxygen reduction in cytochrome c oxidase via proton coupled electron transfers. A unified view of the oxygenic life from four electron oxidation-reduction reactions

Basic concepts and theoretical foundations of broken symmetry (BS) and post BS methods for strongly correlated electron systems (SCES) such as electron-transfer (ET) diradical, multi-center polyradicals with spin frustration are described systematically to elucidate structures, bonding and reactivity of the high-valent transition metal oxo bonds in metalloenzymes: photosystem II (PSII) and cytochrome c oxidase (CcO). BS hybrid DFT (HDFT) and DLPNO coupled-cluster (CC) SD(T0) computations are performed to elucidate electronic and spin states of CaMn4Ox cluster in the key step for oxygen evolution, namely S4 [S3 with Mn(IV) = O + Tyr161-O radical] state of PSII and PM [Fe(IV) = O + HO-Cu(II) + Tyr161-O radical] step for oxygen reduction in CcO. The cycle of water oxidation catalyzed by the CaMn4Ox cluster in PSII and the cycle of oxygen reduction catalyzed by the CuA-Fea-Fea3-CuB cluster in CcO are examined on the theoretical grounds, elucidating similar concerted and/or stepwise proton transfer coupled electron transfer (PT-ET) processes for the four-electron oxidation in PSII and four-electron reduction in CcO. Interplay between theory and experiments have revealed that three electrons in the metal sites and one electron in tyrosine radical site are characteristic for PT-ET in these biological redox reaction systems, indicating no necessity of harmful Mn(V) = O and Fe(V) = O bonds with strong oxyl-radical character. Implications of the computational results are discussed in relation to design of artificial systems consisted of earth abundant transition metals for water oxidation.

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average error of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Finally, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

Orbital entanglement and the double d-shell effect in binary transition metal molecules

Accurate modeling of transition metal-containing compounds is of great interest due to their wide-ranging and significant applications. These systems present several challenges from an electronic structure perspective, including significant multi-reference characters and many chemically relevant orbitals. A further complication arises from the so-called double d-shell effect, which is known to cause a myriad of issues in the treatment of first-row transition metals with both single- and multi-reference methods. While this effect has been well documented for several decades, a comprehensive understanding of its consequences and underlying causes is still evolving. Here, we characterize the second d-shell effect by analyzing the information entropy of correlated wavefunctions in a periodic series of 3d and 4d transition metal molecular hydrides and oxides. These quantum information techniques provide unique insight into the nuanced electronic structure of these species and are powerful tools for the study of weak and strong correlations in the transition metal d manifold.

Gas-Phase Characterization of Redox-Active Manganese Complexes

Manganese complexes exhibit a rich redox chemistry, usually accompanied by structural reorganization during the redox processes often followed by ligand dissociation or association. The push-pull ligand 2,6-diguanidylpyridine (dgpy) stabilizes manganese in the oxidation states +II, +III, and + IV in the complexes [Mn(dgpy)2]n+ (n = 2-4) without change in the coordination sphere in the condensed phase [Heinze et al., Inorganic Chemistry, 2022, 61, 14616]. In the condensed phase, the manganese(IV) complex is a very strong oxidant. In the present work, we investigate the stability and redox activity of the MnIV complex and its counterion (PF6-) adducts in the gas phase, using two modified 3D Paul ion trap mass spectrometers. Six different cationic species of the type [Mnx(dgpy)2(PF6)y]n+ (x = II, III, IV, y = 0-3, n = 1-3) could be observed for the three oxidation states MnIV, MnIII, and MnII, of which one observed complex also contains a reduced dgpy ligand. MnII species showed the highest relative stability in collision induced dissociation and UV/vis photo dissociation experiments. The lowest stability is observed in the presence of one or more counterions, which correlates to a lower total charge n+. Gas phase UV/vis spectra show similar features as the condensed phase spectra only differing in relative band intensities. The strongly oxidizing MnIV complex reacts with triethylamine (NEt3) in the gas phase to give MnIII, while MnIII species show little reactivity toward NEt3.

Ground and Excited Electronic Structure Analysis of FeH with Correlated Wave Function Theory and Density Functional Approximations

FeH is one of the most challenging diatomic molecules to study under electronic structure theory. Here, we have successfully studied 22 electronic states of FeH using ab initio multireference configuration interaction (MRCI), Davidson-corrected MRCI (MRCI+Q), and coupled cluster singles, doubles, and perturbative triples [CCSD(T)] levels of theory. We report their potential energy curves (PECs), excitation energies, dissociation energies, equilibrium electronic configurations, and a series of spectroscopic constants with the use of augmented triple-ζ, quadruple-ζ, and quintuple-ζ quality correlation consistent basis sets. The scalar relativistic effects and active space and core electron correlation contribution on the properties of FeH are also explored. The use of a large CASSCF active space that includes 4s, 4p, 3d, and 4d orbitals of Fe and the 1s of H is critical for producing accurate full PECs with proper dissociations and predicting the exact order of the electronic states. Our findings are in harmony with the experimental results available in the literature and will serve as reference values for future studies of FeH. Furthermore, with the use of PECs, the total internal partition function sum (TIPS) of FeH was calculated across a range of temperatures. Finally, we exploited the single-reference nature of the a6Δ of FeH and its ionized product FeH+ (X5Δ) to evaluate the associated density functional theory (DFT) errors on their dissociation energies and spectroscopic parameters.

Quantum algorithms for scientific computing

Quantum computing promises to provide the next step up in computational power for diverse application areas. In this review, we examine the science behind the quantum hype, and the breakthroughs required to achieve true quantum advantage in real world applications. Areas that are likely to have the greatest impact on high performance computing (HPC) include simulation of quantum systems, optimization, and machine learning. We draw our examples from electronic structure calculations and computational fluid dynamics which account for a large fraction of current scientific and engineering use of HPC. Potential challenges include encoding and decoding classical data for quantum devices, and mismatched clock speeds between classical and quantum processors. Even a modest quantum enhancement to current classical techniques would have far-reaching impacts in areas such as weather forecasting, aerospace engineering, and the design of 'green' materials for sustainable development. This requires significant effort from the computational science, engineering and quantum computing communities working together.

Metallicious: Automated Force-Field Parameterization of Covalently Bound Metals for Supramolecular Structures

Metal ions play a central, functional, and structural role in many molecular structures, from small catalysts to metal-organic frameworks (MOFs) and proteins. Computational studies of these systems typically employ classical or quantum mechanical approaches or a combination of both. Among classical models, only the covalent metal model reproduces both geometries and charge transfer effects but requires time-consuming parameterization, especially for supramolecular systems containing repetitive units. To streamline this process, we introduce metallicious, a Python tool designed for efficient force-field parameterization of supramolecular structures. Metallicious has been tested on diverse systems including supramolecular cages, knots, and MOFs. Our benchmarks demonstrate that parameters accurately reproduce the reference properties obtained from quantum calculations and crystal structures. Molecular dynamics simulations of the generated structures consistently yield stable simulations in explicit solvent, in contrast to similar simulations performed with nonbonded and cationic dummy models. Overall, metallicious facilitates the atomistic modeling of supramolecular systems, key for understanding their dynamic properties and host-guest interactions. The tool is freely available on GitHub (https://github.com/duartegroup/metallicious).

Machine learning-guided strategies for reaction conditions design and optimization

This review surveys the recent advances and challenges in predicting and optimizing reaction conditions using machine learning techniques. The paper emphasizes the importance of acquiring and processing large and diverse datasets of chemical reactions, and the use of both global and local models to guide the design of synthetic processes. Global models exploit the information from comprehensive databases to suggest general reaction conditions for new reactions, while local models fine-tune the specific parameters for a given reaction family to improve yield and selectivity. The paper also identifies the current limitations and opportunities in this field, such as the data quality and availability, and the integration of high-throughput experimentation. The paper demonstrates how the combination of chemical engineering, data science, and ML algorithms can enhance the efficiency and effectiveness of reaction conditions design, and enable novel discoveries in synthetic chemistry.

Wavefunction theory and density functional theory analysis of ground and excited electronic states of TaB and WB

Several low-lying electronic states of TaB and WB molecules were studied using ab initio multireference configuration interaction (MRCI), Davidson corrected MRCI (MRCI+Q), and coupled cluster singles doubles and perturbative triples [CCSD(T)] methods. Their full potential energy curves (PECs), equilibrium electron configurations, equilibrium bond distances (res), dissociation energies (Des), excitation energies (Tes), harmonic vibrational frequencies (ωes), and anharmonicities (ωexes) are reported. The MRCI dipole moment curves (DMCs) of the first 5 electronic states of both TaB and WB are also reported and the equilibrium dipole moment (μ) values are compared with the CCSD(T) μ values. The most stable 13Π (1σ22σ23σ11π3) and 15Δ (1σ22σ23σ11π21δ1) electronic states of TaB lie close in energy with ∼62 kcal mol-1De with respect to the Ta(4F) + B(2P) asymptote. However, spin-orbit coupling effects make the 15Δ0+ state the true ground state of TaB. The ground electronic state of WB (16Π) has the 1σ22σ13σ11π31δ2 electron configuration and is followed by the excited 16Σ+ and 14Δ states. Finally, the MRCI De, re, ωe, and ωexe values of the 13Π state of TaB and 16Π and 14Δ states of WB are used to assess the density functional theory (DFT) errors on a series of exchange-correlation functionals that span multiple-rungs of the Jacob's ladder of density functional approximations (DFA).

Ligand engineering of tetra N-heterocyclic carbenes for boosting catalytic aziridination

A comprehensive computational study on the underlying reactivity of iron tetra-NHC complexes for C2 + N1 aziridination catalysis is presented. A library of 18 unique iron tetra-NHC complexes was constructed, and a computational screening was performed on the reaction barriers associated with the rate-determining step (formation of an open chain radical intermediate). Thermodynamic barriers were computed along with a variety of steric and electronic properties, including the percentage of buried volume, orbital energies and ETS-NOCV analysis, which were used to identify key characteristics related to reactivity. The analysis performed in this study successfully identified key differences in tetracarbenes, such as linking groups (BMe2 or CH2) and the identity of the NHC groups (imidazole, imidazoline or benzimidazole) in terms of sterics, electronics and thermodynamics. Additionally, we have proposed two reaction pathways based on electronic structure arguments for the formation of the key open-chain radical intermediate. The first reaction pathway proceeds through a σ-hole channel where the Fe(IV)-imide intermediate evolves into Fe(III)-imidyl radical through electron donation into the antibonding σ* orbital, while the second involves a Fe(III)-imidyl radical formed through a π-hole channel (donation into π*). These pathways are consistent with the isoelectronic iron(IV)-oxo species for hydrogen atom abstraction mechanisms and they can be used as descriptors of the rate-determining step of the aziridination reaction.

A Caged Neutral 17-Valence-Electron Iron(I) Radical [Fe(CO)2(Cl)(P((CH2)10)3P)]•: Synthetic, Structural, Spectroscopic, Redox, and Computational Studies

UV irradiation of yellow CH2Cl2 solutions of trans-Fe(CO)3(P((CH2)10)3P) (2a) and PMe3 (10 equiv) gives, in addition to the previously reported dibridgehead diphosphine P((CH2)10)3P (46%), a green paramagnetic complex that crystallography shows to be the trigonal-bipyramidal iron(I) radical trans-[Fe(CO)2(Cl)(P((CH2)10)3P)]• (1a•; 31% after workup). This is a rare example of an isolable species of the formula [Fe(CO)4-n(L)n(X)]• (n = 0-3, L = two-electron-donor ligand; X = one-electron-donor ligand). Analogous precursors with longer P(CH2)nP segments (n = 12, 14, 16, 18) give only the demetalated diphosphines, and a rationale is proposed. The magnetic susceptibility of 1a•, assayed by Evans' method and SQUID measurements, indicates a spin (S) of 1/2. Cyclic voltammetry shows that 1a• undergoes a partially reversible one-electron oxidation, but no facile reduction. The UV-visible, EPR, and 57Fe Mössbauer spectra are analyzed in detail. Complex 2a is similarly studied, and, despite the extra valence electron, exhibits a comparable oxidation potential (ΔE1/2 ≤ 0.04 V). The crystal structure shows a cage conformation, solvation level, disorder motif, and unit cell parameters essentially identical to those of 1a•. DFT calculations provide much insight regarding the structural, redox, and spectroscopic properties.

Catalytic mechanisms for As(III) oxidation by H2O2 over TiO2 surfaces, and effects of support, vacancy and photoirradiation

As(III) is much more toxic than As(V) while shows apparently lower affinity at minerals surfaces. Oxidation of As(III) to As(V) by H2O2 over anatase surface provides an attractive avenue for pollution control, and the chemocatalytic and photocatalytic mechanisms are unraveled by means of the DFT + D3 approach. Impacts of anatase as support, O2c/O3c vacancy, photoirradiation are addressed as well. As(III) oxidation under various reaction conditions leads to As(V) through dual electron transfers, while energy barriers differ substantially and decline as 1.80 (direct oxidation) > 1.35 (anatase as support) > 1.24 (O3c vacancy) > 0.50 (chemocatalysis) > 0.28 (photocatalysis) ≥ 0.26 (O2c vacancy) eV. Anatase as support promotes the reaction through bonding with H2O2/As(OH)3 and electron transfers, and its close participation during chemocatalysis produces the TiOOH active site that causes As(III) oxidation to proceed facilely under ambient circumstances. TiOOH exists in two forms (monodentate and bidentate mononuclear) and is critical for chemocatalysis, while its destruction for O3c vacancy exhibits strongly adverse effects to As(III) oxidation. Photoirradiation readily generates the OH• radicals, and corresponding mechanism is plausible while less preferred than the newly posed mechanism based on the Ti(H2O2) active site. Synergism among a number of surface atoms conduces to the superior activity for O2c vacancy and photocatalysis. Results provide a comprehensive understanding for As(III) oxidation to As(V) by H2O2, and facilitate catalysts design for As(III) oxidation that alleviates environmental pollution.

How Doping Regulates As(III) Adsorption at TiO2 Surfaces: A DFT + U Study

The efficient adsorption and removal of As(III), which is highly toxic, remains difficult. TiO2 shows promise in this field, though the process needs improvement. Herein, how doping regulates As(OH)3 adsorption over TiO2 surfaces is comprehensively investigated by means of the DFT + D3 approach. Doping creates the bidentate mononuclear (Ce doping at the Ti5c site), tridentate (N, S doping at the O2c site), and other new adsorption structures. The extent of structural perturbation correlates with the atomic radius when doping the Ti site (Ce >> Fe, Mn, V >> B), while it correlates with the likelihood of forming more bonds when doping the O site (N > S > F). Doping the O2c, O3c rather than the Ti5c site is more effective in enhancing As(OH)3 adsorption and also causes more structural perturbation and diversity. Similar to the scenario of pristine surfaces, the bidentate binuclear complexes with two Ti-OAs bonds are often the most preferred, except for B doping at the Ti5c site, S doping at the O2c site, and B doping at the O3c site of rutile (110) and Ce, B doping at the Ti5c site, N, S doping at the O2c site, and N, S, B doping at the O3c site of anatase (101). Doping significantly regulates the As(OH)3 adsorption efficacy, and the adsorption energies reach -4.17, -4.13, and -4.67 eV for Mn doping at the Ti5c site and N doping at the O2c and O3c sites of rutile (110) and -1.99, -2.29, and -2.24 eV for Ce doping at the Ti5c site and N doping at the O2c and O3c sites of anatase (101), respectively. As(OH)3 adsorption and removal are crystal-dependent and become apparently more efficient for rutile vs. anatase, whether doped at the Ti5c, O2c, or O3c site. The auto-oxidation of As(III) occurs when the As centers interact directly with the TiO2 surface, and this occurs more frequently for rutile rather than anatase. The multidentate adsorption of As(OH)3 causes electron back-donation and As(V) re-reduction to As(IV). The regulatory effects of doping during As(III) adsorption and the critical roles played by crystal control are further unraveled at the molecular level. Significant insights are provided for As(III) pollution management via the adsorption and rational design of efficient scavengers.

Metal complexes featuring a quinazoline schiff base ligand and glycine: synthesis, characterization, DFT and molecular docking analyses revealing their potent antibacterial, anti-helicobacter pylori, and Anti-COVID-19 activities

Mixed ligand complexes of manganese(II), cobalt(II), copper(II), and cadmium(II)with an innovative Schiff base ligand denoted as (L1), 4-(2-((1E,2E)-1-(2-(p-tolyl)hydrazineylidene)propan-2-ylidene)hydrazineyl), served as the principal ligand, while glycine (L2) was employed as secondary ligand were successfully effectively characterized through a comprehensive set of analyses, including Elemental analysis, UV-Visible, FT-IR, Mass spectra, and conductometric measurements. Density functional theory (DFT) computations were executed to discern the enduring electronic arrangement, the energy gap, dipole moment and chemical hardness of the hybrid ligand assemblies. The proposed geometry for the complexes is a distorted octahedral structure. The antimicrobial efficacy of these compounds was assessed against a range of bacterial and fungal strains. Notably, these complexes exhibited promising antimicrobial activities, with the cadmium (II) complex demonstrating superior efficacy towards all tested organisms. These compounds were also examined for their antibiotic properties against H. pylori to explore their broader medical potential. The Schiff base ligand and its corresponding metal complexes displayed substantial potential as an antibiotic against H. pylori. Additionally, the antitumor potential of the synthesized complexes was assessed against MCF-7 (Breast carcinoma) cells-the Cu (II) complex demonstrated superior activity with the lowest IC50 value compared to cisplatin. Moreover, it exhibited reduced cytotoxicity towards normal cells (VERO cells) compared to cisplatin, establishing it as the most potent compound in the study. Furthermore, molecular docking was explored of the Schiff base ligand and its corresponding cadmium(II) complex. The analysis of the docking study yielded valuable structural insights that can be effectively utilized in conducting inhibition studies for example against COVID-19. This comprehensive study highlights these synthesized compounds' multifaceted applications and promising bioactive properties.

The Use of Effective Core Potentials with Multiconfiguration Pair-Density Functional Theory

The reliable and accurate prediction of chemical properties is a key goal in quantum chemistry. Transition-metal-containing complexes can often pose difficulties to quantum mechanical methods for multiple reasons, including many electron configurations contributing to the overall electronic description of the system and the large number of electrons significantly increasing the amount of computational resources required. Often, multiconfigurational electronic structure methods are employed for such systems, and the cost of these calculations can be reduced by the use of an effective core potential (ECP). In this work, we explore both theoretical considerations and performances of ECPs applied in the context of multiconfiguration pair-density functional theory (MC-PDFT). A mixed-basis set approach is used, using ECP basis sets for transition metals and all-electron basis sets for nonmetal atoms. We illustrate the effects that an ECP has on the key parameters used in the computation of MC-PDFT energies, and we explore how ECPs affect the prediction of physical observables for chemical systems. The dissociation curve for a metal dimer was explored, and ionization energies for transition metal-containing diatomic systems were computed and compared to experimental values. In general, we find that ECP approaches employed with MC-PDFT are able to predict ionization energies with improved accuracy compared to traditional Kohn-Sham density functional theory approaches.

Automated prediction of ground state spin for transition metal complexes

Exploiting crystallographic data repositories for large-scale quantum chemical computations requires the rapid and accurate extraction of the molecular structure, charge and spin from the crystallographic information file. Here, we develop a general approach to assign the ground state spin of transition metal complexes, in complement to our previous efforts on determining metal oxidation states and bond order within the cell2mol software. Starting from a database of 31k transition metal complexes extracted from the Cambridge Structural Database with cell2mol, we construct the TM-GSspin dataset, which contains 2063 mononuclear first row transition metal complexes and their computed ground state spins. TM-GSspin is highly diverse in terms of metals, metal oxidation states, coordination geometries, and coordination sphere compositions. Based on TM-GSspin, we identify correlations between structural and electronic features of the complexes and their ground state spins to develop a rule-based spin state assignment model. Leveraging this knowledge, we construct interpretable descriptors and build a statistical model achieving 98% cross-validated accuracy in predicting the ground state spin across the board. Our approach provides a practical way to determine the ground state spin of transition metal complexes directly from crystal structures without additional computations, thus enabling the automated use of crystallographic data for large-scale computations involving transition metal complexes.

Overcoming the Pitfalls of Computing Reaction Selectivity from Ensembles of Transition States

The prediction of reaction selectivity is a challenging task for computational chemistry, not only because many molecules adopt multiple conformations but also due to the exponential relationship between effective activation energies and rate constants. To account for molecular flexibility, an increasing number of methods exist that generate conformational ensembles of transition state (TS) structures. Typically, these TS ensembles are Boltzmann weighted and used to compute selectivity assuming Curtin-Hammett conditions. This strategy, however, can lead to erroneous predictions if the appropriate filtering of the conformer ensembles is not conducted. Here, we demonstrate how any possible selectivity can be obtained by processing the same sets of TS ensembles for a model reaction. To address the burdensome filtering task in a consistent and automated way, we introduce marc, a tool for the modular analysis of representative conformers that aids in avoiding human errors while minimizing the number of reoptimization computations needed to obtain correct reaction selectivity.

Cobalt-Mediated Photochemical C-H Arylation of Pyrroles

Precious metal complexes remain ubiquitous in photoredox catalysis (PRC) despite concerted efforts to find more earth-abundant catalysts and replacements based on 3d metals in particular. Most otherwise plausible 3d metal complexes are assumed to be unsuitable due to short-lived excited states, which has led researchers to prioritize the pursuit of longer excited-state lifetimes through careful molecular design. However, we report herein that the C-H arylation of pyrroles and related substrates (which are benchmark reactions for assessing the efficacy of photoredox catalysts) can be achieved using a simple and readily accessible octahedral bis(diiminopyridine) cobalt complex, [1-Co](PF6)2. Notably, [1-Co]2+ efficiently functionalizes both chloro- and bromoarene substrates despite the short excited-state lifetime of the key photoexcited intermediate *[1-Co]2+ (8 ps). We present herein the scope of this C-H arylation protocol and provide mechanistic insights derived from detailed spectroscopic and computational studies. These indicate that, despite its transient existence, reduction of *[1-Co]2+ is facilitated via pre-assembly with the NEt3 reductant, highlighting an alternative strategy for the future development of 3d metal-catalyzed PRC.

Coupled-cluster treatment of complex open-shell systems: the case of single-molecule magnets

We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

Understanding Cooperativity in Homo- and Heterometallic Complexes: From Basic Concepts to Design

Cooperative effects have attracted considerable attention in recent years. These effects are ubiquitous in chemistry and biology and can govern interactions of proteins with other biomolecules, mechanisms of supramolecular recognition and polymerization, catalysis, assembly of compounds on surfaces, and physical properties such as magnetic, electronic or optical properties, e. g. Consequently, the understanding of cooperative effects can lead to a structure-property relation that can pave the way to future applications in various research areas; however, with regard to cooperative effects in homo- and heterometallic complexes, we still are at the beginning of understanding. Nevertheless, concepts to describe cooperativity of metal centers as well as methodologies to investigate and model these effects have emerged over the last years. This concept article gives an overview of these existing concepts, approaches, and strategies to understand cooperative effects in homo- and heterometallic complexes. Special emphasis is put on concepts to define cooperative effects, their quantification, as well as methods to investigate cooperative effects.

DFT Struggles to Predict the Energy Landscape for Iron Pyridine Diimine-Catalyzed [2 + 2] Cycloaddition of Alkenes: Insights into the Problem and Alternative Solutions

In recent years, noninnocent pyridine diimine (PDI) complexes featuring first-row transition metals have emerged as prominent catalysts, demonstrating efficacy in a diverse range of vital organometallic transformations. However, the inherent complexity of the fundamental reactivity paradigm in these systems arises from the presence of a noninnocent ligand and the multispin feasibility of 3d metals. While density functional theory (DFT) has been widely used to unravel mechanistic insights, its limitations as a single-reference method can potentially misrepresent spin-state energetics, compromising our understanding of these intricate systems. In this study, we employ extensive high-level ab initio state averaged-complete active space self-consistent field/N-electron valence state perturbation theory (SA-CASSCF/NEVPT2) calculations in combination with DFT to investigate an iron-PDI-catalyzed [2 + 2] cycloaddition reaction of alkenes. The transformation proceeds through two major steps: oxidative cyclization and reductive elimination. Contrary to the predictions of DFT calculations, which suggest two-state reactivity in the reaction and identify reductive elimination as the turnover-limiting step, SA-CASSCF/NEVPT2-corrected results unequivocally establish a single-state reactivity scenario with oxidative cyclization as the turnover-limiting step. SA-CASSCF/NEVPT2-based insights into electronic ground states and electron distribution elucidate the intriguing interactions between the PDI ligand and the iron center, revealing the highly multiconfigurational nature of these species and providing a precise depiction of metal-ligand cooperativity throughout the transformation. A comparative assessment of several widely recognized DFT functionals against SA-CASSCF/NEVPT2-corrected data indicates that single-point energy calculations using the modern density functional MN15 on TPSSh geometries offer the most reliable density functional methodology, in scenarios where SA-CASSCF/NEVPT2 computational cost is a consideration.

Explainable chemical artificial intelligence from accurate machine learning of real-space chemical descriptors

Machine-learned computational chemistry has led to a paradoxical situation in which molecular properties can be accurately predicted, but they are difficult to interpret. Explainable AI (XAI) tools can be used to analyze complex models, but they are highly dependent on the AI technique and the origin of the reference data. Alternatively, interpretable real-space tools can be employed directly, but they are often expensive to compute. To address this dilemma between explainability and accuracy, we developed SchNet4AIM, a SchNet-based architecture capable of dealing with local one-body (atomic) and two-body (interatomic) descriptors. The performance of SchNet4AIM is tested by predicting a wide collection of real-space quantities ranging from atomic charges and delocalization indices to pairwise interaction energies. The accuracy and speed of SchNet4AIM breaks the bottleneck that has prevented the use of real-space chemical descriptors in complex systems. We show that the group delocalization indices, arising from our physically rigorous atomistic predictions, provide reliable indicators of supramolecular binding events, thus contributing to the development of Explainable Chemical Artificial Intelligence (XCAI) models.

In-silico-assisted derivatization of triarylboranes for the catalytic reductive functionalization of aniline-derived amino acids and peptides with H2

Cheminformatics-based machine learning (ML) has been employed to determine optimal reaction conditions, including catalyst structures, in the field of synthetic chemistry. However, such ML-focused strategies have remained largely unexplored in the context of catalytic molecular transformations using Lewis-acidic main-group elements, probably due to the absence of a candidate library and effective guidelines (parameters) for the prediction of the activity of main-group elements. Here, the construction of a triarylborane library and its application to an ML-assisted approach for the catalytic reductive alkylation of aniline-derived amino acids and C-terminal-protected peptides with aldehydes and H2 is reported. A combined theoretical and experimental approach identified the optimal borane, i.e., B(2,3,5,6-Cl4-C6H)(2,6-F2-3,5-(CF3)2-C6H)2, which exhibits remarkable functional-group compatibility toward aniline derivatives in the presence of 4-methyltetrahydropyran. The present catalytic system generates H2O as the sole byproduct.

A human-machine interface for automatic exploration of chemical reaction networks

Autonomous reaction network exploration algorithms offer a systematic approach to explore mechanisms of complex chemical processes. However, the resulting reaction networks are so vast that an exploration of all potentially accessible intermediates is computationally too demanding. This renders brute-force explorations unfeasible, while explorations with completely pre-defined intermediates or hard-wired chemical constraints, such as element-specific coordination numbers, are not flexible enough for complex chemical systems. Here, we introduce a STEERING WHEEL to guide an otherwise unbiased automated exploration. The STEERING WHEEL algorithm is intuitive, generally applicable, and enables one to focus on specific regions of an emerging network. It also allows for guiding automated data generation in the context of mechanism exploration, catalyst design, and other chemical optimization challenges. The algorithm is demonstrated for reaction mechanism elucidation of transition metal catalysts. We highlight how to explore catalytic cycles in a systematic and reproducible way. The exploration objectives are fully adjustable, allowing one to harness the STEERING WHEEL for both structure-specific (accurate) calculations as well as for broad high-throughput screening of possible reaction intermediates.

A review on transition metal oxides based photocatalysts for degradation of synthetic organic pollutants

This review provides insight into the current research trend in transition metal oxides (TMOs)-based photocatalysis in removing the organic colouring matters from water. For easy understanding, the research progress has been presented in four generations according to the catalyst composition and mode of application, viz: single component TMOs (the first-generation), doped TMOs/binary TMOs/doped binary TMOs (the second-generation), inactive/active support-immobilized TMOs (the third-generation), and ternary/quaternary compositions (the fourth-generation). The first two generations represent suspended catalysts, the third generation is supported catalysts, and the fourth generation can be suspended or supported. The review provides an elaborated comparison between suspended and supported catalysts, their general/specific requirements, key factors controlling degradation, and the methodologies for performance evaluation. All the plausible fundamental and advanced dye degradation mechanisms involved in each generation of catalysts were demonstrated. The existing challenges in TMOs-based photocatalysis and how the researchers approach the hitch to resolve it effectively are discussed. Future research trends are also presented.

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.

Ligand Field Sensitive Spin Acceleration in the Iron-Catalyzed [2 + 2] Cycloaddition of Unactivated Alkenes and Dienes

Redox-active pyridine(diimine) (PDI) iron catalysts promote the reversible [2 + 2] cycloaddition of alkenes and dienes to cyclobutane derivatives that have applications ranging from fuels to chemically recyclable polymers. Metallacycles were identified as key intermediates, and spin crossover from the singlet to the triplet surface was calculated to facilitate the reductive coupling step responsible for the formation of the four-membered ring. In this work, a series of sterically and electronically differentiated PDI ligands was studied for the [2 + 2] cycloaddition of ethylene and butadiene to vinylcyclobutane. Kinetic studies revealed that the fastest and slowest turnover were observed with equally electron-deficient supporting ligands that either feature phenyl-substituted imine carbon atoms (MeBPDI) or a pyrazine core (MePZDI). While the oxidative cyclization was comparatively slow for both catalysts, the rate of reductive coupling─determined by stoichiometric 13C2H4 labeling studies─correlated with the turnover frequencies. Two-state density functional theory studies and the distinct electronic structures of related (iPrBPDI) and (iPrPZDI) iron methyl complexes revealed significantly different ligand field strengths due to either diminished ligand σ-donation (MeBPDI) or promoted metal π-backbonding (MePZDI). Spin acceleration, leading to fast reductive coupling and catalytic turnover, was promoted in the case of the weaker ligand field and depends on both the nature and position of the electron-withdrawing group. This study provides strong evidence for the role of two-state reactivity in C(sp3)-C(sp3) bond formation and insights on how ligand design either promotes or inhibits spin acceleration in earth-abundant metal catalysis.

A Synergistic Bimetallic Ti/Co-Catalyzed Isomerization of Epoxides to Allylic Alcohols Enabled by Two-State Reactivity

Isomerization of epoxides into versatile allylic alcohols is an atom-economical synthetic method to afford vicinal bifunctional groups. Comprehensive density functional theory (DFT) calculations were carried out to elucidate the complex mechanism of a bimetallic Ti/Co-catalyzed selective isomerization of epoxides to allyl alcohols by examining several possible pathways. Our results suggest a possible mechanism involving (1) radical-type epoxide ring opening catalyzed by Cp2Ti(III)Cl leading to a Ti(IV)-bound β-alkyl radical, (2) hydrogen-atom transfer (HAT) catalyzed by the Co(II) catalyst to form the Ti(IV)-enolate and Co(III)-H intermediate, (3) protonation to give the alcohols, and (4) proton abstraction to form the Co(I) species followed by electron transfer to regenerate the active Co(II) and Ti(III) species. Moreover, bimetallic catalysis and two-state reactivity enable the key rate-determining HAT step. Furthermore, a subtle balance between dispersion-driven bimetallic processes and entropy-driven monometallic processes determines the most favorable pathway, among which the monometallic process is energetically more favorable in all steps except the vital hydrogen-atom transfer step. Our study should provide an in-depth mechanistic understanding of bimetallic catalysis.

Electronic Structure and Ligand Effects on the Activation and Cleavage of N2 on a Molybdenum Center

Dinitrogen fixation under ambient conditions remains a challenge in the field of catalytic chemistry due to the inertness of N2. Nitrogenases and heterogeneous solid catalysts have displayed remarkable performance in the catalytic conversion of dinitrogen to ammonia. By introduction of molybdenum centers in molecular complexes, one of the most azophilic metals of the transitional metal series, moderate ammonia yields have been attained. Here, we present a combined multiconfigurational/density functional theory study that addresses how ligand fields of different strengths affect the binding and activation of dinitrogen on molybdenum atoms. First, we explored with MRCI computations the diatomic Mo-N and triatomic Mo-N2 molecular systems. Then, we performed a systematic examination on the stabilization effects introduced by external NH3 ligands, before we explore model neutral and charged complexes with different types of ligands (H2O, NH3, and PH3) and their consequences on the N2 binding and activation.

Computationally Guided Ligand Discovery from Compound Libraries and Discovery of a New Class of Ligands for Ni-Catalyzed Cross-Electrophile Coupling of Challenging Quinoline Halides

Although screening technology has heavily impacted the fields of metal catalysis and drug discovery, its application to the discovery of new catalyst classes has been limited. The diversity of on- and off-cycle pathways, combined with incomplete mechanistic understanding, means that screens of potential new ligands have thus far been guided by intuitive analysis of the metal binding potential. This has resulted in the discovery of new classes of ligands, but the low hit rates have limited the use of this strategy because large screens require considerable cost and effort. Here, we demonstrate a method to identify promising screening directions via simple and scalable computational and linear regression tools that leads to a substantial improvement in hit rate, enabling the use of smaller screens to find new ligands. The application of this approach to a particular example of Ni-catalyzed cross-electrophile coupling of aryl halides with alkyl halides revealed a previously overlooked trend: reactions with more electron-poor amidine ligands result in a higher yield. Focused screens utilizing this trend were more successful than serendipity-based screening and led to the discovery of two new types of ligands, pyridyl oxadiazoles and pyridyl oximes. These ligands are especially effective for couplings of bromo- and chloroquinolines and isoquinolines, where they are now the state of the art. The simplicity of these models with parameters derived from metal-free ligand structures should make this approach scalable and widely accessible.

In Situ UV-Vis-NIR Absorption Spectroscopy and Catalysis

This review highlights in situ UV-vis-NIR range absorption spectroscopy in catalysis. A variety of experimental techniques identifying reaction mechanisms, kinetics, and structural properties are discussed. Stopped flow techniques, use of laser pulses, and use of experimental perturbations are demonstrated for in situ studies of enzymatic, homogeneous, heterogeneous, and photocatalysis. They access different time scales and are applicable to different reaction systems and catalyst types. In photocatalysis, femto- and nanosecond resolved measurements through transient absorption are discussed for tracking excited states. UV-vis-NIR absorption spectroscopies for structural characterization are demonstrated especially for Cu and Fe exchanged zeolites and metalloenzymes. This requires combining different spectroscopies. Combining magnetic circular dichroism and resonance Raman spectroscopy is especially powerful. A multitude of phenomena can be tracked on transition metal catalysts on various supports, including changes in oxidation state, adsorptions, reactions, support interactions, surface plasmon resonances, and band gaps. Measurements of oxidation states, oxygen vacancies, and band gaps are shown on heterogeneous catalysts, especially for electrocatalysis. UV-vis-NIR absorption is burdened by broad absorption bands. Advanced analysis techniques enable the tracking of coking reactions on acid zeolites despite convoluted spectra. The value of UV-vis-NIR absorption spectroscopy to catalyst characterization and mechanistic investigation is clear but could be expanded.

MACE: Automated Assessment of Stereochemistry of Transition Metal Complexes and Its Applications in Computational Catalysis

Computational chemistry pipelines typically commence with geometry generation, well-established for organic compounds but presenting a considerable challenge for transition metal complexes. This paper introduces MACE, an automated computational workflow for converting chemist SMILES/MOL representations of the ligands and the metal center to 3D coordinates for all feasible stereochemical configurations for mononuclear octahedral and square planar complexes directly suitable for quantum chemical computations and implementation in high-throughput computational chemistry workflows. The workflow is validated through a structural screening of a data set of transition metal complexes extracted from the Cambridge Structural Database. To further illustrate the power and capabilities of MACE, we present the results of a model DFT study on the hemilability of pincer ligands in Ru, Fe, and Mn complexes, which highlights the utility of the workflow for both focused mechanistic studies and larger-scale high-throughput pipelines.

High-spin Co3+ in cobalt oxyhydroxide for efficient water oxidation

Cobalt oxyhydroxide (CoOOH) is a promising catalytic material for oxygen evolution reaction (OER). In the traditional CoOOH structure, Co3+ exhibits a low-spin state configuration ([Formula: see text]), with electron transfer occurring in face-to-face [Formula: see text] orbitals. In this work, we report the successful synthesis of high-spin state Co3+ CoOOH structure, by introducing coordinatively unsaturated Co atoms. As compared to the low-spin state CoOOH, electron transfer in the high-spin state CoOOH occurs in apex-to-apex [Formula: see text] orbitals, which exhibits faster electron transfer ability. As a result, the high-spin state CoOOH performs superior OER activity with an overpotential of 226 mV at 10 mA cm-2, which is 148 mV lower than that of the low-spin state CoOOH. This work emphasizes the effect of the spin state of Co3+ on OER activity of CoOOH based electrocatalysts for water splitting, and thus provides a new strategy for designing highly efficient electrocatalysts.

Spin-Adapted Externally Contracted Multireference Configuration Interaction Method Based on Selected Reference Configurations

As one kind of approximation of the full configuration interaction solution, the selected configuration interaction (sCI) methods have been shown to be valuable for large active spaces. However, the inclusion of dynamic correlation beyond large active spaces is necessary for more quantitative results. Since the sCI wave function can provide a compact reference for multireference methods, previously, we proposed an externally contracted multireference configuration interaction method using the sCI reference reconstructed from the density matrix renormalization group wave function [J. Chem. Theory Comput. 2018, 14, 4747-4755]. The DMRG2sCI-EC-MRCI method is promising for dealing with more than 30 active orbitals and large basis sets. However, it suffers from two drawbacks: spin contamination and low efficiency when using Slater determinant bases. To solve these problems, in this work, we adopt configuration state function bases and introduce a new algorithm based on the hybrid of tree structure for convenient configuration space management and the graphical unitary group approach for efficient matrix element calculation. The test calculation of naphthalene shows that the spin-adapted version could achieve a speed-up of 6.0 compared with the previous version based on the Slater determinant. Examples of dinuclear copper(II) compound as well as Ln(III) and An(III) complexes show that the sCI-EC-MRCI can give quantitatively accurate results by including dynamic correlation over sCI for systems with large active spaces and basis sets.

Spin-Crossing in the (Z)-Selective Alkyne Semihydrogenation Mechanism Catalyzed by Mo3S4 Clusters: A Density Functional Theory Exploration

Semihydrogenation of internal alkynes catalyzed by the air-stable imidazolyl amino [Mo3S4Cl3(ImNH2)3]+ cluster selectively affords the (Z)-alkene under soft conditions in excellent yields. Experimental results suggest a sulfur-based mechanism with the formation of a dithiolene adduct through interaction of the alkyne with the bridging sulfur atoms. However, computational studies indicate that this mechanism is unable to explain the experimental outcome: mild reaction conditions, excellent selectivity toward the (Z)-isomer, and complete deuteration of the vinylic positions in the presence of CD3OD and CH3OD. An alternative mechanism that explains the experimental results is proposed. The reaction begins with the hydrogenation of two of the Mo3(μ3-S)(μ-S)3 bridging sulfurs to yield a bis(hydrosulfide) intermediate that performs two sequential hydrogen atom transfers (HAT) from the S-H groups to the alkyne. The first HAT occurs with a spin change from singlet to triplet. After the second HAT, the singlet state is recovered. Although the dithiolene adduct is more stable than the hydrosulfide species, the large energy required for the subsequent H2 addition makes the system evolve via the second alternative pathway to selectively render the (Z)-alkene with a lower overall activation barrier.

A DFT perspective on organometallic lanthanide chemistry

Computational studies of the coordination chemistry and bonding of lanthanides have grown in recent decades as the need for understanding the distinct physical, optical, and magnetic properties of these compounds increased. Density functional theory (DFT) methods offer a favorable balance of computational cost and accuracy in lanthanide chemistry and have helped to advance the discovery of novel oxidation states and electronic configurations. This Frontier article examines the scope and limitations of DFT in interpreting structural and spectroscopic data of low-valent lanthanide complexes, elucidating periodic trends, and predicting their properties and reactivity, presented through selected examples.

DFT-Guided Mechanistic Insights into Chemodivergence: A Mixed Explicit-Implicit Solvent Description to Dictate the Chemoselectivity

Herein we report a density functional theory (DFT)-guided mechanistic investigation of the nitrile reduction reaction, which exhibits a solvent-dependent chemodivergence. This study reveals an interesting mechanistic picture, highlighting the exact role of a protic solvent, isopropanol, in regulating the reaction outcome. The explicit solvent effect involving polar protic isopropanol favors imine metathesis by proton hopping through stepwise addition and elimination steps and thus produces a secondary amine as the final product. In contrast, the nonpolar solvent n-hexane is incapable of facilitating the proton migration and stops the solvent-assisted imine metathesis. As a result, only primary amines are obtained as the final product. This DFT study provides a recipe for the choice of solvents that can dictate chemoselectivity in product formation.

First-row transition metal for isocyanide-involving multicomponent reactions (IMCR)

First-row transition metal catalyzed transformations that are able to construct complex molecules from simple, readily obtainable feedstocks have become a keystone of modern synthetic organic chemistry. Particularly, the multicomponent reaction (MCR) involving carbon-carbon (C-C) as well as carbon-heteroatom (C-X) bond formation plays an essential role in many chemical conversions, and insurgencies in these reactions powerfully improve the overall synthetic efficiency. Recently, MCRs emerges rapidly because of its greener sides like eco-friendly nature, swift and straightforward execution, high atom/step economy, and construction of aimed product with lowest or no by-product, usually in quantitative yield. Curiously, the exceptional divalent carbon atoms of isocyanides make them predominantly useful components in multicomponent reactions. As a result of widespread research over the past few decades, numerous well-designed and effective procedures for the first-row TM-catalyzed MCR to afford the various entities have been reported. These aspects are summarized in this review article. A particular focus on comparative discussion of various first-row transition-metal catalyzed isocyanide-based multicomponent reactions through mechanistic details included in the review article.

Electrochemical interfacial catalysis in Co-based battery electrodes involving spin-polarized electron transfer

Interfacial catalysis occurs ubiquitously in electrochemical systems, such as batteries, fuel cells, and photocatalytic devices. Frequently, in such a system, the electrode material evolves dynamically at different operating voltages, and this electrochemically driven transformation usually dictates the catalytic reactivity of the material and ultimately the electrochemical performance of the device. Despite the importance of the process, comprehension of the underlying structural and compositional evolutions of the electrode material with direct visualization and quantification is still a significant challenge. In this work, we demonstrate a protocol for studying the dynamic evolution of the electrode material under electrochemical processes by integrating microscopic and spectroscopic analyses, operando magnetometry techniques, and density functional theory calculations. The presented methodology provides a real-time picture of the chemical, physical, and electronic structures of the material and its link to the electrochemical performance. Using Co(OH)2 as a prototype battery electrode and by monitoring the Co metal center under different applied voltages, we show that before a well-known catalytic reaction proceeds, an interfacial storage process occurs at the metallic Co nanoparticles/LiOH interface due to injection of spin-polarized electrons. Subsequently, the metallic Co nanoparticles act as catalytic activation centers and promote LiOH decomposition by transferring these interfacially residing electrons. Most intriguingly, at the LiOH decomposition potential, electronic structure of the metallic Co nanoparticles involving spin-polarized electrons transfer has been shown to exhibit a dynamic variation. This work illustrates a viable approach to access key information inside interfacial catalytic processes and provides useful insights in controlling complex interfaces for wide-ranging electrochemical systems.

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.

Enhanced immobilization of Arsenic(III) and Auto-oxidation to Arsenic(V) by titanium oxide (TiO2), due to Single-Atom vacancies and oxyanion formation

Pollution control of As(III), a naturally occurring carcinogen, has recently gained a global attention, while due to the dominance of neutral H3AsO3 over a wide pH range, As(III) immobilization by most minerals is not efficient as As(V) immobilization. TiO2 shows promise for controlling As(III) pollution, and herein, a comprehensive study about As(III) adsorption by TiO2 and oxyanion formation is conducted by means of DFT + D3 methods. Both anatase and rutile are effective for As(III) adsorption, while As(III) adsorption affinities differ significantly and are -1.48 and -3.79 eV for pristine surfaces, ascend to -3.85 and -5.08 eV for O vacancies, and further to -5.37 and -5.26 eV for Ti vacancies, respectively. The bidentate binuclear complexes dominate for pristine surfaces, and O vacancies prefer OAs insertion into TiO2 lattice, while for Ti vacancies, all As(III) centers are auto-oxidized to As(V). Ti-3d, O-2p or/and As-4p rather than other orbitals contribute significantly to As adsorption, and O and Ti vacancies promote adsorption through stronger orbital hybridization. The superior adsorption for Ti vacancies originates from As(V) formation instead of bonding interactions. The formation of As oxyanions, which may occur spontaneously at pristine surfaces and is greatly promoted by O and Ti vacancies, enhances As(III) adsorption pronouncedly and becomes a viable strategy for As(III) immobilization. H2AsO3- and HAsO32- dominate for pristine surfaces and O vacancies, and for Ti vacancies, H2AsO4- and HAsO42- dominate over anatase whereas AsO43- also makes an important contribution over rutile. Results rationalize experimental observations available, and provide significantly new insights about the migration, bioavailability and fate of As(III) over TiO2 surfaces that facilitate the exploration of scavengers for As and other pollutants.

Cobalt-based tripodal complexes as molecular catalysts for photocatalytic CO2 reduction

Construction of artificial photosynthetic systems including CO2 reduction is a promising pathway to produce carbon-neutral fuels and mitigate the greenhouse effect concurrently. However, the exploitation of earth-abundant catalysts for photocatalytic CO2 reduction remains a fundamental challenge, which can be assisted by a systematic summary focusing on a specific catalyst family. Cobalt-based complexes featuring tripodal ligands should merit more insightful discussion and summarization, as they are one of the most examined catalyst families for CO2 photoreduction. In this feature article, the key developments of cobalt-based tripodal complexes as molecular catalysts for light-driven CO2 reduction are discussed to offer an upcoming perspective, analyzing the present progress in electronic/steric tuning through ligand modification and dinuclear design to achieve a synergistic effect, as well as the bottlenecks for further development.

Quantum chemical calculations for reaction prediction in the development of synthetic methodologies

Quantum chemical calculations have been used in the development of synthetic methodologies to analyze the reaction mechanisms of the developed reactions. Their ability to estimate chemical reaction pathways, including transition state energies and connected equilibria, has led researchers to embrace their use in predicting unknown reactions. This perspective highlights strategies that leverage quantum chemical calculations for the prediction of reactions in the discovery of new methodologies. Selected examples demonstrate how computation has driven the development of unknown reactions, catalyst design, and the exploration of synthetic routes to complex molecules prior to often laborious, costly, and time-consuming experimental investigations.

How does multi-reference computation change the catalysis chemistry? DFT and CASPT2 studies of the Cu-catalysed coupling reactions between aryl iodides and β-diketones

The molecular mechanism of a Cu-catalysed coupling reaction was theoretically studied using density functional theory (DFT) and the complete active space self-consistent field method followed by the second-order perturbation theory (CASSCF/CASPT2) to investigate the effects of the strong electron correlation of the Cu centre on the reaction profile. Both DFT and CASSCF/CASPT2 calculations showed that the catalytic cycle proceeds via an oxidative addition (OA) reaction, followed by a reductive elimination (RE) reaction, where OA is the rate-determining step. Although the DFT-calculated activation energies of the OA and RE steps are highly dependent on the choice of functionals, the CASSCF/CASPT2 results are less affected by the choice of DFT-optimised geometries. Therefore, with a careful assessment based on the CASSCF/CASPT2 single-point energy evaluation, an optimal choice of the DFT geometry is of good qualitative use for energetics at the CASPT2 level of theory. Based on the changes in the electron populations of the 3d orbitals during the OA and RE steps, the characteristic features of the DFT-calculated electronic structure were qualitatively consistent with those calculated using the CASSCF method. Further electronic structure analysis by the natural orbital occupancy of the CASSCF wavefunction showed that the ground state is almost single-reference in this system and the strong electron correlation effect of the Cu centre can be dealt with using the MP2 or CCSD method, too. However, the slightly smaller occupation numbers of the 3dπ orbital in the course of reactions suggested that the electron correlation effect of the Cu(III) centre appears through the interaction between the 3dπ orbital and the C-I antibonding σ* orbital in the OA step, and between the 3dπ orbital and the Cu-C antibonding σ* orbital in the RE step.

Gas-phase and solid-state electronic structure analysis and DFT benchmarking of HfCO

Ab initio multi-reference configuration interaction (MRCI) and coupled cluster singles doubles and perturbative triples [CCSD(T)] levels of theory were used to study ground and excited electronic states of HfCO. We report potential energy curves, dissociation energies (De), excitation energies, harmonic vibrational frequencies, and chemical bonding patterns of HfCO. The 3Σ- ground state of HfCO has an 1σ22σ21π2 electron configuration and a ∼30 kcal mol-1 dissociation energy with respect to its lowest-energy fragments Hf(3F) + CO(X1Σ+). We further evaluated the De of its isovalent HfCX (X = S, Se, Te, Po) series and observed that they increase linearly from the lighter HfCO to the heavier HfCPo with the dipole moment of the CX ligand. The same linear relationship was observed for TiCX and ZrCX. We utilized the CCSD(T) benchmark values of De, excitation energy, and ionization energy (IE) values to evaluate density functional theory (DFT) errors with 23 exchange-correlation functionals spanning GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, and double-hybrid functional families. The global GGA hybrid B3LYP and range-separated hybrid ωB97X performed well at representing the ground state properties of HfCO (i.e., De and IE). Finally, we extended our DFT analysis to the interaction of a CO molecule with a Hf surface and observed that the surface chemisorption energy and the gas-phase molecular dissociation energy are very similar for some DFAs but not others, suggesting moderate transferability of the benchmarks on these molecules to the solid state.

Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC

Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.

The mechanism of visible light-induced C-C cross-coupling by Csp3-H bond activation

Csp3-C cross-coupling by activating Csp3-H bonds is a dream reaction for the chemical community, and visible light-induced transition metal-catalysis under mild reaction conditions is considered a powerful tool to achieve it. Advancement of this research area is still in its infancy because of the chemical and technical complexity of this catalysis. Mechanistic studies illuminating the operative reaction pathways can rationalize the increasing amount of experimental catalysis data and provide the knowledge allowing faster and rational advances in the field. This goal requires complementary experimental and theoretical mechanistic studies, as each of them is unfit to clarify the operative mechanisms alone. In this tutorial review we summarize representative experimental and computational mechanistic studies, highlighting weaknesses, strengths, and synergies between the two approaches.