Computational approach to the study of thermal spin crossover phenomena

Improving the Accuracy in the Prediction of Transition-Metal Spin-State Energetics Using a Robust Variation-Based Approach: Density Functional Theory, CASPT2 and MC-PDFT Applied to the Case Study of Tris-Diimine Fe(II) Complexes

Designing ligands for transition metal complexes with a specified low-spin, high-spin or spin-crossover behavior is challenging. A major advance was recently made by Phan et al. [J. Am. Chem. Soc. 2017, 139, 6437-6447] who showed that the spin state of a homoleptic tris-diimine Fe(II) complex can be predicted from the N-N distance in the free diimine. They could thus predict the change in magnetic behavior on passing from the complexes of 2,2'-bipyridine (bpy), 2,2'-biimidazole (bim) and 2,2'-bis-2-imidazoline (bimz) ligands to those obtained with the modified analogs 4,5-diazafluoren-9-one (dafo), 1,1'-(α,α'-o-xylyl)-2,2'-bisimidazole (xbim) and 2,3,5,6,8,9-hexahydrodiimidazo[1,2-a:2', 1'-c]pyrazine (etbimz), respectively. Theoretically, the challenge lies in the accurate determination of the HS-LS zero-point energy difference ΔEHL°. The issue can be circumvented by using a variation-based approach, wherein ΔEHL° is not directly evaluated but obtained from the estimate of its variation Δ(ΔEHL°) in series of related systems, which include one whose ΔEHL° is accurately known [Phys. Chem. Chem. Phys. 2013, 15, 3752-3763; J. Phys. Chem. A 2022, 126, 6221-6235]. In this study, density functional theory (DFT), second-order multireference perturbation theory in its CASPT2 formulation, multiconfigurational pair DFT (MC-PDFT) and its hybrid formulation (HMC-PDFT) have been applied to the determination of Δ(ΔEHL°) in the pairs of complexes ( [ F e ( b p y ) 3 ] 2 + , [ F e ( d a f o ) 3 ] 2 + ) , ( [ F e ( b i m ) 3 ] 2 + , [ F e ( x b i m ) 3 ] 2 + ) and ( [ F e ( b i m z ) 3 ] 2 + , [ F e ( e t b i m z ) 3 ] 2 + ) . In DFT, we used several semilocal functionals and their global hybrids, as well as their D2, D3, D3BJ and D4 dispersion-corrected forms; and in MC-PDFT, different translated and fully translated functionals. The results are consistent with one another and in very good agreement with experiments. They show small to vanishing dependence on key details of the methods used: namely, the exact-exchange contribution to global hybrids; the ionization potential-electron affinity shift and basis sets used in the CASPT2 calculations; and the active spaces employed for the CASSCF wave functions used in the MC-PDFT and HMC-PDFT calculations. Insights into the change in the spin-state energetics accompanying the ligand exchanges were gained through a complexation energy analysis. Using the accurate CCSD(T) estimate of the HS-LS adiabatic energy difference in [ F e ( N C H ) 6 ] 2 + [J. Chem. Theory Comput. 2012, 8, 4216-4231], the Δ(ΔEHL°)-approach has been applied to the determination of ΔEHL° in the diimine complexes. The CASPT2 and DFT-D2 methods only give results in agreement with experiments. This suggests for the other methods a limitation in their treatment of dispersion which prevents them from accurately describing the spin-state energetics change accompanying the passing from [ F e ( N C H ) 6 ] 2 + with the tetragonal arrangement of its nitrile ligands to the tris-diimine complexes with the trigonal packing of their bulkier ligands.

Investigating the influence of oriented external electric fields on modulating spin-transition temperatures in Fe(II) SCO complexes: a theoretical perspective

Spin-crossover complexes, valued for their bistability, are extensively studied due to their numerous potential applications. A primary challenge in this molecular class is to identify effective methods to adjust the spin-transition temperature, which frequently falls outside the desired temperature range. This typically necessitates intricate chemical design and synthesis or the use of stimuli such as light or pressure, each introducing its own set of challenges for integrating these molecules into end-user applications. In this work, we aim to address this challenge using an oriented external electric field (OEEF) as one stimulus to modulate the spin-transition temperatures. For this purpose, we have employed both periodic and non-periodic calculations on three well-characterised Fe(II) SCO complexes, namely [Fe(phen)2(NCS)2] (1, phen = 1,10-phenanthroline), [Fe(bt)2(NCS)2] (2, bt = 2,2'-bi-2-thiazoline) and [Fe(py)2phen(NCS)2] (3, py = pyridine) possessing a similar structural motif of {FeN4N'2}. To begin with, DFT calculations employing the TPSSh functional were performed on complexes 1 to 3, and the estimated low-spin (LS) and high-spin (HS) gaps are 24.6, 15.3 and 15.4 kJ mol-1, and these are in the range expected for Fe(II) SCO complexes. In the next step, an OEEF was applied in the molecule along the pseudo-C2 axis that bisects two coordinated -NCS groups. Application of an OEEF was found to increase the Fe-ligand bond length and found to affect the spin-transition at the particular applied OEEF. While the HS state of 1 becomes the ground state at an applied field of 0.514 V Å-1, the LS state lies at a higher energy of 1.3 kJ mol-1. Similarly, complexes 2 and 3 also show the HS ground state at an applied field of 0.514 V Å-1, where the LS state stays at higher energies of 6.13 and 11.62 kJ mol-1, respectively. It is found that the overall change in enthalpy (ΔHHL) and entropy (ΔSHL) for the spin transition in the presence of OEEFs decreases upon increasing the strength of the applied field. The computed spin-transition temperature (T1/2) using DFT was found to be in close agreement with the experimentally reported values. It is estimated that on increasing the strength of the applied electric field, the T1/2 increases significantly. While the DFT computed T1/2 values for the optimised geometry of 1, 2 and 3 were found to be 134.6 K, 159.9 K and 111.4 K respectively, at the applied field of 0.6425 V Å-1T1/2 increases up to 187.3 K, 211.0 K and 184.4 K respectively, unveiling an hitherto unknown strategy to tune the T1/2 values. A limited benchmarking was performed with five additional exchange-correlation functionals: PBE, BLYP, B3LYP*, B3LYP, and PBE0. These functionals were found to be unsuitable for predicting the correct SCO behaviour for complex 2, and their behaviour under various electric fields did not improve. This emphasises the importance of choosing the correct functional at zero OEEF prior to testing them under various electric fields. Furthermore, calculations were performed with complex 1 adsorbed on the Au(111) surface. The formation of an Au-S bond during adsorption significantly stabilises the low-spin (LS) state, hindering the observation of spin-crossover (SCO) behaviour. Nonetheless, the application of an OEEF reduces this gap and brings the T1/2 value closer to the desired temperature. This offers a novel post-fabrication strategy for attaining SCO properties at the interface.

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.

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.

Transition Temperature for Spin-Crossover Materials with the Mean Value Ensemble Hubbard-U Correction

Calculation of transition temperatures T1/2 for thermally driven spin-crossover in condensed phases is challenging, even with sophisticated state-of-the-art density functional approximations. The first issue is the accuracy of the adiabatic crossover energy difference ΔEHL between the low- and high-spin states of the bistable metal-organic complexes. The other is the proper inclusion of entropic contributions to the Gibbs free energy from the electronic and vibrational degrees of freedom. We discuss the effects of treatments of both contributions upon the calculation of thermochemical properties for a set of 20 spin-crossover materials using a Hubbard-U correction obtained from a reference ensemble spin-state. The U values obtained from a simplest bimolecular representation may overcorrect, somewhat, the ΔEHL values, hence giving somewhat excessive reduction of the T1/2 results with respect to their U = 0 values in the crystalline phase. We discuss the origins of the discrepancies by analyzing different sources of uncertainties. By use of a first-coordination-sphere approximation and the assumption that vibrational contributions from the outermost atoms in a metal-organic complex are similar in both low- and high-spin states, we achieve T1/2 results with the low-cost, widely used PBE generalized gradient density functional approximation comparable to those from the more costly, more sophisticated r2SCAN meta-generalized gradient approximation. The procedure is promising for use in high-throughput materials screening, because it combines rather low computational effort requirements with freedom from user manipulation of parameters.

Surface-Induced Electronic and Vibrational Level Shifting of [Fe(py)2bpym(NCS)2] on Al(100)

It is essential that one understands how the surface degrees of freedom influence molecular spin switching to successfully integrate spin crossover (SCO) molecules into devices. This study uses density functional theory calculations to investigate how spin state energetics and molecular vibrations change in a Fe(II) SCO compound named [Fe(py)2bpym(NCS)2] when deposited on an Al(100) surface. The calculations consider an environment-dependent U to assess the local Coulomb correlation of 3d electrons. The results show that the adsorption configurations heavily affect the spin state splitting, which increases by 10-40 kJmol-1 on the surface, and this is detrimental to spin conversion. This effect is due to the surface binding energy variation across the spin transition. The preference for the low-spin state originates partly from the strong correlation effect. Furthermore, the surface environment constrains the vibrational entropy difference, which decreases by 8-17 Jmol-1K-1 (at 300 K) and leads to higher critical temperatures. These results suggest that the electronic energy splitting and vibrational level shifting are suitable features for characterizing the spin transition process on surfaces, and they can provide access to high-throughput screening of spin crossover devices.

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).

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.

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

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

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.

Benchmarking CASPT3 Vertical Excitation Energies

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

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

Panchromatic Absorption and Oxidation of an Iron(II) Spin Crossover Complex

In order to expand and exploit the useful properties of d⁶-iron(II) and d⁵-iron(III) complexes in potential magnetic, photophysical, or magnetooptical applications, crucial ligand-controlled parameters are the ligand field strength in a given coordination mode and the availability of suitable metal and ligand frontier orbitals for charge-transfer processes. The push-pull ligand 2,6-diguanidylpyridine (dgpy) features low-energy π* orbitals at the pyridine site and strongly electron-donating guanidinyl donors combined with the ability to form six-membered chelate rings for optimal metal-ligand orbital overlap. The electronic ground states of the pseudo-octahedral d⁶- and d⁵-complexes mer-[Fe(dgpy)₂]²⁺, cis-fac-[Fe(dgpy)₂]²⁺, and mer-[Fe(dgpy)₂]³⁺ as well as their charge-transfer (CT) and metal-centered (MC) excited states are probed by variable temperature UV/vis absorption, NMR, EPR, and Mössbauer spectroscopy, magnetic susceptibility measurements at variable temperature as well as quantum chemical calculations.

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.

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.

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.

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.

Quantitative and Chemically Intuitive Evaluation of the Nature of M-L Bonds in Paramagnetic Compounds: Application of EDA-NOCV Theory to Spin Crossover Complexes

To improve understanding of M-L bonds in 3d transition metal complexes, analysis by energy decomposition analysis and natural orbital for chemical valence model (EDA-NOCV) is desirable as it provides a full, quantitative and chemically intuitive ab initio description of the M-L interactions. In this study, a generally applicable fragmentation and computational protocol was established and validated by using octahedral spin crossover (SCO) complexes, as the transition temperature (T1/2 ) is sensitive to subtle changes in M-L bonding. Specifically, EDA-NOCV analysis of Fe-N bonds in five [FeII (Lazine )2 (NCBH3 )2 ], in both low-spin (LS) and paramagnetic high-spin (HS) states led to: 1) development of a general, widely applicable, corrected M+L6 fragmentation, tested against a family of five LS [FeII (Lazine )3 ](BF4 )2 complexes; this confirmed that three Lazine are stronger ligands (ΔEorb,σ+π =-370 kcal mol-1 ) than 2 Lazine +2 NCBH3 (=-335 kcal mol-1 ), as observed. 2) Analysis of Fe-L bonding on LS→HS, reveals more ionic (ΔEelstat ) and less covalent (ΔEorb ) character (ΔEelstat :ΔEorb 55:45 LS→64:36 HS), mostly due to a big drop in σ (ΔEorb,σ ↓50 %; -310→-145 kcal mol-1 ), and a drop in π contributions (ΔEorb,π ↓90 %; -30→-3 kcal mol-1 ). 3) Strong correlation of observed T1/2 and ΔEorb,σ+π , for both LS and HS families (R2 =0.99 LS, R2 =0.95 HS), but no correlation of T1/2 and ΔΔEorb,σ+π (LS-HS) (R2 =0.11). Overall, this study has established and validated an EDA-NOCV protocol for M-L bonding analysis of any diamagnetic or paramagnetic, homoleptic or heteroleptic, octahedral transition metal complex. This new and widely applicable EDA-NOCV protocol holds great promise as a predictive tool.

Spin-phonon coupling and dynamic zero-field splitting contributions to spin conversion processes in iron(II) complexes

Magnetization dynamics of transition metal complexes manifest in properties and phenomena of fundamental and applied interest [e.g., slow magnetic relaxation in single molecule magnets, quantum coherence in quantum bits (qubits), and intersystem crossing (ISC) rates in photophysics]. While spin-phonon coupling is recognized as an important determinant of these dynamics, additional fundamental studies are required to unravel the nature of the coupling and, thus, leverage it in molecular engineering approaches. To this end, we describe here a combined ligand field theory and multireference ab initio model to define spin-phonon coupling terms in S = 2 transition metal complexes and demonstrate how couplings originate from both the static and dynamic properties of ground and excited states. By extending concepts to spin conversion processes, ligand field dynamics manifest in the evolution of the excited state origins of zero-field splitting (ZFS) along specific normal mode potential energy surfaces. Dynamic ZFSs provide a powerful means to independently evaluate contributions from spin-allowed and/or spin-forbidden excited states to spin-phonon coupling terms. Furthermore, ratios between various intramolecular coupling terms for a given mode drive spin conversion processes in transition metal complexes and can be used to analyze the mechanisms of ISC. Variations in geometric structure strongly influence the relative intramolecular linear spin-phonon coupling terms and will define the overall spin state dynamics. While the findings of this study are of general importance for understanding magnetization dynamics, they also link the phenomenon of spin-phonon coupling across fields of single molecule magnetism, quantum materials/qubits, and transition metal photophysics.

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

Light-Induced Spin Crossover in an Intermediate-Spin Penta-Coordinated Iron(III) Complex

(PMe₃)₂FeCl₃ is an Fe(III) complex that exists in the intermediate-spin ground state in a distorted trigonal bipyramidal geometry. An electronic state with high-spin configuration lies in close vicinity to the ground state, making it a potential spin crossover candidate. A mechanistic account of the spin crossover from the lowest quartet state (Q₀) to the lowest sextet state (S₁) of this complex is provided by exploring both thermal and light-induced pathways. The presence of a large barrier between the two spin states suggests a possible thermal spin crossover at a rather high temperature. The light-induced spin crossover is investigated by employing complete active space self-consistent field calculations together with dynamic correlation and spin-orbit coupling for the lowest seven quartet and lowest five sextet states. The system in the Q₀ state upon light absorption is excited to the optically bright Q₄ LMCT state. By following minimum energy pathways along the electronic states, two light-induced pathways for spin crossover are identified. From the Q₄ state, the system can photo-regenerate the ground intermediate-spin state (Q₀) through an internal conversion of Q₄/Q₃ followed by Q₃/S₁ and S₁/Q₀ intersystem crossings. In an alternate route, through Q₄/S₂ intersystem crossing followed by S₂/S₁ internal conversion, the system can complete the spin crossover from the Q₀ to S₁ state.

The Role of Vibrational Anharmonicity in the Computational Study of Thermal Spin Crossover

Spin crossover in transition metal complexes can be studied in great detail with computational chemistry. Over the years, the understanding has grown that the relative stability of high-spin (HS) versus low-spin (LS) states is a subtle balance of many factors that all need to be taken into account for a reliable description. Among the different contributions, the zero-point energy (ZPE) and the entropy play key roles. These quantities are usually calculated assuming a harmonic oscillator model for the molecular vibrations. We investigated the impact of including anharmonic corrections on the ZPE and the entropy and indirectly on the critical temperature of spin crossover. As test systems, we used a set of ten Fe(II) complexes and one Fe(III) complex, covering different coordination modes (mono-, bi-, and tri-dentate ligands), decreasing coordination number upon spin crossover, coordination by second- and third-row atoms, and changes in the oxidation state. The results show that the anharmonicity has a measurable effect, but it is in general rather small, and tendencies are not easily recognized. As a conclusion, we put forward that for high precision results, one should be aware of the anharmonic effects, but as long as computational chemistry is still struggling with other larger factors like the influence of the environment and the accurate determination of the electronic energy difference between HS and LS, the anharmonicity of the vibrational modes is a minor concern.

Observation of the First Spin Crossover in an Iron(II) Complex with an S6 Coordination Environment: Tris[bis(N,N-diethylamino)carbeniumdithiocarboxylato]iron(II) Hexafluorophosphate

For the first time, the spin-crossover (SCO) phenomenon has been observed in an Fe^II -S₆ system in a tris(chelate)-type iron(II) complex with a zwitterionic sulfur donor bidentate, bis(N,N-diethylamino)carbeniumdithiocarboxylate (EtL), [Fe^II (EtL)₃ ](PF₆ )₂ (1), as synthesized by the reaction of a precursor complex [Fe^II (CH₃ CN)₆ ](PF₆ )₂ with EtL. In the solid state, the high-spin (HS) d⁶ state at ambient temperature and the low-spin (LS) d⁶ state at temperatures lower than approximately 240 K were evidenced by magnetic measurements with SQUID and Mössbauer spectra in the temperature range 4-290 K. X-ray analyses of the crystals at various temperatures disclosed that the distorted trigonal prismatic coordination environments essentially do not change; however, contraction of Fe-S distances by approximately 10 % (0.22 Å), ordering of alkyl groups in EtL and PF₆ ^- counteranions, and formation of significant intermolecular S⋅⋅⋅S interactions between adjacent molecules (average distances of 3.59 Å) take place during the transition from the HS to the LS state. A large decrease in the volume of the formula unit (78.1 ų ) might be responsible for the large activation barrier, thereby resulting in a slow phase transition upon cooling.

The IPEA dilemma in CASPT2

Multi-configurational second order perturbation theory (CASPT2) has become a very popular method for describing excited-state properties since its development in 1990. To account for systematic errors found in the calculation of dissociation energies, an empirical correction applied to the zeroth-order Hamiltonian, called the IPEA shift, was introduced in 2004. The errors were attributed to an unbalanced description of open-shell versus closed-shell electronic states and is believed to also lead to an underestimation of excitation energies. Here we show that the use of the IPEA shift is not justified and the IPEA should not be used to calculate excited states, at least for organic chromophores. This conclusion is the result of three extensive analyses. Firstly, we survey the literature for excitation energies of organic molecules that have been calculated with the unmodified CASPT2 method. We find that the excitation energies of 356 reference values are negligibly underestimated by 0.02 eV. This value is an order of magnitude smaller than the expected error based on the calculation of dissociation energies. Secondly, we perform benchmark full configuration interaction calculations on 137 states of 13 di- and triatomic molecules and compare the results with CASPT2. Also in this case, the excited states are underestimated by only 0.05 eV. Finally, we perform CASPT2 calculations with different IPEA shift values on 309 excited states of 28 organic small and medium-sized organic chromophores. We demonstrate that the size of the IPEA correction scales with the amount of dynamical correlation energy (and thus with the size of the system), and gets immoderate already for the molecules considered here, leading to an overestimation of the excitation energies. It is also found that the IPEA correction strongly depends on the size of the basis set. The dependency on both the size of the system and of the basis set, contradicts the idea of a universal IPEA shift which is able to compensate for systematic CASPT2 errors in the calculation of excited states.

Diffusion Monte Carlo Perspective on the Spin-State Energetics of [Fe(NCH)6](2.)

The energy difference between the high spin and the low spin state of the model compound [Fe(NCH)6](2+) is investigated by means of Diffusion Monte Carlo (DMC), where special attention is dedicated to analyzing the effect of the fix node approximation on the accuracy of the results. For this purpose, we compare several Slater-Jastrow and multireference Slater-Jastrow trial wave functions. We found that a Slater-Jastrow trial wave function constructed with the generalized Kohn-Sham orbitals from hybrid DFT represents the optimal choice. This is understood by observing that hybrid functionals account for the subtle balance between exchange and correlation effects and the respective orbitals accurately describe the ligand-metal hybridization as well as the charge reorganization accompanying the spin transition. Finally the DMC results are compared with those obtained by Hartree-Fock, DFT, CASSCF, and CASPT2. While there is no clear reference value for the high spin-low spin energy difference, DMC and high level CCSD(T) calculations agree within around 0.3 eV.

Surface induces different crystal structures in a room temperature switchable spin crossover compound

We investigated the influence of surfaces in the formation of different crystal structures of a spin crossover compound, namely [Fe(L)2] (LH: (2-(pyrazol-1-yl)-6-(1H-tetrazol-5-yl)pyridine), which is a neutral compound thermally switchable around room temperature. We observed that the surface induces the formation of two different crystal structures, which exhibit opposite spin transitions, i.e. on heating them up to the transition temperature, one polymorph switches from high spin to low spin and the second polymorph switches irreversibly from low spin to high spin. We attributed this inversion to the presence of water molecules H-bonded to the complex tetrazolyl moieties in the crystals. Thin deposits were investigated by means of polarized optical microscopy, atomic force microscopy, X-ray diffraction, X-ray absorption spectroscopy and micro Raman spectroscopy; moreover the analysis of the Raman spectra and the interpretation of spin inversion were supported by DFT calculations.

TD-DFT study of the light-induced spin crossover of Fe(iii) complexes

Two light-induced spin-crossover Fe(iii) compounds have been studied with time-dependent density functional theory (TD-DFT) to investigate the deactivation mechanism and the role of the ligand-field states as intermediates in this process.

Abrupt versus Gradual Spin-Crossover in Fe(II)(phen)2(NCS)2 and Fe(III)(dedtc)3 Compared by X-ray Absorption and Emission Spectroscopy and Quantum-Chemical Calculations

Molecular spin-crossover (SCO) compounds are attractive for information storage and photovoltaic technologies. We compared two prototypic SCO compounds with Fe(II)N6 (1, [Fe(phen)2(NCS)2], with phen = 1,10-phenanthroline) or Fe(III)S6 (2, [Fe(dedtc)3], with dedtc = N,N'-diethyldithiocarbamate) centers, which show abrupt (1) or gradual (2) thermally induced SCO, using K-edge X-ray absorption and Kβ emission spectroscopy (XAS/XES) in a 8-315 K temperature range, single-crystal X-ray diffraction (XRD), and density functional theory (DFT). Core-to-valence and valence-to-core electronic transitions in the XAS/XES spectra and bond lengths change from XRD provided benchmark data, verifying the adequacy of the TPSSh/TZVP DFT approach for the description of low-spin (LS) and high-spin (HS) species. Determination of the spin densities, charge distributions, bonding descriptors, and valence-level configurations, as well as similar experimental and calculated enthalpy changes (ΔH), suggested that the varying metal-ligand bonding properties and deviating electronic structures converge to similar enthalpic contributions to the free-energy change (ΔG) and thus presumably are not decisive for the differing SCO behavior of 1 and 2. Rather, SCO seems to be governed by vibrational contributions to the entropy changes (ΔS) in both complexes. Intra- and intermolecular interactions in crystals of 1 and 2 were identified by atoms-in-molecules analysis. Thermal excitation of individual dedtc ligand vibrations accompanies the gradual SCO in 2. In contrast, extensive inter- and intramolecular phen/NCS vibrational mode coupling may be an important factor in the cooperative SCO behavior of 1.

Stabilization of the Low-Spin State in a Mononuclear Iron(II) Complex and High-Temperature Cooperative Spin Crossover Mediated by Hydrogen Bonding

The tetrapyridyl ligand bbpya (bbpya=N,N-bis(2,2'-bipyrid-6-yl)amine) and its mononuclear coordination compound [Fe(bbpya)(NCS)2 ] (1) were prepared. According to magnetic susceptibility, differential scanning calorimetry fitted to Sorai's domain model, and powder X-ray diffraction measurements, 1 is low-spin at room temperature, and it exhibits spin crossover (SCO) at an exceptionally high transition temperature of T1/2 =418 K. Although the SCO of compound 1 spans a temperature range of more than 150 K, it is characterized by a wide (21 K) and dissymmetric hysteresis cycle, which suggests cooperativity. The crystal structure of the LS phase of compound 1 shows strong NH⋅⋅⋅S intermolecular H-bonding interactions that explain, at least in part, the cooperative SCO behavior observed for complex 1. DFT and CASPT2 calculations under vacuum demonstrate that the bbpya ligand generates a stronger ligand field around the iron(II) core than its analogue bapbpy (N,N'-di(pyrid-2-yl)-2,2'-bipyridine-6,6'-diamine); this stabilizes the LS state and destabilizes the HS state in 1 compared with [Fe(bapbpy)(NCS)2 ] (2). Periodic DFT calculations suggest that crystal-packing effects are significant for compound 2, in which they destabilize the HS state by about 1500 cm(-1) . The much lower transition temperature found for the SCO of 2 compared to 1 appears to be due to the combined effects of the different ligand field strengths and crystal packing.

Towards an accurate and computationally-efficient modelling of Fe(ii)-based spin crossover materials

A theoretical approach is proposed to accurately calculate the LS–HS energy gap of SCO complexes in the solid state.

Spin crossover with thermal hysteresis: practicalities and lessons learnt

Memory applications of spin crossover require bistability: magnetic data must be appropriately collected and reported, and consideration given to lifetimes.