Cooling Rate Dependence and Dynamic Heterogeneity Below the Glass Transition in a Lennard-Jones Glass

We investigate a binary Lennard-Jones mixture with molecular dynamics simulations. We consider first a system cooled linearly in time with the cooling rate gamma. By varying gamma over almost four decades we study the influence of the cooling rate on the glass transition and on the resulting glass. We find for all investigated quantities a cooling rate dependence; with decreasing cooling rate the system falls out of equilibrium at decreasing temperatures, reaches lower enthalpies and obtains increasing local order. Next we study the dynamics of the melting process by investigating the most immobile and most mobile particles in the glass. We find that their spatial distribution is heterogeneous and that the immobile/mobile particles are surrounded by denser/less dense cages than an average particle.