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

Authors: Rui Li; Hui Du; Zepeng Fan; Jianzhong Pei
Journal: Advances in Materials Science and Engineering
 (2016)

---

## Abstract

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.


---

## Body

## 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:(1)Etotal=∑bondEbb+∑angleEθθ+∑dihedralEϕϕ+∑out-of-planeEχχ+∑crossEb,θ,ϕ+Ecoulomb+Evdw.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; b, θ, ϕ, 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 Ecoulomb and van der Waals interaction Evdw, respectively. Ecoulomb and Evdw are calculated by the following: (2)Ecoulomb=∑i>jqiqjrij,(3)Evdw=∑εij2rij0rij9-3rij0rij6.In (2) and (3), i and j represent different atoms; q is atomic charge; r 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 Figure1.Figure 1
Schematic diagram of the interfacial energy.In Figure1, to calculate the interfacial energy between materials M and N, the interface model is constructed as shown in the right of Figure 1. The model is calculated in temperature T and system NVT by molecular dynamics until the system reaches a steady state; the energy of the model is Etotal. The interface models M and N are constructed as shown in the left of Figure 1 and calculated by molecular dynamics; the energy of the two models is EM and EN, respectively. If the area of the combination surface is A, the interfacial energy is calculated as follows:(4)γint=EM+EN-Etotal2A.


## 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:(1)Etotal=∑bondEbb+∑angleEθθ+∑dihedralEϕϕ+∑out-of-planeEχχ+∑crossEb,θ,ϕ+Ecoulomb+Evdw.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; b, θ, ϕ, 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 Ecoulomb and van der Waals interaction Evdw, respectively. Ecoulomb and Evdw are calculated by the following: (2)Ecoulomb=∑i>jqiqjrij,(3)Evdw=∑εij2rij0rij9-3rij0rij6.In (2) and (3), i and j represent different atoms; q is atomic charge; r 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 Figure1.Figure 1
Schematic diagram of the interfacial energy.In Figure1, to calculate the interfacial energy between materials M and N, the interface model is constructed as shown in the right of Figure 1. The model is calculated in temperature T and system NVT by molecular dynamics until the system reaches a steady state; the energy of the model is Etotal. The interface models M and N are constructed as shown in the left of Figure 1 and calculated by molecular dynamics; the energy of the two models is EM and EN, respectively. If the area of the combination surface is A, the interfacial energy is calculated as follows:(4)γint=EM+EN-Etotal2A.


## 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 Figure2.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(M1,M2,…,M10) 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 Figure5.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 Figure7, 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 Figures8–15.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:(5)γint=Epolymer+Esurface-Etotal2A.The four-interface model size is 49.000 × 49.000 × 26.38067 Å. So the interfacial energy equation can be simplified as follows:(6)γint′=Epolymer+Esurface-Etotal.After the molecular dynamics simulation, the four kinds of interfacial energy are shown in Tables1 and 2.Table 1
Asphaltene/SiO2 and asphaltene/CaO.
Temp. (°C)
E
total (kJ/mol)
E
polymer (kJ/mol)
E
surface (kJ/mol)
E
Interface (kJ/mol)
Asphaltene/SiO2
10
9.03
E
+
06
3.63
E
+
02
9.03
E
+
06
2.22
E
+
01
25
9.03
E
+
06
3.83
E
+
02
9.03
E
+
06
2.09
E
+
02
40
9.03
E
+
06
3.73
E
+
02
9.03
E
+
06
3.25
E
+
01
60
9.03
E
+
06
3.72
E
+
02
9.03
E
+
06
3.88
E
+
01
80
9.03
E
+
06
3.83
E
+
02
9.03
E
+
06
1.36
E
+
02
100
9.03
E
+
06
3.85
E
+
02
9.03
E
+
06
1.31
E
+
02
Asphaltene/CaO
10
-
1.57
E
+
06
3.71
E
+
03
-
1.57
E
+
06
8.99
E
+
02
25
-
1.57
E
+
06
3.69
E
+
03
-
1.57
E
+
06
9.63
E
+
02
40
-
1.57
E
+
06
3.68
E
+
03
-
1.57
E
+
06
7.44
E
+
02
60
-
1.57
E
+
06
3.87
E
+
03
-
1.57
E
+
06
8.51
E
+
02
80
-
1.57
E
+
06
3.83
E
+
03
-
1.57
E
+
06
1.00
E
+
03
100
-
1.57
E
+
06
3.81
E
+
03
-
1.57
E
+
06
8.13
E
+
02Table 2
Asphaltene/MgO and asphaltene/Al2O3.
Temp. (°C)
E
total (kJ/mol)
E
polymer (kJ/mol)
E
surface (kJ/mol)
E
Interface (kJ/mol)
Asphaltene/MgO
10
-
2.61
E
+
06
3.71
E
+
03
-
2.61
E
+
06
1.67
E
+
03
25
-
2.61
E
+
06
3.79
E
+
03
-
2.61
E
+
06
1.79
E
+
03
40
-
2.61
E
+
06
3.73
E
+
03
-
2.61
E
+
06
1.61
E
+
03
60
-
2.61
E
+
06
3.76
E
+
03
-
2.61
E
+
06
1.61
E
+
03
80
-
2.61
E
+
06
3.76
E
+
03
-
2.61
E
+
06
1.68
E
+
03
100
-
2.61
E
+
06
3.81
E
+
03
-
2.61
E
+
06
1.63
E
+
03
Asphaltene/Al2O3
10
-
2.01
E
+
07
3.75
E
+
03
-
2.01
E
+
07
4.61
E
+
02
25
-
2.01
E
+
07
3.72
E
+
03
-
2.01
E
+
07
4.32
E
+
02
40
-
2.01
E
+
07
3.92
E
+
03
-
2.01
E
+
07
1.19
E
+
02
60
-
2.01
E
+
07
3.98
E
+
03
-
2.01
E
+
07
3.10
E
+
02
80
-
2.01
E
+
07
3.97
E
+
03
-
2.01
E
+
07
6.43
E
+
02
100
-
2.01
E
+
07
3.87
E
+
03
-
2.01
E
+
07
2.71
E
+
02The interfacial energy of asphaltene/SiO2 varies from 2.22E+01 to 2.09E+02 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 7.44E+02 to 1.00E+03 kJ/mol, 1.61E+03 to 1.79E+03 kJ/mol, and 1.19E+02 to 6.43E+02 kJ/mol, respectively. Further, the interfacial energy was plotted in Figure 16.Figure 16
The interfacial energy with different temperature.As shown in Figure16, 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.


## 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 Figure2.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(M1,M2,…,M10) 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 Figure5.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.2.1. Selection of Representative Oxides

To simulate aggregate, we used 10 chemical compositions of aggregate(M1,M2,…,M10) 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 Figure5.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 Figure7, 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 Figures8–15.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:(5)γint=Epolymer+Esurface-Etotal2A.The four-interface model size is 49.000 × 49.000 × 26.38067 Å. So the interfacial energy equation can be simplified as follows:(6)γint′=Epolymer+Esurface-Etotal.After the molecular dynamics simulation, the four kinds of interfacial energy are shown in Tables1 and 2.Table 1
Asphaltene/SiO2 and asphaltene/CaO.
Temp. (°C)
E
total (kJ/mol)
E
polymer (kJ/mol)
E
surface (kJ/mol)
E
Interface (kJ/mol)
Asphaltene/SiO2
10
9.03
E
+
06
3.63
E
+
02
9.03
E
+
06
2.22
E
+
01
25
9.03
E
+
06
3.83
E
+
02
9.03
E
+
06
2.09
E
+
02
40
9.03
E
+
06
3.73
E
+
02
9.03
E
+
06
3.25
E
+
01
60
9.03
E
+
06
3.72
E
+
02
9.03
E
+
06
3.88
E
+
01
80
9.03
E
+
06
3.83
E
+
02
9.03
E
+
06
1.36
E
+
02
100
9.03
E
+
06
3.85
E
+
02
9.03
E
+
06
1.31
E
+
02
Asphaltene/CaO
10
-
1.57
E
+
06
3.71
E
+
03
-
1.57
E
+
06
8.99
E
+
02
25
-
1.57
E
+
06
3.69
E
+
03
-
1.57
E
+
06
9.63
E
+
02
40
-
1.57
E
+
06
3.68
E
+
03
-
1.57
E
+
06
7.44
E
+
02
60
-
1.57
E
+
06
3.87
E
+
03
-
1.57
E
+
06
8.51
E
+
02
80
-
1.57
E
+
06
3.83
E
+
03
-
1.57
E
+
06
1.00
E
+
03
100
-
1.57
E
+
06
3.81
E
+
03
-
1.57
E
+
06
8.13
E
+
02Table 2
Asphaltene/MgO and asphaltene/Al2O3.
Temp. (°C)
E
total (kJ/mol)
E
polymer (kJ/mol)
E
surface (kJ/mol)
E
Interface (kJ/mol)
Asphaltene/MgO
10
-
2.61
E
+
06
3.71
E
+
03
-
2.61
E
+
06
1.67
E
+
03
25
-
2.61
E
+
06
3.79
E
+
03
-
2.61
E
+
06
1.79
E
+
03
40
-
2.61
E
+
06
3.73
E
+
03
-
2.61
E
+
06
1.61
E
+
03
60
-
2.61
E
+
06
3.76
E
+
03
-
2.61
E
+
06
1.61
E
+
03
80
-
2.61
E
+
06
3.76
E
+
03
-
2.61
E
+
06
1.68
E
+
03
100
-
2.61
E
+
06
3.81
E
+
03
-
2.61
E
+
06
1.63
E
+
03
Asphaltene/Al2O3
10
-
2.01
E
+
07
3.75
E
+
03
-
2.01
E
+
07
4.61
E
+
02
25
-
2.01
E
+
07
3.72
E
+
03
-
2.01
E
+
07
4.32
E
+
02
40
-
2.01
E
+
07
3.92
E
+
03
-
2.01
E
+
07
1.19
E
+
02
60
-
2.01
E
+
07
3.98
E
+
03
-
2.01
E
+
07
3.10
E
+
02
80
-
2.01
E
+
07
3.97
E
+
03
-
2.01
E
+
07
6.43
E
+
02
100
-
2.01
E
+
07
3.87
E
+
03
-
2.01
E
+
07
2.71
E
+
02The interfacial energy of asphaltene/SiO2 varies from 2.22E+01 to 2.09E+02 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 7.44E+02 to 1.00E+03 kJ/mol, 1.61E+03 to 1.79E+03 kJ/mol, and 1.19E+02 to 6.43E+02 kJ/mol, respectively. Further, the interfacial energy was plotted in Figure 16.Figure 16
The interfacial energy with different temperature.As shown in Figure16, 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.


## 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 Figure7, 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 Figures8–15.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:(5)γint=Epolymer+Esurface-Etotal2A.The four-interface model size is 49.000 × 49.000 × 26.38067 Å. So the interfacial energy equation can be simplified as follows:(6)γint′=Epolymer+Esurface-Etotal.After the molecular dynamics simulation, the four kinds of interfacial energy are shown in Tables1 and 2.Table 1
Asphaltene/SiO2 and asphaltene/CaO.
Temp. (°C)
E
total (kJ/mol)
E
polymer (kJ/mol)
E
surface (kJ/mol)
E
Interface (kJ/mol)
Asphaltene/SiO2
10
9.03
E
+
06
3.63
E
+
02
9.03
E
+
06
2.22
E
+
01
25
9.03
E
+
06
3.83
E
+
02
9.03
E
+
06
2.09
E
+
02
40
9.03
E
+
06
3.73
E
+
02
9.03
E
+
06
3.25
E
+
01
60
9.03
E
+
06
3.72
E
+
02
9.03
E
+
06
3.88
E
+
01
80
9.03
E
+
06
3.83
E
+
02
9.03
E
+
06
1.36
E
+
02
100
9.03
E
+
06
3.85
E
+
02
9.03
E
+
06
1.31
E
+
02
Asphaltene/CaO
10
-
1.57
E
+
06
3.71
E
+
03
-
1.57
E
+
06
8.99
E
+
02
25
-
1.57
E
+
06
3.69
E
+
03
-
1.57
E
+
06
9.63
E
+
02
40
-
1.57
E
+
06
3.68
E
+
03
-
1.57
E
+
06
7.44
E
+
02
60
-
1.57
E
+
06
3.87
E
+
03
-
1.57
E
+
06
8.51
E
+
02
80
-
1.57
E
+
06
3.83
E
+
03
-
1.57
E
+
06
1.00
E
+
03
100
-
1.57
E
+
06
3.81
E
+
03
-
1.57
E
+
06
8.13
E
+
02Table 2
Asphaltene/MgO and asphaltene/Al2O3.
Temp. (°C)
E
total (kJ/mol)
E
polymer (kJ/mol)
E
surface (kJ/mol)
E
Interface (kJ/mol)
Asphaltene/MgO
10
-
2.61
E
+
06
3.71
E
+
03
-
2.61
E
+
06
1.67
E
+
03
25
-
2.61
E
+
06
3.79
E
+
03
-
2.61
E
+
06
1.79
E
+
03
40
-
2.61
E
+
06
3.73
E
+
03
-
2.61
E
+
06
1.61
E
+
03
60
-
2.61
E
+
06
3.76
E
+
03
-
2.61
E
+
06
1.61
E
+
03
80
-
2.61
E
+
06
3.76
E
+
03
-
2.61
E
+
06
1.68
E
+
03
100
-
2.61
E
+
06
3.81
E
+
03
-
2.61
E
+
06
1.63
E
+
03
Asphaltene/Al2O3
10
-
2.01
E
+
07
3.75
E
+
03
-
2.01
E
+
07
4.61
E
+
02
25
-
2.01
E
+
07
3.72
E
+
03
-
2.01
E
+
07
4.32
E
+
02
40
-
2.01
E
+
07
3.92
E
+
03
-
2.01
E
+
07
1.19
E
+
02
60
-
2.01
E
+
07
3.98
E
+
03
-
2.01
E
+
07
3.10
E
+
02
80
-
2.01
E
+
07
3.97
E
+
03
-
2.01
E
+
07
6.43
E
+
02
100
-
2.01
E
+
07
3.87
E
+
03
-
2.01
E
+
07
2.71
E
+
02The interfacial energy of asphaltene/SiO2 varies from 2.22E+01 to 2.09E+02 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 7.44E+02 to 1.00E+03 kJ/mol, 1.61E+03 to 1.79E+03 kJ/mol, and 1.19E+02 to 6.43E+02 kJ/mol, respectively. Further, the interfacial energy was plotted in Figure 16.Figure 16
The interfacial energy with different temperature.As shown in Figure16, 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 2.22E+01 kJ/mol to 1.79E+03 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.



---

*Source: 3817123-2016-02-17.xml*