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.

Tests of mode coupling theory in a simple model for two-component miscible polymer blends

We present molecular dynamics simulations on the structural relaxation of a simple bead-spring model for polymer blends. The introduction of a different monomer size induces a large time scale separation for the dynamics of the two components. Simulation results for a large set of observables probing density correlations, Rouse modes, and orientations of bond and chain end-to-end vectors, are analyzed within the framework of the Mode Coupling Theory (MCT). An unusually large value of the exponent parameter is obtained. This feature suggests the possibility of an underlying higher-order MCT scenario for dynamic arrest.Comment: Revised version. Additional figures and citation

Lattice Boltzmann versus Molecular Dynamics simulation of nano-hydrodynamic flows

A fluid flow in a simple dense liquid, passing an obstacle in a two-dimensional thin film geometry, is simulated by Molecular Dynamics (MD) computer simulation and compared to results of Lattice Boltzmann (LB) simulations. By the appropriate mapping of length and time units from LB to MD, the velocity field as obtained from MD is quantitatively reproduced by LB. The implications of this finding for prospective LB-MD multiscale applications are discussed.Comment: 4 pages, 4 figure

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 Gaussian particle displacements.Comment: 24 pages, 14 figures, J. Phys.: Condens. Matter, in pres

Non-exponential kinetic behavior of confined water

We present the results of molecular dynamics simulations of SPC/E water confined in a realistic model of a silica pore. The single-particle dynamics have been studied at ambient temperature for different hydration levels. The confinement near the hydrophilic surface makes the dynamic behaviour of the liquid strongly dependent on the hydration level. Upon decrease of the number of water molecules in the pore we observe the onset of a slow dynamics due to the ``cage effect''. The conventional picture of a stochastic single-particle diffusion process thus looses its validity

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

Wall-liquid and wall-crystal interfacial free energies via thermodynamic integration: A molecular dynamics simulation study

A method is proposed to compute the interfacial free energy of a Lennard-Jones system in contact with a structured wall by molecular dynamics simulation. Both the bulk liquid and bulk face-centered-cubic crystal phase along the (111) orientation are considered. Our approach is based on a thermodynamic integration scheme where first the bulk Lennard-Jones system is reversibly transformed to a state where it interacts with a structureless flat wall. In a second step, the flat structureless wall is reversibly transformed into an atomistic wall with crystalline structure. The dependence of the interfacial free energy on various parameters such as the wall potential, the density and orientation of the wall is investigated. The conditions are indicated under which a Lennard-Jones crystal partially wets a flat wall.Comment: 15 pages, 11 figure

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) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and D"unweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely the effective Stokes friction and long time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.Comment: 37 pages, 12 figure

Tension and stiffness of the hard sphere crystal-fluid interface

A combination of fundamental measure density functional theory and Monte Carlo computer simulation is used to determine the orientation-resolved interfacial tension and stiffness for the equilibrium hard-sphere crystal-fluid interface. Microscopic density functional theory is in quantitative agreement with simulations and predicts a tension of 0.66 kT/\sigma^2 with a small anisotropy of about 0.025 kT and stiffnesses with e.g. 0.53 kT/\sigma^2 for the (001) orientation and 1.03 kT/\sigma^2 for the (111) orientation. Here kT is denoting the thermal energy and \sigma the hard sphere diameter. We compare our results with existing experimental findings

The Influence of Chemical Short Range Order on Atomic Diffusion in Al-Ni Melts

We use inelastic neutron scattering and molecular dynamics (MD) simulation to investigate the chemical short range order (CSRO), visible through prepeaks in the structure factors, and its relation to self diffusion in Al-Ni melts. As a function of composition at 1795K Ni self diffusion coefficients from experiment and simulation exhibit a non-linear dependence with a pronounced increase on the Al-rich side. This comes along with a change in CSRO with increasing Al content that is related to a more dense packing of the atoms in Ni-rich Al-Ni systems.Comment: 11 pages, 4 figure