Time-dependent density functional theory/discrete reaction field spectra of open shell systems: The visual spectrum of [FeIII(PyPepS)2]- in aqueous solution

We report the calculated visible spectrum of FeIII(PyPepS)2]- in aqueous solution. From all-classical molecular dynamics simulations on the solute and 200 water molecules with a polarizable force field, 25 solute/solvent configurations were chosen at random from a 50 ps production run and subjected the systems to calculations using time-dependent density functional theory (TD-DFT) for the solute, combined with a solvation model in which the water molecules carry charges and polarizabilities. In each calculation the first 60 excited states were collected in order to span the experimental spectrum. Since the solute has a doublet ground state several excitations to states are of type "three electrons in three orbitals," each of which gives rise to a manifold of a quartet and two doublet states which cannot properly be represented by single Slater determinants. We applied a tentative scheme to analyze this type of spin contamination in terms of Delta and Delta transitions between the same orbital pairs. Assuming the associated states as pure single determinants obtained from restricted calculations, we construct conformation state functions (CFSs), i.e., eigenfunctions of the Hamiltonian Sz and S2, for the two doublets and the quartet for each Delta,Delta pair, the necessary parameters coming from regular and spin-flip calculations. It appears that the lower final states remain where they were originally calculated, while the higher states move up by some tenths of an eV. In this case filtering out these higher states gives a spectrum that compares very well with experiment, but nevertheless we suggest investigating a possible (re)formulation of TD-DFT in terms of CFSs rather than determinants.