Abstract: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 ([Fe(bpy)3]2+, [Fe(dafo)3]2+), ([Fe(bim)3]2+, [Fe(xbim)3]2+) and ([Fe(bimz)3]2+, [Fe(etbimz)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 [Fe(NCH)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 [Fe(NCH)6]2+ with the tetragonal arrangement of its nitrile ligands to the tris-diimine complexes with the trigonal packing of their bulkier ligands.

Benchmarks and reliable DFT results for spin-crossover complexes

Performance of Quantum Chemistry Methods for Benchmark Set of Spin–State Energetics Derived from Experimental Data of 17 Transition Metal Complexes (SSE17)

DC-DFT for Open Shells: How to Deal with Spin Contamination

Can DFT Accurately Predict Spin Densities? Analysis of Discrepancies in Iron Nitrosyl Complexes

Density Functional Theory Calculations for Spin Crossover Complexes

Spin-state Gaps and Self-Interaction-Corrected Density Functional Approximations: Octahedral Fe(II) Complexes as Case Study

Predicting Spin States of Iron Porphyrins with DFT Methods Including Crystal Packing Effects and Thermodynamic Corrections

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

Conceptual Insights into DFT Spin-State Energetics of Octahedral Transition-Metal Complexes through a Density Difference Analysis.

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state

Benchmark Study of Density Functional Theory Methods in Geometry Optimization of Transition Metal-Dinitrogen Complexes

Improvements in DFT Calculations of Spin–Spin Coupling Constants

Appraising spin-state energetics in transition metal complexes using double-hybrid models: accountability of SOS0-PBESCAN0-2(a) as a promising paradigm

Approaching the Complete Basis Set Limit for Transition Metal Spin–State Energetics

Coupled-Cluster Treatment of Complex Open-Shell Systems: The Case of Single-Molecule Magnets

Low-rank approximations for accelerating plane-wave hybrid functional calculations in unrestricted and noncollinear spin density functional theory

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

Improving results by improving densities: Density-corrected density functional theory