El presente trabajo se examina los patrones de unión entre los grupos de cobre Cun (n = 3 − 20) y compuestos aromáticos (benceno, fenol y benzaldehído) utilizando un enfoque de teoría de densidad funcional (DFT) y la teoría de átomos en moléculas (AIM). En este estudio se efectuó una búsqueda estocástica de distintas estructuras para los cúmulos de cobre mediante recocido molecular, optimizando las estructuras con el funcional PBE y la base LANL2DZ para obtener los cúmulos más estables. A su vez los compuestos aromáticos se optimizaron empleando el mismo funcional y la base 6-311++G**. La reactividad química, tanto de los cúmulos como de los compuestos aromáticos, se analizó en función de los orbitales moleculares de frontera, dureza, diferencia HOMO-LUMO, segunda diferencia de energía y la función de Fukui, los cuales permitieron determinar las zonas adecuadas para realizar la interacción metal-ligante.