Low-CO2 Optimization Design of Quaternary Binder Containing Calcined Clay, Slag, and Limestone

Blended cement is commonly used for producing sustainable concretes. This paper presents an experimental study and an optimization design of a low-CO2 quaternary binder containing calcined clay, slag, and limestone using the response surface method. First, a Box–Behnken design with three influencing factors and three levels was used for the combination design of the quaternary composite cement. The lower limit of the mineral admixtures was 0%. The upper limits of slag, calcined clay, and limestone powder were 30%, 20%, and 10%, respectively. The water-to-binder ratio (water/binder) was 0.5. Experimental works to examine workability and strength (at 3 and 28 days) were performed for the composite cement. The CO2 emissions were calculated considering binder compositions. A second-order polynomial regression was used to evaluate the experimental results. In addition, a low-CO2 optimization design was conducted for the composite cement using a composite desirability function. The objectives of the optimization design were the target 28-day strength (30, 35, 40, and 45 MPa), target workability (160 mm flow), and low CO2 emissions. The trends of the properties of optimal combinations were consistent with those in the test results. In summary, the proposed optimization design can be used for designing composite cement considering strength, workability, and ecological aspects.


Introduction
Clay and limestone are minerals available in many countries worldwide, and granulated blast furnace slag is a byproduct of the steel-making industry.Calcined clay, slag, and limestone are being increasingly used as mineral admixtures in concrete production.Using these mineral admixtures, numerous advantageous properties can be achieved, such as good workability, low hydration heat, and low environmental impact.To rationally use calcined clay, slag, and limestone, both experimental studies and numerical optimization studies of the workability, mechanical, and ecological aspects are necessary [1,2].
Many experimental studies have been conducted on binary or ternary blended concretes containing calcined clay, slag, and limestone.First, regarding binary blended concretes, Kang et al. [3] found that cement hydration can be accelerated by the inclusion of limestone powder.Li et al. [4] reported that the threshold value of chloride ions decreases as the content of limestone increases.Tironi et al. [5] showed that calcination and grinding improve the reaction of raw clay.Lin et al. [6] found that calcined hwangtoh clay has nucleation, dilution, and chemical effects on the hydration reaction of ordinary Portland cement.Kocaba et al. [7] reported that calorimetry calibration and backscattered electron imaging are promising methods to measure the reaction degree of slag.Second, for ternary hybrid concrete, Arora et al. [8] found a synergistic effect between limestone and slag and determined the linear relationship between carbonate consumption and carboaluminate formation.Ramezanianpour and Hooton [9] found that a high content of alumina in slag or metakaolin increases the reaction rate of limestone and its optimum level.Dhandapani and Santhanam [10] showed that cement-based material with a limestone-calcined clay-cement (LC3) binder has much lower permeability than a fly ash hybrid binary concrete.Dhandapani et al. [11] presented results showing that LC3 concretes are suitable for structures in a chloride attack environment.Pillai et al. [12] reported that LC3 concrete has a much lower CO 2 footprint than plain concrete while having similar strength.In summary, previous experimental studies have mainly focused on the properties of ternary or binary blended concretes; there have been few studies on quaternary blended concretes incorporating calcined clay, slag, and limestone filler [13,14].In some countries, such as China, Korea, and India, many mineral admixtures are in use, such as slag, calcined clay, and limestone.Concrete producers and researchers are interested in whether slag, calcined clay, and limestone can be used together for producing quaternary blended cement.In addition, what about the performance of such quaternary-blended-cement-based materials?Are there possible synergistic effects among calcined clay, slag, and limestone?
Several numerical models for blended cement have been proposed to evaluate the properties of blended concretes.First, regarding binary blended concretes, Kolani et al. [15] modeled the kinetic process of hydration of slag-cement binary blends and evaluated the heat of hydration, combined water, and calcium hydroxide.Wang [16] proposed an integrated model for metakaolin-cement hybrid blends and evaluated the degree of reaction of metakaolin and cement, strength, and permeability properties of metakaolin hybrid concretes.Wang [17] evaluated the degree of hydration, strength development, and carbonation durability of limestone-blended concrete considering the nucleation effect and dilution effect of limestone filler.Second, regarding ternary blended concretes, Kunther et al. [18] conducted thermodynamic modeling and evaluated the compositions of reaction products of ternary blends of cement-metakaolin-limestone.Wang and Luan [19] proposed a kinetic hydration and strength model for ternary blends of cement-slag-limestone and determined the optimum combination ratio of slag and limestone for different ages.Carrasco et al. [20] proposed experimental design methods and produced market-oriented cement containing slag and limestone.Yang et al. [21] simulated chloride penetration through LC3 highperformance concretes considering the binder composition and the damage.Avet et al. [22] determined the reacted extent of metakaolin in calcined clay and evaluated the phase assemblage of an LC3 paste by thermodynamic modeling and mass balance.In summary, the current numerical processes are mainly valid for plain, binary, or ternary hybrid concretes, but models that are effective for quaternary blended concretes require further investigation.
Based on this literature study, we found that compared with fundamental studies on plain, binary, or ternary blended concretes, studies on quaternary concretes are limited.To fill this gap, this paper presents experimental studies and numerical modeling of quaternary cement containing slag, calcined clay, and limestone.Tests of the workability and strength were performed, and the CO 2 emissions were calculated.Furthermore, based on the test results, a numerical model was used to determine the optimum combinations with target flow and strength and low CO 2 emissions.
Compared with previous research, the scientific innovations of this article mainly include the following aspects: 1.The research object is four mixed cement-based materials.Previous studies were mostly on two and three mixed cement-based materials; 2. The response surface design method proposed in this paper systematically considers fluidity, strength, and CO 2 emissions, while CO 2 emissions were not considered in previous concrete mix design methods; 3. The design method proposed in this article has certain portability and is suitable for designing low-carbon concretes in different countries and regions.

Materials
Ordinary Portland cement of ASTM type I (ASTM C150) was obtained from Sung Shin Cement [23], Republic of Korea.Hwangtoh is a clay mineral widely distributed in Korea.It contains high levels of calcium potassium chloride and calcium.Clay was ground before calcination, and calcined clay was produced by calcining hwangtoh clay at 800 • C for 60 min [6].Figure 1 shows the XRD analysis results of clay before and after calcination.We can see that the calcination eliminated the kaolinite peak (2θ = 12.37 • ) [6].Limestone powder and granulated blast furnace slag were obtained from Daejung Chemical Co., Ltd.(Sinan, Republic of Korea) and Asia Cement Co., Ltd.(Seoul, Republic of Korea), respectively.Table 1 presents the chemical compositions of the cement, limestone, slag, and calcined clay.The standard deviation for the XRF analysis results in Table 1 is 0.01%.Calcined clay and slag have higher aluminum contents than cement, which can support a synergistic effect between the limestone and slag or calcined clay.Calcined clay presents a higher silicon content than cement, which can contribute to the production of secondary calcium silicate hydrate (CSH).The average particle sizes of the cement, slag, calcined clay, and limestone were 18.Previous studies were mostly on two and three mixed cement-based materials; 2. The response surface design method proposed in this paper systematically considers fluidity, strength, and CO2 emissions, while CO2 emissions were not considered in previous concrete mix design methods; 3. The design method proposed in this article has certain portability and is suitable for designing low-carbon concretes in different countries and regions.

Materials
Ordinary Portland cement of ASTM type I (ASTM C150) was obtained from Sung Shin Cement [23], Republic of Korea.Hwangtoh is a clay mineral widely distributed in Korea.It contains high levels of calcium potassium chloride and calcium.Clay was ground before calcination, and calcined clay was produced by calcining hwangtoh clay at 800 °C for 60 min [6].Figure 1 shows the XRD analysis results of clay before and after calcination.We can see that the calcination eliminated the kaolinite peak (2θ = 12.37°) [6].Limestone powder and granulated blast furnace slag were obtained from Daejung Chemical Co., Ltd.(Sinan, Republic of Korea) and Asia Cement Co., Ltd.(Seoul, Republic of Korea), respectively.Table 1 presents the chemical compositions of the cement, limestone, slag, and calcined clay.The standard deviation for the XRF analysis results in Table 1 is 0.01%.Calcined clay and slag have higher aluminum contents than cement, which can support a synergistic effect between the limestone and slag or calcined clay.Calcined clay presents a higher silicon content than cement, which can contribute to the production of secondary calcium silicate hydrate (CSH).The average particle sizes of the cement, slag, calcined clay, and limestone were 18.3 µm, 12.8 µm, 12.3 µm, and 7.4 µm, respectively.

Experimental Method
Mixture preparation: The cement paste was mixed using a mechanical mixer and immediately put into a steel mold.At the age of 1 day, the steel mold was removed, and after removal, the paste specimens were placed in a curing chamber for sealed curing until they reached test age.The temperature for sealed curing was 20 degrees Celsius.
Test: The flow properties of the slurries were measured based on the mini-slump test method.The dimensions of the mini-slump cone were a height of 50 mm, a top diameter of 70 mm, and a bottom diameter of 100 mm.In addition, at curing ages of 3 days and 28 days, the compressive strength of the paste samples was tested according to ASTM C349 [24].The sample size for compressive strength testing was 50 × 50 × 50 mm.Samples for compressive strength testing were cured via sealed curing.For each mixture, three specimens were prepared.The average of the measured values for the three specimens was taken as the strength value for each group of specimens.However, when the difference between the maximum or minimum value among the three measured values and the middle value exceeded 15% of the middle value, the maximum and minimum values were eliminated, and the middle value was taken as the compressive strength of the group of specimens.

Experimental Design
A Box-Behnken design (BBD) was employed to design the quaternary hybrid cement combinations [25,26].As can be seen from Table 2, the BBD consisted of three independent factors (limestone, calcined clay, and slag).For each factor, there were three values corresponding to codes +1, 0, and −1.The maximum replacement ratios (high-coded +1) of limestone, calcined clay, and slag were 10%, 20%, and 30%, respectively.The maximum replacement ratio by the mineral admixture was 60% (10% + 20% + 30% = 60%).Cases of control, binary blends, ternary blends, and quaternary blends were considered.The minimum ratio (low-coded −1) for limestone, calcined clay, and slag was 0%.The water-tobinder ratio of the specimens was 0.5.The combinations of components of the binders are shown in Figure 2 and listed in Table 3.A total of 21 specimens were prepared, consisting of 17 BBD specimens and 4 additional specimens.The additional specimens were added to consider the lower limit combination (M1), the upper limit combination (M21), and the synergetic effect between limestone and calcined clay (M5 and M14).CO2 emissions are an important indicator of sustainable concrete.Herein, quaternary cement was designed with CO2 emissions as the optimization target.According to the sample mixtures and the specific CO2 emissions of each binder (listed in Table 4) [27], the CO2 emissions per unit volume of each sample can be determined.CO 2 emissions are an important indicator of sustainable concrete.Herein, quaternary cement was designed with CO 2 emissions as the optimization target.According to the sample mixtures and the specific CO 2 emissions of each binder (listed in Table 4) [27], the CO 2 emissions per unit volume of each sample can be determined.The CO 2 emissions of 1 m 3 of composite paste can be determined as follows.First, we determine the individual masses of the paste components: where m i and ρ i are the mass and density (shown in Table 1), respectively, of component i of the composite binder.Because the composite binder consists of five constituents, i.e., cement, calcined clay, slag, limestone powder, and water, the value of subscript i ranges from 1 to 5. Second, from the individual masses and individual CO 2 emissions of the components of paste (shown in Table 4), we can calculate the CO 2 emissions of 1 m 3 of paste.
Here, CO 2i is the individual CO 2 emissions (shown in Table 4) of component i of the paste, and CO 2 is the total CO 2 emissions of 1 m 3 of paste.
A second-order polynomial model was employed to perform a regression of the results (flow, strength, and CO 2 emissions) and examine the effect of the type and content of the mineral admixture [26].The general equation of the second-order polynomial regression model is as follows: where y is an experimental result (such as flow, strength, or CO 2 emissions), x i is an independent factor (three independent factors were employed, i.e., the contents of limestone, calcined clay, and slag), β 0 is the intercept coefficient, β i is the linear term coefficient, β ii is the quadratic term coefficient, and β ij is the coefficient of the interaction term.For different experimental tests, the coefficients (β 0 , β i , β ii , and β ij ) may differ.

Experimental Results
The experimental results (strength, flow) and CO 2 emission calculation results are summarized in Table 5.Table 6 shows the regression coefficients of all the results.As shown in Table 6, the parameters of the equations consist of an intercept, linear term, quadratic terms, and interaction terms.Based on the results for flow, strength, and CO 2 emissions, the coefficients of the equations were calibrated.The proposed equations can be validated by examining the p-value, lack of fit, and R 2 of the regression.The coefficients of determination between the experimental and predicted values were higher than 0.95, which means that the proposed equations for regression of the results are reasonable.The p values of the regression equations were less than 0.0001, which suggests that these regression models are significant.In addition, the lack of fit for each result was not significant relative to the pure error.

Flow of Paste
As shown in Table 5, the overall trend of flow is that after adding SCM, the flow of the paste was reduced.The flow of M21 was the lowest, at 160 mm, mainly because M21 contains a high content of clay and slag.These two substances have high aluminum contents, which can reduce the flow.From the experimental results (shown in Figure 3a), the flow of a paste was regressed as follows: where x 1 , x 2 , and x 3 denote the contents of limestone, calcined clay, and slag, respectively.The predictions of flow agreed with the experimental results.Figure 3b shows the impact of each component on flow.Limestone has a negligible effect on flow.With increasing slag content, the flow rate decreases slightly.This is because of the higher surface area of slag than cement.However, as the content of calcined clay increases, the flow decreases significantly.This is because of the higher aluminum content in calcined clay than in cement.A higher aluminum content makes the mixture prone to forming ettringite, which can lower the flow [28].In addition, the high surface area and specific structure of calcined clays also contribute to a lower slump.Likewise, Lin et al. [6] found that the addition of calcined hwangtoh clay can significantly reduce the slump of concrete.x , and x denote the contents of limestone, calcined clay, and slag, respec- tively.The predictions of flow agreed with the experimental results.Figure 3b shows the impact of each component on flow.Limestone has a negligible effect on flow.With increasing slag content, the flow rate decreases slightly.This is because of the higher surface area of slag than cement.However, as the content of calcined clay increases, the flow decreases significantly.This is because of the higher aluminum content in calcined clay than in cement.A higher aluminum content makes the mixture prone to forming ettringite, which can lower the flow [28].In addition, the high surface area and specific structure of calcined clays also contribute to a lower slump.Likewise, Lin et al. [6] found that the addition of calcined hwangtoh clay can significantly reduce the slump of concrete.

Compressive Strength
Table 5 shows the overall trend of the 3-day strength (shown in Figure 4a): after adding SCM, the strength of the paste was reduced.This is mainly because of the dilution effect of SCM.M21 had the lowest 3-day strength, which is due to the higher SCM content in M21.In addition, M2 showed slightly higher strength than the control specimen M1, which is due to the limestone powder's nucleation effect, which accelerates the early-age hydration rate of Portland cement and improves its strength.From the experimental results, the 3-day strength of a paste was regressed as follows:

Compressive Strength
Table 5 shows the overall trend of the 3-day strength (shown in Figure 4a): after adding SCM, the strength of the paste was reduced.This is mainly because of the dilution effect of SCM.M21 had the lowest 3-day strength, which is due to the higher SCM content in M21.
In addition, M2 showed slightly higher strength than the control specimen M1, which is due to the limestone powder's nucleation effect, which accelerates the early-age hydration rate of Portland cement and improves its strength.From the experimental results, the 3-day strength of a paste was regressed as follows: hydration extent of cement, and can enhance the early-age strength of concrete.However, when the content of limestone powder further increases from 5% to 10% (code 0 to code +1), owing to the limestone's dilution effect, the reduction in the strength becomes remarkable.In addition, as displayed in Figure 4d, at the age of 28 days, the strength decreases with increasing mineral admixture.This is due to the dilution effect and the lower reactivity of the mineral admixture compared with that of the cement.It should be noticed that in Figure 4d, the code +1 for limestone and calcined clay represents 10% and 20%, respectively.In other words, for limestone and calcined clay, when the x-axis represents the replacement ratio, the regression curves will be different.Figure 4e,f show the normalized strengths at the ages of 3 and 28 days, respectively.The paste's normalized strength was calculated as the compressive strength ratio between hybrid paste and plain paste.The proportional reduction line corresponds to the dilution effect.At a curing age of 3 days, the normalized strength values are close to the dilution effect line, whereas, at the age of 28 days, they are much higher than the dilution effect line.This indicates the increasing reaction extent of the mineral admixtures from curing ages of 3 days to 28 days.

CO2 Emissions
Table 5 shows the overall trend of CO2 emissions: after adding SCM, CO2 emissions were reduced.This is mainly because of the dilution effect of SCM reducing the amount of cement used.M21 has the lowest CO2 emissions due to its higher SCM substitution.From the CO2 emission calculation results (shown in Figure 5a), the CO2 emissions of a paste were regressed based on a second-order polynomial as follows: As shown in Table 5, the overall trend of the 28-day strength was that after adding SCM, the strength of the paste was reduced.This is mainly because of the dilution effect of SCM.M21 had the lowest 28-day strength, which is due to the higher SCM content in M21.In addition, the strength of M2 was slightly lower than that of the standard specimen M1, which differs from the trend for the 3-day strength.This is because the dilution effect of limestone is the dominant factor in the later stage (28 days), while the nucleation effect is the dominant factor in the early stage (3 days).The 28-day strength of a paste (shown in Figure 4b) was determined as follows: As shown in Figure 4c, at the curing age of 3 days, when the limestone content ranges from 0% to 5% (code changes from −1 to 0), the decrease in the strength is not significant.This is because limestone has a nucleation effect on cement hydration, can improve the hydration extent of cement, and can enhance the early-age strength of concrete.However, when the content of limestone powder further increases from 5% to 10% (code 0 to code +1), owing to the limestone's dilution effect, the reduction in the strength becomes remarkable.In addition, as displayed in Figure 4d, at the age of 28 days, the strength decreases with increasing mineral admixture.This is due to the dilution effect and the lower reactivity of the mineral admixture compared with that of the cement.It should be noticed that in Figure 4d, the code +1 for limestone and calcined clay represents 10% and 20%, respectively.In other words, for limestone and calcined clay, when the x-axis represents the replacement ratio, the regression curves will be different.Figure 4e,f show the normalized strengths at the ages of 3 and 28 days, respectively.The paste's normalized strength was calculated as the compressive strength ratio between hybrid paste and plain paste.The proportional reduction line corresponds to the dilution effect.At a curing age of 3 days, the normalized strength values are close to the dilution effect line, whereas, at the age of 28 days, they are much higher than the dilution effect line.This indicates the increasing reaction extent of the mineral admixtures from curing ages of 3 days to 28 days.

CO 2 Emissions
Table 5 shows the overall trend of CO 2 emissions: after adding SCM, CO 2 emissions were reduced.This is mainly because of the dilution effect of SCM reducing the amount of cement used.M21 has the lowest CO 2 emissions due to its higher SCM substitution.From the CO 2 emission calculation results (shown in Figure 5a), the CO 2 emissions of a paste were regressed based on a second-order polynomial as follows: CO Figure 5b presents the effect of the components on CO 2 emissions.As the individual replacement ratios of the binder components increase, the CO 2 emissions linearly decrease.It should be noticed in Figure 5b that the code +1 for limestone and calcined clay represents 10% and 20%, respectively.In other words, for limestone and calcined clay, when the x-axis represents the replacement ratio, the regression curves will be different.Figure 5c shows the normalized CO 2 emissions of the blended slurry.The normalized CO 2 emissions are the ratios of the CO 2 emissions of blended paste to those of control paste.Figure 5c shows that, for various mixtures, normalized CO 2 emissions decrease linearly with increasing total substitution rate.Therefore, adding SCMs is an effective way to make low-CO 2 concrete.sents 10% and 20%, respectively.In other words, for limestone and calcined clay, when the x-axis represents the replacement ratio, the regression curves will be different.Figure 5c shows the normalized CO2 emissions of the blended slurry.The normalized CO2 emissions are the ratios of the CO2 emissions of blended paste to those of control paste.Figure 5c shows that, for various mixtures, normalized CO2 emissions decrease linearly with increasing total substitution rate.Therefore, adding SCMs is an effective way to make low-CO2 concrete.

Multi-Objective Optimal Design
The optimized design had multiple goals: target strength and flo CO2 emissions.Desirability functions were used for the optimization.that the expected value is dependent on the response value [29], wi desirability function ranging from zero to one.For multi-objective optimization, based on the individual desirab composite desirability function can be determined using the following ( ) where D , n, and r are the composite desirability function, the total nu tion objections, and the relative importance of each response, respectiv Ecological aspects, such as low CO2 emissions, are an important sign.When the minimum value is set as the aim of optimization, the determined using the following equation [25]:

Multi-Objective Optimal Design
The optimized design had multiple goals: target strength and flow values and low CO 2 emissions.Desirability functions were used for the optimization.Equation (8) shows that the expected value is dependent on the response value [29], with the value of the desirability function ranging from zero to one.
Here, Y i and d i are the response value and the desirability function, respectively.Unity and zero correspond to fully achieved and not achieved optimization objectives, respectively.
For multi-objective optimization, based on the individual desirability functions, the composite desirability function can be determined using the following equation [25,29]: where D, n, and r are the composite desirability function, the total number of optimization objections, and the relative importance of each response, respectively.
Ecological aspects, such as low CO 2 emissions, are an important aim of optimal design.When the minimum value is set as the aim of optimization, the desirability can be determined using the following equation [25]: where L and U denote the lower limit and upper limit of the response value, respectively, and S denotes the weight factor.When the value of S is equal to 1, Equation ( 10) is a linear equation.We choose S > 1 when closeness to the target value is more important; we choose 0 < S < 1 when it is less important.Equation (10) shows that desirability equals 1 when the value of the response is below the lower bound.Desirability equals 0 when the value of the response is above the upper limit.
As shown in Table 5, through the combination of the four binder components, different responses, such as flow and 28-day strength, can be achieved.This can be achieved by optimizing the design.When a range is set as the aim of optimization, the desirability can be determined using the following equation: This equation shows that the desirability function is equal to 1 when the response is between the lower bound and upper bound, and it is 0 when the response is outside these bounds.
The optimization goals are listed in Table 7.The independent variable-the content of limestone, calcined clay, and slag-should be within the lower and upper limits.The variables of the response are as follows: (1) The 28-day strength should not be less than the target strength.In this study, the 28-day target strength values were set to 30, 35, 40, and 45 MPa.(2) The flow should be not less than the target flow (160 mm).( 3) The CO 2 emissions should be as low as possible for purposes of ecofriendliness.In this study, Design Expert software (version 12) was used for the optimization design [29].The relative importance r and weighting factor S of all responses were set to 1 [26].In summary, we consider four optimal design cases with different design strengths.Each design case considers basic requirements (strength and workability) and an ecolog-ical factor (CO 2 emissions).Multi-objective optimization was implemented based on a composite desirability function, which was calculated using the responses and individual desirability.

Results of Optimal Design
The results of the design case combinations are identified and summarized in Table 8.Mix-30, Mix-35, Mix-40, and Mix-45 have 28-day strengths of 30, 35, 40, and 45 MPa, respectively.Table 9 shows that the values of strength and flow for Mix-30, Mix-35, Mix-40, and Mix-45 meet the design requirements.9 show the independent variables (limestone, calcined clay, and slag) and responses (strength, flow, and CO 2 emissions) of Mix-30 to Mix-45, respectively.In Figure 6a, Figure 7a, Figure 8a, and Figure 9a, the first row represents limestone, clay, and slag contents, and the second row represents the results of multi-objective optimization (such as flow, strength, and CO 2 emissions).The optimized flow and strength meet the constraints of the target range, and the optimized CO 2 emissions are the minimum value.From the minimum value of 501.59 to the maximum value of 1052.04, the CO 2 optimization decreases with the falling line.This means that the closer the CO 2 emissions are to the minimum value of 501.59, the higher the satisfaction degree of the optimization; on the contrary, a lower optimization satisfaction degree occurs when the CO 2 emissions are closer to the maximum value of 1052.04. Figure 6a shows that the limestone, calcined clay, and slag contents of Mix-30 are 8.81%, 20%, and 30%, respectively.Solid red circles in the first row of Figure 6a show the component values of limestone, clay, and slag.The content of the mineral admixtures is close to the upper limit.Solid blue circles in the second row of Figure 6a show the response values for strength, flow, and CO 2 emissions.The strength and flow rate both meet the design requirements.The CO 2 emissions are close to a minimum.Figure 6b shows the desirability of the flow, 28-day strength, and CO 2 emissions, as well as the composite desirability.The desirability is 1 for both the 28-day strength and flow.The composite desirability value for Mix-30 is 0.977.
8.81%, 20%, and 30%, respectively.Solid red circles in the first row of Figure 6a show the component values of limestone, clay, and slag.The content of the mineral admixtures is close to the upper limit.Solid blue circles in the second row of Figure 6a show the response values for strength, flow, and CO2 emissions.The strength and flow rate both meet the design requirements.The CO2 emissions are close to a minimum.Figure 6b shows the desirability of the flow, 28-day strength, and CO2 emissions, as well as the composite desirability.The desirability is 1 for both the 28-day strength and flow.The composite desirability value for Mix-30 is 0.977.Figure 7 shows the desirability of Mix-35.The limestone, calcined clay, and slag contents of Mix-35 are 3.75%, 19.98%, and 30.00%, respectively (see the first row in Figure 7a).Since the replacement rate of mineral admixture in Mix-35 is smaller than that in Mix-30, Mix-35 has higher CO2 emissions than Mix-30.The desirability is 1 for both the strength and flow, but the CO2 emissions are less desirable than those for Mix-30.This is because as the strength of the concrete increases, its carbon dioxide emissions also increase, resulting in a decrease in the desirability function corresponding to CO2 emissions.In addition, Mix-35 has a lower composite desirability than Mix-30 (as shown in Figures 6b and 7b). Figure 7 shows the desirability of Mix-35.The limestone, calcined clay, and slag contents of Mix-35 are 3.75%, 19.98%, and 30.00%, respectively (see the first row in Figure 7a).Since the replacement rate of mineral admixture in Mix-35 is smaller than that in Mix-30, Mix-35 has higher CO2 emissions than Mix-30.The desirability is 1 for both the strength and flow, but the CO2 emissions are less desirable than those for Mix-30.This is because as the strength of the concrete increases, its carbon dioxide emissions also increase, resulting in a decrease in the desirability function corresponding to CO2 emissions.In addition, Mix-35 has a lower composite desirability than Mix-30 (as shown in Figures 6b and 7b).The desirability values of Mix-40 and Mix-45 are shown in Figures 8 and 9, respectively.Compared to Mix-30, Mix-35, and Mix-40, the total substitution rate of mineral admixture in Mix-45 is much lower.Therefore, Mix-45 has lower CO2 emissions than the other cases (as shown in Figure 8b).Consequently, the composite desirability of CO2 emissions for Mix-45 is also lower than those for the other cases.
(a) The desirability values of Mix-40 and Mix-45 are shown in Figures 8 and 9, respectively.Compared to Mix-30, Mix-35, and Mix-40, the total substitution rate of mineral admixture in Mix-45 is much lower.Therefore, Mix-45 has lower CO2 emissions than the other cases (as shown in Figure 8b).Consequently, the composite desirability of CO2 emissions for Mix-45 is also lower than those for the other cases.

Optimized Design of CO2 Emissions Based on 1 kg Binder
In the optimal design method shown in Section 4.2, the CO2 emission is based on each cubic meter of cement paste.In actual engineering, CO2 is often calculated based on each kilogram of binder.Based on the mix ratio and the CO2 emissions of binder components in Table 4, the CO2 emissions per kilogram of binder can be calculated (shown in Table 10).Moreover, the CO2 emission of 1 kg binder can be determined using a linear equation as follows: CO 2 of binder = 0.86 − 0.00852 *  1 -0.0059* 2 -0.0077* 3 (12) Following a similar method in Section 4.2, we conducted optimized designs for different strength levels, and the calculation results are shown in Figure 10.The combination of the mix ratio and performance of optimized designs are shown in Tables 11 and 12, respectively.According to these calculation results, we can find: First, as shown in Table 11, as the design strength increases, the cement content increases.As shown in Table 12, as the design strength increases, the CO2 emissions also increase.These trends are consistent with that in Section 4.2.Second, for Mix-30, Mix-35, Mix-40, and Mix-45, the overall satisfaction degree is 0.977, 0.877, 0.757, and 0.610, respectively.Comparing the overall satisfaction degree results in Section 4.2 and the satisfaction degree in this section, we can find that the satisfaction degree values are similar for CO2 emissions per cubic meter of paste and CO2 emissions per kilogram of binder.Overall, the design method proposed in this article has certain versatility and can be suitable for different optimization objectives.Figure 7 shows the desirability of Mix-35.The limestone, calcined clay, and slag contents of Mix-35 are 3.75%, 19.98%, and 30.00%, respectively (see the first row in Figure 7a).Since the replacement rate of mineral admixture in Mix-35 is smaller than that in Mix-30, Mix-35 has higher CO 2 emissions than Mix-30.The desirability is 1 for both the strength and flow, but the CO 2 emissions are less desirable than those for Mix-30.This is because as the strength of the concrete increases, its carbon dioxide emissions also increase, resulting in a decrease in the desirability function corresponding to CO 2 emissions.In addition, Mix-35 has a lower composite desirability than Mix-30 (as shown in Figures 6b and 7b).
The desirability values of Mix-40 and Mix-45 are shown in Figures 8 and 9, respectively.Compared to Mix-30, Mix-35, and Mix-40, the total substitution rate of mineral admixture in Mix-45 is much lower.Therefore, Mix-45 has lower CO 2 emissions than the other cases (as shown in Figure 8b).Consequently, the composite desirability of CO 2 emissions for Mix-45 is also lower than those for the other cases.

Optimized Design of CO 2 Emissions Based on 1 kg Binder
In the optimal design method shown in Section 4.2, the CO 2 emission is based on each cubic meter of cement paste.In actual engineering, CO 2 is often calculated based on each kilogram of binder.Based on the mix ratio and the CO 2 emissions of binder components in Table 4, the CO 2 emissions per kilogram of binder can be calculated (shown in Table 10).Moreover, the CO 2 emission of 1 kg binder can be determined using a linear equation as follows: CO 2 of binder = 0.86 − 0.00852 * x 1 −0.0059 * x 2 −0.0077 * x 3 (12) Following a similar method in Section 4.2, we conducted optimized designs for different strength levels, and the calculation results are shown in Figure 10.The combination of the mix ratio and performance of optimized designs are shown in Tables 11 and 12, respectively.According to these calculation results, we can find: First, as shown in Table 11, as the design strength increases, the cement content increases.As shown in Table 12, as the design strength increases, the CO 2 emissions also increase.These trends are consistent with that in Section 4.2.Second, for Mix-30, Mix-35, Mix-40, and Mix-45, the overall satisfaction degree is 0.977, 0.877, 0.757, and 0.610, respectively.Comparing the overall satisfaction degree results in Section 4.2 and the satisfaction degree in this section, we can find that the satisfaction degree values are similar for CO 2 emissions per cubic meter of paste and CO 2 emissions per kilogram of binder.Overall, the design method proposed in this article has certain versatility and can be suitable for different optimization objectives.

Discussion of the Optimal Design
In traditional design processes, concrete mix proportional design considers only strength and workability [17,30,31].However, with the development of concrete techniques, customers demand not only strength and workability but also better ecological performance, such as with regard to CO2 emissions [32,33].The model presented in this study considers more objectives than the conventional mix design approaches, which is important for performance-based design.
In traditional response surface optimization design, all optimization goals should be measured based on experimental results.However, there have been some new developments, and some researchers have adopted the response surface method to optimize concrete materials.In prior studies, the mechanical properties of the materials were measured experimentally, and the cost and CO2 emissions of concrete materials were calculated based on the mix ratio [33,34].Iman and Camoes [34] measured the strength and workability of ultra-high-performance concrete (UHPC), calculated the cost and CO2 emissions of UHPC using concrete components, and produced cost-efficient or CO2-efficient optimal designs using the response surface method.Similarly, Mohammad and Ozgur [33] measured the strength of UHPC, calculated the cost of UHPC, and determined optimal designs with a low cost.The rationality of these studies was verified by experiments.Inspired by these studies, we adopted a similar method.The strength and flow were measured experimentally, and the CO2 emissions were calculated using a formula.

Discussion of the Optimal Design
In traditional design processes, concrete mix proportional design considers only strength and workability [17,30,31].However, with the development of concrete techniques, customers demand not only strength and workability but also better ecological performance, such as with regard to CO 2 emissions [32,33].The model presented in this study considers more objectives than the conventional mix design approaches, which is important for performance-based design.
In traditional response surface optimization design, all optimization goals should be measured based on experimental results.However, there have been some new developments, and some researchers have adopted the response surface method to optimize concrete materials.In prior studies, the mechanical properties of the materials were measured experimentally, and the cost and CO 2 emissions of concrete materials were calculated based on the mix ratio [33,34].Iman and Camoes [34] measured the strength and workability of ultra-high-performance concrete (UHPC), calculated the cost and CO 2 emissions of UHPC using concrete components, and produced cost-efficient or CO 2 -efficient optimal designs using the response surface method.Similarly, Mohammad and Ozgur [33] measured the strength of UHPC, calculated the cost of UHPC, and determined optimal designs with a low cost.The rationality of these studies was verified by experiments.Inspired by these studies, we adopted a similar method.The strength and flow were measured experimentally, and the CO 2 emissions were calculated using a formula.
In addition, in the experimental design of this paper, the water used was distilled water.More recently, wastewater has also begun to be used in concrete materials.It is necessary to pay attention to the compositional differences between wastewater and distilled water and the resulting impact on the physical and chemical properties of concrete.The response surface design method proposed in this paper is suitable not only for distilled water concrete but also for wastewater concrete.

Conclusions
This study presents experimental works and multi-objective optimization design regarding a quaternary hybrid cement containing calcined clay, slag, and limestone.A total of 21 mixtures were prepared, of which 17 were BBD samples and 4 were additional samples.The flow, compressive strength (at 3 and 28 days), and CO 2 emissions of the various mixtures were determined.The following conclusions were obtained based on the test results: (1) Flow, strength, and CO 2 emissions decrease with increasing mineral admixture content.The flow decreases significantly with increasing calcined clay content.At 3 days of age, the normalized strength value is close to the dilution effect line, but at 28 days of age, the normalized strength value is much higher than the dilution effect line.Normalized CO 2 emissions decrease linearly as the substitution rate increases.(2) Four design cases (Mix-30, Mix-35, Mix-40, and Mix-45) with different 28-day design strengths (30,35,40, and 45 MPa) were considered.Each design case considered various aspects, namely basic requirements (strength and workability) and ecological aspects (CO 2 emissions).Multi-objective optimization was implemented based on a composite desirability function that was calculated using the individual responses and desirability for each objective.(3) From Mix-30 to Mix-45, as the content of supplementary cementing materials decreased, CO 2 emissions increased.Because the aim of CO 2 emissions optimization is to reach low CO 2 emissions, increasing CO 2 emissions lowered the individual desirability of the cases.In addition, the individual desirability values of strength and flow were 1 for each case.Consequently, the composite desirability decreased from 0.977 to 0.609 as the design strength increased from 30 to 45 MPa.The performance trend of the best combination was consistent with the experimental results.(4) The results presented in this paper can be used to guide a general method for designing low-carbon concrete.Adopting this method requires two steps.The first step is to produce an experimental design using the response surface method, conduct experimental research on the strength and fluidity, and calculate the CO 2 emissions.The second step is to optimize the design according to the required strength and flow level and choose the right combination of low-CO 2 cementitious materials.

Figure 1 .
Figure 1.XRD analysis of calcined clay and raw clay.

Figure 1 .
Figure 1.XRD analysis of calcined clay and raw clay.

Figure 3 .
Figure 3. Results of flow analysis.(a) Experimental versus predicted results.(b) Perturbations of flow.

Figure 3 .
Figure 3. Results of flow analysis.(a) Experimental versus predicted results.(b) Perturbations of flow.

Figure 4 .
Figure 4. Results of compressive strength analysis.(a) Experimental versus predicted results of 3day strength.(b) Experimental versus predicted results of 28-day strength.(c) Perturbations of variables of 3-day strength.(d) Perturbations of variables of 28-day strength.(e) Normalized 3-day strength.(f) Normalized 28-day strength.

Figure 4 .
Figure 4. Results of compressive strength analysis.(a) Experimental versus predicted results of 3-day strength.(b) Experimental versus predicted results of 28-day strength.(c) Perturbations of variables of 3-day strength.(d) Perturbations of variables of 28-day strength.(e) Normalized 3-day strength.(f) Normalized 28-day strength.

d
are the response value and the desirability fun Unity and zero correspond to fully achieved and not achieved optim respectively.

Figure 5 .
Figure 5. Results of CO 2 emissions analysis.(a) Calculated versus regressed results of CO 2 emissions.(b) Perturbations of variables of CO 2 emissions.(c) CO 2 emissions versus replacement ratio.
a Loss on ignition.

Table 2 .
Replacement ranges of mineral admixtures.

Table 5 .
Test results and CO 2 emissions.

Table 6 .
Coefficients of regression of the equations.

Table 7 .
Objectives of optimization.

Table 8 .
Results of optimal mixtures.

Table 9 .
Performance of optimal mixtures.

Table 11 .
Optimization design results (CO 2 emissions based on 1 kg binder).

Table 12 .
Performance of optimal mixtures (CO 2 emissions based on 1 kg binder).