skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Effects of system size and cooling rate on the structure and properties of sodium borosilicate glasses from molecular dynamics simulations

ORCiD logo [1]; ORCiD logo [1]
  1. Department of Materials Science and Engineering, University of North Texas, Denton, Texas 76203, USA
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 148; Journal Issue: 2; Related Information: CHORUS Timestamp: 2018-02-15 00:10:56; Journal ID: ISSN 0021-9606
American Institute of Physics
Country of Publication:
United States

Citation Formats

Deng, Lu, and Du, Jincheng. Effects of system size and cooling rate on the structure and properties of sodium borosilicate glasses from molecular dynamics simulations. United States: N. p., 2018. Web. doi:10.1063/1.5007083.
Deng, Lu, & Du, Jincheng. Effects of system size and cooling rate on the structure and properties of sodium borosilicate glasses from molecular dynamics simulations. United States. doi:10.1063/1.5007083.
Deng, Lu, and Du, Jincheng. Sun . "Effects of system size and cooling rate on the structure and properties of sodium borosilicate glasses from molecular dynamics simulations". United States. doi:10.1063/1.5007083.
title = {Effects of system size and cooling rate on the structure and properties of sodium borosilicate glasses from molecular dynamics simulations},
author = {Deng, Lu and Du, Jincheng},
abstractNote = {},
doi = {10.1063/1.5007083},
journal = {Journal of Chemical Physics},
number = 2,
volume = 148,
place = {United States},
year = {Sun Jan 14 00:00:00 EST 2018},
month = {Sun Jan 14 00:00:00 EST 2018}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on January 10, 2019
Publisher's Accepted Manuscript

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • 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 issuemore » 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{sup 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.« less
  • We have carried out classical molecular dynamics simulations in order to get insight into the atomistic mechanisms of the deformation during nanoindentation of the pristine and irradiated forms of a sodium borosilicate glass. In terms of the glass hardness, we have found that the primary factor affecting the decrease of hardness after irradiation is depolymerization rather than free volume, and we argue that this is a general trend applicable to other borosilicate glasses with similar compositions. We have analyzed the changes of the short- and medium-range structures under deformation and found that the creation of oxygen triclusters is an importantmore » mechanism in order to describe the deformation of highly polymerized borosilicate glasses and is essential in the understanding of the folding of large rings under stress. We have equally found that the less polymerized glasses present a higher amount of relative densification, while the analysis of bond-breaking during the nanoindentation has showed that shear flow is more likely to appear around sodium atoms. The results provided in this study can be proven to be useful in the interpretation of experimental results.« less
  • /sup 9/Be and /sup 19/F NMR techniques are used to gain information about the short-range order in two beryllium fluoride glasses: A simple BeF/sub 2/ glass and a NaFx2BeF/sub 2/ glass. The magnitude of the average quadrupole coupling constant for /sup 9/Be sites is deduced from the location of the first-order satellites and from low-frequency second-order broadening of the central transition of the /sup 9/Be spectrum. The average quadrupole coupling constant for the binary glass was a factor of 4 larger than for BeF/sub 2/ glass. It is argued that this is consistent with the existence of non-tetrahedrally-coordinated Be ionsmore » in the sodium fluoroberyllate glass. The central transition of the /sup 9/Be NMR spectrum for the nominally pure BeF/sub 2/ glass exhibited an anomalous narrow line which may be due to the presence of Be metallic clusters. The linewidth of the /sup 19/F resonance and the /sup 19/F spin-lattice relaxation time, T/sub 1/, of the two glasses were also studied as a function of temperature. These measurements indicate that there is some motion of the fluorines in the binary sample slightly above its transition temperature of 388 K. The two glasses were simulated by methods of molecular dynamics. It was found that whereas in simple BeF/sub 2/ glass virtually all of the Be ions were tetrahedrally coordinated by four F ions, in multicomponent glasses a large number of Be ions have five F ions in the first coordination shell. This is qualitatively consistent with the experimental observations. Meaningful quantitative predictions of the effective /sup 9/Be quadrupole coupling constants for the two glasses failed, however, because the simulated glasses are more disordered than glasses prepared in the laboratory. There is a broad distribution of F--Be--F angles in the simulated BeF/sub 2/ glass.« less
  • We have performed classical molecular dynamics simulations in order to study the changes under compression in the local and medium range structural properties of three sodium borosilicate glasses with varying sodium content. These glasses have been isostatically compressed up to 20 GPa and then decompressed in order to analyze the different mechanisms that affect densification, alongside with the permanent modifications of the structure after a full compression/decompression cycle. The results show that the atomic packing is the prominent characteristic that governs the amount of densification in the glass, as well as the setup of the permanent densification. During compression, themore » bulk modulus increases linearly up to approximately 15 GPa and more rapidly for higher pressures, a behavior which is reflected on the rate of increase of the average coordination for B and Na. Radial distribution functions at different pressures during the cycle help to quantify the amount of distortions in the elementary structural units, with a pronounced shortening of the Na–Na and Na–O bond lengths during compression. A subsequent decomposition of the glassy matrix into elementary Voronoi volumes verifies the high compressibility of Na-rich regions.« less