The amorphous aluminium silicate (Al2O3)2(SiO2) [AS2] is investigated by means of large scale molecular dynamics computer simulations. We consider fully equilibrated melts in the temperature range 6100K >= T >= 2300K as well as glass configurations that were obtained from cooling runs from T=2300K to 300K with a cooling rate of about 10^12K/s. Already at temperatures as high as 4000K, most of the Al and Si atoms are four-fold coordinated by oxygen atoms. Thus, the structure of AS2 is that of a disordered tetrahedral network. The packing of AlO4 tetrahedra is very different from that of SiO4 tetrahedra in that Al is involved with a relatively high probability in small-membered rings and in triclusters in which an O atom is surrounded by four cations. We find as typical configurations two-membered rings with two Al atoms in which the shared O atoms form a tricluster. On larger length scales, the system shows a microphase separation in which the Al-rich network structure percolates through the SiO2 network. The latter structure gives rise to a prepeak in the static structure factor at a wavenumber q=0.5\\AA^{-1}. The comparison of experimental X-ray data with the results from the simulation shows a good agreement for the structure function. The diffusion dynamics in AS2 is found to be much faster than in SiO2. We show that the self-diffusion constants for O and Al are very similar and that they are by a factor of 2-3 larger than the one for Si.

We present the results of molecular dynamics computer simulations of a binary Lennard-Jones liquid confined between two parallel rough walls. These walls are realized by frozen amorphous configurations of the same liquid and therefore the structural properties of the confined fluid are identical to the ones of the bulk system. Hence this setup allows us to study how the relaxation dynamics is affected by the pure effect of confinement, i.e. if structural changes are completely avoided. We find that the local relaxation dynamics is a strong function of z, the distance of the particles from the wall, and that close to the surface the typical relaxation times are orders of magnitude larger than the ones in the bulk. Because of the cooperative nature of the particle dynamics, the slow dynamics also affects the dynamics of the particles for large values of z. Using various empirical laws, we are able to parameterize accurately the z-dependence of the generalized incoherent intermediate scattering function F_s(q,z,t) and also the spatial dependence of structural relaxation times. These laws allow us to determine various dynamical length scales and we find that their temperature dependence is compatible with an Arrhenius law. Furthermore, we find that at low temperatures time and space dependent correlation function fulfill a generalized factorization property similar to the one predicted by mode-coupling theory for bulk systems. For thin films and/or at sufficiently low temperatures, we find that the relaxation dynamics is influenced by the two walls in a strongly non-linear way in that the slowing down is much stronger than the one expected from the presence of only one confining wall. ....

We present extensive molecular dynamics simulations of the motion of a single linear rigid molecule in a two-dimensional random array of fixed overlapping disklike obstacles. The diffusion constants for the center of mass translation, D-CM, and for rotation, D-R, are calculated for a wide range of the molecular length, L, and the density of obstacles, rho. The obtained results follow a master curve Drho(mu)similar to(L(2)rho)(-nu) with an exponent mu=-3/4 and 1/4 for D-R and D-CM, respectively, that can be deduced from simple scaling and kinematic arguments. The nontrivial positive exponent nu shows an abrupt crossover at L(2)rho=zeta(1). For D-CM we find a second crossover at L(2)rho=zeta(2). The values of zeta(1) and zeta(2) correspond to the average minor and major axis of the elliptic holes that characterize the random configuration of the obstacles. A violation of the Stokes-Einstein-Debye relation is observed for L(2)rho>zeta(1), in analogy with the phenomenon of enhanced translational diffusion observed in supercooled liquids close to the glass transition temperature. (C) 2004 American Institute of Physics.

We present molecular dynamics simulations of the motion of a single rod in a two-dimensional random static array of discs. For long rods the mean-squared displacement of the centre of mass shows a cage effect similar to that observed in supercooled liquids or dense colloidal systems. We have determined the time-dependence of the non-Gaussian parameter for different rod lengths. It is found that the long-time regime is strongly non-Gaussian even on length scales of the order of 10-15 times the rod length, thus showing the heterogeneity of the dynamics on such length scales.

We use inelastic neutron scattering and molecular dynamics simulation to investigate the interplay between the structure and the fast sodium ion diffusion in various sodium silicates. With increasing temperature and decreasing density the structure factors exhibit an emerging prepeak around 0.9 Angstrom(-1). We show that this prepeak has its origin in the formation of sodium rich channels in the static structure. The channels serve as preferential ion conducting pathways in the relative immobile Si-O matrix. On cooling below the glass transition this intermediate range order is frozen in.