Abstract

Metakaolin (MK) is an aluminosilicate pozzolan material and can contribute to property development of concrete due to the pozzolanic reaction. Limestone (LS) powder presents the dilution effect, nucleation effect, and chemical effect on hydration of cement. When metakaolin and limestone are used together, due to the additional chemical reaction between the aluminum phase in MK and limestone, the synergetic benefit can be achieved. This study presents a hydration model for cement-limestone-metakaolin ternary blends. Individual reactions of cement, metakaolin, and limestone are simulated separately, and the interactions among cement hydration, limestone reaction, and metakaolin reaction are considered through the contents of calcium hydroxide and capillary water. The hydration model considers the pozzolanic reaction of metakaolin, chemical and physics effects of limestone, and synergetic effect between metakaolin and limestone. Furthermore, the gel-space ratio of hydrating concrete is calculated using reaction degrees of binders and concrete mixtures. The strength development of ternary blends is evaluated using the gel-space ratio. Based on parameter analysis, the synergetic effect on strength development is shown and the optimal combinations of cement-limestone-metakaolin ternary blends are determined.

1. Introduction

Metakaolin is increasingly used in modern concrete industry. The addition of metakaolin can make many advantages on the performance of concrete. Metakaolin can enhance workability and finishing ability, increase compressive and flexural strength, and reduce chloride permeability. However, metakaolin-blended binary concrete has some drawbacks. As the replacement level of cement by metakaolin increases, an increased amount of superplasticizer is necessary to reach the required consistency. The material cost of the metakaolin is higher than that of cement. To avoid these drawbacks, the substitution level of cement by metakaolin is generally lower than 25%. Summarily, the benefits of metakaolin, such as mechanical performance and extended service life, may overcome its negative effect. Metakaolin is a very promising supplementary cementitious material for the concrete industry [1, 2].

Limestone powder can improve the workability, reduce the bleeding, and reduce the amount of CO2 emission of the concrete industry. The price of limestone is much lower than that of metakaolin. The addition of limestone lowers the late-age strength of concrete. When limestone and metakaolin are used together, the pozzolanic reaction of metakaolin can contribute to concrete late-age strength. In addition, metakaolin has a high aluminum content which can react with limestone, form carboaluminate phases, and increase solid volume and strength of concrete. This is the synergetic effect of ternary concrete. While metakaolin or limestone is added individually, the synergetic effect cannot be achieved. In summary, when using metakaolin- and limestone ternary-blended concrete, the economic benefit and strength benefit can be achieved [1, 2].

Many experimental studies have been done about workability, mechanical, and durability aspects of cement-limestone-metakaolin ternary blends. Vance et al. [3] reported that, for cement-limestone-metakaolin ternary blends, the yield stress reduces as the limestone content increases. This is because of the particle packing, water demand, and interparticle spacing and contacts. Vance et al. [4] found that the synergetic effect of limestone and MK incorporation can improve early-age properties and maintain late-age properties of concrete. Alvarez et al. [5] presented that combined mixtures of limestone and MK enhance compressive strength compared with 100% Portland cement concrete. Ramezanianpour and Hooton [6] presented that carboaluminate hydrates are formed for cement-limestone-MK ternary blends, and there is an optimum level of limestone in terms of the maximum strength and minimum porosity. Perlot et al. [7] reported that ternary-blended mix presents a real benefit for carbonation durability performances due to the refinement of the pore structure. Sotiriadis et al. [8] reported that the use of MK in the limestone cement concrete retards and inhibits deterioration due to the sulfate attack and improves its durability.

Compared with abundant experimental studies, the theoretical models about cement-limestone-metakaolin ternary blends are very limited. Antoni et al. [9] made the thermodynamic model for cement-limestone-metakaolin ternary blends and present phase assemblage for different combinations of binders. Shi et al. [10, 11] presented the thermodynamic model for carbonation and chloride ingression of cement-limestone-metakaolin ternary blends. Changes in phase assemblages and total porosities due to carbonation and chloride exposure are evaluated [10, 11]. The thermodynamic models [911] mainly focus on the chemical aspect of ternary blends, such as the phase assemblages of hydrated concrete and reaction products between concrete and ingression ions. However, limited works was done about mechanical aspects, such as the evaluation of strength development and optimal combinations of binders of ternary blends [911].

Optimal combinations of binders of ternary blends are an interesting topic for concrete manufacture and construction company. This study presents a blended hydration model for ternary blends considering the synergetic effect between limestone and MK. The strength development of ternary blends is evaluated using reaction degrees of binders and gel-space ratios. Based on parameter analysis, the optimal combinations of cement-limestone-metakaolin ternary blends are determined.

2. Hydration Model and Strength Model for Ternary Blends

2.1. Hydration Model

For MK- and limestone-blended concrete, the hydration of cement, the reaction of MK, and reaction of limestone coexist. In this study, we simulate the hydration of cement, MK, and limestone, respectively. Meanwhile, the interactions among cement hydration, metakaolin reaction, and limestone reaction are considered through the contents of capillary water and calcium hydroxide in hydrating blends.

2.1.1. Cement Hydration Model

The hydration of cement can be described using a kinetic model shown in our former studies [12]. The degree of hydration can be calculated as , where t is the time and is the rate of hydration, which can be determined as follows:where and are the rate-determining coefficients in the initial dormant period, is the reaction parameter in the diffusion-controlled period, is the reaction parameter in the phase boundary-controlled period, considers the reduction of the hydration rate due to development of microstructure, and considers the reduction of hydration rate due to the consumption of capillary water (, where is the amount of capillary water and is the amount of water in concrete mixtures).

The kinetic processes involved in cement hydration, such as the initial dormant process, phase boundary reaction process, and diffusion process, are considered in the cement hydration model. The input variables of the cement hydration model are cement compound compositions, Blaine surface of cement, concrete mixing proportions, and curing conditions. The values of hydration parameters , , , and can be determined using cement compound compositions. Furthermore, the reaction degree of cement can be calculated automatically using hydration parameters , , , and . The effect of curing temperature on cement hydration is considered using Arrhenius law [12]. For high-strength concrete, the water-to-cement ratio is low, and the hydration rate is significantly lowered due to the reduction of capillary water [13, 14]. This effect is considered using capillary water concentration . Summarily, the proposed cement hydration model is valid for concrete with different strength levels, different types of Portland cement, and different curing conditions [12].

2.1.2. MK Reaction Model

The reaction process of MK also consists of the initial dormant process, phase boundary reaction process, and diffusion process, which is similar with the processes involved in hydration of cement [15]. On the contrary, MK is a pozzolanic material. The reaction rate of the pozzolanic reaction is dependent on the amount of calcium hydroxide in blends [16, 17]. Considering the kinetic reaction processes and the essence of MK as the pozzolanic material, Wang [15] proposed that the reaction equation of MK can be written as follows:where is the reaction degree of MK, is the reaction rate of MK, and are reaction parameters of MK in the dormant period, is the reaction parameter of MK in the diffusion-controlled period, is the reaction parameter of MK in the phase boundary-controlled period, is the content of calcium hydroxide in blends, and is the content of MK in concrete mixtures. The verifications of the MK reaction model are available in our former study [15]. An integrated hydration-strength-durability model for MK-blended concrete is proposed for evaluating the reaction degrees of binders, strength development, and chloride penetrability [15].

2.1.3. Limestone Reaction Model

The addition of limestone powder presents the dilution effect, nucleation effect, and chemical effect on hydration of cement. In this study, dilution effect is considered through the amount of capillary water, nucleation effect is considered through a nucleation effect indicator, and chemical effect is considered through a logarithm function with multiple modification factors [18, 19].

The hydration products of cement can form on the surface of limestone powder. This is called as the nucleation effect. The nucleation effect indicator of limestone powder can be written as follows [2]:where is the limestone nucleation effect indicator; and C0 are the mass of limestone and cement in mixing proportions, respectively; and and are the Blaine surface area of limestone powder and cement, respectively.

In our former study [2], based on the experimental results of hydration degree of cement in cement-limestone binary blends, Wang and Luan [2] proposed that the nucleation effect of limestone powder can be described as follows:where is the updated phase boundary reaction coefficient in cement-limestone blends, 1.2 is enhancing coefficients of [2], is the updated diffusion coefficient in cement-limestone blends, and is enhancing coefficients of [2].

Until now, the experimental results about the reaction degree of limestone are very limit. Tentatively, Wang and Luan [2] proposed an empirical model with multimodification factors for analyzing the reaction degree of limestone. The empirical model considers the effects of various factors on the chemical reaction of limestone, such as limestone replacement ratios, mineral admixtures additions, limestone fineness, cement fineness, water-to-binder ratio, and curing temperature. The empirical model for the limestone reaction is shown as follows:where is the reaction degree of limestone in a reference mixture. This reference mixture is Portland cement, and limestone binary blends with a water-to-binder ratio 0.5 and 20% limestone addition cured at 20°C. considers the effect of limestone replacement ratios on the reaction degree of limestone, considers the effect of limestone fineness, considers the effect of cement fineness, considers the effect of MK addition, considers the effect of the water-to-binder ratio, and considers the effect of the curing temperature. Table 1 shows the summary of influencing factors of the limestone reaction. As limestone replacement ratio increases, reaction degree of LS decreases. While fineness of limestone, fineness of cement, MK addition, and water to binder ratio increase, the reaction degree of LS increases. Especially, for modification factor , where is the aluminum content in MK, is the aluminum content in cement, in the numerator is reacted aluminum content from the MK reaction, and in denominator is the reacted aluminum content from the cement reaction. Because the aluminum content in MK is much higher than that in cement, the addition of MK can significantly improve the reactivity of limestone. The factor m4 considers the synergetic effect of limestone and MK. A higher aluminum content and a higher reactivity of MK are effective to enhance the reactivity of limestone.

Summarily, this study considers the dilution effect, nucleation effect, and chemical effect of limestone additions. The enhancement of limestone reactivity due to the addition of metakaolin is considered through a modification factor. The influences of other factors, such as limestone replacement ratio, fineness of binders, and water-to-binder ratios, are also considered in the limestone reaction model.

2.1.4. Interaction Model among Cement, Metakaolin, and Limestone

In this study, the interactions among cement hydration, metakaolin reaction, and limestone reaction are considered through the contents of capillary water and calcium hydroxide. Maekawa et al. [13] proposed that, as 1 g cement hydrates, 0.4 g capillary water will be consumed. Dunster et al. [20] proposed that, as 1 g metakaolin reacts, 0.55 g capillary water will be consumed. Bentz [21] proposed that, as 1 g limestone reacts, 1.62 g capillary water will be consumed. For hydrating cement-metakaolin-limestone ternary blends, the content of capillary water can be determined as follows:where , , and are the contents of consumed water from cement hydration, metakaolin reaction, and limestone reaction, respectively [13, 20, 21]. The consumption of capillary water from the 1 g limestone reaction is much higher than those of cement and metakaolin. This is because the reaction products of limestone powder are monocarboaluminate and ettringite which contains abundant water.

For hydrating cement-metakaolin-limestone ternary blends, the content of calcium hydroxide can be determined as follows:where means the mass of CH produced from the hydration of 1 unit mass of cement and means the mass of CH consumed from the reaction of 1 unit mass of metakaolin [15]. is the mass of CH produced from cement hydration. is the mass of CH consumed from the metakaolin reaction.

Summarily, the effect of cement hydration, metakaolin reaction, and limestone reaction on the contents of capillary water and calcium hydroxide is considered. The capillary water content can be used for the cement hydration model (equation (1)), and the calcium hydroxide content can be used for the metakaolin reaction model (equation (2)). In addition, because the blended hydration model has considered the interactions among cement, metakaolin, and limestone reactions, the coefficients of the hydration model do not vary with different mixtures. When the mixtures change from one to another, the coefficients of the hydration model are constant.

2.2. Strength Development Model

The gel-space ratio denotes the ratio of the volume of binder hydration products to the sum of volume of hydrated binders and capillary pore. For cement-metakaolin-limestone blends, 1 ml hydrated cement, 1 ml reacted metakaolin, and 1 ml reacted limestone occupy 2.06 ml of space [16, 22], 2.52 ml space [16, 22], and 4.1 ml space, respectively. Reacted products of 1 ml limestone can occupy much higher space that those of cement (4.1 vs. 2.06). This is due to the development of ettringite and monocarboaluminate from the limestone reaction. Considering the reactions of cement, metakaolin, and limestone, the gel-space ratio of cement-metakaolin-limestone ternary-blended cement can be determined as follows:where , , and are densities of cement, metakaolin, and limestone powder, respectively.

According to Powers’ strength theory, the compressive strength of hydrating concrete can be evaluated using the gel-space ratio as follows:where is the concrete compressive strength, is the intrinsic strength of concrete, and n is the strength exponent.

For cement-metakaolin-limestone blends, cement, metakaolin, and limestone will affect the intrinsic strength of concrete and strength exponent. We assume that the intrinsic strength of concrete and strength exponent n is proportional to the weight fractions of cement, metakaolin, and limestone in the mixing proportion as follows:where coefficients , , and in equation (12) represent the contributions of cement, metakaolin, and limestone to the intrinsic strength of concrete, respectively, and the units of a1, a2, and a3 are MPa; the coefficients , b2, and b3 in equation (13) represent the contributions of cement, metakaolin, and limestone to the strength exponent, respectively. For neat Portland cement concrete without limestone or metakaolin, the strength of concrete only pertains to a1 and b1. For metakaolin-blended binary concrete without limestone, the strength of concrete pertains to coefficients a1, a2, b1, and b2. For ternary-blended concrete, the strength of concrete pertains to coefficients a1, a2, a3, b1, b2, and b3. These coefficients a1, a2, a3, b1, b2, and b3 do not change for various mixing proportions of concrete.

The flowchart of calculation is proven in Figure 1. Each time step, the response levels of cement, metakaolin, and limestone powder are calculated by utilizing ternary-blended hydration model. The quantity of CH and capillary water are based on using reaction levels of binders and concrete mixtures. In addition, the gel-space ratio of hydrating concrete is decided, thinking about the contributions from reactions of cement, metakaolin, and limestone reactions. By utilizing Powers’ strength theory, the compressive strength of hardening concrete is calculated.

3. Verifications of Proposed Models

3.1. Verification of Hydration Model

Experimental results from Antoni et al. [9] are used to verify the proposed blended hydration model and strength development model. Antoni et al. [9] measured the reaction degrees of binder and compressive strength of cement-MK-LS ternary-blended concrete. The chemical compositions of cement, metakaolin, and limestone are shown in Table 2. The mixing proportions are shown in Table 3. Paste specimens with a water-to-binder ratio of 0.4 were used for measuring reaction degree of binders. Mortar specimens with a water-to-binder ratio of 0.5 were used for measuring compressive strength. For cement-limestone binary blends, the replacement ratio of limestone was 15%, while for cement-metakaolin binary blends, the replacement ratio of metakaolin was 30%. For ternary-blended specimens, the sum of limestone and metakaolin ranged from 15% to 60%, and the mass ratio of metakaolin to limestone was fixed as 2. The reaction degrees and strength were measured at the ages of 1, 7, 28, and 90 days.

The input parameters of the ternary-blended cement hydration model are concrete mixtures, curing temperature, and compound compositions and Blaine surface areas of binders. By using the blended cement hydration model, the reaction degree of MK and LS is calculated and shown in Figure 2.

As shown in Figure 2(a), the sequences of the reaction degree of MK from higher to lower are B15 > B30 > MK30 > B45 > B60. This can be explained using the MK reaction model (equation (2)). As shown in equation (2), the reaction degree of MK mainly depends on the mass ratio of cement to MK. As the ratio of cement to MK increases, the activation effect from cement hydration is enhanced, and reaction degree of MK increases. The mass ratios of cement to MK were 8.5, 3.5, 2.33, 1.83, and 1 in the mixtures of B15, B30, MK30, B45, and B60, respectively. The orders of reaction degree of MK are consistent with the mass ratios of cement to MK.

As shown in Figure 2(b), the sequences of reaction degree of limestone from higher to lower are B15 > B30 > B45 > B60. The proposed ternary-blended hydration model can reflect this trend of reaction degree of LS. In this study, the mass ratio of MK to LS in ternary blends is constant, and the difference of the reaction degree of LS is mainly due to the variations of the cement-to-limestone ratio. The mass ratios of cement to LS were 17, 7, 3.66, and 2 in the mixtures of B15, B30, B45, and B60, respectively. As the mass ratio of cement to LS decreases, the reaction degree of LS also decreases (parameter m1 of equation (7)). The trend of the reaction degree of LS is consistent with the mass ratio of cement to LS. In addition, the reaction degree of LS at 1 day is almost zero. This also agrees with our analysis. As shown in equation (6), we assumed that the reaction of LS starts after 21 hours. Furthermore, the reactivity of LS is very low. At the age of 90 days, the reaction degree of LS for B15 is 12% which is much lower than cement.

Figure 3 shows the parameter analysis of the hydration model. Figure 3(a) shows the reaction degree of LS in cement-LS binary blends. As limestone replacement ratio increases, the reaction degree of limestone decreases. Similar to the contents shown in Figure 3(a), Aqel and Panesar [23] also found the reactivity of limestone will be lower with the increasing of limestone content.

Figure 3(b) shows the reaction degree of MK in cement-MK binary blends. As MK replacement ratio increases, the activation effect from cement hydration becomes weaker, and the reaction degree of MK decreases. Similar to the contents shown in Figure 3(b), Poon et al. [24] also found similar results that the reaction degree of MK will be lower as the content of MK increases.

Figure 3(c) shows the effect of MK additions on the reaction degree of LS. The addition of MK presents the twofold effect on reaction of LS. First, when MK is added to replace partial cement in the mixtures, the mass ratio of cement to LS decreases which will lower the reaction degree of LS (this is considered through parameter m1 of equation (7)). However, the aluminum content in MK (46%) is about ten times of the aluminum content in cement (4.6%). The addition of MK will enhance the reaction of LS (this is considered through parameter m4 of equation (7)). Because the enhancing effect is much more significant than the lowering effect, the addition of MK can increase the reaction degree of limestone (shown in Figure 3(c)). Similar to the contents shown in Figure 3(c), many researchers [4, 9] also experimentally found that reactivity of limestone can be improved due to MK addition.

Figure 3(d) shows the effect of metakaolin and limestone contents on the reaction degree of cement. When metakaolin and limestone are used to replace partial cement, the reaction degree of cement is improved due to the dilution effect and nucleation effect (the dilution effect is considered through parameter in equation (1), and the nucleation effect is considered through equations (4) and (5)). Similar to the contents shown in Figure 3(d), Lam et al. [25] also found that the addition of mineral admixtures can improve the reaction degree of cement.

3.2. Verification of Strength Development Model

By using the cement-MK-LS ternary-blended hydration model, the gel-space ratio of hydrating concrete can be calculated (equation (10)). Furthermore, based on the strength of concrete at different ages, the values of strength coefficients of a1, a2, and a3 and b1, b2, and b3 can be calibrated (a1 = 140 MPa, a2 = 258 MPa, a3 = 120 MPa, b1 = 3.85, b2 = 1.13, and b3 = 1.34). These coefficients do not vary with concrete mixtures. The values of a1 and b1 relate to cement hydration, the values of a2 and b2 relate to the metakaolin reaction, and the values of a3 and b3 relate to the limestone reaction. For cement-metakaolin binary blends, the development of strength relates to a1, a2, b1, and b2. For cement-limestone binary blends, the development of strength relates to a1, a3, b1, and b3. For cement-metakaolin-limestone ternary blends, the development of strength relates to a1, a2, a3, b1, b2, and b3. The analyzed results of compressive strength are shown in Figure 4. The analysis results generally agree with experimental results. At the ages of 28 days, the B15 concrete (cement 85% + metakaolin 10% + limestone 5%) has a highest strength than other mixtures. This may be because of the synergetic effect of metakaolin and limestone.

Because the strength coefficients of strength evaluation equation are constants for different concrete mixtures, we can make parameter analysis for different concrete mixtures. Figure 5(a) shows the strength development of cement-limestone binary blends. At the early age, due to the nucleation effect, the strength of limestone blends concrete shows higher strength than control concrete. While at late ages, due to the dilution effect, the strength of limestone blends concrete is lower than control concrete. As the contents of limestone increases from 10% to 20%, the late-age strength decreases. The trend shown in Figure 5(a) agrees with Bonavetti et al.’s [19] studies about strength development of limestone-blended concrete.

Figure 5(b) shows the strength development of cement-metakaolin binary blends. Metakaolin-blended concrete has a higher strength than control concrete. As the contents of metakaolin increases from 5% to 10%, the strength also increases. The trend shown in Figure 5(b) agrees with Poon et al.’s [24] studies about strength development of metakaolin-blended concrete.

Damidot et al. [26] studies the strength development of 70% cement + 30% clay-limestone ternary blends. The sum of clay and limestone was fixed as 30%, and the weight fraction of clay/(clay + limestone) ranges from 0 to 100% [26]. Damidot et al. [26] found that, at the age of 28 days, the mix with 70% metakaolin has the highest strength than other mixes. This is because of the synergetic effect of limestone and metakaolin [26]. Based on the proposed strength development in this study, we make parameter analysis of strength development for 70% cement + 30% clay-limestone ternary blends. In our analysis, the sum of clay and limestone is also fixed as 30%, the weight fractions of clay/(clay + limestone) are given as 0, 25%, 50%, 75%, and 100%, and the ages of parameter analysis are 1.5 days, 3 days, 28 days, and 90 days, respectively. The results of parameter analysis are shown in Figures 6(a)6(d). As shown in Figure 6(a), at the age of 1.5 days, the strength of blended concrete is higher than base Portland cement. This is because of the nucleation effect of limestone. While as shown in Figures 6(b)6(d), at the age of 3 days, 28 days, and 90 days, when the mk/(mk + limestone) equals to zero (the content of metakaolin is zero, and binder consists of 30% limestone and 70% cement), the strength of limestone-blended concrete is lower than base Portland cement. This is because of the dilution effect of limestone. While for other mk/(mk + limestone) ratios of 25%, 50%, 75%, and 100%, because the reaction of metakaolin can contribute to the strength, the strength of blended concrete is higher than base Portland cement.

At the age of 1.5 days, 3 days, 28 days, and 90 days, the optimum weight fractions of clay/(clay + limestone) are 25%, 50%, 75%, and 75%, respectively (shown in Figures 6(e)). Our analysis result about optimum weight fraction of clay/(clay + limestone) is similar to that of Damidot et al. [26] studies. In addition, our analysis shows that, at the ages of 1.5 days, 3 days, 28 days, and 90 days, the optimum weight fraction of limestone/(clay + limestone) is 75%, 50%, 25%, and 25%, respectively. It means that, at the early age, limestone is effective to improve the strength of concrete (this is because of the nucleation effect of limestone), and at the late age, metakaolin is effective to improve the strength of concrete (this is because of the pozzolanic reaction of metakaolin).

In Figure 6, the sum of metakaolin and limestone is fixed as 30%. To find the optimum combinations of cement, metakaolin, and limestone, we make much wider parameter analysis. In this wider parameter analysis, the sum of metakaolin and limestone is not a fixed value. The contents of metakaolin vary from 0 to 30%, and the contents of limestone vary from 0 to 20%. The analysis results of isoline of strength are shown in Figures 7(a)7(d). At the early age 1.5 days, the concrete with a higher limestone content and a lower metakaolin content has highest strength (shown in Figure 7(a)), while at the late age of 90 days, the concrete with a higher metakaolin content and a lower limestone content has highest strength (shown in Figure 7(d)). In other words, to achieve the highest strength of cement-metakaolin-limestone ternary blends, the optimum combination of metakaolin and limestone is dependent on ages. From the early age to late age, the optimum combinations change shift from high limestone-low metakaolin zone to low limestone-high metakaolin zone (shown in Figure 7(e)).

4. Conclusions

This study presents an integrated hydration-strength model for cement-limestone-metakaolin ternary blends.

First, a cement hydration model, a metakaolin reaction model, and a limestone reaction model are proposed in the ternary-blended hydration model. Pozzolanic reaction of metakaolin, chemical and physics effect of limestone, and synergetic effect between metakaolin and limestone are detailed considered in the ternary-blended hydration model. Moreover, the interactions among cement hydration, limestone reaction, and metakaolin reaction are considered through the contents of calcium hydroxide and capillary water. The coefficients of the hydration model do not change for various concrete mixtures.

Second, based on the hydration model, the gel-space ratio of hydrating blends is calculated considering the contributions from the reactions of cement, metakaolin, and limestone. Furthermore, the strength development of ternary blends is evaluated using the gel-space ratio. The coefficients of the strength model do not change for various concrete mixtures. Based on parameter analysis, the synergetic effect on strength development is shown, and the optimal combinations of cement-limestone-metakaolin ternary blends are determined. From the early age to late age, the optimum combinations of ternary blends shift from high limestone-low metakaolin zone to low limestone-high metakaolin zone.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (no. 2015R1A5A1037548) and a NRF grant (NRF-2017R1C1B1010076).