#### Abstract

The behaviors of viscosity and thermal expansion for different compositions of bioactive glasses have been studied. The effect of phosphorous pentoxide as a second glass former in addition to silica was investigated. Consequently, the nonlinear behaviors of viscosity and thermal expansion with respect to the oxide composition have been modeled. The modeling uses published data on bioactive glass compositions with viscosity and thermal expansion. -regression optimization technique has been utilized for analysis. Linear and nonlinear relations are shown to establish the viscosity and thermal expansion coefficients associated with oxide components of the glasses under study. The modeling allows the calculation of viscosity for a given temperature and, accordingly, the fusion temperature of these glasses along with the coefficient of thermal expansion. The established model relations also suggest first- and second-order phosphorus-alkali and alkaline earth oxides interaction which is reflected on the model coefficient that calculates viscosity and thermal expansion.

#### 1. Introduction

In a review of the importance of comprehensive information on glass and melt properties both for the glass science and technology by Mazurin [1], it is shown that the number of publications per year concerning glass property measurements have decreased quite considerably after 1975. Moreover, the frequency of some quite important but labor-consuming property measurements decreases steadily. Besides that, these measurements are quite time and money consuming; thus, the specialist on glass science and technology needs to predict the glass properties via calculation of it. This task must be based on modeling methods which use experimental property data and calculate glass properties including newly designed glasses.

The changes in the composition area of glass often lead to the change of the reliable modeling, that is, different models may be needed. Finding the most reliable model is possible only by comparing the results of property calculations with experimental data for glasses belonging to a system of interest. Next, the model with a minimal root-mean-square error is selected. After that, the model is used for future calculations of glass properties. At present, this is the best way to obtain a property value that is as near to the true one as possible. This explains why modeling of glass properties is an enduring active field in glass science and technology [2].

Viscosity at given temperature and the coefficient of thermal expansion CTE for bioactive glasses are the glass properties which are focused on in this study. Bioactive glasses have the property of forming an apatite layer; thus, when the implants is coated with bioactive glass, it can chemically bond to bone [3]. However, if there is a large difference in coefficient of thermal expansions between the coating and the implant, the bioactive glass will crack resulting in incomplete coating. A second problem is that the fusion temperature is needed to be well estimated for efficient cast of the glass coating to the implant. The fusion temperature should be low enough to avoid phase transformation and degradation of the mechanical properties of the implant. Therefore, scientific efforts must be carried out to achieve glass compositions that are suitable in view of fusion temperature and coefficient of thermal expansion [4]. The glass designer intends to assemble bioactive glass with fusion temperature below 900°C to avoid possible phase transformation of for example, Ti-alloys; also, the value CTE should be around which should match that of Ti-alloys [5]. In addition, the trend in bioactive glass compositions is to set P_{2}O_{5} content to 6% with high CaO content to ensure the formation of apatite layer after implantation that is responsible for the biocompatibility of the glass.

In general, viscosity and thermal expansion for bioactive glass shows nonlinear behavior with glass composition. Therefore, any proposed model should express effective composition-related terms that describe these properties for bioactive glass. Consequently, the model calculation can be used in design technology in such glass category [6].

The topic of predicting the viscosity and thermal expansion of bioactive glass is still vigorous for different sets of bioactive glass compositions [6, 7], and the modeling study based on formerly published data is still on focus in the applied ceramics community [8–10]. This work presents modeling methodology, results, and discussion of the modeled bioactive class along with the comparison on the effect of CaO and MgO content on the properties of the bioactive glass.

#### 2. Composition-Viscosity and Thermal Expansion Model

The variation of viscosity with the absolute temperature is expressed as the following [11]: The parameter is constants, is the gas constant, and is the activation energy of the viscous process.

The natural logarithm of (1) readsas the following: The above equation shows that and − are linked linearly. Accordingly, values of ln can be calculated for a given set of viscosity at temperature and for a range of glass compositions.

is as function of glass composition as well as the coefficient of thermal expansion CTE [12]:

is the weight percentages wt% of the th component composing the glass and or is the coefficient associated with the th glass component. The terms () or () can also include nonlinear combination of more than one glass component. Finally, *n* represents the total number of terms.

The final step is to explicitly express the composition-viscosity or thermal expansion model that is specific to the bioactive glass composition set under study. In this study, the effort to present such a model is enlightened by the discussion of Kuppinger and Shelby [13], that is, inspecting possible interaction between ions composing the glass which may lead to the nonlinear relations of the bioactive glass composition with and CTE. In this study, phosphorous-alkali and alkaline earth effect is found. So, the suggested composition-glass property model is stated as below:

, , , and are first- and second-order functions of phosphorous-alkali and phosphorous-alkaline earth content:

Equations (4a)–(6b) are the model equations presented in this study that calculate and CTE as function of the glass compositions.

#### 3. Modeling Technique

Table 1 shows published bioactive glass composition data with their thermal properties [5]. In this table, Ts represent the temperature at constant viscosity level of *μ* = 10^{6.6} Pa·s, which corresponds to the Littleton softening temperature. Tg represents another constant viscosity level of *μ* = 10^{12} Pa·s, which is in the range of the expected glass transition temperature (*μ* = 10^{11}10^{12.3} Pa·s).

Utilizing (2), the two unknowns ( and ln) can be found by substitution of Ts and Tg in Table 1. Figure 1 is a plot of against –ln which shows linear relation connecting them. The figure also shows the equation of the least-squares line passing through the data points. The reference values of and ln found utilizing Table 1 is used together with the glass compositions in -regression [14] calculations. Accordingly, calculated glass properties ( and CTEc) are found as function of their composition.

The assumption of linear additive behavior is firstly examined. This assumption yields to large differences between reference experimental and calculated glass properties. These differences remain at large values even when several combinations of the wt% of oxide constituents are taken into account. The differences between the reference experimental and calculated glass properties are reduced only when phosphorous-alkali and alkaline earth effect as shown in (4a)–(6b) is taken into account. The calculations astonishingly show that the nonadditive relations are equally needed for modeling of both and CTE. The final root-mean-square error for and CTE calculations are 26.47 and 0.23, respectively. These are only around 2% of the reference experimental value for and CTE.

Table 2 shows the reference experimental , reference experimental ln , calculated via -regression, calculated ln (utilizing calculated and the least-squares relation shown in Figure 1), reference experimental CTE, and calculated CTE via -regression. The final coefficients and that obtained by -regression to calculate and CTE for bioactive glasses are shown in Table 3.

To find the temperature at a given viscosity value (an inverse problem), it needs to calculate and ln by the above methodology, then (2) can be solved for . This can be accomplished using any numerical technique such as the “finite step method.” Finally, setting to 20 Pa·sec; the fusion temperature for a bioactive glass can be calculated [15].

To conclude, Figure 2 is a way to show how and CTE respond to CaO and Na_{2}O content. and CTE are calculated by setting K_{2}O, MgO, and Na_{2}O to 3, 6, and 9 wt%, respectively, and let SiO_{2} varies with CaO to make total contents to 100 wt%. The same is done when varying Na_{2}O by setting K_{2}O, MgO, and CaO to 3, 6, and 11 wt%, respectively. The figure is presented to show an example of the nonlinear behavior of and CTE with CaO and Na_{2}O. Yet, and CTE are linear with other oxide contents as shown in (4a)–(6b)

**(a)**

**(b)**

**(c)**

**(d)**

#### 4. Discussion

Table 3 shows that SiO_{2} has a trend to increase and to decrease CTE. At the same time, P_{2}O_{5} tends to decrease both and CTE. MgO shows similar behavior with SiO_{2}; also, K_{2}O behaves similar to P_{2}O_{5} with different values of their coefficients. Figure 2 describes the behavior of CaO and Na_{2}O more obviously. Na_{2}O shows a trend to increase both and CTE in nonlinear manner. The effect of Na_{2}O is distinguished from that of K_{2}O although they are both alkali oxides. Similarly, Figure 2 shows that CaO is decreasing and increasing CTE in nonlinear fashion. This effect of CaO is also distinguished from that of MgO although they are both alkaline earth oxides.

The aforementioned different type of contribution for both Na_{2}O and CaO to the viscosity and thermal expansion compared with that of K_{2}O and MgO, respectively, may be understood in view of that Na2O and CaO content is generally higher than that of K_{2}O and MgO in the bioactive glass as seen in Table 1. Consequently, large ion concentration of Na and Ca give a chance for P-Na and P-Ca interactions. These interactions are seen in a way that P_{2}O_{5} is modifying the role of Na_{2}O and CaO to that observed in Figure 2 and reflected formally in (9)–(12). Nevertheless, those interactions on a relatively low-concentration level of K_{2}O and MgO do not influence the properties apparently. This explanation is illuminated by the discussion of Kuppinger and Shelby [13] that suggest short range ion-ion interaction in the glass medium that yields nonlinear relations of glass viscosity and thermal expansion with the oxide content. They study sodium potassium borate glasses with high sodium and potassium oxides contents which led to ion-ion interaction between the alkali ions. Nevertheless, in this study the high sodium and calcium oxide contents do not led to mutual interaction between them, but phosphorous-alkali and alkaline earth interaction takes place. The acidic nature of P_{2}O_{5} may explain the interaction with alkali and alkaline earth oxides of the basic nature.

It is worth to mention that multilayer coating is one of the experimental experiences to avoid crack generation which originates from CTE mismatch [7]. However, multilayer coating may lead to other problems such as thick coating. As a direct application of the modeling presented in this study, a bioactive glass can be designed with suitable fusion temperature and CTE passing up lengthy experiments and cost.

#### 5. Conclusions

(1)Model relations of viscosity and thermal expansion with bioactive glass compositions have been obtained. The model allows the predictions of viscosity and thermal expansion for bioactive glass necessary in the field of glass science and technology.(2)The model relations are discussed for the nonlinear behavior of viscosity and thermal expansion (deviation from additivity) with Na_{2}O and CaO content. It is suggested that the deviation from additivity originates from short-range interaction between alkali and alkaline earth ions with phosphorus as formally presented in the model.