Simulation of short-chain polymer collapse with an explicit solvent

Coil-globule transition in two-dimensional polymer chains in an explicit solvent

The structure of two-dimensional polymer chains in a solvent at different temperatures is still far from being fully understood. Computer simulations of high-density macromolecular systems require the use of appropriate algorithms, and therefore the simulations were carried out using the Cooperative Motion Algorithm. The polymer model studied was exactly two-dimensional, coarse-grained and based on a triangular lattice. The theta temperature and temperature of coil-to-globule transition, and critical exponents were determined. The differences between the structure of such a disk and that of a chain in a dense polymer liquid were shown.

Dynamic Feature of Incipient Polymer Collapse below the Theta Point

We study the dynamics of a polymer chain with gradually changing the solvent quality from good to poor by dissipative particle dynamics simulation. We find several spectral modes related to internal motions of intrachain interaction. Approaching the coil-to-globule transition point, all fast modes of spectrum ω > 1 (ns)^-1 disappear. There is only a slow mode at ω ≈ 0.66 (ns)^-1. Moreover, the spectral density at this slow mode reaches a maximum value at the transition point. We suggest that, at the transition point, the chain conformation relaxes to the most probable distribution only by the slow mode. There is a critical slowing down of internal motion with passing through the transition point.

Kinetics of polymer collapse: effect of temperature on cluster growth and aging

Using state of the art Monte Carlo simulations of a bead-spring model we investigate both the equilibrium and the nonequilibrium behavior of the homopolymer collapse. The equilibrium properties obtained via multicanonical sampling recover the well-known finite-size scaling behavior of collapse for our model polymer. For the nonequilibrium dynamics we study the collapse by quenching the homopolymer from an expanded coiled state into the globular phase. The sequence of events observed during the collapse is independent of the quench depth. In particular, we focus on finding out universal scaling behaviors related to the growth or coarsening of clusters of monomers, by drawing phenomenological analogies with ordering kinetics. We distinguish the cluster coarsening stage from the initial stage of primary cluster formation. By successful application of a nonequilibrium finite-size scaling analysis we show that at all quench temperatures, during the coarsening stage, the cluster growth is roughly linear and can be characterised by a universal finite-size scaling function. In addition, we provide evidence of aging by constructing a suitable autocorrelation function and its corresponding dynamical power-law scaling with respect to the growing cluster sizes. The predicted theoretical bound for the exponent governing such scaling is strictly obeyed by the numerical data irrespective of the quench temperature. The results and methods presented here in general should find application in similar phenomena such as the collapse of a protein molecule preceding its folding.

Conformations of Low-Molecular-Weight Lignin Polymers in Water

Low-molecular-weight lignin binds to cellulose during the thermochemical pretreatment of biomass for biofuel production, which prevents the efficient hydrolysis of the cellulose to sugars. The binding properties of lignin are influenced strongly by the conformations it adopts. Here, we use molecular dynamics simulations in aqueous solution to investigate the dependence of the shape of lignin polymers on chain length and temperature. Lignin is found to adopt collapsed conformations in water at 300 and 500 K. However, at 300 K, a discontinuous transition is found in the shape of the polymer as a function of the chain length. Below a critical degree of polymerization, Nc =15, the polymer adopts less spherical conformations than above Nc. The transition disappears at high temperatures (500 K) at which only spherical shapes are adopted. An implication relevant to cellulosic biofuel production is that lignin will self-aggregate even at high pretreatment temperatures.

Directed Assembly of Soft Colloids through Rapid Solvent Exchange

We studied the directed assembly of soft nanoparticles through rapid micromixing of polymers in solution with a nonsolvent. Both experiments and computer simulations were performed to elucidate the underlying physics and to investigate the role of various process parameters. In particular, we discovered that no external stabilizing agents or charged end groups are required to keep the colloids separated from each other when water is used as the nonsolvent. Furthermore, the size of the nanoparticles can be reliably tuned through the mixing rate and the ratio between polymer solution and nonsolvent. Our results demonstrate that this mechanism is highly promising for the mass fabrication of uniformly sized colloidal particles, using a wide variety of polymeric feed materials.

Effect of confinement on the collapsing mechanism of a flexible polymer chain

In this paper, Brownian dynamics simulation (BDS) studies are executed to demonstrate the distinctive influences of the extent of confinement on the collapsing mechanism and kinetics of a flexible hydrophobic polymer chain in a poor solvent. The collapsing behavior is quantified by the time of collapse, which below a critical dimension of the confinement (h(c)), encounters a drastic reduction with a further strengthening in the degree of confinement. For dimensions greater than this critical one, the collapse occurs through the well-known hydrodynamic interaction (HI) controlled multiple-globule-mediated mechanisms. However, for channel dimensions less than this critical one, the collapse mechanism is drastically altered. Under such circumstances, the collapse gets predominantly controlled by the confinement effects (with negligible contribution of the HIs) and occurs via the formation of a single central globule. This central globule rapidly engulfs the noncondensed polymer segments, and in the process largely hastens up the collapsing event. Under such circumstances, the collapse time is found to decrease linearly with decrements in the channel height. On the contrary, for channel heights greater than h(c), the multiple-globule-mediated collapse is characterized by a collapse time that shows an exponential dependence on the channel height, rapidly attaining a state in which the confinement effect becomes inconsequential and HIs dictate the entire collapsing behavior. We further propose detailed arguments based on physical reasoning as well as free energy estimations to conclusively support the qualitative and quantitative nature of influences of the confinement on the polymer collapse.

Probing solvation decay length in order to characterize hydrophobicity-induced bead-bead attractive interactions in polymer chains

In this paper, we quantitatively demonstrate that exponentially decaying attractive potentials can effectively mimic strong hydrophobic interactions between monomer units of a polymer chain dissolved in aqueous solvent. Classical approaches to modeling hydrophobic solvation interactions are based on invariant attractive length scales. However, we demonstrate here that the solvation interaction decay length may need to be posed as a function of the relative separation distances and the sizes of the interacting species (or beads or monomers) to replicate the necessary physical interactions. As an illustrative example, we derive a universal scaling relationship for a given solute-solvent combination between the solvation decay length, the bead radius, and the distance between the interacting beads. With our formalism, the hydrophobic component of the net attractive interaction between monomer units can be synergistically accounted for within the unified framework of a simple exponentially decaying potential law, where the characteristic decay length incorporates the distinctive and critical physical features of the underlying interaction. The present formalism, even in a mesoscopic computational framework, is capable of incorporating the essential physics of the appropriate solute-size dependence and solvent-interaction dependence in the hydrophobic force estimation, without explicitly resolving the underlying molecular level details.

Brownian dynamics simulation of polymer collapse in a poor solvent: influence of implicit hydrodynamic interactions

The effect of solvent on the collapse dynamics of homopolymers is investigated with Brownian dynamics simulations of a non-linear bead-spring chain model incorporating implicit hydrodynamic interactions. Our simulations suggest that the polymer collapse takes place via a three-stage mechanism, namely, formation of pearls, coarsening of pearls and the formation of a compact globule. The collapse pathways from a good solvent state to a poor solvent state are found to be independent of hydrodynamic interactions. On the other hand, hydrodynamic interaction is found to speed up the collapse rate. At a large quench depth (the depth of the Lennard-Jones potential), independent of the presence of hydrodynamic interaction, polymer molecules are found to be trapped in metastable states for long periods before acquiring their native globular state. The exponents characterizing the decay of various properties such as the radius of gyration are determined and compared with the values reported in the literature.

Solvation potentials for flexible chain molecules in solution: on the validity of a pairwise decomposition

The effects of a solvent on the conformation of a flexible n-site solute molecule can be described formally in terms of an n-body solvation potential. Given the practical difficulty in computing such multibody potentials, it is common to carry out a pairwise decomposition in which the n-body potential is approximated by a sum of two-body potentials. Here we investigate the validity of this two-site approximation for short interaction-site chain-in-solvent systems. Using exact expressions for the conformation of an isolated chain, we construct a mapping between the full chain-in-solvent system and its solvation potential representation. We present results for both hard-sphere and square-well systems with n=5 that show that the two-site approximation is sufficient to completely capture the effects of an explicit solvent on chain conformation for a wide range of conditions (which include varying the solvent diameter in the hard-sphere system and varying the chain-solvent coupling in the square-well system). In all cases, a set of two-site potentials (one for each distinct site-site pair) is required. We also show that these two-site solvation potentials can be used to accurately compute a multisite intramolecular correlation function.

Simulation study of the coil-globule transition of a polymer in solvent

Molecular dynamics simulations are used to study the coil-globule transition for a system composed of a bead-spring polymer immersed in an explicitly modeled solvent. Two different versions of the model are used, which are differentiated by the nature of monomer-solvent, solvent-solvent, and nonbonded monomer-monomer interactions. For each case, a model parameter lambda determines the degree of hydrophobicity of the monomers by controlling the degree of energy mismatch between the monomers and solvent particles. We consider a lambda-driven coil-globule transition at constant temperature. The simulations are used to calculate average static structure factors, which are then used to determine the scaling exponents of the system in order to determine the theta-point values lambdatheta separating the coil from the globule states. For each model we construct coil-globule phase diagrams in terms of lambda and the particle density rho. The results are analyzed in terms of a simple Flory-type theory of the collapse transition. The ratio of lambdatheta for the two models converges in the high density limit exactly to the value predicted by the theory in the random mixing approximation. Generally, the predicted values of lambdatheta are in reasonable agreement with the measured values at high rho, though the accuracy improves if the average chain size is calculated using the full probability distribution associated with the polymer-solvent free energy, rather than merely using the value obtained from the minimum of the free energy.

Discontinuous molecular dynamics simulation study of polymer collapse

Discontinuous molecular dynamics simulations were used to study the coil-globule transition of a polymer in an explicit solvent. Two different versions of the model were employed, which are differentiated by the nature of monomer-solvent, solvent-solvent, and nonbonded monomer-monomer interactions. For each case, a model parameter lambda determines the degree of hydrophobicity of the monomers by controlling the degree of energy mismatch between the monomers and solvent particles. We consider a lambda-driven coil-globule transition at constant temperature. The simulations are used to calculate average static structure factors, which are then used to determine the scaling exponents of the system in order to determine the theta-point values lambda(theta) separating the coil from the globule state. For each model we construct coil-globule phase diagrams in terms of lambda and the particle density rho. Additionally, we explore for each model the effects of varying the range of the attractive interactions on the phase boundary separating the coil and globule phases. The results are analyzed in terms of a simple Flory-type theory of the collapse transition.

Detailed Structures and Mechanism of Polymer Solvation

In the present study, we simulated a model system, PE in biphenyl, to gain the insight into the detailed solvation structures and the molecular mechanism of polymer chain solvation. Using atomistic molecular dynamics (MD) simulation, it was found that when the biphenyl is far from PE chain or in the bulk, the dihedral angle of the two rings in the solvent molecule are approximately 32 degrees. But, the dihedral angel is about 27 degrees when the biphenyls are very close to the PE chain. In the first solvation shell, the orientation angle of the biphenyl long axis to the chain segment backbone was found to be enhanced around two values: approximately 0 and approximately 60 degrees. The detailed solvation structures found here include all dyad conformations (TT, TG, TG', GT, GG, GG', G'T, G'G, and G'G') and vary as a function of the distance between PE chain and biphenyls in the first solvation shell. The closer the the solvent molecule to the PE segment, the higher the TT conformation fraction response is. The other dyad conformations such as TG, GG', etc. undergo different decreases, respectively. This study shows that the solvation even in the Theta condition makes the overall size expansion or the chain stretched. Such a cooperative change was examined here and found not due to generating or losing a conformational state but due to a change in conformational distribution. This change occurs in the middle location of the chain instead of the chain end locations.

Mesoscopic description of solvent effects on polymer dynamics

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.

Conformation of a polymer chain in solution: An exact density expansion approach

The conformation of a polymer chain in solution is intrinsically coupled to the thermodynamic and structural properties of the solvent. Here we study such solvent effects in a system consisting of a flexible interaction-site n-mer chain immersed in a monomeric solvent. Chain conformation is described with a set of intramolecular site-site probability functions. We derive an exact density expansion for these intramolecular probability functions and give a diagrammatic representation of the terms contributing at each order of the expansion. The expansion is tested for a short hard-sphere chain (n=3 or 4) with site diameter sigma in a hard-sphere solvent with solvent diameter D. In comparison with Monte Carlo simulation results for 0.2< or =D/sigma< or =100, the expansion (taken to second order) is found to be quantitatively accurate for low to moderate solvent volume fractions for all size ratios. Average chain dimensions are predicted accurately up to liquidlike solvent densities. The hard-sphere chains are compressed with both increasing solvent density and decreasing solvent size. For small solvent (D<sigma), depletion effects are found and the chain structure is strongly perturbed even at low solvent volume fractions.

Kinetics of a polysoap collapse

We study the dynamics of collapse of a polysoap by means of large-scale molecular dynamics simulation and scaling arguments. A polysoap consists of a hydrophilic backbone and hydrophobic side chains attached at regular intervals along the backbone. In selective solvent conditions, the hydrophobic components aggregate, forcing the hydrophilic backbone to form loops anchored at the surface of the core, ultimately forming a micelle. The kinetics of polysoap collapse includes two major mechanisms: (1) early aggregation of the hydrophobic side chains controlled by first-order kinetics whose rate constant is given by a contact probability and (2) coalescence into larger clusters which requires activation to overcome energy barriers due to excluded volume repulsions between intermediate micelle coronas. In the late stage, the energy barrier is increasing as p(3/2), with p the number of aggregated side chains in an intermediate micelle. The corresponding late-stage rate constant decays exponentially as approximately exp(-p(3/2)).

Molecular description of the collapse of hydrophobic polymer chains in water

We propose a self-consistent molecular theory of conformational properties of flexible polymers in solution. It is applied to the collapse of a hydrophobic polymer chain in water, and can be readily generalized to any polymer-solvent system (e.g., copolymers with high complexity). We stress the potential of this method for a variety of problems, such as protein folding.