Uno de los cúmulos de metales de transición más estudiados teórica y experimentalmente son los de hierro (Fen). El creciente interés en estudiarlos se debe a sus propiedades químicas que permite su aplicación como catalizador en síntesis orgánica. Una manera de estudiar la reactividad del cúmulo en casos prácticos es mediante la adsorción de moléculas orgánicas, analizar las interacciones Fen‒ligante. En el presente estudio, se determinaron las propiedades estructurales y electrónicas del cúmulo de Fe6 y de compuestos N‒heterocíclicos piridinas y quinolinas, utilizando el funcional de intercambio y correlación PBE y la base orbital 6–311++G** para los átomos de los compuestos orgánicos, y el potencial de core efectivo (ECP) LanL2DZ para el átomo de Fe, se estudiaron tanto las propiedades del cúmulo de Fe6 y los compuestos orgánicos aislados como su interacción. La reactividad química del cúmulo y de los compuestos orgánicos se evaluó mediante diversas propiedades electrónicas que permitieron identificar factores que influyen en la interacción metal‒ligante, en función al tipo de sustituyente presente. Se aplicó la Teoría Cuántica de Átomos en Moléculas (QTAIM) para estudiar la naturaleza de la interacción metal–ligante revelando los puntos críticos de enlace y las interacciones débiles. Esto permitió comprender la naturaleza de unión metal‒ligante como parcialmente covalentes y parcialmente electrostáticas, formadas por los orbitales d del hierro y los enlaces π en los N‒heterociclos, de forma que se caracterizaron las interacciones entre el cúmulo de Fe6 y derivados de piridina como interacciones tipo η^1‒η^3 y en las quinolinas como η^2.