Um método teórico geral em duas etapas para estudar redistribuições eletrônicas em processos catalíticos é apresentado. Na primeira etapa, teoria do funcional da densidade (DFT) é usada para otimizar completamente duas geometrias: o aglomerado que representa o catalisador e o sistema aglomerado mais molécula adsorvida. Na segunda etapa, a densidade eletrônica molecular convergida é dividida em multipolos centrados em sítios atômicos obtidos segundo uma análise de multipolos distribuídos, que fornece informação topológica detalhada da redistribuição de carga do catalisador e da molécula, antes e depois da adsorção. Esta abordagem é aplicada para a adsorção de tiofeno sobre a borda metálica 10 -10 dos catalisadores Ni(Co)MoS e comparada com a mesma reação em MoS 2 não-substituído. Energias de adsorção, geometrias e a análise de multipolos calculada indicam uma fraca adsorção do tiofeno em ambas as situações. Um modelo coulombiano para a ligação química mostra que a magnitude da ligação metal-enxofre na superfície dos catalisadores de Ni(Co)MoS promovidos é consideravelmente menor do que em MoS 2 não substituído, desta forma confirmando a origem do aumento da atividade de hidrodessulfurização (HDS) nestes catalisadores.A general two-step theoretical approach to study electronic redistributions in catalytic processes is presented. In the first step, density functional theory (DFT) is used to fully optimize two geometries: the cluster representing the catalyst and the cluster plus adsorbed molecule system. In the second step, the converged electron density is divided into multipoles centered on atomic sites according to a distributed multipole analysis which provides detailed topological information on the charge redistribution of catalyst and molecule before and after adsorption. This approach is applied to thiophene adsorption on the 10 -10 metal edge of Ni(Co)MoS catalysts and compared to the same reaction on bare MoS 2 . Calculated adsorption energies, geometries and multipole analysis indicate weak thiophene chemisorption on both cases. A Coulombic bond model showed that surface metal-sulfur bond strengths in Ni(Co)MoS promoted catalysts are considerably smaller than in bare MoS 2 , thus confirming the origin of the enhancement of hydrodesulfurization (HDS) activity in these catalysts.