Based on their connectivity, network glasses can be classified as flexible, stressed–rigid, or isostatic, if the number of topological constraints is lower, larger, or equal to the number of atomic degrees of freedom, respectively. Thanks to the absence of any stoichiometric requirement, the rigidity of glasses can be continuously tuned (e.g., from flexible to stressed–rigid) by changing their chemical composition. Interestingly, optimally-constrained isostatic glasses have been noted to exhibit unusual properties (e.g., nearly-reversible glass transition, low relaxation, desirable mechanical properties, etc.). Especially, the non-aging intermediate state features an almost vanished endotherm between the first and second heat scan across glass transition, providing a pathway for phase change material optimization in the application of non-volatile rewriteable media. However, the physical origin of the unusual behaviors and properties of isostatic glasses remain unclear.This thesis begins with investigating how the flexible-to-rigid transition in network glasses is encoded in their energy landscape based on molecular dynamics simulations. To this end, we introduce a simplified, yet realistic model of network glasses with varying connectivity. We characterize the topography of these glasses by adopting the activation-relaxation technique (ART), which enables a systematic search of saddle points and transition pathways in the energy landscape surface. We then demonstrate that the flexible-to-rigid transition arises from an interplay between low-energy saddle points (in flexible glasses) and topological frustration (in stressed–rigid glasses). Also, by utilizing the ring structure, we expand the transition correlation with ring size distribution. Meanwhile, we highlight the local heterogeneity with all the energy landscape features by dicing the model into small cubes. Comparing within a single glass helps exclude the effect of different configurations, further consolidating our conclusion on the physical origin of rigidity transition.
Finally, to explore the role of chemistry effect in rigidity transition, we compare the behavior of the simple connectivity model with a realistic GexSe1-x model. With the similar shape of enthalpic differences, the realistic model could reveal the effect of glass-forming ability with experimental results where the simple model fails. Overall, we have a clear pathway towards understanding the physical origin of rigidity transition of GexSe1-x glass.