Next Article in Journal
Structural, Optical, Charge-Transport, and Dielectric Properties of Double-Perovskite La2Co1−zFezMnO6 (z = 0, 0.2–1.0)
Next Article in Special Issue
Three-Dimensional-Printed Molds from Water-Soluble Sulfate Ceramics for Biocomposite Formation through Low-Pressure Injection Molding
Previous Article in Journal
Effect of the Preparation Conditions on the Catalytic Properties of CoPt for Highly Efficient 4-Nitrophenol Reduction
Previous Article in Special Issue
Mn-Containing Bioactive Glass-Ceramics: BMP-2-Mimetic Peptide Covalent Grafting Boosts Human-Osteoblast Proliferation and Mineral Deposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of Multicomponent Fluoridated Borate Glasses through a Design of Mixtures Approach

1
IR Scientific Inc., Halifax, NS B3H 0A8, Canada
2
Department of Applied Oral Sciences, Dalhousie University, Halifax, NS B3H 4R2, Canada
*
Author to whom correspondence should be addressed.
Materials 2022, 15(18), 6247; https://doi.org/10.3390/ma15186247
Submission received: 22 July 2022 / Revised: 1 September 2022 / Accepted: 2 September 2022 / Published: 8 September 2022
(This article belongs to the Special Issue Bioceramics and Bioactive Glass-Based Composites)

Abstract

:
Due to their enhanced dissolution, solubility and reaction speed, borate glasses offer potential advantages for the design and development of therapeutic ion-release systems. However, the field remains poorly understood relative to traditional phosphosilicate and silicate bioglasses. The increased structural complexity and relative lack of published data relating to borates, particularly borofluorates, also decreases the accuracy of artificial intelligence models, which are used to predict glass properties. To develop predictive models for borofluorate networks, this paper uses a design of mixtures approach for rapid screening of composition–property relationships, including the development of polynomial equations that comprehensively establish the predictive capabilities for glass transition, density, mass loss and fluoride release. A broad range of glass compositions, extending through the boron anomaly range, were investigated, with the inclusion of 45 to 95 mol% B2O3 along with 1–50 mol% MgO, CaO and Na2O as well as 1–30% KF and NaF. This design space allows for the investigation of the impact of fluorine as well as mixed alkali–alkaline earth effects. Glass formation was found to extend past 30 mol% KF or NaF without a negative impact on glass degradation in contrast to the trends observed in phosphosilicates. The data demonstrates that fluoroborate materials offer an exceptional base for the development of fluoride-releasing materials.

1. Introduction

The principal of using therapeutic ion release from bioactive glasses enables the engineering of versatile and tailorable biomaterial platforms to address therapeutic needs at the local tissue level. Such glassy materials react with an aqueous environment to release ionic constituents capable of directing local cellular activity, while concurrently enabling improved processing, sterilization and shelf life relative to traditional drug loading approaches [1,2]. In this regard, antibacterial [3], anti-inflammatory, angiogenic [4] and osteogenic [5,6] materials have all been developed through the inclusion of therapeutic ions into glassy materials [7]; though, to date, this has largely been associated with silicate and phosphosilicate glass chemistries. Increasingly, interest in borate glass chemistries as the basis for controlled release bioactive glasses has become the focus of attention in the literature.
Commercially, the use of borate glasses (in fiber form) for the management of skin ulcers has proven the utility of borate-based bioactive glasses as therapeutic agents to promote tissue healing and repair [8]. Clearly, the selection of therapeutic ions to be incorporated into a soluble glass network depends on the intended use and indications for use. Two therapeutic inorganic ions of interest for hard tissue engineering in dentistry are magnesium and fluorine.
Magnesium is known to play a role in the regulation of bone mineral density; both through direct effects on bone cells and indirectly through its impact on bone signaling pathways [9,10]. Along with its ability to alter cell interactions, magnesium ions are capable of substituting into the mineral phase of bone and teeth, forming a magnesium substituted hydroxy apatite and altering the crystal lattice [11,12].
Magnesium substitution has been demonstrated to result in smaller crystallite size, imparting both a whiter appearance due to increased light scattering and resulting in an increase in microhardness of enamel [13,14]. Fluoride ions, in turn, are arguably the most well recognized therapeutic ion, with a long history of use in the prevention of tooth decay and have also shown beneficial effects for the maintenance of bone mineral density. For the design of dental materials, optimization of the loading and release of fluoride ions from bioactive glasses has gained both research and commercial interest, both for the prevention of dental decay and acid erosion and for the treatment of dentin hypersensitivity [15,16,17,18,19].
When considering the design of glasses for the release of therapeutic inorganic ions, borate glass networks are an appealing choice due to their high solubility and increased rate of ion release [20]. The substitution of boron oxide into the original bioglass formulation (often referred to as 45S5) in place of the original phosphosilicate glass network for example has been reported to result in a five-fold increase in reaction rate [21]. Furthermore, glasses based on boron networks have been reported to undergo bulk dissolution, proving them a useful base of the design of completely soluble biomaterials [22].
The design of borate glasses to meet specific material characteristics, however, is less straightforward than those of their silicate counterparts due to the ability of boron to switch between three- and four-fold coordination [23]. In contrast to the behavior of silicate glasses, the impact of modifier oxide addition to a borate glass impacts glass properties, such as the density, glass transition and network connectivity in a non-monotonic function, with a localized maxima/minima around ~30% modifier addition [24].
As the reactivity of bioactive glasses is commonly thought to be attributable to their network connectivity, with the best bioactivity around a network connectivity of two [25], the nonlinear variations in structural changes have a significant impact on glass dissolution and bioactivity. Further complicating the design of borate glasses, the mixed alkali and mixed alkali earth effect results in changes to the location of this minima/maxima when the glass complexity is increased beyond a single modifier binary glass [26,27].
As such, while borates are promising materials for the development of ion releasing bioactive glasses, they present an increased challenge for the prediction of composition–structure–property relationships. Difficulties in the prediction of structural drivers of glass reactivity in borate glasses are further complicated when considering magnesium and fluorine-containing borate glasses. Binary magnesium borate glasses have a narrow range of glass formation (45–55% MgO) when formed using traditional melt quench techniques [24,28]. This narrow glass-forming window has been attributed to magnesium’s high field strength and low ionic radius.
The literature reports are in conflict regarding the role that magnesium plays in a borate glass network, identifying its oxygen bonding as ionic [29] (indicative of a network modifier) and covalent [28] (indicative of network formation). Fluoride addition, in turn, is known to cause the formation of fluoride rich cation domains, removing glass modifiers from the network.
As the network connectivity of borate glasses varies nonlinearly with oxide modifier addition, it can be expected that the removal of modifying cations into fluoride clustering domains would also impart a nonlinear change in glass structure and properties, making glass behaviors less simple to predict. Some understanding of the impact of fluoride salt addition on the properties of borate glasses can be gained from the field of optical glasses, where fluoroborate glasses have been investigated due to their improved ionic conductivity relative to oxide glasses [30,31,32,33,34].
While some structural information can be gained from this work, these studies did not investigate the potential for fluoride release from the materials, or the reactivity of the glass. Furthermore, as there is little consensus on the impact of borate glass structure as a modulating factor for dissolution, extrapolation of composition–structure–property information from existing publications is not feasible.
For example, investigations on the effect of fluoride salt inclusion on the bioactivity of glasses in a B2O3–CaO–Na2O–CaF–NaF system has been conducted, however the boron content of the glasses was kept constant, limiting our understanding of the impact of the fraction of glass network formers to modifiers [35,36]. The impact of fluoride addition on the degradation, and thus fluoride release, or a borate oxyfluoride glass system is an underexplored research area hindered by our lack of a priori and a posteriori knowledge in the larger field of borate glass networks.
Compounding these issues, borate glasses have historically received little research interest, with less than 2% of published glass science literature from 2007 to 2013 focused on borates [37]. This historical under investigation has left the field at a significant disadvantage in terms of collective knowledge. This disadvantage grows exponentially in the age of machine learning and “materials by design” where the historical knowledge of glass materials can be compiled and interpolated to allow for more refined and rapid identification of useful compositions [38,39].
Such a disadvantage is again compounded with the considerations of less studied glass modifiers, such as magnesium and fluoride, where multiple structural roles are possible, increasing the complexity of the structural drivers of material properties. To exemplify this challenge, when using artificial intelligence techniques to model fluoride glass density, Ahmmad et al. found that, while a coefficient of determination (R2) of 0.980 was achievable for their dataset of 1258 fluoride glass compositions, it decreased to 0.792 for high-borate glasses [40].
A similar study aiming to model the density of oxide glasses found R2 values to decrease from 0.84 when modeling binary systems to 0.65 for multicomponent systems (defined as 5 or more components), emphasizing the increased difficulty in predicting multicomponent glass properties [41]. As such, alternative glass design approaches must be used to investigate borate glasses, especially in the context of less-studied glass components, such as fluoride or magnesium.
To better understand the complex relationship(s) governing the composition–property relationships in the understudied field of multicomponent fluoridated borate glasses, a design of mixtures approach is presented in this work. This approach investigates the impact of variations in boron oxide, alkali earth oxide, alkali oxide and fluoride salt loading in a multifactorial design of mixtures methodology. The design of mixtures is a form of experimental design in which all variable components sum to unity (such as the composition of a glass) and can be expressed as relative weights; this feature allows for a mixture of n components to be simplified down to n − 1 variables.
The compensatory changes required to balance the mixture to unity complicate the use of intuitive, step wise experimental design (such as the replacement of component 1 with component 2, then 3 and then 4). For these reasons, the selection of experimental compositions is typically made though a statistics software program to allow optimal spacing for the generation and assessment of composition-property surface response models. This allows for the development of robust and objectively verifiable composition-property relationships over a broad range of compositions, while minimizing the number of compositions to be studied.
The compositional range to be investigated covers boron oxide contents of 45–95%, allowing for a broad understanding of the impacts of modifier additions up into the inverted glass range. Minimum inclusion levels of 1% were chosen for all compositions to maintain model focus on multicomponent mixtures (i.e., to prevent the inclusion of binary of ternary compositions), requiring a 95% maximum for B2O3.
Oxide modifier ranges were selected between 1 and 50% for oxides magnesium, calcium and sodium, while sodium and potassium fluoride inclusions were limited to 1–30% due to experimental difficulty with high fluoride compositions (unpublished laboratory data) This study aims to serve as an initial screening tool to allow for the identification of regions of interest that merit more in-depth study for their potential use in biomaterials design and development for a variety of applications.

2. Materials and Methods

2.1. Experimental Design

Composition–property relationships were investigated for borate oxyfluoride glasses in the region of 45–95% B2O3, 1–50% MgO, 1–50% CaO, 1–50% Na2O, 1–30% NaF and 1–30% KF using a design of mixtures approach. Design Expert Software (Stat-Ease, Minneapolis, Minneapolis, MN, USA, Version-12.0.1.0, https://www.statease.com/software/design-expert/ accessed on 6 August 2022) was used to develop an I-optimal quadratic model based on a 6-component mixture, with the design constraints as shown in Table 1. From the model, 31 randomised “runs” corresponding to a glass compositions (Table 2) were identified for investigation.

2.2. Glass Synthesis

The 31 glass compositions as specified in Table 2 were made through traditional melt quench methods in the molar ratios shown above. Briefly, appropriate weights of reagent grade modifiers (sodium carbonate, calcium carbonate, magnesium oxide, sodium fluoride and potassium fluoride) were combined with appropriate weights of 98% pure boron oxide (Alfa Aesa), sealed and mixed to homogeneity in HDPE containers using a rotary dry mixer for 1 h (Glas-Col model 099A RD9912). Reagents were then packed in 100 mL 90% Pt/10% Rh crucibles (XRF Scientific) and placed in a Carbolite Gero RHF 14/3 furnace at 600 °C.
The furnace program was set to dwell for 60 min and then ramp at 20 °C/min to 1200 °C and hold for a further 60 min. Following the final dwelling period, glasses were quenched between two stainless steel plates and stored in desiccated conditions until further processing. Glasses were ground using a two-step process, first using a manual stainless steel glass crusher before further milling with a planetary micro mill using zirconia accessories (Pulverisetter 6 Classic Line, Fritsch, Idar-Oberstein, Germany). The ground glass powder was then sieved with an ASTM E-11 compliant sieve to retain particles of <25 um. All materials were stored in desiccated conditions until the time of further analysis.

2.3. X-ray Diffraction

X-ray diffraction was performed using a Bruker D2 diffractometer with a Cu source and Lynx-Eye XE linear array detector. Glass powder samples were analyzed in the region of 10° ≤ 2 θ ≤ 60° with a step width of 0.03 and a step time of 2 s. XRD spectra were analysed using Bruker Diffrac.Eva to assess the relative fraction of crystalline material in the samples (as the ratio of crystalline peaks to amorphous halo in the scan).

2.4. Helium Pycnometry

The glass density was measured using an Accupyc 1340 Helium pycnometer (Micromeritics, Norcross, GA, USA) with a 1 cm3 sample insert. The pycnometer chamber volume was calibrated with a traceable reference standard as per the manufacturer’s instruction prior to each use. For each glass sample, approximately 0.9 g of ground glass powder was firmly packed into the sample cup and analysed as the average of 10 fill and purge cycles using high purity helium gas (99.995% He, Praxair, Canada).

2.5. Differential Scanning Calorimetry

Differential scanning calorimetry was performed using a Pegasus F404 DSC (Netzsch, Selb, Germany), with a silicon carbide furnace in an argon environment. Approximately 25 mg of ground glass was analyzed in a platinum rhodium crucible, fitted with a lid. Glass samples were heated from 50 to 900 °C at a rate of 10 °C/min, with 10 data points collected for each degree (100 data points/min). Thermal curves were assessed to identify the glass transition temperature (as Tonset, Tinflection and Tfinal) and first crystallization temperature using the Netzsch Proteus software.

2.6. Mass Loss Assessment

The mass loss assessment was performed on all glass compositions in triplicate at timepoints of 1, 12 and 24 h. For each glass sample, 0.1 g of glass powder was placed in a pre weighed 15 mL HDPE centrifuge tube with 10 mL of TRIS-buffered saline solution (BioUltra, Sigma Aldrich, Oakville, ON, Canada), capped and placed horizontally on a shaking incubator at 37 °C and agitated at 2 Hz until the desired time point.
Following incubation, the samples were centrifuged for 15 min and 300 RCF to remove residual glass from the solution. Tubes were carefully decanted to remove the supernatant, which was capped and stored at 4 °C until the time of fluoride release analysis. Test tubes with glass residuals were dried in an oven at 70 °C until constant weight was achieved (approximately 2 weeks) to determine residual glass mass and calculate the % of weight loss.

2.7. Fluoride Release

The fluoride release was assessed using the decanted supernatants from the mass loss assessment. The concentration of F released into solution was measured using an Accumet AB250 pH/ion selective electrode meter equipped with a fluoride-ion-selective electrode (Accumet AB250). The ion-selective electrode was calibrated using a series of NaF standard solutions ranging from 0.1 to 10,000 ppm (NaF 0.1 M F, Sigma Aldrich, Oakville, ON, Canada). All samples were mixed with TISAB III as per the manufacturer’s instruction prior to analysis.

3. Results

3.1. Glass Synthesis

All compositions were successfully melted, with the exception of composition #20, which did not form a liquid melt at 1200 °C. Composition 20 was excluded from further analysis, and no attempt to remelt the composition at a higher temperature was made due to the potential impacts of variable thermal processing on glass properties. Multiple compositions did not form homogeneous clear glasses, evidenced both through nonuniform and cloudy appearance (denoted by an asterisk next to composition ID in Table 3).

3.2. XRD Analysis

XRD analysis of the 30 successfully melted oxide blends revealed varying degrees of crystallinity, ranging from 1.4 to 88% (Figure 1, Table 3). A reduced quadratic model was developed to understand the influence of glass components on melt crystallization, as shown in Table 4. While the R2 of 0.88 suggests a strong correlation between the modeled and true crystallization values, a large difference between the adjusted R2 and predicted R2 suggests an impact of outliers on the model, identified as composition 1.
The removal of outliers did not result in an improvement in the model’s predictive ability; however, the adequate precision result of 20.10 suggests that the model is useful in providing an estimate of the crystallinity of materials in this composition space. The factors that had a positive effect on crystallinity (favoured crystallization) in ranking of degree of impact were MgO > B2O3 > CaO > KF, factors that had a negative effect on crystallinity (favoured amorphous glass formation) in ranking of degree of impact were Na2O > NaF > B2O3 × MgO > MgO × KF > MgO × CaO > MgO × NaF > MgO × Na2O, as represented in the 3D surface plot of Figure 2a.
To better visualise the impact of single glass components with opposing impacts as combination factors a graphical correlation table is presented in Figure 2. For example, it can be visualised that the impact of KF additions imparts a positive linear effect, along with a negative second order effect for MgO × KF, resulting in an overall negative correlation to crystallinity (sixth row, seventh column of Figure 3).

3.3. Helium Pycnometry

The densities measured ranged from 1.88 to 2.70 g/cm3 as shown in Table 3 (with the exceptions of 1, 6 and 13, which yielded insufficient material for analysis). The density data was fitted to a reduced quadratic model as presented in Table 4. All glass component terms had a positive effect on glass density, with the effects ranked from most to least significant as CaO > NaF > KF > MgO > Na2O > B2O3 > B2O3 × Na2O > B2O3 × CaO. The only equation term with a negative impact on the glass density was CaO × Na2O, as visualised in Figure 1b. The molar volume ranged from 22.64 cm3/mol to 36.77 cm3/mol, as shown in Table 3.
The molar volume was fitted to a quadratic model as presented in Table 4. All individual glass components had a positive effect on molar volume, along with the Na2O × CaO interaction term. The relative magnitude of the terms positively impacting molar volume was B2O3 > CaO > Na2O > KF > MgO > NaF > CaO × Na2O. The terms that had negative effects on molar volume (caused a contraction of the glass network) were all interaction terms, their ranking from greatest negative impact to least impact is B2O3 × Ca2O > B2O3 × MgO > B2O3 × Na2O > B2O3 × NaF > B2O3 × KF, as visualised in Figure 2c. To better visualise overall impact of glass components, a graphical correlation table is presented in Figure 3.

3.4. Differential Scanning Calorimetry

The Tg onset values for the glasses investigated ranged from 346.3 to 607.8 °C (Figure 4, Table 3) (exceptions of compositions 1, 13, 17 and 20, which yielded insufficient material for processing and analysis) and was fitted to a linear model as shown in Table 4. All glass components had a positive effect on the Tg onset (increased Tg) with the effects ranked from most to least significant as CaO > MgO > B2O3 > NaF > Na2O > KF as visualised in Figure 2d. Tg inflection ranged from 357.1 to 618.4 °C (Table 3) and was best fitted to a reduced quadratic model. The glass components with a positive effect on Tg inflection, ranked from greatest to least impact were B2O3 > NaF > MgO > KF > CaO > B2O3 × CaO > B2O3 × Na2O > CaO × NaF > B2O3 × MgO > MgO × KF, factors with a negative impact on Tg inflection ranked from greatest to least impact were Na2O > CaO × Na2O > MgO × Na2O, as visualised in Figure 2e.
The glass transition as measured by Tg final ranged from 354.7 to 617.4 °C and was best fitted by a reduced quadratic equation. The factors that favoured an increase in Tg final, in ranked form greatest to least impact were MgO > B2O3 > NaF > KF > CaO > B2O3 × CaO > B2O3 × Na2O> CaO × NaF. The factors that favoured a decrease in Tg final, in order of greatest to least impact were Na2O > MgO × Na2O > CaO × Na2O > MgO × NaF > MgO × CaO as visualised in Figure 2f. To better visualise the overall impact of the glass components, a graphical correlation table is presented in Figure 3.

3.5. Mass Loss Assessment

The mass loss assessment was performed on all glasses with the exceptions of 1, 6, 9, 13, 17 and 20 due to insufficient material yield for processing and analysis. The mass loss was from 40.5% to 99.8%, 55% to 99.5% and 59% to 99.4% following 1, 12 and 24 h, respectively (Table 5). The mass loss for all timepoints was best modelled following a logit transformation with lower and upper limits of 0 and 100% (Table 4). The mass loss following 1 h was best fitted to a reduced quadratic model, and the factors favouring increased mass loss ranked from greatest to least impact were CaO > MgO > B2O3 > Na2O.
The factors that favoured a decrease in mass loss were NaF > KF > B2O3 × CaO > B2O3 × MgO > CaO × Na2O as visualised in Figure 5a. The mass loss following 12 h was best fitted to a linear model, with the factor favouring increased mass loss ranked from greatest to least impact being Na2O > B2O3 > NaF. The factors that resulted in a decrease in mass loss were CaO > MgO > KF as visualised in Figure 5c. The mass loss at 24 h was best modeled using a reduced quadratic model, factors that had a positive impact on, ranked from greatest to least impact were B2O3 > Na2O > NaF > KF, and the factors that favoured a decrease in mass loss, ranked from greatest to least impact were CaO > MgO as visualised in Figure 5e.

3.6. Fluoride Release

The fluoride release was assessed for all glasses with the exceptions of 1, 6, 9, 13 and 20 due to insufficient materials for analysis. The fluoride release was measured at timepoints of 1, 12 and 24 h and ranged from 14.4 to 1918.14 ppm at one hour, 16.6 to 946 ppm at 12 h and 16.7 to 1064.5 ppm after 24 h (Table 5). The fluoride release was best modelled following a logit transformation with upper and lower bounds of 0 and 2000 ppm (Table 4). The glass components favouring an increase in fluoride release at 1 h, in order of greatest to least impactful were Na2O > NaF > KF> Na2O × KF, components favouring a decrease in fluoride release were B2O3 > CaO > MgO > CaO × NaF as visualised in Figure 5b.
The fluoride release following 12 h of incubation in TRIS-buffered saline was fitted to a similar reduced quadratic model following a logit transformation with upper and lower bounds of 0 and 2000 ppm. The factors that favoured an increased in F release following 12 h in order of greatest to least impact were Na2O > KF > B2O3 × Na2O > MgO × Na2O > Na2O × NaF, factors that negatively impacted fluoride release at 12 h were, in order of greatest to least impact NaF > B2O3 > CaO > MgO > NaF × KF > B2O3 × Na2O as visualised in Figure 5d.
The fluoride release following 24 h of incubation in TRIS-buffered saline was fitted to a similar reduced quadratic model following a logit transformation with upper and lower bounds of 0 and 2000 ppm. The factors that favoured an increased in F release following 12 h in order of greatest to least impact were KF > Na2O > B2O3 × NaF > Na2O × NaF, and the factors that favoured a decrease in F release at 24 h, from greatest to least impact were NaF > B2O3 > CaO > MgO > CaO × KF > MgO × KF as visualised in Figure 5f. To better visualise overall impact of glass components, a graphical correlation table is presented in Figure 3.

4. Discussion

The prediction of bioactive glass reactivity is frequently correlated to the network connectivity in phosphosilicate glasses. Stemming from Hill’s alternative view on bioactive glass dissolution, it is proposed that a network connectivity around 2 allows the glass to behave as a tangled linear polymer, allowing for easy release of glass modifiers and the formation of the silica gel layer as a precursor to calcium phosphate deposition [25]. Similar correlations between borate glass network structure and glass dissolution have been found, suggesting that network connectivity is a strong predictor of the rate of glass dissolution in borate systems as well [22].
As such, structural modifications are a natural design consideration for the development of novel bioactive glasses. Additions of modifying oxides, provide for a linear decrease in network connectivity, which can be empirically calculated in silicate glasses, providing strong predictors of glass properties. The addition of fluoride to the network of a phosphosilicate has the opposing effect, removing network modifiers from the glass structure and increasing network connectivity, again with a linear effect. For borate glasses, direct manipulation of the glass network structure through modifier additions is less straightforward than in silicate glasses.
In borate glasses, network modifiers may act to generate a transformation from three- to four-fold coordination of boron, or to generate non-bridging oxygens. The balance between these two structural roles is dependent on the modifier and the amount added to the formulation chemistry. While some rules governing the formation of tetrahedrally-coordinated boron and non-bridging oxygen have been generated, they do not offer a consensus, and they cannot be extrapolated across all borate glasses [42].
For example, while many these models treat glass modifiers equally, it is recognized that the fraction of tetrahedrally-coordinated boron is dependent also on the modifier chosen, as demonstrated by the comparison of binary lithium and caesium borate glasses, where the lithium glasses showed strong agreement to models previously developed, but the caesium glasses had substantially fewer tetrahedrally-coordinated boron (46% vs. 43% tetrahedrally-coordinated boron at 41% modifier addition for lithium and caesium borates, respectively) [43].
Following from this, the effects of fluoride additions cannot be easily predicted in borate glasses. If, as suggested in some literature, fluoride ions serve to remove cationic modifiers from the glass network, they can serve both to decrease the number of non-bridging oxygens or to decrease the number of tetrahedral coordinated boron, having opposing effects on the overall network connectivity of the glass. For this reason, the selection of a region of interest for investigation in a borate glass network is much less easily identified than their silicate counter parts (where centering the network connectivity around 2 would be a reasonable guess to begin investigation).
In magnesium oxyfluoride borate glasses, yet another level of structural complexity is added, which further obscures predictions of composition-structure property relationships. Magnesium oxide within a borate glass network has been shown to act both as a glass modifier [29] (decreasing network connectivity) and a network former [28] (increasing network connectivity). Furthermore, when forming cation-anion clusters with fluoride, magnesium may exist in a four- or six-fold coordinated state [44].
Due to the multiple confounding structural variables in these systems, more structured and guided investigations for the identification of useful compositions is needed. While artificial intelligence models have been used for this purpose in silicate glasses, allowing from the predictions of glass properties across broad ranges of glass compositions, they have a lower modeling accuracy for more complex networks, such as multicomponent, fluoridated and borate glasses. As such, rapid screening across a broad range of compositions through design of mixtures experiments remains a relevant method of maximizing research efficacy, allowing for the selection of specific regions of interest for further investigation.
Most materials studied in this work were predominantly amorphous, with limited crystalline phases. Limits of inclusion to minimize crystallinity (defined in this case as below the 3% limit of accuracy of the methods used) were successfully identified using the polynomial equations in Table 4 and found to be: 82% for B2O3, 27% for MgO, 41% for CaO, 43% for Na2O. Upper limits for NaF and KF inclusion in this system were not identifiable (i.e., above 30% maximum investigated). The crystallinity models demonstrated a decreased tendency to crystalize with increased fluoride addition.
These findings counter previous investigations into the effect of fluoride addition in borate glasses. Reports of the effects of CaF2 addition on a SrO-B2O3 glass ceramics at either 2 or 5 mol% was also found to increase nucleation, improving the ease of formation of nanoceramics [45]. Glass forming ability also became increasingly difficult as fluoride was added when investigating the effect of fluoride addition to lithium borate glasses [46]. Such reports are consistent in the literature for bioactive glasses formed from phosphate and silicate networks.
Phase separation and crystallization was reported when CaF2 content was increased beyond 10 mol% in a phosphate glass [47]. Glass formation has also been reported to be hindered by CaF2 inclusion beyond 24% in phosphosilicate bioactive glasses [48,49], however it has been suggested that improper glass network design is responsible due to the increase in network connectivity observed when replacing CaO with CaF2 [49]. Some synthesis techniques allow for high rates of fluoride inclusion in borate glasses; notably, heavy metal borate glasses have been produced with up to 70% fluoride salts, however large cation radii of the heavy metal have been proposed to stabilize fluoride inclusion [30,31,32,33].
The impact of glass modifier oxides on glass density increased with increasing molar mass, suggesting that their impact is in part a result of atomic mass, as would be expected. This trend was not followed by NaF and KF, which imparted greater increases in density than would be expected from their atomic mass alone. Models developed for molar volume in this study found that B2O3 had the greatest positive effect on molar volume, that is to say, all glass modifiers investigated resulted in decreased molar volume relative to B2O3.
This finding is consistent with knowledge on borate glass structure, where glass modifier addition is known to increase boron coordination, increasing the packing efficiency. It is noteworthy that this trend is continuous in this glass series, with no reversal when studying borate contents from 95 to 45%, without evidence of a boron anomaly. While further structural investigation is required to better understand this finding, it appears that tetrahedral borate coordination continues to be the main function of modifier oxide addition through the glass series.
While cation radius has previously been used predict molar density in borate glass networks [50,51], no such trend was observed in this study. When compared to oxide glasses, molar volumes for fluoroborate glasses show less consistent trends. El-Egili et al. found an initial decrease of molar volume with increasing BaF2 a binary borate glass, but found this trend to reverse after a local minima at 40% BaF2, a finding attributed to the formation of BaF4 tetrahedra [31].
In further investigation of tertiary borate glasses, the substitution of PbF2 for CdF2 or BaF2 for PbF2 increased molar volume, a finding attributed to the larger size of the structural units formed [30,33]. The strong impact of NaF and KF on density may be explained by their impact on the molar volume. When looking at the relative effects of Na2O and NaF on molar volume, their coefficients in the predictive model are 0.31 and 0.18, respectively, indicating the NaF addition results in significantly smaller molar volume than Na2O.
The exact nature for this decrease in molar volume with NaF addition however is unclear. If fluoride ions segregate from the glass network into cation clusters, the addition of fluoride salts would be expected to decrease the fraction of tetrahedrally-coordinated boron, which would cause the glass network to expand. Alternatively, fluoride addition may favor the removal of non-bridging oxygen from the network, however the linear trend in molar volume and boron content suggests that tetrahedrally-coordinated boron’s likely dominate the structure.
Previous work suggests that molar volume may be predicted from the unit volumes of structural units formed within the glass, which can be garnered from binary borate glasses formed of each of their components [30,31,32,33,34,52]. This modeling however requires accurate structural knowledge of the glass, which is difficult to obtain for multicomponent borate glasses. Machine learning algorithms have been applied to predict glass density from large data sets, on both oxide glasses and fluoridated glasses [40,41,53]. While these studies have been relatively successful, the algorithms developed are skewed by the data sets on which they are trained.
They therefore can present accurate predictions for glasses that fall within the pallet of previously studied glasses but may (and do) fail to accurately model glasses with more novel compositions. Such an issue is highlighted in the works by Ahmmad et al. who used artificial intelligence techniques to predict the density of fluoride glasses from their chemical composition and ionic radii, but found that modeling of borate rich glasses was significantly less accurate [40].
A similar decrease in model accuracy was observed when attempting to predict the densities of multicomponent glasses when compared to binary, ternary, or quaternary materials [41]. Two factors may be contributing to the failure of the predictive models in both of these cases, (1) a relative lack of data for model fitting in borate and multicomponent glasses and (2) the added structural complexity of the glasses themselves, including the “boron anomaly” and mixed alkali/alkali earth effects.
Artificial intelligence algorithms are complex mathematic models that find patterns in the data on which they are “trained”, which when developed appropriately can be used to predict outcomes (in this case glass density) from inputs (in this case glass compositions). As such the failing of this approach is not surprising in the field of borates, which have been described as “enigmatic materials”, where discussion of structural changes center on the “borate anomaly” [24].
Similar difficulties can be expected with the prediction of magnesium containing glasses, where controversy lies in its bonding coordination and trends followed by other alkali earth oxides are broken [28,29,54]. Magnesium oxide, unlike oxides of Ca, Sr and Ba, is often identified as an intermediate glass modifier, capable of both tetrahedral coordination with oxygen as a glass network former and octahedral coordination with oxygen as a network modifier, further impacting glass network density [55].
As such, while the models are helpful in providing a prediction in density trends, their accuracy decreases for more unusual or novel glasses. Simply put, if artificial intelligence models are designed to predict expected results and can not be relied on to identify unusual behaviors. As such, the models presented in this work, offer a unique insight to navigate multicomponent borate glasses and allow for the identification of regions of interest for further investigation.
Glass transitions onset ranged from 346 to 607 °C. These temperatures are considerably lower than glass transition temperatures observed for binary borate glasses made from alkali and alkali earth metal oxides, which range from approximately 560 to 925 °C [28]. Alkali earth oxides had the greatest impact on glass transition, which could be explained by the formation of strongly bonded oxygen centers, increasing the network stiffness [26,27].
Similarly, low glass transition temperatures ranging from 382 to 660 °C were reported by Manchester et al. in the investigation of mixed alkali-alkali earth oxide effects in a K2O-SrO-B2O3 system, attributed to the lower fraction of B4 coordinated boron caused by the mixed alkali-alkali earth oxide effect [26]. This effect, thought to be of a similar nature to the mixed alkali effect, likely occurs due to the increased stability of adjacent B3 coordinated boron modified by the polarization effect caused by dissimilar metal ion modifiers [56].
In borate glasses, glass transition temperature is recognized to be closely linked to the fraction of B4 in the glass network, which increase rigidity. As such, glass transition temperature varies nonlinearly with modifier oxide addition similarly to the formation of B4 tetrahedra, reaching a maximum at approximately 30% oxide addition. Magnesium oxide is unique in that the glass-forming window is limited to 45–55% MgO, and the maximum B4 fraction is observed at ~50% MgO [42].
While rapid cooling techniques have been used to expand the glass forming region, they have resulted in significant decreases in glass transition (approximately 600 °C in the region of 50% and below 300 °C in the region of 0–30% MgO addition) [51]. The low glass transition temperatures observed in this study may be related to the low glass transition temperatures observed following low-level MgO inclusion in binary glasses, however the overall impact of magnesium addition in this work suggests it serves to increase the Tg.
This positive correlation (as visualised in Figure 1 and Figure 2) could be explained by the formation of covalent bonds as a network former, as well as by the stiffening effect of divalent cation bridges, as observed with calcium addition. While these finding suggest that magnesium in likely to act as a glass modifier in a multicomponent borate glass, such as the one studied (as its impact on glass transition is similar in magnitude and direction to that of calcium—a known network modifier), further structural analysis is required to fully understand the interaction between the simultaneous inclusion of magnesium among other glass components.
Fluoride addition offers another alternative explanation for the decreased glass transitions observed in this study. A decrease in Tg with increasing fluoride content was also observed in investigation of fluoridated lithium borate glasses [46]. Similar decreases in Tg with the addition of BaF2 in a binary glass has been attributed to the formation of terminal B-F bonds at the expense of B-O bonds [31]. Similar trends are found in silicate and phosphosilicate glasses [49,57], where it is attributed to the binding of fluorine ions to calcium in the glass network, preventing the calcium from acting as an ionic bridge between two adjacent non bridging oxygens and thus decreasing network rigidity [49].
A study on the impact of alkali earth fluorides (MgF2, CaF2, SrF2, BaF2) on the electrical properties of the MeF2-Na2B4O7 system similarly suggests that fluoride occupies some direct boron bonding sites, in the form of Na+[F-BO3/2] [44]. The MgF2 containing glasses were noted to retain more fluoride than their Ca, Sr or Ba counterparts during glass melting, due to the formation of octahedrally coordinated MgF6/2 centers [44]. Magnesium fluoride centers may decrease the impact of magnesium addition on the fraction of B4 however further investigation would be required to demonstrate if the additions of NaF or KF along with MgO has the same potential of forming MgF6/2.
Direct comparison of mass loss between the compounds studied and previously materials reported in the literature is difficult to achieve due to the lack of standardization in test methods, including the selection of elution media, particle size and mass (or surface) to volume ratios. Effort to standardize these variables using ISO 10993-14 has provided beneficial insight, however the methods are developed for more resilient glass ceramic systems and are not designed to study rapidly degrading materials.
With that said, some general trends may be compared to gain insight into the changes observed in these glasses. The glasses in this study demonstrated fast and extensive mass loss ranging from 41% to 100% after a single hour and rising to 60–100% following 24 h. This rapid dissolution could be beneficial for the use of glasses as ion delivery vehicles in oral care application, where daily exposure to the glasses could present issues related to the accumulation of unreacted glass particles [58].
While the high reactivity of borate glasses with low modifier inclusion is not unexpected, the ability of this glass system to create glasses of high calcium and fluoride inclusions while maintaining a high degree of reactivity is promising. While binary magnesium borate glasses show unusually high rates of dissolution, calcium oxide has been reported to have the greatest impact on supressing glass dissolution in alkali and alkali earth oxide borates [22]. This effect was similarly observed with the addition of calcium in this work, however extensive dissolution was maintained, providing a platform for controlled dissolution in a relatively high calcium content glass.
The strongest predictor of high rates of mass loss was associated with the B content in the glass, which is consistent with previous reports of boron being a soluble component of the glass [20]. Similarly, Na2O, NaF and KF served to increase mass loss when added to the glass network. Notably, KF increased mass loss at early time points, whereas B2O3 has a greater effect at extended time points. The finding of increased glass degradation with the addition of fluoride salts counters the findings in the study of phosphosilicate glasses, where fluoride addition (typically in the form of CaF2) decreases the rate of silica release into solution [59,60].
Two factors have been identified in the literature to explain these finding: changes to the glass network structure and changes to the pH of the extraction solution. Increases in network connectivity have been proposed to explain the decrease in reactivity in glass series where CaF2 is added into the glass at the expense of CaO or Na2O [61]. Similar findings are also observed in studies where the network connectivity is kept constant, suggesting this may be only partially responsible for the observations [59,60]. The decrease in reactivity for glasses in which the network connectivity is not altered may be explained by the ability of F– ions to exchange with OH– ions in solution, thereby decreasing the pH of the solution and causing the formation of a protective silica gel layer [59].
In contrast with phosphosilicate glasses, borate glass dissolution is not limited by the formation on a pacifying gel layer and results in a lesser increase in solution pH during dissolution. As such, borate glass networks can allow for a greater addition of fluoride to the glass network without compromising reactivity. This is a desirable feature of the design of a fluoride-releasing material, as full dissolution allows for greater release of fluoride into a bioavailable state (i.e., not retained in unreacted glass).
When investigating the weight loss in a fluoridated borate glass based off of full substitution of B2O3 for SiO2 in the original 45S5 composition with 2 wt% of fluoride salt weight loss was reported to be in the order of LiF > NaF > CaF2 > ZnF2, however no fluoride free composition was provided for comparison [45]. It is notable that due to the choice to use an equal weight percentage for the addition of fluoride salts, the amount of fluoride incorporated into the system would follow the same trend.
To best visualize the effects of fluoride addition on mass loss in their work, the 3D surface responses at both high and low calcium concentration can be observed as in Figure 6. At low calcium concentration (exemplified as 10% CaO), Na2O glasses demonstrate greater mass loss than their NaF counterparts. At high calcium concentration, (exemplified as 35% CaO), NaF glasses demonstrate greater mass loss than their Na2O counterparts.
This finding is counter intuitive to finding that glass dissolution is dependent on the solubilities of the dissolution products in solution, where the combination of Ca and F would be expected to decrease the solubility of the glass. In turn, structural effect, namely the decrease in the B:O ratio, and a decrease in the fraction of B4 coordinated fraction is likely to contribute to the change in dissolution. For a better understanding of the underlying mechanism through which fluoride addition alters glass dissolution structural analysis of the glasses is needed.

5. Conclusions

A design of mixtures approach allowed for the development of strong predictive models to assess the composition–property relationships in multicomponent fluoroborate glasses. Complex nonlinear relationships demonstrated the increased complexity in the borate glass systems relative to their silicate glass counterparts. High rates of fluoride inclusion, without negatively impacting glass formation or mass loss, demonstrated that mixed alkali–alkaline earth fluoroborate glasses offer an exceptional base for the design of fluoride-releasing glasses.

Author Contributions

Conceptualization, K.M. and D.B.; methodology, K.M. and D.B.; formal analysis, K.M.; writing—original draft preparation, K.M.; writing—review and editing, K.M. and D.B.; visualization, K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study is the property of IR Scientific and is available on request from the corresponding author.

Acknowledgments

The authors would like to thank Effie Boudreau for assistance in the synthesis and analysis of the materials studied.

Conflicts of Interest

The authors declare the following conflict of interest, Daniel Boyd and Kathleen MacDonald are shareholders and/or option holders in IR Scientific.

References

  1. Mouriño, V.; Cattalini, J.P.; Boccaccini, A.R. Metallic ions as therapeutic agents in tissue engineering scaffolds: An overview of their biological applications and strategies for new developments. J. R. Soc. Interface 2012, 9, 401–419. [Google Scholar] [CrossRef] [PubMed]
  2. Habibovic, P.; Barralet, J.E. Bioinorganics and biomaterials: Bone repair. Acta Biomater. 2011, 7, 3013–3026. [Google Scholar] [CrossRef]
  3. Fernandes, J.S.; Martins, M.; Neves, N.M.; Fernandes, M.H.V.; Reis, R.L.; Pires, R.A. Intrinsic antibacterial borosilicate glasses for bone tissue engineering applications. ACS Biomater. Sci. Eng. 2016, 2, 1143–1150. [Google Scholar] [CrossRef]
  4. Day, R.M. Bioactive glass stimulates the secretion of angiogenic growth factors and angiogenesis in vitro. Tissue Eng. 2005, 11, 768–777. [Google Scholar] [CrossRef]
  5. Fernandes, J.S.; Gentile, P.; Crawford, A.; Pires, R.A.; Hatton, P.V.; Reis, R.L. Substituted borosilicate glasses with improved osteogenic capacity for bone tissue engineering. Tissue Eng. Part A 2017, 23, 1331–1342. [Google Scholar] [CrossRef] [PubMed]
  6. Schmitz, S.; Widholz, B.; Essers, C.; Becker, M.; Tulyaganov, D.; Moghaddam, A.; de Juan, I.G.; Westhauser, F. Superior biocompatibility and comparable osteoinductive properties: Sodium-reduced fluoride-containing bioactive glass belonging to the CaO–MgO–SiO2 system as a promising alternative to 45S5 bioactive glass. Bioact. Mater. 2020, 5, 55–65. [Google Scholar] [CrossRef]
  7. Jones, J.R. Review of bioactive glass: From Hench to hybrids. Acta Biomater. 2013, 9, 4457–4486. [Google Scholar] [CrossRef]
  8. Wray, P. Cotton candy that heals? Borate glass nanofibers look promising. Am. Ceram. Soc. Bull. 2011, 90, 25–29. [Google Scholar]
  9. Castiglioni, S.; Cazzaniga, A.; Albisetti, W.; Maier, J.A.M. Magnesium and osteoporosis: Current state of knowledge and future research directions. Nutrients 2013, 5, 3022–3033. [Google Scholar] [CrossRef]
  10. Rondanelli, M.; Faliva, M.A.; Tartara, A.; Gasparri, C.; Perna, S.; Infantino, V.; Riva, A.; Petrangolini, G.; Peroni, G. An update on magnesium and bone health. Biometals 2021, 34, 715–736. [Google Scholar] [CrossRef]
  11. Wopenka, B.; Pasteris, J.D. A mineralogical perspective on the apatite in bone. Mater. Sci. Eng. C 2005, 25, 131–143. [Google Scholar] [CrossRef]
  12. Bigi, A.; Falini, G.; Foresti, E.; Gazzano, M.; Ripmonti, A.; Roveri, N. Rietveld structure refinements of calcium hydroxylapatite containing magnesium. Acta Crystallogr. Sect. B Struct. Sci. 1996, 52, 87–92. [Google Scholar] [CrossRef]
  13. Abdallah, M.-N.; Eimar, H.; Bassett, D.C.; Schnabel, M.; Ciobanu, O.; Nelea, V.; McKee, M.D.; Cerruti, M.; Tamimi, F. Diagenesis-inspired reaction of magnesium ions with surface enamel mineral modifies properties of human teeth. Acta Biomater. 2016, 37, 174–183. [Google Scholar] [CrossRef] [PubMed]
  14. Abdallah, M.N. Surface Reactivity of Tooth Enamel with Dyes Oxidizing Agents and Magnesium Ions and Its Effect on Tooth Color. Master’s Thesis, McGill University, Montreal, QC, Canada, 2012. [Google Scholar]
  15. Shaikh, K.; Pereira, R.; Gillam, D.G.; Phad, S. Comparative evaluation of desensitizing dentifrices containing BioMin®, Novamin® and fluoride on dentinal tubule occlusion before and after a citric acid challenge—A scanning electron microscope in-vitro study. J. Odontol. 2018, 2, 1000105. [Google Scholar]
  16. Reddy, G.V.; Surakanti, J.R.; Vemisetty, H.; Doranala, S.; Hanumanpally, J.R.; Malgikar, S. Comparative assessment of effectiveness of Biomin, NovaMin, herbal, and potassium nitrate desensitizing agents in the treatment of hypersensitive teeth: A clinical study. J. NTR Univ. Health Sci. 2019, 8, 24–28. [Google Scholar]
  17. Tirapelli, C.; Panzeri, H.; Soares, R.G.; Peitl, O.; Zanotto, E.D. A novel bioactive glass-ceramic for treating dentin hypersensitivity. Braz. Oral Res. 2010, 24, 381–387. [Google Scholar] [CrossRef]
  18. De Caluwe, T.; Vercruysse, C.W.J.; Declercq, H.A.; Schaubroeck, D.; Verbeeck, R.M.H.; Martens, L.C. Bioactivity and biocompatibility of two fluoride containing bioactive glasses for dental applications. Dent. Mater. 2016, 32, 1414–1428. [Google Scholar] [CrossRef]
  19. MacDonald, K.; Boudreau, E.; Thomas, G.V.; Badrock, T.C.; Davies, L.J.; Lloyd, M.J.; Spradbery, P.S.; Turner-Cahill, S.; Boyd, D. In vitro evaluation of Sensi-IP®: A soluble and mineralizing sensitivity solution. Heliyon 2021, 8, e08672. [Google Scholar] [CrossRef]
  20. Balasubramanian, P.; Büttner, T.; Pacheco, V.M.; Boccaccini, A.R. Boron-containing bioactive glasses in bone and soft tissue engineering. J. Eur. Ceram. Soc. 2018, 38, 855–869. [Google Scholar] [CrossRef]
  21. Liang, W.; Rüssel, C.; Day, D.E.; Völksch, G. Bioactive comparison of a borate, phosphate and silicate glass. J. Mater. Res. 2006, 21, 125–131. [Google Scholar] [CrossRef]
  22. Goetschius, K.L.; Beuerlein, M.A.; Bischoff, C.M.; Brow, R.K. Dissolution behavior of ternary alkali–alkaline earth-borate glasses in water. J. Non-Cryst. Solids 2018, 487, 12–18. [Google Scholar] [CrossRef]
  23. Wright, A.C.; Dalba, G. Borate versus silicate glasses: Why are they so different? Phys. Chem. Glasses Eur. J. Glass Sci. Technol. B 2010, 5, 233–265. [Google Scholar]
  24. Wright, A.C. My borate life: An enigmatic journey. Int. J. Appl. Glass Sci. 2015, 6, 45–63. [Google Scholar] [CrossRef]
  25. Hill, R. An alternative view of the degradation of bioglass. J. Mater. Sci. Lett. 1996, 15, 1122–1125. [Google Scholar] [CrossRef]
  26. Manchester, R.A.; Todorova, T.Z.; Werner-Zwanziger, U.; Boyd, D. Mixture designs to investigate the role of alkali and alkaline earth cations on composition–structure–property relationships in ternary borate glass networks. J. Non-Cryst. Solids 2021, 569, 120982. [Google Scholar] [CrossRef]
  27. MacDonald, K.; Hanson, M.A.; Boyd, D. Modulation of strontium release from a tertiary borate glass through substitution of alkali for alkali earth oxide. J. Non-Cryst. Solids 2016, 443, 184–191. [Google Scholar] [CrossRef]
  28. Yiannopoulos, Y.D.; Chryssikos, G.D.; Kamitsos, E.I. Structure and properties of alkaline earth borate glasses. Phys. Chem. Glasses 2001, 42, 164–172. [Google Scholar]
  29. Kamitsos, E.I.; Chryssikos, G.D.; Karakassides, M.A. Vibrational spectra of magnesium-sodium-borate glasses. 1. Far-infrared investigation of the cation-site interactions. J. Phys. Chem. 1987, 91, 1067–1073. [Google Scholar] [CrossRef]
  30. Doweidar, H.; El-Egili, K.; Ramadan, R.; Khalil, E. Structural species in mixed-fluoride PbF2–CdF2–B2O3 borate glasses; FTIR investigation. Vib. Spectrosc. 2019, 102, 24–30. [Google Scholar] [CrossRef]
  31. El-Egili, K.; Doweidar, H.; Ramadan, R.; Altawaf, A. Role of F− ions in the structure and properties of BaF2B2O3 glasses. J. Non-Cryst. Solids 2016, 449, 83–93. [Google Scholar] [CrossRef]
  32. Doweidar, H.; El-Egili, K.; Ramadan, R.; Khalil, E. Structural studies and properties of CdF2–B2O3 glasses. J. Non-Cryst. Solids 2018, 481, 494–502. [Google Scholar] [CrossRef]
  33. Doweidar, H.; El-Egili, K.; Altawaf, A. Structural units and properties of BaF2–PbF2–B2O3 glasses. J. Non-Cryst. Solids 2017, 464, 73–80. [Google Scholar] [CrossRef]
  34. Doweidar, H.; El-Damrawi, G.; Abdelghany, M. Structure and properties of CaF2–B2O3 glasses. J. Mater. Sci. 2012, 47, 4028–4035. [Google Scholar] [CrossRef]
  35. ElBatal, F.H.; Ouis, M.A.; ElBatal, H.A. Comparative studies on the bioactivity of some borate glasses and glass–ceramics from the two systems: Na2O–CaO–B2O3 and NaF–CaF2–B2O3. Ceram. Int. 2016, 42, 8247–8256. [Google Scholar] [CrossRef]
  36. Ouis, M.A.; Abdelghany, A.M.; Elbatal, H.A. Corrosion mechanism and bioactivity of borate glasses analogue to Hench’s bioglass. Process. Appl. Ceram. 2012, 6, 141–149. [Google Scholar] [CrossRef]
  37. Mauro, J.C.; Philip, C.S.; Vaughn, D.J.; Pambianchi, M.S. Glass science in the United States: Current status and future directions. Int. J. Appl. Glass Sci. 2014, 5, 2–15. [Google Scholar] [CrossRef]
  38. Han, T.; Stone-Weiss, N.; Huang, J.; Goel, A.; Kumar, A. Machine learning as a tool to design glasses with controlled dissolution for healthcare applications. Acta Biomater. 2020, 107, 286–298. [Google Scholar] [CrossRef]
  39. Krishnan, N.A.; Mangalathu, S.; Smedskjaer, M.M.; Tandia, A.; Burton, H.; Bauchy, M. Predicting the dissolution kinetics of silicate glasses using machine learning. J. Non-Cryst. Solids 2018, 487, 37–45. [Google Scholar] [CrossRef]
  40. Ahmmad, S.K.; Jabeen, N.; Ahmed, S.T.U.; Hussainy, S.F.; Ahmed, B. Density of fluoride glasses through artificial intelligence techniques. Ceram. Int. 2021, 47, 30172–30177. [Google Scholar] [CrossRef]
  41. Ahmmad, S.K.; Jabeen, N.; Ahmed, S.T.U.; Ahmed, S.A.; Rahman, S. Artificial intelligence density model for oxide glasses. Ceram. Int. 2021, 47, 7946–7956. [Google Scholar] [CrossRef]
  42. Wright, A.C. Borate structures: Crystalline and vitreous. Phys. Chem. Glasses Eur. J. Glass Sci. Technol. 2010, 51, 1–39. [Google Scholar]
  43. Aguiar, P.M.; Kroeker, S. Boron speciation and non-bridging oxygens in high-alkali borate glasses. J. Non-Cryst. Solids 2007, 353, 1834–1839. [Google Scholar] [CrossRef]
  44. Sokolov, I.A.; Naraev, V.N.; Nosakin, A.N.; Pronkin, A.A. Influence of MeF2 (Me = Mg, Ca, Sr, and Ba) on the electrical properties of glasses in the MeF2-Na2B4O7 system. Glass Phys. Chem. 2000, 26, 383–389. [Google Scholar] [CrossRef]
  45. Abdelghany, A.M.; Kamal, H. Spectroscopic investigation of synergetic bioactivity behavior of some ternary borate glasses containing fluoride anions. Ceram. Int. 2014, 40, 8003–8011. [Google Scholar] [CrossRef]
  46. Chowdari, B.V.R.; Rong, Z. Study of the fluorinated lithium borate glasses. Solid State Ion. 1995, 78, 133–142. [Google Scholar] [CrossRef]
  47. Shaharyar, Y.; Wein, E.; Kim, J.-J.; Youngman, R.E.; Muñoz, F.; Kim, H.-W.; Tilocca, A.; Goel, A. Structure-solubility relationships in fluoride-containing phosphate based bioactive glasses. J. Mater. Chem. B 2015, 3, 9360–9373. [Google Scholar] [CrossRef]
  48. Lusvardi, G.; Malavasi, G.; Cortada, M.; Menabue, L.; Menziani, M.C.; Pedone, A.; Segre, U. Elucidation of the structural role of fluorine in potentially bioactive glasses by experimental and computational investigation. J. Phys. Chem. B 2008, 112, 12730–12739. [Google Scholar] [CrossRef]
  49. Brauer, D.S.; Karpukhina, N.; Law, R.V.; Hill, R.G. Structure of fluoride-containing bioactive glasses. J. Mater. Chem. 2009, 19, 5629–5636. [Google Scholar] [CrossRef]
  50. Gaafar, M.S.; Marzouk, S.Y.; Zayed, H.A.; Soliman, L.I.; El-Deen, A.H.S. Structural studies and mechanical properties of some borate glasses doped with different alkali and cobalt oxides. Curr. Appl. Phys. 2013, 13, 152–158. [Google Scholar] [CrossRef]
  51. Lower, N.P.; McRae, J.L.; Feller, H.A.; Betzen, A.R.; Kapoor, S.; Affatigato, M.; Feller, S.A. Physical properties of alkaline-earth and alkali borate glasses prepared over an extended range of compositions. J. Non-Cryst. Solids 2001, 293, 669–675. [Google Scholar] [CrossRef]
  52. Doweidar, H.; El-Damrawi, G.; Mansour, E.; Fetouh, R.E. Structural role of MgO and PbO in MgO–PbO–B2O3 glasses as revealed by FTIR; A new approach. J. Non-Cryst. Solids 2012, 358, 941–946. [Google Scholar] [CrossRef]
  53. Deng, B. Machine learning on density and elastic property of oxide glasses driven by large dataset. J. Non-Cryst. Solids 2020, 529, 119768. [Google Scholar] [CrossRef]
  54. Kamitsos, E.I. Modifying role of alkali-metal cations in borate glass networks. J. Phys. Chem. 1989, 93, 1604–1611. [Google Scholar] [CrossRef]
  55. Varshneya, A.K. Fundamentals of Inorganic Glasses; Academic Press, Inc.: San Diego, CA, USA, 2013. [Google Scholar]
  56. Zhong, J.; Bray, P.J. Change in boron coordination in alkali borate glasses, and mixed alkali effects, as elucidated by NMR. JNCS 1988, 111, 67–76. [Google Scholar] [CrossRef]
  57. Bachar, A.; Mercier, C.; Tricoteaux, A.; Hampshire, S.; Leriche, A.; Follet, C. Effect of nitrogen and fluorine on mechanical properties and bioactivity in two series of bioactive glasses. J. Mech. Behav. Biomed. 2013, 23, 133–148. [Google Scholar] [CrossRef] [PubMed]
  58. Huang, W.; Day, D.E.; Kittiratanapiboon, K.; Rahaman, M.N. Kinetics and mechanisms of the conversion of silicate (45S5), borate, and borosilicate glasses to hydroxyapatite in dilute phosphate solutions. J. Mater. Sci. Mater. Med. 2006, 17, 583–596. [Google Scholar] [CrossRef]
  59. Brauer, D.S.; Karpukhina, N.; O’Donnell, M.D.; Law, R.V.; Hill, R.G. Fluoride-containing bioactive glasses: Effect of glass design and structure on degradation, pH and apatite formation in simulated body fluid. Acta Biomater. 2010, 6, 3275–3282. [Google Scholar] [CrossRef]
  60. Brauer, D.S.; Mneimne, M.; Hill, R.G. Fluoride-containing bioactive glasses: Fluoride loss during melting and ion release in tris buffer solution. J. Non-Cryst. Solids 2011, 357, 3328–3333. [Google Scholar] [CrossRef]
  61. Shah, F.A. Fluoride-containing bioactive glasses: Glass design, structure, bioactivity, cellular interactions, and recent developments. Mater. Sci. Eng. C 2016, 58, 1279–1289. [Google Scholar] [CrossRef]
Figure 1. Representative XRD spectrum collected from ground powder from quenched glass of (a) composition 9 and (b) composition 10.
Figure 1. Representative XRD spectrum collected from ground powder from quenched glass of (a) composition 9 and (b) composition 10.
Materials 15 06247 g001
Figure 2. 3D surface plots demonstrating the nonlinear relationships between the glass composition and (a) crystallinity, (b) density, (c) molar volume, (d) Tg onset, (e) Tf inflection and (f) Tg final.
Figure 2. 3D surface plots demonstrating the nonlinear relationships between the glass composition and (a) crystallinity, (b) density, (c) molar volume, (d) Tg onset, (e) Tf inflection and (f) Tg final.
Materials 15 06247 g002
Figure 3. Correlation table visualising the impact of glass components on the measured glass properties with correlations displayed from blue for negative correlations to red for positive correlations.
Figure 3. Correlation table visualising the impact of glass components on the measured glass properties with correlations displayed from blue for negative correlations to red for positive correlations.
Materials 15 06247 g003
Figure 4. The representative DSC trace collected from the powder of the quenched glass of composition 9.
Figure 4. The representative DSC trace collected from the powder of the quenched glass of composition 9.
Materials 15 06247 g004
Figure 5. 3D surface plots demonstrating the nonlinear relationships between the glass composition and (a) 1 h mass loss, (b) 1 h fluoride release, (c) 12 h mass loss, (d) 12 h fluoride release, (e) 24 h mass loss and (f) 24 h fluoride release.
Figure 5. 3D surface plots demonstrating the nonlinear relationships between the glass composition and (a) 1 h mass loss, (b) 1 h fluoride release, (c) 12 h mass loss, (d) 12 h fluoride release, (e) 24 h mass loss and (f) 24 h fluoride release.
Materials 15 06247 g005
Figure 6. Visualization of the effects of Na2O and NF inclusion on mass loss at 1 h at (a) low and (b) high CaO loadings.
Figure 6. Visualization of the effects of Na2O and NF inclusion on mass loss at 1 h at (a) low and (b) high CaO loadings.
Materials 15 06247 g006
Table 1. Design of the mixture components.
Table 1. Design of the mixture components.
ComponentMinimum (mol%)Maximum (mol%)Coded LowCoded HighMean StdDev.
B2O34595+0 ↔ 45+1 ↔ 9554.1413.12
MgO150+0 ↔ 1+0.98 ↔ 509.3312.33
CaCO3150+0 ↔ 1+0.98 ↔ 509.7712.66
Na2CO3150+0 ↔ 1+0.98 ↔ 5011.1812.71
NaF130+0 ↔ 1+0.58 ↔ 307.958.80
KF130+0 ↔ 1+0.58 ↔ 307.638.64
Total = 100%L Pseudo Coding
Table 2. Compositions of the 31 investigational glasses, denoted in molar % as identified by the design of mixtures outlined in Table 1 (repeated compositions demoted by matching Greek letters).
Table 2. Compositions of the 31 investigational glasses, denoted in molar % as identified by the design of mixtures outlined in Table 1 (repeated compositions demoted by matching Greek letters).
B2O3MgOCaCO3Na2CO3NaFKF
146.00001.000001.0000050.000001.000001.00000
267.885201.321781.8248023.484401.000004.48374
354.2486011.601909.071792.3575821.720101.00000
445.000005.304831.0000017.6952030.000001.00000
5 α45.0000023.094201.0000022.938103.047754.91997
6 β68.646701.0000023.119501.000002.954583.27929
772.525601.000001.000006.8845217.589901.00000
848.200904.0099113.572809.3802511.0428013.79340
9 α45.0000023.094201.0000022.938103.047754.91997
1045.000003.993035.033661.0212514.9521030.00000
1145.000001.000001.0000022.000001.0000030.00000
1252.190401.000001.000001.0000030.000014.80960
13 β68.646701.0000023.119501.000002.954583.27929
14 γ45.8861023.5111024.032802.007223.562741.00000
1545.0000027.527001.000001.0000018.738306.73472
16 δ45.000001.0000023.6165023.344804.330262.70839
1795.000001.000001.000001.000001.000001.00000
1846.249201.0000026.880101.0000020.223204.64742
19 ε45.000001.000001.0000025.2943014.0757013.63000
2045.0000050.000001.000002.000001.000001.00000
2145.000004.0348329.276702.659461.0000018.02900
2272.343901.000001.000001.000004.8365719.81960
2354.174901.000003.3621833.290407.172491.00000
24 ε45.000001.000001.0000025.2943014.0757013.63000
2546.000001.0000050.000001.000001.000001.00000
2648.4771028.602702.032781.000001.0000018.88740
2759.1174013.302001.0000011.269801.0000014.31080
2868.3276023.232401.185482.064534.189961.00000
2978.388507.961525.197006.452991.000001.00000
30 γ45.8861023.5111024.032802.007223.562741.00000
31 δ45.000001.0000023.6165023.344804.330262.70839
Table 3. The material characterization results.
Table 3. The material characterization results.
% CrystallinityDensity
(g/cm3)
Molar Volume
(cm3)
Tg Onset
(°C)
Tg Inflection
(°C)
Tg Final
(°C)
1 *88.0-----
29.02.2629.45472.4489.4479.9
32.02.4024.49505.7521.5513.0
43.02.4523.77392.4404.7403.1
57.92.4224.66422.9434.4431.6
6 *28.0--589.2604.9599.4
7 *25.32.1829.26433.1453.1444.9
8 *17.92.5024.61429.1444.1435.2
93.52.4424.41420.4433.0429.1
10 *19.02.4125.02381.3400.6394.4
11 *13.12.3627.08345.3357.1354.7
122.32.3025.77384.2399.4395.6
13 *19.9-----
142.02.5722.68566.9583.5575.6
155.02.4522.65480.9494.9490.4
164.92.4625.56430.2441.0440.3
17 *42.91.8836.77---
185.02.5023.90555.1569.8566.5
191.92.4025.85358.9369.9367.1
20------
211.92.5124.75489.8503.2501.9
22 *18.42.1230.95387.5400.2400.6
233.62.4126.77419.4433.1427.0
241.42.4225.59357.7369.9367.8
25 *38.12.7023.03607.8618.4617.4
264.22.3425.02509.1523.8521.9
271.72.3227.11457.6471.7468.1
282.42.2727.01538.8558.8556.4
292.52.1730.32476.0501.5492.1
302.32.5722.70568.6582.4576.7
317.12.4625.64582.4439.8435.2
* Compositions found to be cloudy or coloured upon quenching. Note: Missing values are not reported for compositions that yielded insufficient material for analysis.
Table 4. Regression model equations (actual equations) generated through design expert for the prediction of glass properties, along with the model fit statistics.
Table 4. Regression model equations (actual equations) generated through design expert for the prediction of glass properties, along with the model fit statistics.
ResponseRegression ModelR2R2 AdjustedR2 PredictedAdequate Precision
Crystallinity0.40 × B2O3 + 6.42 × MgO + 0.19 × CaCO3 − 0.45 × Na2CO3 − 0.27 × NaF + 0.12 × KF − 0.10 × B2O3 × MgO − 0.09 × MgO × CaCO3 − 0.05 × MgO × Na2CO3 − 0.06 × MgO × NaF − 0.10 × MgO × KF0.880.860.4420.10
Density0.018 × B2O3 + 0.023 × MgO + 0.035 × CaCO3 + 0.021 × Na2CO3 + 0.029 × NaF + 0.027 × KF + 0.00015 × B2O3 × MgO + 0.00018 × B2O3 × Na2CO3 − 0.00017 × CaCO3 × Na2CO30.980.970.9545.83
Molar Volume0.39 × B2O3 + 0.23 × MgO + 0.34 × CaCO3 + 0.31 × Na2CO3 + 0.18 × NaF + 0.26 × KF − 0.0033 × B2O3 × MgO + -0.0054 × B2O3 × CaCO3 − 0.0031 × B2O3 × Na2CO3 − 0.0019 × B2O3 × NaF − 0.0016 × B2O3 × KF + 0.0015 × CaCO3 × Na2CO30.990.990.9766.66
Tg Onset5.11 × B2O3 + 6.12 × MgO + 7.69 × CaCO3 + 2.41 × Na2CO3 + 2.98 × NaF + 1.25 × KF0.840.800.7417.83
Tg Inflection4.17 × B2O3 + 3.06 × MgO + 0.99 × CaCO3 − 2.96 × Na2CO3 + 3.87 × NaF + 2.38 × KF + 0.092 × B2O3 × MgO + 0.16 × B2O3 × CaCO3 + 0.14 × B2O3 × Na2CO3 − 0.067 × MgO × Na2CO3 + 0.066 × MgO × KF − 0.10 × CaCO3 × Na2CO3 + 0.10× CaCO3 × NaF0.990.980.8335.24
Tg Final4.35 × B2O3 + 9.32 × MgO + 1.83 × CaCO3 − 1.89 × Na2CO3 + 3.75 × NaF + 2.11 × KF + 0.14 × B2O3 × CaCO3 + 0.11 × B2O3 × Na2CO3 − 0.087 × MgO × CaCO3 − 0.14 × MgO × Na2CO3 − 0.090 × MgO × NaF − 0.11 × CaCO3 × Na2CO3 + 0.10 × CaCO3 × NaF0.990.980.8833.42
1 h Mass Loss(100 × ey)/(1 + ey)
where y = 0.10 × B2O3 + 0.17 × MgO + 0.88 × CaCO3 + 0.010 × Na2CO3 − 0.061 × NaF − 0.036 × KF − 0.0059 × B2O3 × MgO − 0.021 × B2O3 × CaCO3 − 0.0032 × CaCO3 × Na2CO3
0.940.920.8618.12
12 h Mass Loss(100 × ey)/(1 + ey)
where y = 0.045 × B2O3 − 0.026 × MgO − 0.055 × CaCO3 + 0.046 × Na2CO3 + 0.0024 × NaF − 0.00077 × KF
0.760.700.5912.32
24 h Mass Loss(100 × ey)/(1 + ey)
where y = 0.045 × B2O3 − 0.027 × MgO + −0.054 × CaCO3 + 0.043 × Na2CO3 + 0.0022 × NaF + 0.0010 × KF
0.740.670.5511.51
1 h Fluoride Release(2000 × ey)/(1 + ey)
where y = −0.057 × B2O3 − 0.015 × MgO − 0.045 × CaO+ 0.18 × Na2O
+ 0.031 × KF + 0.14 × NaF − 0.0045 × CaO × NaF + 0.0039 × Na2O × KF.
0.850.790.7011.68
12 h Fluoride Release(2000 × ey)/(1 + ey)
where y = −0.053 × B2O3 − 0.027 × MgO − 0.040 × CaCO3 + 0.080 × Na2CO3 -0.20 × NaF + 0.036 × KF − 0.0019 × B2O3 × Na2CO3 + 0.0050× B2O3 × NaF + 0.0031 × MgO × Na2CO3 + 0.0029 × Na2CO3 × NaF − 0.0026 × NaF × KF
0.950.900.7618.31
24 h Fluoride Release(2000 × ey)/(1 + ey)
where y = −0.059 × B2O3 − 0.0029 × MgO − 0.0384 × CaCO3 + 0.014 × Na2CO3 − 0.14 × NaF + 0.080 × KF + 0.0035 × B2O3 × NaF − 0.0016 × MgO × KF − 0.0020 × CaCO3 × KF + 0.0021 × Na2CO3 × NaF
0.980.960.8429.5
Table 5. The mass loss and fluoride release values following incubation of glass powders in TRIS-buffered saline, compositions without reported values yielded insufficient materials for processing and analysis.
Table 5. The mass loss and fluoride release values following incubation of glass powders in TRIS-buffered saline, compositions without reported values yielded insufficient materials for processing and analysis.
Mass Loss (%)Fluoride Release (ppm)
1 h12 h24 h1 h12 h24 h
1------
299.898.498.4168.425.274.6
353.579.579.5360.6318.6152.5
493.791.991.91918.1946.9582.5
582.879.179.1714.6459.1191.0
6------
798.693.293.2670.4463.5201.0
861.973.073.0500448.0165.0
9------
1081.682.382.31342.9177.7483.0
1197.892.492.41725.7394.71064.5
1296.095.295.21216.8273.9766.5
13------
1447.961.161.139.646.749.2
1555.181.981.9674.898.2202.0
1672.573.562.7112111.594.7
17------
1856.960.864.771.981.179.5
1996.899.299.21915.9418.9854.5
20------
2165.559.764.382.966.071.5
2299.795.995.9119.4113.7223.0
2398.099.499.4149.2144.7329.0
2492.696.096.0441.2414.5808.0
2540.555.059.316.921.321.2
2667.885.985.9127.5115.2187.0
2792.191.993.0137.9134.5117.9
2864.493.893.850.564.564.7
2976.397.297.214.416.616.7
3048.362.862.341.650.752.1
3166.967.869.6110.4103.495.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

MacDonald, K.; Boyd, D. Investigation of Multicomponent Fluoridated Borate Glasses through a Design of Mixtures Approach. Materials 2022, 15, 6247. https://doi.org/10.3390/ma15186247

AMA Style

MacDonald K, Boyd D. Investigation of Multicomponent Fluoridated Borate Glasses through a Design of Mixtures Approach. Materials. 2022; 15(18):6247. https://doi.org/10.3390/ma15186247

Chicago/Turabian Style

MacDonald, Kathleen, and Daniel Boyd. 2022. "Investigation of Multicomponent Fluoridated Borate Glasses through a Design of Mixtures Approach" Materials 15, no. 18: 6247. https://doi.org/10.3390/ma15186247

APA Style

MacDonald, K., & Boyd, D. (2022). Investigation of Multicomponent Fluoridated Borate Glasses through a Design of Mixtures Approach. Materials, 15(18), 6247. https://doi.org/10.3390/ma15186247

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop