Unveiling the role of heteroatom compounds in heavy oil viscosity is pivotal for finding targeted viscosity reduction methods to improve oil recovery. This research investigates the impact of heteroatoms in asphaltene molecules by utilizing quantum chemical calculations and molecular dynamics simulations to analyze their electrostatic potential characteristics, pairwise interactions, and dynamic behavior within realistic reservoirs. Heteroatom compounds can influence the molecular-level properties of asphaltenes and thus impact the macroscopic behavior of heavy oils. Research results suggest that the presence of ketone and aromatic rings in asphaltene molecules leads to the unrestricted movement of pi electrons due to their collective electronegativity. Two distinct configurations of asphaltene dimers, face-to-face, and edge-to-face, were observed. Intermolecular interactions were predominantly governed by van der Waals forces, highlighting their significant role in stabilizing asphaltene aggregates. The distribution of asphaltene molecules in the oil phase can be summarized as the "rebar-cement" theory. In the heteroatom-free system, the face-to-face peaks in the radial distribution function exhibit significantly reduced magnitudes compared to those in the heteroatom-containing system. This emphasizes the pivotal function of heteroatoms in connecting molecular components to form a more compact asphaltene structure, which may result in a higher viscosity of heavy oil. These findings give insight into the significance of heteroatoms in bridging molecular components and shaping the intricate structure of asphaltene and advance our understanding of heavy oil viscosity properties.