Table of Contents Author Guidelines Submit a Manuscript
Advances in Materials Science and Engineering
Volume 2016, Article ID 3817123, 10 pages
Research Article

Molecular Dynamics Simulation to Investigate the Interaction of Asphaltene and Oxide in Aggregate

1School of Highway, Chang’an University, Xi’an 710064, China
2School of Transportation Engineering, Southeast University, Nanjing 210096, China

Received 4 July 2015; Accepted 4 January 2016

Academic Editor: Amritendu Roy

Copyright © 2016 Rui Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The asphalt-aggregate interface interaction (AAI) plays a significant role in the overall performances of asphalt mixture, which is caused due to the complicated physicochemical processes and is influenced by various factors, including the acid-base property of aggregates. In order to analyze the effects of the chemical constitution of aggregate on the AAI, the average structure C65H74N2S2 is selected to represent the asphaltene in asphalt and magnesium oxide (MgO), calcium oxide (CaO), aluminium sesquioxide (Al2O3), and silicon dioxide (SiO2) are selected to represent the major oxides in aggregate. The molecular models are established for asphaltene and the four oxides, respectively, and the molecular dynamics (MD) simulation was conducted for the four kinds of asphaltene-oxide system at different temperatures. The interfacial energy in MD simulation is calculated to evaluate the AAI, and higher value means better interaction. The results show that interfacial energy between asphaltene and oxide reaches the maximum value at 25°C and 80°C and the minimum value at 40°C. In addition, the interfacial energy between asphaltene and MgO was found to be the greatest, followed by CaO, Al2O3, and SiO2, which demonstrates that the AAI between asphalt and alkaline aggregates is better than acidic aggregates.

1. Introduction

Studies have shown that asphalt pavement damage phenomenon has a strong dependency upon asphalt-aggregate interfacial nature. If the bond between asphalt and aggregate is weak, then the interface between asphalt and aggregate is prone to failure, which results in the early damage of pavement. Therefore, many methods were used to study the interface between asphalt and aggregate in different perspectives. Elphingstone [1] proposed that the surface free energy is an effective approach to predict the fatigue and moisture properties of mixture, and adhesion process of aggregate plays an important role. Tan et al. [2] tested the phase angle and complex shear modulus of matrix asphalt and asphalt mastic by dynamic shear rheometer, and the asphalt-aggregate interaction capabilities were evaluated by K.Ziegel-B parameter. Cheng et al. [3] measured the surface free energy of different asphalt and aggregates and calculated the adhesion work of asphalt and reactive aggregates with and without water. Bhasin et al. [4] proposed quantifying the free energy of the surface of asphalt-aggregate adhesion energy, and the results have shown that the same type of aggregate and bitumen adhesion between the energies has significant difference. To sum up, on the asphalt-aggregate interface studies have been largely analyzed from the experimental point of view. However, the experimental results were affected by many factors, in the dynamic shear rheological experiments; powder-binder ratio had greater impact on the test results, and how to determine the best powder-binder ratio is a problem. While aggregate composition is complex, powder, and mixed into the asphalt, mastic asphalt does not characterize asphalt-aggregate interface behavior. In the surface energy calculation, measurement results on the final great influence, different methods may result in different surface energy, resulting in a larger difference in the results.

In recent years, with the improvement of computer hardware and the perfection of related theory, molecular dynamics (MD) have become the effective means to study a complex system from the molecular level [5], the effective methods for the deep understanding of asphaltene adsorption characteristics, and the interaction mechanism of the interface from the microscopic view. At present, the domestic and foreign scholars mostly from the angle of test studied asphalt with aggregate interface area properties, and the computer simulation research is less. As a result, the MD method asphalt-aggregate interface behavior of asphalt mixture’s research has the vital significance.

Based on the molecular dynamics calculation, we can simulate the asphalt-aggregate interfacial behavior. In four components of asphalt, asphaltene molecular mass is bigger and has stronger polarity and the asphaltenes content has a great impact on the properties of asphalt [6]. Asphaltene is fused aromatic ring system as the core, with a number of different alkyl chains linked structure, with N, S, and O heteroatoms, not soluble in nonpolar molecules of the corresponding normal alkanes soluble in benzene and toluene matter [7]. The composition, structure, and morphology of the asphaltene are very complex. In the study, the average molecular model is often used to replace the chemical structure [8]. Rogel [9] proposed a representative average asphaltene molecular structure model for the study of this paper. According to the composition of aggregate analysis, aggregates content of more oxides was analyzed, the remaining content of less neglected ones. The asphalt-oxide interface model based on dynamics simulation of the interfacial energy and adsorption configuration, which can quantitatively point out the ability of the main oxides in asphalt and aggregate interface interaction, for further understanding the mechanism of adhesion between asphalt and aggregate, is of great significance.

2. Computer Simulation Theory

2.1. Basic Theory

Calculations were conducted with Accelrys Materials Studio software package in this study and this method relies mainly on Newtonian mechanics to simulate the movement of molecular systems, the samples are taken in different states of the molecular system, and the configuration of the system is integrated. The thermodynamic quantities and other macroscopic properties are further calculated based on the results of the integration configuration.

COMPASS is selected for force field simulation, and the potential energy function form is shown as follows:

In (1), the preceding five letters are bonding interaction of which the potential energy function contains bond lengths, bond angles, dihedral angles, plane, and cross-term; , , , and are expressed as bond lengths, bond angles, dihedral angles, and off-plane vibration; the next two letters are nonbonding interactions which are electrostatic interaction and van der Waals interaction , respectively. and are calculated by the following:

In (2) and (3), and represent different atoms; is atomic charge; is the distance of different atomic; is potential well depth.

The interface model is composed of three layers. The first layer is oxide, the second layer is asphaltene, and the third layer is vacuum, eliminating the impact of the cyclical interaction between the particles resulting from double counting.

The molecular dynamics simulation of the canonical ensemble (NVT) is completed by Discover module. According to the set temperature, the temperature is controlled by Andersen thermostat, the initial rate of each molecule is randomly generated by the Maxwell-Boltzmann distribution, Newton’s equation of motion is calculated by Velocity Verlet algorithm, and the van der Waals interaction and Coulomb interaction are calculated by Charge Group method.

2.2. Interfacial Energy Calculation

Interfacial energy is also known as the interaction energy, the calculation method is similar to the surface energy, and the difference is that the interfacial energy reflects the bonding strength between the two materials. The interface may be between the same materials or different materials; the interfacial energy of the same material can be theoretically equal to the surface energy. The interfacial energy can be calculated as shown in Figure 1.

Figure 1: Schematic diagram of the interfacial energy.

In Figure 1, to calculate the interfacial energy between materials and , the interface model is constructed as shown in the right of Figure 1. The model is calculated in temperature and system NVT by molecular dynamics until the system reaches a steady state; the energy of the model is . The interface models and are constructed as shown in the left of Figure 1 and calculated by molecular dynamics; the energy of the two models is and , respectively. If the area of the combination surface is , the interfacial energy is calculated as follows:

3. Results and Discussion

3.1. Asphaltene Molecular Structure Model

Since the composition and structure of the asphaltene are very complex, Rogel et al. built average molecular structure, which was standardized samples studied in the Strategic Highway Research Program (SHRP). The average structure is shown in Figure 2.

Figure 2: Average molecular structure of asphaltene.

Because asphaltenes are highly polydisperse molecules, the construction is highly complex. So in the research, we used one of the average molecules. For the simulation, we chose the COMPASS force field. For the asphaltene, the target density was set to 1 g/cm3, constructing the configuration of asphaltene system. In order to analyze the temperature dependence of the asphaltene system, six different temperatures were chosen: 10°C (283.15 K), 25°C (298.15 K), 40°C (313.15 K), 60°C (333.15 K), 80°C (353.15 K), and 100°C (373.15 K), and the results are shown in Figure 3.

Figure 3: The configuration of asphaltene molecules in different temperature.
3.2. Oxide Molecular Structure Model
3.2.1. Selection of Representative Oxides

To simulate aggregate, we used 10 chemical compositions of aggregate suggested by Ozkahraman and Işik [10] (as seen in Figure 4). The chemical composition of aggregate is SiO2 > CaO > Al2O3 > MgO > Fe2O3 > NaO > K2O.

Figure 4: The chemical composition of aggregates.

Therefore, the simulation is mainly focused on the four oxides which are SiO2, CaO, Al2O3, and MgO.

3.2.2. Building Oxide Molecular Cell

By the molecular dynamics software that comes with the unit cell model of oxide molecules, molecular cell structure of oxide is formed after importing the unit cell model, as shown in Figure 5.

Figure 5: Four kinds of oxide molecular cell structure.
3.2.3. Building Oxide Supercell

Based on the molecular cell above, the cell is intercepted by section 0 0 1. Secondly, the force field is assigned to each atom of cell. In the paper, the force field is COMPASS. The energy of cell is minimized by the Minimize module of Discover. Finally, the oxide supercell structure is built by the supercell module of Build which is composed of 5 × 5 matrix. For example, the process of building supercell structure of SiO2 is shown in Figure 6.

Figure 6: The process of building supercell structure.

3.3. The Interaction of Asphaltene and Oxide
3.3.1. Creating of Asphaltene/Oxide Interface

Interface model is composed of three layers. Layer 1 is oxide, layer 2 is asphaltene, and the last layer is 30 Å thick vacuum layer. According the depiction earlier, the three-layer interface structure was formed. For example, the asphaltene/SiO2 interface model of 25°C is shown in Figure 7.

Figure 7: Asphaltene/SiO2 interface model.
3.3.2. Dynamic Simulation Interface of Asphaltene/Oxide

As shown in the interface model in Figure 7, molecular dynamic simulation was then conducted in order to equilibrate systems further. The temperature is set and maintained constant by Andersen thermostat; the initial rate of molecular cell structure is randomly generated by the Maxwell-Boltzmann distribution. The Velocity Verlet algorithm is used to solve Newton’s equation of motion. The van der Waals interaction and the Coulomb interaction are calculated by Charge Group method. The intermolecular interaction is corrected by an average density approximation. The system was simulated under NVT conditions with velocity rescaling, which was used for 2000 ps with a time step of 1.0 fs to dissipate high energy interactions between molecules.

MD simulation of four-group interface model was performed for temperatures of 10°C (283.15 K), 25°C (298.15 K), 40°C (313.15 K), 60°C (333.15 K), 80°C (353.15 K), and 100°C (373.15 K). The results are shown in Figures 815.

Figure 8: The interface configuration of asphaltene/SiO2 (10–40°C).
Figure 9: The interface configuration of asphaltene/SiO2 (60–100°C).
Figure 10: The interface configuration of asphaltene/CaO (10–40°C).
Figure 11: The interface configuration of asphaltene/CaO (60–100°C).
Figure 12: The interface configuration of asphaltene/MgO (10–40°C).
Figure 13: The interface configuration of asphaltene/MgO (60–100°C).
Figure 14: The interface configuration of asphaltene/Al2O3 (10–40°C).
Figure 15: The interface configuration of asphaltene/Al2O3 (60–100°C).

As shown in the four interface configurations, the asphaltene molecular structures change in varying degree after the simulation. The asphaltene molecular structure moves closer to the oxide surface. Comparing the four kinds of interface configuration, it is difficult to evaluate the interaction between asphaltene and oxides accurately. For the fixed interface configuration, it is almost impossible to distinguish whether the interaction changes with temperature. Therefore, the interfacial energy is introduced to evaluate the interaction more precisely.

3.3.3. The Calculation of Interfacial Energy

The interfacial energy can be calculated according to the following:

The four-interface model size is 49.000 × 49.000 × 26.38067 Å. So the interfacial energy equation can be simplified as follows:

After the molecular dynamics simulation, the four kinds of interfacial energy are shown in Tables 1 and 2.

Table 1: Asphaltene/SiO2 and asphaltene/CaO.
Table 2: Asphaltene/MgO and asphaltene/Al2O3.

The interfacial energy of asphaltene/SiO2 varies from to  kJ/mol within the temperature range of 10°C to 100°C; the calculated variation ranges of asphaltene/CaO, asphaltene/MgO, and asphaltene/Al2O3 within the same temperature range are to  kJ/mol, to  kJ/mol, and to  kJ/mol, respectively. Further, the interfacial energy was plotted in Figure 16.

Figure 16: The interfacial energy with different temperature.

As shown in Figure 16, the four kinds of interfacial energy reach the extreme twice at 25°C and 80°C which indicates stronger interaction between asphaltene and oxides. Within the temperature range of 10°C to 100°C, the interaction capability ranking of the four oxides to asphaltene is as follows: MgO > CaO > Al2O3 > SiO2.

4. Conclusion

The interaction and the interfacial energy between asphaltene and four kinds of oxides (MgO, CaO, Al2O3, and SiO2) under temperature range from 10°C to 100°C were analyzed by virtue of molecular dynamic simulations and the following conclusions can be made: the interaction between asphaltene and oxides reaches strongest level at the temperature of 25°C and 80°C and reaches weakest level at the temperature of 40°C. It could be concluded according to the calculation results of the interfacial energy that the interfacial energy between asphaltene and four kinds of oxides ranges from  kJ/mol to  kJ/mol within the temperature range of 10°C to 100°C. In this temperature range, the interaction capability ranking of the four oxides to asphaltene is as follows: MgO > CaO > Al2O3 > SiO2.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work is supported by the National Natural Science Foundation of China (Grant nos. 51408048 and 51378037), the Key Program for International Science and Technology Cooperation Projects of Shaanxi province (Grant nos. 2012KW-09 and 2014KW10-03), the China Postdoctoral Science Foundation (Grant nos. 2015M570805 and 2015M582592), the Fundamental Research Funds for the Central Universities of Chang’an University (Grant nos. 310821153502 and 310821152005), and the Research Fund for the Doctoral Program of Higher Education of China (Grant no. 20120205120010).


  1. G. M. Elphingstone, Adhesion and cohesion in asphalt-aggregate systems [Ph.D. thesis], Texas A&M University, College Station, Tex, USA, 1997.
  2. Y.-Q. Tan, J.-T. Wu, X.-M. Li, and L. Ji, “Evaluation indexes for interaction capability of asphalt and aggregate,” Journal of Harbin Institute of Technology, vol. 41, no. 7, pp. 81–84, 2009. View at Google Scholar · View at Scopus
  3. D. Cheng, D. N. Little, R. L. Lytton, and J. C. Holste, “Use of surface free energy properties of the asphalt-aggregate system to predict moisture damage potential,” Journal of the Association of Asphalt Paving Technologists, vol. 71, pp. 59–88, 2002. View at Google Scholar
  4. A. Bhasin, E. Masad, D. Little, and R. Lytton, “Limits on adhesive bond energy for improved resistance of hot-mix asphalt to moisture damage,” Transportation Research Record, vol. 1970, pp. 3–13, 2006. View at Publisher · View at Google Scholar
  5. R. Notman and T. R. Walsh, “Molecular dynamics studies of the interactions of water and amino acid analogues with quartz surfaces,” Langmuir, vol. 25, no. 3, pp. 1638–1644, 2009. View at Publisher · View at Google Scholar · View at Scopus
  6. K. J. Leontaritis, J. O. Amaefule, and R. E. Charles, “Systematic approach for the prevention and treatment of formation damage caused by asphaltene deposition,” in Proceedings of the International Symposium on Formation Damage Control, pp. 383–395, Lafayette, La, USA, February 1992. View at Scopus
  7. Y.-F. Li, G.-W. Lu, W. Sun et al., “Study on the molecular dynamics of petroleum-derived asphaltene aggregate,” Acta Petrolei Sinica (Petroleum Processing Section), vol. 23, no. 4, pp. 25–31, 2007. View at Google Scholar
  8. P. Srivastava, W. G. Chapman, and P. E. Laibinis, “Molecular dynamics simulation of oxygen transport through n-alkanethiolate self-assembled monolayers on gold and copper,” The Journal of Physical Chemistry B, vol. 113, no. 2, pp. 456–464, 2009. View at Publisher · View at Google Scholar · View at Scopus
  9. E. Rogel, “Studies on asphaltene aggregation via computational chemistry,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 104, no. 1, pp. 85–93, 1995. View at Publisher · View at Google Scholar · View at Scopus
  10. H. T. Ozkahraman and E. C. Işik, “The effect of chemical and mineralogical composition of aggregates on tensile adhesion strength of tiles,” Construction and Building Materials, vol. 19, no. 4, pp. 251–255, 2005. View at Publisher · View at Google Scholar · View at Scopus