In this study, in order to explain the solvent and spin state effects on the molecular structure of catechol-Fe complex [Fe(cat)3]n﹣ where n = 2 and 3, Hartree Fock (HF)-Density Functional Theory (DFT) hybrid calculations are performed at the B3LYP/6-311g(d,p) level of theory. The binding energies of Fe2+ and Fe3+ in high-spin state are higher than intermediate and low-spin states which show that the complex formation in a high spin state is more favorable. The calculated binding energies at different solvents indicate that the binding energies in polar solvents are lower than non-polar solvents. Furthermore, spectroscopic studies including FTIR and Raman spectrum in various solvents reveal that the formation of intermolecular bonds between the oxygen atom of carbonyl group and the hydrogen atom of solvent causes a spectral red shift. The calculated FTIR and geometry parameters are in good agreement with previous experimental data. Donor-acceptor interaction energies are evaluated due to the importance of the charge transfer in the complex formation. It is observed that the free electrons of oxygen atom interact with the antibonding orbitals of the iron. Finally, some correlations between the quantum chemical reactivity indices of the complexes and solvent polarity are considered. The study indicates a linear correlation between chemical hardness and binding energies of [Fe(cat)3]3﹣ complex.