Computational Investigation of Ionization Energies and Absorption Spectra of Metallocenes: Effect of the Metal Center on Electronic Properties

Since the synthesis of ferrocene in 1951, metallocenes have attracted attention, making the accurate prediction of their electronic structure and ionization energy crucial for understanding their photophysical and electrochemical behavior in materials and in biological systems. Here, we combined Density Functional Theory (DFT), Complete Active Space Self-Consistent Field (CASSCF), NEVPT2 (N-Electron Valence State Perturbation Theory) and Coupled Cluster approaches (CCSD, DLPNO-CCSD(T)) to study the electronic structure, ionization energies (IEs) and absorption spectra of metallocene and metallocenium complexes in the gas phase and in THF implicit solvent. DFT IEs agree closely with NEVPT2 and DLPNO-CCSD(T) values and with experiment values (deviations 0.02–0.3 eV). For CASSCF and NEVPT2, the minimal active space of the d electrons at six orbitals is not enough for the accurate prediction of the IEs, while an extended active space incorporating all 3d metal electrons plus four ligand valence electrons into 15 orbitals improves the calculated IE values. In solution, computed oxidation energies (OEs) in THF reproduce experimental values and follow the Fe > Ni > Co ordering. Substitution of metallocene complexes with chromophore units results in similar OEs. Overall, the substitution effects remain modest: the effect of substitution on OE values results in differences up to 0.2 eV. These results clarify the effect of the metal center on IE and OE values and UV–vis absorption behavior.