en
Scientific article
Open access
English

Comparison of density functionals for energy and structural differences between the high- [5T2g:(t2g)4(eg)2] and low- [1A1g:(t2g)6(eg)0] spin states of iron(II) coordination compounds. II. More functionals and the hexaminoferrous cation, [Fe(NH3)6]2+

Published inThe Journal of chemical physics, vol. 122, no. 4, 44110
Publication date2005
Abstract

The ability of different density functionals to describe the structural and energy differences between the high- [5T2g:(t2g)4(eg)2] and low- [1A1g:(t2g)6(eg)0] spin states of small octahedral ferrous compounds is studied. This work is an extension of our previous study of the hexaquoferrous cation, [Fe(H2O)6]2+, [J. Chem. Phys. 120, 9473 (2004)] to include a second compound—namely, the hexaminoferrous cation, [Fe(NH3)6]2+—and several additional functionals. In particular, the present study includes the highly parametrized generalized gradient approximations (GGAs) known as HCTH and the meta-GGA VSXC [which together we refer to as highly parametrized density functionals (HPDFs)], now readily available in the GAUSSIAN03 program, as well as the hybrid functional PBE0. Since there are very few experimental results for these molecules with which to compare, comparison is made with best estimates obtained from second-order perturbation theory-corrected complete active space self-consistent field (CASPT2) calculations, with spectroscopy oriented configuration interaction (SORCI) calculations, and with ligand field theory (LFT) estimations. While CASPT2 and SORCI are among the most reliable ab initio methods available for this type of problem, LFT embodies many decades of empirical experience. These three methods are found to give coherent results and provide best estimates of the adiabatic low-spin–high-spin energy difference, DeltaELHadia, of 12 000–13 000 cm–1 for [Fe(H2O)6]2+ and 9 000–11 000 cm–1 for [Fe(NH3)6]2+. All functionals beyond the purely local approximation produce reasonably good geometries, so long as adequate basis sets are used. In contrast, the energy splitting, DeltaELHadia, is much more sensitive to the choice of functional. The local density approximation severely over stabilizes the low-spin state with respect to the high-spin state. This "density functional theory (DFT) spin pairing-energy problem" persists, but is reduced, for traditional GGAs. In contrast the hybrid functional B3LYP underestimates DeltaELHadia by a few thousands of wave numbers. The RPBE GGA of Hammer, Hansen, and Nørskov gives good results for DeltaELHadia as do the HPDFs, especially the VSXC functional. Surprisingly the HCTH functionals actually over correct the DFT spin pairing-energy problem, destabilizing the low-spin state relative to the high-spin state. Best agreement is found for the hybrid functional PBE0.

Keywords
  • Density functional theory
  • Perturbation theory
  • SCF calculations
  • Molecular electronic states
  • Configuration interactions
  • Ab initio calculations
  • Iron compounds
  • Positive ions
  • Bonds (chemical)
  • Molecular configurations
Research group
Citation (ISO format)
FOUQUEAU, Antony et al. Comparison of density functionals for energy and structural differences between the high- [5T2g:(t2g)4(eg)2] and low- [1A1g:(t2g)6(eg)0] spin states of iron(II) coordination compounds. II. More functionals and the hexaminoferrous cation, [Fe(NH3)6]2+. In: The Journal of chemical physics, 2005, vol. 122, n° 4, p. 44110. doi: 10.1063/1.1839854
Main files (1)
Article (Published version)
accessLevelPublic
Identifiers
ISSN of the journal0021-9606
770views
485downloads

Technical informations

Creation09/21/2009 4:22:37 PM
First validation09/21/2009 4:22:37 PM
Update time03/14/2023 3:14:08 PM
Status update03/14/2023 3:14:08 PM
Last indexation01/15/2024 6:57:40 PM
All rights reserved by Archive ouverte UNIGE and the University of GenevaunigeBlack