Sodium effect on static mechanical behavior of MD-modeled sodium silicate glasses

The structural and elastic mechanical properties of xNa(2)O-(100-x)SiO2 sodium silicate glasses were computed and analyzed at different scales, using atomistic simulations and coarse-grain methods based on physical principles. The numerical simulations were performed on large samples (similar to 100(3) angstrom(3) box size with similar to 70 000 atoms), and the results were compared to experimental measurements. It was shown that the cutoff in the non-Coulombic part of the empirical interactions affects the pressure/density relations. Therefore, this value was tuned to achieve the experimental density at ambient pressure. As a result we obtained realistic mechanical and structural properties as well. With this model, we analyzed the elastic response of the samples for different sodium content. We showed, that experimentally measured elastic moduli result from a succession of micro-plastic rearrangements that must be taken into account when calculating microscopic elastic moduli. Moreover, we investigated the size dependence of the elastic moduli, and we showed a strong connection between small scale heterogeneous elasticity and sodium repartition. The transition from small scale to large scale description of elasticity should involve an accurate description of the spatial organization of sodium ions inside the silica network.