Showing papers in "arXiv: Soft Condensed Matter in 1998"

How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm

Abstract: We construct an accurate estimate for the root mean square force error of the particle-particle-particle-mesh (P3M) algorithm by extending a single particle pair error measure which has been given by Hockney and Eastwood. We also derive an easy-to-use analytic approximation to the error formula. This allows a straightforward and precise determination of the optimal splitting parameter (as a function of system specifications and P3M parameters) and hence knowledge of the force accuracy prior to the actual simulation. The high quality of the estimate is demonstrated in several examples.

Influence of wetting properties on the hydrodynamic boundary condition at a fluid-solid interface

Abstract: It is well known that, at a macroscopic level, the boundary condition for a viscous fluid at a solid wall is one of "no-slip". The liquid velocity field vanishes at a fixed solid boundary. In this paper, we consider the special case of a liquid that partially wets the solid, i.e. a drop of liquid in equilibrium with its vapor on the solid substrate has a finite contact angle. Using extensive Non-Equilibrium Molecular Dynamics (NEMD) simulations, we show that when the contact angle is large enough, the boundary condition can drastically differ (at a microscopic level) from a "no-slip" condition. Slipping lengths exceeding 30 molecular diameters are obtained for a contact angle of 140 degrees, characteristic of Mercury on glass. On the basis of a Kubo expression for $\delta$, we derive an expression for the slipping length in terms of equilibrium quantities of the system. The predicted behaviour is in very good agreement with the numerical results for the slipping length obtained in the NEMD simulations. The existence of large slipping lentgh may have important implications for the transport properties in nanoporous media under such "nonwetting" conditions.

The Johnson-Segalman model with a diffusion term in Couette flow

Abstract: We study the Johnson-Segalman (JS) model as a paradigm for some complex fluids which are observed to phase separate, or ``shear-band'' in flow. We analyze the behavior of this model in cylindrical Couette flow and demonstrate the history dependence inherent in the local JS model. We add a simple gradient term to the stress dynamics and demonstrate how this term breaks the degeneracy of the local model and prescribes a much smaller (discrete, rather than continuous) set of banded steady state solutions. We investigate some of the effects of the curvature of Couette flow on the observable steady state behavior and kinetics, and discuss some of the implications for metastability.

Fluid membranes in hydrodynamic flow fields: Formalism and an application to fluctuating quasi-spherical vesicles in shear flow

Abstract: The dynamics of a single fluid bilayer membrane in an external hydrodynamic flow field is considered. The deterministic equation of motion for the configuration is derived taking into account both viscous dissipation in the surrounding liquid and local incompressibility of the membrane. For quasi-spherical vesicles in shear flow, thermal fluctuations can be incorporated in a Langevin-type equation of motion for the deformation amplitudes. The solution to this equation shows an overdamped oscillatory approach to a stationary tanktreading shape. Inclination angle and ellipiticity of the contour are determined as a function of excess area and shear rate. Comparisons to numerical results and experiments are discussed.

Large frequency range of negligible transmission in 1D photonic quantum well structures

Abstract: We show that it is possible to enlarge the range of low transmission in 1D photonic crystals by using photonic quantum well structures. If a defect is introduced in the photonic quantum well structures, defect modes with a very high quality factor may appear. The transmission of the defect mode is due to the coupling between the eigenmodes of the defect and those at the band edges of the constituent photonic crystals.

Glass Transition Temperature and Dynamics of $\alpha$-Process in Thin Polymer Films

Abstract: The glass transition temperature $T_g$ and the temperature $T_{\alpha}$ corresponding to the peak in the dielectric loss due to the $\alpha$-process have been simultaneously determined as functions of film thickness $d$ through dielectric measurements for thin films of polystyrene. A decrease of $T_g$ was observed with decreasing film thickness, while $T_{\alpha}$ was found to remain almost constant for $d > d_c$ and decrease drastically for $d < d_c$. Here, $d_c$ is a critical thickness dependent on molecular weight. The thickness dependence of $T_g$ is related to the distribution of relaxation times of the $\alpha$-process, not to the relaxation time itself.

Quantitative study of laterally inhomogeneous wetting films

Abstract: Based on a microscopic density functional theory we calculate the internal structure of the three-phase contact line between liquid, vapor, and a confining wall as well as the morphology of liquid wetting films on a substrate exhibiting a chemical step. We present a refined numerical analysis of the nonlocal density functional which describes the interface morphologies and the corresponding line tensions. These results are compared with those predicted by a more simple phenomenological interface displacement model. Except for the case that the interface exhibits large curvatures, we find that the interface displacement model provides a quantitatively reliable description of the interfacial structures.

Modeling Market Mechanism with Evolutionary Games

Abstract: This is an essay solicited by Europhysics News, published in its March/April 1998 issue with slight modifications. We outline some highlights of the econophysics models, especially the so-called Minority model of competition and evolution. Even without the usual math, this essay offers an analytical solution to the Minority model, revealing some key features of the solution.

Molecular-Dynamics Simulation of a Glassy Polymer Melt: Incoherent Scattering Function

Abstract: We present simulation results for a model polymer melt, consisting of short, nonentangled chains, in the supercooled state. The analysis focuses on the monomer dynamics, which is monitored by the incoherent intermediate scattering function. The scattering function is recorded over six decades in time and for many different wave-vectors. The lowest temperatures studied are slightly above the critical temperature of mode-coupling theory (MCT), which was determined from a quantitative analysis of the beta- and alpha-relaxations. We find evidence for the space-time factorization theorem in the beta-relaxation regime, and for the time-temperature superposition principle in the alpha-regime, if the temperature is not too close to the critical temperature. The wave-vector dependence of the nonergodicity parameter, of the critical amplitude, and the alpha-relaxation time are in qualitative agreement with calculations for hard spheres. For wave-vectors larger than the maximum of the structure factor the alpha-relaxation time already agrees fairly well with the asymptotic MCT-prediction. The behavior of the relaxation time at small wave-vectors can be rationalized by the validity of the Gaussian approximation and the value of the Kohlrausch stretching exponent.

Simple Model for Attraction between Like-Charged Polyions

Abstract: We present a simple model for the possible mechanism of appearance of attraction between like charged polyions inside a polyelectrolyte solution. The attraction is found to be short ranged, and exists only in presence of multivalent counterions. The attraction is produced by the correlations in the condensed layers of counterions surrounding each polyion, and appears only if the number of condensed counterions exceeds the threshold, $ n > Z/2 \alpha $, where $\alpha$ is the valence of counterions and $Z$ is the polyion charge.

Competition of mixing and segregation in rotating cylinders

Abstract: Using discrete element methods, we study numerically the dynamics of the size segregation process of binary particle mixtures in three-dimensional rotating drums, operated in the continuous flow regime. Particle rotations are included and we focus on different volume filling fractions of the drum to study the interplay between the competing phenomena of mixing and segregation. It is found that segregation is best for a more than half-filled drum due to the non-zero width of the fluidized layer. For different particle size ratios, it is found that radial segregation occurs for any arbitrary small particle size difference and the final amount of segregation shows a linear dependence on the size ratio of the two particle species. To quantify the interplay between segregation and mixing, we investigate the dynamics of the center of mass positions for each particle component. Starting with initially separated particle groups we find that no mixing of the component is necessary in order to obtain a radially segregated core.

Topological Inclusions in 2D Smectic-C Films

Abstract: In thin films of smectic-C liquid crystals, localized regions containing additional smectic layers form circular inclusions that carry a topological charge. Such inclusions nucleate a companion topological defect. These inclusion-defect pairs are modeled as topological dipoles within the context of a one-coupling constant approximation to the 2D Frank free energy. Deviations of the dipole direction from a preferred orientation cause the dipoles to acquire a logarithmic charge. Thermal fluctuations of the dipole direction are calculated and found to be large, scaling as the logarithm of the system size. In addition to dipole-dipole interactions arising from the topological charges, we also find that the thermal fluctuations of the dipole directions are coupled through a preference for global charge neutrality of the logarithmic charges.

Anomalous structure factor of dense star polymer solutions

Abstract: The core-core structure factor of dense star polymer solutions in a good solvent is shown theoretically to exhibit an unusual behaviour above the overlap concentration. Unlike usual liquids, these solutions display a structure factor whose first peak decreases by increasing density while the second peak grows. The scenario repeats itself with the subsequent peaks as the density is further enhanced. For low enough arm numbers $f$ ($f \leq 32$), various different considerations lead to the conclusion that the system remains fluid at all concentrations.

Random Polyelectrolytes and Polyampholytes in Solution

Abstract: The behavior of polyelectrolytes and polyampholytes in semi-dilute solutions is investigated theoretically. Various statistical charge distributions along the polyelectrolyte chains are considered: smeared, annealed, permuted and quenched. Annealed polyampholytes are also considered. Path integral formulation was used to derive mean field free energies for the different models. Self-consistent field equation is obtained for the polymer order parameter and a Poisson-Boltzmann like equation for the electrostatic potential. The random phase approximation is used to calculate the monomer-monomer structure factor S(q) for the different statistical charge distribution models. We show that in the annealed model, fluctuations of the the monomer charges contribute to the electrostatic screening in addition to the free ions in the solution. The strength of this screening depends on the variance of the monomer charge distribution and is especially important for polyampholytes in bad solvent conditions where the mesophase separation is enhanced. The ratio between the variance and the net average charge determines whether polyampholytes behave as polyelectrolytes or as neutral chains.

Investigating the influence of different thermodynamic paths on the structural relaxation in a glass forming polymer melt

Abstract: We present results from Molecular Dynamics simulations of the thermal glass transition in a dense polymer melt. In previous work we compared the simulation data with the idealized version of mode coupling theory (MCT) and found that the theory provides a good description of the dynamics above the critical temperature. In order to investigate the influence of different thermodynamic paths on the structural relaxation (alpha-process), we performed simulations for three different pressures and are thus able to give a sketch of the critical line of MCT in the pressure-temperature-plane [(p,T)-plane], where, according to the idealised version of MCT, an ergodic-nonergodic transition should occur. Furthermore, by cooling our system along two different paths (an isobar and an isochor), with the same impact point on the critical line, we demonstrate that neither the critical temperature nor the exponent gamma depend on the chosen path.

Solvation effects on kinetics of methylene chloride reactions in sub- and supercritical water: theory, experiment, and ab initio calculations

Abstract: The nature of the CH2Cl2 neutral/acidic hydrolysis reaction from ambient to supercritical conditions (25 C to 600 C at 246 bar) is explored. Of primary interest is the effect of the changing dielectric behavior of the water solvent over this temperature range on this reaction. Experiments reveal that significant CH2Cl2 hydrolysis occurs under subcritical temperatures, while relatively little hydrolysis occurs under supercritical conditions. These trends cannot be explain by simple Arrhenius behavior. A combination of Kirkwood theory and ab initio modeling provides a means of successfully accounting for this behavior both qualitatively and quantitatively. The results show that increases in the activation energy and a changing reaction profile with a decreasing dielectric constant provide a mechanism for a slowing of the reaction at higher temperatures by as much as three orders of magnitude. These solvent effects are captured quantitatively in a correction factor to the Arrenius form of the rate constant, which is incorporated into a global rate expression proposed for CH2Cl2 hydrolysis that provides good predictions of the experimental data.

Dynamics of spherical particles on a surface: About collision induced sliding and other effects

Abstract: We present a model for the motion of hard spherical particles on a two dimensional surface. The model includes both the interaction between the particles via collisions, as well as the interaction of the particles with the substrate. We analyze in details the effects of sliding and rolling friction, that are usually overlooked. It is found that the properties of this particulate system are influenced significantly by the substrate-particle interactions. In particular, sliding of the particles relative to the substrate after a collision leads to considerable energy loss for common experimental conditions. The presented results provide a basis, that can be used to realistically model the dynamical properties of the system, and provide further insight into density fluctuations and related phenomena of clustering, structure formations, and inelastic collapse.

Crystallization of a polymer on a surface

Abstract: We have studied the structure and free energy landscape of a semi-flexible lattice-polymer in the presence of a surface. At low temperatures coexistence of two-dimensional integer-folded crystals is observed. As the temperature is increased there is a transition from these crystalline configurations to a disordered coil adsorbed onto the surface. The polymer then gradually develops three-dimensional character at higher temperatures. We compute the free energy as a function of increasing crystallinity and compare with the free energy profiles assumed by the Lauritizen-Hoffman surface nucleation theory of polymer crystallization. Our free energy profiles exhibit a `sawtooth' structure associated with the successive formation of chain folds. However, in the early stages of crystallization our profiles significantly deviate from those assumed by surface nucleation theory because the initial nucleus is not a single stem but two incomplete stems connected by a fold. This finding has significant implications for the theoretical description of polymer crystallization.

DNA-Protein Cooperative Binding through Long-Range Elastic Coupling

Abstract: Cooperativity plays an important role in the action of proteins bound to DNA. A simple, mechanical mechanism for cooperativity, in the form of a tension-mediated interaction between proteins bound to DNA at two different locations is proposed. These proteins are not in direct physical contact. DNA segments intercalating bound proteins are modeled as a Worm-Like Chain, which is free to deform in two dimensions. The tension-controlled protein-protein interaction is the consequence of two effects produced by the protein binding. The first is the introduction of a bend in the host DNA and the second is the modification of the bending modulus of the DNA in the immediate vicinity of the bound protein. The interaction between two bound proteins may be either attractive or repulsive, depending on their relative orientation on the DNA. Applied tension controls both the strength and the range of protein-protein interactions in this model. Properties of the cooperative interaction are discussed, along with experimental implications.

Geodesics around line defects in elastic solids

Abstract: Topological defects in solids, usually described by complicated boundary conditions in elastic theory, may be described more simply as sources of a gravity- like deformation field in the geometric approach of Katanaev and Volovich. This way, the deformation field is described by non-Euclidean metric that incorporates the boundary imposed by the defects. A possible way of gaining some insight into the motion of particles in a medium with topological defects (e.g., electrons in a dislocated metal) is to look at the geodesics of the medium around the defect. In this work, we find the exact solution for the geodesic equation for elastic medium with a generic line defect, the dispiration, that can either be a screw dislocation or a wedge disclination for particular choices of its parameters.

Debye-Hueckel-Bjerrum theory for charged colloids

Abstract: We formulate an extension of the Debye-Hueckel-Bjerrum theory [M. E. Fisher and Y. Levin, Phys. Rev. Lett. 71, 3826 (1993)] to the fluid state of a highly asymmetric charged colloid. Allowing for the formation of clusters consisting of one polyion and n condensed counterions, the total Helmholtz free energy of the colloidal suspension is constructed. The thermodynamic properties, such as the cluster-density distribution and the pressure, are obtained by the minimization of the free energy under the constraints of fixed number of polyions and counterions. In agreement with the current experimental and Monte Carlo results, no evidence of any phase transition is encountered.

Monte Carlo simulations of copolymers at homopolymer interfaces: Interfacial structure as a function of the copolymer density

Abstract: By means of extensive Monte Carlo simulations of the bond fluctuation model, we study the effect of adding AB diblock copolymers on the properties of an interface between demixed homopolymer phases. The parameters are chosen such that the homopolymers are strongly segregated, and the whole range of copolymer concentrations in the two phase coexistence region is scanned. We compare the ``mushroom'' regime, in which copolymers are diluted and do not interact with each other, with the ``wet brush'' regime, where copolymers overlap and stretch, but are still swollen by the homopolymers. A ``dry brush'' regime is never entered for our choice of chain lengths. ``Intrinsic'' profiles are calculated using a block analysis method introduced by us in earlier work. We discuss density profiles, orientational profiles and contact number profiles. In general, the features of the profiles are similar at all copolymer concentrations, however, the profiles in the concentrated regime are much broader than in the dilute regime. The results compare well with self-consistent field calculations.

Monte Carlo and Molecular Dynamics Simulation of the Glass Transition of Polymers

Abstract: Two coarse-grained models for polymer chains in dense glass-forming polymer melts are studied by computer simulation: the bond-fluctuation model on a simple cubic lattice, where a bond-length potential favors long bonds, is treated by dynamic Monte Carlo methods, and a bead-spring model in the continuum with a Lennard-Jones potential between the beads is treated by Molecular Dynamics. While the dynamics of both models differ for short length scales and associated time scales, on mesoscopic spatial and temporal scales both models behave similarly. In particular, the mode coupling theory of the glass transition can be used to interpret the slowing down of the undercooled polymer melt. For the off-lattice model, the approach to the critical point of mode coupling is both studied for constant pressure and for constant volume. The lattice model allows a test of the Gibbs-Di Marzio entropy theory of the glass transition, and our finding is that although the entropy does decrease significantly, there is no ``entropy catastrophe'', where the fluid entropy would turn negative. Finally, an outlook on confinement effects on the glass transition in thin film geometry is given.

Theory and Phenomenology of Mixed Amphiphilic Aggregates

Abstract: We give a short overview of existing approaches describing shapes and energetics of amphiphilic aggregates. In particular, we consider recent experimental data and theory in relation to mixed aggregates. We point out the outstanding questions deserving further investigations such as stability of single-component vesicles and size growth of mixed vesicles induced by micelle-forming surfactants.

Determination of Interaction Potentials of Amino Acids from Native Protein Structures: Test on Simple Lattice Models

Abstract: We propose a novel method for the determination of the effective interaction potential between the amino acids of a protein. The strategy is based on the combination of a new optimization procedure and a geometrical argument, which also uncovers the shortcomings of any optimization procedure. The strategy can be applied on any data set of native structures such as those available from the Protein Data Bank (PDB). In this work, however, we explain and test our approach on simple lattice models, where the true interactions are known a priori. Excellent agreement is obtained between the extracted and the true potentials even for modest numbers of protein structures in the PDB. Comparisons with other methods are also discussed.

Fundamental measure theory for mixtures of parallel hard cubes. II. Phase behavior of the one-component fluid and of the binary mixture

Abstract: A previously developed fundamental measure fucntional [J. Chem. Phys. vol.107, 6379 (1997)] is used to study the phase behavior of a system of parallel hard cubes. The single-component fluid exhibits a continuous transition to a solid with an anomalously large density of vacancies. The binary mixture has a demixing transition for edge-length ratios below 0.1. Freezing in this mixture reveals that at least the phase rich in large cubes lies in the region where the uniform fluid is unstable, hence suggesting a fluid-solid phase separation. A method is develop to study very asymmetric binary mixtures by taking the limit of zero size ratio (scaling the density and fugacity of the solvent as appropriate) in the semi-grand ensemble where the chemical potential of the solvent is fixed. With this procedure the mixture is exactly mapped onto a one-component fluid of parallel adhesive hard cubes. At any density and solvent fugacity the large cubes are shown to collapse into a close-packed solid. Nevertheless the phase diagram contains a large metastability region with fluid and solid phases. Upon introduction of a slight polydispersity in the large cubes the system shows the typical phase diagram of a fluid with an isostructural solid-solid transition (with the exception of a continuous freezing). Consequences about the phase behavior of binary mixtures of hard core particles are then drawn.

Hydrodynamics of bacterial motion

Abstract: In this paper we present a hydrodynamic approach to describe the motion of migrating bacteria as a special class of self-propelled systems. Analytical and numerical calculations has been performed to study the behavior of our model in the turbulent-like regime and to show that a phase transition occurs as a function of noise strength. Our results can explain previous experimental observations as well as results of numerical simulations.

Laser induced quasicrystalline order in charge stabilised colloidal systems

Abstract: We have studied the ordering of a two dimensional charge stabilised colloidal system in the presence of a stationary one dimensionally modulated laser field formed by the superposition of two modulations with wavevectors $q_0 \tau$ and $q_0/\tau$, where $q_0$ is the wavevector corresponding to the first peak of the direct correlation function of the unperturbed liquid, and $\tau$ is the golden mean. In the framework of the Landau-Alexander-McTague theory we find that a decagonal quasicrystalline phase is stabler than the liquid or the triangular lattice in certain regions of the phase diagram. Our study also shows a reentrant melting phenomena for larger laser field strengths. We find that the transition from the modulated liquid to the quasicrystalline phase is continuous in contrast to a first order transition from the modulated liquid to the triangular crystalline phase.

Protein Adsorption on Lipid Monolayers at their Coexistence Region

Abstract: We investigate theoretically the behavior of proteins as well as other large macromolecules which are incorporated into amphiphilic monolayers at the air-water interface. We assume the monolayer to be in the coexistence region of the ``main'' transition, where domains of the liquid condensed phase coexist with the liquid expanded background. Using a simple mean-field free energy accounting for the interactions between proteins and amphiphilic molecules, we obtain the spatial protein distribution with the following characteristics. When the proteins preferentially interact with either the liquid condensed or liquid expanded domains, they will be dissolved in the respective phase. When the proteins are energetically rather indifferent to the density of the amphiphiles, they will be localized at the line boundary between the (two-dimensional) liquid expanded and condensed phases. In between these two limiting cases, a delocalization transition of the proteins takes place. This transition is accessible by changing the temperature or the amount of incorporated protein. These findings are in agreement with recent fluorescence microscopy experiments. Our results also apply to lipid multicomponent membranes showing coexistence of distinct fluid phases.

Strong-Segregation Theory of Bicontinuous Phases in Block Copolymers

Abstract: We compute phase diagrams for $A_nB_m$ starblock copolymers in the strong-segregation regime as a function of volume fraction $\phi$, including bicontinuous phases related to minimal surfaces (G, D, and P surfaces) as candidate structures. We present the details of a general method to compute free energies in the strong segregation limit, and demonstrate that the gyroid G phase is the most nearly stable among the bicontinuous phases considered. We explore some effects of conformational asymmetry on the topology of the phase diagram.