Computer Simulation of Adsorbed Polymer Chains with a Different Molecular Architecture

Simulations of simple models of polymer chains were carried out by the means of the dynamic Monte Carlo method. The model chains were confined to a simple cubic lattice. Three different chain architectures were studied: linear, star‐branched and ring chains. The polymer model chain interacted with an impenetrable surface with a simple contact attractive potential. It was found that size parameters of all these polymers obey scaling laws. The temperatures of the transitions from weakly to strongly adsorbed chain were determined. It was shown for weakly adsorbed chains that ring polymers are always ca. 50% more adsorbed than linear and star‐branched ones. The properties of adsorbed linear and star‐branched polymers are very similar in the length of chain and the strength of adsorption studied. Strongly adsorbed ring polymers are still more adsorbed but differences between all kinds of chains become less pronounced.     

Structure of Adsorbed Polymer Chains: A Monte Carlo Study

The structure of adsorbed polymer chains was studied using simplified lattice models. The model chains were adsorbed on an impenetrable surface with an attractive potential. The dynamic Monte Carlo simulations based in the Metropolis scheme were carried out using these models. The influence of the internal chain architecture (linear, star‐branched and ring chains) and the degree of adsorption on the chain's structure was studied. It was shown that for weakly adsorbed chain regime the ring polymers which exhibit an almost twice as high degree of adsorption compared to linear and star chains have a higher number of adsorbed parts of chain (trains). But the length of such train remains almost the same for all types of a polymer chain. Star‐branched chains exhibit a slightly different change in number and the mean length of trains, loops and tails with the temperature and the chain total length compared to two other types of chain.     

Conformational behavior of polymers adsorbed on nanotubes

The importance of hydrophobic interactions in determining polymer adsorption and wrapping of carbon nanotubes is still under debate. In this work, we concentrate on the effect of short-ranged weakly attractive hydrophobic interactions between polymers and nanotubes (modeled as an infinitely long and smooth cylindrical surface), neglecting all other interactions apart for chain flexibility. Using coarse-grained Monte Carlo simulation of such simplified systems, we find that uniform adsorption and wrapping of the nanotube occur for all degrees of chain flexibility for tubes with sufficiently large outer radii. However, the adsorbed conformations depend on chain stiffness, ranging from randomly adsorbed conformations of the flexible chain to perfect helical or multihelical conformations (in the case of more concentrated solutions) of the rigid chains. Adsorption appears to occur in a sequential manner, wrapping the nanotube nearly one monomer at a time from the point of contact. Once adsorbed, the chains travel on the surface of the cylinder, retaining their helical conformations for the semiflexible and rigid chains. Our findings may provide additional insight to experimentally observed ordered polymer wrapping of carbon nanotubes.     

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.     

Properties of Star-Branched Polymer Chains Near a Surface

We considered two model systems of star-branched polymers near an impenetrable surface. The model chains were constructed on a simple cubic lattice. Each star polymer consisted of f = 3 arms of equal length and the total number of segments was up to 799. The excluded volume effect was included into these models only and therefore the system was studied at good solvent conditions. In the first model system polymer chain was terminally attached with one arm to the surface. The grafted arm could slide along the surface. In the second system the star-branched chain was adsorbed on the surface and the strength of adsorption was were varied. The simulations were performed using the dynamic Monte Carlo method with local changes of chain conformations. The internal and local structures of a polymer layer were determined. The lateral diffusion and internal mobility of star-branched chains were studied as a function of strength of adsorption and the chain length. The lateral diffusion and internal mobility of star-branched chains were studied as a function of strength of adsorption and the chain length. It was shown that the behavior of grafted and weakly adsorbed chains was similar to that of a free three-dimensional polymer, while the strongly adsorbed chains behave as a two-dimensional system.     

Transitions of tethered polymer chains: a simulation study with the bond fluctuation lattice model

A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensional state space. The simulations' density-of-states results have been evaluated for interaction parameters spanning the range from good- to poor-solvent conditions and from repulsive to strongly attractive surfaces. In this work, we describe the simulation method and present results for the overall phase behavior and for some of the transitions. For adsorption in good solvent, we compare with Metropolis Monte Carlo data for the same model and find good agreement between the results. For the collapse transition, which occurs when the solvent quality changes from good to poor, we consider two situations corresponding to three-dimensional (hard surface) and two-dimensional (very attractive surface) chain conformations, respectively. For the hard surface, we compare tethered chains with free chains and find very similar behavior for both types of chains. For the very attractive surface, we find the two-dimensional chain collapse to be a two-step transition with the same sequence of transitions that is observed for three-dimensional chains: a coil-globule transition that changes the overall chain size is followed by a local rearrangement of chain segments.     

Star-branched polymers in an adsorbing slit: a Monte Carlo study

A coarse-grained model of star-branched polymer chains confined in a slit was studied. The slit was formed by two parallel impenetrable surfaces, which were attractive for polymer beads. The polymer chains were flexible homopolymers built of identical united atoms whose positions in space were restricted to the vertices of a simple cubic lattice. The chains were regular star polymers consisted of f = 3 branches of equal length. The chains were modeled in good solvent conditions and, thus, there were no long-range specific interactions between the polymer beads-only the excluded volume was present. Monte Carlo simulations were carried out using the algorithm based on a chain's local changes of conformation. The influence of the chain length, the distances between the confining surfaces, and the strength of the adsorption on the properties of the star-branched polymers was studied. It was shown that the universal behavior found previously for the dimension of chains was not valid for some dynamic properties. The strongly adsorbed chains can change their position so that they swap between both surfaces with frequency depending on the size of the slit and on the temperature only.     

Motion of a branched polymer chain in confinement: a Monte Carlo study

The aim of the study was a theoretical investigation of the polymer molecules located between two parallel and impenetrable surfaces which were also attractive for polymer segments. The chains were constructed of identical segments and were restricted to knots of a simple cubic lattice. Since the chains were at good solvent conditions the only interactions between the segments of the chain were the excluded volume. The properties of the model chains were determined by means of Monte Carlo simulations with a sampling algorithm based on the chain's local changes of conformation. The differences and similarities in the structure for different adsorption regimes and the size of the slit were shown and discussed. It was observed that at certain conditions the polymer chain was adsorbed at one of the confining surfaces, and then after a certain period of time it detached from this surface and approached the opposite wall; this switch was repeated many times. The influence of the strength of the adsorption, the size of the slit, and the chain's length on the frequency of these jumps were determined. The mechanism of the chain's motion during the switch was also shown.     

Adsorption of comb copolymers on weakly attractive solid surfaces

In this work continuum and lattice Monte Carlo simulation methods are used to study the adsorption of linear and comb polymers on flat surfaces. Selected polymer segments, located at the tips of the side chains in comb polymers or equally spaced along the linear polymers, are attracted to each other and to the surface via square-well potentials. The rest of the polymer segments are modeled as tangent hard spheres in the continuum model and as self-avoiding random walks in the lattice model. Results are presented in terms of segment-density profiles, distribution functions, and radii of gyration of the adsorbed polymers. At infinite dilution the presence of short side chains promotes the adsorption of polymers favoring both a decrease in the depletion-layer thickness and a spreading of the polymer molecule on the surface. The presence of long side chains favors the adsorption of polymers on the surface, but does not permit the spreading of the polymers. At finite concentration linear polymers and comb polymers with long side chains readily adsorb on the solid surface, while comb polymers with short side chains are unlikely to adsorb. The simple models of comb copolymers with short side chains used here show properties similar to those of associating polymers and of globular proteins in aqueous solutions, and can be used as a first approximation to investigate the mechanism of adsorption of proteins onto hydrophobic surfaces.     

A Monte Carlo study of the coil-to-globule transition of model polymer chains near an attractive surface

When a polymer chain in solution interacts with an atomically smooth solid substrate, its conformational properties are strongly modified and deviate substantially from those of chains in bulk. In this work, the interplay of two competing transitions that affect the conformations of polymer chains near an energetically attractive surface is studied by means of Monte Carlo simulations on a cubic lattice. The transition from an extended to a compact conformation of a polymer chain near an attractive wall, as solubility deteriorates, exhibits characteristics akin to the “coil-to-globule” transition in bulk. An effective θ-temperature is determined. Its role as the transition point is confirmed in a variety of ways. The nature of the coil-to-compact transition is not qualitatively different from that in the bulk. Adsorbed polymer chains may assume “globular” or “pancake” configurations depending on the competition among adsorption strength, cohesive energy, and entropy. In a very relevant range of conditions, the dependence of the adsorbate thickness on chain-length is intermediate between that of 3-d (“semidroplets”) and 2-d (“pancake”) objects. The focus of this study is on rather long polymer chains. Several crucial features of the transitions of the adsorbed chains are N-dependent and various aspects of the adsorption and “dissolution” process are manifested clearly only at the “long chain” limit. © 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 47: 2462–2476, 2009     

Nanoscale colloids in a freely adsorbing polymer solution: a Monte Carlo simulation study

A key issue in nanoscale materials and chemical processing is the need for thermodynamic and kinetic models covering colloid-polymer systems over the mesoscopic length scale (approximately 1-100 nm). We have applied Monte Carlo simulations to attractive nanoscale colloid-polymer mixtures toward developing a molecular basis for models of these complex systems. The expanded ensemble Monte Carlo simulation method is applied to calculate colloid chemical potentials (micro(c)) and polymer adsorption (gamma) in the presence of freely adsorbing Lennard-Jones (LJ) homopolymers (surface modifiers). gamma and micro(c) are studied as a function of nanoparticle diameter (sigma(c)), modifier chain length (n) and concentration, and colloid-polymer attractive strength over 0.3 < Rg/sigma(c) < 6 (Rg is the polymer radius of gyration). In the attractive regime, nanocolloid chemical potential decreases and adsorbed amount increases as sigma(c), or n is increased. The scaling of gamma with n from the simulations agrees with the theory of Aubouy and Raphael (Macromolecules 1998, 31, 4357) in the extreme limits of Rg/sigma(c). When Rg/sigma(c) is large, the "colloid" approaches a molecular size and interacts only locally with a few polymer segments and gamma approximately n. When Rg/sigma(c) is small, the system approaches the conventional colloid-polymer size regime where multiple chains interact with a single particle, and gamma approximately sigma(c)2, independent of n. In contrast, adsorption in the mesoscopic range of Rg/sigma(c) investigated here is represented well by a power law gamma approximately n(p), with 0 < p < 1 depending on concentration and LJ attractive strength. Likewise, the chemical potential from our results is fitted well with micro(c) approximately n(q)sigma(c)3, where the cubic term results from the sigma(c) dependence of particle surface area (approximately sigma(c)2) and LJ attractive magnitude (approximately sigma(c)). The q-exponent for micro(c) (micro(c) approximately n(q)) varies with composition and LJ attractive strength but is always very close to the power exponent for gamma (gamma approximately n(p)). This result leads to the conclusion that in attractive systems, polymer adsorption (and thus polymer-colloid attraction) dominates the micro(c) dependence on n, providing a molecular interpretation of the effect of adsorbed organic layers on nanoparticle stability and self-assembly.     

A multiple chain Monte Carlo method for atomistic simulation of high molecular weight polymer melts

A new Monte Carlo method is proposed for the simulation of bulk systems of atomistically detailed polymers. Each move consists of a configurational rearrangement of the atoms in a specified region of the material, rather than a specified molecule. Thus atoms within different chains may be displaced cooperatively in each Monte Carlo move. Here, the method is implemented for the case of melts of linear chains, where the bond lengths and bond angles are held constant during the move. The performance of the algorithm is examined for linear polyethylene systems with chain lengths of 100 and 1000 backbone atoms, under a range of conditions. The method shows a considerable potential as a very general and flexible tool for simulating realistic polymer materials, subject to a number of performances limiting factors which are described in detail.     

Structure and thermodynamics of protein-polymer solutions: effects of spatially distributed hydrophobic surface residues

Protein-polymer association in solution driven by a short-range attraction has been investigated using a simple coarse-grain model solved by Monte Carlo simulations. The effect of the spatial distribution of the hydrophobic surface residues of the protein on the adsorption of weakly hydrophobic polymers at variable polymer concentration, polymer length, and polymer stiffness has been considered. Structural data on the adsorbed polymer layer and thermodynamic properties, such as the free energy, energy, and entropy, related to the protein-polymer interaction were calculated. It was found that a more heterogeneous distribution of the surface residues promotes adsorption and that this also applies for different polymer concentrations, polymer chain lengths, and polymer flexibilities. Furthermore, the polymer adsorption onto proteins with more homogeneous surface distributions displayed larger sensitivity to polymer properties such as chain length and flexibility. Finally, a simple relation between the adsorption probability and the change in the free energy was found and rationalized by a simple two-state adsorption model.     

Monte Carlo simulation of mixed lennard-jones nonionic surfactant adsorption at the liquid/vapor interface

New Monte Carlo simulations are presented for nonionic surfactant adsorption at the liquid/vapor interface of a monatomic solvent specifically investigating the roles of tail attraction and binary mixtures of different tail lengths. Surfactant molecules consist of an amphiphilic chain with a solvophilic head and a solvophobic tail. All molecules in the system, solvent and surfactant, are characterized by the Lennard-Jones (LJ) potential. Adjacent atoms along the surfactant chain are connected by finitely extensible harmonic springs. Solvent molecules move via the Metropolis random-walk algorithm, whereas surfactant molecules move according to the continuum configurational bias Monte Carlo (CBMC) method. We generate thermodynamic adsorption and surface-tension isotherms and compare results quantitatively to single-surfactant adsorption (Langmuir, 2007, 23, 1835). Surfactant tail groups with attractive interaction lead to cooperative adsorption at high surface coverage and higher maximum adsorption at the interface than those without. Moreover, adsorption and surface-tension isotherms with and without tail attraction are identical at low concentrations, deviating only near maximum coverage. Simulated binary mixtures of surfactants with differing lengths give intermediate behavior between that of the corresponding single-surfactant adsorption and surface-tension isotherms both with and without tail attraction. We successfully predict simulated mixture results with the thermodynamically consistent ideal adsorbed solution (IAS) theory for binary mixtures of unequal-sized surfactants using only the simulations from the single surfactants. Ultimately, we establish that a coarse-grained LJ surfactant system is useful for understanding actual surfactant systems when tail attraction is important and for unequal-sized mixtures of amphiphiles.     

A simple model of stiff and flexible polymer chain adsorption: the influence of the internal chain architecture

In this study, we investigated the process of random sequential adsorption of stiff and flexible polymer chains on a two-dimensional square lattice. The polymer chains were represented by sequence of lattice points forming needles, T shapes, and crosses as well as flexible linear chains and star-branched chains consisted of three and four arms. The Monte Carlo method was employed to generate the model systems. The percolation threshold and the jamming threshold were determined for all systems under consideration. The influence of the chain length and the chain architecture on both thresholds was calculated and discussed. The changes in the ordering of the system were also studied.     

The Structure of Star-Branched Chains in a Confined Space

Summary. We studied the properties of a simplified model of star-branched polymers confined in a slit formed by two parallel and impenetrable surfaces. The chains were built of identical united atoms (segments) whose positions were restricted to vertices of a simple cubic lattice. The polymer excluded volume and polymer segment-surface contact interactions were also introduced into the model. The properties of the model chains were determined by means of Monte Carlo simulations with a Metropolis-type sampling algorithm based on local changes of chain’s conformation. The structure of star-branched chains was investigated and the influence of the confinement and the temperature on the chain dimensions and structure was studied. It was shown that for chains in the adsorbing slits their sizes do not exhibit a universal behavior contrary to confined athermal polymers. The polymers in narrow slits at higher temperatures still exhibited features of a three-dimensional chain. It was also shown that chains in small slits and at low temperatures were fully adsorbed at one of the surfaces but could also switch the surface rapidly.     

Monte Carlo simulations of the exchange kinetics of polymers adsorbed on a solid surface

The exchange kinetics of polymers adsorbing on a solid surface is extensively studied by dynamic Monte Carlo simulations. A model employed simulates a semidilute polymer solution placed in contact with a solid surface that attracts polymer segments by the adsorption interaction (χ_s). The exchange process of polymer chains, between the solution and the adsorbed polymer layer, is examined under various conditions. The exchange kinetics shows two characteristic regimes with increasing chain length. One is the diffusion‐controlled regime found with a small χ_s , and the other the detachment‐controlled regime with a large χ_s . These two regimes are well described by a kinetic theory. Various dynamic quantities show that the diffusion‐controlled regime is not due to sluggish dynamics near the surface, but rather to bulk diffusion of chains. The diffusion‐controlled regime found in this study is considered to appear at the high temperature limit.     

Semiflexible polymers grafted to a solid planar substrate: changing the structure from polymer brush to "polymer bristle"

Monte Carlo simulations are presented for a coarse-grained model of polymer brushes with polymers having a varying degree of stiffness. Both linear chains and ring polymers grafted to a flat structureless non-adsorbing substrate surface are considered. Applying good solvent conditions, it is shown that with growing polymer stiffness the brush height increases significantly. The monomer density profiles for the case of ring polymers (chain length N(R) = 64) are very similar to the case of corresponding linear chains (N(L) = 32, grafting density larger by a factor of two) in the case of flexible polymers, while slight differences appear with increasing stiffness. Evidence is obtained that the chain dynamics in brushes is slowed down dramatically with increasing stiffness. Very short stiff rings (N(R) ≤ 16) behave like disks, grafted to the substrate such that the vector, perpendicular to the disk plane, is oriented parallel to the substrate surface. It is suggested that such systems can undergo phase transitions to states with liquid crystalline order.     

Dynamics of star branched polymers in a matrix of linear chains — a Monte Carlo study

A simple cubic lattice model of the melt of 3-arm star-branched polymers of various length dissolved in a matrix of long linear chains (n₁ = 800 beads) is studied using a dynamic Monte Carlo method. The total polymer volume fraction is equal to 0,5, while the volume fraction of the star polymers is about ten times smaller. The static and dynamic properties of these systems are compared with the corresponding model systems of isolated star-branched polymers and with the melt of linear chains. It has been found that the number of dynamic entanglements for the star polymers with arm length up to 400 segments is too small for the onset of the arm retraction mechanism of polymer relaxation. In this regime dynamics of star-branched polymers is close to the dynamics of linear polymers at corresponding concentration and with equivalent chain length. The entanglement length for star polymers appears to be somewhat larger compared with linear chains.     

Simulation of interaction forces between nanoparticles in the presence of Lennard-Jones polymers: freely adsorbing homopolymer modifiers

The force between two nanoscale colloidal particles dispersed in a solution of freely adsorbing Lennard-Jones homopolymer modifiers is calculated using the expanded grand canonical Monte Carlo simulation method. We investigate the effect of polymer chain length (N), nanoparticle diameter (sigma(c)), and colloid-polymer interaction energy (epsilon(cp)) on polymer adsorption (Gamma) and polymer-induced forces (F(P)(r)) between nanoparticles in the full thermodynamic equilibrium condition. There is a strong correlation between polymer adsorption and the polymer-mediated nanoparticle forces. When the polymer adsorption is weak, as in the case of smaller diameters and short polymer chain lengths (sigma(c) = 5, N = 10), the polymers do not have any significant effect on the bare nanoparticle interactions. The adsorbed amount increases with increasing particle diameter, polymer chain length, and colloid-polymer interaction energy. In general, for strong polymer-particle adsorption the polymer-governed force profiles between nanoparticles show short-range repulsion and long-ranged attraction, suggesting that homopolymers would not be ideal for achieving stabilization in nanoparticle dispersions. The attraction is likely due to bridging, as well as polymer segment-segment interactions. The location and magnitude of attractive minimum in the force profile can be controlled by varying N and epsilon(cp). The results show partial agreement and some marked differences with previous theoretical and experimental studies of forces in the limit of flat walls in an adsorbing polymer solution. The difference could be attributed to incorporation of long-ranged colloid-polymer potential in our simulations and the influence of the curvature of the nanoparticles.