Relaxation in glassforming liquids and amorphous solids

The field of viscousliquid and glassysolid dynamics is reviewed by a process of posing the key questions that need to be answered, and then providing the best answers available to the authors and their advisors at this time. The subject is divided into four parts, three of them dealing with behavior in different domains of temperature with respect to the glass transition temperature, Tg,and a fourth dealing with “short time processes.” The first part tackles the high temperature regime T\u3eTg, in which the system is ergodic and the evolution of the viscousliquid toward the condition at Tg is in focus. The second part deals with the regime T∼Tg, where the system is nonergodic except for very long annealing times, hence has time-dependent properties (aging and annealing). The third part discusses behavior when the system is completely frozen with respect to the primary relaxation process but in which secondary processes, particularly those responsible for “superionic” conductivity, and dopart mobility in amorphous silicon, remain active. In the fourth part we focus on the behavior of the system at the crossover between the low frequency vibrational components of the molecular motion and its high frequency relaxational components, paying particular attention to very recent developments in the short time dielectric response and the high Qmechanical response

New fitting scheme to obtain effective potential from Car-Parrinello molecular dynamics simulations: Application to silica

A fitting scheme is proposed to obtain effective potentials from Car-Parrinello molecular dynamics (CPMD) simulations. It is used to parameterize a new pair potential for silica. MD simulations with this new potential are done to determine structural and dynamic properties and to compare these properties to those obtained from CPMD and a MD simulation using the so-called BKS potential. The new potential reproduces accurately the liquid structure generated by the CPMD trajectories, the experimental activation energies for the self-diffusion constants and the experimental density of amorphous silica. Also lattice parameters and elastic constants of alpha-quartz are well-reproduced, showing the transferability of the new potential.Comment: 6 pages, 5 figure

Frequency dependent specific heat of viscous silica

We apply the Mori-Zwanzig projection operator formalism to obtain an expression for the frequency dependent specific heat c(z) of a liquid. By using an exact transformation formula due to Lebowitz et al., we derive a relation between c(z) and K(t), the autocorrelation function of temperature fluctuations in the microcanonical ensemble. This connection thus allows to determine c(z) from computer simulations in equilibrium, i.e. without an external perturbation. By considering the generalization of K(t) to finite wave-vectors, we derive an expression to determine the thermal conductivity \lambda from such simulations. We present the results of extensive computer simulations in which we use the derived relations to determine c(z) over eight decades in frequency, as well as \lambda. The system investigated is a simple but realistic model for amorphous silica. We find that at high frequencies the real part of c(z) has the value of an ideal gas. c'(\omega) increases quickly at those frequencies which correspond to the vibrational excitations of the system. At low temperatures c'(\omega) shows a second step. The frequency at which this step is observed is comparable to the one at which the \alpha-relaxation peak is observed in the intermediate scattering function. Also the temperature dependence of the location of this second step is the same as the one of the $\alpha-$peak, thus showing that these quantities are intimately connected to each other. From c'(\omega) we estimate the temperature dependence of the vibrational and configurational part of the specific heat. We find that the static value of c(z) as well as \lambda are in good agreement with experimental data.Comment: 27 pages of Latex, 8 figure

Static and Dynamic Properties of a Viscous Silica Melt Molecular Dynamics Computer Simulations

We present the results of a large scale molecular dynamics computer simulation in which we investigated the static and dynamic properties of a silica melt in the temperature range in which the viscosity of the system changes from O(10^-2) Poise to O(10^2) Poise. We show that even at temperatures as high as 4000 K the structure of this system is very similar to the random tetrahedral network found in silica at lower temperatures. The temperature dependence of the concentration of the defects in this network shows an Arrhenius law. From the partial structure factors we calculate the neutron scattering function and find that it agrees very well with experimental neutron scattering data. At low temperatures the temperature dependence of the diffusion constants $D$ shows an Arrhenius law with activation energies which are in very good agreement with the experimental values. With increasing temperature we find that this dependence shows a cross-over to one which can be described well by a power-law, D\propto (T-T_c)^gamma. The critical temperature T_c is 3330 K and the exponent gamma is close to 2.1. Since we find a similar cross-over in the viscosity we have evidence that the relaxation dynamics of the system changes from a flow-like motion of the particles, as described by the ideal version of mode-coupling theory, to a hopping like motion. We show that such a change of the transport mechanism is also observed in the product of the diffusion constant and the life time of a Si-O bond, or the space and time dependence of the van Hove correlation functions.Comment: 30 pages of Latex, 14 figure