An experimental and theoretical approach on stability towards hydrolysis of triethyl phosphate and its effects on the microstructure of sol-gel-derived bioactive silicate glass

The sol-gel method is versatile and one of the well-established synthetic approaches for preparing bioactive glass with improved microstructure. In a successful approach, alkoxide precursors undergo rapid hydrolysis, followed by immediate condensation leading to the formation of three-dimensional gels. On the other hand, a slow kinetics rate for hydrolysis of one or more alkoxide precursors generates a mismatch in the progression of the consecutive reactions of the sol-gel process, which makes it difficult to form homogeneous multicomponent glass products. The amorphous phase separation (APS) into the gel is thermodynamically unstable and tends to transform into a crystalline form during the calcination step of xerogel. In the present study, we report a combined experimental and theoretical method to investigate the stability towards hydrolysis of triethyl phosphate (TEP) and its effects on the mechanism leading to phase separation in 58S bioactive glass obtained via sol-gel route. A multitechnical approach for the experimental characterization combined with calculations of functional density theory (DFT) suggest that TEP should not undergo hydrolysis by water under acidic conditions during the formation of the sol or even in the gel phase. The activation energy barrier (ΔG^‡) showed a height of about 20 kcal·mol^-1 for the three stages of hydrolysis and the reaction rates calculated for each stage of TEP hydrolysis were k_FHR = 7.0 × 10^-3s^-1, k_SHR = 6.8 × 10^-3s^-1 and k_THR = 3.5 × 10^-3s^-1. These results show that TEP remains in the non-hydrolyzed form segregated within the xerogel matrix until its thermal decomposition in the calcination step, when P species preferentially associate with calcium ions (labile species) and other phosphate groups present nearby, forming crystalline domains of calcium pyrophosphates permeated by the silica-rich glass matrix. Together, our data expand the knowledge about the synthesis by the sol-gel method of bioactive glass and establishes a mechanism that explains the role played by the precursor source of phosphorus (TEP) in the phase separation, an event commonly observed for these biomaterials.

Understanding the sodium ion transport properties, deintercalation mechanism, and phase evolution of a Na2Mn2Si2O7 cathode by atomistic simulation

Molecular dynamics (MD) together with the first principles method (DFT) reveal that Na+ is capable of migrating three dimensionally in a Na2Mn2Si2O7 cathode material. Migration along the a-axis and c-axis have the same mechanism, that is, alternating between the Na1 and Na2 route with a similar local environment and distance. Long-distance hopping between two Na2 atoms or between Na1 and Na2 atoms is crucial for continuous migration along the b-axis. Also, the anti-site phenomenon is identified, and it facilitates the migration of the Na ions. Four intermediate phases are determined according to the formation energy curve and, as a result, the voltage profile is predicted accurately. The state of charge (SOC) dependency of the Na+ energy shows that the mobility of Na+ is highly inhibited in the fully discharged state. Upon the deintercalation of sodium ions, Na+ is activated immediately. A maximal DNa+ value of 3.6 × 10-9 cm2 s-1 and a low energy barrier of ca. 0.26 eV at the deintercalation level of x = 0.25 are observed. Because of the scarcity of Na+, DNa+ experiences a sharp decrease at the end of deintercalation. Despite the low level of Na+ mobility in the range of 0.25 < x < 1, Na2Mn2Si2O7 is still a potential cathode material for use in sodium ion batteries (SIBs).

B2O3/SiO2 substitution effect on structure and properties of Na2O–CaO–SrO–P2O5–SiO2 bioactive glasses from molecular dynamics simulations

The effect of B₂O₃/SiO₂ substitution in SrO-containing 55S4.3 bioactive glasses on glass structure and properties, such as ionic diffusion and glass transition temperature, was investigated by combining experiments and molecular dynamics simulations with newly developed potentials.

Structure-solubility relationships in fluoride-containing phosphate based bioactive glasses

Structural role of fluoride on chemical dissolution behavior of bioactive phosphate glasses has been studied.

Cooling rate and size effects on the medium-range structure of multicomponent oxide glasses simulated by molecular dynamics

A set of molecular dynamics simulations were performed to investigate the effect of cooling rate and system size on the medium-range structure of melt-derived multicomponent silicate glasses, represented by the quaternary 45S5 Bioglass composition. Given the significant impact of the glass degradation on applications of these materials in biomedicine and nuclear waste disposal, bulk structural features which directly affect the glass dissolution process are of particular interest. Connectivity of the silicate matrix, ion clustering and nanosegregation, distribution of ring and chain structural patterns represent critical features in this context, which can be directly extracted from the models. A key issue is represented by the effect of the computational approach on the corresponding glass models, especially in light of recent indications questioning the suitability of conventional MD approaches (that is, involving melt-and-quench of systems containing ~10(3) atoms at cooling rates of 5-10 K/ps) when applied to model these glasses. The analysis presented here compares MD models obtained with conventional and nonconventional cooling rates and system sizes, highlighting the trend and range of convergence of specific structural features in the medium range. The present results show that time-consuming computational approaches involving much lower cooling rates and/or significantly larger system sizes are in most cases not necessary in order to obtain a reliable description of the medium-range structure of multicomponent glasses. We identify the convergence range for specific properties and use them to discuss models of several glass compositions for which a possible influence of cooling-rate or size effects had been previously hypothesized. The trends highlighted here represent an important reference to obtain reliable models of multicomponent glasses and extract converged medium-range structural features which affect the glass degradation and thus their application in different fields. In addition, as a first application of the present findings, the fully converged structure of the 45S5 glass was further analyzed to shed new light on several dissolution-related features whose interpretation has been rather controversial in the past.

On the structure of biomedical silver-doped phosphate-based glasses from molecular dynamics simulations

First-principles and classical molecular dynamics simulations have been carried out on undoped and silver-doped phosphate-based glasses with 50 mol% P₂O₅, 0–20 mol% Ag₂O, and varying amounts of Na₂O and CaO.

Bioactive glasses as potential radioisotope vectors for in situ cancer therapy: investigating the structural effects of yttrium

The incorporation of yttrium in bioactive glasses (BGs) could lead to a new generation of radionuclide vectors for cancer therapy, with high biocompatibility, controlled biodegradability and the ability to enhance the growth of new healthy tissues after the treatment with radionuclides. It is essential to assess whether and to what extent yttrium incorporation affects the favourable properties of the BG matrix: ideally, one would like to combine the high surface reactivity typical of BGs with a slow release of radioactive yttrium. Molecular Dynamics simulations show that, compared to a BG composition with the same silica fraction, incorporation of yttrium results in two opposing effects on the glass durability: a more fragmented silicate network (leading to lower durability) and a stronger yttrium-mediated association between separate silicate fragments (leading to higher durability). The simulations also highlight a high site-selectivity and some clustering of yttrium cations, which are likely linked to the observed slow rate of yttrium released from related Y-BG compositions. Optimisation of yttrium BG compositions for radiotherapy applications thus depends on the delicate balance between these effects.

Surface signatures of bioactivity: MD simulations of 45S and 65S silicate glasses

The surface of a bioactive (45S) and a bioinactive (65S) glass composition has been modeled using shell-model classical molecular dynamics simulations. Direct comparison of the two structures allowed us to identify the potential role of specific surface features in the processes leading to integration of a bioglass implant with the host tissues, focusing in particular on the initial dissolution of the glass network. The simulations highlight the critical role of network fragmentation and sodium enrichment of the surface in determining the rapid hydrolysis and release of silica fragments in solution, characteristic of highly bioactive compositions. On the other hand, no correlation has been found between the surface density of small (two- and three-membered) rings and bioactivity, thus suggesting that additional factors need to be taken into account to fully understand the role of these sites in the mechanism leading to calcium phosphate deposition on the glass surface.

Modeling the water-bioglass interface by ab initio molecular dynamics simulations

The hydration of the surface of a highly bioactive silicate glass was modeled using ab initio (Car-Parrinello) molecular dynamics (CPMD) simulations, focusing on the structural and chemical modifications taking place at the glass-water interface immediately after contact and on the way in which they can affect the bioactivity of these materials. The adsorption of a water dimer and trimer on the dry surface was studied first, followed by the extended interface between the glass and liquid water. The CPMD trajectories provide atomistic insight into the initial stages relevant to the biological activity of these materials: following contact of the glass with an aqueous (physiological) medium, the initial enrichment of the surface region in Na+ cations establishes dominant Na+-water interactions at the surface, which allow water molecules to penetrate into the open glass network and start its partial dissolution. The model of a Na/H-exchanged interface shows that Ca2+-water interactions are mainly established after the dominant fraction of Na is leached into the solution. Another critical role of modifier cations was highlighted: they provide the Lewis acidity necessary to neutralize OH(-) produced by water dissociation and protonation of nonbridging oxygen (NBO) surface sites. The CPMD simulations also highlighted an alternative, proton-hopping mechanism by which the same process can take place in the liquid water film. The main features of the bioactive glass surface immediately after contact with an aqueous medium, as emerged from the simulations, are (a) silanol groups formed by either water dissociation at undercoordinated Si sites or direct protonation of NBOs, (b) OH(-) groups generally stabilized by modifier cations and coupled with the protonated NBOs, and (c) small rings, relatively stable and unopened even after exposure to liquid water. The possible role and effect of these sites in the bioactive process are discussed.

Structural models of bioactive glasses from molecular dynamics simulations

The bioactive mechanism, by which living tissues attach to and integrate with an artificial implant through stable chemical bonds, is at the core of many current medical applications of biomaterials, as well as of novel promising applications in tissue engineering. Having been employed in these applications for almost 40 years, soda-lime phosphosilicate glasses such as 45S5 represent today the paradigm of bioactive materials. Despite their strategical importance in the field, the relationship between the structure and the activity of a glass composition in a biological environment has not been studied in detail. This fundamental gap negatively affects further progress, for instance, to improve the chemical durability and tailor the biodegradability of these materials for specific applications. This paper reviews recent advances in computer modelling of bioactive glasses based on molecular dynamics simulations, which are starting to unveil key structural features of these materials, thus contributing to improve our fundamental understanding of how bioactive materials work.