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.     

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.     

Interactions between colloidal particles in polymer solutions: A density functional theory study

We present a density functional theory study of colloidal interactions in a concentrated polymer solution. The colloids are modeled as hard spheres and polymers are modeled as freely jointed tangent hard sphere chains. Our theoretical results for the polymer-mediated mean force between two dilute colloids are compared with recent simulation data for this model. Theory is shown to be in good agreement with simulation. We compute the colloid-colloid potential of mean force and the second virial coefficient, and analyze the behavior of these quantities as a function of the polymer solution density, the polymer chain length, and the colloid/polymer bead size ratio.     

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.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

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.     

Steric and bridging interactions between two plates induced by grafted polyelectrolytes

If colloidal particles are grafted with a polymer, then the grafted chains can provide steric repulsion between them. If some of the grafted polymer chains are also adsorbed to a second particle, then a bridging force is generated as well. For uncharged plates and polymer, the following contributions to the free energy of the system have been taken into account in the calculation of the interaction force: (i) the Flory-Huggins expression for the mixing free energy of the grafted chains with the liquid; (ii) the entropy loss due to the connectivity of the polymeric segments; (iii) the van der Waals interactions between the segments and the plates; and (iv) the free energy of adsorption of the polymer segments of the grafted chains on the other plate. For charged plates, the electrostatic free energy as well as the free energy of the electrolyte are included in the total free energy of the system. By minimizing the free energy with respect to the segment concentration and, when it is the case, with respect to the electrical potential, equations for the segment number density distribution and for the electrical potential are obtained, on the basis of which the interactions between two plates grafted with polymer chains that can be also adsorbed on the other plate were calculated. The interaction thus obtained includes steric and bridging forces.     

Probing chain ends in adsorbed polymer layers

We study theoretically the pull-out of polymer chains from an adsorbed polymer layer by sticking of the chain ends on an opposing surface using scaling arguments and a mean field theory. When only one chain is pulled out from the layer, we extend previous results obtained for a single adsorbed chain and calculate the force necessary to extract the chain from the layer. We then discuss end adsorption from an adsorbed layer of polymers bearing specific end groups onto a second surface. Two bridging regimes are predicted: a diffuse layer regime at weak separations (or/and weak interaction) and a large separation strong interaction regime where the bridges stretch into a brush like structure. Bridging fractions and force profiles are displayed that could be compared to atomic force microscope or surface force apparatus experiments.     

Lennard-Jones parameters for the combined QM/MM method using the B3LYP/6-31G*/AMBER potential

A combined DFT quantum mechanical and AMBER molecular mechanical potential (QM/MM) is presented for use in molecular modeling and molecular simulations of large biological systems. In our approach we evaluate Lennard-Jones parameters describing the interaction between the quantum mechanical (QM) part of a system, which is described at the B3LYP/6-31+G* level of theory, and the molecular mechanical (MM) part of the system, described by the AMBER force field. The Lennard-Jones parameters for this potential are obtained by calculating hydrogen bond energies and hydrogen bond geometries for a large set of bimolecular systems, in which one hydrogen bond monomer is described quantum mechanically and the other is treated molecular mechanically. We have investigated more than 100 different bimolecular systems, finding very good agreement between hydrogen bond energies and geometries obtained from the combined QM/MM calculations and results obtained at the QM level of theory, especially with respect to geometry. Therefore, based on the Lennard-Jones parameters obtained in our study, we anticipate that the B3LYP/6-31+G*/AMBER potential will be a precise tool to explore intermolecular interactions inside a protein environment.     

Tethered polymer layers: phase transitions and reduction of protein adsorption

The structural and thermodynamic properties of tethered polymer layers formed by spreading diblock copolymers at a solid surface or at a fluid‐fluid interface are studied using a molecular mean‐field theory. The role of the anchoring block in determining the properties of the tethered polymer layer is studied in detail. It is found that both the anchoring and the tethered blocks are very important in determining the phase behavior of the polymer layer. The structures of the coexisting phases, the phase boundaries and the stability of the layer are found to depend on the ratio of molecular weight between the two blocks, the polymer‐interface (surface) interactions and the strength of the interactions between the two blocks. The different phase transitions found are related to experimental observations. The properties of the polymer layers at coexistence reflect the block that is the dominant driving force for phase separation. The ability of the tethered polymer layers, under different conditions, to control protein adsorption to surfaces is also studied. It is found that the most important factors determining the ability of a polymer layer to reduce the equilibrium amount of proteins adsorbed to a surface are the surface coverage of polymer and the surface‐polymer interactions. The polymer chain length plays only a secondary role. For the kinetic control, however, it is found that the potential of mean‐force, and thus the early stages of adsorption, depends strongly on polymer molecular weight. Further, it is found that the molecular factors determining the ability of the tethered polymer layer to reduce the equilibrium amount of protein adsorption are different than those that control the kinetic behavior. Comparisons with experimental observations are presented. The predictions of the theory are in very good agreement with the measured adsorption isotherms. Guidelines for building optimal surface protection for protein adsorption, both kinetic and thermodynamic, are discussed.     

Interactions between polymer brush-coated spherical nanoparticles: the good solvent case

The interaction between two spherical polymer brushes is studied by molecular dynamics simulation varying both the radius of the spherical particles and their distance, as well as the grafting density and the chain length of the end-grafted flexible polymer chains. A coarse-grained bead-spring model is used to describe the macromolecules, and purely repulsive monomer-monomer interactions are taken throughout, restricting the study to the good solvent limit. Both the potential of mean force between the particles as a function of their distance is computed, for various choices of the parameters mentioned above, and the structural characteristics are discussed (density profiles, average end-to-end distance of the grafted chains, etc.). When the nanoparticles approach very closely, some chains need to be squeezed out into the tangent plane in between the particles, causing a very steep rise of the repulsive interaction energy between the particles. We consider as a complementary method the density functional theory approach. We find that the quantitative accuracy of the density functional theory is limited to large nanoparticle separation and short chain length. A brief comparison to Flory theory and related work on other models also is presented.     

van der Waals interaction of simple, parallel polymers

We study the mutual interactions of simple parallel polymers within the framework of density-functional theory (DFT). As the conventional implementations of DFT do not treat the long-range dispersion [van der Waals (vdW)] interactions, we develop a systematic correction scheme for the nonlocal energy contribution of the polymer interaction at the intermediate to the asymptotic separations. We primarily focus on the three polymers, polyethylene, isotactic polypropylene, and isotactic polyvinylchloride, but the scheme presented applies also more generally to other simple polymers. From first-principle calculations we extract the geometrical and electronic structures of the polymers and the local part of their interaction energy, as well as the static electric response. The dynamic electrodynamic response is modeled on the basis of these static calculations, from which the nonlocal vdW interaction of the polymers is extracted.     

Density functional theory for the recognition of polymer at nanopatterned surface

The recognition of homopolymer at nanopatterned surface has been investigated by density functional theory (DFT). Chain conformation and pattern transfer parameter predicted from the DFT are in good agreement with Monte Carlo simulation results. The theory describes satisfactorily the transition from depletion at low packing fractions to adsorption and double-layer adsorption at high packing fractions and also accounts for the crucial effect of the segment-wall interaction. It is found that homopolymer is better recognized at a low bulk density and a stronger interaction with the surface. The polymer can not only recognize the surface but also invert the surface at high bulk densities. The chain in the solution-wall interface exhibits a typical "brush" conformation with a length approximated by half the length of polymer chain.     

Dissecting the concave–convex π‐π interaction in corannulene and sumanene dimers: SAPT(DFT) analysis and performance of DFT dispersion‐corrected methods

The characteristics of the concave–convex π‐π interactions are evaluated in 32 buckybowl dimers formed by corannulene, sumanene, and two substituted sumanenes (with S and CO groups), using symmetry‐adapted perturbation theory [SAPT(DFT)] and density functional theory (DFT). According to our results, the main stabilizing contribution is dispersion, followed by electrostatics. Regarding the ability of DFT methods to reproduce the results obtained with the most expensive and rigorous methods, TPSS‐D seems to be the best option overall, although its results slightly tend to underestimate the interaction energies and to overestimate the equilibrium distances. The other two tested DFT‐D methods, B97‐D2 and B3LYP‐D, supply rather reasonable results as well. M06‐2X, although it is a good option from a geometrical point of view, leads to too weak interactions, with differences with respect to the reference values amounting to about 4 kcal/mol (25% of the total interaction energy). © 2017 Wiley Periodicals, Inc.     

Estimation of polymer-surface interfacial interaction strength by a contact AFM technique

Atomic force microscopy (AFM) measurements were employed to assess polymer-surface interfacial interaction strength. The main feature of the measurement is the use of contact-mode AFM as a tool to scratch off the polymer monolayer adsorbed on the solid surface. Tapping-mode AFM was used to determine the depth of the scraped recess. Independent determination of the layer thickness obtained from optical phase interference microscopy (OPIM) confirmed the depth of the AFM scratch. The force required for the complete removal of the polymer layer with no apparent damage to the substrate surface was determined. Polypropylene (PP), low-density polyethylene (PE), and PP-grafted-maleic anhydride (PP-g-ma) were scraped off silane-treated glass slabs, and the strength of surface interaction of the polymer layer was determined. In all cases it was determined that the magnitude of surface interaction force is of the order of van der Waals (VDW) interactions. The interaction strength is influenced either by polymer ability to wet the surface (hydrophobic or hydrophilic interactions) or by hydrogen bonding between the polymer and the surface treatment.     

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.     

A density functional theory approach to noncovalent interactions via interacting monomer densities

A recently proposed "DFT + dispersion" treatment (Rajchel et al., Phys. Rev. Lett., 2010, 104, 163001) is described in detail and illustrated by more examples. The formalism derives the dispersion-free density functional theory (DFT) interaction energy and combines it with the dispersion energy from separate DFT calculations. It consists of the self-consistent polarization of DFT monomers restrained by the exclusion principle via the Pauli blockade technique. Within the monomers a complete exchange-correlation potential should be used, but between them only the exact exchange operates. The application to a wide range of molecular complexes from rare-gas dimers to hydrogen-bonds to π-electron interactions shows good agreement with benchmark values.     

Structure, surface excess and effective interactions in polymer nanocomposite melts and concentrated solutions

The Polymer Reference Interaction Site Model (PRISM) theory is employed to investigate structure, effective forces, and thermodynamics in dense polymer-particle mixtures in the one and two particle limit. The influence of particle size, degree of polymerization, and polymer reduced density is established. In the athermal limit, the surface excess is negative implying an entropic dewetting interface. Polymer induced depletion interactions are quantified via the particle-particle pair correlation function and potential of mean force. A transition from (nearly) monotonic decaying, attractive depletion interactions to much stronger repulsive-attractive oscillatory depletion forces occurs at roughly the semidilute-concentrated solution boundary. Under melt conditions, the depletion force is extremely large and attractive at contact, but is proceeded by a high repulsive barrier. For particle diameters larger than roughly five monomer diameters, division of the force by the particle radius results in a nearly universal collapse of the depletion force for all interparticle separations. Molecular dynamics simulations have been employed to determine the depletion force for nanoparticles of a diameter five times the monomer size over a wide range of polymer densities spanning the semidilute, concentrated, and melt regimes. PRISM calculations based on the spatially nonlocal hypernetted chain closure for particle-particle direct correlations capture all the rich features found in the simulations, with quantitative errors for the amplitude of the depletion forces at the level of a factor of 2 or less. The consequences of monomer-particle attractions are briefly explored. Modification of the polymer-particle pair correlations is relatively small, but much larger effects are found for the surface excess including an energetic driven transition to a wetting polymer-particle interface. The particle-particle potential of mean force exhibits multiple qualitatively different behaviors (contact aggregation, steric stabilization, local bridging attraction) depending on the strength and spatial range of the polymer-particle attraction.