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

Structural Causes of Brittleness Changes in Aluminosilicate Glasses with Different Cooling Rates

Numerous sources have already demonstrated that varying annealing rates can result in distinct toughness and brittleness in glass. To determine the underlying mechanisms driving this phenomenon, molecular dynamic (MD) simulations were employed to investigate the microstructure of aluminosilicate glasses under different cooling rates, and then uniaxial stretching was performed on them under controlled conditions. Results indicated that compared with short-range structure, cooling rate has a greater influence on the medium-range structure in glass, and it remarkably affects the volume of voids. Both factors play a crucial role in determining the brittleness of the glass. The former adjusts network connectivity to influence force transmission by manipulating the levels of bridging oxygen (BO) and non-bridging oxygen (NBO), and the latter accomplishes the objective of influencing brittleness by modifying the environmental conditions that affect the changes in BO and NBO content. The variation in the void environment results in differences in the strategies of the changes in BO and NBO content during glass stress. These findings stem from the excellent response of BO and NBO to the characteristic points of stress-strain curves during stretching. This paper holds importance in understanding the reasons behind the effect of cooling rates on glass brittleness and in enhancing our understanding of the ductile/brittle transition (DTB) in glass.

Review: understanding the properties of amorphous materials with high-performance computing methods

Amorphous materials have no long-range order in their atomic structure. This makes much of the formalism for the study of crystalline materials irrelevant, and so elucidating their structure and properties is challenging. The use of computational methods is a powerful complement to experimental studies, and in this paper we review the use of high-performance computing methods in the simulation of amorphous materials. Five case studies are presented to showcase the wide range of materials and computational methods available to practitioners in this field. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Multi-technique Approach to Unravel the (Dis)order in Amorphous Materials

The concept of order in disordered materials is the key to controlling the mechanical, electrical, and chemical properties of amorphous compounds widely exploited in industrial applications and daily life. Rather, it is far from being understood. Here, we propose a multi-technique numerical approach to study the order/disorder of amorphous materials on both the short- and the medium-range scale. We combine the analysis of the disorder level based on chemical and physical features with their geometrical and topological properties, defining a previously unexplored interplay between the different techniques and the different order scales. We applied this scheme to amorphous GeSe and GeSeTe chalcogenides, showing a modulation of the internal disorder as a function of the stoichiometry and composition: Se-rich systems are less ordered than Ge-rich systems at the short- and medium-range length scales. The present approach can be easily applied to more complex systems containing three or more atom types without any a priori knowledge about the system chemical-physical features, giving a deep insight into the understanding of complex systems.

Biogenic Synthesis and Characterization of Antioxidant and Antimicrobial Silver Nanoparticles Using Flower Extract of Couroupita guianensis Aubl

Couroupita guianensis Aubl. is an important medicinal tree. This tree is rich in various phytochemicals, and is therefore used as a potent antioxidant and antibacterial agent. This plant is also used for the treatment of various diseases. Here, we have improved its medicinal usage with the biosynthesis of silver nanoparticles (AgNPs) using Couroupita guianensis Aubl. flower extract as a reducing and capping agent. The biosynthesis of the AgNPs reaction was carried out using 1 mM of silver nitrate and flower extract. The effect of the temperature on the biosynthesis of AgNPs was premeditated by room temperature (25 °C) and 60 °C. The continuous stirring of the reaction mixture at room temperature for approximately one hour resulted in the successful formation of AgNPs. A development of a yellowish brown color confirmed the formation of AgNPs. The efficacious development of AgNPs was confirmed by the characteristic peaks of UV-Vis, X-ray diffraction (XRD) and Fourier transform infrared (FT-IR) spectroscopy spectra. The biosynthesized AgNPs exhibited significant free radical scavenging activity through a DPPH antioxidant assay. These AgNPs also showed potent antibacterial activity against many pathogenic bacterial species. The results of molecular dynamics simulations also proved the average size of NPs and antibacterial potential of the flower extract. The observations clearly recommended that the green biosynthesized AgNPs can serve as effective antioxidants and antibacterial agents over the plant extract.

Statistical topology of bond networks with applications to silica

Whereas knowledge of a crystalline material's unit cell is fundamental to understanding the material's properties and behavior, there are no obvious analogs to unit cells for disordered materials despite the frequent existence of considerable medium-range order. This article views a material's structure as a collection of local atomic environments that are sampled from some underlying probability distribution of such environments, with the advantage of offering a unified description of both ordered and disordered materials. Crystalline materials can then be regarded as special cases where the underlying probability distribution is highly concentrated around the traditional unit cell. The H_{1} barcode is proposed as a descriptor of local atomic environments suitable for disordered bond networks and is applied with three other descriptors to molecular dynamics simulations of silica glasses. Each descriptor reliably distinguishes the structure of glasses produced at different cooling rates, with the H_{1} barcode and coordination profile providing the best separation. The approach is generally applicable to any system that can be represented as a sparse graph.

Structure Prediction of Rare Earth Doped BaO and MgO Containing Aluminosilicate Glasses⁻the Model Case of Gd₂O₃

The medium-range atomic structure of magnesium and barium aluminosilicate glasses doped with Gd₂O₃ as a model rare earth oxide is elucidated using molecular dynamics simulations. Our structure models rationalize the strong dependence of the luminescence properties of the glasses on their chemical composition. The simulation procedure used samples’ atomic configurations, the so-called inherent structures, characterizing configurations of the liquid state slightly above the glass transition temperature. This yields medium-range atomic structures of network former and modifier ions in good agreement with structure predictions using standard simulated annealing procedures. However, the generation of a large set of inherent structures allows a statistical sampling of the medium-range order of Gd³⁺ ions with less computational effort compared to the simulated annealing approach. It is found that the number of Si-bound non-bridging oxygen in the vicinity of Gd³⁺ considerably increases with growing ionic radius and concentration of network-modifier ions. In addition, structure predictions indicate a low driving force for clustering of Gd³⁺, yet no precise correlation between the atomic structure and luminescence lifetimes can be conclusively established. However, the structure models provided in this study can serve as a starting point for future quantum mechanical simulations to shed a light on the relation between the atomic structure and optical properties of rare earth doped aluminosilicate glasses.

Structure–composition trends in multicomponent borosilicate-based glasses deduced from molecular dynamics simulations with improved B–O and P–O force fields

Improved B–O and P–O force fields provide accurate molecular dynamics simulations of multicomponent glasses.

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.

Structures and properties of phosphate-based bioactive glasses from computer simulation: a review

Computer simulations have enabled breakthroughs in understanding the connections between the atomic structure and properties of bioactive phosphate glasses.

Influence of Calcium on the Initial Stages of the Sol-Gel Synthesis of Bioactive Glasses

Understanding how calcium interacts with silica sources and influences their polycondensation in aqueous solutions is of central importance for the development of more effective biomaterials by sol-gel approaches. For this purpose, the atomic-scale evolutions of a calcium-containing precursor solution corresponding to a typical sol-gel bioactive glass and of a corresponding Ca-free solution were compared using reactive molecular dynamics simulations. The simulations highlight a significantly faster rate of condensation when calcium is present in the initial solution, resulting in the formation of large and ramified silica clusters within 5 ns, which are absent in the Ca-free system. This different behavior has been analyzed and interpreted in terms of the Ca-induced nanosegregation in calcium-rich and silica-rich regions, which promotes the condensation reactions within the latter. By identifying a possible mechanism behind the limited incorporation of calcium in the silica nanoclusters formed in the early stages of the sol-gel process, these results could guide further studies aimed at identifying favorable experimental conditions to enhance initial calcium incorporation and thus produce sol-gel biomaterials with improved properties.

Simulations reveal the role of composition into the atomic-level flexibility of bioactive glass cements

Bioactive glass ionomer cements (GICs), the reaction product of a fluoro-alumino-silicate glass and polyacrylic acid, have been in effective use in dentistry for over 40 years and more recently in orthopaedics and medical implantation. Their desirable properties have affirmed GIC's place in the medical materials community, yet are limited to non-load bearing applications due to the brittle nature of the hardened composite cement, thought to arise from the glass component and the interfaces it forms. Towards helping resolve the fundamental bases of the mechanical shortcomings of GICs, we report the 1st ever computational models of a GIC-relevant component. Ab initio molecular dynamics simulations were employed to generate and characterise three fluoro-alumino-silicate glasses of differing compositions with focus on resolving the atomic scale structural and dynamic contributions of aluminium, phosphorous and fluorine. Analyses of the glasses revealed rising F-content leading to the expansion of the glass network, compression of Al-F bonding, angular constraint at Al-pivots, localisation of alumino-phosphates and increased fluorine diffusion. Together, these changes to the structure, speciation and dynamics with raised fluorine content impart an overall rigidifying effect on the glass network, and suggest a predisposition to atomic-level inflexibility, which could manifest in the ionomer cements they form.

Cooling rate and stress relaxation in silica melts and glasses via microsecond molecular dynamics

We have conducted extremely long molecular dynamics simulations of glasses to microsecond times, which close the gap between experimental and atomistic simulation time scales by two to three orders of magnitude. Static, thermal, and structural properties of silica glass are reported for glass cooling rates down to 5×10(9) K/s and viscoelastic response in silica melts and glasses are studied over nine decades of time. We present results from relaxation of hydrostatic compressive stress in silica and show that time-temperature superposition holds in these systems for temperatures from 3500 to 1000 K.

Na/Ca Intermixing around Silicate and Phosphate Groups in Bioactive Phosphosilicate Glasses Revealed by Heteronuclear Solid-State NMR and Molecular Dynamics Simulations

We characterize the intermixing of network-modifying Na(+)/Ca(2+) ions around the silicate (QSi(n)) and phosphate (QP(n)) tetrahedra in a series of 16 Na2O–CaO–SiO2–P2O5 glasses, whose P content and silicate network connectivity were varied independently. The set includes both bioactive and bioinactive compositions and also encompasses two soda-lime-silicate members devoid of P, as well as two CaO–SiO2 glasses and one Na2O–SiO2–P2O5 glass. The various Si/P↔Na/Ca contacts were probed by molecular dynamics (MD) simulations together with heteronuclear magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) experimentation utilizing (23)Na{(31)P} and (23)Na{(29)Si} REDOR, as well as (31)P{ (23)Na} and (29)Si{(23)Na} REAPDOR. We introduce an approach for quantifying the extent of Na(+)/Ca(2+) ordering around a given QP(n) or QSi(n) group, encoded by the preference factor 0⩽ PM ⩽ 1 conveying the relative weights of a random cation intermixing (PM = 0) and complete preference/ordering (PM = 1) for one of the species M, which represents either Na(+) or Ca(2+). The MD-derived preference factors reveal phosphate and silicate species surrounded by Na(+)/Ca(2+) ions intermixed nearly randomly (PM ≲ 0.15), except for the QSi(4) and QSi(1) groups, which manifest more significant cation ordering with preference for Na+ and Ca2+, respectively. The overall weak preferences are essentially independent of the Si and P contents of the glass, whereas PM primarily correlates with the total amount of network modifiers: as the latter is increased, the Na/Ca distribution around the {QP(0), QSi(1), QSi(2)} groups with preference for Ca2(+ )tend to randomize (i.e., PCa decreases), while the PNa-values grow slightly for the {QP(1), QSi(3), QSi(4)} species already preferring coordination of Na. The set of experimental preference factors {PCa} for the orthophosphate (QP(0)) groups extracted from (31)P{(23)Na} REAPDOR NMR-derived M2(P–Na) dipolar second moments agrees well with the MD-generated counterparts. Our results on the Na/Ca intermixing in soda-lime-silicate glasses are discussed in relation to previous reports, highlighting the dependence of the conclusion on the approach to data evaluation.

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.

Atomic-scale models of early-stage alkali depletion and SiO2-rich gel formation in bioactive glasses

Molecular dynamics simulations of Na⁺/H⁺-exchanged 45S5 Bioglass® reveal the co-existence of bonded and non-bonded hydroxyls, suggesting a direct mechanism for forming a silica-rich gel structure upon the initial ion exchange.

Probing the fate of interstitial water in bulk bioactive glass by ab initio simulations

The mechanism and effects of the interaction of a water molecule with different sites found in the bulk of 45S5 bioactive glass have been investigated through ab initio simulations.

Current challenges in atomistic simulations of glasses for biomedical applications

Atomic-scale simulations of bioglasses are being used to tackle several challenging aspects, such as new structural markers of bioactivity, ion migration and nanosized samples.