Next Article in Journal
Enhancement in Crystallizability of Poly(L-Lactide) Using Stereocomplex-Polylactide Powder as a Nucleating Agent
Next Article in Special Issue
Study on the Influence of Nano-OvPOSS on the Compatibility, Molecular Structure, and Properties of SBS Modified Asphalt by Molecular Dynamics Simulation
Previous Article in Journal
Synthesis and Characterization of Quadrupolar-Hydrogen-Bonded Polymeric Ionic Liquids for Potential Self-Healing Electrolytes
Previous Article in Special Issue
Comprehensive Study on the Performance of Waste HDPE and LDPE Modified Asphalt Binders for Construction of Asphalt Pavements Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of the Microscopic Mechanism of Natural Rubber (Cis-1, 4-Polyisoprene, NR)/Polyethylene (PE) Modified Asphalt from the Perspective of Simulation

1
School of Highway, Chang’an University, Xi’an 710064, China
2
College of Civil Engineering, Henan University of Technology, Zhengzhou 450001, China
3
Department of Structural Engineering, College of Civil Engineering, Tongji University, Shanghai 200092, China
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(19), 4087; https://doi.org/10.3390/polym14194087
Submission received: 28 August 2022 / Revised: 16 September 2022 / Accepted: 21 September 2022 / Published: 29 September 2022
(This article belongs to the Special Issue Application of Polymer Materials in Pavement Design)

Abstract

:
This paper aims to study the interaction mechanism of waste tire/plastic modified asphalt from the microscopic perspective of molecules. Based on BIOVIA Materials Studio, a classic four-component asphalt model consisting of asphaltene (C149H177N3O2S2), resin (C59H85NOS), aromatic (C46H50S), and saturate (C22H46) was constructed. Waste tires are represented by natural rubber (NR), which uses cis-1, 4-polyisoprene as a repeating unit. In contrast, waste plastics are characterized by polyethylene (PE), whose optimum degree of polymerization is determined by the difference in solubility parameters. Then, the above molecular models are changed to a stable equilibrium state through the molecular dynamics process. Finally, the interaction process is analyzed and inferred using the indexes of radial distribution function, diffusion coefficient, and concentration distribution; further, the interaction mechanism is revealed. The results show that the optimal degree of polymerization of PE is 12, so the solubility parameter between PE and NR-modified asphalt is the lowest at 0.14 (J/cm3) 1/2. These models are in agreement with the characteristics of amorphous materials with the structures ordered in the short-range and long-range disordered. For NR-modified asphalt, the saturate moves fastest, and its diffusion coefficient reaches 0.0201, followed by that of the aromatic (0.0039). However, the molecule of NR ranks the slowest in the NR-modified asphalt. After the addition of PE, the diffusion coefficient of resin increased most significantly from 0.0020 to 0.0127. NR, PE, and asphaltene have a particular attraction with the lightweight components, thus changing to a more stable spatial structure. Therefore, using NR and PE-modified asphalt can change the interaction between asphalt molecules to form a more stable system. This method not only reduces the large waste disposal task but also provides a reference for the application of polymer materials in modified asphalt.

Graphical Abstract

1. Introduction

Petroleum asphalt is an essential raw material in road engineering, especially in asphalt pavements, which are used widely. At present, petroleum is gradually becoming exhausted, so asphalt, a by-product of petroleum, will also face the problem of exhaustion. Against this background, it is indispensable and crucial to lower the consumption of asphalt in road engineering and decrease dependence on petroleum. In addition, the deteriorating environment and traffic conditions put higher requirements on the performance of asphalt pavements. Therefore, the modification of asphalt to reduce cost and improve performance is one of the leading research directions for researchers, which is very significant. The use of modified asphalt has recently increased, and its application scope and applications are becoming broader [1,2,3,4,5].
There are many kinds of materials that can be used as modifiers, and the current research mainly includes the following aspects. Some researchers use biomaterials as modifiers to prepare bio-asphalt [6,7,8], such as beans [9], waste cooking oil [10,11], biochar [12,13], and so on. Forestry resources and animal excreta are liquefied and separated to prepare modifiers. In addition, in the process of biomass raw materials processing to produce bio-based materials (furfural, furfuryl alcohol, dimethylfuran, biodiesel, etc.), the remaining polymer carbohydrates can also be used to create bio-asphalt. Some researchers use soil as a modifier, such as nano clay [14], diatomite [15,16], organically modified bentonite [17,18] and so on. Take bentonite as an example. Its electron microscope diagram is shown in Figure 1, where (a) represents the original bentonite and (b) represents the nano-bentonite, which shows a clear layered structure and good continuity. Bentonite itself belongs to the octahedral structure. If added to the asphalt, it can fill the pores between the micelle and micelle and further increase the anti-skid resistance of asphalt. Fibers are also adopted as modifiers, such as polyester [19,20], lignin [21,22], basalt fibers [23,24] and so on [25]. This kind of modifier can enhance the cohesion between aggregate particles so that the asphalt pavement can adapt to high wear requirements, increased skid resistance, and other sections. Guo [26,27], Zheng [28], and some other researchers [29,30] used waste tires as raw materials to prepare rubber-modified asphalt and conducted research. Some researchers use plastics as modifiers, such as polyethylene, polyurethane, polypropylene, etc. [4,31,32]. Some nano-modifiers are carbon nanotubes, graphene, and other nanomaterials, but the cost is high.
It should be noted that as asphalt modifiers, waste tires and plastics have a significant negative influence on the ecological environment [33]. Waste tires and plastics are non-biodegradable substances, which brings difficulties to waste management [27]. If incinerated, a lot of smoke and other harmful substances will be produced, which will significantly impact the air, soil, rivers, and other environments. However, waste plastics and rubber have good reuse value as asphalt modifiers and can improve asphalt performance [28]. Generally speaking, plastic polymers can significantly enhance asphalt performance at high temperatures. In contrast, rubber polymers can upgrade the deformation resistance of asphalt at low temperatures, such as styrene-butadiene rubber (SBR). Therefore, the use of the two as asphalt modifiers simultaneously becomes one of the best choices for the recycling and reuse of waste material. At present, only one kind of modified agent is doped separately. Still, there are few studies on adding two types of waste materials simultaneously. Hence, it is not easy to give full play to the synergistic effect of the two and improve the high- and low-temperature performance of asphalt simultaneously.
It is a much more cost-effective and scientific to perform research using molecular dynamics. The traditional test has a long test cycle and a high cost. In addition, many conventional experiments only describe and compare the experimental results, which fails to explain some empirical phenomena in depth. Molecular dynamics, on the other hand, refers to the study of microscopic phenomena from the perspective of molecules or atoms. Currently, many researchers are studying the molecular dynamics of asphalt, such as the modification mechanism of modifiers [34]. For example, Hu et al. [35] studied the interaction mechanism between styrene-butadiene-styrene (SBS) and asphalt to explain the formation and separation mechanism of an SBS-modified asphalt two-phase structure from a microscopic perspective.
The purpose of this paper is to study the molecular interaction mechanism after two modifiers simultaneously using molecular dynamics. In this paper, the original asphalt and rubber-modified asphalt models were built from the perspective of molecular dynamics. Then, plastic was added to create a plastic/rubber asphalt model, and the influence of plastic molecules on rubber asphalt was analyzed. Among the models, the classical four-component molecular model was taken to represent asphalt, and rubber is expressed by natural rubber molecules composed of polyisoprene. Additionally, plastics are represented by polyethylene molecules, the optimal polymerization degree of which is measured by solubility parameters. After the molecular dynamics process, all the above models reach a state of equilibrium and stability. The radial distribution function (RDF) was used to verify the molecular model, and the mean square displacement (MSD) and diffusion coefficient were used to study the diffusion of molecules in modified asphalt. Finally, concentration distribution was taken to study the spatial distribution of different molecules.

2. Simulation

2.1. Model of Asphalt

Bitumen is a highly complex mixture of hydrocarbons of different molecular weights and their derivatives of non-metallic components (oxygen, nitrogen, sulfur, etc.). Petroleum bitumen, a vital type of bitumen, is widely used in road engineering. From a scientific research perspective, the parts of petroleum asphalt with some common characteristics, such as similar chemical composition and physical properties, are divided into the same group, called components. The American Society for Testing and Materials (ASTM) D4124-09 has proposed and recommended a four-component analysis method. This method divides petroleum asphalt into asphaltene, aromatic, resin, and saturate.
The simplified model adopted in this study consists of four representative molecules, and this method has been widely used because of its simplicity, convenience, and conformity with the actual situation [36]. The four representative molecules are asphaltene (C149H177N3O2S2), resin (C59H85NOS), aromatic (C46H50S), and saturate (C22H46), as shown in Figure 2. The proportion of the four representative molecules in this model is calculated according to the mass ratio and molecular weight of each component in asphalt, in line with previous studies [36,37], as shown in Table 1.

2.2. Models of Rubber

Rubber powder mainly comes from waste tires, and its composition is complicated. The three most common types of rubber are NR (natural rubber), BR (cis-polybutadiene rubber), and SBR (styrene-butadiene rubber). Since natural rubber is the most widely used, it was selected as the research object in this study.
Natural rubber mainly comprises cis-1, 4-polyisoprene, accounting for more than 97% and 2–3% of the 3- and 4-bonded structure. Therefore, when building the molecular model of rubber, NR was regarded as a homopolymer, which was polymerized from the isoprene mentioned above. The relative molecular weight range of NRs is relatively wide, generally between 30,000 and 30 million, because of the different degrees of polymerization. The degree of polymerization affects the results of the experiments and leads to different simulation costs, so a proper selection is needed. Referring to previous studies [29,38], the degree of polymerization selected in this paper is 16, and the number of repeating units is set at 16. The molecular chain for NR is shown in Figure 3.

2.3. Models of PE

Polyethylene, short for PE, is a thermoplastic resin prepared by the polymerization of ethylene, including ethylene and a small amount of α-olefin copolymer. The degree of polymerization, a key factor, affects the molecular weight and the amount of simulation, and it is also closely related to the solubility parameters. The independent variables of the polymerization degree for polyethylene used in this article are 6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 72, 78, 84, 90, 96, 102, 108, 114, 120, and there are 20 different degrees of polymerization in total. Firstly, the degree of polymerization is closely related to the amount of simulation calculations because the degree of polymerization is a direct factor affecting the molecular weight; the greater the degree of polymerization, and the greater the molecular weight, the larger the amount of simulation calculations. Secondly, the solubility parameters of polyethylene with different degrees of polymerization are different. If the difference in the solubility parameters is less than 1.3~2.1 (J/cm3) 1/2, the two materials are compatible. The smaller the difference is, the more easily the two are compatible. Therefore, in this paper, the solubility parameter difference is selected to evaluate the solubilization effect of PE on rubber-modified asphalt.

2.4. Models of the Mixture

The simulation software used in this study is the BIOVIA Materials Studio, an integrated computer simulation platform involving quantum mechanics, molecular dynamics, and mesoscopic dynamics. Its function is to reveal the interaction mechanism of substances or predict the properties of materials from the microscopic perspective of molecules or atoms; therefore, it is widely used in various research fields.
The basic principles of molecular dynamics simulation were introduced. A system with N particles and the description of classical mechanics depends on 3N degrees of freedom ( q , p )   { q i , p i | i = 1   ( 1 )   N } . The position q { q i | i = 1   ( 1 )   N } , momentum P { p i | i = 1 ( 1 ) N } .   q i is the position of the particle N o .   i and p i is the linear momentum. ( q , p ) determines the system’s state, and each microscopic state is a phase point in phase space.
In molecular dynamics studies, the laws of Newtonian mechanics are usually adopted, as shown in the related literature [39].
d v i d t = F i   ( t )   m i = 1 m i i U       i = 1 ( 1 ) N
d r i   ( t )   d t = v i   ( t )       i = 1 ( 1 ) N
where, r i and v i are the position and velocity of particle N o .   i , respectively. The bond action, similarly to the quantum effect in the above equation, is reflected in the potential energy U . The potential energy U , which depends on the position of all of the nuclei, can be calculated precisely by quantum chemistry or approximated by force-field methods.
Formula (1) and (2) are a set of 3 N differential equations, which should be solved by the numerical integration method of differential equations. Therefore, the position { r i ( t ) } and velocity { v i ( t ) } of the system at each moment can be solved by the integration algorithm. In this case, the trajectory of all the nuclei of the system is given, that is, the time evolution of the system. Then, the change of all system properties with time can be obtained.
COMPASS is a powerful force field that supports the atomic level simulation of condensed matter materials, short for “Condensed-Phase Optimized Molecular Potential for Atomistic Simulation Study.” It is the first ab initio force field parameterized and verified by condensed matter properties and various ab initio and empirical data of isolated molecules. Therefore, the structure, conformation, vibration, and thermophysical properties of different molecules in isolated or condensed systems could be accurately predicted for a wide range of temperatures and pressures by this field [40,41,42,43,44].
The specific steps of building a stable model are as follows.
  • The establishment of initial polymer model.
  • Geometric optimization.
  • The simulation of annealing.
  • In this procedure, the pressure (1.01 × 10−4 Gpa) and temperature of the system remain constant [45].
  • NVT+NPT dynamics simulation.
  • At ambient temperature and standard atmospheric pressure, 50 ps of MD simulation was performed to bring the model closer to its natural state [46,47]. Taking the original asphalt model as an example, after the above steps were performed, parameters such as the energy, density, and others for the system reached a stable state, as shown in Figure 4. For example, the density of the original asphalt model increased and stabilized at about 1.0, which coincided with the actual situation.
  • The calculation and analysis of parameters.
  • After the molecular dynamic simulation, parameters related to this study were analyzed and calculated, such as the solubility parameters and radial distribution function, etc., according to the stable model generated at the end of the simulation process.
The specific composition of the two modified asphalt models, including NR-modified asphalt and PE-modified asphalt, is shown in Table 2 and Table 3. Both the original asphalt model and the modified asphalt model obey the above steps. In the NR-modified asphalt, one NR molecule was included, and its content was 9.28%. The original asphalt model and NR-modified asphalt model in an equilibrium state are shown in Figure 5.
Figure 4. Dynamic equilibrium process.
Figure 4. Dynamic equilibrium process.
Polymers 14 04087 g004
In PE/NR-modified asphalt, the degree of polymerization of the PE molecule should first be determined. In this study, different molecular models were constructed for PE molecules with varying degrees of polymerization, and the kinetic simulation process was performed according to the steps mentioned above. The equilibrium model was finally obtained, as shown in Figure 5.

3. Results and Discussion

3.1. Solubility Parameters

Molecular models of PE with different degrees of polymerization are shown in Figure 6. The solubility parameters were used to evaluate the compatibility between polyethylene with different degrees of polymerization and NR-modified asphalt. If the solubility parameters for two polymer materials become closer, the blending effect is better [42]. In this study, solubility parameters were adopted to evaluate the compatibility between different polymers, primarily to determine the degree of polymerization of PE molecules. Its value is the square root of the cohesive energy density, as shown in the following Formula (3):
δ = C E D = E V
where Δ E represents cohesive energy and V represents volume.
If the action of the hydrogen bond is ignored, the intermolecular interaction force mainly consists of a van der Waals force and electrostatic force; the C E D is composed of the van der Waals cohesive energy density ( C E D v ) and electrostatic cohesive energy density ( C E D e ).
Then, the related cohesive energy index is obtained and used to estimate the solubility parameter. The solubility parameters of the PE additives with individual degrees of polymerization are shown in Table 4, while that of the NR-modified asphalt is 17.23. The difference between the two is shown in Figure 7. According to the principle of similarity-compatibility, when the degree of polymerization is 12, the difference between the two is the smallest, which is 0.14 < 1.3~2.1 (J/cm3) 1/2. Therefore, the optimum polymerization degree of the PE molecule in this paper is 12.
A PE molecular chain with a polymerization degree of 12 was adopted to establish a modified asphalt model. After the steps of dynamics simulation mentioned above were performed, the balanced model of PE/NR modified asphalt was finally obtained for subsequent analysis and calculation.

3.2. Radial Distribution Function

RDF, an essential indicator in molecular dynamics, can be employed to verify the correctness of molecular models.
The radial distribution function (RDF) can be applied to analyze the interaction and microscopic distribution of particles [39,43,44]. The radial distribution function is the spatial probability distribution of particle B approaching the center of particle A, which is in the r to r + d r shell. It can also be understood as the ratio of the number density of particle B to the mean density of particle B in this region.
g ( r ) = d N ρ V
d N is the number of particles B in the shell r to r + d r ; ρ is the numerical density of particle B.
V is the volume of the spherical shell, centered on atom A, with the radius R and the thickness d r . The formula is as follows:
V = 4 3 π   ( r + d r )   3 4 3 π r 3 = 4 π r 2 d r + 4 π r d r 2 + 4 3 π d r 3 4 π r 2 d r
For the RDF of amorphous crystals, there are a few peaks with different heights and sharpness at close distances, but the peak height decreases rapidly with an increase in space. At long distances, the radial distribution function tends to be evenly distributed, g   ( r ) = 1.
The RDF of the crystal still maintains a sharp peak shape at a long distance, indicating that the crystal has a long-range ordered structure, whereas the amorphous crystal has a short-range ordered structure and a long-range disorder structure.

3.3. Diffusion Coefficient

The diffusion coefficient can represent the molecular movement rate and is employed to study the changes of diverse molecular movements in NR-modified asphalt systems with and without PE. According to Einstein’s law [48], if the simulation time is long, the slope of the curve of mean square displacement ( MSD ) to time is six times the diffusion coefficient; thus, the diffusion coefficient of particles can be gained, as shown in formula (6). In molecular dynamics calculation, molecules move continuously from the initial position, and the position of each moment is different. If r   ( t ) represents the position of particle i at time t , MSD represents the average square of displacement. The formula of the diffusion coefficient is shown in formula (6).
D = lim t M S D   ( t )   6 t = lim t | r ( t ) r ( 0 ) | 2 6 t
To analyze the movement of various molecules in NR-modified asphalt, the MSD results of different molecules were obtained through molecular dynamics calculation, as shown in Figure 8. According to the functional relationship between the diffusion coefficient and MSD, the diffusion coefficient is obtained, as shown in Figure 9. The analysis shows that the saturate has the fastest movement speed, which is much higher than other types of molecules, and its diffusion coefficient reaches 0.0201. This phenomenon is most likely related to its negligible molecular weight. Aromatic ranks second, at about 0.0039.
In addition, the molecular movement velocity of asphaltene, resin, and rubber are roughly similar, showing that their diffusion coefficients are about 0.0018, likely connected with the higher molecular weight and unique shapes to some extent.
For the NR-modified asphalt system, the diffusion coefficient of NR is the lowest. More specifically, the long-chain shape of NR and its poor compatibility with the polar components of asphalt may be the main reasons restricting the movement of NR. Therefore, NR is the slowest molecule in the NR-modified asphalt system, with a diffusion coefficient of 0.00171. However, NR attracts non-polar components in asphalt (aromatic and saturate), which results in a more aromatic and saturated distribution around NR molecules. Some light components are attracted by NR molecules and transferred to the vicinity of NR molecules, so it could be noted that the light components have the highest diffusion coefficients, which are 0.0201 and 0.0039, respectively, in NR-modified asphalt.
As shown in Figure 9, the movement of molecules changed significantly due to the addition of PE. The diffusion coefficient of resin molecules increased dramatically from 0.0020 to 0.0127. This may be because there is a specific force between the PE and resin, thus accelerating the movement of resin molecules. According to the colloid theory for asphalt, resins attach molecules in asphalt to form micelles. In this study, the transfer of resin molecules caused the asphaltene molecules to lose part of the constraint, resulting in a slight increase in the diffusion coefficient of asphaltene from 0.00173 to 0.0043. Because of the rearrangement and combination of the dispersed phase, it is inevitable to promote the corresponding transfer of light components, including saturate and aromatic, that is, the corresponding increase of diffusion coefficient; for example, the diffusion coefficient of saturate increases from 0.0201 to 0.0239, and that of the aromatic component increases from 0.0039 to 0.0058.
In the statistical data analysis, the Pearson correlation coefficient is adopted to reflect the linear correlation of two random variables. For NR-modified asphalt and PE/NR-modified asphalt, M S D increased with time. Linear fitting was performed, and the Pearson correlation coefficient was analyzed. For NR-modified asphalt, the Pearson correlation coefficient of asphaltene and resin was 0.67805 and 0.82105, respectively, as shown in Table 5. Meanwhile, the PE/ NR modified asphalt coefficient is above 90%, 0.94514, and 0.99194, respectively, as shown in Table 6.

3.4. Concentration Distribution

Concentration distribution, an important indicator, can show how concentration varies with position in a region at a given time. This parameter can be acquired through the Forcite Analysis module in Materials Studio. The spatial distribution of each molecule in the NR modified asphalt in three different directions (100,010,001) and the spatial distribution of each molecule in the PE/NR modified asphalt were evaluated, and the interaction between each molecule was further analyzed and predicted. In the concentration distribution curve of a molecule, the peak represents the location at which the molecule is most abundant.
As shown in Figure 10a–c, the four components of asphalt basically follow the colloid structure theory to form a blend with an ordered arrangement. After careful analysis, the following specific rules can be found. Firstly, it might be suggested that there is adsorption between the asphaltene, NR molecules, and the light components. Specifically, the concentration curve peak of asphaltene is surrounded by the concentration curve peak of saturate because the saturate is adsorbed around by asphaltene, and the mutual attraction between the two is considerable. Again, it is surrounded by a peak of aromatic. Both indicate adsorption between the asphaltene and light components. At the same time, the concentration curve peak of NR was also surrounded by the concentration curve peak of the light component, indicating that the light component was also adsorbed around the NR molecule. Second, the peak of the NR concentration curve is far away from that of resin, as shown in Figure 10b, which may be closely related to the incompatibility between the NR and resin, thus preventing them from coexisting in the same position. Similarly, the same rule can be found in the relationship between NR and asphaltene, which is more pronounced. As shown in Figure 10a, the concentration curves of the two have apparent peaks, and the distance between the two peaks is greater, indicating that there is probably a more vigorous competition for light components and incompatibility between NR and asphaltene. This is consistent with a previous study [49]. NR has the highest binding strength to aromatic and saturate, followed by resins, and has the lowest binding capacity to asphaltenes.
Figure 10d–f shows the concentration distribution curve of PE/ NR-modified asphalt. It could be predicted from this group diagram that the addition of PE changes the colloid structure composition of asphalt. As shown in Figure 10d, both the PE concentration curve and that of asphaltene have evident peak values, and the peaks of them are far away. It could be the case that the PE molecule has a specific repulsive force with asphaltene. It can also be seen from Figure 10d that the concentration curve of the PE molecule overlaps with that of the saturate to a certain extent, and the saturation fraction is found near the periphery of the peak, which might be related to the intense attraction between the two. It is likely that the adsorption of the saturate could make the PE molecule move and form a more stable system. Therefore, PE molecules tend to further promote the formation of a new balanced and stable system by changing the interaction between the original system.

4. Conclusions

In this paper, two types of polymer materials, PE (polyethylene) and NR (cis-1, 4-polyisoprene), were used as modifiers to study the interaction mechanism in the two modified asphalts. Specifically, the interaction of molecules in the original asphalt and NR-modified asphalt was analyzed before and after the addition of PE. The change of the interaction can help the modified asphalt to become a more stable in structure and thus have better storage performance, which has a reference value for changing the molecular agglomeration phenomenon. On the other hand, the reuse of PE and NR, two typical wastes, is of great significance for the protection of the ecological environment.
The specific conclusions are as follows.
  • The optimal degree of polymerization of PE in this study model is 12. In this case, the solubility parameter between PE and NR-modified asphalt is the smallest at 0.14 < 1.3~2.1 (J/cm3) 1/2.
  • The three models in the paper are typical amorphous substances with structures of short-range order and long-range disorder. According to the RDF diagram, when r = 1.11, the g   ( r ) functions for the three models all have sharp peaks, which are 11.76, 11.74, and 11.8, respectively.
  • In NR-modified asphalt, the diffusion coefficient of saturate is the largest, at far higher than that of other types of molecules, and its value reaches 0.0201. In addition, the molecular velocity of asphaltenes, resins, and rubbers are about the same, showing that their diffusion coefficients are about 0.0018. NR is the slowest molecule.
  • In PE/NR-modified asphalt, the movement of molecules changed significantly. The diffusion coefficient of resin molecules increased considerably from 0.0020 to 0.0127. The most noticeable feature is a specific adsorption phenomenon between the PE molecule and light components.

Author Contributions

Conceptualization, Y.C. and K.H.; methodology, K.H.; software, Y.C.; validation, Y.C.; formal analysis, X.B.; investigation, Y.C.; resources, Y.C.; writing—original draft preparation, Y.C.; writing—review and editing, D.Y.; supervision, C.Y.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (2018YFB1600200), Fundamental Research Fund of Natural Science Foundation of China (52208434), China Postdoctoral Science Foundation (2022M711079), Science and Technology Project of Henan Province (192102310229), Natural Science Foundation of Henan Province (222300420142).

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank all the editors and reviewers for correcting the manuscript very carefully.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Porto, M.; Caputo, P.; Loise, V.; Eskandarsefat, S.; Teltayev, B.; Rossi, C.O. Bitumen and bitumen modification: A review on latest advances. Appl. Sci. 2019, 9, 742. [Google Scholar] [CrossRef]
  2. Presti, D.L. Recycled tyre rubber modified bitumens for road asphalt mixtures: A literature review. Constr. Build. Mater. 2013, 49, 863–881. [Google Scholar] [CrossRef]
  3. Jiang, W.; Yuan, D.D.; Shan, J.H.; Ye, W.L.; Lu, H.H.; Sha, A.M. Experimental study of the performance of porous ultra-thin asphalt overlay. Int. J. Pavement Eng. 2022, 23, 2049–2061. [Google Scholar] [CrossRef]
  4. Wu, S.; Montalvo, L. Repurposing waste plastics into cleaner asphalt pavement materials: A critical literature review. J. Clean. Prod. 2021, 280, 124355. [Google Scholar] [CrossRef]
  5. Jiang, W.; Xiao, J.J.; Yuan, D.D.; Lu, H.H.; Xu, S.D.; Huang, Y. Design and experiment of thermoelectric asphalt pavements with power-generation and temperature-reduction functions. Energy Build. 2018, 169, 39–47. [Google Scholar] [CrossRef]
  6. Fini, E.H.; Al-Qadi, I.L.; You, Z.P.; Zada, B.; Mills-Beale, J. Partial replacement of asphalt binder with bio-binder: Characterisation and modification. Int. J. Pavement Eng. 2012, 13, 515–522. [Google Scholar] [CrossRef]
  7. Sihombing, A.V.R.; Subagio, B.S.; Hariyadi, E.S.; Yamin, A. Mechanical Properties of Bio-Asphalt on Recycled Asphalt Pavement Binder. In Proceedings of the 9th International Conference on Maintenance and Rehabilitation of Pavements—Mairepav9, Zurich, Switzerland, 1–3 July 2020; Springer: Cham, Switzerland, 2020; pp. 529–538. [Google Scholar]
  8. Ren, J.L.; Zang, G.Y.; Wang, S.Y.; Shi, J.; Wang, Y.Y. Investigating the pavement performance and aging resistance of modified bio-asphalt with nano-particles. PLoS ONE 2020, 15, e0238817. [Google Scholar]
  9. Yan, K.Z.; Zhang, M.; You, L.Y.; Wu, S.H.; Ji, H.Y. Performance and optimization of castor beans-based bio-asphalt and European rock-asphalt modified asphalt binder. Constr. Build. Mater. 2020, 240, 117951. [Google Scholar] [CrossRef]
  10. Elahi, Z.; Mohd Jakarni, F.; Muniandy, R.; Hassim, S.; Ab Razak, M.S.; Ansari, A.H.; Ben Zair, M.M. Waste Cooking Oil as a Sustainable Bio Modifier for Asphalt Modification: A Review. Sustainability 2021, 13, 11506. [Google Scholar] [CrossRef]
  11. Jain, S.; Chandrappa, A.K. Rheological and chemical investigation on asphalt binder incorporating high recycled asphalt with waste cooking oil as rejuvenator. Innov. Infrastruct. Solut. 2022, 7, 1–19. [Google Scholar] [CrossRef]
  12. Zhou, X.X.; Zhao, G.Y.; Wu, S.P.; Tighe, S.; Pickel, D.; Chen, M.Z.; Adhikari, S.; Gao, Y.M. Effects of biochar on the chemical changes and phase separation of bio-asphalt under different aging conditions. J. Clean. Prod. 2020, 263, 121532. [Google Scholar] [CrossRef]
  13. Zhou, X.X.; Zhao, G.Y.; Miljković, M.; Tighe, S.; Chen, M.Z.; Wu, S.P. Crystallization kinetics and morphology of biochar modified bio-asphalt binder. J. Clean. Prod. 2022, 349, 131495. [Google Scholar] [CrossRef]
  14. Joohari, M.I.; Aziz, N.A.; Daud, N.M.; Mansor, S.; Halim, M.A. Performance of porous asphalt pavement using clay brick dust as mineral filler[C]//Journal of Physics: Conference Series. IOP Publ. 2019, 1349, 012098. [Google Scholar]
  15. Abdelsalam, M.; Yue, Y.C.; Khater, A.; Luo, D.; Musanyufu, J.; Qin, X.L. Laboratory study on the performance of asphalt mixes modified with a novel composite of diatomite powder and lignin fiber. Appl. Sci. 2020, 10, 5517. [Google Scholar] [CrossRef]
  16. Du, T.; Song, P.; Liu, L. Experimental Study on Activated Diatomite Modified Asphalt Pavement in Deep Loess Area. Processes 2022, 10, 1227. [Google Scholar] [CrossRef]
  17. Chen, X.L.; Sun, Y.S.; Han, Y.X.; Tian, Z.F. Microstructure and DSC of organic bentonite modified bitumen. J. Northeast. Univer. 2012, 33, 743. [Google Scholar]
  18. Chen, S.; Zhang, B.; He, X.Y.; Su, Y.; Liu, Q.; Xu, H. Research on mechanical-activated nanoscale bentonite and surface aging behavior of its modified asphalt. Constr. Build. Mater. 2022, 321, 126356. [Google Scholar] [CrossRef]
  19. Du, D.C.; Li, X.L.; Zheng, G.Y.; Zhang, L.Q. Study on the Performance of Polyester Fiber Modified Asphalt Mastic. Plastics 2012, 41, 42–45. [Google Scholar]
  20. Wu, S.; Ye, Q.; Li, N. Investigation of rheological and fatigue properties of asphalt mixtures containing polyester fibers. Constr. Build. Mater. 2008, 22, 2111–2115. [Google Scholar] [CrossRef]
  21. Cheng, C.; Sun, W.W.; Hu, B.; Tao, G.X.; Peng, C.; Tian, Y.J.; Wu, S.J. Analysis of the mechanism and effectiveness of lignin in improving the high-temperature thermal stability of asphalt. J. Renew. Mater. 2020, 8, 1243–1255. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Wang, X.C.; Ji, G.Y.; Fan, Z.Y.; Guo, Y.C.; Gao, W.Z.; Xin, L. Mechanical performance characterization of lignin-modified asphalt mixture. Appl. Sci. 2020, 10, 3324. [Google Scholar] [CrossRef]
  23. Zheng, Y.X.; Cai, Y.C.; Zhang, G.H.; Fang, H.Y. Fatigue property of basalt fiber-modified asphalt mixture under complicated environment. J. Wuhan Univ. Technol. Mater. Sci. Ed. 2014, 29, 996–1004. [Google Scholar] [CrossRef]
  24. Wang, W.; Cheng, Y.; Ma, G.; Tan, G.; Sun, X.; Yang, S. Further investigation on damage model of eco-friendly basalt fiber modified asphalt mixture under freeze-thaw cycles. Appl. Sci. 2018, 9, 60. [Google Scholar] [CrossRef]
  25. Arabani, M.; Shabani, A. Evaluation of the ceramic fiber modified asphalt binder. Constr. Build. Mater. 2019, 205, 377–386. [Google Scholar] [CrossRef]
  26. Guo, F.C.; Zhang, J.P.; Pei, J.Z.; Zhou, B.C.; Falchetto, A.C.; Hu, Z. Investigating the interaction behavior between asphalt binder and rubber in rubber asphalt by molecular dynamics simulation. Constr. Build. Mater. 2020, 252, 118956. [Google Scholar] [CrossRef]
  27. Guo, F.C.; Zhang, J.P.; Pei, J.Z.; Ma, W.S.; Hu, Z.; Guan, Y.S. Evaluation of the compatibility between rubber and asphalt based on molecular dynamics simulation. Front. Struct. Civ. Eng. 2020, 14, 435–445. [Google Scholar] [CrossRef]
  28. Zheng, W.H.; Wang, H.N.; Chen, Y.; Ji, J.; You, Z.P.; Zhang, Y.Q. A review on compatibility between crumb rubber and asphalt binder. Constr. Build. Mater. 2021, 297, 123820. [Google Scholar] [CrossRef]
  29. Wang, W.S.; Cheng, Y.C.; Chen, H.P.; Tan, G.J.; Lv, Z.H.; Bai, Y.S. Study on the performances of waste crumb rubber modified asphalt mixture with eco-friendly diatomite and basalt fiber. Sustainability 2019, 11, 5282. [Google Scholar] [CrossRef]
  30. Li, J.; Chen, Z.X.; Xiao, F.P.; Amirkhanian, S.N. Surface activation of scrap tire crumb rubber to improve compatibility of rubberized asphalt. Resour. Conserv. Recycl. 2021, 169, 105518. [Google Scholar] [CrossRef]
  31. Adiyanto, O.; Mohamad, E.; Abd Razak, J. Systematic Review of Plastic Waste as Eco-Friendly Aggregate for Sustainable Construction. Int. J. Sustain. Constr. Eng. Technol. 2022, 13, 243–257. [Google Scholar]
  32. Ben Zair, M.M.; Jakarni, F.M.; Muniandy, R.; Hassim, S. A brief review: Application of recycled polyethylene terephthalate in asphalt pavement reinforcement. Sustainability 2021, 13, 1303. [Google Scholar] [CrossRef]
  33. Leng, Z.; Padhan, R.K.; Sreeram, A. Production of a sustainable paving material through chemical recycling of waste PET into crumb rubber modified asphalt. J. Clean. Prod. 2018, 180, 682–688. [Google Scholar] [CrossRef]
  34. Qu, X.; Wang, D.W.; Wang, L.B.; Huang, Y.C.; Hou, Y.; Oeser, M. The state-of-the-art review on molecular dynamics simulation of asphalt binder. Adv. Civ. Eng. 2018, 2018, 4546191. [Google Scholar] [CrossRef]
  35. Yu, C.H.; Hu, K.; Chen, Y.J.; Li, Q.; Zhang, T.L.; Chen, G.X. The Influence of Styrene/butadiene Ratio on the Self-assembly Behavior of SBS Modified Asphalt[C]//IOP Conference Series: Earth and Environmental Science. IOP Publ. 2021, 668, 012060. [Google Scholar]
  36. Chu, L.; Luo, L.; Fwa, T.F. Effects of aggregate mineral surface anisotropy on asphalt-aggregate interfacial bonding using molecular dynamics (MD) simulation. Constr. Build. Mater. 2019, 225, 1–12. [Google Scholar] [CrossRef]
  37. Hu, B.; Huang, W.K.; Yu, J.L.; Xiao, Z.C.; Wu, K.H. Study on the adhesion performance of asphalt-calcium silicate hydrate gel interface in semi-flexible pavement materials based on molecular dynamics. Materials 2021, 14, 4406. [Google Scholar] [CrossRef]
  38. Hu, D.L.; Pei, J.Z.; Li, R.; Zhang, J.P.; Jia, Y.S.; Fan, Z.P. Using thermodynamic parameters to study self-healing and interface properties of crumb rubber modified asphalt based on molecular dynamics simulation. Front. Struct. Civ. Eng. 2020, 14, 109–122. [Google Scholar] [CrossRef]
  39. Yan, L.M.; Zhu, S.H. Theory and Practice of Molecular Dynamics Simulation; Science Press: Beijing, China, 2013. [Google Scholar]
  40. Sun, H.; Rigby, D. Polysiloxanes: Ab initio force field and structural, conformational and thermophysical properties. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 1997, 53, 1301–1323. [Google Scholar] [CrossRef]
  41. Sun, H. An ab initio force-field optimized for condensed-phase applications overview with details on alkane and benzene compounds. J. Phys. Chem. 1998, 102, 7338–7364. [Google Scholar] [CrossRef]
  42. Sun, H.; Ren, P.; Fried, J.R. The COMPASS force field: Parameterization and validation for phosphazenes. Comput. Theor. Polym. Sci. 1998, 8, 229–246. [Google Scholar] [CrossRef]
  43. Bunte, S.W.; Sun, H. Molecular modeling of energetic materials: The parameterization and validation of nitrate esters in the COMPASS force field. J. Phys. Chem. B 2000, 104, 2477–2489. [Google Scholar] [CrossRef]
  44. Sun, H.; Jin, Z.; Yang, C.W.; Akkermans, R.L.C.; Robertson, S.H.; Spenley, N.A.; Miller, S.; Todd, S.M. COMPASS II: Extended coverage for polymer and drug-like molecule databases. J. Mol. Model. 2016, 22, 1–10. [Google Scholar] [CrossRef] [PubMed]
  45. Su, M.M.; Si, C.D.; Zhang, Z.P.; Zhang, H.L. Molecular dynamics study on influence of Nano-ZnO/SBS on physical properties and molecular structure of asphalt binder. Fuel 2020, 263, 116777. [Google Scholar] [CrossRef]
  46. Zhang, Y. Fundamentals of Computational Materials Science; Beijing University of Aeronautics and Astronautics Press: Beijing, China, 2007. [Google Scholar]
  47. Lan, Y.H.; Li, D.H.; Yang, R.J.; Liang, W.S.; Zhou, L.X.; Chen, Z.W. Computer simulation study on the compatibility of cyclotriphosphazene containing aminopropylsilicone functional group in flame retarded polypropylene/ammonium polyphosphate composites. Compos. Sci. Technol. 2013, 88, 9–15. [Google Scholar] [CrossRef]
  48. D’hooge, D.R.; Van Steenberge, P.H.M.; Reyniers, M.F.; Marin, G.B. The strength of multi-scale modeling to unveil the complexity of radical polymerization. Prog. Polym. Sci. 2016, 58, 59–89. [Google Scholar] [CrossRef]
  49. Hu, K.; Yu, C.H.; Chen, Y.J.; Li, W.; Wang, D.D.; Zhang, W.G. Multiscale mechanisms of asphalt performance enhancement by crumbed waste tire rubber: Insight from molecular dynamics simulation. J. Mol. Model. 2021, 27, 1–14. [Google Scholar] [CrossRef]
Figure 1. The SEM images of (a) BT and (b) OBT.
Figure 1. The SEM images of (a) BT and (b) OBT.
Polymers 14 04087 g001
Figure 2. Models of representative molecules in bitumen. (a) asphaltene; (b) resin; (c) saturate; (d) aromatic.
Figure 2. Models of representative molecules in bitumen. (a) asphaltene; (b) resin; (c) saturate; (d) aromatic.
Polymers 14 04087 g002
Figure 3. The molecular chain for NR.
Figure 3. The molecular chain for NR.
Polymers 14 04087 g003
Figure 5. Models (a) original asphalt; (b) NR-modified asphalt; (c) PE/NR-modified asphalt.
Figure 5. Models (a) original asphalt; (b) NR-modified asphalt; (c) PE/NR-modified asphalt.
Polymers 14 04087 g005
Figure 6. PE molecular model with different degrees of polymerization (a) 6 (b) 12 (c) 24 (d) 30 (e) 42 (f) 78.
Figure 6. PE molecular model with different degrees of polymerization (a) 6 (b) 12 (c) 24 (d) 30 (e) 42 (f) 78.
Polymers 14 04087 g006
Figure 7. The difference in the solubility parameter.
Figure 7. The difference in the solubility parameter.
Polymers 14 04087 g007
Figure 8. MSD for modified asphalt (a) NR-modified asphalt; (b) PE/NR-modified asphalt.
Figure 8. MSD for modified asphalt (a) NR-modified asphalt; (b) PE/NR-modified asphalt.
Polymers 14 04087 g008
Figure 9. Diffusion coefficient of each component in NR modified asphalt and PE/NR modified asphalt.
Figure 9. Diffusion coefficient of each component in NR modified asphalt and PE/NR modified asphalt.
Polymers 14 04087 g009
Figure 10. Concentration distribution of each component for NR-modified asphalt and PE/NR-modified asphalt. (The relative concentration of NR-modified asphalt in different directions, including (a) 1 0 0, (b) 0 1 0, and (c) 0 0 1; The relative concentration of PE/NR-modified asphalt in different directions includes (d) 1 0 0, (e) 0 1 0, and (f) 0 0 1).
Figure 10. Concentration distribution of each component for NR-modified asphalt and PE/NR-modified asphalt. (The relative concentration of NR-modified asphalt in different directions, including (a) 1 0 0, (b) 0 1 0, and (c) 0 0 1; The relative concentration of PE/NR-modified asphalt in different directions includes (d) 1 0 0, (e) 0 1 0, and (f) 0 0 1).
Polymers 14 04087 g010
Table 1. Specific information about the asphalt model.
Table 1. Specific information about the asphalt model.
ComponentsMolecular FormulaNumber RatioRelative Molecular Mass Mass Fraction (%) NumAtoms
AsphalteneC149H177N3O2S212106.1919.734333
ResinC59H85NOS3856.39524.071147
SaturateC22H465310.61014.55168
AromaticC46H50S7634.96641.64597
Table 2. NR-modified asphalt model.
Table 2. NR-modified asphalt model.
Serial NumberRepresentative MoleculeChemical CompositionMolecule NumberRelative Molecular Mass Mass Fraction
1AsphalteneC149H177N3O2S212106.19017.90%
2ResinC59H85NOS3856.39521.84%
3SaturateC22H465310.61013.20%
4AromaticC46H50S7634.96637.78%
5NRC80 H13011091.9209.28%
Table 3. PE/NR-modified asphalt model.
Table 3. PE/NR-modified asphalt model.
Serial NumberRepresentative MoleculeChemical CompositionMolecule Number Relative Molecular Mass Mass Fraction
1AsphalteneC149H177N3O2S212106.19017.40%
2ResinC59H85NOS3856.39521.23%
3SaturateC22H465310.61012.83%
4AromaticC46H50S7634.96636.72%
5NRC80 H13011091.9209.02%
6PEC24H501338.6642.80%
Table 4. Solubility parameters of different PE molecules.
Table 4. Solubility parameters of different PE molecules.
Number1234567891011121314
Polymerization degree612182430364248546066727884
Solubility parameter17.8517.0916.6615.6614.3315.5614.9915.7115.4716.0814.7015.9215.7915.02
Table 5. Data analysis of MSD for NR-modified asphalt.
Table 5. Data analysis of MSD for NR-modified asphalt.
Equation y = a + b x
PlotasphalteneresinsaturatearomaticNR
WeightNo Weighting
Intercept1.68015 ± 0.148441.55818 ± 0.10932.44376 ± 0.156381.38694 ± 0.027992.007 ± 0.06071
Slope0.01035 ± 0.00290.01188 ± 0.002130.12047 ± 0.003050.02361 ± 5.46305E-40.01028 ± 0.00118
Residual Sum of Squares1.284240.696311.425370.045660.21482
Pearson’s r0.678050.821050.995220.996010.91308
R-Square (COD) 0.459750.674120.990460.992040.83371
Adj. R-Square0.423730.652390.989830.99150.82263
Table 6. MSD data analysis for PE/NR modified asphalt.
Table 6. MSD data analysis for PE/NR modified asphalt.
Equation y = a + b x
PlotasphalteneresinsaturatearomaticNR
WeightNo Weighting
Intercept1.61155 ± 0.118771.75603 ± 0.128322.74225 ± 0.258821.53794 ± 0.052432.10403 ± 0.06517
Slope0.02598 ± 0.002320.07596 ± 0.00250.14313 ± 0.005050.03484 ± 0.001020.0455 ± 0.00127
Residual Sum of Squares0.822140.959783.904580.16020.24752
Pearson’s r0.945140.991940.990790.993590.99419
R-Square (COD) 0.893290.983950.981660.987220.98842
Adj. R-Square0.886180.982880.980440.986370.98764
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Y.; Hu, K.; Yu, C.; Yuan, D.; Ban, X. Study of the Microscopic Mechanism of Natural Rubber (Cis-1, 4-Polyisoprene, NR)/Polyethylene (PE) Modified Asphalt from the Perspective of Simulation. Polymers 2022, 14, 4087. https://doi.org/10.3390/polym14194087

AMA Style

Chen Y, Hu K, Yu C, Yuan D, Ban X. Study of the Microscopic Mechanism of Natural Rubber (Cis-1, 4-Polyisoprene, NR)/Polyethylene (PE) Modified Asphalt from the Perspective of Simulation. Polymers. 2022; 14(19):4087. https://doi.org/10.3390/polym14194087

Chicago/Turabian Style

Chen, Yujing, Kui Hu, Caihua Yu, Dongdong Yuan, and Xiaoyi Ban. 2022. "Study of the Microscopic Mechanism of Natural Rubber (Cis-1, 4-Polyisoprene, NR)/Polyethylene (PE) Modified Asphalt from the Perspective of Simulation" Polymers 14, no. 19: 4087. https://doi.org/10.3390/polym14194087

APA Style

Chen, Y., Hu, K., Yu, C., Yuan, D., & Ban, X. (2022). Study of the Microscopic Mechanism of Natural Rubber (Cis-1, 4-Polyisoprene, NR)/Polyethylene (PE) Modified Asphalt from the Perspective of Simulation. Polymers, 14(19), 4087. https://doi.org/10.3390/polym14194087

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop