Abstract

Mg-35Zn-3Ca (wt.%) alloy containing nontoxic and biocompatible Zn and Ca as alloying elements was prepared and subjected to heat treatment and artificial aging for different duration of time to reduce its rate of degradation. Solution heat treatment was performed at 310°C while artificial aging was performed at 170°C for 0, 2.5, 5.0, 7.5, and 10.0 h and they were designated as AT0, AT1, AT2, AT3, and AT4, respectively. The finest and most homogenous reticulum was observed on the surface of the AT2 group. The result of immersion test in Hank’s balanced salt solution (HBSS) showed that the corrosion rate of the AT2 group was 2.32 mg/(cm2 day), which was significantly lower as compared to other groups . The hemolysis value was ≤5% in all groups, indicating no toxicity during short-term blood reaction.

1. Introduction

Implants are widely used for the replacement and treatment of abnormally functioning bone tissues due to diseases or injuries. Among the various types of implants that are currently being used, ceramic- or polymer-based implants have unsatisfactory mechanical properties while metal-based ones are considered suitable for load bearing implant applications due to their superior mechanical strength and durability. Stainless steel, titanium and its alloys, and Co-Cr alloys are some of the commonly used materials for fabricating metallic implants. However, these implants pose a high risk of release of toxic heavy metals or generating fine metallic debris particles induced by corrosion or wear or by tribocorrosion, which could cause inflammation in the surrounding tissues [15]. In addition, the difference in Young’s modulus between these implants and bone tissues could cause stress-shielding effect, resulting in stress concentration on the contact plane between the natural bone and the implant. This leads to a reduction of the stimuli required for the new bone growth, thereby lowering the success rate of the implant [6].

Mg and its alloys possess a low density, high mechanical strength, and fracture toughness. In addition, they are biocompatible and biodegradable [7]. Young’s modulus of Mg and its alloys are closer to the natural bone than stainless steels, titanium and its alloys, and Co-Cr alloys. Mg is the fourth abundant ion in the human body and it is an essential component for human metabolism. The daily Mg intake of a normal adult is 300–400 mg and the excess Mg ions are excreted via urine [812]. Mg, which is a natural component of bone tissues, helps bone growth and increases bone strength [13]. These attributes form the basis that Mg and its alloys can be explored towards the development of degradable biomaterials. For load-bearing application, biodegradable implants must be retained in the body for at least 12–18 weeks, a time period which is usually required for natural growth of bone tissues after treatment or surgery, while maintaining sufficient mechanical strength [14, 15]. Unfortunately, Mg and its alloys undergo rapid corrosion in chloride containing solutions such as body fluid, resulting in the generation of a huge volume of hydrogen gas with alkalinization of the medium [16]. The hydrogen gas thus generated could cavern under the subcutaneous tissue and under a worst scenario, the hydrogen gas bubbles could block the blood flow, causing death of the patient [8]. Thus, studies on purification, alloying, and surface coating have been explored to resolve the rapid biodegradation of Mg and its alloys.

Among the various alloys evaluated for biodegradable implant applications, Al-Zn based Mg alloys such as AZ91 and AZ31, as well as rare earth based Mg alloys such as LAE442 and WE43, have received considerable attention. However, when Al-Zn based Mg alloys are used as degradable implants, the accumulation of aluminum in the body could cause neurogenic disorders such as Alzheimer’s disease, particularly dementia, and senile dementia [17]. In case of rare earth based Mg alloys, cerium, praseodymium, and yttrium show hepatotoxicity and hence they are considered inappropriate for use as degradable implants [18]. Calcium is an important component of the natural bone. As a biomaterial, it has advantages for being essential in cellular signal transduction and for having density similar to that of the bone [19]. When Ca is added to Mg alloy, Mg2Ca with high melting point is formed and a protective oxide layer is also formed to improve corrosion resistance, thereby increasing hardness and creep resistance [20]. Zinc is one of the typical inorganic matters that constitutes somatic cells and regulates physiological functions. Its deficiency in the body could cause various anomalies such as low birth weight during delivery and growth and developmental problems. If Zn and Ca are added to Mg, a stable intermetallic compound Ca2Mg6Zn3 is formed that could help in age-hardening [21] and in increasing the tensile strength and ductility of the alloy [22].

In the present study, to reduce the biodegradation rate of Mg and to improve its biocompatibility in the body, a ternary Mg alloy using Zn and Ca as nontoxic and biocompatible alloying elements was prepared. The as-prepared ternary Mg-Zn-Ca alloy was subjected to heat treatment followed by artificial aging for different duration of time to alter the microstructure so as to improve the corrosion resistance, biocompatibility, and hemocompatibility.

2. Materials and Methods

2.1. Materials

Mg-35Zn-3Ca (wt.%) alloys were prepared using magnesium (99.93%, Sincere East Co., China), zinc (99.90%, Sincere East Co., China), and calcium (99.50%, Junsei Chemical Co., Ltd., Japan). The experimental alloys were melted in vacuum arc remelting furnace (Ace Vacuum Co., Korea) equipped with water-cooling Cu hearth. Vacuum arc remelting was performed in presence of argon gas (99.99%) at 5.0 × 10−5 Torr, under a current and voltage of 300 A and 12 V, respectively. Melting was performed five or more times to prevent segregation and remelting 10 more times while flipping over to obtain homogeneity in Mg-Ca master alloy.

Mg-35Zn-3Ca alloy was manufactured using a vacuum induction melting furnace (Ace Vacuum Co., Korea). Zinc was loaded on the bottom of alumina crucible, and the Mg-Ca master alloy and Mg element were cross-loaded on top of the loaded zinc. The alumina crucible was preheated for 200 s, followed by melting and agitation for 10 min by increasing the power in an interval of 200 s. The solution was poured into the alumina cast to produce Mg-35Zn-3Ca alloy.

For the microstructure observation and assessment of the corrosion resistance of heat treated cast alloy, a 16.4 mm diameter Mg stick was sliced into 2 mm thickness before use.

2.2. Methods
2.2.1. Heat Treatment and Artificial Aging

In order to determine the suitable temperature for solution heat treatment, the Mg-35Zn-3Ca alloy specimens were analyzed using Differential Scanning Calorimetry (DSC) (TA Q20, TA Instruments Ltd., USA). Both the Mg-Zn-Ca alloy specimen and reference materials underwent heating/cooling at a constant rate under the same condition. The temperature difference that occurs between the Mg-Zn-Ca alloy specimen and reference material either by exothermic or endothermic reactions was adjusted by supplying the required amount of energy and the transformation temperature was determined by measuring the amount of the supplied energy.

For solution heat treatment, the Mg-Zn-Ca alloy specimens were placed inside an electric arc furnace (Ajeon Industrial Co. Ltd., Korea) and the temperature of the furnace was raised to 310°C at a heating rate of 10°C/min and the samples were soaked for 48 h followed by quenching in distilled water at room temperature. Some of the specimens were subjected to artificial aging following solution heat treatment. For artificial aging, the specimens were placed in an electric arc furnace and the temperature of the furnace was raised to 170°C at a heating rate of 10°C/min. Then the specimens were soaked for various time intervals from 2.5 to 10.0 h. Specimens which were subjected to only heat treatment at 310°C were designated as AT0 whereas those subjected to an additional artificial aging for 2.5, 5.0, 7.5, and 10.0 h were designated as AT1, AT2, AT3, and AT4, respectively. The details of the experimental conditions employed for heat treatment and artificial aging are compiled in Table 1.

2.2.2. Surface Examination and Analysis

The heat treated and artificially aged Mg alloy specimens were sequentially polished by SiC coated abrasive papers (grit size #800 to #1600), followed by polishing using 0.3 μm Al2O3 powder and ultrasonic cleaning in distilled water for 5 min. The polished specimens were carefully etched in 2% Nital etchant (ethanol : HNO3 = 98 mL : 2 mL) for 5 s with a periodic assessment of surface microstructure for any sign of acidic corrosion. The microstructure of the Mg-Zn-Ca alloy specimens was examined by scanning electron microscopy (SEM) (JSM-5900, JEOL, Japan) (low magnification images) and by field emission scanning electron microscopy (FE-SEM) (SUPRA 40VP, Carl Zeiss Co., Ltd., Germany) (high magnification images). In addition, the concentration of the elements present in the alloy and their distribution were determined using energy dispersive X-ray analysis (EDX) (Oxford, England).

2.2.3. Potentiodynamic Polarization Test

Potentiodynamic polarization test was used to assess the electrochemical corrosion behavior of Mg-Zn-Ca alloy specimens. An electrochemical cell consisting of Ag/AgCl, KCl (sat’d) as the reference electrode, a platinum plate as the counter electrode, and Mg-Zn-Ca alloy specimens subjected to heat treatment and artificial aging for different duration of time as the working electrode, connected to a potentiostat/galvanostat (273A, EG&G, PARC, USA), was used for the study. The experiments were performed with exposing the area 1 cm2 in 150 mL of Hanks solution, whose chemical composition is presented in Table 2. The potentiodynamic polarization test was performed in the potential range −0.5 V to +0.5 V at a scan rate of 3 mV/s at room temperature. The corrosion potential () and current density () were calculated from the measured polarization curves using Tafel extrapolation method, followed by assessment of extent of corrosion the specimens.

2.2.4. Immersion Test

The initial pH of the HBSS was adjusted to 7.2 with the incremental additions of either NaOH (0.1 mol/L) or HCl (0.1 mol/L) using a pH meter (HD-PH/P, Davis Co. Ltd., USA). The Mg-Zn-Ca alloy specimens subjected to heat treatment and artificial aging for different duration of time were immersed in 45 mL of HBSS and incubated at 37°C for 72 h. The change in pH of the HBSS was measured as a function of immersion time.

2.2.5. Measurement of Corrosion Rate

The Mg-Zn-Ca alloy specimens subjected to heat treatment and artificial aging for different duration of time were immersed in 45 mL of HBSS and incubated at 37°C for 5 days. The HBSS was replaced every 24 h to suppress the change in concentration of the solution during the test. After 5 days of immersion, the Mg alloy specimens were rinsed in deionized water and ultrasonically cleaned using isopropanol for 5 min to remove the corrosion products attached to the specimen. The specimens were washed with distilled water and then dried. To examine weight changes in the specimen caused by corrosion, the specimen weight was measured before and 5 days after deposition. The corrosion rate was calculated using the following equation [23]: where is corrosion rate, is weight loss, is surface area of the specimen exposed to the etching reagent, and is exposure time.

2.2.6. Hemolysis Assessment

The test solution used for hemolysis assessment was blood samples collected from the coronary artery of rats, which were treated with heparin and diluted with saline solution in a 4 : 5 ratio (vol.%). Specimens having a dimensions of 16.4 mm diameter and 2 mm thickness from each experimental group were treated with 10 mL of saline solution for 30 min and the elute was used as the experimental group. 10 mL saline solution containing 0.1 wt.% of Na2CO3 and 10 mL saline solution were used as the positive and negative control groups, respectively. Each 10 mL solution of each group (experimental group and the positive and negative control group) was mixed with 0.2 mL of the prepared test solution and incubated at 37°C for 60 min. Subsequently, each mixture was centrifuged at 680 G for 3 min, and 200 μL of the supernatant liquid was collected. Then absorbance was measured at 540 nm using a medical enzyme analyzer (Emax, Molecular Devices, USA). The measured absorbance was used to assess the percentage hemolysis by the following equation [24]:

3. Results

The SEM image of the microstructure of the as-cast Mg-35Zn-3Ca alloy etched using Nital for 5 s and its chemical composition analyzed by EDX analysis are shown in Figure 1. The microstructure of the alloy is homogeneous and reveals no casting defects. The chemical composition of the alloy indicates the presence of Mg (63.82 wt.%), Zn (32.63 wt.%), and Ca (3.55 wt.%) as the major elements (Table 3), suggesting that the alloy conforms to the chemical composition requirements of Mg-35Zn-3Ca alloy.

The DSC trace of the as-cast Mg-35Zn-3Ca alloy specimen is shown in Figure 2. The liquidus and solidus temperatures are 361.84°C and 340.11°C, respectively. Based on this result, the suitable temperature for solution heat treatment is set as 310°C, which is slightly lower than the solidus temperature of the alloy. Accordingly, the as-cast Mg-35Zn-3Ca alloy is homogenized at 310°C for 48 h and quenched in deionized water at room temperature to prevent any microsegregation and residual stress during precipitation (designated as AT0). Some of the specimens were subjected to artificial aging heat treatment at 170°C for 2.5, 5.0, 7.5, and 10.0 h and designated as AT1, AT2, AT3, and AT4, respectively.

The SEM images of the surface microstructure of the heat-treated (AT0) and artificially aged (AT1, AT2, AT3, and AT4) Mg-35Zn-3Ca alloy specimens, etched in Nital for 5 s, are shown in Figure 3. Among all the specimens, the AT2 group exhibits the finest and most homogenous distribution of α and β phases. It is also evident from Figure 3 that, with increase in aging time, the size of α phase is increased with the disappearance of the boundary of β phase, resulting in a change in microstructure from a homogeneous to an inhomogeneous state. The X-ray elemental mapping of Mg, Zn, and Ca along with the secondary electron image of the Mg-35Zn-3Ca alloy is shown in Figure 4. It is evident that the darker regions of the microstructure are mostly the α-Mg matrix whereas the lighter regions are mainly the β-phase precipitate, which are relatively lower in Mg and largely consist of zinc and calcium.

The potentiodynamic polarization curves of heat-treated (AT0) and artificially aged (AT1, AT2, AT3, and AT4) Mg-35Zn-3Ca alloy specimens in HBSS are shown in Figure 5. The and values of all the specimens, calculated using Tafel extrapolation method, are compiled in Table 4. There is not much variation in the of all the specimens tested and all of them lie in the range of −1.20 to −1.29 V versus Ag/AgCl/KCl (sat’d). The values of the specimens that are heat-treated and artificially aged for different duration of time (AT1 to AT4) are relatively higher than those subjected only for heat treatment (AT0). Among the heat treated and artificially aged samples (AT1 to AT4), specimens of the AT2 group exhibit a lower . Based on the values, the corrosion resistance of the Mg-35Zn-3Ca alloys specimens subjected to only heat treatment (AT0) and heat treatment followed by artificial aging for different duration of time (AT1 to AT4) can be ranked as follows: AT0 > AT2 > AT3 > AT4 > AT1.

The change in pH of the HBSS in which the Mg-35Zn-3Ca alloys subjected to only heat treatment (AT0) and heat treatment followed by artificial aging for different durations (AT1 to AT4) are immersed, measured for a duration of 72 h, is shown in Figure 6. In general, all the specimens exhibit a rapid increase in pH during the first 6 h, beyond which the extent of increase in pH is decreased. Beyond 48 h of immersion the pH of HBSS reached a stable value for all the specimens studied. During the first 6 h of immersion, the Mg-35Zn-3Ca alloy specimens subjected to heat treatment alone (AT0) showed a rapid increase in pH when compared to those subjected to heat treatment followed by artificial aging for different duration of time (AT1 to AT4). Among the AT1 to AT4 group, specimens of the AT2 group showed a relatively less increase in pH, suggesting the lower extent of dissolution of the Mg alloy.

The corrosion rate of the Mg-35Zn-3Ca alloy specimens subjected to only heat treatment (AT0) and heat treatment followed by artificial aging for different duration of time (AT1 to AT4), calculated using the loss in weight after immersion in HBSS at 37°C for 5 days, is shown in Figure 7. It is evident that specimens subjected to heat treatment alone (AT0) show a higher corrosion rate (3.92 mg/(cm2 day)) when compared to those subjected to heat treatment followed by artificial aging for different durations (AT1 to AT4) (ranging from 2.32 to 3.73 mg/(cm2 day)). Among the AT1 to AT4 group, specimens of the AT2 group showed a lower corrosion rate of 2.32 mg/(cm2 day), which is significantly lower than those of other groups ().

In vitro hemolysis test is a well-known method to assess the biotoxicity in a short duration of time. The hemolysis ability of Mg-35Zn-3Ca alloy specimens subjected to only heat treatment (AT0) and, heat treatment followed by artificial aging for different duration of time (AT1 to AT4), is measured using the positive (10 mL saline solution containing 0.1 wt.% of Na2CO3) and negative (10 mL saline solution) as controls. The measured optical density values and the % hemolysis are compiled in Table 5. The variation in the % hemolysis of all the samples tested along with the negative and positive controls is plotted in Figure 8 for a better comparison and to understand the statistical difference among them. It is evident from Table 5 and Figure 8 that among all the specimens tested those belonging to the AT2 group showed the lowest value of % hemolysis (0.147%). However, when compared to other groups, the difference in % hemolysis is statistically insignificant ().

4. Discussion

Mg undergoes rapid corrosion in electrolytes containing chlorides, such as body fluid. Hence, it will be difficult to use Mg for implant applications. When Mg is exposed to air, the oxide layer formed on its surface offers only a partial protection. In aqueous solution, the formation of porous Mg(OH)2 layer could not offer much protection (equation (3)). In presence of chloride ions, it gets dissolved as MgCl2 while the aggressive nature of the medium causes severe corrosion of the Mg (equations (4) and (5)) [25]. The reaction sequences of Mg in aqueous medium in the absence (equation (3)) and presence (equations (4) and (5)) of chloride ions can be represented as follows [8]:

The corrosion resistance of Mg could be improved through alloying, formation of composites, surface treatment, and heat treatment [2531]. However, the alloying elements, the second phase materials used to prepare the composites, and the surface treated layer should not be harmful to the body when used as a biomaterial. Some of the commonly used Mg-based alloys contain aluminum that is harmful to the nervous system and osteoblasts while the RE-based Mg alloys could cause blood toxicity [17, 18]. In the present study, Mg-35Zn-3Ca alloy is prepared by adding zinc and calcium as alloying elements, both of which are nontoxic. Alloying magnesium with zinc enables precipitation hardening while alloying with calcium could help in grain refinement. The microstructure and corrosion resistance of the Mg-35Zn-3Ca alloy are evaluated after subjecting them to only heat treatment (AT0) and heat treatment and artificial aging for different duration of time (AT1 to AT4).

The microstructural examination after solution heat treatment and artificial aging heat treatment performed over a specific period of time showed that the α and β phases were finest and most homogenous in the AT2 group, which underwent an artificial aging heat treatment at 170°C for 5 h. Moreover, the grain boundary was broken and β phase was inhomogeneous in the AT3 group, which underwent artificial aging heat treatment for 7.5 h. Results also showed that the size of crystal grain became larger and the continuity of the crystal margin has disappeared in the AT4 group, which underwent artificial aging heat treatment for 10 h (Figure 3). Solution heat treatment is a process to prepare the supersaturated solid solution of unstable state; materials are generally softened by solution heat treatment. Previous studies have reported that artificial aging heat treatment following solution heat treatment of Mg alloy improved strength by refinement of the precipitate [29]. Meanwhile, it was reported that heat treatment changed the microstructure of alloy thereby affecting the corrosion behavior [32, 33].

In this study, potentiodynamic polarization studies, change in pH of HBSS, and loss in weight after immersion in HBSS were used to assess the corrosion resistance of the Mg-35Zn-3Ca alloy, both in as-prepared and heat-treated conditions. The results of the potentiodynamic polarization test, conducted to assess electrochemical corrosion behavior, showed the lowest corrosion current density of 1.011 μA/cm2 for the AT0 group, which underwent only the solution heat treatment. Among the groups that underwent artificial aging following solution heat treatment, the lowest corrosion current density of 3.641 μA/cm2 is observed for the AT2 group. The corrosion current density, which is closely associated with corrosion rate, is an intersection point at which the negative Tafel slope is extrapolated to corrosion potential by the Tafel extrapolation method [34]. This result showed that corrosion resistance of the Mg-35Zn-3Ca alloy is much improved for the group that underwent solution heat treatment alone than in groups that underwent artificial aging heat treatment following solution heat treatment. When immersed in HBSS, Mg undergoes anodic dissolution with the release of H2 gas and generation of OH ions, resulting in an increase in pH of the medium [35]. These reaction sequences can be represented as follows [36]:

Thus, the rate of corrosion of Mg/Mg alloys could be estimated from the change in pH of the corrosive medium (HBSS in this study) as a function of time. The change in pH of the HBSS measured as a function of time (Figure 6) showed that it is slowest for specimens in the T2 group when compared to other groups. This inference indicates that the corrosion resistance is increased as the microstructure of the surface gets finer and denser. The loss in weight of the specimens measured after 5 days of immersion in HBSS is 2.32 mg/(cm2 day) for the AT2 group, which is significantly lower than that of the other groups (). The loss in weight is 3.92 mg/(cm2 day) for the AT0 group, which is the highest among all the specimens tested. However, the difference between the specimens is insignificant when compared with the specimens in other groups except with the AT2 group (). Hence, it is clear that the specimens of the AT2 group have offered an excellent corrosion resistance, which is consistent with the results of change in pH. However, the results are inconsistent with those of the potentiodynamic polarization test. The fineness and homogeneity of the β phase largely determine the corrosion behavior, which is reflected in the pH variation and weight loss. These results supported that corrosion behavior of magnesium alloys with multiphase structure is changed according to corrosion potential difference, distribution, and volume fraction between micr-constituent phases [16, 37, 38]. Since these measurements are made for a long term, the ability of the β phase to act as a physical barrier for corrosion becomes effective [33, 39]. However, during the electrochemical study, since the test is performed over a short duration, the effect is not reflected [29].

The hemolysis test, which is conducted to assess the toxicity to the body, uses the principle that red blood cells are divided or necrotized when the hemoglobin is destroyed. The concentration of free hemoglobin is proportional to the number of the injured RBC and the volume of the fluid [24, 40]. The results of hemolysis test showed that the percentage of hemolysis (0.147%) in AT2 group is the lowest among all the samples tested and when compared with other groups, no statistical significance difference could be observed (Figure 8, Table 5) (). In general, if the hemolysis value is 5% or less, then it is considered as nonhemolytic according to ISO 10993-4:2002 [41]. Accordingly, the composition of Mg-35Zn-3Ca has low hemolysis and analogies good biocompatibility. In other words, homogenization and aging process of the Mg-35Zn-3Ca alloy can improve the physical properties and corrosion resistance of the material and it does not have a significant effect on the hemolytic toxicity.

5. Conclusions

The preparation of Mg-35Zn-3Ca alloy containing zinc and calcium, rather than aluminum or RE elements, as nontoxic and biocompatible alloying elements, followed by heat treatment and artificial aging of the alloy to control its rate of degradation, is reported. The study leads to the following conclusions.(1)DSC studies showed that the liquidus and solidus temperatures of the Mg-35Zn-3Ca alloy were 361.84°C and 340.11°C, respectively.(2)The microstructural features of the Mg-35Zn-3Ca alloy indicate the presence of α-Mg matrix and β phase after heat treatment and artificial aging for various duration of time. Among them, the AT2 group (heat treatment at 170°C for 5 h followed by artificial aging at 310°C for 48 h) exhibits the finest and most homogenous distribution of α and β phases.(3)Potentiodynamic polarization studies of the Mg-35Zn-3Ca alloy specimens indicate that the values of specimens that are heat-treated and artificially aged for different duration of time (AT1 to AT4) are relatively higher than those subjected only to heat treatment (AT0). Among heat-treated and artificially aged samples (AT1 to AT4), specimens of the AT2 group exhibit a lower . Based on the values, the specimens can be ranked as follows: AT0 > AT2 > AT3 > AT4 > AT1.(4)The change in pH of the HBSS measured for 72 h indicates that all the specimens exhibit a rapid increase in pH during the first 6 h, during which the Mg-35Zn-3Ca alloy specimens subjected to heat treatment alone (AT0) showed a rapid increase in pH. Beyond 48 h of immersion the pH of HBSS reached a stable value for all the specimens studied. Among the AT1 to AT4 groups, specimens of the AT2 group showed a relatively less increase in pH, suggesting the lower extent of dissolution of the Mg alloy.(5)Mg-35Zn-3Ca alloy specimens subjected to heat treatment alone (AT0) showed a corrosion rate of 3.92 mg/(cm2 day) than those subjected to heat treatment followed by artificial aging for different durations (AT1 to AT4), which ranges between 2.32 and 3.73 mg/(cm2 day). Among the AT1 to AT4 groups, specimens of the AT2 group showed a lower corrosion rate of 2.32 mg/(cm2 day), which is significantly lower than those of other groups ().(6)In vitro hemolysis test showed the lowest hemolysis percentage value of 0.147% for the AT2 group, but the differences were statistically insignificant when compared with that of other groups (). The hemolysis was 5% or less for all groups, which indicates no significant hemoglobin destruction, which suggests that Mg-35Zn-3Ca alloy has no cytotoxicity.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Il-Song Park contributed to this work and is considered as joint first author.

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) Grant funded by the Korea government (MSIP) (no. 2011-0028709 and no. 2014R1A4A1005309) and Regional Strategic Industry project (2013-R0002274).