Synergetic Mechanism of Multiple Industrial Solid Waste-Based Geopolymer Binder for Soil Stabilization: Optimization Using D-Optimal Mixture Design

: In order to improve the comprehensive utilization rate of industrial solid waste and the road quality, a novel low-carbon and environmental friendly soil stabilizer is proposed. In this study, steel slag (SS), carbide slag (CS), blast furnace slag (BFS), fly ash (FA), and desulfurized gypsum (DG) were used as raw materials to develop a multiple industrial solid waste-based soil stabilizer (MSWSS). The optimal mix ratio of the raw materials determined by D-optimal design was as follows: 5% SS, 50% CS, 15% BFS, 15% DG, and 15% FA. The 7-day unconfined compressive strength (UCS) of MSWSS-stabilized soil was 1.7 MPa, which was 36% higher than stabilization with ordinary portland cement (OPC) and met the construction requirements of highways. After 7 days of curing, the UCS of MSWSS-stabilized soil was significantly higher than that in the OPC group. X-ray powder diffraction (XRD), thermogravimetric analysis (TGA), and scanning electron microscopy (SEM) analysis indicated that the prominent hydration products were ettringite (AFt) and C-S-H gel. The results showed that an amount of AFt and C-S-H were formed in the initial stage of curing, resulting in a rapid improvement in early UCS. As the curing proceeded, the content of AFt and C-S-H increased constantly and grew intertwined with each other, which lead to the denser microstructure of stabilized soil and better mechanical strength.


Introduction
Rapid urbanization has led to an exponential increase in the amount of solid waste that is generated.According to the report of the World Bank Group, global cities generate 2.01 billion tons of solid waste each year, including industrial waste, agricultural waste, and domestic waste.It is expected to increase to 3.40 billion tons per year by 2050 [1].According to the statistics, only about 13.5% of solid waste is recycled globally [1].Landfills and incineration are the two main solid waste disposal methods.However, these conventional treatments waste land resources, release harmful gases, and cause a negative impact on the economy, environment, and human health [2][3][4].
On the other hand, the development of urban transportation has brought challenges to the road construction industry [5,6].The quality of pavement depends on the mechanical properties of the underlying soil, including strength, compressibility, permeability, etc. [7][8][9].One of the effective ways to improve road quality is to use soil stabilizers.Traditional soil stabilizers are cement, lime, and coal fly ash, which are characterized by lower costs and being less time consuming [8,10,11].However, although the soil solidified using cement exhibits high Processes 2024, 12, 436 2 of 15 early strength, its later strength development is limited.In contrast, soil solidified with lime and coal fly ash shows lower early compressive strength [12,13].Additionally, the production process of these materials generates a large amount of carbon dioxide, accelerating global warming.The carbon dioxide emissions from cement production account for approximately 8% of global carbon dioxide emissions [14][15][16][17].Therefore, environmentally friendly, efficient, and sustainable soil stabilizers should be explored as alternatives for pavement construction.
Recently, non-traditional stabilizers, especially cementitious materials based on multiple industrial solid wastes, have received more attention due to their low carbon emissions, low cost, low energy consumption, and better sustainability [18][19][20][21][22].
Steel slag (SS) and blast furnace slag (BFS) are the industrial solid waste generated in the process of steelmaking, and their major components are CaO, SiO 2 , and Al 2 O 3 [23][24][25].In 2021, the SS output in China exceeded 120 million tons, but its comprehensive utilization was only about 20% [26].Carbide slag (CS) is the predominant by-product resulting from the hydrolysis of calcium carbide during acetylene production, which is mainly composed of CaOH [23].It is predicted that the annual production of CS exceeds 50 million tons in China, with a utilization rate is less than 30% [27].Desulfurized gypsum (DG) is an industrial waste generated from smelters, large-scale enterprise boilers, and the desulfurization of flue gas in thermal power plants [28].The output of DG in 2019 was approximately 130 million tons [29].The main component of DG is CaSO 4 •2H 2 O, which also contains impurities such as CaCO 3 and CaSO 3 .Fly ash (FA) is a major solid waste emitted from thermal power plants [30,31].The global annual production is estimated at 360 million tons [32], while China's annual output is about 800-900 million tons [33].FA contains many toxic elements including Pb, As, Se, and Cr, and is considered to be a hazardous waste to the environment and human health [34,35].
From a chemical composition point of view, these industrial solid wastes (including SS, BFS, CS, DG, and FA) exhibit higher pozzolanic activity and alkalinity and can be used as raw materials for the preparation of geopolymers [21,36].This is especially the case with CS, which, as an alkaline industrial solid waste, has the same main components and can be considered as a substitute for lime [37].A geopolymer can interact with soil particles to form a cementitious gel through a pozzolanic reaction.The addition of CS, as an alkaline resource, provides OH − and Ca 2+ .Under alkaline conditions, the dissolution of aluminum and silica mineral components can be activated and react with Ca(OH) 2 to generate C-S-H, which can reduce pores and increase the compressive strength and stiffness of soil [38].Several studies have been conducted on the utilization of industrial solid waste as soil stabilizers [19,29,[39][40][41][42][43].However, few studies investigated the action mechanisms of multiple industrial solid wastes as soil solidifiers, e.g., five different industrial solid wastes.The advantages of various solid wastes can be complemented to improve the utilization rate effectively and reduce costs.
In the present study, a novel low-carbon and environmentally friendly soil stabilizer is proposed in order to improve the comprehensive utilization efficiency of industrial solid waste and reduce environmental pollution.Steel slag (SS), carbide slag (CS), blast furnace slag (BFS), fly ash (FA), and desulfurized gypsum (DG) were used as raw materials to formulate a multiple industrial solid waste-based soil stabilizer (MSWSS).The D-optimal mixture design approach was utilized to determine the optimal formulation of the various raw materials for the MSWSS, which can reduce the number of experiments and costs [44].The mechanical characteristics of the stabilized soil were evaluated using the unconfined compressive strength (UCS) test and compared with conventional soil stabilizers.The micromorphology of solidified soil and hydration mechanisms were investigated by using X-ray diffraction (XRD), scanning electron microscopy (SEM), and thermogravimetric analysis (TGA).

Materials
The test soil used in this study was collected from Cangzhou, Hebei Province, China, and had such characteristics as a small plasticity index, poor stability when exposed to water, low strength, and the inability to be directly used as base for the construction of pavement.The characteristics of the test soil are provided in Table 1.The raw materials used to produce the MSWSS included CS, BFS, SS, FA, and DG, which were collected from the following companies: Hebei Jinniu Chemical Co., Ltd., Xingtai, China; Cangzhou Zhongtie Equipment Manufacturing Material Co., Ltd., Cangzhou, China; and Shenhua Group Guohua Power Plant, Yongzhou, China.The chemical compositions of the test soil and the raw materials are listed in Table 2, while the XRD patterns of the raw materials are shown in Figure 1.It can be seen from the results that SiO 2 (52.93%) comprised the largest proportion of the test soil, and the main phases were quartz, plagioclase, muscovite, and calcite.SS, CS, and FA had a high content of CaO (39.92%, 93.83%, 41.82%), while BFS had a high content of SiO 2 (39.83%).Additionally, SS contained 27.26% Fe 2 O 3 .The major components of DG were CaO (47.43%) and SO 3 (43.32%).As shown in Figure 1, the main phases of FA were quartz, mullite, CaO, and corundum.BFS is an amorphous phase structure.The dominant phases of SS were tetracalcium aluminoferrite C 4 AF (4CaO

Materials
The test soil used in this study was collected from Cangzhou, Hebei Province, China, and had such characteristics as a small plasticity index, poor stability when exposed to water, low strength, and the inability to be directly used as base for the construction of pavement.The characteristics of the test soil are provided in Table 1.The raw materials used to produce the MSWSS included CS, BFS, SS, FA, and DG, which were collected from the following companies: Hebei Jinniu Chemical Co., Ltd., Xingtai, China; Cangzhou Zhongtie Equipment Manufacturing Material Co., Ltd., Cangzhou, China; and Shenhua Group Guohua Power Plant, Yongzhou, China.The chemical compositions of the test soil and the raw materials are listed in Table 2, while the XRD patterns of the raw materials are shown in Figure 1.It can be seen from the results that SiO2 (52.93%) comprised the largest proportion of the test soil, and the main phases were quartz, plagioclase, muscovite, and calcite.SS, CS, and FA had a high content of CaO (39.92%, 93.83%, 41.82%), while BFS had a high content of SiO2 (39.83%).Additionally, SS contained 27.26% Fe2O3.The major components of DG were CaO (47.43%) and SO3 (43.32%).As shown in Figure 1, the main phases of FA were quartz, mullite, CaO, and corundum.BFS is an amorphous phase structure.The dominant phases of SS were tetracalcium aluminoferrite C4AF (4CaO•Al2O3•Fe2O3), dicalcium silicate C2S (2CaO•SiO2), and tricalcium silicate C3S (3CaO•SiO2).

Sample Preparation
The test soil was dried at a temperature of 55 • C and ground to a specific surface area reached 450-500 cm 2 /g.The MSWSS was produced by mixing CS, BFS, SS, FA, and DG in a certain proportion obtained using the D-optimal mixture design approach.The amount of the addition of MSWSS was fixed at 10% of the total mass of the test soil.Then, a certain amount of water was added, the MSWSS and test soil were mixed with the water, and the mixture moisture content was 17.86%.Afterwards, the samples were prepared in accordance with the "Test Methods of Materials Stabilized with Inorganic Binders for Highway Engineering" (JTGE51-2009).The specimens were placed within a cylindrical mold with a diameter and height of 50mm and shaped by applying a stable pressure using a TYA-2000 digital display pressure testing machine at a loading rate of 1 mm/min for 2 min.Afterwards, the samples were weighed, enclosed within sealed bags, and moved to a standard environment with a relative humidity of over 95% and a temperature maintained at 20 ± 2 • C for the specified curing time.After reaching the set time, the specimens were tested for compressive strength and microstructure characterization.

D-Optimal Mixture Design
The D-optimal mixture design is commonly used in experimental design to optimize multiple variables simultaneously.In this study, the optimal formulation of raw materials was determined using D-optimal design.The steel slag (SS), carbide slag (CS), blast furnace slag (BFS), desulfurized gypsum (DG), and fly ash (FA) were labelled as A, B, C, D, and E. The lower and upper limits of each component were proposed based on preliminary experiments: A (5-20 wt.%), B (20-50 wt.%), C (5-15 wt.%), D (5-20 wt.%), and E (10-25 wt.%).A 35-run D-optimal mixture design with multiple constraints on the component proportions was implemented by using Design-Expert software version 13.0 (Table 3).The combination of the components and the experimental results of 7d UCS are listed in Table 4.
Table 3.The variables and their constraints for D-optimal mixture design.

Compressive Strength Test
The unconfined compressive strengths (UCS) were tested according to Chinese standard JTG E51-2009 [45] using an electro-hydraulic universal testing machine (CTS-E200) with a loading rate of 1mm/min.The measurements were performed on at least three representative samples to obtain an average value.

X-ray Powder Diffraction (XRD)
The phase compositions of the samples were detected using the X-ray powder diffraction (XRD) method using a Bruker Advance D8 X-ray diffractometer with Cu Kα radiation (λ = 0.1542 nm) operating at 20 kV and 200 mA.The scanning angle ranged from 5 • to 80 • , with a scanning speed of 5 • /min and a step size of 0.02 • .

Thermogravimetric Analysis (TGA)
Thermal analyses were performed on a thermogravimetric differential scanning calorimeter (TG-DSC, NETZSCH STA-409, Selb, Germany) with a nitrogen protective gas.A 30 mg sample was used for each measurement, with a heating rate of 10 • C/min in the temperature range of 25 to 1000 • C.

Scanning Electron Microscopy (SEM)
Microstructural analysis of the specimens was performed using scanning electron microscopy (SEM) alongside an energy dispersive spectra (EDS) analyzer.EDS was utilized to determine the elemental composition of hydration products.SEM-EDS was conducted using a Zeiss SUPRA 55 field emission-scanning electron microscope (FE-SEM) outfitted with a LinkNA1000 energy spectrometer.

Fitting Model and Response Surface Analysis
The interaction between the five components was mathematically modeled against the response to 7d UCS.A comparison and evaluation of different fitting models was performed using the Design-Expert 13 program.The assessment results of the accuracy of the four fitting models are shown in Table 5.The model selection criteria were a higher F-value, the adjusted and predicted R 2 value, and a low p-value (<0.05).According to the result, the special cubic model was selected.The p-value of the special cubic model was below 0.0001, suggesting that the model adequately describes the response.With a coefficient of determination of R 2 = 0.9971, it can be inferred that the experimental data closely align with the predicted values, affirming the model's high significance.Analyses of variance (ANOVAs) were conducted on the response (7d UCS) using the Design-Expert software.The results of the ANOVAs for 7d UCS are shown in Table 6.The p-values and F-values represent the significance of variable factors.In general, smaller p -values (p-value > 0.05 insignificant, p-value < 0.05 significant, p-value < 0.01 extremely significant), along with larger F-values, denote the more significant impact of the variable factors on the response [46].It can be seen from Table 6 that except for the interaction terms CD, CE, DE, ACD, and BCE, the p-values of other interaction terms were all < 0.01, which means that they had a notable influence on the 7d UCS.In contrast, interactive items CD, CE, DE, ACD, and BCE had no significant effects on 7d UCS because their p-values were greater than 0.05.Multiple quadratic regression equations were used to fit the functional relationship between the influencing variable factors and the responses.Equation (1) reflects the mathematical relations between 7d UCS and factors A, B, C, D, and E. The coefficients of the equation could be used to describe the effect of these variables on the response (7d UCS).A positive coefficient in the equation signifies a positive effect, while a negative coefficient indicates a negative effect.According to the results, the interactive terms AB, AD, AE, BC, BE, ABC, ABD, ABE, ADE, BCD, and BDE had an extremely significant effect (p-value < 0.01) on the 7d UCS.Table 7 presents the findings of the model reliability test analysis.R 2 represents the correlation coefficient between the predicted and actual values.A higher R 2 value indicates a stronger correlation.If the values of the adjusted R 2 and predicted R 2 are high and the difference between them is close (adjusted R 2 -predicted R 2 < 0.2), then the regression model can further prove this conclusion.As shown in Table 7, the values of R 2 , adjusted R 2 , predicted R 2 , and adequate precision were 0.9971, 0.9893, 0.8363, and 43.2860, respectively.This indicated that the model equation can represent the relationship between the components and the 7d UCS well.The value of C.V. % was 1.48, indicating that the experiment shows high reliability and accuracy.The results showed that this model was able to predict the 7d UCS relatively accurately.The response surface interaction analysis was used to evaluate the effect of each raw material on the 7d UCS of the solidified soil.Figure 2 shows the 3D surface graphs and contour illustrating the impact of the independent variables on the 7d UCS.Generally, when the response surface shows a curved form, it indicates that there is a degree of interaction among the raw material components [46].As seen in Figure 2, an increase in the dosage of desulfurized gypsum and blast furnace slag resulted in an increase in the 7d UCS, whereas an increase in the dosage of steel slag, carbide slag, and fly ash had the opposite effect.

Optimization of Mixture Composition
The maximum unconfined compressive strength value for 7 days (7d UCS) is used as the target optimization value.Based on regression model analysis combined with the 3D response surface plots and contour plots, the optimal ratio of the soil solidifier can be obtained.According to the predicted result, the optimal raw material proportion for preparing the MSWSS was: 5% steel slag, 50% carbide slag, 15% blast furnace slag, 15% desulfurized gypsum and 15% fly ash, respectively.The predicted value of the 7d UCS was 1.65 MPa.To confirm the precision of the predicted value, experiments on the compressive strength for 7 days were carried out.The outcome is displayed in Table 8.The experimental value closely aligned with the predicted value, which demonstrated the reliability of the model.The value of the absolute relative deviations was 1.81% (<5%), suggesting that the prediction model exhibits high accuracy.

Optimization of Mixture Composition
The maximum unconfined compressive strength value for 7 days (7d UCS) is used as the target optimization value.Based on regression model analysis combined with the 3D response surface plots and contour plots, the optimal ratio of the soil solidifier can be obtained.According to the predicted result, the optimal raw material proportion for

Compressive Strength Test
MSWSS was prepared from multiple solid waste according to the optimal ratio obtained using the D-optimal mixture design method (SS:CS:BFS:DG:FA = 5:50:15:15:15).To compare the solidification effects of multiple solid waste-based stabilizer and ordinary Processes 2024, 12, 436 9 of 15 portland cement (OPC) on the test soil, the unconfined compressive strength (UCS) of these two stabilized soils was tested and analyzed.Standard specimens were prepared using the test soil plus 10% MSWSS, with a moisture content of 17.65%.In addition, 10% P.O32.5 OPC was used instead of MSWSS as a comparative experiment.
The unconfined compressive strength (UCS) of the stabilized soil at various curing ages is presented in Figure 3.The UCS of the stabilized soil showed a progressive increase with the curing time.The soil solidified using MSWSS has a UCS ranging from 0.35 MPa to 3.27 MPa.At the same time, soil solidified using OPC has a UCS ranging from 0.47 MPa to 2.89 MPa.The 7d UCS results of the MSWSS group and of the OPC group were 1.7 MPa and 1.25 MPa, respectively, which met the standard requirements of Chinese technical standards for the Technical Guidelines for Construction of Highway Roadbases.The 7d UCS of MSWSS-stabilized soil is 36% higher than that of the OPC group.It can be seen from Figure 3 that the UCS of OPC stabilized soil was higher than that of MSWSS-stabilized soil when the curing age is in the first 3 days.Afterwards, with the prolongation of curing time, the UCS of soil stabilized using MSWSS was significantly higher than that of the OPC group.After curing for 180 days, the compressive strength exhibited a notable increase to 3.27 MPa and 2.89 MPa, respectively.
MSWSS was prepared from multiple solid waste according to the optimal r tained using the D-optimal mixture design method (SS:CS:BFS:DG:FA = 5:50:15:15 compare the solidification effects of multiple solid waste-based stabilizer and o portland cement (OPC) on the test soil, the unconfined compressive strength ( these two stabilized soils was tested and analyzed.Standard specimens were p using the test soil plus 10% MSWSS, with a moisture content of 17.65%.In additi P.O32.5 OPC was used instead of MSWSS as a comparative experiment. The unconfined compressive strength (UCS) of the stabilized soil at variou ages is presented in Figure 3.The UCS of the stabilized soil showed a progressive with the curing time.The soil solidified using MSWSS has a UCS ranging from 0 to 3.27 MPa.At the same time, soil solidified using OPC has a UCS ranging from 0 to 2.89 MPa.The 7d UCS results of the MSWSS group and of the OPC group were and 1.25 MPa, respectively, which met the standard requirements of Chinese t standards for the Technical Guidelines for Construction of Highway Roadbases UCS of MSWSS-stabilized soil is 36% higher than that of the OPC group.It can from Figure 3 that the UCS of OPC stabilized soil was higher than that of MSWS lized soil when the curing age is in the first 3 days.Afterwards, with the prolong curing time, the UCS of soil stabilized using MSWSS was significantly higher tha the OPC group.After curing for 180 days, the compressive strength exhibited a increase to 3.27 MPa and 2.89 MPa, respectively.

XRD Analysis
The XRD patterns of the solidified test soil at curing times of 7, 28, and 90 days are presented in Figure 4, and the primary phase compositions have been identified and labeled.Based on the results, the main phases were quartz, calcite, and plagioclase, and the primary hydration products were ettringite (AFt) and calcium silicate hydrate (C-S-H).This is consistent with the results of the SEM analysis (Figure 5).In addition, it also contained small amounts of chlorite and microcline.After curing for 7 days, ettringite and C-S-H can be observed from XRD pattern.As the curing time progressed, the intensity of AFt increased gradually, and the intensity of C-S-H also showed a weak increase, indicating that the hydration reaction was constantly developing.Meanwhile, a small amount of C 2 S was detected, and the peak intensity decreased from 7 to 28 days of curing time, consistent with previous research [4,47].Moreover, no distinct peaks of C 2 S were observed after 90 days of curing.This showed that that C 2 S is involved in the hydration process.Generally, it was assumed that the hydration reaction of C 2 S generated C-S-H.
AFt increased gradually, and the intensity of C-S-H also showed a weak increase, indicating that the hydration reaction was constantly developing.Meanwhile, a small amount of C2S was detected, and the peak intensity decreased from 7 to 28 days of curing time, consistent with previous research [4,47].Moreover, no distinct peaks of C2S were observed after 90 days of curing.This showed that that C2S is involved in the hydration process.Generally, it was assumed that the hydration reaction of C2S generated C-S-H.ing that the hydration reaction was constantly developing.Meanwhile, a small amount of C2S was detected, and the peak intensity decreased from 7 to 28 days of curing time, consistent with previous research [4,47].Moreover, no distinct peaks of C2S were observed after 90 days of curing.This showed that that C2S is involved in the hydration process.Generally, it was assumed that the hydration reaction of C2S generated C-S-H.

TG-DSC Analysis
The TG-DSC curves of the samples after curing for 7 and 28 days are displayed in Figure 5.The endothermic peaks at 50-400 • C could be primarily attributed to the thermal decomposition of C-S-H gel and AFt [24].The endothermic peaks at 88 • C and 150 • C of the DSC curves were due to the release of free water [48].The mass losses at 88 • C and 150 • C of the stabilized soil at 7 days and 28 days were 0.49% and 1.33%, and 0.76% and 1.78%, respectively.The endothermic peak at 362 • C of the DSC curves at a curing age of 28 days can be attributed to the dehydration of the crystal water of C-S-H and ettringite.The endothermic peaks at 680 • C and 800 • C were due to the decarburization and decomposition of calcite [25].The results of the thermal analysis showed that the total mass loss of the samples after curing for 7 days and 28 days was 11.32% and 12.86%, respectively.This indicated that the hydration reaction continues as the curing age increases.The XRD analysis results also showed that the aggregate mass of the hydration products increases with the increase in curing time.

SEM-EDS Analysis
Figure 6 illustrates the microstructure of the solidified soils after 7 and 28 days of hydration.Figure 6a,b represent the samples at 7 curing days and 28 curing days.After curing for 7days, the sample presented a loose arrangement and showed the typical needlerod crystal structures of AFt (Figure 6a).Simultaneously, C-S-H gel was detected, which covered the needle-rod crystals of the AFt.The AFt crystals and C-S-H grew intertwined with each other, and this kind of structure could effectively improve the early compressive strength of the stabilized soil (Figure 3).

SEM-EDS Analysis
Figure 6 illustrates the microstructure of the solidified soils after 7 and 28 days of hydration.Figure 6a,b represent the samples at 7 curing days and 28 curing days.After curing for 7days, the sample presented a loose arrangement and showed the typical needle-rod crystal structures of AFt (Figure 6a).Simultaneously, C-S-H gel was detected, which covered the needle-rod crystals of the AFt.The AFt crystals and C-S-H grew intertwined with each other, and this kind of structure could effectively improve the early compressive strength of the stabilized soil (Figure 3).
As the curing time progressed to 28 days (Figure 6b), the content of the needle-rod structures representing ettringite was increased, which corresponds with the findings of the XRD analysis.Meanwhile, the loosely structured solidified soil became denser, which might be attributed to the gradually increasing output of the C-S-H gel filling the pore space.At this stage, the AFt is surrounded by C-S-H gel, resulting in a significant enhancement in the compressive strength.

Discussion
According to characterization of the solidified soil, the primary hydration products were C-S-H gel (Ca5Si6O16(OH)•4H2O) and AFt (3CaO•Al2O3•3CaSO4•32H2O). Carbide slag (CS) plays a crucial role in the hydration process because it provides Ca 2+ ions and an alkaline environment.At the initial phase of hydration, the C2S (2CaO•SiO2) and C3S (3CaO•SiO2) in SS were hydrolyzed to form C-S-H gel and calcium hydroxide (Ca(OH)2) (Equations ( 2) and (3)), which improved the early compressive strength (Figure 3).Meanwhile, a notable quantity of Ca(OH)2 was produced, raising the alkalinity of the binder.As the curing time progressed to 28 days (Figure 6b), the content of the needle-rod structures representing ettringite was increased, which corresponds with the findings of the XRD analysis.Meanwhile, the loosely structured solidified soil became denser, which might be attributed to the gradually increasing output of the C-S-H gel filling the pore space.At this stage, the AFt is surrounded by C-S-H gel, resulting in a significant enhancement in the compressive strength.

Discussion
According to characterization of the solidified soil, the primary hydration products were C-S-H gel (Ca 5 Si 6 O 16 (OH)•4H 2 O) and AFt (3CaO•Al 2 O 3 •3CaSO 4 •32H 2 O).Carbide slag (CS) plays a crucial role in the hydration process because it provides Ca 2+ ions and an alkaline environment.At the initial phase of hydration, the C 2 S (2CaO•SiO 2 ) and C 3 S (3CaO•SiO 2 ) in SS were hydrolyzed to form C-S-H gel and calcium hydroxide (Ca(OH) 2 ) (Equations ( 2) and (3)), which improved the early compressive strength (Figure 3).Meanwhile, a notable quantity of Ca(OH) 2 was produced, raising the alkalinity of the binder.4).The chemical reaction can be expressed in Equations ( 4)- (7).Consequently, large amounts of C-S-H and AFt grew interlacedly, filling the gaps between soil particles and improving the compactness and strength of the solidified soil.The hydration mechanism of the MSWSS is shown in Figure 7.
tions ( 4)- (7).Consequently, large amounts of C-S-H and AFt grew interlacedly, filling the gaps between soil particles and improving the compactness and strength of the solidified soil.The hydration mechanism of the MSWSS is shown in Figure 7.

Conclusions
This study proposed the preparation of a multiple industrial solid waste-based soil stabilizer (MSWSS) using steel slag (SS), carbide slag (CS), blast furnace slag (BFS), fly ash (FA), and desulfurized gypsum (DG).The D-optimal mixture approach was utilized in the experimental design to obtain the optimal formulation of raw materials.The unconfined compressive strength (UCS) of the solidified test soil was evaluated, and the hydration mechanism was analyzed using the XRD, TG-DSC, and SEM methods.In summary, the following conclusions can be drawn: (1) Interaction between the five components was modeled mathematically against the responses of 7d UCS.According to the result of the ANOVA and the model reliability analysis, the prediction model showed a lower p-value (<0.0001) and lower absolute

Conclusions
This study proposed the preparation of a multiple industrial solid waste-based soil stabilizer (MSWSS) using steel slag (SS), carbide slag (CS), blast furnace slag (BFS), fly ash (FA), and desulfurized gypsum (DG).The D-optimal mixture approach was utilized in the experimental design to obtain the optimal formulation of raw materials.The unconfined compressive strength (UCS) of the solidified test soil was evaluated, and the hydration mechanism was analyzed using the XRD, TG-DSC, and SEM methods.In summary, the following conclusions can be drawn: (1) Interaction between the five components was modeled mathematically against the responses of 7d UCS.According to the result of the ANOVA and the model reliability analysis, the prediction model showed a lower p-value (<0.0001) and lower absolute relative deviations (1.81%), and a higher F-value (406.92),R2 value (0.9971), and adjusted R2 value (0.9893).This indicated that the model had high reliability and was able to predict 7d UCS relatively accurately.(2) The D-optimal mixture approach indicated that the optimal ratio of MSWSS should be 5% steel slag (SS), 50% carbide slag (CS), 15% blast furnace slag (BFS), 15% desulfurized gypsum (DG), and 15% fly ash (FA).The UCS of the stabilized soil increased with the curing time.The 7d UCS of soil stabilized using the multiple industrial solid waste-based soil stabilizer was 1.7 MPa and the 180d UCS increased to 3.27 MPa, which all met the requirements for the "Technical Guidelines for Construction of Highway Roadbases".The UCS of the MSWSS stabilized test soil group was significantly higher than that of the OPC group.(3) XRD, TG-DSC, and SEM results revealed that the primary hydration products were AFt and C-S-H gel.As the curing progressed, the content of the AFt and C-S-H increased and the microstructure of stabilized soil exhibited a denser structure, leading to better mechanical strength.
The research proposed a novel low-carbon and environmentally friendly soil stabilizer, MSWSS, based on multiple industrial solid waste products, which could achieve better mechanical properties than traditional OPC stabilizers and meet the requirements of highway

Figure 1 .
Figure 1.XRD patterns of the test soil and the raw materials.Figure 1. XRD patterns of the test soil and the raw materials.

Figure 1 .
Figure 1.XRD patterns of the test soil and the raw materials.Figure 1. XRD patterns of the test soil and the raw materials.

Figure 2 .
Figure 2. 3D surface graphs and contour plots for the effects of the independent variables on 7d UCS (A: Steel slag, B: Carbide slag, C: Blast furnace slag, D: Desulfurized gypsum, E: Fly ash): (a) the effect of A, D, E on 7d UCS; (b) the effect of A, C, E on 7d UCS; (c) the effect of A, C, D on 7d UCS; (d) the effect of B, C, E on 7d UCS.

Figure 2 .
Figure 2. 3D surface graphs and contour plots for the effects of the independent variables on 7d UCS (A: Steel slag, B: Carbide slag, C: Blast furnace slag, D: Desulfurized gypsum, E: Fly ash): (a) the effect of A, D, E on 7d UCS; (b) the effect of A, C, E on 7d UCS; (c) the effect of A, C, D on 7d UCS; (d) the effect of B, C, E on 7d UCS.

Figure 3 .
Figure 3. Compressive strength of stabilized soil at different curing times.

Figure 3 .
Figure 3. Compressive strength of stabilized soil at different curing times.

Figure 5 .
Figure 5. TG-DSC diagram of solidified soil at different curing times.

Figure 5 .
Figure 5. TG-DSC diagram of solidified soil at different curing times.Figure 5. TG-DSC diagram of solidified soil at different curing times.

Figure 5 .
Figure 5. TG-DSC diagram of solidified soil at different curing times.Figure 5. TG-DSC diagram of solidified soil at different curing times.

Figure 6 .
Figure 6.SEM morphology of the solidified soil for (a) 7 days and (b) 28 days curing times.

Figure 6 .
Figure 6.SEM morphology of the solidified soil for (a) 7 days and (b) 28 days curing times.
In alkaline environments, OH − reacts with the amphoteric oxide Al 2 O 3 and acidic oxide SiO 2 which can accelerate the hydrolysis of the glass phase, destroying the Ca-O, Si-O, and Al-O bonds and form [Al(OH) 6 ] 3− and [Si(OH) 3 ] − [49].Subsequently, with the participation of the Ca 2+ and SO 4 2+ produced by the hydrolysis of DG, it reacts with [Al(OH) 6 ] 3− and [Si(OH) 3 ] − to generate AFt.At the same time, [Si(OH) 3 ] − could react with Ca 2+ to generate C-S-H gel.With the progression of the hydration reaction, the content of C-S-H and AFt increased significantly (Figure

Figure 7 .
Figure 7. Hydration mechanism of the multiple industrial solid waste-based soil stabilizer (MSWSS).

Table 1 .
Physical and mechanical properties of the test soil.

Table 2 .
Chemical composition of the test soil and the raw materials (wt.%).

Table 1 .
Physical and mechanical properties of the test soil.

Table 2 .
Chemical composition of the test soil and the raw materials (wt.%).

Table 4 .
Experimental design of the five components and the results for 7-day unconfined compressive strength (7d UCS).

Table 5 .
Accuracy evaluation of different fitting models.

Table 6 .
ANOVA table of experimental results.

Table 7 .
Model reliability test analysis.

Table 8 .
Comparison of the observed values and the predicted values after optimization.