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. 相似文献
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. 相似文献
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. 相似文献
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. 相似文献
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. 相似文献
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. 相似文献
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. 相似文献
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 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. 相似文献
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. 相似文献
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. 相似文献
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.
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 相似文献
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. 相似文献
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. 相似文献
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. 相似文献