Interactions between nanocolloidal particles in polymer solutions: effect of attractive interactions

We present a density-functional theory study of nanoparticle interactions in a concentrated polymer solution. The polymers are modeled as freely jointed tangent chains; all nonbonded interactions between polymer segments and nanoparticles are described by Lennard-Jones potentials. We test several recently proposed methods of treating attractive interactions within the density-functional theory framework by comparing theoretical results with recent simulation data. We find that the simple van der Waals approach provides the most accurate results for the polymer-mediated potential of mean force between two dilute nanoparticles. We employ this approach to study nanoparticle interactions as a function of nanoparticle-segment interaction strength and polymer solution density and temperature.     

Demixing in athermal mixtures of colloids and excluded-volume polymers from density functional theory

We study the structure and interfacial properties of model athermal mixtures of colloids and excluded volume polymers. The colloid particles are modeled as hard spheres whereas the polymer coils are modeled as chains formed from tangentially bonded hard spheres. Within the framework of the nonlocal density functional theory we study the influence of the chain length on the surface tension and the interfacial width. We find that the interfacial tension of the colloid-interacting polymer mixtures increases with the chain length and is significantly smaller than that of the ideal polymers. For certain parameters we find oscillations on the colloid-rich parts of the density profiles of both colloids and polymers with the oscillation period of the order of the colloid diameter. The interfacial width is few colloid diameters wide and also increases with the chain length. We find the interfacial width for the end segments to be larger than that for the middle segments and this effect is more pronounced for longer chains.     

Two- and three-body interactions among nanoparticles in a polymer melt

We perform direct three-dimensional density functional theory (DFT) calculations of two- and three-body interactions in polymer nanocomposites. The nanoparticles are modeled as hard spheres, immersed in a hard-sphere homopolymer melt of freely jointed chains. The two-particle potential of mean force obtained from the DFT is in near quantitative agreement with the potential of mean force obtained from self-consistent polymer reference interaction site model theory. Three-body interactions among three nanoparticles are found to be significant, such that it is not possible to describe these systems with a polymer-mediated two-body interaction calculated from the potential of mean force.     

Effect of polymer size and chain length on depletion interactions between two colloids

A density functional theory based on the weighted density has been developed to investigate the depletion interactions between two colloids immersed in a bath of the binary polymer mixtures, where the colloids are modeled as hard spheres and the polymers as freely jointed tangent hard-sphere chain mixtures. The theoretical calculations for the depletion forces between two colloids induced by the polymer are in good agreement with the computer simulations. The effects of polymer packing fraction, degree of polymerization, polymer/polymer size ratio, colloid/polymer size ratio on the depletion interactions, and colloid-colloid second virial coefficient B2 due to polymer-mediated interactions have been studied. With increasing the polymer packing fraction, the depletion interaction becomes more long ranged and the attractive interaction near the colloid becomes deeper. The effect of degree polymerization shows that the long chain gives a more stable dispersion for colloids rather than the short chain. The strong effective colloid-colloid attraction appears for the large colloid/polymer and polymer/polymer size ratio. The location of maximum repulsion Rmax is found to appear Rmax approximately sigmac+Rg2 for the low polymer packing fraction and this is shifted to smaller separation Rmax approximately sigmac+sigmap2 with increasing the polymer packing fraction, where sigmap2 and Rg2 are the small-particle diameter and the radius of gyration of the polymer with the small-particle diameter, respectively.     

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.     

Microstructure and depletion forces in polymer-colloid mixtures from an interfacial statistical associating fluid theory

By using a classical density functional theory (interfacial statistical associating fluid theory), we investigate the structure and effective forces in nonadsorbing polymer-colloid mixtures. The theory is tested under a wide range of conditions and performs very well in comparison to simulation data. A comprehensive study is conducted characterizing the role of polymer concentration, particle/polymer-segment size ratio, and polymer chain length on the structure, polymer induced depletion forces, and the colloid-colloid osmotic second virial coefficient. The theory correctly captures a depletion layer on two different length scales, one on the order of the segment diameter (semidilute regime) and the other on the order of the polymer radius of gyration (dilute regime). The particle/polymer-segment size ratio is demonstrated to play a significant role on the polymer structure near the particle surface at low polymer concentrations, but this effect diminishes at higher polymer concentrations. Results for the polymer-mediated mean force between colloidal particles show that increasing the concentration of the polymer solution encourages particle-particle attraction, while decreasing the range of depletion attraction. At intermediate to high concentrations, depletion attraction can be coupled to a midrange repulsion, especially for colloids in solutions of short chains. Colloid-colloid second virial coefficient calculations indicate that the net repulsion between colloids at low polymer densities gives way to net attraction at higher densities, in agreement with available simulation data. Furthermore, the results indicate a higher tendency toward colloidal aggregation for larger colloids in solutions of longer chains.     

Interaction between colloids with grafted diblock polyampholytes

The interaction between composite colloidal particles composed of a spherical core and grafted AB-diblock polyampholytes (diblock copolymers with oppositely charged blocks) are investigated by using a coarse-grained model solved with Monte Carlo simulations. The B block is end-grafted onto the core of the colloid and its linear charge density is varied, whereas the linear charge density of the A block is fixed. The brush structure of a single colloid, the mean force between two colloids, and the structure of solutions of such colloids have been determined for different linear charge densities of the B blocks and block lengths. Many features of the present system are controlled by the charge of the B blocks. In the limit of uncharged B blocks, (i) the grafted chains are stretched and form an extended polyelectrolyte brush, (ii) a strong repulsive force is operating between two colloids, (iii) and the solution is thermodynamic stable and displays strong spatial correlation among the colloids. In the limit where the charges of the two types of blocks exactly compensate each other, (i) the chains are collapsed and form a polyelectrolyte complex surrounding the cores, (ii) an attractive force appears between two colloids, and (iii) strong colloid clustering appears in the solution. These features become more pronounced as the length of the polymer blocks is increased, and a phase instability occurs at sufficiently long chains. A comparison with properties for other related colloidal particles is also provided.     

Sterically stabilized lock and key colloids: a self-consistent field theory study

A self-consistent field theory study of lock and key type interactions between sterically stabilized colloids in polymer solution is performed. Both the key particle and the lock cavity are assumed to have cylindrical shape and their surfaces are uniformly grafted with polymer chains. The lock-key potential of mean force is computed for various model parameters, such as length of free and grafted chains, lock and key size matching, free chain volume fraction, grafting density, and various enthalpic interactions present in the system. The lock-key interaction is found to be highly tunable, which is important in the rapidly developing field of particle self-assembly.     

An atomic force microscopy investigation on self-assembled peptide nucleic acid structures on gold(111) surface

We studied interaction of hydrophilic polymer chain and hydrophilic silica nanoparticles in a dilute aqueous system using an idealized model system comprised of a well characterized polyvinyl alcohol of 100 Å R_g and hard spherical LUDOX® silica of 80 Å radii. Interaction among the polymer chains forming polymer clusters with collective polymer structure factor induced by the polymer-mediated potentials of mean force between the nanoparticles, was observed. However, Gaussian nature of individual polymer chain remains unaltered. The dilute system of polymer with low silica volume fraction has the scattering form which was appropriately modeled as the sum of the individual profiles of spherical silica particles and polymer cluster of interchain packing. With increasing silica volume fraction in the dilute solution, the spatial range parameter between the particles is reduced; hence there is a net increase in the mean potential force and consequently to stronger interaction between the silica and polymer. In the dilute systems of high silica with low polymer volume fraction, the polymer chain apparently attracted closer to the silica and concurrently absorbed to the silica hard surface and their scattering data were excellently fit with a model form factor as comprising of one unit forming the core of the spherical silica particles and the interacting polymer as the corona. This result of severe change in polymer interchain conformation in a dilute system corroborated with reduced polymer viscosity observed.     

Alignment of particles in sheared viscoelastic fluids

We investigate the shear-induced structure formation of colloidal particles dissolved in non-Newtonian fluids by means of computer simulations. The two investigated visco-elastic fluids are a semi-dilute polymer solution and a worm-like micellar solution. Both shear-thinning fluids contain long flexible chains whose entanglements appear and disappear continually as a result of Brownian motion and the applied shear flow. To reach sufficiently large time and length scales in three-dimensional simulations with up to 96 spherical colloids, we employ the responsive particle dynamics simulation method of modeling each chain as a single soft Brownian particle with slowly evolving inter-particle degrees of freedom accounting for the entanglements. Parameters in the model are chosen such that the simulated rheological properties of the fluids, i.e., the storage and loss moduli and the shear viscosities, are in reasonable agreement with experimental values. Spherical colloids dispersed in both quiescent fluids mix homogeneously. Under shear flow, however, the colloids in the micellar solution align to form strings in the flow direction, whereas the colloids in the polymer solution remain randomly distributed. These observations agree with recent experimental studies of colloids in the bulk of these two liquids.     

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.     

Potential of mean force in confined colloids: integral equations with fundamental measure bridge functions

The potential of mean force for uncharged macroparticles suspended in a fluid confined by a wall or a narrow pore is computed for solvent-wall and solvent-macroparticle interactions with attractive forces. Bridge functions taken from Rosenfeld's density-functional theory are used in the reference hypernetted chain closure of the Ornstein-Zernike integral equations. The quality of this closure is assessed by comparison with simulation. As an illustration, the role of solvation forces is investigated. When the "residual" attractive tails are given a range appropriate to "hard sphere-like" colloids, the unexpected role of solvation forces previously observed in bulk colloids is confirmed in the confinement situation.     

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.     

环状高分子对胶体稳定性影响的理论研究

将改进的基本度量理论与热力学微扰理论相结合,提出了胶体/高分子系统排空相互作用的IRDFT理论。该理论解释了在实际胶体系统中起重要作用的排斥体积效应和高分子链内相关性的竞争机制。应用该IRDFT,分别以线性链状和环状结构的高分子及其单体为排空元,计算了不同排空元条件下的胶体间排空相互作用。研究表明:对于絮凝破坏,自由环状高分子具有更大的优势。     

Theoretical aspects and computer simulations of flexible charged oligomers in salt-free solutions

The structural and thermodynamic properties of a model solution containing flexible charged oligomers and an equivalent number of counterions were studied by means of the canonical Monte Carlo simulation and integral equation theory. The oligomers were represented as freely jointed chains of charged hard spheres. In accordance with the primitive model of electrolyte solutions, the counterions were modeled as charged hard spheres and the solvent as a dielectric continuum. Simulations were performed for a set of model parameters, independently varying the chain length and concentration of the oligomers. Structural properties in the form of pair distribution functions were calculated as functions of model parameters. In addition, thermodynamic properties such as the excess energy of solution and the excess chemical potential of counterions were obtained. These properties were correlated with the conformational averages of oligomers as reflected in the end-to-end distances and radii of gyration obtained from the simulations. The relation with the experimental data for heats of dilution and for the activity coefficient is discussed. Finally, theories based on Wertheim's integral equation approach (product reactant Ornstein-Zernike approach) [J. Stat. Phys. 42, 477 (1986)] in the so-called polymer mean spherical and polymer hypernetted chain approximations were tested against the new and existing computer simulations. For the values of parameters examined in this study, the integral equation theory yields semiquantitative agreement with computer simulations.     

Modeling of drying in films of colloidal particles

The process of film formation on a solid substrate from polymer colloid dispersion during solvent evaporation has been investigated by means of the Monte Carlo simulation method. Colloid particles are modeled as hard spheres. Time evolution of the colloid density distribution and coverage of the solid substrate are studied. Both density and structure of colloid film is shown to depend strongly on the evaporation rate. At a low evaporation rate, the coexistence of hexagonal and tetragonal domains of dried colloid monolayer has been observed. The results of monolayer structure are in good agreement with the confocal scanning laser microscopy observations of Dullens et al. (2004).     

Interactions between spherical colloids mediated by a liquid crystal: a molecular simulation and mesoscale study

Monte Carlo simulations and dynamic field theory (DyFT) are used to study the interactions between dilute spherical particles, dispersed in nematic and isotropic phases of a liquid crystal. A recently developed simulation method (expanded ensemble density of states) was used to determine the potential of mean force (PMF) between the two spheres as a function of their separation and size. The PMF was also calculated by a dynamic field theory that describes the evolution of the local tensor order parameter. Both methods reveal an overall attraction between the colloids in the nematic phase; in the isotropic phase, the overall attraction between the colloids is much weaker, whereas the repulsion at short range is stronger. In addition, both methods predict a new topology of the disclination lines, which arises when the particles approach each other. The theory is found to describe the results of simulations remarkably well, down to length scales comparable to the size of the molecules. At separations corresponding to the width of individual molecular layers on the particles' surface, the two methods yield different defect structures. We attribute this difference to the neglect of density inhomogeneities in the DyFT. We also investigate the effects of the size of spherical colloids on their interactions.     

Phase behavior of polymer/nanoparticle blends near a substrate

We use the recent fluids density functional theory of Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005); J. Chem. Phys. 122, 094506 (2005)] to investigate the phase behavior of athermal polymer/nanoparticle blends near a substrate. The blends are modeled as a mixture of hard spheres and freely jointed hard chains, near a hard wall. There is a first order phase transition present in these blends in which the nanoparticles expel the polymer from the surface to form a monolayer at a certain nanoparticle concentration. The nanoparticle transition density depends on the length of the polymer, the nanoparticle diameter, and the overall bulk density of the system. The phase transition is due to both packing entropy effects related to size asymmetry between the components and to the polymer configurational entropy, justifying the so-called "entropic push" observed in experiments. In addition, a layered state is found at higher densities which resembles that in colloidal crystals, in which the polymer and nanoparticles form alternating discrete layers. We show that this laminar state has nearly the same free energy as the homogeneously mixed fluid in the bulk and is nucleated by the surface.     

Coarse graining of star-polymer--colloid nanocomposites

We consider mixtures of self-avoiding multiarm star polymers with hard colloids that are smaller than the star polymer size. By employing computer simulations, and by extending previous theoretical approaches, developed for the opposite limit of small star polymers [A. Jusufi et al., J. Phys.: Condens. Matter 13, 6177 (2001)], we coarse-grain the mixture by deriving an effective cross-interaction between the unlike species. The excellent agreement between theory and simulation for all size ratios examined demonstrates that the theoretical approaches developed for the colloidal limit can be successfully modified to maintain their validity also for the present case of the protein limit, in contrast to the situation for mixtures of colloids and linear polymers. We further analyze, on the basis of the derived interactions, the non-additivity parameter of the mixture as a function of size ratio and star functionality and delineate the regions in which we expect mixing as opposed to demixing behavior. Our results are relevant for the study of star-colloid nanocomposites and pave the way for further investigations of the structure and thermodynamics of the same.