Compressible or incompressible blend of interacting monodisperse star and linear polymers near a surface

We consider a lattice model of a mixture of repulsive, attractive, or neutral monodisperse star (species A) and linear (species B) polymers with a third monomeric species C, which may represent free volume. The mixture is next to a hard, infinite plate whose interactions with A and C can be attractive, repulsive, or neutral. These two interactions are the only parameters necessary to specify the effect of the surface on all three components. We numerically study monomer density profiles using the method of Gujrati and Chhajer that has already been previously applied to study polydisperse and monodisperse linear-linear blends next to surfaces. The resulting density profiles always show an enrichment of linear polymers in the immediate vicinity of the surface due to entropic repulsion of the star core. However, the integrated surface excess of star monomers is sometimes positive, indicating an overall enrichment of stars. This excess increases with the number of star arms only up to a certain critical number and decreases thereafter. The critical arm number increases with compressibility (bulk concentration of C). The method of Gujrati and Chhajer is computationally ultrafast and can be carried out on a personal computer (PC), even in the incompressible case, when simulations are unfeasible. Calculations of density profiles usually take less than 20 min on PCs.     

Structure of inhomogeneous attractive and repulsive hard-core yukawa fluid: grand canonical Monte Carlo simulation and density functional theory study

A density functional theory is proposed for an inhomogeneous hard-core Yukawa (HCY) fluid based on Rosenfeld's perturbative method. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the hard-core repulsion and a quadratic functional Taylor expansion for the long-ranged attractive or repulsive interactions. To test the established theory, grand canonical ensemble Monte Carlo simulations are carried out to simulate the density profiles of attractive and repulsive HCY fluid near a wall. Comparison with the results from the Monte Carlo simulations shows that the present density functional theory gives accurate density profiles for both attractive and repulsive HCY fluid near a wall. Both the present theory and simulations suggest that there is depletion for attractive HCY fluid at low temperature, but no depletion is found for repulsive HCY fluid. The calculated results indicate that the present density functional theory is better than those of the modified version of the Lovett-Mou-Buff-Wertheim and other density functional theories. The present theory is simple in form and computationally efficient. It predicts accurate radial distribution functions of both attractive and repulsive HCY fluid except for the repulsive case at high density, where the theory overestimates the radial distribution function in the vicinity of contact.     

Analytic model for surface-induced molecular weight segregation in thermoplastic composites

An analytic model based on a random-walk formulation for absorbing barriers is presented. The range of applicability of the model is established from results derived from Monte Carlo simulations. The model takes into account the influences of neutral, repulsive, and attractive polymer-surface interactions on the spatial variation of species concentrations in the vicinity of a surface. The model predicts that the lower molecular weight species of a polydisperse mixture will accumulate in the vicinity of a surface to create an interphase region. The extent of this surface-perturbed interphase region is on the order of three times the radius of gyration of the longest species of the mixture. The degree of molecular weight segregation is dependent upon the breadth of the molecular weight distribution. © 1994 John Wiley & Sons. Inc.     

Comparison of random-walk density functional theory to simulation for bead-spring homopolymer melts

Density profiles for a homopolymer melt near a surface are calculated using a random-walk polymeric density functional theory, and compared to results from molecular dynamics simulations. All interactions are of a Lennard-Jones form, for both monomer-monomer interactions and surface-monomer interactions, rather than the hard core interactions which have been most investigated in the literature. For repulsive systems, the theory somewhat overpredicts the density oscillations near a surface. Nevertheless, near quantitative agreement with simulation can be obtained with an empirical scaling of the direct correlation function. Use of the random phase approximation to treat attractive interactions between polymer chains gives reasonable agreement with simulation of dense liquids near neutral and attractive surfaces.     

Density profiles and solvation forces for a Yukawa fluid in a slit pore

The effect of varying wall-particle and particle-particle interactions on the density profiles near a single wall and the solvation forces between two walls immersed in a fluid of particles is investigated by grand canonical Monte Carlo simulations. Attractive and repulsive particle-particle and particle-wall interactions are modeled by a versatile hard-core Yukawa form. These simulation results are compared to theoretical calculations using the hypernetted chain integral equation technique, as well as with fundamental measure density functional theory (DFT), where particle-particle interactions are either treated as a first order perturbation using the radial distribution function or else with a DFT based on the direct-correlation function. All three theoretical approaches reproduce the main trends fairly well, but exhibit inconsistent accuracy, particularly for attractive particle-particle interactions. We show that the wall-particle and particle-particle attractions can couple together to induce a nonlinear enhancement of the adsorption and a related "repulsion through attraction" effect for the effective wall-wall forces. We also investigate the phenomenon of bridging, where an attractive wall-particle interaction induces strongly attractive solvation forces.     

Modeling of adsorption in pores with strongly heterogeneous walls: parametric lattice-site wall model

We present results of grand canonical Monte Carlo simulations of adsorption in cylindrical pores with rough surface modeled by a parametric lattice-site approach. The sites are randomly distributed over the pore walls. They could be attractive, neutral or repulsive with respect to the smooth pore model. Each site is characterized by two amplitudes (structural and energetic) which modify locally the structure and energetic properties of the surface. The results presented here show how different parameters of the model affect the mechanism of adsorption and, consequently, the form of the isotherm.     

Monte Carlo study of polyelectrolyte gels

The discontinuous transition between dense and dilute phases in polyelectrolyte gels is observed in Bond-Fluctuation Method Monte Carlo simulations of gels. The transition is driven by the competition between local attractive interactions of a poor-quality solvent and global repulsive interactions from counter-ion pressure. A procedure is introduced that prevents local attractive interactions from destroying ergodicity. Under good solvent conditions, lengths and volumes of gels are found to follow self-avoiding random walk scaling. © 1996 John Wiley & Sons, Inc.     

Influence of perpendicular external magnetic field on microstructures of monolayer composed of ferromagnetic particles: analysis by means of quasi-two-dimensional Monte Carlo simulation

We investigated the influences of the magnetic field strength and particle areal density on the microstructure of a quasi-two-dimensional monolayer composed of ferromagnetic particles by means of a Monte Carlo simulation. The magnetic field was applied along a direction perpendicular to the plane of the monolayer. Microstructures of the monolayer obtained in the simulations were analyzed in terms of radial distribution and orientational distribution functions. Formation of the microstructures is discussed from the perspective of particle-particle interaction energy and the perpendicular magnetic susceptibility of the monolayer was calculated from simulated magnetization curves. The obtained results are summarized as follows. For small areal density of particles, formation of chain-like structures is prevented by the repulsive magnetic interaction between particles due to orientations of the magnetic moments in the particles along the magnetic field direction. For intermediate areal density of particles, the chain-like structures remain even when a relatively strong magnetic field is applied, because contributions of the attractive magnetic interactions increase. For large areal density of particles, mixtures of chain-like and locally ordered structures appear due to the anisotropic attractive magnetic interactions in the absence of the magnetic field. However, when a sufficiently strong magnetic field is applied, the magnetic interactions between particles change to isotropic repulsive interactions, which results in the short-range repulsive steric interactions between particles becoming dominant with the appearance of hexagonal close packed structures.     

Transitions of tethered polymer chains: a simulation study with the bond fluctuation lattice model

A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensional state space. The simulations' density-of-states results have been evaluated for interaction parameters spanning the range from good- to poor-solvent conditions and from repulsive to strongly attractive surfaces. In this work, we describe the simulation method and present results for the overall phase behavior and for some of the transitions. For adsorption in good solvent, we compare with Metropolis Monte Carlo data for the same model and find good agreement between the results. For the collapse transition, which occurs when the solvent quality changes from good to poor, we consider two situations corresponding to three-dimensional (hard surface) and two-dimensional (very attractive surface) chain conformations, respectively. For the hard surface, we compare tethered chains with free chains and find very similar behavior for both types of chains. For the very attractive surface, we find the two-dimensional chain collapse to be a two-step transition with the same sequence of transitions that is observed for three-dimensional chains: a coil-globule transition that changes the overall chain size is followed by a local rearrangement of chain segments.     

基于密度泛函理论研究二元排斥Yukawa流体的表面结构性质

基于自由能密度泛函理论(DFT)考察了二元排斥Yukawa (HCRY)流体在不同外场下的密度分布. 基于微扰理论, 体系的Helmholtz自由能泛函采用硬球排斥部分和长程色散部分贡献之和, 其中Kierlik和Rosinberg的加权密度近似(WDA)被用来计算硬球排斥部分, 而色散部分采用平均场理论(MFT)进行描述. 为了验证DFT计算结果的合理性, 研究中采用巨正则Monte Carlo(GCMC)模拟计算了在不同主体相密度、硬核直径和位能参数比的条件下二元HCRY混合流体的密度分布. 结果表明, 该DFT计算结果与GCMC模拟值吻合良好.     

Simple model for grafted polymer brushes

The first theories of grafted polymer brushes assumed a step profile for the monomer density. Later, the real density profile was obtained from Monte Carlo or molecular dynamics simulations and calculated numerically using a self-consistent field theory. The analytical approximations of the solutions of the self-consistent field equations provided a parabolic dependence of the self-consistent field, which in turn led to a parabolic distribution for the monomer density in neutral brushes. As shown by numerical simulations, this model is not accurate for dense polymer brushes, with highly stretched polymers. In addition, the scaling laws obtained from the analytical approximations of the self-consistent field theory are identical to those derived from the earlier step-profile-approximation and predict a vanishing thickness of the brush at low graft densities, and a thickness exceeding the length of the polymer chains at high graft densities. Here a simple model is suggested to calculate the monomer density and the interaction between surfaces with grafted polymer brushes, based on an approximate calculation of the partition function of the polymer chains. The present model can be employed for both good and poor solvents, is compatible with a parabolic-like profile at moderate graft densities, and leads to an almost steplike density for highly stretched brushes. While the thickness of the brush depends strongly on solvent quality, it is a continuous function in the vicinity of the temperature. In good and moderately poor solvents, the interactions between surfaces with grafted polymer brushes are always repulsive, whereas in poor solvents the interactions are repulsive at small separations and become attractive at intermediate separation distances, in agreement with experiment. At large separations, a very weak repulsion is predicted.     

A minimal proteinlike lattice model: an alpha-helix motif

A simple protein model of a four-helix bundle motif on a face-centered cubic lattice has been studied. Total energy of a conformation includes attractive interactions between hydrophobic residues, repulsive interactions between hydrophobic and polar residues, and a potential that favors helical turns. Using replica exchange Monte Carlo simulations we have estimated a set of parameters for which the native structure is a global minimum of conformational energy. Then we have shown that all the above types of interactions are necessary to guarantee the cooperativity of folding transition and to satisfy the thermodynamic hypothesis.     

Excluded volume effects in macromolecular forces and ion-interface interactions

A charged Yukawa liquid confined in a slit nanopore is studied in order to understand excluded volume effects in the interaction force between the pore walls. A previously developed self-consistent scheme [S. Buyukdagli, C. V. Achim, and T. Ala-Nissila, J. Stat. Mech. 2011, P05033] and a new simpler variational procedure that self-consistently couple image forces, surface charge induced electric field, and pore modified core interactions are used to this aim. For neutral pores, it is shown that with increasing pore size, the theory predicts a transition of the interplate pressure from an attractive to a strongly repulsive regime associated with an ionic packing state, an effect observed in previous Monte Carlo simulations for hard core charges. We also establish the mean-field theory of the model and show that for dielectrically homogeneous pores, the mean-field regime of the interaction between the walls corresponds to large pores of size d > 4 ?. The role of the range of core interactions in the ionic rejection and interplate pressure is thoroughly analyzed. We show that the physics of the system can be split into two screening regimes. The ionic packing effect takes place in the regime of moderately screened core interactions characterized with the bare screening parameter of the Yukawa potential b ? 3/l(B), where l(B) is the Bjerrum length. In the second regime of strongly screened core interactions b ? 3/l(B), solvation forces associated with these interactions positively contribute to the ionic rejection driven by electrostatic forces and enhance the magnitude of the attractive pressure. For weakly charged pores without a dielectric discontinuity, core interactions make a net repulsive contribution to the interplate force and also result in oscillatory pressure curves, whereas for intermediate surface charges, these interactions exclusively strengthen the external pressure, thereby reducing the magnitude of the net repulsive interplate force. The pronounced dependence of the interplate pressure and ionic partition coefficients on the magnitude and the range of core interactions indicates excluded volume effects as an important ion specificity and a non-negligible ingredient for the stability of macromolecules in electrolyte solutions.     

Steric stabilization of spherical colloidal particles: Implicit and explicit solvent

We present the results of Monte Carlo simulations and density functional theory treatment of interactions between spherical colloidal brushes both in implicit (good) solvent and in an explicit polymeric solution. Overall, theory is seen to be in good agreement with simulations. We find that interactions between hard-sphere particles grafted with hard-sphere chains are always repulsive in implicit solvent. The range and steepness of the repulsive interaction is sensitive to the grafting density and the length of the grafted chains. When the brushes are immersed in an explicit solvent of hard-sphere chains, a weak mid-range attraction arises, provided the length of the free chains exceeds that of the grafted chains.     

Monte Carlo Simulations of Polymers in Nanoslits

Summary: Monte Carlo computer simulations have been performed for model polymers confined in slits of thickness comparable to the transverse diameter of the chains. The density of polymer within the slits is allowed to vary with the slit thickness in such a way that the content of the slits is always in equilibrium with a large reservoir of bulk polymer. The calculations reveal the presence of polymer‐mediated attractive or repulsive interactions between the slit plates, oscillating with the slit thickness in good agreement with experimental results.

The base cell used in the simulations.     

Collapse transitions of a supersized neutral chain due to irreversibly adsorbed small colloidal particles

 We performed Monte Carlo simulations to study the destabilization processes of large neutral and flexible polymer chains due to irreversibly adsorbed colloidal particles attached to the chains like beads on a necklace. The particles are modeled as charged spherical units which interact with each other via repulsive electrostatic and attractive van der Waals (vdW) potentials. The usual Monte Carlo search procedure is extended and carefully checked to completely sample the chain conformational space and achieve dense conformations in the limit of both strong attractive and repulsive interaction potentials. Configurational properties, such as the radius of gyration, the end-to-end length, and the Kuhn length, are calculated as a function of the intensity of the vdW interactions and ionic strength values. It is observed that chains exhibit a new range of possible conformations compared to the classical random walk and self avoiding walk chains or polyelectrolytes. In the limit of low salt concentration, by gradually increasing vdW interactions, chains undergo a cascade of transitions from extended structures to dumbbells, from dumbbells to pearl necklaces, and from pearl necklaces to collapsed coils. Because of strong competition between the vdW and electrostatic forces, the distance along the chain between the interacting particles, and the sampling limitations, these transitions are found to sample metastable domains and to depend on the initial conformations. To gain insight into the spatial organization of the collapsed conformations, the pair correlation functions of both monomers and particles are calculated. It is shown that collapsed conformations which are the result of strong particle–particle interactions exhibit two distinct parts: a hard core mainly composed of particles and a surrounding polymeric shell composed of loops and tails. Possible effects of such a collapsed transition on the kinetics of flocculation of a mixture containing large flexible chains and small adsorbing colloidal particles are discussed. Received: 26 July 1999 Accepted in revised form: 9 November 1999     

Simulation of phase equilibria in lamellar surfactant systems

The coexistence of two lamellar liquid crystalline phases has been investigated by means of Monte Carlo simulations. The surfaces of the negatively charged bilayers formed by the surfactant molecules are modeled as planar infinite walls with a uniform surface charge density. Water is treated as a dielectric continuum, and only electrostatic interactions are considered. The counterions are mono- and divalent point ions, and their ratio is allowed to vary. Monovalent counterions lead to a repulsive osmotic pressure at all separations, while an attractive region exists when the counterions are divalent. In the latter case, one would expect a phase separation to take place, although it is not observed experimentally due to the limited stability of the lamellar phase at high water content. In a system with mixed counterions, however, the osmotic pressure exhibits a van der Waals loop under such conditions that two phases can coexist. A phase diagram is constructed, and the agreement with experimental data is excellent.     

Effect of discrete macroion charge distributions in solutions of like-charged macroions

The effect of replacing the conventional uniform macroion surface charge density with discrete macroion charge distributions on structural properties of aqueous solutions of like-charged macroions has been investigated by Monte Carlo simulations. Two discrete charge distributions have been considered: point charges localized on the macroion surface and finite-sized charges protruding into the solution. Both discrete charge distributions have been examined with fixed and mobile macroion charges. Different boundary conditions have been applied to examine various properties. With point charges localized on the macroion surface, counterions become stronger accumulated to the macroion and the effect increases with counterion valence. As a consequence, with mono- and divalent counterions the potential of mean force between two macroions becomes less repulsive and with trivalent counterions more attractive. With protruding charges, the excluded volume effect dominates over the increased correlation ability; hence the counterions are less accumulated near the macroions and the potential of mean force between two macroions becomes more repulsive/less attractive.     

Nonequilibrium effects in diffusion of interacting particles on vicinal surfaces

We study the influence of nonequilibrium conditions on the collective diffusion of interacting particles on vicinal surfaces. To this end, we perform Monte Carlo simulations of a lattice-gas model of an ideal stepped surface, where adatoms have nearest-neighbor attractive or repulsive interactions. Applying the Boltzmann-Matano method to spreading density profiles of the adatoms allows the definition of an effective, time-dependent collective diffusion coefficient D(C) (t)(theta) for all coverages theta. In the case of diffusion across the steps and strong binding at lower step edges we observe three stages in the behavior of the corresponding D(xx,C) (t)(theta). At early times when the adatoms have not yet crossed the steps, D(xx,C) (t)(theta) is influenced by the presence of steps only weakly. At intermediate times, where the adatoms have crossed several steps, there are sharp peaks at coverages theta<1L and theta>1-1L, where L is the terrace width. These peaks are due to different rates of relaxation of the density at successive terraces. At late stages of spreading, these peaks vanish and D(xx,C) (t)(theta) crosses over to its equilibrium value, where for strong step edge binding there is a maximum at theta=1L. In the case of diffusion in direction along the steps the nonequilibrium effects in D(yy,C) (t)(theta) are much weaker, and are apparent only when diffusion along ledges is strongly suppressed or enhanced.