Graphene-modified asphalt nanocomposite has shown an advanced self-healing ability compared to pure asphalt. In this protocol, molecular dynamics simulations have been applied in order to understand the role of graphene in the self-healing process and to explore the self-healing mechanism of asphalt components from the atomistic level.
Graphene can improve the self-healing properties of asphalt with high durability. However, the self-healing behaviors of graphene-modified asphalt nanocomposite and the role of incorporated graphene are still unclear at this stage. In this study, the self-healing properties of pure asphalt and graphene-modified asphalt are investigated through molecular dynamics simulations. Asphalt bulks with two crack widths and locations for graphene are introduced, and the molecular interactions among asphalt components and the graphene sheet are analyzed. The results show that the location of graphene significantly affects the self-healing behaviors of asphalt. Graphene near the crack surface can greatly accelerate the self-healing process by interacting with the aromatic molecules through π-π stacking, while graphene at the top area of the crack tip has a minor impact on the process. The self-healing process of asphalt goes through the reorientation of asphaltene, polar aromatic, and naphthene aromatic molecules, and the bridging of saturate molecules between crack surfaces. This in-depth understanding of the self-healing mechanism contributes to the knowledge of the enhancement for self-healing properties, which will help to develop durable asphalt pavements.
Deterioration under daily vehicle loadings and variant environmental conditions, and the aging of asphalt during service result in degradation or even structural failures, i.e., cracking and raveling, which can further weaken the durability of asphalt pavements. The inherent response of asphalt to repair micro-cracks and voids automatically helps it recover from damages and restore strength1. This self-healing capability can considerably extend the service life of asphalt, save costs on maintenance, and reduce the emission of greenhouse gases2,3. The self-healing behavior of asphalt generally depends on several influencing factors, including its chemical composition, the degree of damage, and environmental conditions4. The improved self-healing capability of asphalt that can fully heal damage within a short period is desired; this has attracted extensive research interest in better mechanical performance and durability for asphalt pavements within civil engineering.
Novel methods to improve the self-healing capability of asphalt mainly include three approaches – inducing heating, encapsulation healing, and incorporating nanomaterials – which can be applied individually or simultaneously5,6. Inducing heating can significantly improve the mobility of asphalt and activate its self-healing for recovery7. The self-healing technology of asphalt by inducing heating can be ascribed to the assisted self-healing technique, which indicates that the self-healing properties of asphalt are improved by external stimuli. The objective of adding the steel wool fibers is to enhance the electrical conductivity so as to increase the healing capacity of the asphalt binder8. The approach to induce heat is to expose these electrically conductive fibers to the high-frequency alternating electromagnetic field, which can induce eddy currents, and the heat energy can diffuse into the asphalt binder by the conductive fibers9. The steel wool fibers enhance not just the electrical conductivity but also the thermal conductivity, both of which can positively affect the self-healing properties of asphalt. However, it is challenging to select the proper mixing time for fibers10. The length of fibers decreases with increased mixing time and influences thermal conductivity, while the decreased mixing time leads to clusters of fibers and impedes the mechanical properties of asphalt9. The encapsulation method can supply light components of aged asphalt such as aromatics and saturates and refresh the self-healing capability of asphalt11,12. However, this is a once-only treatment, and the healing materials cannot be replenished after the release. With the development of nanotechnology, nanomaterials have become promising modifiers for enhancing asphalt-based materials. Asphalt binders incorporated with nanomaterials present better thermal conductivity and mechanical properties13. Graphene with excellent mechanical performance and high thermal performance is regarded as an excellent candidate to improve the self-healing ability of asphalt14,15,16,17. The increased healing properties of graphene-modified asphalt can be attributed to the fact that graphene increases the capacity of the asphalt binder to be heated and produce heat transfer inside the asphalt binder, which means that graphene-modified asphalt can be heated more rapidly and reach up to higher temperature than pure asphalt18. The generated heat can be transferred throughout the graphene-modified asphalt at a faster speed than that through pure asphalt. The crack region of the asphalt binder can be influenced easily and healed faster by the heat flow with higher temperature and higher heating capacity. The self-healing reaction will begin if the energy that is equal to or larger than the healing activation energy exists at the crack surface of the asphalt19. Graphene can improve the thermal activation healing performance and accelerate the healing rate of asphalt19,20. Besides, graphene can save heating energy up to 50% during the healing process, which can benefit energy efficiency and reduce maintenance costs21. As a microwave-absorbent material, graphene is reported to improve the healing ability of asphalt during the rest period of microwave heating22. It is expected that the addition of graphene into asphalt will improve not only the mechanical performance but also the self-healing and energy-saving capacity, which requires in-depth knowledge of the self-healing mechanism.
Self-healing at the nanoscale is mainly due to the wetting and diffusion of asphalt molecules at the fractured faces23. As asphalt consists of various polar and non-polar molecules, its self-healing capability is strongly related to molecular interactions and movements of asphalt molecules of different components1. However, current research mainly relies on experimental techniques to quantify macroscopic mechanical properties, which causes missing information in the change of microstructures and the interactions among asphalt molecules when trying to understand the healing mechanism. The reinforcing mechanism of graphene in the self-healing capability of asphalt is also unclear at this stage. Molecular dynamics (MD) simulations play an influential role in investigating molecular interactions and motions of nanocomposite systems, and link microstructural deformation with molecular interactions and movements24,25,26,27,28,29,30,31. MD simulations have become more and more popular for analyzing material behaviors that cannot be accessed easily by experiments32,33. Existing studies have shown the feasibility and availability of MD simulations in asphalt systems; the cohesion, adhesion, aging, and thermomechanical properties of asphalt and asphalt composites can be explored by MD simulations34,35,36,37. The self-healing behaviors of asphalt can also be predicted by MD simulations38,39,40. Therefore, it is believed that the investigation using MD simulations is an effective way to understand both the self-healing and reinforcing mechanisms.
The objectives of this study are to investigate the self-healing behaviors of pure asphalt and graphene-modified asphalt nanocomposites and to understand the role of graphene in improving the healing capacity of asphalt through MD simulations. The self-healing simulations of pure asphalt and graphene-modified asphalt composites are carried out by introducing cracks into the initial structures. The self-healing capabilities are characterized by the contour of atom numbers, the reorientation and entanglement of molecules at the fractured face, and the mobility of asphalt components during the self-healing processes. By investigating the healing efficiency of graphene at different sites, the reinforcing mechanism of graphene contributing to the self-healing abilities of asphalt is unveiled, which can help with the monitoring of nanofillers in an optimal way and thus enable the life extension of asphalt pavements. An investigation of the self-healing capacity at the atomistic scale can provide an efficient way to develop advanced asphalt-based materials for future research.
According to asphalt chemistry, asphalt consists of various types of hydrocarbons and non-hydrocarbons with different polarity and shapes, which can mainly be divided into the four components of asphaltene, polar aromatics, naphthene aromatics, and saturates41,42. Asphaltene molecules are relatively larger and heavier than other molecules in asphalt, with a mean atomic mass of roughly 750 g/mol and a molecular diameter in the range of 10-20 Å. It has been widely accepted that asphaltene is composed of large aromatic cores that contain heteroatoms and are surrounded by different lengths of alkyl groups43. A modified asphaltene molecule is constructed, as shown in Figure 1a. The molecules of polar aromatics and naphthene aromatics are constructed based on the polarity and the element ratio of asphalt molecules, with benzobisbenzothiophene (C18H10S2) representing the polar aromatic molecule and 1,7-dimethylnaphthalene (C12H12) chosen as the representative naphthene aromatic molecule, as shown in Figure 1b–c. N-docosane (n-C22H46) is constructed as shown in Figure 1d. The parameters listed in Table 1 for asphalt molecules are selected and used to meet the desired criteria, including the elemental mass fraction, the atom ratio, and the aromatic/aliphatic ratio, of real asphalt from experiments41. The same mass ratio has been defined in our previous studies, and the other thermomechanical properties like density, glass transition temperature, and viscosity are in good agreement with experimental data of real asphalt36. The molecular structure of graphene applied in this study is shown in Figure 1e. The adopted graphene sheet in this study has no defect and no fold compared to that of the real case, while the real graphene sheet usually has several defects such as atomic vacancies and Stone-Wales defects44, and some of the graphene sheets can be folded during the mixing process in the asphalt matrix45. These imperfect situations are not considered in this study, since we focus on the effect of the site of the graphene sheet on the self-healing properties and choose it as the only variable. The variables of graphene sheets in terms of the defects and folded cases will be the focus of our future studies. The mass ratio of graphene to asphalt in this study is 4.75%, which is the normal situation (<5%) for graphene modified asphalt in the experiment46,47.
Figure 1: Chemical structure. The atomistic models of (a) asphaltene molecule (C53H55NOS), (b) naphthene aromatic molecule (C12H12), (c) polar aromatic molecule (C18H10S2), (d) saturate molecule (C22H46), (e) graphene, and (f) pure asphalt. For the atomistic asphalt model, the carbon, oxygen, nitrogen, sulfur, and hydrogen atoms are shown in gray, red, blue, yellow, and white, respectively. Please click here to view a larger version of this figure.
Asphalt model | Mass (g/mol) | Chemical formula | Numbers of molecules | Total mass (g/mol) | Mass fraction (%) |
Asphaltene | 754.04 | C53H55NOS | 43 | 32423.72 | 26 |
Naphthene aromatic | 156.22 | C12H12 | 65 | 10154.3 | 8 |
Polar aromatic | 290.38 | C18H10S2 | 74 | 21485.16 | 17 |
Saturate | 310.59 | C22H46 | 205 | 63670.95 | 49 |
Asphalt binder | 387 | 127734.13 | 100 | ||
Graphene | 6369.28 | C525H63 | 1 | 6369.28 |
Table 1: Overall components of pure asphalt model and graphene-modified asphalt model.
With respect to the protocol described below, two types of wedge-like cracks with different sizes are inserted into the middle of the asphalt model with a blunt crack tip and two parallel crack surfaces, while the middle-top area of the asphalt bulk remains intact. Two crack widths are chosen as 15 Å and 35 Å, as shown in Figure 2a–b. The reason for selecting 15 Å is that the crack width should be wider than the cutoff of 12 Å to avoid the early self-healing of asphalt molecules during the equilibrium process while investigating an extreme case for a small crack. The reason for selecting 35 Å is that the crack width should be wider than the length of the saturate molecules of 34 Å in order to prevent the bridging effect. The height of the crack is 35 Å, the same as the box width, and the depth of the crack is 70 Å, the same as the box length. In the real situation, the observed micro-crack sizes can be varied in the range from several micrometers to several millimeters, which is far larger than the length scale we are modeling here. Normally, the length scale in MD simulation is limited to the scale of 100 nm, which is still several orders of magnitude smaller than the real crack size. However, the cracks initiate at the nanoscale and grow into macroscale cracks with continuous deformation48. The understanding of the self-healing mechanism at the nanoscale can help to prevent the growth and further propagation of the crack at the macroscale. Even though the selected crack sizes are in the range of nanometer, the results can still be influential and applicable to explore the self-healing behaviors of asphalt molecules. There are two locations for the graphene sheets in the crack areas: one is on top of the crack tip and the other is perpendicular to the left crack surface. It has been found that these are the most common positions for graphene in graphene-modified nanocomposites with cracks49.
Figure 2: The self-healing schemes for pure asphalt and graphene-modified asphalt. The self-healing model of pure asphalt with a crack width of (a) 15 Å and (b) 35 Å. The self-healing model of graphene-modified asphalt with the graphene sheet is located (c) at the top of the crack tip and (d) perpendicular to the crack surface. Please click here to view a larger version of this figure.
In MD simulations, the intramolecular and intermolecular interactions in the asphalt nanocomposites are described by the Consistent Valence Forcefield (CVFF)50, which works well with asphalt and graphene-based materials. The functional form of CVFF is expressed as the following expression:
1
Here, the total energy Etotal is composed of the bonded energy terms and the non-bonded energy terms. The bonded interactions consist of the covalent bond stretching, the bond angle bending energy, the torsion angle rotation, and the improper energies as expressed in the first four terms. The non-bonded energy includes an LJ-12-6 function for the van der Waals (vdW) term and a Coulombic function for the electrostatic interactions. CVFF has been widely employed in simulating asphalt materials51,52. The simulated physical and mechanical properties such as density, viscosity, and bulk modulus are in good agreement with the experimental data, which demonstrates the reliability of CVFF51. CVFF is not only suitable for inorganic materials, but it has also been successfully employed in structures consisting of organic and inorganic phases such as asphalt-silica52 and the system of epoxy-graphene53. In addition, the interfacial interactions between graphene and asphalt can be characterized by CVFF36,54. Since the major part in selecting forcefield is to determine the asphalt-graphene interface, the non-bonded interactions described by CVFF are more reliable, which is also considered in our previous study36. Overall, the forcefield CVFF is adopted in this study. The partial charges for different kinds of atoms are calculated by the forcefield-assigned method.
1. Build the atomistic models
2. Perform the simulations
3. Postprocessing
The contour of atom number
The contours of the atom number of pure asphalt and graphene-modified asphalt models in the yz plane are shown in Figure 3, where the color bar from blue to red exhibits atom numbers varying from 0 to 28. Figure 3a–c illustrates the contour of the atom number of the structures with 15 Å crack width in pure asphalt and asphalt nanocomposites modified by graphene at the crack tip and the crack surface. For pure asphalt, complete healing occurs after about 300 ps. The self-healing behavior starts from the area of the crack tip, as the area around the crack tip becomes a melted blunt shape with a blue color after 50 ps, and several asphalt molecules bridge the two crack surfaces at the middle of the crack tip. The green color in the contour presents the bulk asphalt, which is the stage of the crack zone getting full self-healing. At around 100 ps, the crack zone is almost closed with a small void left, and the colors of the initial crack surfaces change to green, which indicates that the self-healing process is finished in these areas; however, there are still some blue and white areas remaining to be self-healed. After about 300 ps, most of the color of the crack zone has changed to green, which is the same as that of the asphalt bulk, indicating that the self-healing process is complete. As shown in Figure 3b, the self-healing process is not significantly changed after adding the graphene sheet on the top of the crack. The self-healing process takes about 500 ps to complete, and the crack zone is sharply decreased at 50 ps and almost disappears at 200 ps. The graphene sheet on the top of the crack tip seems to have little influence on the self-healing process of the crack surface. However, inserting the graphene on the left of the crack surface can significantly accelerate the self-healing process, as shown in Figure 3c, where the red line in the contour is the graphene sheet. The self-healing period is shortened to around 200 ps, which is half as much as that of pure asphalt. The crack width is significantly decreased at 20 ps, and the asphalt molecules from the bulk tend to move to the graphene area and fill up the crack area. The crack zone almost disappears at around 150 ps, although some of the areas at the bottom remain blue. After another 50 ps of the self-healing process, the crack area is full of blue color, which indicates the end of the process.
The self-healing process of the models with 35 Å crack width takes almost twice as long as that of the models with 15 Å crack width, while the self-healing process of pure asphalt lasts around 1,000 ps. The self-healing behavior starts in the crack tip area, and the crack shape becomes shrunken and irregular at 100 ps. Most of the crack zone is healed by 500 ps, with a small void left in the middle of the crack zone. After performing the self-healing process for another 500 ps, the crack zone is filled with asphalt molecules until the self-healing process is completed. The graphene sheet is located on the top of the crack tip, as shown by the red line of Figure 3e. The self-healing period is about 1,100 ps, which is close to that of pure asphalt. However, the crack shape changes differently. There are some asphalt molecules that bridge the crack area at around 400 ps, which can advance the self-healing process. As shown in Figure 3f, the self-healing behaviors can be significantly improved when the graphene sheet is located at the left crack surface. A phenomenon can be observed similar to the model with 15 Å crack width: some of the asphalt molecules in the asphalt bulk tend to move to the graphene area and wrap around the graphene sheet, which can decrease the crack area significantly and help the self-healing process. The width of the crack is decreased to around half of the initial crack width by only 50 ps, and most of the crack area is healed at around 300 ps. The whole self-healing process lasts about 600 ps and most of the crack zone disappears; this takes only half the length of time taken by pure asphalt.
Figure 3: The contour of the atom number during the self-healing process. The contour of the atom number during the self-healing process for the models with 15 Å crack width for (a) pure asphalt, (b) graphene on the crack tip, and (c) graphene at the left surface of the crack tip, and the models with 35 Å crack width for (d) pure asphalt, (e) graphene on the crack tip, and (f) graphene at the left surface of the crack tip. The dotted black boxes refer to the locations of graphene. The color bar from blue to red stands for the atom numbers varying from 0 to 28 in the contour. Please click here to view a larger version of this figure.
Molecular interactions
To explore the difference in self-healing behaviors between pure asphalt and graphene-modified asphalt composites, the molecular interactions and movement during the self-healing process are captured and analyzed, as shown in Figure 4. From Figure 4a, it can be observed that aromatic molecules like asphaltene, polar aromatics and naphthene aromatics are attracted by the graphene sheet through π-π stacking when graphene is placed on the top area of the crack tip. These asphalt molecules are captured tightly by the graphene sheet and cannot readily diffuse into the neighborhood of the crack zone or fill up the crack, which hinders the self-healing process to a certain extent. However, the self-healing behaviors mainly arise from the asphalt molecules near the crack surface, and the influence of these molecules in the top area needs further exploration. From Figure 4b, it is observed that the polar aromatic molecule at the crack surface is attracted by the graphene sheet at the other crack surface, which can further increase the nearby naphthene aromatic molecule's likelihood of moving onto the crack area. The gathered asphalt molecules attracted by the graphene sheet can pack the crack zone with a higher speed than that of pure asphalt and the self-healing ability of the graphene-modified asphalt nanocomposite can be improved. The self-healing process of the model with 35 Å crack width modified by graphene at the left crack surface is shown in Figure 4c. The polar aromatic molecules are attracted by the graphene sheet through π-π stacking when the self-healing begins, and these asphalt molecules can rapidly wrap around the graphene sheet and reduce the space of the crack zone, as shown in Figure 3f. This indicates that graphene plays an important role in the initial stage of self-healing when it is located around the crack surface. A self-healing snapshot of pure asphalt with 15 Å crack width is shown in Figure 4d. It can be clearly observed that the chain structure of saturate is important for the self-healing process since the molecules can become entangled with each other and bridge the crack surface. This bridging effect between saturate molecules and the side chains of asphaltene molecules can significantly increase the packing efficiency and decrease the self-healing period. It is also observed that asphalt molecules with polyaromatic rings, such as asphaltene, polar aromatics, and naphthene aromatics, reorientate themselves at the crack surface by π-π stacking. This reorientation allows the asphalt molecules to move in a parallel direction and contributes to crack wetting, which further closes the crack surfaces.
Figure 4: Details of non-bond interaction of pure asphalt and graphene-modified asphalt nanocomposites during the self-healing process. For the model with (a) 15 Å crack width and graphene located on the top of the crack tip, aromatic molecules in asphalt are attracted by the graphene sheet through π-π stacking. For the model with (b) 15 Å crack width and graphene at the left side of the crack surface, the polar aromatic molecules at the other crack surface move onto the graphene surface due to strong aromatic interactions. For the model with (c) 35 Å crack width and graphene at the left side of the crack surface, the polar aromatic molecules are attracted by the graphene sheet and thus protrude from the crack surface. For the model with (d) 15 Å crack width and pure asphalt, there is a reorientation of aromatic molecules at the crack surface and a chain bridging and entanglement of saturate molecules during the self-healing process. The blue dotted boxes and purple dotted boxes in the figure indicate the π-π stacking and reorientation behaviors, respectively. Please click here to view a larger version of this figure.
The reorientation of aromatic molecules including asphaltene, polar aromatics, and naphthene aromatics around the crack surface during the self-healing process, is shown in Figure 5. Figure 5a shows that the tracked molecules before self-healing are almost perpendicular between naphthene aromatics and polar aromatics and between asphaltene and polar aromatics. The distance between asphaltene and the other two aromatic molecules is 13.3 Å, which is larger than the distance between the aromatic molecules. After self-healing for 40 ps, the naphthene aromatic molecule diffuses to the space between asphaltene and polar aromatics and plays an important role in interacting with the other two molecules. In Figure 5b, it can be observed that the distance and angle between the polar aromatic molecule and naphthene aromatic molecule are 4.6 Å and 89°, which indicates a T-shaped π-π stacking interaction between the two aromatic molecules. The angle and distance between the naphthene aromatic and asphaltene decrease to 32° and 4.6 Å, respectively. This indicates that the non-bond interactions between naphthene aromatics and asphaltene cause them to rotate and adjust the orientation gradually, contributing to the wetting of the crack surface. The orientations of the three molecules are almost parallel after 50 ps, as the angles between them are 26° and 35°, as shown in Figure 5c. The distance between them declines below 4.0 Å, which indicates that the π-π stacking facilitates the parallel structure and brings the aromatic molecules closer together. Overall, the reorientation at the crack surface promotes the interaction of asphalt molecules, which shortens the intermolecular distance and increases the attraction among them. The reorientation and diffusion of asphalt molecules further aid the filling up of the crack zone and speed up the self-healing process.
Figure 5: The reorientation of asphalt molecules during the self-healing process. Angles and distances among asphalt molecules (a) before self-healing, (b) after 40 ps, and (c) at 50 ps. Please click here to view a larger version of this figure.
Mobility of asphalt components
To understand the roles of different components quantitatively in the self-healing behaviors of asphalt, the MSD of the center mass for asphalt composite is calculated to represent the transitional mobility during the self-healing process, which is expressed by:
where ri(t) is the position vector of particle i at time t, and the angular bracket indicates the average value of the traveling distance. The MSD values of pure asphalt and graphene-modified asphalt are tracked and shown in Figure 6. Figure 6a–c shows the MSD of asphalt composite with a 15 Å width crack, while those with a 35 Å width crack are shown in Figure 6d–f. It can be observed that saturates are the most active component in the self-healing behaviors of asphalt, while asphaltene is the least active. There are two possible reasons: one relates to the molecular mass, as asphaltenes have the highest molecular mass in asphalt, and they are less able to move and fill the crack zone. The other is the chain-like structure of saturates, which have higher mobility than other components and are more likely to become entangled and stretch out at the crack surface. The mobility of polar aromatics is higher than that of naphthene aromatics; this is because of the higher molecular mass and polarity of polar aromatics. The polar atoms on polar aromatics, such as sulfur atoms, can form an H-bond with asphaltenes, and the mobility can be impeded. The MSD figures for graphene-modified asphalt on the top area of the crack tip and the left crack surface are shown in Figure 6b and Figure 6c. It can be seen in Figure 6b that the MSD of graphene is lower than those of asphalt components since graphene occupies the largest volume and has the highest molecular mass in the asphalt nanocomposites. The MSD values of asphalt components are relatively lower than those of pure asphalt; this is because the interactions among these molecules and graphene hinder the mobility of the asphalt molecules and slow the self-healing process. However, when graphene is placed at the left crack surface, the mobilities of polar aromatics, naphthene aromatics, and graphene improve significantly compared to that of pure asphalt. This indicates that graphene plays an important role in the self-healing process and that its interactions with aromatic molecules in asphalt contribute to the self-healing process of asphalt. For the 35 Å width crack case in Figure 6d, the MSD of pure asphalt follows a similar trend to that of the case with 15 Å crack width, as the MSD of asphaltene, polar aromatics, naphthene aromatics, and saturates vary in an increasing manner. When inserting graphene into the top area of the crack tip, the MSD of saturates decreases by about 15 Å2. The existence of graphene sheets in asphalt bulk vertically influences the mobile space of the saturate molecules and prevents the routes of self-healing. From Figure 6f, it can be observed that the MSD values of asphaltene, polar aromatics and naphthene aromatics are all improved compared to pure asphalt, while the MSD of saturates decreases slightly. Graphene is highly responsible for the improvement of the self-healing process, especially with molecules containing aromatics. The π-π stacking interactions among graphene and asphaltene, polar aromatics, and naphthene aromatics improve the mobility of these asphalt molecules and help to form a stable packing structure in the crack zone, which accelerates the self-healing process of asphalt.
Figure 6: MSD of pure asphalt and graphene-modified asphalt molecules during the self-healing process. For the models with 15 Å crack width, MSD of (a) pure asphalt and graphene-modified asphalt on (b) the top of the crack tip and (c) the left surface is presented. For the models with 35 Å crack width, MSD of (d) pure asphalt and graphene-modified asphalt on (e) the top of the crack tip and (f) the left crack surface is presented. The X-axis represents the time of the simulation and the Y-axis represents the MSD values of asphalt components and the graphene molecule during the self-healing process. Please click here to view a larger version of this figure.
Molecular locations after self-healing
To explore the relative locations between graphene and asphalt molecules during the self-healing process, the radial distribution functions between graphene and aromatic molecules in asphalt are calculated and shown in Figure 7. Figure 7a–c shows the RDF of the model with 15 Å crack width before and after the self-healing process. It can be seen that the aromatic molecules in asphalt move closer to the graphene sheet following the self-healing process, especially the polar aromatic molecules and naphthene aromatic molecules. As indicated in Figure 4, there are strong π-π stacking interactions between graphene and aromatic molecules such as asphaltene, polar aromatics, and naphthene aromatics, which cause the graphene sheet to attract these molecules toward the crack surface. However, the difference in the g(r) values of asphaltene before and after self-healing is not as significant as those of polar aromatics and naphthene aromatics. This is because the asphaltene molecules obtain a higher molecular mass and volume than the polar aromatic and naphthene aromatic molecules, making it harder for them to rotate and diffuse to the graphene area and fill the crack zone. The increased g(r) values between graphene and polar aromatic or naphthene aromatic molecules within 4.0 Å are within the typical interaction distance for π-π stacking, and the increased g(r) values beyond 4.0 Å are due to the combination of molecular interactions and the elimination of the crack zone. An RDF of the model with 35 Å crack width before and after the self-healing process is shown in Figure 7d–f. The g(r) values between graphene and asphaltene beyond 4.0 Å through the self-healing process are more obvious than those of the 15 Å crack width; this is because asphaltene has more space to diffuse and move toward the graphene in the larger crack zone. The g(r) values within 4.0 Å are more significant for naphthene aromatics than those for polar aromatics; this is due to the smaller molecular mass and better diffusion ability of naphthene aromatic molecules.
Figure 7: The RDF values between graphene at the left crack surface and the asphalt components. The RDF values between graphene at the left crack surface and the asphalt components of (a) asphaltene, (b) polar aromatics, (c) naphthene aromatics from the model with 15 Å crack width, (d) asphaltene, (e) polar aromatics, and (f) naphthene aromatics from the models with 35 Å crack width. The X-axis represents the distance of the two molecules and the Y-axis represents the RDF values. Please click here to view a larger version of this figure.
The critical steps within the Protocol part are as follows: step 1.4 – Build and pack the four types of asphalt molecules; step 1.5 – Build the asphalt structure with the crack; step 2.3 – Achieve the equilibrium; step 2.4 – Perform the self-healing process. These steps indicate the most cohesive and important contents of the protocol. To create the desired shapes of the inserted crack, the packing process is modified compared to the normal packing in Materials Studio. The crack shape is created and filled inside the simulation box, and then the asphalt molecules are packed into the other part of the simulation box. After that, the redundant asphalt molecules are deleted around the created crack contour. The limitation of MD simulations is that the time scale and length scale are relatively small in the order of nanosecond and nanometer compared to traditional methods such as the finite element method, in which the simulations can be analyzed up to seconds and meters57. The significance of this method is that it can reveal the self-healing mechanism of asphalt and graphene-modified asphalt at the atomistic level by capturing nanostructure evolution, molecular interactions, and motions, which are hard to be accessed by traditional approaches58. The self-healing mechanism can help researchers and engineers to apply nanomaterials at the appropriate site and improve asphalt in an efficient way. The future application of this technique is that it can monitor the molecular structure in a decent manner and help to investigate the effect of other variables of nanomaterials such as defects, folded structures, and functional groups. This technique can also be combined with other approaches to observe the self-healing behaviors of asphalt nanocomposites from a multiscale aspect. The self-healing properties of asphalt can be understood thoroughly and get improved significantly in the future.
Graphene is critical in the changes and migration of the interface and components during the self-healing process. Without inserting the graphene sheet, saturate plays an important role in the self-healing process, since the chain structure of saturates can entangle with each other and bridge the crack surface. The bridging effect between saturate molecules and the side chains of asphaltene molecules can strongly increase the packing density and decrease the time of the self-healing process. Besides, asphalt molecules with polyaromatic rings, such as asphaltene, polar aromatics, and naphthene aromatics, reorientate themselves at the crack surface by π-π stacking, which makes the asphalt molecules move in a parallel direction and contributes to crack wetting and closes the crack surfaces. With the insertion of graphene, the polar aromatic molecules at one side of the crack surface are attracted by the graphene sheet at the other side of the crack surface, which can further increase the possibility of nearby naphthene aromatic molecules moving into the crack area. The gathered asphalt molecules attracted by the graphene sheet can fill in the crack zone with a higher speed than that in pure asphalt, and the self-healing ability is significantly improved in the graphene-modified asphalt nanocomposite. Asphaltene molecules have higher molecular mass and volume in the asphalt matrix, which makes them hard to diffuse to the graphene part and fill in the crack zone. Naphthene aromatics have faster movement than polar aromatics, which is due to the smaller molecular mass and better diffusion ability of naphthene aromatic molecules39.
In this study, the self-healing properties of pure asphalt and graphene-modified asphalt nanocomposites are investigated in consideration of different crack widths and graphene locations using MD simulations. It is observed that the self-healing behavior starts from the crack tip area, with the sharp tip becoming blunt and fuzzy. Asphalt molecules at the crack boundary can diffuse to decrease the width of the crack and continue filling the gap. The complete self-healing process is confirmed when the atom density of the crack area is the same as that of the asphalt bulk. MD simulations can help to reveal the molecular interactions and the chain movement in the asphalt matrix during the self-healing process. The entanglement and reorientation of asphalt molecules play an important role in self-healing behaviors. The self-healing rate with the incorporation of the graphene sheet is determined by its location. For the graphene sheet located at the crack tip area, the movement of the asphalt molecules is hindered and cannot easily diffuse into the crack zone. For the graphene sheet at the side of the crack zone, the asphalt molecules are attracted by the graphene sheet due to π-π stacking interaction and easily gather at the crack zone, indicating an increasing self-healing rate. The simulation results show that the modification of asphalt by nanomaterials can improve both thermomechanical and self-healing properties, which has great potential for the development of smart asphalt pavements. The fundamental understanding of the self-healing mechanism in asphalt nanocomposites based on MD simulations can facilitate efficient manipulation of nanomaterials at the optimal site, which is beneficial to the advanced design of asphalt nanocomposites with desired properties and functions.
The authors have nothing to disclose.
The authors are grateful for the support from City University of Hong Kong Strategic Research Grant with the Project No. 7005547, the support from the Research Grants Council (RGC) of the Hong Kong Special Administrative Region, China, with the Project No. R5007-18, and the support from Shenzhen Science and Technology Innovation Committee under the grant JCYJ20170818103206501.
Atomistic models of asphalt and graphene/Materials Studio | BIOVIA | Materials Studio 8.0 | The atomistic models are built for molecular dynamics simulations. |
Large-scale Atomic/Molecular Massively Parallel Simulator Package | Sandia National Laboratories | lammps-stable20 | The equilibrium is achieved under NPT ensemble, and the atomistic models get self-healed. |
OVITO | Materials Science Department of Technische Universität Darmstadt, Germany | ovito-basic-3.1.0-win64 | The self-healing behaviors of the atomistic models are visualized. |
Origin | OriginLab | Origin 2018 64Bit | The contours of the atom numbers of the trajectory are drawn and analyzed. |