En este trabajo se describe el estudio y reactividad de cúmulos de oro Aun (n=3−11) hacia compuestos derivados de anisol y tioanisol utilizando el funcional VWN y la base def2-TZVP para la optimización de las geometrías y el funcional PBE y la base def2-TZVP para el cálculo de los parámetros energéticos. En este estudio se efectuó la optimización local de distintas estructuras posibles para los cúmulos de oro Aun (n=3−11) y se seleccionaron los cúmulos más estables para su interacción con los compuestos orgánicos. Para la optimización de los compuestos orgánicos derivados de anisol y tioanisol se empleó el funcional PBE y la base 6−311G**. La reactividad química fue analizada en función de los orbitales moleculares de frontera, la dureza, la diferencia de energía LUMO-HOMO y la función de Fukui, los cuales permitieron seleccionar las moléculas orgánicas y los sitios adecuados para realizar la interacción con los cúmulos de oro.