Theoretical study of the Fe(phen)(2)(NCS)(2) spin-crossover complex with reparametrized density functionals

Accurate and interpretable representation of correlated electronic structure via Tensor Product Selected CI

The task of computing wavefunctions that are accurate, yet simple enough mathematical objects to use for reasoning, has long been a challenge in quantum chemistry. The difficulty in drawing physical conclusions from a wavefunction is often related to the generally large number of configurations with similar weights. In Tensor Product Selected Configuration Interaction (TPSCI), we use a locally correlated tensor product state basis, which has the effect of concentrating the weight of a state onto a smaller number of physically interpretable degrees of freedom. In this paper, we apply TPSCI to a series of three molecular systems ranging in separability, one of which is the first application of TPSCI to an open-shell bimetallic system. For each of these systems, we obtain accurate solutions to large active spaces, and analyze the resulting wavefunctions through a series of different approaches including (i) direct inspection of the TPS basis coefficients, (ii) construction of Bloch effective Hamiltonians, and (iii) computation of cluster correlation functions.

Computational Screening of Putative Catalyst Transition Metal Complexes as Guests in a Ga4L612- Nanocage

Metal-organic cages form well-defined microenvironments that can enhance the catalytic proficiency of encapsulated transition metal complexes (TMCs). We introduce a screening protocol to efficiently identify TMCs that are promising candidates for encapsulation in the Ga4L612- nanocage. We obtain TMCs from the Cambridge Structural Database with geometric and electronic characteristics amenable to encapsulation and mine the text of associated manuscripts to curate TMCs with documented catalytic functionality. By docking candidate TMCs inside the nanocage cavity and carrying out electronic structure calculations, we identify a subset of successfully optimized candidates (TMC-34) and observe that encapsulated guests occupy an average of 60% of the cavity volume, in line with previous observations. Notably, some guests occupy as much as 72% of the cavity as a result of linker rotation. Encapsulation has a universal effect on the electrostatic potential (ESP), systematically decreasing the ESP at the metal center of each TMC in the TMC-34 data set, while minimally altering TMC metal partial charges. Collectively these observations support geometry-based screening of potential guests and suggest that encapsulation in Ga4L612- cages could electrostatically stabilize diverse cationic or electropositive intermediates. We highlight candidate guests with associated known reactivity and solubility most amenable for encapsulation in experimental follow-up studies.

Creation of spin switching in graphene oxide-based hybrid film materials with an anionic Fe(III) 5Cl-salicyaldehyde-thiosemicarbazone complex

The present article describes the synthesis of hybrid composite film materials formed during the self-assembly process through non-covalent interactions of graphene oxide (GO) nanosheets with salt 1, represented by an anionic spin-crossover complex [FeIII(5Cl-thsa)2]- (5Cl-thsa - 5-chlorosalicylaldehyde thiosemicarbazone) and the organic tetraethylammonium cation [Et4N]+. The insertion of the salt 1 molecules into the interlayer space of GO nanosheets with the subsequent formation of a hybrid material GO-1 was observed. The film of the hybrid material GO-1 was characterized by scanning electron and confocal laser microscopy, EDX and XPS analysis, IR, Raman and 57Fe Mössbauer spectroscopy, dc magnetic measurements, and powder X-ray diffraction. Comparison of the magnetic properties of salt 1 and a film of the hybrid material GO-1 demonstrated a significant influence of the GO nanosheets matrix on the completeness of spin transition and showed a slight shift of the hysteresis loop by 1 K in the temperature range of 200-230 K. DFT calculations showed an important role of the organic cation [Et4N]+ in the process of adsorption of the spin-crossover anion [FeIII(5Cl-thsa)2]- on the GO nanosheet surface.

Ab initio calculation of magnetic anisotropy and thermal spin transition in the variable temperature crystal conformations of [Co(terpy)2]2

The structure-property correlation of [Co(terpy)2]2+, which shows a spin crossover at 270 K, has been computationally investigated based on its variable temperature crystal structures. Among the employed DFT functionals, only the re-parametrized hybrid B3LYP* functional could describe the correct spin transition temperature. Explicit consideration of metal-ligand sigma bonding with dynamic electron correlation is found to be necessary for an accurate determination of the SCO temperature with multi-reference calculations. The metal-ligand axial bond distances are found to be the most significant internal coordinates in deciding SCO. A small structural change along the axial distance causes a change in the t2g orbital splitting pattern and a reorientation of the magnetization axes at the SCO temperature. The complex shows an unusual triaxial magnetic anisotropy, with an easy axis of magnetization developing at higher temperatures. The strong coupling of low-frequency wagging motion of the two terpyridine ligands with the spin states of the complex provides an effective pathway for the relaxation of magnetization, resulting in a small magnetic anisotropy barrier.

Fano resonance in molecular junctions of spin crossover complexes

In this paper, we introduce a novel molecular switch paradigm that integrates spin crossover complexes with the Fano resonance effect. Specifically, by performing density-functional theory calculations, the feasibility of achieving Fano resonance using spin crossover complexes is demonstrated in our designed molecular junctions using the complex {Fe[H2B(pz)2]2[Bp(bipy)]} [pz = 1-pyrazolyl, Bp(bipy) = bis(phenylethynyl)(2,2'-bipyridine)]. It is further revealed that the Fano resonance, particularly the Fano dip, is most prominent in the junction with cobalt tips among all the schemes, together with the spin-filtering effect. Most importantly, this junction of cobalt tips is able to exhibit three distinct conductance states, which are controlled by the modulation of Fano resonance due to the spin-state transition of the complex and the applied gate voltage. Such a molecular switch paradigm holds potential for applications in logic gates, memory units, sensors, thermoelectrics, and beyond.

Modeling Fe(II) Complexes Using Neural Networks

We report a Fe(II) data set of more than 23000 conformers in both low-spin (LS) and high-spin (HS) states. This data set was generated to develop a neural network model that is capable of predicting the energy and the energy splitting as a function of the conformation of a Fe(II) organometallic complex. In order to achieve this, we propose a type of scaled electronic embedding to cover the long-range interactions implicitly in our neural network describing the Fe(II) organometallic complexes. For the total energy prediction, the lowest MAE is 0.037 eV, while the lowest MAE of the splitting energy is 0.030 eV. Compared to baseline models, which only incorporate short-range interactions, our scaled electronic embeddings improve the accuracy by over 70% for the prediction of the total energy and the splitting energy. With regard to semiempirical methods, our proposed models reduce the MAE, with respect to these methods, by 2 orders of magnitude.

Machine Learning Prediction of the Experimental Transition Temperature of Fe(II) Spin-Crossover Complexes

Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. We leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T1/2) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T1/2 values. We then compare ML T1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T1/2 values with a Pearson's r of 0.82. In contrast, DFA-predicted T1/2 values have large errors and show no correlation to estimated T1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.

Theoretical and DFT Study of Atypical Pentanuclear [(iPr3P)Ni]5Hn (n = 4, 6, 8) Clusters: What are the Rules?

The structure, bonding, and properties of a series of atypical pentanuclear nickel hydride clusters supported by electron-rich iPr3P of the type [(iPr3P)Ni]5Hn (n = 4, 6, 8; H4, H6, H8) and their anionic models where iPr3P are substituted by H- (H4', H6', H8') were investigated by density functional theory (DFT) calculations. All clusters were calculated to adopt a similar square pyramidal core geometry. Calculations indicate singlet ground states with small singlet-triplet gaps for H4 and H6, similar to previously reported experimental values. Molecular orbital theory description clusters were investigated using the simplified model complexes [HNi]5Hn5- (n = 4, 6, 8; H4', H6', H8'). The results show that there are three skeletal electron pairs (SEPs) in H4'. The addition of two molecules of H2 to form H6' and H8' results in the partial or full occupation of two degenerate MOs (e* set) that give two SEPs and one SEP, respectively. Indeed, the occupation of these low-lying weakly antibonding orbitals governs the multielectron chemistry available for these clusters and plays a role in their unique reactivity.

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.

Spin state and magnetic coupling in polynuclear Ni(II) complexes from density functional theory: is there an optimal amount of Fock exchange?

Reliable prediction of the ground-state spin and magnetic coupling constants in transition-metal complexes is a well-known challenge for density functional theory (DFT). One popular strategy for addressing this long-standing issue involves the modification of the fraction of Fock exchange in a hybrid functional. Here we explore the viability of this approach using three polynuclear metal-organic complexes based on a Ni4O4 cubane motif, having different ground state spin values (S = 0, 2, 4) owing to the use of different ligands. We systematically search for an optimum fraction of Fock exchange, across various global, range-separated, and double hybrid functionals. We find that for all functionals tested, at best there only exists a very narrow range of Fock exchange fractions which results in a correct prediction of the ground-state spin for all three complexes. The useful range is functional dependent, but general trends can be identified. Typically, at least two similar systems must be used in order to determine both an upper and lower limit of the optimal range. This is likely owing to conflicting demands of minimizing delocalization errors, which typically requires a higher percentage of Fock exchange, and addressing static correlation, which typically requires a lower one. Furthermore, we find that within the optimal range of Fock exchange, the sign and relative magnitude of Ni-Ni magnetic coupling constants are reasonably well reproduced, but there is still room for quantitative improvement in the prediction. Thus, the prediction of spin state and magnetic coupling in polynuclear complexes remains an ongoing challenge for DFT.

How complex-surface interactions modulate the spin transition of Fe(II) SCO complexes supported on metallic surfaces?

The deposition of a prototypical spin-crossover [Fe(phen)2(NCS)2] complex on Au(111), Cu(111) and Ag(111) surfaces has been investigated by means of periodic DFT+U calculations, with the aim of understanding how different metallic surfaces affect the spin state switching. Our results show that adsorption is metal- and spin-dependent, with different preferred adsorption sites for the different surfaces and spin states. For the three considered surfaces adsorption energies are larger in the LS state than in the HS one, which increases the transition enthalpy by 58.7 kJ mol-1 for Cu(111), 14.6 kJ mol-1 for Au(111) and 9.6 kJ mol-1 for Ag(111) with respect to the free molecule. There is a clear correlation between this effect and the extent of the complex-surface interaction, which can be established from adsorption energies, surface-complex distances and charge density difference plots as: Cu(111) > Au(111) > Ag(111). Therefore, a stronger interaction with the surface produces a larger energy difference between two spin states, making the spin transition less probable to occur. Finally, our calculations show that it would be possible to probe the spin-state of the deposited molecules from the STM images, in line with the recent experimental results.

Is the surface of Hofmann-like spin-crossover {Fe(pz)[Pt(CN)4]} the same as its bulk?

Temperature dependent X-ray photoemission spectroscopy (XPS) has been employed to examine the spin-crossover (SCO) transition in the nanocrystals of 3D Hoffman-like {Fe(pz)[Pt(CN)4]}. Consistent with the existing literature, the temperature-dependent variations in the Fe 2p core-level spectrum provide unambiguous evidence of the spin-state transition in this SCO complex. One of the many possible reasons behind a lack of discernible temperature-driven shifts in the binding energies of both the N 1s core-level components could be the immunity of its HS electronic configuration to thermal fluctuations. The high-spin fraction versus temperature plot, extrapolated from the XPS measurements, reveals that the surface of the nanocrystals of {Fe(pz)[Pt(CN)4]} is in the high-spin state at room temperature, rendering it promising for room-temperature spintronics and quantum information science applications.

Spin-state energetics and magnetic anisotropy in penta-coordinated Fe(III) complexes with different axial and equatorial ligand environments

The penta-coordinated trigonal-bi-pyramidal (TBP) Fe(III) complex (PMe₂Ph)₂FeCl₃ shows a reduced magnetic anisotropy in its intermediate-spin (IS) state as compared to its methyl-analog (PMe₃)₂Fe(III)Cl₃. In this work, the ligand environment in (PMe₂Ph)₂FeCl₃ is systematically altered by replacing the axial -P with -N and -As, the equatorial -Cl with other halides, and the axial methyl group with an acetyl group. This has resulted in a series of Fe(III) TBP complexes modelled in their IS and high-spin (HS) states. Lighter ligands -N and -F stabilize the complex in the HS state, while the magnetically anisotropic IS state is stabilized by -P and -As at the axial site, and -Cl, -Br, and -I at the equatorial site. Larger magnetic anisotropies appear for complexes with nearly degenerate ground electronic states that are well separated from the higher excited states. This requirement, largely controlled by the d-orbital splitting pattern due to the changing ligand field, is achieved with a certain combination of axial and equatorial ligands, such as -P and -Br, -As and -Br, and -As and -I. In most cases, the acetyl group at the axial site enhances the magnetic anisotropy compared to its methyl counterpart. In contrast, the presence of -I at the equatorial site compromises the uniaxial type of anisotropy of the Fe(III) complex leading to an enhanced rate of quantum tunneling of magnetization.

Electronic Structure and Geometry of the Lowest 2LMCT State of Fe(III) Potential Fluorescent Emitters†

Metal complexes with a 3d⁶ electron count are emerging as an alternative to 4d⁶-based photosensitizers, emitters, or photoredox catalysts. In recent years, several Fe(II) potential emitters have been proposed, based on strongly donating ligand sets. Those tend to facilitate oxidation to their 3d⁵ species, whose photophysics is based on low-lying ligand-to-metal charge-transfer (LMCT) states. The geometry and electronic structure of ²LMCT states are unveiled in this work.

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew-Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. The quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

First crystal structure of an Fe(III) anionic complex based on a pyruvic acid thiosemicarbazone ligand with Li+: synthesis, features of magnetic behavior and theoretical analysis

The iron(III) anionic complex based on a pyruvic acid thiosemicarbazone ligand with the lithium cation Li[Fe^III(thpy)₂]·3H₂O (1) has been synthesized and characterized by FTIR spectroscopy, powder and single crystal X-ray diffraction, direct current magnetic susceptibility measurements, and ⁵⁷Fe Mössbauer spectroscopy. Moreover, the molecular structure of the [Fe(thpy)₂]^- anion has been determined for the first time. The [Fe(thpy)₂]^- units in the triclinic P1̄ lattice of 1 are assembled into layers parallel to the bc plane. The Li⁺ cations and water molecules are located between the layers and the structure is stabilized by hydrogen bonding. The [Fe(thpy)₂]^- anions form interconnected dimer pairs through hydrogen bonds and short contacts with Fe⋯Fe separation of 6.7861(4) Å. According to dc magnetic measurements, compound 1 demonstrates an incipient spin-crossover transition from the LS (S = 1/2) to the HS (S = 5/2) state above 250 K. The Bleaney-Bowers equation for a model of an isolated LS dimer with a mean-field correction was applied to fit the experimental data of magnetic susceptibility dependence on temperature in the temperature range of 2-250 K. The intra-dimer J₁ = -1.79(1) K and inter-dimer J₂ = -0.24(3) K antiferromagnetic coupling constants were defined. The analysis of the ⁵⁷Fe Mössbauer spectra at 80 K and 296 K confirms the presence of the shortened distances between the iron nuclei. Moreover, the influence of the lithium cation on the stabilization of the LS state was shown for the [Fe(thpy)₂]^- anion. BS-DFT calculations for the optimized structure of two isolated [Fe(thpy)₂]^- anions also correctly predict a weak exchange J₁(calc) = -0.92 K. DFT calculations revealed the OPBE (GGA-type) functional that correctly predicts the spin-crossover transition for the iron(III) thpy compounds. Besides, the effect of the N₂O₄, N₂S₂O₂, and N₂Se₂O₂ coordination environments on the energy stabilization of the LS state of iron(III) anionic thpy complexes was noted as well.

Computational demonstration of isomer- and spin-state-dependent charge transport in molecular junctions composed of charge-neutral iron(II) spin-crossover complexes

Chemistry offers a multitude of opportunities towards harnessing functional molecular materials with application propensity. One emerging area of interest is molecular spintronics, in which charge and spin degrees of freedom have been used to achieve power-efficient device architectures. Herein, we show that, with the aid of state-of-the-art quantum chemical calculations on designer molecular junctions, the conductance and spin filtering capabilities are molecular structure-dependent. As inferred from the calculations, structural control over the transport can be achieved by changing the position of the thiomethyl (SMe) anchoring groups for Au(111) electrodes in a set of isomeric 2,2'-bipyridine-based metal coordinating ligand entities L1 and L2. The computational studies on heteroleptic iron(II) coordination complexes (1 and 2) composed of L1 and L2 reveal that switching the spin-state of the iron(II) centers, from the low-spin (LS) to high-spin (HS) state, by means of an external electric field stimulus, could, in theory, be performed. Such switching, known as spin-crossover (SCO), renders charge transport through single-molecule junctions of 1 and 2 spin-state-dependent, and the HS junctions are more conductive than the LS junctions for both complexes. Additionally, the LS and HS junctions based on complex 1 are more conductive than those featuring complex 2. Moreover, it is predicted that the spin filtering efficiency (SFE) of the HS junctions strongly depends on the bridging complex geometry, with 1 showing a voltage-dependent SFE, whereas 2 exhibits an SFE of practically 100% over all the studied voltage range. To be pragmatic towards applications, the ligands L1 and L2 and complex 1 have been successfully synthesized, and the spin-state switching propensity of 1 in the bulk state has been elucidated. The results shown in this study might lead to the synthesis and characterization of isomeric SCO complexes with tuneable spin-state switching and charge transport properties.

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.

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.

Push and Pull Effect of Methoxy and Nitro Groups Modifies the Spin-State Switching Temperature in Fe(III) Complexes

We have explored the impact of electron-donating (methoxy) and electron-withdrawing (nitro) substituents on SalEen ligand based spin crossover (SCO) behavior of Fe(III) complexes. Thus, 3-X-substituted SalEen ligands were employed to prepare [Fe(3-X-SalEen)₂]·NCSe, where X = OMe (1), H (2), and NO₂ (3) (3-X-SalEen is the condensation product of 3-substituted salicylaldehyde and N-ethylethylenediamine). The characteristic spin transition temperature (T _1/2) is shown to shift to a lower temperature when an electron-donating substituent (OMe) is used and to a higher temperature when an electron-withdrawing substituent (NO₂) is used. We used experimental and theoretical methods to determine the reasons for this behavior. The solid-state magnetic data revealed the transition temperatures for complexes 1, 2, and 3 to be 219, 251, and 366 K, respectively. The solution-state magnetic data also support this trend in T _1/2 values. UV-vis spectra analysis indicates that there is greater delocalization in the π-manifold of the ligand when the nitro group is the substituent. Theoretical studies through density functional theory methods suggest the methoxy substituent decreases the energy gap between the t_2g and e_g orbitals (explaining the lower T _1/2 value), while the nitro substituent increases the energy gap between the t_2g and e_g orbitals and thus increases the T _1/2 value.

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.

Synthesis and characterisation of group 8 tris(1-benzyl-1,2,3-triazol-4-yl)-p-anisolylmethane complexes

The tris(1,2,3-triazol-4-yl)methane framework offers a highly versatile architecture for ligand design, yet the coordination chemistry of this class of ligand remains largely unexplored. We report here the synthesis and characterisation of the homoleptic complexes [M(ttzm)₂](PF₆)₂ (ttzm = tris(1-benzyl-1,2,3-triazol-4-yl)-p-anisolylmethane; M = Fe (Fe), Ru (Ru), Os (Os)). Initial attempts to prepare Ru by reaction of [Ru(p-cymene)Cl₂]₂ and ttzm also led to the isolation of the heteroleptic complex [Ru(p-cymene)(ttzm)](PF₆)₂. The structures of [Ru(p-cymene)(ttzm)](PF₆)₂, [Fe(ttzm)₂]²⁺ (as its BPh₄^- salt) and Os were solved by X-ray diffraction. The homoleptic Fe(II) and Os(II) containing cations adopt distorted octahedral geometries due to the steric interactions between the ansiole and triazole rings of the ttzm ligands. The homoleptic complexes all adopt a low-spin d⁶ configuration and exhibit reversible M(II)/M(III) processes (+0.35 to +0.72 V vs. Fc/Fc⁺). These oxidation processes are cathodically shifted relative to those of related hexatriazole donor based complexes with density functional theory (DFT) calculations showing the metal d-orbitals are destabilised through a π-donor contribution from the triazole rings. The complexes all show prominent UV-visible absorption bands between 350 and 450 nm assigned to transitions of ¹MLCT character. Whilst none of the homoleptic complexes are emissive in room temperature fluid solutions, Os is emissive at 77 K in an EtOH/MeOH glass (λ_max 472 nm).

Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks

Fueled by advances in hardware and algorithm design, large-scale automated explorations of chemical reaction space have become possible. Here, we present our approach to an open-source, extensible framework for explorations of chemical reaction mechanisms based on the first-principles of quantum mechanics. It is intended to facilitate reaction network explorations for diverse chemical problems with a wide range of goals such as mechanism elucidation, reaction path optimization, retrosynthetic path validation, reagent design, and microkinetic modeling. The stringent first-principles basis of all algorithms in our framework is key for the general applicability that avoids any restrictions to specific chemical systems. Such an agile framework requires multiple specialized software components of which we present three modules in this work. The key module, Chemoton, drives the exploration of reaction networks. For the exploration itself, we introduce two new algorithms for elementary-step searches that are based on Newton trajectories. The performance of these algorithms is assessed for a variety of reactions characterized by a broad chemical diversity in terms of bonding patterns and chemical elements. Chemoton successfully recovers the vast majority of these. We provide the resulting data, including large numbers of reactions that were not included in our reference set, to be used as a starting point for further explorations and for future reference.

The Apparently Unreactive Substrate Facilitates the Electron Transfer for Dioxygen Activation in Rieske Dioxygenases

Rieske dioxygenases belong to the non-heme iron family of oxygenases and catalyze important cis-dihydroxylation as well as O-/N-dealkylation and oxidative cyclization reactions for a wide range of substrates. The lack of substrate coordination at the non-heme ferrous iron center, however, makes it particularly challenging to delineate the role of the substrate for productive O 2 activation. Here, we studied the role of the substrate in the key elementary reaction leading to O 2 activation from a theoretical perspective by systematically considering (i) the 6-coordinate to 5-coordinate conversion of the non-heme FeII upon abstraction of a water ligand, (ii) binding of O 2 , and (iii) transfer of an electron from the Rieske cluster. We systematically evaluated the spin-state-dependent reaction energies and structural effects at the active site for all combinations of the three elementary processes in the presence and absence of substrate using naphthalene dioxygenase as a prototypical Rieske dioxygenase. We find that reaction energies for the generation of a coordination vacancy at the non-heme FeII center through thermoneutral H2 O reorientation and exothermic O 2 binding prior to Rieske cluster oxidation are largely insensitive to the presence of naphthalene and do not lead to formation of any of the known reactive Fe-oxygen species. By contrast, the role of the substrate becomes evident after Rieske cluster oxidation and exclusively for the 6-coordinate non-heme FeII sites in that the additional electron is found at the substrate instead of at the iron and oxygen atoms. Our results imply an allosteric control of the substrate on Rieske dioxygenase reactivity to happen prior to changes at the non-heme FeII in agreement with a strategy that avoids unproductive O 2 activation.

Toward Simulation of Fe(II) Low-Spin → High-Spin Photoswitching by Synergistic Spin-Vibronic Dynamics

A new theoretical approach is presented and applied for the simulation of Fe(II) low-spin (LS, singlet, t_2g⁶e_g⁰) → high-spin (HS, quintet, t_2g⁴e_g²) photoswitching dynamics of the octahedral model complex [Fe(NCH)₆]²⁺. The utilized synergistic methodology heavily exploits the strengths of complementary electronic structure and spin-vibronic dynamics methods. Specifically, we perform 3D quantum dynamics (QD) and full-dimensional trajectory surface hopping (TSH, in conjunction with a linear vibronic coupling model), with the modes for QD selected by TSH. We follow a hybrid approach which is based on the application of time-dependent density functional theory (TD-DFT) excited-state potential energy surfaces (PESs) and multiconfigurational second-order perturbation theory (CASPT2) spin–orbit couplings (SOCs). Our method delivers accurate singlet–triplet–quintet intersystem crossing (ISC) dynamics, as assessed by comparison to our recent high-level ab initio simulations and related time-resolved experimental data. Furthermore, we investigate the capability of our simulations to identify the location of ISCs. Finally, we assess the approximation of constant SOCs (calculated at the Franck–Condon geometry), whose validity has central importance for the combination of TD-DFT PESs and CASPT2 SOCs. This efficient methodology will have a key role in simulating LS → HS dynamics for more complicated cases, involving higher density of states and varying electronic character, as well as the analysis of ultrafast experiments.

CASPT2 molecular geometries of Fe(II) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.

Decoding Regioselective Reaction Mechanism of the Gentisic Acid Catalyzed by Gentisate 1,2-Dioxygenase Enzyme

Gentisate 1,2-dioxygenase (GDO), a ring-fission non-heme dioxygenase enzyme, displays a unique regioselective reaction of gentisic acid (GTQ) in the presence of molecular oxygen. GTQ is an important intermediate in the...

Accurate calculation of spin-state energy gaps in Fe(III) spin-crossover systems using density functional methods

Fe(III) complexes are receiving ever-increasing attention as spin crossover (SCO) systems because they are usually air stable, as opposed to Fe(II) complexes, which are prone to oxidation. Here, we present the first systematic study exclusively devoted to assess the accuracy of several exchange-correlation functionals when it comes to predicting the energy gap between the high-spin (S = 5/2) and the low-spin (S = 1/2) states of Fe(III) complexes. Using a dataset of 24 different Fe(III) hexacoordinated complexes, it is demonstrated that the B3LYP* functional is an excellent choice not only for predicting spin-state energy gaps for Fe(III) complexes undergoing spin-transitions but also for discriminating Fe(III) complexes that are either low- or high-spin in the whole range of temperatures. Our benchmark study has led to the identification of a very versatile Fe(III) compound whose SCO properties can be engineered upon changing a single axial ligand. Overall, this work demonstrates that B3LYP* is a reliable functional for screening new spin-crossover systems with tailored properties.

Peculiar Spin-Crossover Behavior in the 2D Polymer K[FeIII(5Cl-thsa)2]

A potassium salt of the N₂S₂O₂-coordination Fe(III) anion K[Fe(5Cl-thsa)₂] (1) (5Cl-thsa - 5-chlorosalicylaldehyde thiosemicarbazone) is synthesized and characterized structurally and magnetically over a wide temperature range. Two polymorphs of salt 1 characterized by the common 2D polymer nature and assigned to the same orthorhombic Pbcn space group have been identified. The molecular structure of the minor polymorph of 1 was solved and refined at 100, 250, and 300 K is shown to correspond to the LS configuration. The dominant polymorph of 1 features K⁺ cations disordered over a few crystallographic sites, while the minor polymorph includes fully ordered K⁺ cations. The major polymorph exhibits a complete three-step cooperative spin-crossover transition both in the heating and cooling modes: The first step occurs in a temperature range from 2 to 50 K; the second abrupt hysteretic step occurs from 200 to 250 K with T_1/2 = 230 K and a 6 K hysteresis loop. The third gradual step occurs from 250 to 440 K. According to ⁵⁷Fe Mössbauer, XRPD, and EXAFS data, the spin-crossover transition for the dominant polymorph is quite peculiar. Indeed, the increase in the HS concentration by 57% at the second step does not result in the expected significant increase in the iron(III)-ligand bond lengths. In addition, the final step of the spin conversion (Δγ_HS = 26%) is associated with a structural phase transition with a symmetry lowering from the orthorhombic (Pbcn) to the monoclinic (P2₁/n) space group. This nontrivial phenomenon was investigated in detail by applying magnetization measurements, electron spin resonance, ⁵⁷Fe Mössbauer spectroscopy, and DFT calculations. These results provide a new platform for understanding the multistep spin-crossover character in the Fe(III) thsa-complexes and related compounds.

Halogenation affects driving forces, reorganization energies and "rocking" motions in strained [Fe(tpy)2]2+ complexes

Controlling the energetics of spin crossover (SCO) in Fe(II)-polypyridine complexes is critical for designing new multifunctional materials or tuning the excited-state lifetimes of iron-based photosensitizers. It is well established that the Fe-N "breathing" mode is important for intersystem crossing from the singlet to the quintet state, but this does not preclude other, less obvious, structural distortions from affecting SCO. Previous work has shown that halogenation at the 6 and 6'' positions of tpy (tpy = 2,2';6',2''-terpyridine) in [Fe(tpy)₂]²⁺ dramatically increased the lifetime of the excited MLCT state and also had a large impact on the ground state spin-state energetics. To gain insight into the origins of these effects, we used density functional theory calculations to explore how halogenation impacts spin-state energetics and molecular structure in this system. Based on previous work we focused on the ligand "rocking" motion associated with SCO in [Fe(tpy)₂]²⁺ by constructing one-dimensional potential energy surfaces (PESs) along the tpy rocking angle for various spin states. It was found that halogenation has a clear and predictable impact on ligand rocking and spin-state energetics. The rocking is correlated to numerous other geometrical distortions, all of which likely affect the reorganization energies for spin-state changes. We have quantified trends in reorganization energy and also driving force for various spin-state changes and used them to interpret the experimentally measured excited-state lifetimes.

Photoinduced Low-Spin → High-Spin Mechanism of an Octahedral Fe(II) Complex Revealed by Synergistic Spin-Vibronic Dynamics

The Fe(II) low-spin (LS; ¹A_1g, t_2g⁶e_g⁰) → high-spin (HS; ⁵T_2g, t_2g⁴e_g²) light-induced excited spin state trapping (LIESST) mechanism solely involving metal-centered states is revealed by synergistic spin-vibronic dynamics simulations. For the octahedral [Fe(NCH)₆]²⁺ complex, we identify an initial ∼100 fs ¹T_1g → ³T_2g intersystem crossing, controlled by vibronic coupling to antisymmetric Fe-N stretching motion. Subsequently, population branching into ³T_1g, ⁵T_2g (HS), and ¹A_1g (LS) is observed on a subpicosecond time scale, with the dynamics dominated by coherent Fe-N breathing wavepackets. These findings are consistent with ultrafast experiments, methodologically establish a new state of the art, and will give a strong impetus for further intriguing dynamical studies on LS → HS photoswitching.

Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning

Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.

H2 activation by hydrogenase-inspired NiFe catalyst using frustrated Lewis pair: effect of buffer and halide ion in the heterolytic H-H bond cleavage

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production. We investigate the reaction mechanism of H₂ activation in an aqueous solution by the recently developed NiFe complex (Ogo et al. Sci. Adv. 2020, 6, eaaz8181) using density functional theory (DFT) calculation. Our computational results showed that H₂ is activated using frustrated Lewis pair. That is, H₂ binds to the Fe site of the NiFe complex, acting as a Lewis acid, while the added buffer, acting as Lewis base, abstracts protons to form a hydride complex. Furthermore, the higher basicity in the proton abstraction reaction characterises reaction more exergonic and lowers the reaction barrier. In addition, in the proton abstraction by the water molecule, the reaction barrier was lowered when anion such as Cl⁻ is in the vicinity of the water. Understanding the chemical species that contribute to the catalytic process in cooperation with the metal catalyst at the atomic level should help to maximise the function of the catalyst.

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production.

Role of Triplet States in the Photodynamics of Aniline

The dynamics of excited heteroaromatic molecules is a key to understanding the photoprotective properties of many biologically relevant chromophores that dissipate their excitation energy nonreactively and thereby prevent the detrimental effects of ultraviolet radiation. Despite their structural variability, most substituted aromatic compounds share a common feature of a repulsive ¹πσ* potential energy surface. This surface can lead to photoproducts, and it can also facilitate the population transfer back to the ground electronic state by means of a ¹πσ*/S₀ conical intersection. Here, we explore a hidden relaxation route involving the triplet electronic state of aniline, which has recently been discovered by means of time-selected photofragment translational spectroscopy [J. Chem. Phys. 2019, 151, 141101]. By using the recently available analytical gradients for multiconfiguration pair-density functional theory, it is now possible to locate the minimum-energy crossing points between states of different spin and therefore compute the intersystem crossing rates with a multireference method, rather than with the less reliable single-reference methods. Using such calculations, we demonstrate that the population loss of aniline in the T₁(³ππ*) state is dominated by C₆H₅NH₂ → C₆H₅NH· + H· dissociation, and we explain the long nonradiative lifetimes of the T₁(³ππ*) state at the excitation wavelengths of 294-264 nm.

First-Principles Investigations of Magnetic Anisotropy and Spin-Crossover Behavior of Fe(III)-TBP Complexes

With the ongoing effort to obtain mononuclear 3d-transition-metal complexes that manifest slow relaxation of magnetization and, hence, can behave as single-molecule magnets (SMMs), we have modeled 14 Fe(III) complexes based on an experimentally synthesized (PMe₃)₂FeCl₃ complex [J. Am. Chem. Soc. 2017, 139 (46), 16474-16477], by varying the axial ligands with group XV elements (N, P, and As) and equatorial halide ligands from F, Cl, Br, and I. Out of these, nine complexes possess large zero field splitting (ZFS) parameter D in the range of -40 to -60 cm^-1. The first-principles investigation of the ground-spin state applying density functional theory (DFT) and wave function-based multiconfigurations methods, e.g., SA-CASSCF/NEVPT2, are found to be quite consistent except for few delicate cases with near-degenerate spin states. In such cases, the hybrid B3LYP functional is found to be biased toward high-spin (HS) state. Altering the percentage of exact exchange admixed in the B3LYP functional leads to intermediate-spin (IS) ground state consistent with the multireference calculations. The origin of large zero field splitting (ZFS) in the Fe(III)-based trigonal bipyramidal (TBP) complexes is investigated. Furthermore, a number of complexes are identified with very small ΔG_HS-IS^adia. values indicating the possible spin-crossover phenomenon between the bistable spin states.

Theoretical inspection of the spin-crossover [Fe(tzpy)2(NCS)2] complex on Au(100) surface

We explore the deposition of the spin-crossover [Fe(tzpy)₂(NCS)₂] complex on the Au(100) surface by means of density functional theory (DFT) based calculations. Two different routes have been employed: low-cost finite cluster-based calculations, where both the Fe complex and the surface are maintained fixed while the molecule approaches the surface; and periodic DFT plane-wave calculations, where the surface is represented by a four-layer slab and both the molecule and surface are relaxed. Our results show that the bridge adsorption site is preferred over the on-top and fourfold hollow ones for both spin states, although they are energetically close. The LS molecule is stabilized by the surface, and the HS-LS energy difference is enhanced by about 15%-25% once deposited. The different Fe ligand field for LS and HS molecules manifests on the composition and energy of the low-lying bands. Our simulated STM images indicate that it is possible to distinguish the spin state of the deposited molecules by tuning the bias voltage of the STM tip. Finally, it should be noted that the use of a reduced size cluster to simulate the Au(100) surface proves to be a low-cost and reliable strategy, providing results in good agreement with those resulting from state-of-the-art periodic calculations for this system.

On the Spin-State Dependence of Redox Potentials of Spin Crossover Complexes

Resistance switching properties of nanoscale junctions of spin crossover molecules have received recently much interest. In many cases, this property has been traced back to the variation of molecular orbital energies upon spin transition. However, one can also expect a substantial reorganization of the molecular structure due to charge localization, which calls for a better understanding of the relationship between the redox potential and the spin state of the molecule. To investigate this issue, we carried out a detailed density functional theory and variable temperature cyclic voltammetry investigation of the benchmark compound [Fe(HB(1,2,4-triazol-1-yl)₃)₂] in solution. We show that, for a correct thermodynamical picture, it is necessary to take into account the charge transfer-induced electronic and structural reorganization as well as spin equilibria in the oxidized and reduced species.

Quantum-chemistry-aided ligand engineering for potential molecular switches: changing barriers to tune excited state lifetimes

Substitution of terpyridine at the 4' position with electron withdrawing and donating groups is used to tune the quintet lifetime of its iron(ii) complex. DFT calculations suggest that the energy barrier between the quintet and singlet states can be altered significantly upon substitution, inducing a large variation of the lifetime of the photoexcited quintet state. This prediction was experimentally verified by transient optical absorption spectroscopy and good agreement with the trend expected from the calculations was found. This demonstrates that the potential energy landscape can indeed be rationally tailored by relevant modifications based on DFT predictions. This result should pave the way to advancing efficient theory-based ligand engineering of functional molecules to a wide range of applications.

Mechanical force-induced manipulation of electronic conductance in a spin-crossover complex: a simple approach to molecular electronics

The atomic-scale technological sophistication from the last half-decade provides new avenues for the atom-by-atom fabrication of nanostructures with extraordinary precision. This urges the appraisal of the fabrication scheme layout for a modular nanoelectronic device based on an individual molecular complex. The mechanical force-induced distortion to the metal coordination sphere triggers a low-spin (LS) to high-spin (HS) electronic transition in the complex. The controlled structural distortions (relative to a specific bond-angle) are deemed to be the switching parameter for the observed spin-transitions. Mechanical stretching is the key to engineering a spin-state switch in the proposed molecular device. The spin-dependent reversible variation in the electronic conductance concurrent to the unique spin-states can be understood from the state-of-the-art Nonequilibrium Green's Function (NEGF) calculations. Combined with NEGF calculations, the DFT study further provides a qualitative perception of the electronic conductance in the two-terminal device architecture. From the transport calculations, there is also evidence of considerable fluctuation in the spin-dependent electronic conductance at the molecular junction with relative variations in the scattering limit. Subsequently, the present study shows significant advances in the transmission probabilities for the high-spin state of the Fe(ii) complex. The results empower the progress of nanoelectronics at the single molecule level.

Modulation of Jahn-Teller distortion and electromechanical response in a Mn3+spin crossover complex

Structural, magnetic and electromechanical changes resulting from spin crossover between the spin quintet and spin triplet forms of a mononuclear Mn3+complex embedded in six lattices with different charge balancing counterions are reported. Isostructural ClO4-and BF4-salts (1) and (2) each have two unique Mn3+sites which follow different thermal evolution pathways resulting in a crossover from the spin quintet form at room temperature to a 1:1 spin triplet:quintet ratio below 150 K. The PF6-(3) and NO3-(4) salts which each have one unique Mn3+site show a complete conversion from spin quintet to spin triplet over the same temperature range. A complete two step spin crossover is observed in the CF3SO3-lattice (5) with a 1:1 ratio of spin quintet and spin triplet forms at intermediate temperature, while the BPh4-lattice (6) stabilizes the spin triplet form over most of the temperature range with gradual and incomplete spin state switching above 250 K. An electromechanical piezoresponse was detected in NO3-complex4despite crystallization in a centrosymmetric space group. The role of deformations associated with stress-induced spin triplet-spin quintet switching in breaking the local symmetry are discussed and computational analysis is used to estimate the energy gap between the two spin states.

Using electronegativity and hardness to test density functionals

Density functional theory (DFT) is used in thousands of papers each year, yet lack of universality reduces DFT's predictive capacity, and functionals may produce energy-density imbalances. The absolute electronegativity (χ) and hardness (η) directly reflect the energy-density relationship via the chemical potential ∂E/∂N and we thus hypothesized that they probe universality. We studied χ and η for atoms Z = 1-36 using 50 diverse functionals covering all major classes. Very few functionals describe both χ and η well. η benefits from error cancellation, whereas χ is marred by error propagation from IP and EA; thus, almost all standard GGA and hybrid functionals display a plateau in the MAE at ∼0.2 eV-0.3 eV for η. In contrast, variable performance for χ indicates problems in describing the chemical potential by DFT. The accuracy and precision of a functional is far from linearly related, yet for a universal functional, we expect linearity. Popular functionals such as B3LYP, PBE, and revPBE perform poorly for both properties. Density sensitivity calculations indicate large density-derived errors as occupation of degenerate p- and d-orbitals causes "non-universality" and large dependency on exact exchange. Thus, we argue that performance for χ for the same systems is a hallmark of an important aspect of universality by probing ∂E/∂N. With this metric, B98, B97-1, PW6B95D3, MN-15, rev-TPSS, HSE06, and APFD are the most "universal" among the tested functionals. B98 and B97-1 are accurate for very diverse metal-ligand bonds, supporting that a balanced description of ∂E/∂N and ∂E²/∂N², via χ and η, is probably a first simple probe of universality.

Macroscopic Polarization Change via Electron Transfer in a Valence Tautomeric Cobalt Complex

Polarization change induced by directional electron transfer attracts considerable attention owing to its fast switching rate and potential light control. Here, we investigate electronic pyroelectricity in the crystal of a mononuclear complex, [Co(phendiox)(rac-cth)](ClO₄)·0.5EtOH (1·0.5EtOH, H₂phendiox = 9, 10-dihydroxyphenanthrene, rac-cth = racemic 5, 5, 7, 12, 12, 14-hexamethyl-1, 4, 8, 11-tetraazacyclotetradecane), which undergoes a two-step valence tautomerism (VT). Correspondingly, pyroelectric current exhibits double peaks in the same temperature domain with the polarization change consistent with the change in dipole moments during the VT process. Time-resolved Infrared (IR) spectroscopy shows that the photo-induced metastable state can be generated within 150 ps at 190 K. Such state can be trapped for tens of minutes at 7 K, showing that photo-induced polarization change can be realized in this system. These results directly demonstrate that a change in the molecular dipole moments induced by intramolecular electron transfer can introduce a macroscopic polarization change in VT compounds.

Polarization change from directional electron transfer attracts considerable attention owing to its fast switching rate and potential light control. Here, the authors provide a proof-of-concept of electronic pyroelectricity induced by intramolecular electron transfer in the single crystal of a valence tautomeric compound.

Seeing Is Believing: Experimental Spin States from Machine Learning Model Structure Predictions

Determination of ground-state spins of open-shell transition-metal complexes is critical to understanding catalytic and materials properties but also challenging with approximate electronic structure methods. As an alternative approach, we demonstrate how structure alone can be used to guide assignment of ground-state spin from experimentally determined crystal structures of transition-metal complexes. We first identify the limits of distance-based heuristics from distributions of metal–ligand bond lengths of over 2000 unique mononuclear Fe(II)/Fe(III) transition-metal complexes. To overcome these limits, we employ artificial neural networks (ANNs) to predict spin-state-dependent metal–ligand bond lengths and classify experimental ground-state spins based on agreement of experimental structures with the ANN predictions. Although the ANN is trained on hybrid density functional theory data, we exploit the method-insensitivity of geometric properties to enable assignment of ground states for the majority (ca. 80–90%) of structures. We demonstrate the utility of the ANN by data-mining the literature for spin-crossover (SCO) complexes, which have experimentally observed temperature-dependent geometric structure changes, by correctly assigning almost all (>95%) spin states in the 46 Fe(II) SCO complex set. This approach represents a promising complement to more conventional energy-based spin-state assignment from electronic structure theory at the low cost of a machine learning model.