Assessing the performance of approximate density functional theory on 95 experimentally characterized Fe(II) spin crossover complexes
Name
024120_1_5.0157187.pdf
Description
Published version
Size
7.15 MB
Format
Adobe PDF
Checksum (MD5)
754c596dcf3adb7f5ce403b15ee788c9
Author(s)
Vennelakanti, Vyshnavi
Taylor, Michael G
Nandy, Aditya
Duan, Chenru
Kulik, Heather J
Date Issued
July 14, 2023
Journal
The Journal of Chemical Physics
Publisher
AIP Publishing
Citation
Vyshnavi Vennelakanti, Michael G. Taylor, Aditya Nandy, Chenru Duan, Heather J. Kulik; Assessing the performance of approximate density functional theory on 95 experimentally characterized Fe(II) spin crossover complexes. J. Chem. Phys. 14 July 2023; 159 (2): 024120.
Version
Final published version
Abstract
Spin crossover (SCO) complexes, which exhibit changes in spin state in response to external stimuli, have applications in molecular electronics and are challenging materials for computational design. We curate a dataset of 95 Fe(II) SCO complexes (SCO-95) from the Cambridge Structural Database that have available low- and high-temperature crystal structures and, in most cases, confirmed experimental spin transition temperatures (T1/2). We study these complexes using density functional theory (DFT) with 30 functionals spanning across multiple rungs of “Jacob’s ladder” to understand the effect of exchange–correlation functional on electronic and Gibbs free energies associated with spin crossover. We specifically assess the effect of varying the Hartree–Fock exchange fraction (aHF) in structures and properties within the B3LYP family of functionals. We identify three best-performing functionals, a modified version of B3LYP (aHF = 0.10), M06-L, and TPSSh, that accurately predict SCO behavior for the majority of the complexes. While M06-L performs well, MN15-L, a more recently developed Minnesota functional, fails to predict SCO behavior for all complexes, which could be the result of differences in datasets used for parametrization of M06-L and MN15-L and also the increased number of parameters for MN15-L. Contrary to observations from prior studies, double-hybrids with higher aHF values are found to strongly stabilize high-spin states and therefore exhibit poor performance in predicting SCO behavior. Computationally predicted T1/2 values are consistent among the three functionals but show limited correlation to experimentally reported T1/2 values. These failures are attributed to the lack of crystal packing effects and counter-anions in the DFT calculations that would be needed to account for phenomena such as hysteresis and two-step SCO behavior. The SCO-95 set thus presents opportunities for method development, both in terms of increasing model complexity and method fidelity.
MIT Department
Massachusetts Institute of Technology. Department of Chemical Engineering
Massachusetts Institute of Technology. Department of Chemistry
Terms of Use
Creative Commons Attribution
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1063/5.0157187
