“…In this context, insights provided by molecular modeling and simulations have all their attractiveness and importance. Traditionally, computational methods have been used to predict structures, describe spin-crossover phenomena, [77][78][79] predict electrochemical properties, [80,81] rationalize optical properties, [60,[82][83][84][85][86] and elucidate the nature and relative alignment of MLCT and MC states and the barriers for their mutual interconversion. [87,88] The large size of the iron-organic complexes often imposes the use of DFT and TD-DFT as methods of choice for tackling this problem, although the development and use of other methodologies derived from multiconfigurational approaches (sometimes in combination with DFT, as in hybrid DFT/wavefunction formalisms) is growing in the last years.…”