A direct numerical simulation method for complex modulus of particle dispersions

We report an extension of the smoothed profile method (SPM)[Y. Nakayama, K. Kim, and R. Yamamoto, Eur. Phys. J. E {\bf 26}, 361(2008)], a direct numerical simulation method for calculating the complex modulus of the dispersion of particles, in which we introduce a temporally oscillatory external force into the system. The validity of the method was examined by evaluating the storage $G'(\omega)$ and loss $G"(\omega)$ moduli of a system composed of identical spherical particles dispersed in an incompressible Newtonian host fluid at volume fractions of $\Phi=0$, 0.41, and 0.51. The moduli were evaluated at several frequencies of shear flow; the shear flow used here has a zigzag profile, as is consistent with the usual periodic boundary conditions

A Smooth Interface Method for Simulating Liquid Crystal Colloid Dispersions

A new method is presented for mesoscopic simulations of particle dispersions in liquid crystal solvents. It allows efficient first-principle simulations of the dispersions involving many particles with many-body interactions mediated by the solvents. Demonstrations have been performed for the aggregation of colloid dispersions in two-dimensional nematic and smectic-C* solvents neglecting hydrodynamic effects, which will be taken into account in the near future.Comment: 13 pages, 4 figure

Simulating (electro)hydrodynamic effects in colloidal dispersions: smoothed profile method

Previously, we have proposed a direct simulation scheme for colloidal dispersions in a Newtonian solvent [Phys.Rev.E 71,036707 (2005)]. An improved formulation called the ``Smoothed Profile (SP) method'' is presented here in which simultaneous time-marching is used for the host fluid and colloids. The SP method is a direct numerical simulation of particulate flows and provides a coupling scheme between the continuum fluid dynamics and rigid-body dynamics through utilization of a smoothed profile for the colloidal particles. Moreover, the improved formulation includes an extension to incorporate multi-component fluids, allowing systems such as charged colloids in electrolyte solutions to be studied. The dynamics of the colloidal dispersions are solved with the same computational cost as required for solving non-particulate flows. Numerical results which assess the hydrodynamic interactions of colloidal dispersions are presented to validate the SP method. The SP method is not restricted to particular constitutive models of the host fluids and can hence be applied to colloidal dispersions in complex fluids

Statistical Mechanical Calculation of Anisotropic Step Stiffness of a Two-Dimensional Hexagonal Lattice Gas Model with Next-Nearest-Neighbor Interactions: Application to Si(111) Surface

We study a two-dimensional honeycomb lattice gas model with both nearest- and next-nearest-neighbor interactions in a staggered field, which describes the surface of stoichiometrically binary crystal. We calculate anisotropic step tension, step stiffness, and equilibrium island shape, by an extended random walk method. We apply the results to Si(111) 7$\times$7 reconstructed surface and high-temperature Si(111) 1$\times$1 surface. We also calculate inter-step interaction coefficient.Comment: revised on May 29 1999: RevTeX v3.1, 10 pages with 9 figures (one figure added

Vicinal Surface with Langmuir Adsorption: A Decorated Restricted Solid-on-solid Model

We study the vicinal surface of the restricted solid-on-solid model coupled with the Langmuir adsorbates which we regard as two-dimensional lattice gas without lateral interaction. The effect of the vapor pressure of the adsorbates in the environmental phase is taken into consideration through the chemical potential. We calculate the surface free energy $f$, the adsorption coverage $\Theta$, the step tension $\gamma$, and the step stiffness $\tilde{\gamma}$ by the transfer matrix method combined with the density-matrix algorithm. Detailed step-density-dependence of $f$ and $\Theta$ is obtained. We draw the roughening transition curve in the plane of the temperature and the chemical potential of adsorbates. We find the multi-reentrant roughening transition accompanying the inverse roughening phenomena. We also find quasi-reentrant behavior in the step tension.Comment: 7 pages, 12 figures (png format), RevTeX 3.1, submitted to Phys. Rev.