Molecular Dynamics Computer Simulation of the Dynamics of Supercooled Silica

We present the results of a large scale computer simulation of supercooled silica. We find that at high temperatures the diffusion constants show a non-Arrhenius temperature dependence whereas at low temperature this dependence is also compatible with an Arrhenius law. We demonstrate that at low temperatures the intermediate scattering function shows a two-step relaxation behavior and that it obeys the time temperature superposition principle. We also discuss the wave-vector dependence of the nonergodicity parameter and the time and temperature dependence of the non-Gaussian parameter.Comment: 5 pages, Latex, 6 postscript figure

Role of saddles in mean-field dynamics above the glass transition

Recent numerical developments in the study of glassy systems have shown that it is possible to give a purely geometric interpretation of the dynamic glass transition by considering the properties of unstable saddle points of the energy. Here we further develop this program in the context of a mean-field model, by analytically studying the properties of the closest saddle point to an equilibrium configuration of the system. We prove that when the glass transition is approached the energy of the closest saddle goes to the threshold energy, defined as the energy level below which the degree of instability of the typical stationary points vanishes. Moreover, we show that the distance between a typical equilibrium configuration and the closest saddle is always very small and that, surprisingly, it is almost independent of the temperature

Effect of dynamic stall on the aerodynamics of vertical-axis wind turbines

Accurate simulations of the aerodynamic performance of vertical-axis wind turbines pose a significant challenge for computational fluid dynamics methods. The aerodynamic interaction between the blades of the rotor and the wake that is produced by the blades requires a high-fidelity representation of the convection of vorticity within the wake. In addition, the cyclic motion of the blades induces large variations in the angle of attack on the blades that can manifest as dynamic stall. The present paper describes the application of a numerical model that is based on the vorticity transport formulation of the Navier–Stokes equations, to the prediction of the aerodynamics of a verticalaxis wind turbine that consists of three curved rotor blades that are twisted helically around the rotational axis of the rotor. The predicted variation of the power coefficient with tip speed ratio compares very favorably with experimental measurements. It is demonstrated that helical blade twist reduces the oscillation of the power coefficient that is an inherent feature of turbines with non-twisted blade configurations

Inhomogeneous elastic response of silica glass

Using large scale molecular dynamics simulations we investigate the properties of the {\em non-affine} displacement field induced by macroscopic uniaxial deformation of amorphous silica,a strong glass according to Angell's classification. We demonstrate the existence of a length scale $\xi$ characterizing the correlations of this field (corresponding to a volume of about 1000 atoms), and compare its structure to the one observed in a standard fragile model glass. The "Boson-peak'' anomaly of the density of states can be traced back in both cases to elastic inhomogeneities on wavelengths smaller than $\xi$, where classical continuum elasticity becomes simply unapplicable

Molecular structural order and anomalies in liquid silica

The present investigation examines the relationship between structural order, diffusivity anomalies, and density anomalies in liquid silica by means of molecular dynamics simulations. We use previously defined orientational and translational order parameters to quantify local structural order in atomic configurations. Extensive simulations are performed at different state points to measure structural order, diffusivity, and thermodynamic properties. It is found that silica shares many trends recently reported for water [J. R. Errington and P. G. Debenedetti, Nature 409, 318 (2001)]. At intermediate densities, the distribution of local orientational order is bimodal. At fixed temperature, order parameter extrema occur upon compression: a maximum in orientational order followed by a minimum in translational order. Unlike water, however, silica's translational order parameter minimum is broad, and there is no range of thermodynamic conditions where both parameters are strictly coupled. Furthermore, the temperature-density regime where both structural order parameters decrease upon isothermal compression (the structurally anomalous regime) does not encompass the region of diffusivity anomalies, as was the case for water.Comment: 30 pages, 8 figure

Temperature dependence of spatially heterogeneous dynamics in a model of viscous silica

Molecular dynamics simulations are performed to study spatially heterogeneous dynamics in a model of viscous silica above and below the critical temperature of the mode coupling theory, $T_{MCT}$. Specifically, we follow the evolution of the dynamic heterogeneity as the temperature dependence of the transport coefficients shows a crossover from non-Arrhenius to Arrhenius behavior when the melt is cooled. It is demonstrated that, on intermediate time scales, a small fraction of oxygen and silicon atoms are more mobile than expected from a Gaussian approximation. These highly mobile particles form transient clusters larger than that resulting from random statistics, indicating that dynamics are spatially heterogeneous. An analysis of the clusters reveals that the mean cluster size is maximum at times intermediate between ballistic and diffusive motion, and the maximum size increases with decreasing temperature. In particular, the growth of the clusters continues when the transport coefficients follow an Arrhenius law. These findings imply that the structural relaxation in silica cannot be understood as a statistical bond breaking process. Though the mean cluster sizes for silica are at the lower end of the spectrum of values reported in the literature, we find that spatially heterogeneous dynamics in strong and fragile glass formers are similar on a qualitative level. However, different from results for fragile liquids, we show that correlated particle motion along quasi one-dimensional, string-like paths is of little importance for the structural relaxation in this model of silica, suggesting that string-like motion is suppressed by the presence of covalent bonds.Comment: 13 pages, 11 figure

What does the potential energy landscape tell us about the dynamics of supercooled liquids and glasses?

For a model glass-former we demonstrate via computer simulations how macroscopic dynamic quantities can be inferred from a PEL analysis. The essential step is to consider whole superstructures of many PEL minima, called metabasins, rather than single minima. We show that two types of metabasins exist: some allowing for quasi-free motion on the PEL (liquid-like), the others acting as traps (solid-like). The activated, multi-step escapes from the latter metabasins are found to dictate the slowing down of dynamics upon cooling over a much broader temperature range than is currently assumed

Gaussian excitations model for glass-former dynamics and thermodynamics

We describe a model for the thermodynamics and dynamics of glass-forming liquids in terms of excitations from an ideal glass state to a Gaussian manifold of configurationally excited states. The quantitative fit of this three parameter model to the experimental data on excess entropy and heat capacity shows that ``fragile'' behavior, indicated by a sharply rising excess heat capacity as the glass transition is approached from above, occurs in anticipation of a first-order transition -- usually hidden below the glass transition -- to a ``strong'' liquid state of low excess entropy. The dynamic model relates relaxation to a hierarchical sequence of excitation events each involving the probability of accumulating sufficient kinetic energy on a separate excitable unit. Super-Arrhenius behavior of the relaxation rates, and the known correlation of kinetic with thermodynamic fragility, both follow from the way the rugged landscape induces fluctuations in the partitioning of energy between vibrational and configurational manifolds. A relation is derived in which the configurational heat capacity, rather than the configurational entropy of the Adam Gibbs equation, controls the temperature dependence of the relaxation times, and this gives a comparable account of the experimental observations.Comment: 21 pp., 17 fig

Relaxation dynamics and ionic conductivity in a fragile plastic crystal

We report a thorough characterization of the dielectric relaxation behavior and the ionic conductivity in the plastic-crystalline mixture of 60% succinonitrile and 40% glutaronitrile. The plastic phase can be easily supercooled and the relaxational behavior is investigated by broadband dielectric spectroscopy in the liquid, plastic crystalline, and glassy crystal phases. The very pronounced alpha-relaxation found in the spectra is characterized in detail. From the temperature dependence of the alpha-relaxation time, a fragility parameter of 62 was determined making this material one of the most fragile plastic-crystalline glass formers. A well-pronounced secondary and faint indications for a third relaxation process were found, the latter most likely being of Johari-Goldstein type. In addition, relatively strong conductivity contributions were detected in the spectra exhibiting the typical features of ionic charge transport.Comment: 8 pages, 7 figure