Differential scanning calorimetry (DSC) is one of the most versatile probes for silicate glasses, allowing determination of, e.g., transition temperatures (glass, crystallization, melting) and the temperature dependence of heat capacity. However, complications arise for glasses featuring overlapping transitions and low sensitivity, e.g., arising from SiO2-rich compositions with small change in heat capacity during glass transition or the low sensitivity of thermocouples at high temperature. These challenges might be overcome using temperature-modulated DSC (TM-DSC), which enables separation of overlapping signals and improved sensitivity at the expense of increased measurement duration.

The existence of nanoscale ductility during the fracture of silicate glasses remains controversial. Here, based on molecular dynamics simulations coupled with topological constraint theory, we show that nano-ductility arises from the spatial heterogeneity of the atomic network's rigidity. Specifically, we report that localized floppy modes of deformation in under-constrained regions of the glass enable plastic deformations of the network, resulting in permanent change in bond configurations. Ultimately, these heterogeneous plastic events percolate, thereby resulting in a non-brittle mode of fracture. This suggests that nano-ductility is intrinsic to multi-component silicate glasses having nanoscale heterogeneities.

Although molecular dynamics (MD) simulations are commonly used to predict the structure and properties of glasses, they are intrinsically limited to short time scales, necessitating the use of fast cooling rates. It is therefore challenging to compare results from MD simulations to experimental results for glasses cooled on typical laboratory time scales. Based on MD simulations of a sodium silicate glass with varying cooling rate (from 0.01 to 100 K/ps), here we show that thermal history primarily affects the medium-range order structure, while the short-range order is largely unaffected over the range of cooling rates simulated. This results in a decoupling between the enthalpy and volume relaxation functions, where the enthalpy quickly plateaus as the cooling rate decreases, whereas density exhibits a slower relaxation. Finally, we show that, using the proper extrapolation method, the outcomes of MD simulations can be meaningfully compared to experimental values when extrapolated to slower cooling rates.