We use large-scale simulations to investigate the dynamic fracture of silica and sodium-silicate glasses under uniaxial tension. The stress-strain curves demonstrate that silica glass is brittle whereas the glasses rich in Na show pronounced ductility. A strong composition dependence is also seen in the crack velocity which is on the order of 1800 m/s for glasses with low Na concentration and decreases to 700 m/s if the concentration is high. We find that during the fracture of Na-rich glasses very irregular cavities as large as 3-4 nm form ahead of the crack front, indicating the presence of nanoductility in these glasses. Before fracture occurs, the local composition, structure, and mechanical properties are heterogeneous in space and show a strong dependence on the applied strain. Further analysis of the correlations between these local properties allows to obtain a better microscopic understanding of the deformation and fracture of glasses and how the local heating close to the crack tip, up to several hundred degrees, permits the structure to relax.

Quasi-two-dimensional (quasi-2D) colloidal hard-sphere suspensions confined in a slit geometry are widely used as two dimensional (2D) model systems in experiments that probe the glassy relaxation dynamics of 2D systems. However, the question to what extent these quasi-2D systems indeed represent 2D systems is rarely brought up. Here, we use computer simulations that take into account hydrodynamic interactions to show that dense quasi-2D colloidal bi-disperse hard-sphere suspensions exhibit much more rapid diffusion and relaxation than their 2D counterparts at the same area fraction. This difference is induced by the additional vertical space in the quasi-2D samples in which the small colloids can move out of the 2D plane, therefore allowing overlap between particles in the projected trajectories. Surprisingly, this difference in the dynamics can be accounted for if, instead of using the surface density, one characterizes the systems by means of a suitable structural quantity related to the radial distribution function. This implies that in the two geometries the relevant physics for glass-formation is essentially identical. Our results provide not only practical implications on 2D colloidal experiments but also interesting insights into the 3D-to-2D crossover in glass-forming systems.

For small tension the response of a solid to an applied stress is given by Hooke’s law. Outside this linear regime the relation between stress and strain is no longer universal and at present there is no satisfactory insight on how to connect for disordered materials the stress-strain relation to the microscopic properties of the system. Here we use atomistic computer simulations to establish this connection for the case of silicate glasses containing alkali modifiers. By probing how in the highly non-linear regime the stress-strain curve depends on composition, we are able to identify the microscopic mechanisms that are responsible for the complex dependence of stress on strain in these systems, notably the presence of an unexpected quasi-plateau in the tangent modulus. We trace back this dependence to the mobility of the modifiers which, without leaving their cage or modifying the topology of the network, are able to relieve the local stresses. Since the identified mechanism is general, the results obtained in this study will also be helpful for understanding the mechanical response of other disordered materials.