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.     

Friction and normal interaction forces between irreversibly attached weakly charged polymer brushes

Polyelectrolyte brushes were built on mica by anchoring polystyrene-poly(acrylic acid) (PS-b-PAA) diblock copolymers at a controlled surface density in a polystyrene monolayer covalently attached to OH-activated mica surfaces. Compared to physisorbed polymer brushes, these irreversibly attached charged brushes allow the polymer grafting density to remain constant upon changes in environmental conditions (e.g., pH, salt concentration, compression, and shear). The normal interaction and friction forces as a function of surface separation distance and at different concentrations of added salt (NaCl) were investigated using a surface forces apparatus. The interaction force profiles were completely reversible both on loading and receding and were purely repulsive. For a constant polymer grafting density, the influence of the polyelectrolyte charges and the Debye screening effect on the overall interaction forces was investigated. The experimental interaction force profiles agree very well with scaling models developed for neutral and charged polymer brushes. The variation of the friction force between two PAA brushes in motion with respect to each other as a function of surface separation distance appeared to be similar to that observed with neutral brushes. This similarity suggests that the increase in friction is associated with an increase in mutual interpenetration upon compression as observed with neutral polymers. The effect of the PAA charges and added ions was more significant on the repulsive normal forces than on the friction forces. The reversible characteristics of the normal force profiles and friction measurements confirmed the strong attachment of the PAA brushes to the mica substrate. High friction coefficients (ca 0.3) were measured at relatively high pressures (40 atm) with no surface damage or polymer removal.     

Density functional theory for rod-coil polymers with different size segments

A polymer density functional theory (PDFT) for rod-coil copolymers with different size segments is proposed, in which the PDFT approach combines a modified fundamental measure theory for the excluded-volume effects, Wertheim's first-order thermodynamics perturbation theory for the chain connectivity and the mean field approximation for van der Waals attraction. First, for testing the PDFT derived, we compare the density profiles from present theory to simulation data, and find that the present theory successfully reproduces the simulation data. Therefore, we use the PDFT to further investigate the local density and solvation forces of rod-coils with different size (A(5)D(3)) and the same size (A(5)B(3)) segments. Results indicate that the excluded volume effect from the coil part determines the solvation force profiles of two rod-coil brushes at strong surface energy. In addition, owing to the vacuum effect, the weak attraction around the classical contact of the rod-coil brushes is also observed. In short, the present theory can be easily applied to the other architecture polymers containing different size segments. It is expected that the calculation results in this work could provide useful reference to select the rod-coils as stabilizer for the protection of surfaces or the colloidal stabilization.     

Optimization of functionalized polymer layers for specific targeting of mobile receptors on cell surfaces

The reversible binding between a planar polymer layer functionalized by ligands and a planar cell surface containing different densities of mobile receptors has been studied by Monte Carlo simulations. Using the acceptance-ratio method, the distance-dependent profiles for the average number of ligands bound to receptors, the total free energy for the polymer layer-cell surface interaction and the interaction force were obtained. Four main design parameters for the polymer layer were considered: the degree of functionalization, chain degree of polymerization, polymer grafting density and the binding energy for the ligand-receptor interaction. We found that an increase in the degree of functionalization or in the absolute energy of ligand-receptor binding results in a larger number of ligands bound to the receptors, lower free energy, and stronger attractive force. Polymer layers composed of shorter chains were found to exhibit a deeper and narrower free energy profile and a larger attractive force, while longer tethers can interact with the cell surface at a larger and broader range of separation distances, in agreement with experimental observations. Our simulation results show that the increase in polymer grafting density from the mushroom to brush regime enhances the ligand availability and results in a stronger attractive force, increases the maximum binding distance, but exhibits a shallower free energy minimum due to the smaller tolerance to compression for polymer layers with high grafting density. We used two measures of the polymer layer binding affinity to the cell surface: the free energy minimum, related to the equilibrium binding constant and the fraction of bound ligands. We found that the polymer layers with a smaller chain length and grafting density, larger degree of functionalization, and larger absolute binding energy exhibit both a larger equilibrium binding constant to the cell surface and a larger average number of bound ligands, except for high binding energies when the maximum level of binding is reached independently of polymer length and grafting density. We showed that high binding specificity can be achieved by the polymer layers with intermediate ligand-receptor binding energies or an intermediate number of ligands, as a larger binding energy or number of ligands ensures a high binding affinity but lacks specificity while a smaller binding energy or number of ligands provides inadequate affinity. We found that the results for polymer layers with different properties follow a similar pattern when both high binding affinity to cells with high receptor density and high binding specificity are considered. As a result, the optimal design of the polymer layers can be achieved by using several different strategies, which are discussed.     

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.     

Grafting acrylic polymers from flat nickel and copper surfaces by surface-initiated atom transfer radical polymerization

Acrylic polymers, including poly(methyl methacrylate), poly(2,2,2-trifluoroethyl methacrylate), poly( N,N'-dimethyaminoethyl methacrylate), and poly(2-hydroxyethyl methacrylate) were grafted from flat nickel and copper surfaces through surface-initiated atom transfer radical polymerization (ATRP). For the nickel system, there was a linear relationship between polymer layer thickness and monomer conversion or molecular weight of "free" polymers. The thickness of the polymer brush films was greater than 80 nm after 6 h of reaction time. The grafting density was estimated to be 0.40 chains/nm2. The "living" chain ends of grafted polymers were still active and initiated the growth of a second block of polymer. Block copolymer brushes with different block sequences were successfully prepared. The experimental surface chemical compositions as measured by X-ray photoelectron spectroscopy agreed very well with their theoretical values. Water contact angle measurements further confirmed the successful grafting of polymers from nickel and copper surfaces. The surface morphologies of all samples were studied by atomic force microscopy. This study provided a novel approach to prepare stable functional polymer coatings on reactive metal surfaces.     

Combining electrochemistry and direct force measurements: from the control of surface properties towards applications

Direct force measurements contributed in the last years much to our understanding of the diffuse double layer of charged interfaces in electrolyte solutions. Such measurements have been performed with the atomic force microscope or the surface force apparatus. This review gives an overview over the recent studies based on force measurements with electrode surfaces. Not only bare metal electrodes but also electrodes modified by different organic layers, including electroactive films, have been studied by these techniques. Direct force measurements indicate that further effects besides classical Gouy–Chapman–Stern theory have to be taken into consideration in order to describe the force profiles. In addition to the long-range forces also the adhesion between surfaces can be tuned by potentiostatic control. New single-molecule techniques based on the atomic force microscope allow to probe the extension of polymer strands or their desorption from solid interfaces. In combination with electrochemistry, it became now possible to tune the desorption behavior of polymer strands or to measure the electromechanical coupling of motors from single strands of electroactive polymers.     

Use of numerical methods in the prediction of polymer adsorption

Average polymer segment densities and thermodynamic properties of polymer adsorbed at liquid–solid interfaces were computed by extension of the polymer adsorption theory of Forsman and Hughes. Expressions were derived for the total free energy of adsorbed polymer chains by using the Flory-Huggins theory to represent free energy of mixing. A square-well potential was used to represent segment–surface interaction, and configurational entropy was calculated from the probability density function for the radius of gyration of random-flight chains. For each specified amount of surface coverage the free energy of the adsorbed polymer was minimized by varying the density of segments normal to the adsorbing surface and using a modified gradient search algorithm on a digital computer. Two different segment densities were considered, and they both gave qualitatively the same results. The two densities were (1) the sum of two Gaussian distributions and (2) a two-step density distribution. Isotherms were then calculated by equating the partial molal free energy of polymer at the surface to that of polymer in bulk solution for each specified amount of surface coverage. The results showed that for the initial region of the isotherms the distribution of polymer segments normal to the surface consisted of a high-density layer adjacent to the surface and a low-density “tail” extending far out into the solution. At higher amounts of adsorbed polymer, i.e., in the general concentration range of the pseudo-plateau, the tail of the polymer density distribution was predicted to thicken, and a single Gaussian distribution best described the segment density. Predicted adsorptions were in good agreement with reported experimental values.     

Use of numerical methods in the prediction of polymer adsorption

Average polymer segment densities and thermodynamic properties of polymer adsorbed at liquid–solid interfaces were computed by extension of the polymer adsorption theory of Forsman and Hughes. Expressions were derived for the total free energy of adsorbed polymer chains by using the Flory-Huggins theory to represent free energy of mixing. A square-well potential was used to represent segment–surface interaction, and configurational entropy was calculated from the probability density function for the radius of gyration of random-flight chains. For each specified amount of surface coverage the free energy of the adsorbed polymer was minimized by varying the density of segments normal to the adsorbing surface and using a modified gradient search algorithm on a digital computer. Two different segment densities were considered, and they both gave qualitatively the same results. The two densities were (1) the sum of two Gaussian distributions and (2) a two-step density distribution. Isotherms were then calculated by equating the partial molal free energy of polymer at the surface to that of polymer in bulk solution for each specified amount of surface coverage. The results showed that for the initial region of the isotherms the distribution of polymer segments normal to the surface consisted of a high-density layer adjacent to the surface and a low-density “tail” extending far out into the solution. At higher amounts of adsorbed polymer, i.e., in the general concentration range of the pseudo-plateau, the tail of the polymer density distribution was predicted to thicken, and a single Gaussian distribution best described the segment density. Predicted adsorptions were in good agreement with reported experimental values.     

Combined depletion and electrostatic forces in polymer-induced membrane adhesion: a theoretical model

We develop a semi-quantitative analytical theory to describe adhesion between two identical planar charged surfaces embedded in a polymer-containing electrolyte solution. Polymer chains are uncharged and differ from the solvent by their lower dielectric permittivity. The solution mimics physiological fluids: It contains 0.1 M of monovalent ions and a small number of divalent cations that form tight bonds with the headgroups of charged lipids. The components have heterogeneous spatial distributions. The model was derived self-consistently by combining: (a) a Poisson-Boltzmann like equation for the charge densities, (b) a continuum mean-field theory for the polymer profile, (c) a solvation energy forcing the ions toward the polymer-poor regions, and (d) surface interactions of polymers and electrolytes. We validated the theory via extensive coarse-grained Molecular Dynamics (MD) simulations. The results confirm our analytical model and reveal interesting details not detected by the theory. At high surface charges, polymer chains are mainly excluded from the gap region, while the concentration of ions increases. The model shows a strong coupling between osmotic forces, surface potential and salting-out effects of the slightly polar polymer chains. It highlights some of the key differences in the behaviour of monomeric and polymeric mixed solvents and their responses to Coulomb interactions. Our main findings are: (a) the onset of long-ranged ion-induced polymer depletion force that increases with surface charge density and (b) a polymer-modified repulsive Coulomb force that increases with surface charge density. Overall, the system exhibits homeostatic behaviour, resulting in robustness against variations in the amount of charges. Applications and extensions of the model are briefly discussed.     

Weak polyelectrolytes tethered to surfaces: Effect of geometry,acid–base equilibrium and electrical permittivity

The structural and thermodynamical properties of weak polyelectrolytes end-tethered to surfaces of arbitrary geometry are studied using a molecular theory. The theory is based on writing down the free energy functional of the system including all the basic interactions and the explicit acid–base equilibrium for the chargeable groups of the polymer. The theory explicitly includes the size, shape, conformations, and charge distribution of all the molecular species. The electrostatic interactions include a density-dependent dielectric function, modeled with the Maxwell–Garnett mixing formula, to account for the composition-dependent permittivity. The minimization of the free energy leads to the distribution of all molecular species and their dependence on bulk pH and salt concentration. We apply the theory to polymer chains end-tethered to planar, cylindrical, and spherical surfaces. The radius of the curved surfaces is small to enhance the curvature effect. We find that when the grafting surfaces are uncharged, the approximation of a constant dielectric function works very well for both structural and thermodynamic properties. The structure of weak polyelectrolytes tethered on cylindrical and spherical surfaces is different from that of polymers tethered on planar surfaces due to the available volume as a function of the distance from the surface. Specifically, the degree of dissociation increases with increasing curvature of the surface. This is a manifestation of the coupling between the local density of protons, counterions, and polymer segments. The results can be interpreted in terms of the local Le Chatelier principle for the acid–base equilibrium, with proper account of the three local contributions: counterions, protons, and chargeable groups. We find that one can achieve local changes of pH between one to two units within 1–2 nm. The thickness of the tethered layers as a function of bulk pH shows a large increase when the pH is equal to the bulk pK. However, the variation with salt concentration is different for the different geometries. The largest swelling is found for cylindrical surfaces. The predictions from scaling theories of a maximum in the thickness of the film as a function of salt concentration is found for planar films, but not for curved surfaces. Finally, the interactions between cylinders with tethered polyelectrolytes is very different from the equivalent planar surfaces. These results are important for the interpretation of force measurements with nanoscale AFM tips. The implications of the results for the rational design of responsive tethered polymer layers is discussed together with the limitations of the theoretical approach. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2638–2662, 2006     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. 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 on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.     

A Constant Potential Molecular Dynamics Simulation Study of the Atomic‐Scale Structure of Water Surfaces Near Electrodes

Novel and technologically important processes and phenomena arise at water surfaces in the presence of electric fields. However, experimental measurements on water surfaces are challenging, and the results are scarce and inconclusive. In this work, the constant potential molecular dynamics method, in which the electrode charges are allowed to fluctuate to keep the electric potential fixed, was implemented in the study of a near‐electrode water surface systems. This simulation system was set up with a vapor/liquid‐water/vapor slab and two electrodes under different sets of applied electrostatic potential, yielding i) a detailed characterization of the external E‐field dependent electrostatic potential/density/dipole moment density profiles, and ii) the relationship between the water surface width and the applied electrode voltage differences which has been rarely reported. The adjustments in the number density profiles in the vicinity of water surfaces due to external E‐fields were observed, while the capillary interfacial widths for the surfaces near both cathode and anode were found with different increment rates under increasing E‐fields. By examining dipole density profiles across the water surfaces, we found that external E‐field induced polarization occurs in both bulk and surface regimes, yet the surface polarization densities vary asymmetrically with respect to the increasing E‐fields. Detailed discussions were carried out to explain the correlation between water surface tension and surface widths, as well as the interplay between the surface polarization densities and the hydrogen bond network structure. We conclude that the mechanical and structural properties of the water surfaces could be tuned by both magnitude and direction of the strong external E‐fields. We also recognize that more surface properties with application value, such as dielectric permittivity tensor or surface potential, could also be regulated by the external E‐fields.     

Polymer depletion profiles around nonspherical colloidal particles

We study the effect of chain self-avoidance on the polymer density profiles that are induced by a single colloidal particle of nonspherical shape such as an ellipsoid, a dumbbell, or a lens in a solution of nonadsorbing polymers. For colloid sizes sigma much smaller than the size R(x) of the polymers, we observe a pronounced difference between ideal and self-avoiding chains. In the case of ideal polymers, the surfaces of constant density always have the same character as the surface of the particle, e.g., are oblate for an oblate ellipsoid. In the self-avoiding case, however, the character changes with increasing distance r from the particle, and an oblate particle induces prolate surfaces of constant density if sigma相似文献   

Static and dynamic properties of the interface between a polymer brush and a melt of identical chains

Molecular-dynamics simulations of a short-chain polymer melt between two brush-covered surfaces under shear have been performed. The end-grafted polymers which constitute the brush have the same chemical properties as the free chains in the melt and provide a soft deformable substrate. Polymer chains are described by a coarse-grained bead-spring model, which includes excluded volume and backbone connectivity of the chains. The grafting density of the brush layer offers a way of controlling the behavior of the surface without altering the molecular interactions. We perform equilibrium and nonequilibrium molecular-dynamics simulations at constant temperature and volume using the dissipative particle dynamics thermostat. The equilibrium density profiles and the behavior under shear are studied as well as the interdigitation of the melt into the brush, the orientation on different length scales (bond vectors, radius of gyration, and end-to-end vector) of free and grafted chains, and velocity profiles. The obtained boundary conditions and slip length show a rich behavior as a function of grafting density and shear velocity.     

Ligand-receptor interactions between surfaces: the role of binary polymer spacers

The interactions between a receptor-modified planar surface and a surface grafted with a bimodal polymer layer, where one of the polymer species is ligand functionalized, are studied using a molecular theory. The effects of changing the binding energy of the ligand-receptor pair, the polymer surface coverage, the composition, and molecular weight of both the unfunctionalized and ligand functionalized polymers on the interactions between the surfaces are investigated. Our findings show that bridging exists between the surfaces including when the molecular weight of the ligand-bearing polymer is smaller than that of the unfunctionalized polymer, even though the ligand is initially buried within the polymer layer. The distance at which the surfaces bind depends only on the molecular weight of the ligand-modified polymer, while the strength of the interaction at a given surface separation can be tuned by changing the molecular weight of the polymers, the total polymer surface coverage, and the fraction of ligated polymers. The composition of the bimodal layer alters the structure of the polymer layer, thereby influencing the strength of the steric repulsions between the surfaces. Our theoretical results show good agreement with experimental data. The present theoretical study can be used as guidelines for the design of surfaces with tailored abilities for tunning the binding strength and surface-ligand separation distances for polymer-grafted surfaces bearing specific targeting ligands.     

Choose sides: differential polymer adhesion

AFM-based single molecule desorption measurements were performed on surface end-grafted poly(acrylic acid) monolayers as a function of the pH of the aqueous buffer to study the adhesion properties of polymers that bridge two surfaces. These properties were found to depend on the adhesion forces of both surfaces in a differential manner, which is explained with a simple model in analogy to the Bell-Evans formalism used in dynamic force spectroscopy. The measured interaction forces between the poly(acrylic acid) chains and silicon nitride AFM tips depend on the grafting density of the polymer monolayers as well as on the contour length of the polymer chains. This study demonstrates that the stability of polymer bridges is determined by the adhesion strengths on both surfaces, which can be tuned by using pH-dependent polyelectrolyte monolayers.     

Simulating equilibrium surface forces in polymer solutions using a canonical grid method

A new simulation method for nonuniform polymer solutions between planar surfaces at full chemical equilibrium is described. The technique uses a grid of points in a two-dimensional thermodynamic space, labeled by surface area and surface separations. Free energy differences between these points are determined via Bennett's optimized rates method in the canonical ensemble. Subsequently, loci of constant chemical potential are determined within the grid via simple numerical interpolation. In this way, a series of free energy versus separation curves are determined for a number of different chemical potentials. The method is applied to the case of hard sphere polymers between attractive surfaces, and its veracity is confirmed via comparisons with established alternative simulation techniques, namely, the grand canonical ensemble and isotension ensemble methods. The former method is shown to fail when the degree of polymerization is too large. An interesting interplay between repulsive steric interactions and attractive bridging forces occurs as the surface attraction and bulk monomer density are varied. This behavior is further explored using polymer density functional theory, which is shown to be in good agreement with the simulations. Our results are also discussed in light of recent self-consistent field calculations which correct the original deGennes results for infinitely long polymers. In particular, we look at the role of chain ends by investigating the behavior of ring polymers.     

Surface characterization of poly(methyl methacrylate-co-n-butyl acrylate-co-cyclopentylstyryl-polyhedral oligomeric silsesquioxane) by inverse gas chromatography

The surface properties of poly(methyl methacrylate-co-n-butyl acrylate-co-cyclopentylstyryl polyhedral oligomeric silsesquioxane) (poly(MMA-co-BA-co-styryl-POSS)) were studied by means of inverse gas chromatography (IGC) using 10 non-polar and polar solvents as the probes. Thermodynamic parameters of adsorption, e.g., specific retention volume, the dispersive component of the surface free energy, the specific interaction contribution to the free energy of adsorption and the acid/base constants were obtained to investigate the interactions between the surfaces of the copolymers and different solvents. It was found that incorporation of styryl-POSS into polymer resulted in increasing interactions between polymers and solvents, dispersive component of surface free energy of polymer and acidity of the surfaces of the polymers. The more the styryl-POSS were embedded, the stronger the interaction between the polymer surface and solvent, the dispersive component of the surface free energy and the acidity of the polymer surface were.