How do the properties of a glass depend on the cooling rate? A computer simulation study of a Lennard-Jones system
- 15 September 1996
- journal article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 105 (11) , 4714-4728
- https://doi.org/10.1063/1.472326
Abstract
Using molecular dynamics computer simulations we investigate how the glass transition and the properties of the resulting glass depend on the cooling rate with which the sample has been quenched. The system we study is a two component Lennard‐Jones model which is coupled to a heat bath whose temperature is decreased from a high temperature, where the system is a liquid, to zero temperature, where the system is a glass. The temperature Tb of this heat bath is decreased linearly in time, i.e. Tb=Ti−γt, where γ is the cooling rate, and we study the cooling rate dependence by varying γ over several orders of magnitude. In accordance with simple theoretical arguments and with experimental observations we find that the glass transition, as observed in the specific heat and the thermal expansion coefficient, becomes sharper when γ is decreased. A decrease of the cooling rate also leads to a decrease of the glass transition temperature Tg and we show that the dependence of Tg on γ can be rationalized by assuming that the temperature dependence of the relaxation times of the system is given by either a Vogel–Fulcher law or a power law. By investigating the structural properties of the glass, such as the radial distribution functions, the coordination numbers and the angles between three neighbor‐sharing particles, we show how the local order of the glass increases with decreasing cooling rate. The enthalpy H and the density ρ of the glass decrease and increase, respectively, with decreasing γ. By investigating the γ dependence of clusters of nearest neighbors, we show how the cooling rate dependence of H and ρ can be understood from a microscopic point of view. Furthermore we demonstrate that the frequency of icosahedral‐like structures is decreasing with decreasing cooling rate. We also show that the spectrum of the glass, as computed from the dynamical matrix, shows a shift towards higher frequencies when γ is decreased. All these effects show that there is a significant dependence of the properties of glasses on the cooling rate with which the glass is produced.Keywords
All Related Versions
This publication has 38 references indexed in Scilit:
- The role of localization in glasses and supercooled liquidsThe Journal of Chemical Physics, 1996
- A normal-mode study of a polymer glass containing a chromophore impurityThe Journal of Chemical Physics, 1996
- Characterization of free volume in atomic models of metallic glassesPhysical Review B, 1995
- Structure and dynamics of selenium chain melts: A molecular dynamics studyThe Journal of Chemical Physics, 1993
- Low-frequency excitations in glassy selenium: A comparison of neutron-scattering and molecular-dynamics resultsPhysical Review B, 1993
- Relaxation processes in supercooled liquidsReports on Progress in Physics, 1992
- Non-universal freezing in an alternating-bond Glauber chainJournal of Physics A: General Physics, 1991
- Freezing of nonequilibrium domain structures in a kinetic Ising modelJournal of Statistical Physics, 1991
- Residual Energies after Slow Cooling of Disordered SystemsPhysical Review Letters, 1986
- Short time structural relaxation processes in liquids: Comparison of experimental and computer simulation glass transitions on picosecond time scalesThe Journal of Chemical Physics, 1983