No description available.

Owing to their excellent infrared transmittance and good rheological properties, selenium based chalcogenide glasses have been materials of choice for a number of technological applications. However, chalcogenide glasses can undergo substantial structural relaxation even at room temperature due to their low glass transition temperatures. The origins of these dynamic processes and their correlation to the glass structure is therefore of fundamental and practical interest. In particular, a deep understanding of the dynamic response near the glass transition region could help elucidate the mechanism of these structural relaxation processes. The correlation between structure and dynamic properties of selenium based glass systems were therefore investigated. NMR and Raman spectroscopy measurements reveal that the structure of AsₓSe₁₋ₓ glass follow the chain crossing model in selenium-rich glasses but contain increasing amounts of cage molecules in arsenic-rich compositions. This structural pattern leads to systematic extrema in physical properties at the stoichiometric composition As₄₀Se₆₀.The dynamic response of AsₓSe₁₋ₓ glasses investigated by heat capacity spectroscopy shows two minima in melt fragility as a function of composition which correlate well with the dimensionality of the glassy network. The structure evolves from 2D to 3D during crosslinking of selenium chains by arsenic but reduces into a 2D layer-like structure at the stoichiometric composition. Upon precipitation of arsenic-rich cages the network first reverts back to 3D and eventually becomes a mix of 2D and 0D structural units. The presence of molecular clusters in the network is evidenced by a strong bimodal dynamic response at high arsenic contents. NMR and Raman spectroscopy measurements of GeₓSe₁₋ₓ glasses suggest a structure composed of aggregated tetrahedral units and long selenium chains with little or no connectivity. Distinct dynamic responses of these two separated structural motifs are revealed by heat capacity spectroscopy. A non-Gaussian distribution of the imaginary heat capacity peak provided further evidence for the structural heterogeneity. This behavior is consistent with high temperature NMR measurements which show that the dynamic response of floppy selenium chains is distinct from that of rigid tetrahedral units. Finally, heat capacity spectroscopy applied to pure selenium provides strong evidence for the microscopic origin of the non-exponential structural relaxation, a universal feature of fragile glasses. The evolution of the imaginary heat capacity peak shape during annealing shows a non-monotonic trend which remarkably matches model predictions based on the enthalpy landscape. These results indicate that the non-exponential character of the relaxation process is linked to density fluctuations in the glass.

No description available.

No description available.

No description available.

No description available.

An increase in the co-operativity in the motion of particles and a growth of a suitably defined dynamical correlation length seem to be generic features exhibited by all liquids upon supercooling. These features have been observed both in experiments and in numerical simulations of glass-forming liquids. Specially designed NMR experiments have estimated that the rough magnitude of this correlation length is of the order of a few nanometers near the glass transition. Simulations also predict that there are regions in the system which are more liquid-like than other regions. A complete theoretical understanding of this behaviour is not available at present. In recent calculations, Berthier, Biroli and coworkers [1, 2] extended the simple mode coupling theory (MCT) to incorporate the effects of dynamic heterogeneity and predicted the existence of a growing dynamical correlation length associated with the cooperativity of the dynamics. MCT also predicts a power law divergence of different dynamical quantities at the mode coupling temperature and at temperatures somewhat higher than the mode coupling temperature, these predictions are found to be consistent with experimental and simulation results. The system size dependence of these quantities should exhibit finite size scaling (FSS) similar to that observed near a continuous phase transition in the temperature range where they show power law growth. Hence we have used the method of finite size scaling in the context of the dynamics of supercooled liquids. In chapter 2, we present the results of extensive molecular dynamics simulations of a model glass forming liquid and extract a dynamical correlation length ξ associated with dynamic heterogeneity by performing a detailed finite size scaling analysis of a four-point dynamic susceptibility χ4(t) [3] and the associated Binder cumulant. We find that although these quantities show the “normal” finite size scaling behaviour expected for a system with a growing correlation length, the relaxation time τ does not. Thus glassy dynamics can not be fully understood in terms of “standard” critical phenomena. Inspired by the success of the empirical Adam-Gibbs relation [4] which relates dynamics with the configurational entropy, we have calculated the configurational entropy for different system sizes and temperatures to explain the nontrivial scaling behaviour of the relaxation time. We find that the behaviour of the relaxation time τ can be explained in terms of the Adam-Gibbs relation [4] for all temperatures and system sizes. This observation raises serious questions about the validity of the mode coupling theory which does not include the effects of the potential energy (or free energy) landscape on the dynamics. On the other hand, in the “random first order transition” theory (RFOT), introduced by Wolynes and coworkers [5], the configurational entropy plays a central role in determining the dynamics. So we also tried to explain our simulation results in terms of RFOT. However, this interpretation has the drawback that the value of one of the exponents of this theory extracted from our numerical results does not satisfy an expected physical bound, and there is no clear explanation for the obtained values of other exponents. Thus we find puzzling values for the exponents relevant to the applicability of RFOT, which are in need of explanation. This can be due to the fact that RFOT focuses only near the glass transition, while all our simulation results are for temperatures far above the glass transition temperature (actually, above the mode coupling temperature). Interestingly, results similar to ours were obtained in a recent analysis [6] of experimental data near the laboratory glass transition, on a large class of glass-forming materials. Thus right now we do not have any theory which can explain our simulation data consistently from all perspectives. There have been some attempts to extend the RFOT analysis to temperatures above the mode coupling temperature [7, 8] and to estimate a length scale associated with the configurational entropy at such temperatures. We compare our results with the predictions arising from these analyses. In chapter 3, we present simulation results that suggest that finite size scaling analysis is probably the only feasible method for obtaining reliable estimates of the dynamical correlation length for supercooled liquids. As mentioned before, although there exists a growing correlation length, the behaviour of all measured quantities (specifically, the relaxation time) is not in accordance with the behaviour expected in “standard” critical phenomena. So one might suspect the results for the correlation length extracted from the scaling analysis. To find out whether the results obtained by doing finite size scaling are correct, we have done simulations of very large system sizes for the same model glass forming liquid. In earlier studies, the correlation length has been extracted from the wave vector dependence of the dynamic susceptibility in the limit of zero wave vector, but to estimate the correlation length with reasonable accuracy one needs data in the small wave vector range. This implies that one needs to simulate very large systems. But as far as we know, in all previous studies typical system sizes of the order of 10, 000 particles have been used to do this analysis. In this chapter we show by comparing results for systems of 28, 000 and 350, 000 particles that these previous estimates are not reliable. We also show that one needs to simulate systems with at least a million particles to estimate the correlation length correctly near the mode coupling temperature and this size increases with decreasing temperature. We compare the correlation length obtained by analyzing the wave vector dependence of the dynamic susceptibility for a 350, 000particle system with the results obtained from the finite size scaling analysis. We were only able to compare the results in the high temperature range due to obvious reasons. However the agreement in the high temperature range shows that the finite size scaling analysis is robust and also establishes the fact that finite size scaling is the only practical method to extract reliable correlation lengths in supercooled liquids. In chapter 4, we present a free energy landscape analysis of dynamic heterogeneity for a monodisperse hard sphere system. The importance of the potential energy landscape for particles interacting with soft potentials is well known in the glass community from the work of Sastry et al. [9] and others, but the hard sphere system which does not have any well defined potential energy landscape also exhibits similar slow dynamics in the high density limit. Thus it is not clear how to treat the hard sphere systems within the same energy landscape formalism. Dasgupta et al. [10, 11, 12, 13, 14, 15] showed that one can explain the slow dynamics of these hard core systems in term of a free energy landscape picture. They and other researchers showed that these system have many aperiodic local minima in its free energy landscape, with free energy lower than that of the liquid. Using the Ramkrishnan-Yussouff free energy functional, we have performed multi parameter variational minimizations to map out the detailed density distribution of glassy free energy minima. We found that the distribution of the widths of local density peaks at glassy minima is spatially heterogeneous. By performing hard sphere event driven molecular dynamics simulation, we show that there exists strong correlation between these density inhomogeneity and the local Debye-Waller factor which provides a measure of the dynamic heterogeneity observed in simulations. This result unifies the system of hard core particles with the other soft core particles in terms of a landscapebased description of dynamic heterogeneity. In chapter 5, we extend the same free energy analysis to a polydisperse system and show that there is a critical polydispersity beyond which the crystal state is not stable and glassy states are thermodynamically stable. We also found a reentrant behaviour in the liquid-solid phase transition within this free-energy based formalism. These results are in qualitative agreement with experimental observations for colloidal systems.