Structural and dynamical properties of sodium silicate melts: An investigation by molecular dynamics computer simulation

We present the results of large scale computer simulations in which we investigate the static and dynamic properties of sodium disilicate and sodium trisilicate melts. We study in detail the static properties of these systems, namely the coordination numbers, the temperature dependence of the Q^(n) species and the static structure factor, and compare them with experiments. We show that the structure is described by a partially destroyed tetrahedral SiO_4 network and the homogeneously distributed sodium atoms which are surrounded on average by 16 silicon and other sodium atoms as nearest neighbors. We compare the diffusion of the ions in the sodium silicate systems with that in pure silica a…

Simulations of Glassforming Network Fluids: Classical Molecular Dynamics versus Car-Parrinello Molecular Dynamics

Abstract Static and dynamic Properties of molten germanium dioxide are studied by two simulation methods, classical Molecular Dynamics (MD) using the Oeffner-Elliott (OE) potential, and “ab initio” Car-Parrinello Molecular Dynamics (CPMD). While CPMD provides a (presumably) more accurate description of the local structure and the forces, it severely suffers from finite size effects when the structure beyond the first neighbor shells is considered. For glassforming fluids, the demanding equilibrium needs are a further reason, why simply MD is still preferable, when a “good” effective potential is available.

The build-up and relaxation of stresses in a glass-forming soft-sphere mixture under shear: A computer simulation study

Molecular-dynamics computer simulations in conjunction with Lees-Edwards boundary conditions are used to investigate a glass-forming binary Yukawa fluid under shear. The transition from the elastic response to plastic flow is elucidated by studying the stress relaxation after switching off the shear. We find a slow stress relaxation starting from states in the elastic regime and a fast one starting from states in the plastic-flow regime. We show that these relaxation patterns are related to a different distribution of local microscopic stresses in both cases.

Amorphous silica modeled with truncated and screened Coulomb interactions: A molecular dynamics simulation study

We show that finite-range alternatives to the standard long-range BKS pair potential for silica might be used in molecular dynamics simulations. We study two such models that can be efficiently simulated since no Ewald summation is required. We first consider the Wolf method, where the Coulomb interactions are truncated at a cutoff distance r_c such that the requirement of charge neutrality holds. Various static and dynamic quantities are computed and compared to results from simulations using Ewald summations. We find very good agreement for r_c ~ 10 Angstroms. For lower values of r_c, the long--range structure is affected which is accompanied by a slight acceleration of dynamic properties…

Electrophoretic properties of charged colloidal suspensions: Application of a hybrid MD/LB method

Abstract Electrophoretic properties of charged colloidal suspensions are investigated using a hybrid simulation method. In this method, the colloidal particles are propagated via Newton’s equations of motion using molecular dynamics (MD), while they are coupled to a structureless solvent that is modelled by the Lattice-Boltzmann (LB) method.

Molecular Dynamics Simulations

A tutorial introduction to the technique of Molecular Dynamics (MD) is given, and some characteristic examples of applications are described. The purpose and scope of these simulations and the relation to other simulation methods is discussed, and the basic MD algorithms are described. The sampling of intensive variables (temperature T, pressure p) in runs carried out in the microcanonical (NVE) ensemble (N= particle number, V = volume, E = energy) is discussed, as well as the realization of other ensembles (e.g. the NVT ensemble). For a typical application example, molten SiO2, the estimation of various transport coefficients (self-diffusion constants, viscosity, thermal conductivity) is d…

Residual Stresses in Glasses

The history dependence of the glasses formed from flow-melted steady states by a sudden cessation of the shear rate $\dot\gamma$ is studied in colloidal suspensions, by molecular dynamics simulations, and mode-coupling theory. In an ideal glass, stresses relax only partially, leaving behind a finite persistent residual stress. For intermediate times, relaxation curves scale as a function of $\dot\gamma t$, even though no flow is present. The macroscopic stress evolution is connected to a length scale of residual liquefaction displayed by microscopic mean-squared displacements. The theory describes this history dependence of glasses sharing the same thermodynamic state variables, but differi…

Capillary Waves in a Colloid-Polymer Interface

The structure and the statistical fluctuations of interfaces between coexisting phases in the Asakura-Oosawa (AO) model for a colloid--polymer mixture are analyzed by extensive Monte Carlo simulations. We make use of a recently developed grand canonical cluster move with an additional constraint stabilizing the existence of two interfaces in the (rectangular) box that is simulated. Choosing very large systems, of size LxLxD with L=60 and D=120, measured in units of the colloid radius, the spectrum of capillary wave-type interfacial excitations is analyzed in detail. The local position of the interface is defined in terms of a (local) Gibbs surface concept. For small wavevectors capillary wa…

Specific Heat of Amorphous Silica within the Harmonic Approximation

We investigate to what extent the specific heat of amorphous silica can be calculated within the harmonic approximation. For this we use molecular dynamics computer simulations to calculate, for a simple silica model (the BKS potential), the velocity autocorrelation function and hence an effective density of states g(ν). We find that the harmonic approximation is valid for temperatures below 300 K but starts to break down at higher temperatures. We show that, to obtain a reliable description of the low-frequency part of g(ν), i.e., where the boson peak is observed, it is essential to use large systems for the simulations and small cooling rates to quench the samples. We find that the calcul…

Computer Simulations of the Dynamics of Amorphous Silica

We present the results of a large scale computer simulation we performed to investigate the dynamical properties of supercooled silica. We show that parallel supercomputers such as the CRAY-T3E are very well suited to solve these type of problems. We find that at low temperatures the transport properties such as the diffusion constants and the viscosity agree well with the experimental data. At high temperatures this simulation predicts that in the transport quantities significant deviations from the Arrhenius law should be observed. Finally we show that such types of simulations can be used to investigate also complex dynamical quantities, such as the dynamical structure factor, and that t…

ChemInform Abstract: Specific Heat of Amorphous Silica within the Harmonic Approximation.

We investigate to what extent the specific heat of amorphous silica can be calculated within the harmonic approximation. For this we use molecular dynamics computer simulations to calculate, for a simple silica model (the BKS potential), the velocity autocorrelation function and hence an effective density of states g(ν). We find that the harmonic approximation is valid for temperatures below 300 K but starts to break down at higher temperatures. We show that, to obtain a reliable description of the low-frequency part of g(ν), i.e., where the boson peak is observed, it is essential to use large systems for the simulations and small cooling rates to quench the samples. We find that the calcul…

Force-induced diffusion in microrheology

We investigate the force-induced diffusive motion of a tracer particle inside a glass-forming suspension when a strong external force is applied to the probe (active nonlinear microrheology). A schematic model of mode-coupling theory introduced recently is extended to describe the transient dynamics of the probe particle, and used to analyze recent molecular-dynamics simulation data. The model describes non-trivial transient displacements of the probe before a steady-state velocity is reached. The external force also induces diffusive motion in the direction perpendicular to its axis. We address the relation between the transverse diffusion coefficient D(perpendicular) and the force-depende…

From capillary condensation to interface localization transitions in colloid-polymer mixtures confined in thin-film geometry.

Monte Carlo simulations of the Asakura-Oosawa (AO) model for colloid-polymer mixtures confined between two parallel repulsive structureless walls are presented and analyzed in the light of current theories on capillary condensation and interface localization transitions. Choosing a polymer to colloid size ratio of q=0.8 and studying ultrathin films in the range of D=3 to D=10 colloid diameters thickness, grand canonical Monte Carlo methods are used; phase transitions are analyzed via finite size scaling, as in previous work on bulk systems and under confinement between identical types of walls. Unlike the latter work, inequivalent walls are used here: while the left wall has a hard-core rep…

Phase Behavior and Microscopic Transport Processes in Binary Metallic Alloys: Computer Simulation Studies

In a binary liquid mixture, different kinds of phase transitions can occur that are associated with various mass transport phenomena in the liquid. First, there is the possibility that the liquid undergoes a liquid-liquid demixing transition [1]. Near the critical point of this transition, a slowing down of dynamic properties is observed which is characterized, e.g., by a vanishing interdiffusion coefficient at the critical point [2, 3]. Another possible phase transition is a first-order transition of the liquid into a crystalline structure. In this case, crystal nucleation and growth are limited by the diffusive transport in the liquid [1, 4]. In a binary liquid, crystal nucleation process…

Onset of flow in a confined colloidal glass under an imposed shear stress.

A confined colloidal glass, under the imposition of a uniform shear stress, is investigated using numerical simulations. Both at macro- and microscales, the consequent dynamics during the onset of flow is studied. When the imposed stress is gradually decreased, the time scale for the onset of steady flow diverges, associated with long-lived spatial heterogeneities. Near this yield-stress regime, persistent creep in the form of shear-banded structures is observed.

Transport of Mobile Particles in an Immobile Environment: Computer Simulations of Sodium Silicates

Molecular dynamics (MD) simulations of various sodium silicate melts, (Na2O)x(SiO2) with x=2, 3, 20, are presented. In these systems, the mobility of sodium ions is much higher, often by orders of magnitude, than that of the silicon and oxygen atoms forming a tetrahedral network structure. We show that the high mobility of sodium is intimately related to the chemical ordering in sodium silicates. A network of percolating sodium-rich channels is formed in the static structure that serve as diffusion channels for the sodium ions. This channel network is revealed in static structure factors by a prepeak at the wavenumber q=0.95 A-1. Inelastic neutron scattering experiments of sodium silicate m…

Structure and transport properties of amorphous aluminium silicates: computer simulation studies

The structure and transport properties of SiO2-Al2O3 melts containing 13 mol% and 47 mol% Al2O3 are investigated by means of large scale molecular dynamics computer simulations. The interactions between the atoms are modelled by a pair potential which is a modified version of the one proposed by Kramer et al. [J. Am. Chem. Soc. 64, 6435 (1991)]. Fully equilibrated melts in the temperature range 6000 K >= T > 2000 K are considered as well as glass configurations, that were obtained by a rapid quench from the lowest melt temperatures. Each system is simulated at two different densities in order to study the effect of pressure on structural and dynamic properties. We find that the Al ato…

Simulation of Transport in Partially Miscible Binary Fluids: Combination of Semigrandcanonical Monte Carlo and Molecular Dynamics Methods

Binary Fluids that exhibit a miscibility gap are ubiquitous in nature (glass melts, polymer solutions and blends, mixtures of molten metals, etc.) and exhibit a delicate interplay between static and dynamic properties. This is exemplified for a simple model system, the symmetrical AB Lennard-Jones mixture. It is shown how semigrandcanonical Monte Carlo methods, that include A→B (B→A) identity switches as Monte Carlo moves, can yield the phase diagram, the interfacial tension between coexisting phases, and various pair correlation functions and structure factors. In addition to the build-up of long-ranged concentration correlations near the critical point, unmixing is also accompanied by the…

Nonlinear effects in charge stabilized colloidal suspensions

Molecular Dynamics simulations are used to study the effective interactions in charged stabilized colloidal suspensions. For not too high macroion charges and sufficiently large screening, the concept of the potential of mean force is known to work well. In the present work, we focus on highly charged macroions in the limit of low salt concentrations. Within this regime, nonlinear corrections to the celebrated DLVO theory [B. Derjaguin and L. Landau, Acta Physicochem. USSR {\bf 14}, 633 (1941); E.J.W. Verwey and J.T.G. Overbeck, {\em Theory of the Stability of Lyotropic Colloids} (Elsevier, Amsterdam, 1948)] have to be considered. For non--bulklike systems, such as isolated pairs or triples…

Electrokinetic Phenomena Revisited: A Lattice—Boltzmann Approach

The Lattice-Boltzmann method (LBM) is an efficient tool to solve the Navier-Stokes equations. Based on this method we have developed a scheme to investigate electrokinetic phenomena in charged colloidal suspensions. The equations of motion that are solved are the so-called electrokinetic equations, i.e. a set of partial differential equations that couple the gradient of the electrostatic potential to the hydrodynamic flow by means of a mean field theory. These equations have been extensively used to study electroviscous phenomena for the limit of a weakly charged sphere in an unbounded electrolyte. We demonstrate that our method can be applied beyond these limit. As an example we discuss th…

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 th…

Monte Carlo simulations of the solid-liquid transition in hard spheres and colloid-polymer mixtures

Monte Carlo simulations at constant pressure are performed to study coexistence and interfacial properties of the liquid-solid transition in hard spheres and in colloid-polymer mixtures. The latter system is described as a one-component Asakura-Oosawa (AO) model where the polymer's degrees of freedom are incorporated via an attractive part in the effective potential for the colloid-colloid interactions. For the considered AO model, the polymer reservoir packing fraction is eta_p^r=0.1 and the colloid-polymer size ratio is q=sigma_p/\sigma=0.15 (with sigma_p and sigma the diameter of polymers and colloids, respectively). Inhomogeneous solid-liquid systems are prepared by placing the solid fc…

Multiscale Computer Simulations in Physics, Chemistry, and Biology: The Example Of Silica

We show to what extent molecular dynamics simulations (MD) can explore struc-tural and dynamic properties of atomic systems whereby the system under consideration is amorphous silica (SiO2). Two studies are presented: (i) a large scale simulation of the dynam-ics of a SiO2 melt and (ii) the investigation of free silica surfaces where a mixture of a classical MD and a Car-Parrinello molecular dynamics is used.

Computer simulations of SiO2 and GeO2

Classical Molecular Dynamics (MD) simulations are used to study structural and dynamic properties of amorphous germania (GeO2) in comparison to those of silica (SiO2). The total structure factor, as obtained from these simulations, is in very good agreement with that of neutron scattering experiments, both for germania and silica. The tetrahedral network structure in silica and germania leads to a prepeak in the structure factor that appears at slightly smaller wavenumbers in GeO2 than in SiO2. At high temperatures the diffusion constants are very similar in both systems whereas at low temperatures diffusion is significantly faster in germania than in silica. We also outline the strategy fo…

Computer Simulation of Molten and Glassy Silica and its Mixtures with Sodium Oxide and Aluminium Oxide

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 fo…

Atomistic Simulation of Transport Phenomena in Simple and Complex Fluids and Fluid Mixtures

Computer simulations of fluids in thermal equilibrium can yield information on transport coefficients such as self—diffusion and interdiffusion coefficients, viscosity, and thermal conductivity. While the estimation of self—diffusion coefficients from the mean square displacements of the respective particles is rather straightforward, the estimation of other transport coefficients is less straightforward, and can be based on either an analysis of time correlation functions of the appropriate collective variables, or on nonequilibrium techniques where the linear response to appropriate perturbations is measured.

Critical phenomena in colloid-polymer mixtures: interfacial tension, order parameter, susceptibility, and coexistence diameter.

The critical behavior of a model colloid-polymer mixture, the so-called AO model, is studied using computer simulations and finite size scaling techniques. Investigated are the interfacial tension, the order parameter, the susceptibility and the coexistence diameter. Our results clearly show that the interfacial tension vanishes at the critical point with exponent 2\nu ~ 1.26. This is in good agreement with the 3D Ising exponent. Also calculated are critical amplitude ratios, which are shown to be compatible with the corresponding 3D Ising values. We additionally identify a number of subtleties that are encountered when finite size scaling is applied to the AO model. In particular, we find …

LATTICE–BOLTZMANN SIMULATION OF DENSE NANOFLOWS: A COMPARISON WITH MOLECULAR DYNAMICS AND NAVIER–STOKES SOLUTIONS

In a recent work, a dense fluid flow across a nanoscopic thin plate was simulated by means of Molecular Dynamics (MD) and Lattice Boltzmann (LB) methods. It was found that in order to recover quantitative agreement with MD results, the LB simulation must be pushed down to sub–nanoscopic scales, i.e. fractions of the range of molecular interactions. In this work, we point out that in this sub–nanoscopic regime, the LB method works outside the hydrodynamic limit at the level of a single cell spacing. A quantitative comparison with the Navier–Stokes (NS) solution shows however that LB and NS results are quite similar, thereby indicating that, apart for a small region past the plate, this nano…

Critical Dynamics in a Binary Fluid: Simulations and Finite-Size Scaling

We report comprehensive simulations of the critical dynamics of a symmetric binary Lennard-Jones mixture near its consolute point. The self-diffusion coefficient exhibits no detectable anomaly. The data for the shear viscosity and the mutual-diffusion coefficient are fully consistent with the asymptotic power laws and amplitudes predicted by renormalization-group and mode-coupling theories {\it provided} finite-size effects and the background contribution to the relevant Onsager coefficient are suitably accounted for. This resolves a controversy raised by recent molecular simulations.

Kinetics of phase separation in thin films: simulations for the diffusive case.

We study the diffusion-driven kinetics of phase separation of a symmetric binary mixture (AB), confined in a thin-film geometry between two parallel walls. We consider cases where (a) both walls preferentially attract the same component (A), and (b) one wall attracts A and the other wall attracts B (with the same strength). We focus on the interplay of phase separation and wetting at the walls, which is referred to as {\it surface-directed spinodal decomposition} (SDSD). The formation of SDSD waves at the two surfaces, with wave-vectors oriented perpendicular to them, often results in a metastable layered state (also referred to as ``stratified morphology''). This state is reminiscent of th…

Computer Simulations of Undercooled Fluids and Glasses

An introduction to the Molecular Dynamics (MD) simulation of chemically realistic models for undercooled fluids and glasses is given, emphasizing silicatic materials such as molten silicon dioxide and its mixtures with sodium oxide and aluminium oxide, and comparing the simulation results to experimental data whenever possible.

Monte Carlo simulations of phase transitions of systems in nanoscopic confinement

Abstract When simple or complex fluids are confined to ultrathin films or channels or other cavities of nanoscopic linear dimensions, the interplay of finite size and surface controls the phase behavior, and may lead to phase transitions rather different from the corresponding phenomena in the bulk. Monte Carlo simulation is a very suitable tool to clarify the complex behavior of such systems, since the boundary conditions providing the confinement can be controlled and arbitrarily varied, and detailed structural information on the inhomogeneous states of the considered systems is available. Examples used to illustrate these concepts include simple Ising models in pores and double-pyramid-s…

Structural relaxation in a binary metallic melt: Molecular dynamics computer simulation of undercooledAl80Ni20

Molecular dynamics computer simulations are performed to study structure and structural relaxation in the glassforming metallic alloy ${\text{Al}}_{80}{\text{Ni}}_{20}$. The interactions between the particles are modeled by an effective potential of the embedded atom type. Our model of ${\text{Al}}_{80}{\text{Ni}}_{20}$ exhibits chemical short-range order (CSRO) that is reflected in a broad prepeak around a wave number of $1.8\text{ }{\text{\AA{}}}^{\ensuremath{-}1}$ in the partial static structure factor for the Ni-Ni correlations. The CSRO is due to the preference of Ni atoms to have Al rather than Ni atoms as nearest neighbors. By analyzing incoherent and coherent intermediate scattering…

The relaxation dynamics of a viscous silica melt: II The intermediate scattering functions

We use molecular dynamics computer simulations to study the relaxation dynamics of a viscous melt of silica. The coherent and incoherent intermediate scattering functions, F_d(q,t) and F_s(q,t), show a crossover from a nearly exponential decay at high temperatures to a two-step relaxation at low temperatures. Close to the critical temperature of mode-coupling theory (MCT) the correlators obey in the alpha-regime the time temperature superposition principle (TTSP) and show a weak stretching. We determine the wave-vector dependence of the stretching parameter and find that for F_d(q,t) it shows oscillations which are in phase with the static structure factor. The temperature dependence of the…

A combined molecular dynamics and Monte Carlo study of the approach towards phase separation in colloid-polymer mixtures.

A coarse-grained model for colloid-polymer mixtures is investigated where both colloids and polymer coils are represented as point-like particles interacting with spherically symmetric effective potentials. Colloid-colloid and colloid-polymer interactions are described by Weeks-Chandler-Andersen potentials, while the polymer-polymer interaction is very soft, of strength k(B)T/2 for maximum polymer-polymer overlap. This model can be efficiently simulated both by Monte Carlo and molecular dynamics methods, and its phase diagram closely resembles that of the well-known Asakura-Oosawa model. The static and dynamic properties of the model are presented for systems at critical colloid density, va…

The Importance of Intermediate Range Order in Silicates: Molecular Dynamics Simulation Studies

We present the results of large scale computer simulations in which we investigate the structural and dynamic properties of silicate melts with the compositions (Na2O)2(SiO2) and (Al2O3)2(Si02). In order to treat such systems on a time scale of several nanoseconds and for system sizes of several thousand atoms it is necessary to use parallel supercomputers like the CRAY T3E. We show that the silicates under consideration exhibit additional intermediate range order as compared to silica (SiO2) where the characteristic intermediate length scales stem from the tetrahedral network structure. For the sodium silicate system it is demonstrated that the latter structural features are intimately con…

Simulation of surface-controlled phase separation in slit pores: Diffusive Ginzburg-Landau kinetics versus Molecular Dynamics

The phase separation kinetics of binary fluids in constrained geometry is a challenge for computer simulation, since nontrivial structure formation occurs extending from the atomic scale up to mesoscopic scales, and a very large range of time needs to be considered. One line of attack to this problem is to try nevertheless standard Molecular Dynamics (MD), another approach is to coarse-grain the model to apply a time-dependent nonlinear Ginzburg–Landau equation that is numerically integrated. For a symmetric binary mixture confined between two parallel walls that prefer one species, both approaches are applied and compared to each other. There occurs a nontrivial interplay between the forma…

O Triclusters Revisited:  Classical MD and Quantum Cluster Results for Glasses of Composition (Al2O3)2(SiO2)

The (17)O NMR spectrum of CaAl(2)Si(2)O(8) glass shows two types of O sites that are not present in the crystalline material. One of these, with (17)O NMR parameters C(Q) = 2.3 MHz and delta = +20 ppm, has been assigned to a "tricluster" O, a local geometry in which the O is coordinated to three tetrahedrally coordinated atoms, either Al or Si. For crystalline CaAl(4)O(7), a tricluster site (with three Al linkages to O, i.e., OAl(3)) has been characterized experimentally, with a C(Q) of 2.5 MHz and a delta of about +40 ppm. Thus, a C(Q) value of 2.5 MHz or less seems to be a characteristic of such sites, although they may show a range of delta values. However, several different quantum chem…

Channel Formation and Intermediate Range Order in Sodium Silicate Melts and Glasses

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 A^-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.

Negative thermal expansion of quartz glass at low temperatures: An ab initio simulation study

Abstract Using a mixed classical Molecular dynamics (MD)/ab initio simulation scheme combined with a quasi-harmonic approximation, we calculate the linear thermal expansion coefficient αL(T) in vitreous silica glasses. The systems are first cooled down by classical MD simulations. Then they are structurally relaxed by ab initio DFT calculations. The vibrational properties are calculated employing the frozen phonon method, and these results are finally used to calculate the Helmholtz free energy as a function of volume. In agreement with experiments, our simulations predict that αL(T) is negative at low temperatures up to T ≈ 150 K. In this low-temperature regime, the simulation results are …

Amorphous Silica at Surfaces and Interfaces: Simulation Studies

The structure of surfaces and interfaces of silica (SiO2) is investigated by large scale molecular dynamics computer simulations. In the case of a free silica surface, the results of a classical molecular dynamics simulation are compared to those of an ab initio method, the Car—Parrinello molecular dynamics. This comparative study allows to check the accuracy of the model potential that underlies the classical simulation. By means of a pure classical MD, the interface between amorphous and crystalline SiO2 is investigated, and as a third example the structure of a silica melt between walls is studied in equilibrium and under shear. We show that in the latter three examples important structu…

Surface-directed spinodal decomposition: Lattice model versus Ginzburg-Landau theory

When a binary mixture is quenched into the unstable region of the phase diagram, phase separation starts by spontaneous growth of long-wavelength concentration fluctuations ("spinodal decomposition"). In the presence of surfaces, the latter provide nontrivial boundary conditions for this growth. These boundary conditions can be derived from lattice models by suitable continuum approximations. But the lattice models can also be simulated directly, and thus used to clarify the conditions under which the Ginzburg–Landau type theory is valid. This comparison shows that the latter is accurate only in the immediate vicinity of the bulk critical point, if thermal fluctuations can also be neglecte…

Statics and dynamics of colloid-polymer mixtures near their critical point of phase separation: A computer simulation study of a continuous Asakura–Oosawa model

We propose a new coarse-grained model for the description of liquid-vapor phase separation of colloid-polymer mixtures. The hard-sphere repulsion between colloids and between colloids and polymers, which is used in the well-known Asakura-Oosawa (AO) model, is replaced by Weeks-Chandler-Anderson potentials. Similarly, a soft potential of height comparable to thermal energy is used for the polymer-polymer interaction, rather than treating polymers as ideal gas particles. It is shown by grand-canonical Monte Carlo simulations that this model leads to a coexistence curve that almost coincides with that of the AO model and the Ising critical behavior of static quantities is reproduced. Then the …

Molecular dynamics study of phase separation kinetics in thin films.

We use molecular dynamics to simulate experiments where a symmetric binary fluid mixture (AB), confined between walls that preferentially attract one component (A), is quenched from the one-phase region into the miscibility gap. Surface enrichment occurs during the early stages, yielding a B-rich mixture in the film center with well-defined A-rich droplets. The droplet size grows with time as l(t) proportional t(2/3) after a transient regime. The present atomistic model is also compared to mesoscopic coarse-grained models for this problem.

Kinetics of phase separation in thin films: Lattice versus continuum models for solid binary mixtures

A description of phase separation kinetics for solid binary (A,B) mixtures in thin film geometry based on the Kawasaki spin-exchange kinetic Ising model is presented in a discrete lattice molecular field formulation. It is shown that the model describes the interplay of wetting layer formation and lateral phase separation, which leads to a characteristic domain size $\ell(t)$ in the directions parallel to the confining walls that grows according to the Lifshitz-Slyozov $t^{1/3}$ law with time $t$ after the quench. Near the critical point of the model, the description is shown to be equivalent to the standard treatments based on Ginzburg-Landau models. Unlike the latter, the present treatmen…

The dynamics of sodium in sodium disilicate: Channel relaxation and sodium diffusion

We use molecular dynamics computer simulations to study the dynamics of amorphous (Na_2O)2(SiO_2). We find that the Na ions move in channels embedded in a SiO_2 matrix. The characteristic distance between these channels gives rise to a prepeak in the structure factor at around q=0.95 A^-1. The dynamics of sodium is given by a fast process which can be seen in the incoherent scattering function and a slow process which is seen in the coherent function. The relaxation time of the latter coincides with the alpha-relaxation time of the matrix. The Kohlrausch exponent of the fast process for q>1.6 A^1 is the same as the von Schweidler exponent for the slow one, demonstrating that the two proc…

Phase diagram and structure of colloid-polymer mixtures confined between walls

The influence of confinement, due to flat parallel structureless walls, on phase separation in colloid-polymer mixtures, is investigated by means of grand-canonical Monte Carlo simulations. Ultra-thin films, with thicknesses between $D=3-10$ colloid diameters, are studied. The Asakura-Oosawa model [J. Chem. Phys. 22, 1255 (1954)] is used to describe the particle interactions. To simulate efficiently, a ``cluster move'' [J. Chem. Phys. 121, 3253 (2004)] is used in conjunction with successive umbrella sampling [J. Chem. Phys. 120, 10925 (2004)]. These techniques, when combined with finite size scaling, enable an accurate determination of the unmixing binodal. Our results show that the critica…

Structure and dynamics of B2O3 melts and glasses: From ab initio to classical molecular dynamics simulations

Abstract Boron oxide (B2O3) is investigated by a combination of ab initio (DFT-based) molecular dynamics (MD) simulations and classical MD simulations. From the trajectories of the ab initio MD simulation, we derive a three-body interaction potential which is used in classical MD simulations to study various structural and dynamic properties on larger time and length scales than possible in the ab initio simulations. Differences and similarities to the structure and dynamics of other network glass formers such as SiO2 and GeO2 are discussed. Moreover, various properties as obtained from the simulations are compared to those from experiments of B2O3.

Finite-size scaling analysis of the anisotropic critical behavior of the two-dimensional Ising model under shear

The critical behavior of the two-dimensional Ising Model with non-conserved order parameter in steady-state shear is studied by large-scale Monte Carlo simulations. Studying the structure factor S(qx,qy) in the disordered phase, the ratio of correlation length exponents νx/νy in the two lattice directions (x,y) is estimated, and the critical temperature is determined as a function of the shear rate as Tc() − Tc(0) ∝ with ≈0.45. Critical exponents β≈0.37, γ≈1.1, ; ν⊥≈0.46, ν∥≈1.38 are roughly compatible with anisotropic hyperscaling.

Confinement effects on phase behavior of soft matter systems.

When systems that can undergo phase separation between two coexisting phases in the bulk are confined in thin film geometry between parallel walls, the phase behavior can be profoundly modified. These phenomena shall be described and exemplified by computer simulations of the Asakura-Oosawa model for colloid-polymer mixtures, but applications to other soft matter systems (e.g. confined polymer blends) will also be mentioned. Typically a wall will prefer one of the phases, and hence the composition of the system in the direction perpendicular to the walls will not be homogeneous. If both walls are of the same kind, this effect leads to a distortion of the phase diagram of the system in thin …

Molecular-dynamics computer simulation of crystal growth and melting in Al 50 Ni 50

The melting and crystallization of Al50Ni50} are studied by means of molecular dynamics computer simulations, using a potential of the embedded atom type to model the interactions between the particles. Systems in a slab geometry are simulated where the B2 phase of AlNi in the middle of an elongated simulation box is separated by two planar interfaces from the liquid phase, thereby considering the (100) crystal orientation. By determining the temperature dependence of the interface velocity, an accurate estimate of the melting temperature is provided. The value k=0.0025 m/s/K for the kinetic growth coefficient is found. This value is about two orders of magnitude smaller than that found in …

Modeling glass materials

Abstract Structural and dynamic properties of silicate melts and glasses (SiO 2 and its mixtures with Na 2 O and Al 2 O 3 ) are derived from Molecular Dynamics simulations and compared to pertinent experimental data. It is shown that these mixtures exhibit additional intermediate order as compared to pure silica, where the characteristic length scales stem from the tetrahedral network structure. While sodium ions show much faster diffusion through percolating channels than the silicon and oxygen ions forming the surrounding network, aluminium ions are incorporated into the network (leading to tricluster formation) and do not show such an enhanced mobility.

Some Finite Size Effects in Simulations of Glass Dynamics

We present the results of a molecular dynamics computer simulation in which we investigate the dynamics of silica. By considering different system sizes, we show that in simulations of the dynamics of this strong glass former surprisingly large finite size effects are present. In particular we demonstrate that the relaxation times of the incoherent intermediate scattering function and the time dependence of the mean squared displacement are affected by such finite size effects. By compressing the system to high densities, we transform it to a fragile glass former and find that for that system these types of finite size effects are much weaker.

Critical behavior of a colloid-polymer mixture confined between walls

We investigate the influence of confinement on phase separation in colloid-polymer mixtures. To describe the particle interactions, the colloid-polymer model of Asakura and Oosawa [J. Chem. Phys. 22, 1255 (1954)] is used. Grand canonical Monte Carlo simulations are then applied to this model confined between two parallel hard walls, separated by a distance D=5 colloid diameters. We focus on the critical regime of the phase separation and look for signs of crossover from three-dimensional (3D) Ising to two-dimensional (2D) Ising universality. To extract the critical behavior, finite size scaling techniques are used, including the recently proposed algorithm of Kim et al. [Phys. Rev. Lett. 91…

Computer simulation studies of finite-size broadening of solid–liquid interfaces: from hard spheres to nickel

Using Molecular Dynamics (MD) and Monte Carlo (MC) simulations interfacial properties of crystal-fluid interfaces are investigated for the hard sphere system and the one-component metallic system Ni (the latter modeled by a potential of the embedded atom type). Different local order parameters are considered to obtain order parameter profiles for systems where the crystal phase is in coexistence with the fluid phase, separated by interfaces with (100) orientation of the crystal. From these profiles, the mean-squared interfacial width w^2 is extracted as a function of system size. We rationalize the prediction of capillary wave theory that w^2 diverges logarithmically with the lateral size o…

Static and dynamic properties of a viscous silica melt

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}^{\ensuremath{-}2})$ P to ${O(10}^{2})$ P. 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 scatte…

O Triclusters Revisited: Classical MD and Quantum Cluster Results for Glasses of Composition (Al2O3)2(SiO2).

The (17)O NMR spectrum of CaAl(2)Si(2)O(8) glass shows two types of O sites that are not present in the crystalline material. One of these, with (17)O NMR parameters C(Q) = 2.3 MHz and delta = +20 ppm, has been assigned to a "tricluster" O, a local geometry in which the O is coordinated to three tetrahedrally coordinated atoms, either Al or Si. For crystalline CaAl(4)O(7), a tricluster site (with three Al linkages to O, i.e., OAl(3)) has been characterized experimentally, with a C(Q) of 2.5 MHz and a delta of about +40 ppm. Thus, a C(Q) value of 2.5 MHz or less seems to be a characteristic of such sites, although they may show a range of delta values. However, several different quantum chem…

Combining Molecular Dynamics with Lattice-Boltzmann: A Hybrid Method for the Simulation of (Charged) Colloidal Systems

We present a hybrid method for the simulation of colloidal systems, that combines molecular dynamics (MD) with the Lattice-Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point (with a velocity V(r) at any off-lattice position r is interacting with the neighboring eight LB nodes by a frictional force F=\xi_0(V(r)-u(r)) with \xi_0 being a friction force and u(r)…

The structural relaxation of molten sodium disilicate

We use molecular dynamics computer simulations to study the relaxation dynamics of Na2O-2(SiO2) in its molten, highly viscous state. We find that at low temperatures the incoherent intermediate scattering function for Na relaxes about 100 times faster than the one of the Si and O atoms. In contrast to this all coherent functions relax on the same time scale if the wave-vector is around 1AA^-1. This anomalous relaxation dynamics is traced back to the channel-like structure for the Na atoms that have been found for this system. We find that the relaxation dynamics for Si and O as well as the time dependence of the coherent functions for Na can be rationalized well by means of mode-coupling th…

Active nonlinear microrheology in a glass-forming Yukawa fluid.

A molecular dynamics computer simulation of a glass-forming Yukawa mixture is used to study the anisotropic dynamics of a single particle pulled by a constant force. Beyond linear response, a scaling regime is found where a force-temperature superposition principle of a Peclet number holds. In the latter regime, the diffusion dynamics perpendicular to the force can be mapped on the equilibrium dynamics in terms of an effective temperature, whereas parallel to the force a superdiffusive behavior is seen in the long-time limit. This behavior is associated with a hopping motion from cage to cage and can be qualitatively understood by a simple trap model.

The dynamics of melts containing mobile ions: computer simulations of sodium silicates

We present the results of large-scale computer simulations in order to discuss the structural and dynamic properties of sodium silicate melts with the compositions (Na2O)2(SiO2) (NS2) and (Na2O)20(SiO2) (NS20). We show that, compared to silica (SiO2), these systems exhibit additional intermediate range order on intermediate length scales that stem from the tetrahedral network structure. By means of intermediate-scattering functions, we characterize the dynamics of sodium in the system under consideration. Whereas in NS2 the incoherent scattering functions for Na decay much faster to zero than the coherent ones for Na–Na, in NS20 this different behaviour of the incoherent and coherent functi…

Qualitative characterisation of effective interactions of charged spheres on different levels of organisation using Alexander’s renormalised charge as reference

Abstract Effective interactions are conveniently determined from experimental or numerical data by fitting a Debye–Huckel potential with an effective charge Z ∗ and an effective electrolyte concentration c ∗ as free parameters. In this contribution we numerically solved the Poisson–Boltzmann equation to obtain the so-called renormalised charge Z PBC ∗ . For sufficiently large bare charge Z one finds a saturation of Z ∗ which scales as Z ∗ = A a / λ B , where a is the particle radius, λ B the Bjerrum length and A a proportionality factor of order (8–10). The saturation value increases with increased total micro-ion concentration and shows a shallow minimum as a function of packing fraction. …

From equilibrium to steady state: The transient dynamics of colloidal liquids under shear

We investigate stresses and particle motion during the start up of flow in a colloidal dispersion close to arrest into a glassy state. A combination of molecular dynamics simulation, mode coupling theory and confocal microscopy experiment is used to investigate the origins of the widely observed stress overshoot and (previously not reported) super-diffusive motion in the transient dynamics. A link between the macro-rheological stress versus strain curves and the microscopic particle motion is established. Negative correlations in the transient auto-correlation function of the potential stresses are found responsible for both phenomena, and arise even for homogeneous flows and almost Gaussia…

Intermediate Range Order in Silicate Melts and Glasses: Computer Simulation Studies

ABSTRACTWe present the results of large scale computer simulations to discuss the structural and dynamic properties of silicate melts with the compositions (Na2O)(2·SiO2), (Na2O)(20·SiO2) and (Al2O3)(2·SiO2). We show that these systems exhibit additional intermediate range order as compared to silica (SiO2) where the characteristic intermediate length scales stem from the tetrahedral network structure. Furthermore we show that the sodium dynamics in the sodium silicate systems exhibits a very peculiar feature: the long–time decay of the incoherent intermediate scattering function can be described by a Kohlrausch law with a constant exponent β for q > qth whereby qth is smaller than the l…