Structure of short polymers at interfaces: a combined simulation and theoretical study

The structure of polymers confined between surfaces is studied using computer simulation and a density functional approach. The simple model system considers the polymer molecule as a pearl necklace of freely jointed hard spheres, having attractions among the beads, confined between attractive surfaces. This approach uses the universality of the free-energy functional to obtain the self-consistent field required in the single chain simulation. The second-order direct correlation function for the uniform bulk fluid required as input has been calculated from the reference interaction site model integral equation theory using mean spherical approximation. The theoretical results are shown to compare well with the Monte Carlo simulation results for varying densities, chain lengths, and with different attractive interaction parameters. The simulation results on the conformational properties give important indications regarding the behavior of chains as they approach the surfaces.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

Effect of attractive interactions on the structure of polymer melts confined between surfaces: a density-functional approach

A density-functional theory is presented to study the structure of polymers, having attractive interactions, confined between attractive surfaces. The theory treats the ideal-gas free-energy functional exactly and uses weighted density approximation for the hard-chain contribution to the excess free-energy functional. The bulk interactions of freely jointed hard spheres are obtained from generalized Flory equation of state and the attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions are found to be in quite good agreement with the Monte Carlo simulation results for varying densities, chain lengths, and different interaction potentials. The results confirm important implications of using different approximations for the hard-sphere and attractive interactions.     

Structure and phase behavior of a confined nanodroplet composed of the flexible chain molecules

A polymer density functional theory has been employed for investigating the structure and phase behaviors of the chain polymer, which is modelled as the tangentially connected sphere chain with an attractive interaction, inside the nanosized pores. The excess free energy of the chain polymer has been approximated as the modified fundamental measure-theory for the hard spheres, the Wertheim's first-order perturbation for the chain connectivity, and the mean-field approximation for the van der Waals contribution. For the value of the chemical potential corresponding to a stable liquid phase in the bulk system and a metastable vapor phase, the flexible chain molecules undergo the liquid-vapor transition as the pore size is reduced; the vapor is the stable phase at small volume, whereas the liquid is the stable phase at large volume. The wide liquid-vapor coexistence curve, which explains the wide range of metastable liquid-vapor states, is observed at low temperature. The increase of temperature and decrease of pore size result in a narrowing of liquid-vapor coexistence curves. The increase of chain length leads to a shift of the liquid-vapor coexistence curve towards lower values of chemical potential. The coexistence curves for the confined phase diagram are contained within the corresponding bulk liquid-vapor coexistence curve. The equilibrium capillary phase transition occurs at a higher chemical potential than in the bulk phase.     

Confinement free energy of surfaces bearing end-grafted polymers in the mushroom regime and local measurement of the polymer density

End-tethered polymer chains usually adopt mushroomlike structures on the surface when their density is low. The behaviors of these surface-attached hemicoils are described by existing polymer theory. Dolan and Edwards derived the free energy of a single polymer chain confined between two planar surfaces. Their theory was used to approximate the steric interaction free energy, E, of two identical surfaces bearing polymers in the mushroom regime and to compare with experimental data obtained from surface force measurements. However, because of a mislabeled plot in the original paper, experimental force profiles did not seem to fit the free energy approximation satisfactorily. We have correctly relabeled the involved plot and derived a new simple expression for E. In order to verify this expression, we have performed experiments on PEG45 polymers incorporated in lipid bilayers using a surface force apparatus. The measured force profiles are in perfect agreement with the prediction. We show that such measurements can be used to determine the local density of grafted polymer with good precision.     

Density functional theory of homopolymer mixtures confined in a slit

A density functional theory (DFT) is developed for polymer mixtures with shorted-ranged attractive interparticle interactions confined in a slit. Different weighting functions are used separately for the repulsive part and the attractive part of the excess free energy functional by applying the weighted density approximation. The predicted results by DFT are in good agreement with the corresponding simulation data indicating the reliability of the theory. Furthermore, the center-of-mass profiles and the end-to-end distance distributions are obtained by the single chain simulation; the predictions also agree well with simulation data. The results reveal that both the attraction of the slit wall and the temperature has stronger effect on longer chains than on shorter ones because the intrasegment correlation of chains increases with increasing chain length.     

Ordering of amphiphilic Janus particles at planar walls: a density functional study

We investigate the structure formation of amphiphilic molecules at planar walls using density functional theory. The molecules are modeled as (hard) spheres composed of a hydrophilic and hydrophobic part. The orientation of the resulting Janus particles is described as a vector representing an internal degree of freedom. Our density functional approach involves fundamental measure theory combined with a mean-field approximation for the anisotropic interaction. Considering neutral, hydrophilic, and hydrophobic walls, we study the adsorption of the particles, focusing on the competition between the surface field and the interaction-induced ordering phenomena. Finally, we consider systems confined between two planar walls. It is shown that the anisotropic Janus interaction yields pronounced frustration effects at low temperatures.     

Solvent density inhomogeneities and solvation free energies in supercritical diatomic fluids: a density functional approach

Density functional theory is used to explore the solvation properties of a spherical solute immersed in a supercritical diatomic fluid. The solute is modeled as a hard core Yukawa particle surrounded by a diatomic Lennard-Jones fluid represented by two fused tangent spheres using an interaction site approximation. The authors' approach is particularly suitable for thoroughly exploring the effect of different interaction parameters, such as solute-solvent interaction strength and range, solvent-solvent long-range interactions, and particle size, on the local solvent structure and the solvation free energy under supercritical conditions. Their results indicate that the behavior of the local coordination number in homonuclear diatomic fluids follows trends similar to those reported in previous studies for monatomic fluids. The local density augmentation is particularly sensitive to changes in solute size and is affected to a lesser degree by variations in the solute-solvent interaction strength and range. The associated solvation free energies exhibit a nonmonotonous behavior as a function of density for systems with weak solute-solvent interactions. The authors' results suggest that solute-solvent interaction anisotropies have a major influence on the nature and extent of local solvent density inhomogeneities and on the value of the solvation free energies in supercritical solutions of heteronuclear molecules.     

Density functional theory study on the structure and capillary phase transition of a polymer melt in a slitlike pore: effect of attraction

A density functional theory is proposed to investigate the effects of polymer monomer-monomer and monomer-wall attractions on the density profile, chain configuration, and equilibrium capillary phase transition of a freely jointed multi-Yukawa fluid confined in a slitlike pore. The excess Helmholtz energy functional is constructed by using the modified fundamental measure theory, Wertheim's first-order thermodynamic perturbation theory, and Rosenfeld's perturbative method, in which the bulk radial distribution function and direct correlation function of hard-core multi-Yukawa monomers are obtained from the first-order mean spherical approximation. Comparisons of density profiles and bond orientation correlation functions of inhomogeneous chain fluids predicted from the present theory with the simulation data show that the present theory is very accurate, superior to the previous theory. The present theory predicts that the polymer monomer-monomer attraction lowers the strength of oscillations for density profiles and bond orientation correlation functions and makes the excess adsorption more negative. It is interesting to find that the equilibrium capillary phase transition of the polymeric fluid in the hard slitlike pore occurs at a higher chemical potential than in bulk condition, but as the attraction of the pore wall is increased sufficiently, the chemical potential for equilibrium capillary phase transition becomes lower than that for bulk vapor-liquid equilibrium.     

Surface forces in polymer fluids: a comparison between simulations and density functional theory

A polymer density functional theory is evaluated in terms of its ability to predict interactions between large surfaces in a polymer fluid. Comparisons are made with results from simulations in an expanded isotension ensemble. The variation of the net surface-surface interaction with adsorption strength is examined. Cases where the monomers interact via a pure hard sphere potential are investigated, but we have also studied the effect of attractions between the monomers. In all cases, we obtain an almost quantitative agreement between the simulated results and the predictions from the polymer density functional theory.     

Density functional approach to the adsorption of spherical molecules on a surface modified with attached short chains

A density functional and Monte Carlo simulation study of end-grafted polymers immersed by simple fluids is presented. The polymer molecules are modeled as freely jointed tangent hard spheres with the end segments linked to the surface. The authors analyze an influence of the chain length, the grafting density, and a nature of solvent on the brush structure. Adsorption of hard-sphere mixtures on the modified surface is also discussed. The theory precisely approximates simulation data.     

The solvophobic effect in simple fluid mixtures

The “solvophobic” effect is the tendency of solute particles to cluster as the attractive interaction between solvent particles is strengthened. The potentials of mean force in the hypernetted chain approximation have been obtained for a simple fluid mixture consisting of two hard spheres immersed in a hard core/Yukawa tail solvent. The results clearly exhibit features attributable to the “solvophobic” effect.     

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.     

Exact density functional theory for ideal polymer fluids with nearest neighbor bonding constraints

We present a new density functional theory of ideal polymer fluids, assuming nearest-neighbor bonding constraints. The free energy functional is expressed in terms of end site densities of chain segments and thus has a simpler mathematical structure than previously used expressions using multipoint distributions. This work is based on a formalism proposed by Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005)]. Those authors obtain an approximate free energy functional for ideal polymers in terms of monomer site densities. Calculations on both repulsive and attractive surfaces show that their theory is reasonably accurate in some cases, but does differ significantly from the exact result for longer polymers with attractive surfaces. We suggest that segment end site densities, rather than monomer site densities, are the preferred choice of "site functions" for expressing the free energy functional of polymer fluids. We illustrate the application of our theory to derive an expression for the free energy of an ideal fluid of infinitely long polymers.     

Effect of surface grafted polymers on the adsorption of different model proteins

Adsorption of a model protein to a surface with end-grafted polymers was studied by Monte Carlo simulations. In the model the effect on protein adsorption in the presence of end-grafted polymers was evaluated by calculating the change in free energy between an end-grafted surface and a surface without polymers. The change in free energy was calculated using statistical mechanical perturbation theory. Apart from ordinary athermal polymer-polymer and protein-polymer interactions we also study a broad selection of systems by varying the interaction between proteins and polymers and effective polymer-solvent interactions. The interactions between the molecules span an interval from -0.5 to +0.5 kT. Consequently, general features of protein adsorption to end-grafted surfaces is investigated by systematically changing properties like hydrophilicity/hydrophobicity of the polymer, protein and surface as well as grafting density, degree of polymerization and protein size. Increasing grafting density as well as degree of polymerization decreases the adsorption of protein except in systems with attractive polymer-protein interactions, where adsorption increases with increasing chain length and higher grafting density. At a critical polymer-protein interaction neither chain length nor grafting density affects the free energy of adsorption. Hydrophilic polymers were found to prevent adsorption better than hydrophobic polymers. Very small particles with radii comparable to the size of a polymer segment were, however, better excluded from the surface when using hydrophobic than hydrophilic polymers. For systems with attractive polymer-protein interaction not only the volume of the protein was shown to be of importance but also the size of the exposed surface.     

Hybrid density functional theory for homopolymer mixtures confined in a selective nanoslit

By integrating polymer density function theory (DFT) and single-chain molecular simulation, a hybrid DFT is developed for homopolymer mixtures confined in a selective nanoslit. Two weighting functions are adopted separately in the polymer DFT for repulsive and attractive contributions to the excess free energy functional. The theoretical results agree well with simulation data for the density profiles, configurations (tail, loop and train), adsorption amounts, layer thicknesses, and partition coefficients. The polymer-slit interaction is found to have a large effect on the density profiles and partition coefficients but is found to have a small effect on the average sizes and percentages of the configurations. Nearly half of the polymer segments form tails, and the other half form trains. In addition, bridges are observed to form for sufficiently long polymer chains. As the length difference between two polymers increases, the effect of chain connectivity becomes increasingly important.     

Adsorption of short heteropolymers in slitlike pores

Adsorption of short linear heteropolymers in slitlike pores is studied using the density functional theory and Monte Carlo simulations. The molecules are assumed to be freely jointed tangent hard spheres. The segments have different affinity with regard to the walls. Each molecule contains one surface-binding segment that interacts with the walls via Lennard-Jones (3,9) potential and a number of segments interacting with surfaces via the hard-wall potentials. A position of the surface-binding segment in the chain can be arbitrarily chosen. We have studied the influence of the pore width, the chain length and the chemical structure of molecules on adsorption and the microscopic structure of the confined fluid. The theoretical predictions are compared with Monte Carlo simulations carried out for different 'isomeric' pentamers.     

改进的聚合物溶液理论

在高分子溶液理论中引入Gibbs分布 ,用统计物理学方法重新推导出了聚合物溶液的热力学公式 .将高分子溶液的自由能和熵分三部分进行了计算 ,无热平动部分 ,无热构象部分和构象有热部分 .无热平动自由能和无热构象自由能分别等于Flory Huggins混合自由能公式的前两项 ,构象有热部分引入了Gibbs分布 ,考虑了链段 溶剂分子相互作用对高分子构象的影响 .在分子间的相互作用足够小时 ,又回到了FH公式     

Solvation forces in ionic and neutral liquids: A density functional approach

The solvation forces between two planar charged surfaces in ionic solutions, corresponding to charged and neutral hard spheres representing the ions and the solvent, respectively, are studied here using a weighted density functional theory for inhomogeneous Coulomb systems developed by us recently. The hard sphere contributions to the one-particle correlation function are evaluated nonperturbatively using a position-dependent effective density, while the electrical contributions are obtained through a perturbative expansion around this weighted density. The calculated results on the solvation forces between two charged hard walls compare well with available simulation results for ionic systems. For a neutral system, the present results show good agreement with the experimentally observed oscillating forces for two mica surfaces in octamethylcyclotetrasiloxane. The present approach thus provides a direct route to the calculation of interaction energies between colloidal particles.     

Ionic solvation at the solid/electrolyte interface: A statistical-mechanical approach

We have studied studied the influence of the size of ions on their adsorbability at a solid surface in the presence of a molecular solvent. Ions and molecules are represented respectively by charged hard spheres and dipolar hard spheres and the surface is just a neutral hard wall. We have found that the electrostatic interaction between ions and molecules can induce the exclusion of small ions from the surface. A pure MSA (mean spherical approximation) calculation would not give any effect of the solvation on the ionic density profile. The present calculation is limited to the case of infinite ionic dilution.