1.
Introduction
Waste rubber powder modified asphalt (RA) has excellent engineering performance in reducing driving noise, delaying reflective cracking, minimizing pavement water damage and withstanding heavy traffic [1,2]. In addition, RA has a significant role in improving the high temperature fatigue resistance and resistance to low temperature cracking of base asphalt [3]. As a polymer material which is hard to degrade by natural conditions, scrap tires can be used in road engineering after crushing to improve road performance and driving comfort [4,5]. At the same time, the harmless treatment and resource utilization of solid waste are realized to avoid "black pollution, " which is aligned with the concept of green sustainable development and has good application prospects and environmental protection value [6,7]. However, the incompatibility of the rubber and base asphalt make them vulnerable to segregation during storage and transportation, which limits the application and further development of RA in engineering [8,9,10,11].
Nanomaterials are a common class of additives for enhancing the storage stability of RA. Graphene has shown excellent performance as an asphalt modifier. Wang et al. [12] evaluated the viability of adding carbon nanoparticles to asphalt binder. Softening point differences, rheological responses and micromorphology were tested, demonstrating that the carbon nanomaterial-modified asphalt was capable to satisfy the storage stability requirements. In accordance with Chen et al. [13], modified bitumen storage stability could be improved by the incorporation of graphene. Graphene enhances the viscoelasticity and the high- and low-temperature characteristics of RA, and it also increases the adhesion of ore and asphalt. Li et al. [14] discovered that the incorporation of nanoscale graphene to bitumen could remarkably increase the material's high-temperature resilience and temperature sensitivity, leading to improved thermal storage stability of modified bitumen. Singh et al. [15] analyzed the compatibility of graphene oxide (GO) rubber powder composite modified hot mix bitumen together with the structure and morphology of the modified bitumen. The results showed that WRP and GO were completely compatible with asphalt, offering better road performance. Meng et al. [16] discovered that the rheological and mechanical performances of RA incorporated with graphene were improved, with higher elastic recovery and deformation resistance. The phase distribution was also more uniform and continuous, resulting in improved storage stability and compatibility. Liu et al. [17] discovered that the doping of graphene and carbon black composites (GC) greatly enhanced the fatigue damage resistance and thermally activated healing characteristics of RA. Through physical and chemical interactions, the loose GC reacted with asphalt and formed tight connections. Ferrari et al. [18] demonstrated that an important indicator of graphene properties is the number of layers. From the above findings, it is clear that, as a new nanomaterial, graphene could improve the storage stability performance of RA. However, most existing studies mainly emphasized the reinforcement effect of graphene on RA. Although they can to some extent reflect the property changes and microscopic morphology of asphalt, the majority of them are simply phenomenological or qualitative descriptions, making it hard to clearly reveal the mechanism of the interaction between graphene and RA.
The properties of the material itself are impacted by the particle sizes, specific surface areas and other metrics that differ between various sizes and layers of graphene [19]. The characteristics of graphenes with different layers and their respective advantages in terms of performance and applications in emerging fields, such as novel energy materials, composites, biomaterials and sensors, have been summarized [20]. However, in the road materials field, there are insufficient studies about graphene size effect on modified asphalt properties, which results in the inability to optimize the modification effect of graphene. Owing to advances in graphene preparation technology, it is now possible to produce graphene with specific layers and sizes. To better guide the design and development of graphene materials and achieve optimal asphalt performance enhancement when graphene is used as a modifier, it is necessary to explore further the effect of graphene microstructural properties on RA storage stability and the modification mechanism [21].
Molecular dynamics simulation (MD) characterizes the behavior of materials at an atomic or molecular level, especially to describe the structure-property relationships of materials, which could be used to explain the molecular motion behaviors of multiple components in modified asphalt [22]. At present, the modifying effect of graphene on RA storage stability still lacks quantitative explanation. Graphene/rubber composite modified asphalt (GRA) was taken as the research object. The interaction behavior in GRA and the graphene size effect on modified asphalt storage stability were studied by combining multiple-scale tests with MD.
In this paper, the graphene size effect on RA storage stability is discussed in combination with experimental analysis and MD. The influence of different graphene types in modified asphalt on the enhancement of storage stability was investigated from a macroscopic perspective using storage stability tests and dynamic shear rheometry (DSR) tests. The effect of different layers and sizes of graphene in GRA was explained from a microscopic perspective by fluorescence microscopy experiments (FM). Meanwhile, the size effect of graphene and its mechanism on the storage stability of GRA were investigated by MD for the application of graphene in modified asphalt.
2.
Objective
This research uses multiple-scale experimental characterizations and molecular dynamics simulation to study the enhancement effect of graphene on RA storage stability. An emphasis was placed on how the geometric properties of graphene affected the modification effect.
3.
Materials and methodology
3.1. Materials and preparation
3.1.1. Materials
AH-70 road asphalt was selected for this study, and the performance indices of the original asphalt used were examined based on the Chinese standard JTG E20-2011, as shown in Table 1. The waste rubber powder (WRP) was produced from whole automobile radial tires, with a fineness of 40 mesh. Natural rubber (NR), butadiene rubber (BR) and styrene-butadiene rubber (SBR) are the three major ingredients of the chosen WRP. The proportion of each component was NR: BR: SBR = 36.66: 8.88: 54.46% [23]. According to Chinese standards JT/T797-2011, GB/T 14837 and GB/T 3516, the technical indices of the measured WRP are listed in Table 2. Three types of graphene particles, e.g., single-layer graphene (SLG), few-layer graphene (FLG) and multi-layer graphene (MLG), were used in this study, and their performance indices are listed in Table 3.
Using scanning electron microscopy (SEM) at a magnification of 7000, three different types of graphene were observed. Figure 1 depicts the apparent morphologies of these particles. SLG shows a thin sheet layer with a relatively smooth and flat surface. FLG films show a large number of folds on the surface, showing a silk-like appearance. MLG is multiple lamellar layers stacked together, with a rough surface and some protrusions. The results illustrated significant differences in the apparent morphologies of SLG, FLG and MLG.
3.1.2. Preparation of modified asphalt
First, an oven was used to melt the base asphalt. Next, the asphalt was transferred into a heating sleeve with the temperature set to 180 ℃. When the asphalt temperature reached 175 ℃, low speed shearing was started at a shearing speed of 1000 r/min to eliminate air bubbles and prevent swelling after adding the WRP. Then, 20% pre-dried WRP and 5% ultrasonically dispersed graphene were gradually added while being stirred with a glass rod using the external blending method [24]. When 180 ℃ was reached, the shear speed was set to 4000 r/min for 50 minutes, and finally, the modified asphalt specimens were placed in a 175 ℃ oven for 45 minutes of continuous development. The preparation process is illustrated in Figure 2.
3.2. Property Tests
The flow chart of this work is illustrated in Figure 3.
3.2.1. Storage stability test
According to the specification ASTM D7173, single-layer graphene-modified rubber asphalt (SGRA), few-layer graphene-modified rubber asphalt (FGRA), multi-layer graphene-modified rubber asphalt (MGRA) and RA were taken for 48 h storage stability testing. After curing, the material was cut into three equally divided segments, of which the upper and bottom two were taken for subsequent experiments. The 48 h storage stability test simulated the process of a binary mixture of asphalt-rich phase and polymer-rich phase dissociation under high-temperature storage in the mixing plant. The softening points of the specimens of different tube segments were evaluated and compared after sufficient mixing, and storage stability was indicated by softening point test (ASTM D36) and DSR test (ASTM D7175).
In accordance with ASTM D7175 and related literature [25,26], the separation ratios (RS) under high-temperature and low-frequency conditions are reliable metrics to evaluate RA storage stability, followed by the percentage of separation (SP). DSR tests were performed on samples from the upper and bottom segments of modified bitumen tubes after the storage stability experiments. Performing DSR at 58 ℃, frequency sweep experiments were performed from 0.15 to 30 Hz. The dynamic shear modulus G* and phase angle δ of specimens at 0.15 and 1.5 Hz were selected to calculate RS according to Eq (1) [25,27]. Using DSR's LAS-test mode, the oscillatory time sweep tests were performed at 10 rad/s and 58 ℃, and the SP was calculated according to Eq (2) [25,27].
where u represents the upper specimen of test tubes, b represents the bottom specimen of test tubes, and avg indicates the average of the values taken for the upper and bottom segments.
3.2.2. Fluorescence microscopy experiments (FM)
Fluorescence microscopy experiments (FM) were performed using a German Leica DM-40 ortho-mounted fluorescence microscope with excitation wavelength of 470 ± 20 nm and emission sheet wavelength of 515 nm. The effects of different sizes and layers of graphene on the degree of separation from the asphalt phase of the polymer phase were investigated by comparing fluorescence microscopy images for the upper and bottom segments of the specimen after storage stability tests. The preparation process of fluorescence microscopy experimental specimens was the following: The asphalt specimen was heated to flow dynamics, and a drop of asphalt was taken onto the slide. Then, the slide was covered and pressed flat and thin until there was some light transmission. If the sample was too thick, it could be heated in the oven and pressed again until there was light transmission. The cooled specimen could be placed in a fluorescence microscope for observation.
3.3. Simulation models and methods
3.3.1. Modeling method
Molecular dynamics simulations were implemented using the COMPASSII force field of the Materials Studio (MS) software. The asphalt molecular models were constructed according to the chemical four-component method, using the AAA-1 model suggested by Li et al. [28]. Asphalt in the system is composed of 12 molecules. The number of each molecule of asphalt is determined according to the asphaltene, saturate, aromatic and resin contents measured by the four-component test (ASTM D4124). Molecular structures of twelve components of asphalt are illustrated in Figure 4. Numbers of the twelve molecules in the models are presented in Figure 5, and the ratios of four components in the models and experimentally obtained values are shown in Figure 6. Natural rubber (NR), butadiene rubber (BR) and styrene-butadiene rubber (SBR) are the three main ingredients of the WRP used [30]. The monomers cis-1, 4-polybutadiene and cis-1, 4-polyisoprene are polymerized to form the homopolymers NR and BR, respectively. SBR is a random copolymer primarily composed of trans-1, 4-butadiene, 1, 2-butadiene, styrene and cis-1, 4-butadiene [23]. The structures of rubber molecules and their monomers are shown in Figure 7. In this study, graphene molecules with different combinations of layers and sizes were constructed. Selected layer numbers were 1, 3, 5, 7, 9, and the sizes were selected in supercell with three sizes of repeating monomer, u = v = 3, 6, 9, corresponding to edge lengths of 7.38, 14.76 and 22.14 Å, respectively. In Figure 8, the constructed graphene molecular models are shown.
Based on the above work, after Geometry Optimization of each molecule, the RA model and GRA models were constructed in the Amorphous Cell.
3.3.2. Simulation method
Once the modified asphalt system model was completed, the geometry was first optimized in Forcite module until convergence conditions were met, and the Smart algorithm was chosen to find the lowest potential energy conformation. Then, to remove irrational energies from each model, asphalt molecules were annealed in the NPT ensemble (constant-pressure, constant-temperature) at 300 to 500 K and an atmospheric pressure of 1.01 × 10-4 GPa. The model was then sufficiently relaxed to reach the equilibrium state after 100,000 steps of dynamics simulations in the NVT ensemble (canonical ensemble). After completing molecular dynamics simulations for each asphalt system at 298 and 436 K with 0.2 fs time step and 100,000 steps of NPT ensemble, the obtained files were subjected to parametric calculations.
3.3.3. Model Verification
In a system, the normalized probability of a particle appearing radially around a reference particle at a given distance r is described by the radial distribution function (RDF) [29]. The trajectory file of each modified asphalt system model after dynamic simulation was operated by the Analysis function of MS to obtain the intramolecular RDF. In Figure 9, the RDFs of the constructed models show multiple peaks in the range of 0–4Å, which gradually become flat and converge to 1 after r > 5Å. This demonstrates that the constructed asphalt models are non-crystalline structures with remotely disordered and proximally ordered internal masses [30]. The presence of the peaks in RDF also reflects the fact that the type of interaction between non-bonded atoms in the asphalt models is mainly van der Waals forces and hydrogen bonding, which is in line with actual asphalt [30]. Figure 10 displays the evolution of system parameters throughout the MD. After the simulation ran for 25 ps, it was evident that the asphalt model had almost reached equilibrium.
4.
Results and discussion
4.1. Multiple-scale experimental characterizations
4.1.1. Storage stability test results
The storage stabilities of the four modified asphalts were assessed by three indicators: softening point difference, RS value and SP value. After storage stability testing, segments at the upper and bottom parts of the tubes were tested for the softening point, the frequency sweep and the oscillatory time sweep of DSR. The results obtained are summarized in Table 4.
Softening point difference results revealed the values of SGRA, FGRA and MGRA were significantly smaller than that of RA. Compared to the RA, with a softening point difference of 4.2 ℃, the decrease in GRA values ranged from 11.9 to 45.24%. The data revealed that doping with graphene could alleviate the phase separation of RA, with the most obvious improvement with FLG. The improvement effect on asphalt storage stability was FLG > SLG > MLG.
As the RS value approaches the range of 0.8–1.2, the trend of separating the two phases of polymer-modified asphalt is less pronounced [25]. The RS values of RA were 0.21 at 0.15 Hz frequency condition and 0.11 at 1.5 Hz frequency condition, which were apparently lower than the specified range, indicating that more severe phase separation occurred from RA in the thermal storage stability experiments. The SGRA, FGRA and MGRA all had significantly higher RS values than those of RA, which proved that the incorporation of graphene enhanced the bonding of the polymer phase to the asphalt phase under thermal storage conditions. A similar conclusion can be obtained in the SP results. The SP values of SGRA, FGRA and MGRA were significantly lower than that of RA. Lower SP values mean that the modified asphalt is less likely to segregate under thermal storage conditions. FGRA in 0.15 and 1.5 Hz frequency conditions can provide the desired RS value of 0.8–1.2, and the SP value of FGRA is the lowest among the four types of modified asphalt. These results suggest that FGRA is virtually free from phase separation and has significantly improved storage stability compared to RA. According to the analysis of RS values, the improvement effect on asphalt storage stability was FLG > MLG > SLG. However, according to the analysis of SP values, the improvement effect on asphalt storage stability was FLG > SLG > MLG.
From the above laboratory phenomena, FLG showed the best performance in improving RA storage stability. On the other hand, the enhancement efficiencies of SLG and MLG are different when different evaluation indices are adopted, so the graphene size and the number of layers will be further subdivided in subsequent segments.
4.1.2. Fluorescence microscope
Samples of the upper and bottom segments of four modified asphalts were observed by fluorescence microscopy after storage stability tests. The distribution state of rubber powders in asphalt after a long storage time visually reflects the separation characteristics of the polymer phase and asphalt phase. The images are shown in Figure 11. The fluorescent spots in the images obtained by the fluorescence microscope were counted using ImageJ software, and the difference in the number of fluorescent spots between the upper and bottom segments of each specimen was calculated, with the results shown in Figure 12.
The polymer in RA was distributed in particle form in the continuous phase of the asphalt. As shown in Figure 11 (1a), (1b), the rubber powder particles at the upper part of the RA were fewer than at the bottom segments, indicating obvious segregation inside. In the images of SGRA, FGRA and MGRA, the fluorescence spots difference between the upper and bottom segments is small. From Figure 12, the differences in the number of fluorescent spots in the upper and bottom two segments of SGRA, FGRA and MGRA decreased by 78.5, 87.9 and 75.8%, respectively, compared with RA, indicating that the incorporation of graphene obviously reduced the phase separation in GRA.
A possible reason for segregation in RA is that the asphalt viscosity has decreased under thermal storage conditions, and the network of polymer-asphalt linkage partially fails, making the binding effect on the rubber powder in the system smaller [31]. Graphene, with its large surface area, acts as a physical cross-link from rubber to asphalt and promotes bonding between them to limit the movement of asphalt molecules under thermal storage conditions, thus improving the storage stability of RA [32,33]. Graphene particles of different sizes and layers have different physical and surface chemical properties, which have different degrees of enhancing effects on the storage stability performance of RA.
4.2. Molecular dynamics simulation results
4.2.1. Solubility parameter (δ)
According to the theory of mixing heat of polymer blends, the cohesive energy density (CED) is the energy needed to counteract the intermolecular forces for 1 mol of cohesion in a unit volume, and the open square of the CED is the solubility parameter (δ) [34,35]. Solubility parameters are used to gauge the compatibility between materials. The tendency of segregation decreases as the solubility parameter difference (Δδ) between the materials decreases [36]. It is calculated as follows:
where Ecoh is the condensed matter cohesion energy, and V is the molar volume.
The solubility parameters obtained from the simulations are shown in Figure 13. The Δδ between RA and graphene can indicate the compatibility between those two phases, and good compatibility is possible when Δδ ≤ 4.10 (J/cm3)1/2 [36]. As shown in Figure 13, with the graphene size of u = v = 3, the solubility difference falls in the range of 2.4–2.7 (J/cm3)1/2 at 298 K and in the range of -0.2–0.1 (J/cm3)1/2 at 436 K. These values satisfy the condition that Δδ ≤ 4.10 (J/cm3)1/2, indicating that small-sized graphene molecules (u = v = 3) are well compatible with RA. At high temperatures, the solubility difference reduces and almost equals to zero, indicating that the increased temperature could avoid the precipitation of graphene due to incompatibility. If the compatibility between asphalt and graphene is poor, some graphene outside the network structure will fail to provide reinforcement due to stress concentration caused by agglomeration [37,38].
As graphene size increases, the solubility of graphene molecules increases greatly. Therefore, the difference in solubility between graphene and RA shows an increasing trend, indicating the worsening compatibility between RA and graphene. When compared with the influence of the molecular size, the variation in layer number has relatively little impact on solubility parameters. Overall, the results indicate that small-sized graphene is more compatible with RA, and the solubility parameters are more related to the change in graphene size rather than its layers.
4.2.2. Binding energy(∆E)
Based on the theory of weak intermolecular interactions, the final frame of the MD trajectory is selected to calculate the binding energy (ΔE) of graphene molecules with different combinations of layers and sizes in GRA system models [39]. A bound system maintains its integrity by having a system energy that is lower than the total energy of its components. The process of breaking this bound state requires energy absorption, and the more energy that is required, i.e., the larger the absolute value of ΔE is, the more stable the system will be when it is bound [40]. The ΔE values at 298 and 436 K were used to evaluate the association properties of graphene and RA at the interface under normal temperature and thermal storage conditions. The calculation of ∆E is given in Eq (4).
where ∆E represents the binding energy between graphene and RA, E(RA-Graphene) represents the energy of graphene rubber composite modified asphalt, E(RA) represents RA's energy, and E(Graphene) represents graphene's energy.
The ∆E values of different sizes and layers of graphene in the GRA models are listed in Table 5. The binding energy exhibits negative values (i.e., ∆E < 0) in all models. This indicates that the components in graphene and RA are mutually attracted and that the doping of graphene enhances the stability of the combined asphalt phase and polymer phase [40].
According to the data in Table 5, graphene with the same number of layers shows an increasing trend of ∆E as the size increases. This may be because, as the size becomes larger, the contact area with surrounding components increases, with the contact surface being more prone to twisting and folding, thus increasing the roughness. The layered graphene structure forms an intercalation structure with the chain structure of saturated molecules, which has stable physical cross-linking and improves the thermal storage stability between asphalt and polymer [41]. For graphene molecules of size u = v = 6 and u = v = 9, ∆E tends to rise with increasing number of layers. When the size is same, the FLG and MLG molecules can produce slip between the layers to increase the contact area with surrounding molecules [20], and the intercalation structure formed easily with the surrounding molecules improves the intermolecular binding energy and forms a graphene-asphalt-polymer association network [32,33].
However, only small-sized graphene (u = v = 3) can provide the desired compatibility between RA and graphene. Therefore, the binding energy was further characterized for the series with the size of u = v = 3. Regarding the effect of graphene layers, the ∆E of RA and graphene with increasing layers shows an increasing and then decreasing trend. Meanwhile, FLG and MLG provide notably higher binding energies than SLG. For FLG, its small size has good compatibility with asphalt. Also, its interlayer slip and fold increase the contact surface area and roughness, thus enhancing the intermolecular interaction force [20]. The highest binding energy for a layer number of 3 in molecular simulation results is consistent with the macroscopic experiments, where FLG has the best effect on the improvement of the storage stability of RA. The lower binding energy provided by the small-sized graphene with 5–9 layers is consistent with the slightly lower effect of MLG on improving storage stability than FLG in macroscopic experiments. This might be associated with the fact that graphene with 3 to 4 layers is more inclined to undergo interlamellar slippage, which was discovered by Gong et al. [42].
Overall, based on the balance between compatibility and binding energy, the simulation results point out that graphene molecules with size u = v = 3 and layer number of 3 layers have the best modification effect on RA storage stability.
4.2.3. Relative concentration distribution
The relative concentration distribution (RCD) of asphaltenes, light components and rubber components in the model at 436 K was obtained in the Forcite module. Based on the theory of colloidal structure, colloids are composed of micelles formed by resins adsorbed around asphalt as dispersions, dispersed in a dispersion medium consisting of saturates and aromatics, and swollen rubber powders are dispersed in asphalt. Homogeneous dispersion of the asphaltenes leads to a stable colloidal structure through interaction with lighter components [43]. Graphene has a similar morphology to asphaltene, and it can act as a nucleus for the asphaltic micelles when dispersed in the asphaltic colloid [12]. This could enhance the link between the polymer and asphalt phases, reducing the possibility of segregation. Its excellent thermal conductivity also facilitates rapid and homogeneous dispersion of bitumen for heat transfer, thus improving the high temperature rheology of bitumen [37,38]. Based on the above findings, the GRA model of graphene with a layer number of 3 and a size of u = v = 3 was chosen to analyze the RCD of constituents and compared with the results of the RA model. Figure 14 shows the distribution of the constituents in the RA and GRA models. Figure 15 shows the RCDs of asphaltenes, light components and rubber powder molecules in the directions of the X, Y and Z coordinate axes in the RA and GRA models.
As can be seen from Figure 14, asphaltene dispersion of the GRA model is more uniform than that of the RA model. As shown in Figure 15(a), (b), both in the RA and GRA models, there are multiple asphaltene peaks in the triaxial directions. However, more asphaltene peaks appear in the GRA model, and these peaks are lower and uniformly distributed throughout the model cross-section, indicating a more uniform asphaltene dispersion in the system with graphene, which is consistent with the direct observation of the model. From Figure 15(c), (d), the light components in the GRA model have more peaks and fewer peak fluctuations compared to those in the RA model, and they are more uniformly distributed in all three directions. From Figure 15(e), (f), the rubber components have obvious peaks in all three directions in the GRA model, the peaks are higher, and the polymer phase is more aggregated, which is consistent with the molecular distribution characteristics of rubber powder in the GRA model in Figure 14 [39].
Comparing Figure 15(a), (e), it can be found that the asphaltene peaks and the rubber components peaks in the RA model have a slight overlap and close distance in triaxial directions. In contrast, in the GRA model as illustrated in Figure 15(b), (f), the asphaltene peaks do not overlap with the rubber component peaks. Moreover, the asphaltene peaks appear at the trough of the RCD curve of the rubber components. This indicates that both asphaltenes and rubber molecules are dispersed in the light components, which facilitates the spatial distribution of rubber molecules to absorb the light components and dissolve them. This strengthens the weak interface of asphalt and WRP, increasing the cohesion of RA and improving storage stability [39].
5.
Conclusions
To investigate the graphene size effect on the storage stability of RA and explore the better application and development direction of graphene in modified asphalt, this research combined macro and micro experiments and MD to reach the following conclusions:
1) The storage stability test and DSR results indicated that graphene could enhance the storage stability of RA, with FLG being the most effective in improving the phase separation that exists in RA at high temperature storage conditions.
2) The fluorescence microscopy experiment indicated that graphene could promote the bonding of rubber and asphalt and prevent the partial failure of polymer-asphalt connection network during thermal storage, thus improving the stability of RA. In addition, GRA doped with few-layer graphene has a superior microstructure, which is in line with the conclusions obtained from macroscopic experiments.
3) The molecular dynamics simulation indicated that graphene with a layer number of 3 and a size of u = v = 3 was the most effective in enhancing the storage stability of RA. Graphene layer number mainly influences the binding energy, while the graphene size primarily dominates both solubility parameters and binding energy.
4) The Δδ and ΔE results indicated that both Δδ and ΔE between RA and graphene increase as the graphene size increases, indicating a decrease in compatibility and an increase in binding energy. For GRA models with graphene sizes of 6 and 9, ΔE rises with increasing number of layers, while for GRA models with graphene size of 3, an optimal layer number exists to ensure a lower Δδ and a higher ΔE.
5) The results of the relative concentration distributions support that the presence of graphene is conducive to the spatial distribution of asphaltene, rubber powder and light components, as well as the enhancement of the weak interfaces of asphalt and rubber powder.
Acknowledgments
Authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (52278239 and 51978080), Key Research and Development Project of Hunan Provincial Science and Technology Department (2021SK2044) and Changsha University of Science & Technology Graduate Student Research Innovation Program (CX2021SS114).
Conflict of interest
The authors declare there is no conflict of interest.