Structural origin of thermal shrinkage in soda-lime silicate glass below the glass transition temperature: A theoretical investigation by microsecond timescale molecular dynamics simulations.

Microscopic dynamical features in the relaxation of glass structures are one of the most important unsolved problems in condensed matter physics. Although the structural relaxation processes in the vicinity of glass transition temperature are phenomenologically expressed by the Kohlrausch-Williams-Watts function and the relaxation time can be successfully interpreted by Adam-Gibbs theory and/or Narayanaswamy's model, the atomic rearrangement, which is the origin of the volume change, and its driving force have not been elucidated. Using the microsecond time-scale molecular dynamics simulations, this study provides insights to quantitatively determine the origin of the thermal shrinkage below Tg in a soda-lime silicate glass. We found that during annealing below Tg, Na ions penetrate into the six-membered silicate rings, which remedies the acute O-O-O angles of the energetically unstable rings. The ring structure change makes the space to possess the cation inside the rings, but the ring volume is eventually reduced, which results in thermal shrinkage of the soda-lime silica glass. In conclusion, the dynamical structural relaxation due to the cation displacement evokes the overall volume relaxation at low temperature in the glassy material.

Disclosing crystal nucleation mechanism in lithium disilicate glass through molecular dynamics simulations and free-energy calculations.

Unraveling detailed mechanism of crystal nucleation from amorphous materials is challenging for both experimental and theoretical approaches. In this study, we have examined two methods to understand the initial stage of crystal precipitation from lithium disilicate glasses using molecular dynamics simulations. One of the methods is a modified exploring method to find structurally similar crystalline clusters in the glass models, enabling us to find three different embryos, such as Li2Si2O5 (LS2), Li2SiO3 (LS) and Li3PO4 (LP), in the 33Li2O·66SiO2·1P2O5 glass (LS2P1), in which P2O5 is added as a nucleating agent. Interestingly, LS2 and LP crystals were found inside the LS2P1 glass while LS crystal appeared on the glass surface, which agrees with experimental observations. The other method is free energy calculation using a subnano-scale spherical crystal embedded in the glass model. This method, which we called Free-Energy Seeding Method (FESM), allows us to evaluate free energy change as a function of crystal radius and to identify critical size of the crystal precipitation. The free energy profiles for LS and LS2 crystal nuclei in the LS2 glass models possess maximum energy at a critical radius as expected by classical nucleation theory. Furthermore, the critical radius and the energy barrier height agree well with recent experimental investigation, proving the applicability of this method to design glass-ceramics by atomistic modeling.

Structural origins of the Mixed Alkali Effect in Alkali Aluminosilicate Glasses: Molecular Dynamics Study and its Assessment.

The comprehension of the nonlinear effects provided by mixed alkali effect (MAE) in oxide glasses is useful to optimize glass compositions to achieve specific properties that depend on the mobility of ions, such as the chemical durability, glass transition temperature, viscosity and ionic conductivity. Although molecular dynamics (MD) simulations have already been applied to investigate the MAE on silicates, less effort has been devoted to study such phenomenon in mixed alkali aluminosilicate glasses where alkali cations can act both as modifiers, forming non-bridging oxygens and percolation channels, and as charge compensator of the AlO4- units present in the network. Moreover, the ionic conductivity has not been computed yet; thus, the accuracy of the atomistic simulations in reproducing the MAE on the property is still open to question. In this work, we have validated five major interatomic potentials for the classical MD simulations by modelling the structure, density, glass transition temperature and ionic conductivity for three aluminosilicate glasses, (25 - x)Na2O - x(K2O) - 10(Al2O3) - 65(SiO2) (x = 0, 12.5, 25). It was observed that only the core-shell (CS) polarizable force field well reproduces the experimentally measured MAE on Tg and the ionic conductivity as well as the higher conductivity of single sodium aluminosilicate glass at low temperature and the higher conductivity of single potassium aluminosilicate glass at high temperature. The MAE is related to the suppression of jump events of the alkaline ions between dissimilar sites in the percolation channels consisting of both sodium and potassium ions as in the case of alkaline silicates. The superior reproducibility of the CS potential is originated from the larger and the flexible ring structures due to the smaller Si-O-Si inter-tetrahedra angle, creating appropriate percolation channels for ion conductivity. We also report detailed assessments for using the potential models including the CS potential for investigating MAE on aluminosilicates.

Higher Toughness of Metal-nanoparticle-implanted Sodalime Silicate Glass with Increased Ductility.

In this report, we propose a novel framework for toughening brittle oxide glass originated from enhanced ductility by implanting a secondary material comprising different mechanical properties. To do so, copper-metal nanoparticles are implanted into the subsurface layer of commercial soda-lime silica glass by using the electrofloat method. The crack initiation load of the implanted glass is found to be comparable to the glass chemically strengthened in ordinary tempering conditions. By observing crack propagation and stress distribution from cross-section, it is found that the crack propagation stops within the metal nanoparticle implanted layer, due to the stress dissipation or relaxation. The copper-implanted glass shows improved toughness with decreased hardness. The toughening mechanism of the composite glass is theoretically studied using molecular dynamics calculations on an amorphous silica model with copper nanoparticles embedded, and Peridynamics fracture simulations for indentation on a glass sheet model whose surface was implicitly modeled as the copper-implanted oxide glass. The experimentally observed phenomena of intrinsic toughening were well explained by the series of the conducted simulations.

Small nanoparticles, surface geometry and contact forces.

In this molecular dynamics study, we examine the local surface geometric effects of the normal impact force between two approximately spherical nanoparticles that collide in a vacuum. Three types of surface geometries-(i) crystal facets, (ii) sharp edges, and (iii) amorphous surfaces of small nanoparticles with radii R<10 nm-are considered. The impact forces are compared with their macroscopic counterparts described by nonlinear contact forces based on Hertz contact mechanics. In our simulations, edge and amorphous surface contacts with weak surface energy reveal that the average impact forces are in excellent agreement with the Hertz contact force. On the other hand, facet collisions show a linearly increasing force with increasing compression. Our results suggest that the nearly spherical nanoparticles are likely to enable some nonlinear dynamic phenomena, such as breathers and solitary waves observed in granular materials, both originating from the nonlinear contact force.

Influence of material stretchability on the equilibrium shape of a Möbius band.

We use a two-dimensional discrete, lattice-based model to show that Möbius bands made with stretchable materials are less likely to crease or tear. This stems from a delocalization of twisting strain that occurs if stretching is allowed. The associated low-energy configurations provide strategic target shapes for the guided assembly of nanometer and micron scale Möbius bands. To predict macroscopic band shapes for a given material, we establish a connection between stretchability and relevant continuum moduli, leading to insight regarding the practical feasibility of synthesizing Möbius bands from materials with continuum parameters that can be measured experimentally or estimated by upscale averaging.

Rich collision dynamics of soft and sticky crystalline nanoparticles: numerical experiments.

A molecular dynamics study on the collisional dynamics of soft and sticky single face-centered cubic crystal nanoparticles is presented. The softness and stickiness of the nanoparticles are controlled by varying parameters in the Lennard-Jones potential that is used to describe the interatomic interactions. Softening of nanoparticles due to extensive plastic deformations is observed as was previously found in hard nanoparticles. Further, two primary plastic deformation modes, slip and twinning, of the nanoparticles are found to play important roles in the temperature dependence of the coefficient of restitution. Additionally, we observe the effects of surface roughness, facets, and edges in the collisional behaviors of the sticky nanoparticles in low-velocity collisions. Nevertheless, the Johnson-Kendall-Roberts theory for macroscopic adhesive bodies still remains valid in nearly spherical nanoparticles.

Strong plastic deformation and softening of fast colliding nanoparticles.

Nanoparticles, with sizes ranging between 1 and ∼102 nm, show dynamical properties distinctly different than those of bulk materials. Due to their large surface area to volume ratio, their properties often depend on length scales. We investigate the size and the collision velocity (vcoll) dependence of the coefficient of restitution (COR) for nanoparticles made of a face-centered cubic lattice of Lennard-Jones atoms via nonequilibrium molecular dynamics simulations. A sharp crossover between elastic collision and plastic collision occurs when vcoll=vY, where vY is the size-dependent yield velocity. For high-collision velocities the COR ∼vcoll-α, α∼1. This result is in agreement with recent small system simulations and with experiments and is distinct from the elasticity-theory-based result for COR for inelastic collisions which behaves as vcoll-α, with α=14. We find that the size-dependent critical vY approaches the theoretical constant value for macroscopic spheres as our particle sizes grow. Possible insights into the origins of α∼1 and the size dependence of the yield velocity are suggested. The work also suggests that sufficiently fast moving nanoparticles traveling through vacuum could be sticky and hence could be of potential interest in many applications.