Asphaltene aggregation under vacuum at different temperatures was obtained using classical molecular dynamics (MD) simulations under nonperiodic boundary conditions in a monodisperse system of 96 hypothetical asphaltene molecules. Identical asphaltenes were originally set as an array, where the separation between each other was similar to40 Angstrom. Simulations under the canonical ensemble at NVT conditions, using the Verlet numerical method to solve the motion equations, were conducted. Aggregated systems formed by several asphaltene monomers after 100 ps of classical MD simulations were found. The structure of the solution was analyzed using the radial distribution function. Simulations at four different temperatures (273, 312, 342, and 368 K) were accomplished. Another similar MD simulation at a temperature of 310 K for 300 ps was performed, to validate the stability in the previous systems. After this run, good structures for explaining asphaltene interactions were also observed; these structures were never before proposed. The following effects can be observed from the results: (i) aggregates that have different structures indicate different types of interactions; (ii) decreasing aggregation number values with increasing temperature is consistent with experimental reality; (iii) the average molecular weights obtained for different temperatures agree with the expected range of experimental values; and (iv) the minimum of the potential energy well, in the range of 3.5-4.0 Angstrom, is consistent with the Yen model.