This study investigated the thermodynamic and dynamic properties of asphaltene molecules in Toluene and n-Heptane using the OPLS-AA force field in the GROMACS package at 300 K and 1 bar. The results indicate that van der Waals (vdW) interactions have a significant influence on thermodynamics, accounting for around 88–90 % of the total interaction energy between asphaltene molecules, while Coulombic (Coul) interactions contribute approximately 10–12 %. Stable asphaltene aggregates form at concentrations equal to or greater than 8 asphaltene molecules in 4000 molecules of Toluene. The interaction energy between asphaltene chains, which promotes aggregation, is much weaker than the interaction between aromatic rings. Longer chains exhibit more attractive vdW interactions and more repulsive Coul interactions. The presence of heteroatoms (sulfur and oxygen) hinders the parallel and symmetrical interaction of asphaltene molecules. Density profiles and two-dimensional maps were used to analyze the spatial arrangement of asphaltenes. A total interaction energy equal to or more positive than −300 kJ/mol for asphaltene-liquid interactions is introduced as a criterion for classifying a liquid as an asphaltene precipitator. The polarity of the solvent was inferred not to affect asphaltene solubility. In n-Heptane, asphaltene molecules aggregated to form a dense nucleus, while in Toluene, sharp peaks in the density profiles indicated the formation of clusters of varying sizes via asphaltene aggregation. Systems containing 64 and 32 asphaltene molecules had the highest probability of forming clusters comprising 10–16 and 12–16 asphaltene molecules, respectively. In Toluene, the average cluster sizes ranged from 1.4 to 6.1 molecules, depending on asphaltene concentration. Cluster size analysis showed variations in both size and stacking architectures, including face-to-face stacked and T-shaped aggregates. The proportion of face-to-face stacked structures was higher than T-shaped structures in both Toluene and n-Heptane, with the ratio of face-to-face stacked to T-shaped aggregates estimated at ∼2.9. While thermodynamic equilibrium was well attained, simulation times longer than 150 ns are recommended to achieve dynamic equilibrium.
Read full abstract