Colloid-polymer mixtures in the protein limit

Bottlebrush Midblocks Promote Colloidal Bridging of Telechelic Polymers

Telechelic polymers are effective rheological modifiers that bridge between associative constituents to form elastic networks. The performance of linear telechelic chains, however, is controlled by entropic forces and thus suffers from an upper limit on bridge formation. This work overcomes this limitation by utilizing telechelic triblock copolymers containing bottlebrush midblocks. By comparing the rheological properties of emulsions linked by telechelic bottlebrush polymers to those containing linear chains, we determined that telechelic polymers with bottlebrush midblocks form elastic networks more efficiently. These enhanced rheological properties arise from the high stiffness of the bottlebrush midblocks, which offsets the entropic stretching penalty for bridge formation, enabling them to more readily form networks. This molecular-level control over polymer conformation in complex fluids opens avenues for designing highly elastic networks with minimal polymeric additives.

Colloidal adsorption in planar polymeric brushes

The design of nano-functionalised membranes or channels, able to effectively adsorb pollutants in aqueous solutions, is a topic that is gaining a great deal of attention in the materials science community. With this work we explore, through a combination of scaling theories and molecular dynamics simulations, the adsorption of spherical non-deformable colloidal nanoparticles within planar polymeric brushes. Our strategy is twofold: first, we generalise the Alexander-de Gennes theory for planar homopolymeric brushes to the case of diblock copolymer brushes, then we map the adsorbing homopolymeric brushes onto a diblock copolymer system, where the adsorbed colloids and all interacting monomers are considered monomers in bad solvent and we apply the generalised scaling theory to this effective diblock copolymer. This allows the prediction of the average conformation of the grafted substrate, i.e. its average height, as a function of the amount of loaded particles, as well as the introduction of a continuous mapping between a homopolymeric brush, the fraction of loaded particles and the average height of the adsorbing substrate.

Surfactancy in a tadpole model of proteins

We model the environment of eukaryotic nuclei by representing macromolecules by only their entropic properties, with globular molecules represented by spherical colloids and flexible molecules by polymers. We put particular focus on proteins with both globular and intrinsically disordered regions, which we represent with 'tadpole' constructed by grafting single polymers and colloids together. In Monte Carlo simulations, we find these tadpoles support phase separation via depletion flocculation, and demonstrate several surfactant behaviours, including being found preferentially at interfaces and forming micelles in single phase solution. Furthermore, the model parameters can be tuned to give a tadpole a preference for either bulk phase. However, we find entropy too weak to drive these behaviours by itself at likely biological concentrations.

Effects of polymer nonideality on depletion-induced phase behaviour of colloidal disks and rods

Colloidal dispersions composed of either platelets or rods exhibit liquid crystalline phase behaviour that is strongly influenced by the addition of nonadsorbing polymers. In this work we examined how polymer segment-segment interactions affect this phase behaviour as compared to using either penetrable hard spheres (PHS) or ideal ('ghost') chains as depletants. We find that the simplified polymer description predicts the same phase diagram topologies as the more involved polymer descriptions. Therefore the PHS description is still adequate for qualitative predictions. For sufficiently large polymer sizes we find however that the precise polymer description significantly alters the locations of the phase coexistence regions. Especially the stability region of isotropic-isotropic coexistence is affected by the polymer interactions. To illustrate the quantitative effects some examples are presented.

Protein-polymer mixtures in the colloid limit: Aggregation, sedimentation, and crystallization

While proteins have been treated as particles with a spherically symmetric interaction, of course in reality, the situation is rather more complex. A simple step toward higher complexity is to treat the proteins as non-spherical particles and that is the approach we pursue here. We investigate the phase behavior of the enhanced green fluorescent protein (eGFP) under the addition of a non-adsorbing polymer, polyethylene glycol. From small angle x-ray scattering, we infer that the eGFP undergoes dimerization and we treat the dimers as spherocylinders with aspect ratio L/D - 1 = 1.05. Despite the complex nature of the proteins, we find that the phase behavior is similar to that of hard spherocylinders with an ideal polymer depletant, exhibiting aggregation and, in a small region of the phase diagram, crystallization. By comparing our measurements of the onset of aggregation with predictions for hard colloids and ideal polymers [S. V. Savenko and M. Dijkstra, J. Chem. Phys. 124, 234902 (2006) and Lo Verso et al., Phys. Rev. E 73, 061407 (2006)], we find good agreement, which suggests that the behavior of the eGFP is consistent with that of hard spherocylinders and ideal polymers.

Influence of solvent quality on depletion potentials in colloid-polymer mixtures

As first explained by the classic Asakura-Oosawa (AO) model, effective attractive forces between colloidal particles induced by depletion of nonadsorbing polymers can drive demixing of colloid-polymer mixtures into colloid-rich and colloid-poor phases, with practical relevance for purification of water, stability of foods and pharmaceuticals, and macromolecular crowding in biological cells. By idealizing polymer coils as effective penetrable spheres, the AO model qualitatively captures the influence of polymer depletion on thermodynamic phase behavior of colloidal suspensions. In previous work, we extended the AO model to incorporate aspherical polymer conformations and showed that fluctuating shapes of random-walk coils can significantly modify depletion potentials [W. K. Lim and A. R. Denton, Soft Matter 12, 2247 (2016); J. Chem. Phys. 144, 024904 (2016)]. We further demonstrated that the shapes of polymers in crowded environments sensitively depend on solvent quality [W. J. Davis and A. R. Denton, J. Chem. Phys. 149, 124901 (2018)]. Here, we apply Monte Carlo simulation to analyze the influence of solvent quality on depletion potentials in mixtures of hard-sphere colloids and nonadsorbing polymer coils, modeled as ellipsoids whose principal radii fluctuate according to random-walk statistics. We consider both self-avoiding and non-self-avoiding random walks, corresponding to polymers in good and theta solvents, respectively. Our simulation results demonstrate that depletion of polymers of equal molecular weight induces much stronger attraction between colloids in good solvents than in theta solvents and confirm that depletion interactions are significantly influenced by aspherical polymer conformations.

Bacterial Nucleoid: Interplay of DNA Demixing and Supercoiling

This work addresses the question of the interplay of DNA demixing and supercoiling in bacterial cells. Demixing of DNA from other globular macromolecules results from the overall repulsion between all components of the system and leads to the formation of the nucleoid, which is the region of the cell that contains the genomic DNA in a rather compact form. Supercoiling describes the coiling of the axis of the DNA double helix to accommodate the torsional stress injected in the molecule by topoisomerases. Supercoiling is able to induce some compaction of the bacterial DNA, although to a lesser extent than demixing. In this work, we investigate the interplay of these two mechanisms with the goal of determining whether the total compaction ratio of the DNA is the mere sum or some more complex function of the compaction ratios due to each mechanism. To this end, we developed a coarse-grained bead-and-spring model and investigated its properties through Brownian dynamics simulations. This work reveals that there actually exist different regimes, depending on the crowder volume ratio and the DNA superhelical density. In particular, a regime in which the effects of DNA demixing and supercoiling on the compaction of the DNA coil simply add up is shown to exist up to moderate values of the superhelical density. In contrast, the mean radius of the DNA coil no longer decreases above this threshold and may even increase again for sufficiently large crowder concentrations. Finally, the model predicts that the DNA coil may depart from the spherical geometry very close to the jamming threshold as a trade-off between the need to minimize both the bending energy of the stiff plectonemes and the volume of the DNA coil to accommodate demixing.

Confined crowded polymers near attractive surfaces

We present results from molecular dynamics simulations of a spherically confined neutral polymer in the presence of crowding particles, studying polymer shapes and conformations as a function of the strength of the attraction to the confining wall, solvent quality, and the density of crowders. The conformations of the polymer under good solvent conditions are weakly dependent on crowder particle density, even when the polymer is strongly confined. In contrast, under poor solvent conditions, when the polymer assumes a collapsed conformation when unconfined, it can exhibit transitions to two different adsorbed phases, when either the interaction with the wall or the density of crowder particles is changed. One such transition involves a desorbed collapsed phase change to an adsorbed extended phase as the attraction of the polymer towards the confining wall is increased. Such an adsorbed extended phase can exhibit a second transition to an ordered adsorbed collapsed phase as the crowder particle density is increased. The ordered adsorbed collapsed phase of the polymer differs significantly in its structure from the desorbed collapsed phase. We revisit the earlier understanding of the adsorption of confined polymers on attractive surfaces in light of our results.

Soft Interactions Modify the Diffusive Dynamics of Polymer-Grafted Nanoparticles in Solutions of Free Polymer

We examine the dynamics of silica particles grafted with high molecular weight polystyrene suspended in semidilute solutions of chemically similar linear polymer using X-ray photon correlation spectroscopy. The particle dynamics decouple from the bulk viscosity despite their large hydrodynamic size and instead experience an effective viscosity that depends on the molecular weight of the free polymer chains. Unlike for hard-sphere nanoparticles in semidilute polymer solutions, the diffusivities of the polymer-grafted nanoparticles do not collapse onto a master curve solely as a function of normalized length scales. Instead, the diffusivities can be collapsed across two orders of magnitude in free polymer molecular weight and concentration and one order of magnitude in grafted molecular weight by incorporating the ratio of free to grafted polymer molecular weights. These results suggest that the soft interaction potential between polymer-grafted nanoparticles and free polymer allows polymer-grafted nanoparticles to diffuse faster than predicted based on bulk rheology and modifies the coupling between grafted particle dynamics and the relaxations of the surrounding free polymer.

Preferential Localization of the Bacterial Nucleoid

Prokaryotes do not make use of a nucleus membrane to segregate their genetic material from the cytoplasm, so that their nucleoid is potentially free to explore the whole volume of the cell. Nonetheless, high resolution images of bacteria with very compact nucleoids show that such spherical nucleoids are invariably positioned at the center of mononucleoid cells. The present work aims to determine whether such preferential localization results from generic (entropic) interactions between the nucleoid and the cell membrane or instead requires some specific mechanism, like the tethering of DNA at mid-cell or periodic fluctuations of the concentration gradient of given chemical species. To this end, we performed numerical simulations using a coarse-grained model based on the assumption that the formation of the nucleoid results from a segregative phase separation mechanism driven by the de-mixing of the DNA and non-binding globular macromolecules. These simulations show that the abrupt compaction of the DNA coil, which takes place at large crowder density, close to the jamming threshold, is accompanied by the re-localization of the DNA coil close to the regions of the bounding wall with the largest curvature, like the hemispherical caps of rod-like cells, as if the DNA coil were suddenly acquiring the localization properties of a solid sphere. This work therefore supports the hypothesis that the localization of compact nucleoids at regular cell positions involves either some anchoring of the DNA to the cell membrane or some dynamical localization mechanism.

Quantifying the protein-protein association rate in polymer solutions: crowding-induced diffusion and energy modifications

A theoretical framework is developed to study protein-protein association in polymer solutions under diffusion-limited conditions. Starting from the basal association rate, two fundamental aspects concerning macromolecular crowding are particularly taken into account. One is the effect of microviscosity on protein diffusion. The length-scale dependent relations of translational and rotational diffusion coefficients are incorporated. Another is relevant to the crowding-induced effective interaction between the pair of proteins. The resultant energy modifications to the basal association rate are properly introduced, following an instructive classification with increasing crowder size: repulsive interaction dominant, repulsive and attractive interactions competitive, and attractive interaction prevailing. With specific energy modification terms, we are able to investigate the deviations of the association rate from the Stokes-Einstein (SE) behavior (i.e., the linear relationship with respect to the reciprocal of macroviscosity) in a quantitative manner. Our theory is applied to study the association of TEM1-β-lactamase (TEM) and the β-lactamase inhibitor protein (BLIP) in polyethylene glycol (PEG) solutions, with varying concentration and polymerization. We explicitly evaluate the relative association rate constant as a function of the solution macroviscosity. The theoretical results demonstrate very good agreement with the experimental data. Moreover, the complicated non-trivial deviations, either positive (slower than SE) or negative (faster than SE), are systematically rationalized. The precise role of energy modifications and the microviscosity effect on diffusion, in particular on rotational diffusion, are clearly unraveled.

Many-body effects in a binary nano-particle mixture dispersed in ideal polymer solutions

A new mean-field theory is developed to treat a binary mixture of nanoparticles imbedded in a polydisperse polymer solution. The theory is based on a many-body polymer-mediated potential of mean force (PMF) between the particles and remains accurate even in the protein regime, where the particles' diameters cannot necessarily be considered large compared to the polymer radius of gyration. As implemented here, the theory is strictly valid for dilute to semi-dilute polymer solutions near the theta temperature (the so-called theta regime) or when the range of the PMF is strongly affected by the polymer size. For non-adsorbing particles, this is the same regime where the celebrated Asakura-Oosawa (AO) model is often used. Unlike the traditional AO model, however, our approach includes polymer flexibility and is accurate in the protein regime. We use the theory to calculate phase diagrams for a binary mixture of unequal-sized particles, both adsorbing and non-adsorbing. To test the theory, we carry out comparisons with simulations and obtained good quantitative agreement, which gives support to its accuracy. On the other hand, the oft-used approach assuming pairwise-additive potentials of mean force produce quantitatively (and sometime qualitatively) different phase diagrams.

Influence of solvent quality on conformations of crowded polymers

The structure and function of polymers in confined environments, e.g., biopolymers in the cytoplasm of a cell, are strongly affected by macromolecular crowding. To explore the influence of solvent quality on conformations of crowded polymers, we model polymers as penetrable ellipsoids, whose shape fluctuations are governed by the statistics of self-avoiding walks, appropriate for a polymer in a good solvent. Within this coarse-grained model, we perform Monte Carlo simulations of mixtures of polymers and hard-nanosphere crowders, including trial changes in polymer size and shape. Penetration of polymers by crowders is incorporated via a free energy cost predicted by polymer field theory. To analyze the impact of crowding on polymer conformations in different solvents, we compute the average polymer shape distributions, radius of gyration, volume, and asphericity over ranges of the polymer-to-crowder size ratio and crowder volume fraction. The simulation results are accurately predicted by a free-volume theory of polymer crowding. Comparison of results for polymers in good and theta solvents indicates that excluded-volume interactions between polymer segments significantly affect crowding, especially in the limit of crowders much smaller than polymers. Our approach may help to motivate future experimental studies of polymers in crowded environments, with possible relevance for drug delivery and gene therapy.

Perspective: Outstanding theoretical questions in polymer-nanoparticle hybrids

This topical review discusses the theoretical progress made in the field of polymer nanocomposites, i.e., hybrid materials created by mixing (typically inorganic) nanoparticles (NPs) with organic polymers. It primarily focuses on the outstanding issues in this field and is structured around five separate topics: (i) the synthesis of functionalized nanoparticles; (ii) their phase behavior when mixed with a homopolymer matrix and their assembly into well-defined superstructures; (iii) the role of processing on the structures realized by these hybrid materials and the role of the mobilities of the different constituents; (iv) the role of external fields (electric, magnetic) in the active assembly of the NPs; and (v) the engineering properties that result and the factors that control them. While the most is known about topic (ii), we believe that significant progress needs to be made in the other four topics before the practical promise offered by these materials can be realized. This review delineates the most pressing issues on these topics and poses specific questions that we believe need to be addressed in the immediate future.

Self-assembly of polymer-linked nanoparticles and scaling behavior in the assembled phase

An entropic depletion-driven phase separation is known to be observed for mixtures of polymers and nanoparticles. While polymer-linked nanoparticles have been synthesized, their phase behavior has only been predicted for chemically specific interactions. We use integral equation theory to determine the structure and phase behavior of chemically isotropic polymer-linked nanoparticles at high densities. When each end of a linear polymer is grafted to a nanoparticle, we predict an entropy-driven microphase separation of locally segregated polymer-rich and nanoparticle-rich domains. The formation of these self-assembled structures is purely a consequence of the shape of the polymer-linked particle species. The depletion-driven demixing of ungrafted polymer-nanoparticle composites (with small amounts of nanoparticles) is enhanced as particle diameter (D) grows compared to the polymer radius of gyration (R_g). However, this study shows that for polymer-linked nanoparticle systems, the transition from a liquid to microphase separated state shifts to higher densities (i.e. is inhibited) as D/R_g increases. The transition volume fractions attain a unique value (of ∼0.69) at D/R_g ∼ 1.13. The repeating length scale (L*) is 1.4-2.2 times the size of the entire species (D + R_g). Surprisingly, L*/(D + R_g) is a non-monotonic function of the polymer radius of gyration. The repeating length scale also displays a remarkable scaling behavior, as a function of the particle diameter and the polymer density. Additionally, our study implies that two different mechanisms of transitioning to the microphase separated state are possible for these systems, which has important implications for the transition density and the kinds of structures formed.

Tuning the phase diagram of colloid-polymer mixtures via Yukawa interactions

Theory that predicts the phase behavior of interacting Yukawa spheres in a solution containing nonadsorbing polymer is presented. Our approach accounts for multiple overlap of depletion zones. It is found that additional Yukawa interactions beyond hard core interactions strongly affect the location and presence of coexistence regions and phase states. The theoretical phase diagrams are compared with Monte Carlo simulations. The agreement between the two approaches supports the validity of the theoretical approximations made and confirms that, by choosing the parameters of the interaction potentials, tuning of the binodals is possible. The critical end point characterizes the phase diagram topology. It is demonstrated how an additional Yukawa interaction shifts this point with respect to the hard sphere case. Provided a certain depletant-to-colloid size ratio for which a stable colloidal gas-liquid phase coexistence takes place for hard spheres, added direct interactions turn this into a metastable gas-liquid equilibrium. The opposite case, the induction of a stable gas-liquid coexistence where only fluid-solid was present for hard spheres, is also reported.

Contact line of adsorbed colloid-polymer droplets in theory and experiment

The contact line between the colloid-rich bulk liquid and an adsorbed thin film in colloid-polymer mixtures (CPM) is studied by means of an interface displacement model. The interface displacement profiles are compared to laser scanning confocal microscopy (LSCM) images. The mixtures consist of poly(methylmetacrylate) (PMMA) colloids and polystyrene (PS) polymers with polymer-to-colloid size ratio q = 1.18. Based on the experimental parameters, the theoretical model predicts a contact angle for colloid-rich liquid droplets adsorbed on glass of θ∞ = 59°, assuming a contact line with infinite radius, R = ∞. When a contact-line curvature correction and a correction for the protein-limit character of the CPM are taken into account, a modest shift is obtained. The refined theory predicts θ≈ 56°. The contact angle determined visually from the LSCM images is θ≈ 30°. The model predicts a three-phase contact-line tension of τ = -1.2 × 10(-12) N (uncorrected) and τ = -2.3 × 10(-13) N (with protein-limit correction), which is physically sound both in sign and magnitude. The line tension influences the contact angle to a small extent due to the contact line curvature. The predicted width of the transition zone between the thin film and the droplet is about 2 μm. The effect of gravity is noticeable as a deformation near the middle of the droplet.

Size-dependent interaction of silica nanoparticles with lysozyme and bovine serum albumin proteins

The interaction of three different sized (diameter 10, 18, and 28 nm) anionic silica nanoparticles with two model proteins-cationic lysozyme [molecular weight (MW) 14.7 kDa)] and anionic bovine serum albumin (BSA) (MW 66.4 kDa) has been studied by UV-vis spectroscopy, dynamic light scattering (DLS), and small-angle neutron scattering (SANS). The adsorption behavior of proteins on the nanoparticles, measured by UV-vis spectroscopy, is found to be very different for lysozyme and BSA. Lysozyme adsorbs strongly on the nanoparticles and shows exponential behavior as a function of lysozyme concentration irrespective of the nanoparticle size. The total amount of adsorbed lysozyme, as governed by the surface-to-volume ratio, increases on lowering the size of the nanoparticles for a fixed volume fraction of the nanoparticles. On the other hand, BSA does not show any adsorption for all the different sizes of the nanoparticles. Despite having different interactions, both proteins induce similar phase behavior where the nanoparticle-protein system transforms from one phase (clear) to two phase (turbid) as a function of protein concentration. The phase behavior is modified towards the lower concentrations for both proteins with increasing the nanoparticle size. DLS suggests that the phase behavior arises as a result of the nanoparticles' aggregation on the addition of proteins. The size-dependent modifications in the interaction potential, responsible for the phase behavior, have been determined by SANS data as modeled using the two-Yukawa potential accounting for the repulsive and attractive interactions in the systems. The protein-induced interaction between the nanoparticles is found to be short-range attraction for lysozyme and long-range attraction for BSA. The magnitude of attractive interaction irrespective of protein type is enhanced with increase in the size of the nanoparticles. The total (attractive+repulsive) potential leading to two-phase formation is found to be more attractive for larger sized nanoparticles. The nanoparticle aggregates are characterized by mass fractal.

Influence of polymer shape on depletion potentials and crowding in colloid-polymer mixtures

Depletion-induced interactions between colloids in colloid-polymer mixtures depend in range and strength on size, shape, and concentration of depletants. Crowding by colloids in turn affects shapes of polymer coils, such as biopolymers in biological cells. By simulating hard-sphere colloids and random-walk polymers, modeled as fluctuating ellipsoids, we compute depletion-induced potentials and polymer shape distributions. Comparing results with exact density-functional theory calculations, molecular simulations, and experiments, we show that polymer shape fluctuations play an important role in depletion and crowding phenomena.

Depletion, melting and reentrant solidification in mixtures of soft and hard colloids

We present extensive experimental and theoretical investigations on the structure, phase behavior, dynamics and rheology of model soft-hard colloidal mixtures realized with large, multiarm star polymers as the soft component and smaller, compact stars as the hard one. The number and length of the arms in star polymers control their softness, whereas the size ratio, the overall density and the composition are additional parameters varied for the mixtures. A coarse-grained theoretical strategy is employed to predict the structure of the systems as well as their ergodicity properties on the basis of mode coupling theory, for comparison with rheological measurements on the samples. We discovered that dynamically arrested star-polymer solutions recover their ergodicity upon addition of colloidal additives. At the same time the system displays demixing instability, and the binodal of the latter meets the glass line in a way that leads, upon addition of a sufficient amount of colloidal particles, to an arrested phase separation and reentrant solidification. We present evidence for a subsequent solid-to-solid transition well within the region of arrested phase separation, attributed to a hard-sphere-mixture type of glass, due to osmotic shrinkage of the stars at high colloidal particle concentrations. We systematically investigated the interplay of star functionality and size ratio with glass melting and demixing, and rationalized our findings by the depletion of the big stars due to the smaller colloids. This new depletion potential in which, contrary to the classic colloid-polymer case, the hard component depletes the soft one, has unique and novel characteristics and allows the calculation of phase diagrams for such mixtures. This work covers a broad range of soft-hard colloidal mixture compositions in which the soft component exceeds the hard one in size and provides general guidelines for controlling the properties of such complex mixtures.

Polymer-mediated interactions and their effect on the coagulation-fragmentation of nano-colloids: a self-consistent field theory approach

This feature paper reviews our recent efforts to theoretically model the effect of polymer mediated interactions on the coagulation-fragmentation of nano-colloids in different settings encountered in practical systems. The polymer-mediated interactions among nanoparticles play a key role in many biological and technological processes such as red blood cell aggregation, protein crystallization, self-healing of polymer composites, filler reinforcement of rubbers used in tire technology, etc. By developing and making use of the novel potential theory, we investigate several important cases of these interactions acting between nanoparticles in diverse nano-polymer composites. As a demonstration of its practical applicability, we use the developed theory to investigate the effect of polymer mediated interactions on the coagulation-fragmentation of fillers and their kinetic stability in the presence of non-adsorbing and adsorbing polymers. In particular, we use our findings to develop a pragmatic way of evaluating the kinetic stability of nano-filler agglomerates critical for understanding the filler reinforcement of rubbers. Finally, we perform thorough comparison of the present theoretical findings with the available experimental data and simulations.

Phase diagram of mixtures of colloids and polymers in the thermal crossover from good to θ solvent

We determine the phase diagram of mixtures of spherical colloids and neutral nonadsorbing polymers in the thermal crossover region between the θ point and the good-solvent regime. We use the generalized free-volume theory, which takes into account the polymer-concentration dependence of the depletion thickness and of the polymer compressibility. This approach turns out to be quite accurate as long as q = Rg/Rc ≲ 1 (Rg is the radius of gyration of the polymer and Rc is the colloid radius). We find that, close to the θ point, the phase diagram is not very sensitive to solvent quality, while, close to the good-solvent region, changes of the solvent quality modify significantly the position of the critical point and of the binodals. We also analyze the phase behavior of aqueous solutions of charged colloids and polymers, using the approach proposed by Fortini et al. [J. Phys.: Condens. Matter 17, 7783 (2005)].

Theory and simulation studies of effective interactions, phase behavior and morphology in polymer nanocomposites

Polymer nanocomposites are a class of materials that consist of a polymer matrix filled with inorganic/organic nanoscale additives that enhance the inherent macroscopic (mechanical, optical and electronic) properties of the polymer matrix. Over the past few decades such materials have received tremendous attention from experimentalists, theoreticians, and computational scientists. These studies have revealed that the macroscopic properties of polymer nanocomposites depend strongly on the (microscopic) morphology of the constituent nanoscale additives in the polymer matrix. As a consequence, intense research efforts have been directed to understand the relationships between interactions, morphology, and the phase behavior of polymer nanocomposites. Theory and simulations have proven to be useful tools in this regard due to their ability to link molecular level features of the polymer and nanoparticle additives to the resulting morphology within the composite. In this article we review recent theory and simulation studies, presenting briefly the methodological developments underlying PRISM theories, density functional theory, self-consistent field theory approaches, and atomistic and coarse-grained molecular simulations. We first discuss the studies on polymer nanocomposites with bare or un-functionalized nanoparticles as additives, followed by a review of recent work on composites containing polymer grafted or functionalized nanoparticles as additives. We conclude each section with a brief outlook on some potential future directions.

Structure of phase-separated athermal colloid-polymer systems in the protein limit

Structural features of phase-separated athermal colloid-polymer mixtures in the so-called "protein limit," where polymer chain dimensions exceed those of the colloid, are investigated using grand canonical Monte Carlo simulations on a fine lattice. Previous work [N. A. Mahynski et al., Phys. Rev. E 85, 051402 (2012)] has shown that this model accurately captures the phase behavior of experimental systems, and that colloids with sufficiently small diameters, σ(c), relative to that of the monomeric segments, σ(s), phase separate more readily than their large-diameter counterparts. In the present study, we directly connect colloid and polymer structure with their phase behavior by investigating these solutions along their binodal curves; we also explore the role of colloid surface curvature in destabilizing such solutions. Our findings suggest that simple consideration of an additional depletion radius, on the order of the σ(s), leads to a quantitatively accurate prediction of the division between stable and unstable ranges of d=σ(s)/σ(c). We compare these results to continuum models with different bonding potentials between monomer segments in order to elucidate the significance of the lattice model's bond fluctuations and inherently coarse colloid surface. In a number of cases, the continuum models deviate both qualitatively and quantitatively from the lattice results, but the binodals of the continuum models are presently not known, making a strong conclusion about these differences impossible.

Entropy-driven segregation of nanoparticles in polymer melts

We present a Monte Carlo study of the dispersibility of nanoparticles in polymer melts based on the relative entropic contributions of particles and polymer chains. Our model predictions are in terms of the particle diameter d and the chain statistical segment length ℓ. In the limit of small (d=ℓ) particles, a good dispersion is predicted, as the high entropy contribution of the small particles largely offsets the loss in polymer chain entropy. This is not the case, however, at larger d ≈ 3ℓ and particle segregation at the surfaces of thin polymer films is observed. That model prediction is in good agreement with experimental observation for systems with little or no enthalpic interaction between polymer and particles.

Immersion energy and polymer-mediated depletion interactions between nanocolloids as studied by analytic self-consistent field theory

By applying the Edwards self-consistent field theory, we calculate the polymer density and free energy excesses caused by the presence of nanocolloids in the excluded volume polymer system. Using the obtained results, we have calculated the depletion potential U as a function of the separation between colloids, colloid radius, polymer volume fraction, and polymer gyration radius. Upon analyzing the obtained results against the known exact asymptotic small-separation limit and scaling relations, we propose an approximate expression for U that builds upon our exact result. This expression is shown to give excellent agreement with Monte Carlo simulations. Reliability of the derived analytical expression for the depletion potential and its relation to the previous theoretical and simulation work are thoroughly discussed.

Pressure and density scaling for colloid-polymer systems in the protein limit

Grand canonical Monte Carlo and histogram reweighting techniques are used to study the fluid-phase behavior of an athermal system of colloids and nonadsorbing polymers on a fine lattice in the "protein limit," where polymer dimensions exceed those of the colloids. The main parameters are the chains' radius of gyration, R_{g}, the diameter of the colloids, σ_{c}, and the monomer diameter, σ_{s}. The phase behavior is controlled by the macroscopic size ratio, q_{r}=2R_{g}/σ_{c}, and the microscopic size ratio, d=σ_{s}/σ_{c}. The latter ratio is found to play a significant role in determining the critical monomer concentration for q_{r}≲4 and the critical colloid density for all chain lengths. However, the critical (osmotic) pressure is independent of the microscopic size ratio at all macroscopic size ratios studied. Quantitative agreement is observed between our simulation results and experimental data. We scale our results based on the polymer correlation length, which has previously been suggested to universally collapse these binodals [Bolhuis et al., Phys. Rev. Lett. 90, 068304 (2003); Fleer and Tuinier, Phys. Rev. E 76, 041802 (2007)]. While the density binodals exhibit universal characteristics along the low-colloid-density branch, such features are not present in the corresponding high-density phase. However, pressure binodals do collapse nicely under such a scaling, even far from the critical point, which allows us to produce a binodal curve whose shape is independent of either size ratio.