The Structure and Dynamics of Sodium Disilicate

We investigate the structure and dynamics of sodium disilicate by means of molecular dynamics computer simulation. We show that the structure is described by a partially destroyed tetrahedral SiO_4 network and a spherical super structure formed by the silicon and sodium atoms. The static structure factor of our simulation is in very good agreement with one from a neutron scattering experiment. For 1008 particles we find strong finite size effects in the dynamics which are due to the missing of modes contributing to the boson peak.Comment: 7 pages of Latex, 3 figure

Slow Dynamics in Ion-Conducting Sodium Silicate Melts: Simulation and Mode-Coupling Theory

A combination of molecular-dynamics (MD) computer simulation and mode-coupling theory (MCT) is used to elucidate the structure-dynamics relation in sodium-silicate melts (NSx) of varying sodium concentration. Using only the partial static structure factors from the MD as an input, MCT reproduces the large separation in relaxation time scales of the sodium and the silicon/oxygen components. This confirms the idea of sodium diffusion channels which are reflected by a prepeak in the static structure factors around 0.95 A^-1, and shows that it is possible to explain the fast sodium-ion dynamics peculiar to these mixtures using a microscopic theory.Comment: 4 pages, 4 figure

Liquid pair correlations in four spatial dimensions: Theory versus simulation

Using liquid integral equation theory, we calculate the pair correlations of particles that interact via a smooth repulsive pair potential in d = 4 spatial dimensions. We discuss the performance of different closures for the Ornstein-Zernike equation, by comparing the results to computer simulation data. Our results are of relevance to understand crystal and glass formation in high-dimensional systems

A Double-Transition Scenario for Anomalous Diffusion in Glass-Forming Mixtures

We study by molecular dynamics computer simulation a binary soft-sphere mixture that shows a pronounced decoupling of the species' long-time dynamics. Anomalous, power-law-like diffusion of small particles arises, that can be understood as a precursor of a double-transition scenario, combining a glass transition and a separate small-particle localization transition. Switching off small-particle excluded-volume constraints slows down, rather than enhances, small-particle transport. The data are contrasted with results from the mode-coupling theory of the glass transition

The fluid-fluid interface in a model colloid-polymer mixture: Application of grand canonical Monte Carlo to asymmetric binary mixtures

We present a Monte Carlo method to simulate asymmetric binary mixtures in the grand canonical ensemble. The method is used to study the colloid-polymer model of Asakura and Oosawa. We determine the phase diagram of the fluid-fluid unmixing transition and the interfacial tension, both at high polymer density and close to the critical point. We also present density profiles in the two-phase region. The results are compared to predictions of a recent density functional theory.Comment: 4 pages, 4 figure

Diffusion and Interdiffusion in Binary Metallic Melts

We discuss the dependence of self- and interdiffusion coefficients on temperature and composition for two prototypical binary metallic melts, Al-Ni and Zr-Ni, in molecular-dynamics (MD) computer simulations and the mode-coupling theory of the glass transition (MCT). Dynamical processes that are mainly entropic in origin slow down mass transport (as expressed through self diffusion) in the mixture as compared to the ideal-mixing contribution. Interdiffusion of chemical species is a competition of slow kinetic modes with a strong thermodynamic driving force that is caused by non-entropic interactions. The combination of both dynamic and thermodynamic effects causes qualitative differences in the concentration dependence of self-diffusion and interdiffusion coefficients. At high temperatures, the thermodynamic enhancement of interdiffusion prevails, while at low temperatures, kinetic effects dominate the concentration dependence, rationalized within MCT as the approach to its ideal-glass transition temperature $T_c$. The Darken equation relating self- and interdiffusion qualitatively reproduces the concentration-dependence in both Zr-Ni and Al-Ni, but quantitatively, the kinetic contributions to interdiffusion can be slower than the lower bound suggested by the Darken equation. As temperature is decreased, the agreement with Darken's equation improves, due to a strong coupling of all kinetic modes that is a generic feature predicted by MCT.Comment: 16 pages, 12 figure

Spinodal Decomposition in Thin Films: Molecular Dynamics Simulations of a Binary Lennard-Jones Fluid Mixture

We use molecular dynamics (MD) to simulate an unstable homogeneous mixture of binary fluids (AB), confined in a slit pore of width $D$. The pore walls are assumed to be flat and structureless, and attract one component of the mixture (A) with the same strength. The pair-wise interactions between the particles is modeled by the Lennard-Jones potential, with symmetric parameters that lead to a miscibility gap in the bulk. In the thin-film geometry, an interesting interplay occurs between surface enrichment and phase separation. We study the evolution of a mixture with equal amounts of A and B, which is rendered unstable by a temperature quench. We find that A-rich surface enrichment layers form quickly during the early stages of the evolution, causing a depletion of A in the inner regions of the film. These surface-directed concentration profiles propagate from the walls towards the center of the film, resulting in a transient layered structure. This layered state breaks up into a columnar state, which is characterized by the lateral coarsening of cylindrical domains. The qualitative features of this process resemble results from previous studies of diffusive Ginzburg-Landau-type models [S.~K. Das, S. Puri, J. Horbach, and K. Binder, Phys. Rev. E {\bf 72}, 061603 (2005)], but quantitative aspects differ markedly. The relation to spinodal decomposition in a strictly 2-$d$ geometry is also discussed.Comment: 37 pages, 11 figures, to appear in Phys. Rev.

The Dynamics of Silica Melts under High Pressure: Mode-Coupling Theory Results

The high-pressure dynamics of a computer-modeled silica melt is studied in the framework of the mode-coupling theory of the glass transition (MCT) using static-structure input from molecular-dynamics (MD) computer simulation. The theory reproduces the experimentally known viscosity minimum (diffusivity maximum) as a function of density or pressure and explains it in terms of a corresponding minimum in its critical temperature. This minimum arises from a gradual change in the equilibrium static structure which shifts from being dominated by tetrahedral ordering to showing the cageing known from high-density liquids. The theory is in qualitative agreement with computer simulation results.Comment: Presented at ESF EW Glassy Liquids under Pressure, to be published in Journal of Physic