Monte Carlo simulation studies of ring polymers at athermal and theta conditions

By use of an intramolecular criterion, i.e., the direct proportionality between mean square dimension and chain length, theta conditions for linear chains and ring shaped polymers are evaluated for several types of cubic lattice chains (simple cubic, body centered cubic, and face centered cubic). The properties of the rings are evaluated for the same thermodynamic conditions under which they are prepared thus allowing for a natural amount of knots which have been identified by use of Alexander polynomials. For the limit of infinite chain lengths the same theta parameter is found for linear chains and rings. On the contrary, a significant theta point depression occurs due to an additional excluded volume effect if unknots are exclusively regarded. Parameters characteristic of the shape of rings and chains under theta conditions extrapolated to infinite chain length fairly well coincide with respective data for random walks. Mean square dimensions (characteristic of the size) of theta systems are slightly in excess as compared to nonreversal random walks due to the necessity of avoiding overlaps on a local scale. Furthermore athermal systems are studied as well for comparison; mean square dimensions are described by use of scaling relations with proper short chain corrections, shape parameters are given in the limit of infinite chain length.     

Globular structures of a helix-coil copolymer: self-consistent treatment

A self-consistent-field theory was developed in the grand canonical ensemble formulation to study transitions in a helix-coil multiblock globule. Helical and coil parts are treated as stiff rods and self-avoiding walks of variable lengths correspondingly. The resulting field theory takes, in addition to the conventional Zimm-Bragg, [J. Chem. Phys. 31, 526 (1959)] parameters, also three-dimensional interaction terms into account. The appropriate differential equations which determine the self-consistent fields were solved numerically with finite element method. Three different phase states are found: open chain, amorphous globule, and nematic liquid-crystalline (LC) globule. The LC-globule formation is driven by the interplay between the hydrophobic helical segment attraction and the anisotropic globule surface energy of an entropic nature. The full phase diagram of the helix-coil copolymer was calculated and thoroughly discussed. The suggested theory shows a clear interplay between secondary and tertiary structures in globular homopolypeptides.     

The phase diagram of a single polymer chain: New insights from a new simulation method

We present simulation results for the phase behavior of a single chain for a flexible lattice polymer model using the Wang-Landau sampling idea. Applying this new algorithm to the problem of the homopolymer collapse allows us to investigate not only the high temperature coil–globule transition but also an ensuing crystallization at lower temperature. Performing a finite size scaling analysis on the two transitions, we show that they coincide for our model in the thermodynamic limit corresponding to a direct collapse of the random coil into the crystal without intermediate coil–globule transition. As a consequence, also the many chain phase diagram of this model can be predicted to consist only of gas and crystal phase in the limit of infinite chain length. This behavior is in agreement with findings on the phase behavior of hard-sphere systems with a relatively short-ranged attractive square well. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2542–2555, 2006     

Osmotic pressure and polymer size in semidilute polymer solutions under good-solvent conditions

We consider the lattice Domb-Joyce model at a value of the coupling for which scaling corrections approximately vanish and determine the universal scaling functions associated with the osmotic pressure and the polymer size for semidilute polymer solutions (c/c( *)相似文献   

Coil-to-globule transition by dissipative particle dynamics simulation

The dynamics of a collapsing polymer under a temperature quench in dilute solution is investigated by dissipative particles dynamics. Hydrodynamic interactions and many-body interaction are preserved naturally by incorporating explicit solvent particles in this approach. Our simulation suggests a four-stage collapse pathway: localized clusters formation, cluster coarsening in situ, coarsening involving global backbone conformation change into a crumpled globule, and compaction of the globule. For all the quench depths and chain lengths used in our study, collapse proceeds without the chain getting trapped in a metastable "sausage" configuration, as reported in some earlier studies. We obtain the time scales for each of the first three stages, as well as its scaling with the quench depths ξ and chain lengths N. The total collapse time scales as τ(c) ～ ξ(-0.46 ± 0.04)N(0.98 ± 0.09), with the quench depth and degree of polymerization.     

Development of knotting during the collapse transition of polymers

A dynamic Monte Carlo simulation of the collapse transition of polymer chains is presented. The chains are represented as self-avoiding walks on the simple cubic lattice with a nearest-neighbor contact potential to model the effect of solvent quality. The knot state of the chains is determined using the knot group procedure presented in the accompanying paper. The equilibrium knot spectrum and the equilibrium rms radius of gyration as functions of the chain length and the contact potential are reported. The collapse transition was studied following quenches from good-to poor-solvent conditions. Our results confirm the prediction that the newly formed globule is not yet at equilibrium, since it has not yet achieved its equilibrium knot spectrum. For our model system, the relaxation of the knot spectrum is about an order of magnitude slower than that of the radius of gyration. The collapse transition is also studied for a model in which both ends of the chain remain in good-solvent conditions. Over the time scale of these simulations, knot formation is frustrated in this inhomogeneous model, verifying that the mechanism of knotting is the tunneling of chain ends in and out of the globule.     

Shape instabilities in adsorbed polymer condensates

A polymer chain in a poor solvent collapses to a spherical globule. If this globule is subsequently deformed, either by stretching it or if it has some non-zero net electric charge, it will become unstable to sinusoidal perturbations leading to novel final states. When such globules are imaged by direct methods such as atomic force microscopy the substrate on which they adsorb can affect the final morphological state. We therefore investigate strongly adsorbed polymers in poor solvents. We demonstrate that the real-space, Self-Consistent Field method is an ideal numerical tool in predicting equilibrium morphologies. New structures are predicted, which support previous explicit free energy calculations. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 3327–3337, 2007     

Polymers confined between two parallel plane walls

Single three-dimensional polymers confined to a slab, i.e., to the region between two parallel plane walls, are studied by Monte Carlo simulations. They are described by N-step walks on a simple cubic lattice confined to the region 1< or = z < or = D. The simulations cover both regions DRF (where RF approximately Nnu is the Flory radius, with nu approximately 0.587), as well as the cross-over region in between. Chain lengths are up to N=80 000, slab widths up to D=120. In order to test the analysis program and to check for finite size corrections, we actually studied three different models: (a) ordinary random walks (mimicking Theta polymers); (b) self-avoiding walks; and (c) Domb-Joyce walks with the self-repulsion tuned to the point where finite size corrections for free (unrestricted) chains are minimal. For the simulations we employ the pruned-enriched-Rosenbluth method with Markovian anticipation. In addition to the partition sum (which gives us a direct estimate of the forces exerted onto the walls), we measure the density profiles of monomers and of end points transverse to the slab, and the radial extent of the chain parallel to the walls. All scaling laws and some of the universal amplitude ratios are compared to theoretical predictions.     

Conformational relaxation dynamics of a poly(N-isopropylacrylamide) aqueous solution measured using the laser temperature jump transient grating method

We observed phase transition and phase relaxation processes of a poly(N-isopropylacrylamide) (PNIPAM) aqueous solution using the heterodyne transient grating (HD-TG) method combined with the laser temperature jump technique. The sample temperature was instantaneously raised by about 1.0 K after irradiation of a pump pulse to crystal violet (CV) molecules for heating, and the phase transition was induced for the sample with an initial temperature just below the lower critical solution temperature (LCST); the following phase relaxation dynamics was observed. Turbidity relaxation was observed in both the turbidity and HD-TG responses, while another relaxation process was observed only in the HD-TG response, namely via the refractive index change. It is suggested that this response is due to formation of globule molecules or their assemblies since they would have nothing to do with turbidity change but would affect the refractive index, which is dependent on the molar volume of a chemical species. Furthermore, the grating spacing dependence of the HD-TG responses suggests that the response was caused by the counter propagating diffusion of the coil molecules as a reactant species and the globule molecules as a product species and the lifetime of the globule molecules ranged from 1.5 to 5 seconds. Thus, we conclude that the turbidity reflects the dynamics of aggregate conditions, not molecular conditions. The coil and globule sizes were estimated from the obtained diffusion coefficient. The sizes of the coil molecules did not change at the initial temperatures below the LCST but increased sharply as it approaches LCST. We propose that the coil-state molecules associate due to hydrophobic interaction when the initial temperature was higher than LCST minus 0.5 K and that the globule-state molecules generated from the coil-state molecules showed a similar trend in temperature. The phase transition was also induced by heating under a microscope, and the relaxation process was followed using the fluorescence peak shift of a fluorescent molecule-labeled PNIPAM. The result also supports the existence of a globule molecule or its assembly remains for several seconds in the phase relaxation.     

Simulating the collapse transition of a two-dimensional semiflexible lattice polymer

It has been revealed by mean-field theories and computer simulations that the nature of the collapse transition of a polymer is influenced by its bending stiffness epsilon(b). In two dimensions, a recent analytical work demonstrated that the collapse transition of a partially directed lattice polymer is always first order as long as epsilon(b) is positive [H. Zhou et al., Phys. Rev. Lett. 97, 158302 (2006)]. Here we employ Monte Carlo simulation to investigate systematically the effect of bending stiffness on the static properties of a two-dimensional lattice polymer. The system's phase diagram at zero force is obtained. Depending on epsilon(b) and the temperature T, the polymer can be in one of the three phases: crystal, disordered globule, or swollen coil. The crystal-globule transition is discontinuous and the globule-coil transition is continuous. At moderate or high values of epsilon(b) the intermediate globular phase disappears and the polymer has only a discontinuous crystal-coil transition. When an external force is applied, the force-induced collapse transition will either be continuous or discontinuous, depending on whether the polymer is originally in the globular or the crystal phase at zero force. The simulation results also demonstrate an interesting scaling behavior of the polymer at the force-induced globule-coil transition.     

Conformational properties of rigid-chain amphiphilic macromolecules: The phase diagram

The coil-globule transition in rigid-chain amphiphilic macromolecules was studied by means of computer simulation, and the phase diagrams for such molecules in the solvent quality-persistence length coordinates were constructed. It was shown that the type of phase diagram depends to a substantial extent on the degree of polymerization of a macromolecule. Relatively short amphiphilic macromolecules in the poor-solvent region always form a spherical globule, with the transition to this globule involving one or two intermediate conformations. These are the disk globule if the Kuhn segment is relatively large and the string of spherical micelles or the disk globule in the case of relative flexible chains. The phase diagram of a long rodlike amphiphilic chain turned out to be even more complex. Namely, three characteristic regions were distinguished in the region of a poor solvent, depending on the chain rigidity: the region of a cylindrical globule without certain order in the main chain, the region of the cylindrical globule with blobs having the collagen ordering of the chain, and the region of coexistence of collagen-like and toroidal globules. In the intermediate transitional region, not only conformations of strings of spherical micelle beads but also the necklace conformations in which the polymer chain in each bead has collagen ordering can occur in this case.     

Diagram of State of Stiff Amphiphilic Macromolecules

Summary: We studied coil-globule transitions in stiff-chain amphiphilic macromolecules via computer modeling and constructed phase diagrams for such molecules in terms of solvent quality and persistence length. We showed that the shape of the phase diagram essentially depends on the macromolecule degree of polymerization. Relatively short amphiphilic molecules always form a spherical globule in a poor solvent, and the coil-globule transition includes one or two intermediate conformations, depending on the chain's stiffness. These are a disk-like globule in case of high enough Kuhn segment length, and a pearl necklace-like structure of spherical micelles and a disk-like globule in case of relatively flexible chains. The phase diagram of a long stiff amphiphilic chain was found to be more complex still. Thus three specific regions can be distinguished in the poor solvent region, depending on the chain stiffness. These correspond to a cylindrical globule without any specific backbone ordering, a cylindrical globule containing blobs with collagen-like ordering of the chain, and co-existence of collagen-like and toroidal globules. In the intermediate transition region in this case, apart from the pearl necklace-like conformations with spherical micelles, necklace conformations can be also observed where the polymeric chain has collagen-like ordering within each bead.     

激光光散射研究聚(N-异丙基丙烯酰胺)单链及其智能凝胶微球在水中的相变(下)

（接上期）２聚（Ｎ－异丙基丙烯酰胺）微凝胶在水中的体积相变２．１理论部分凝胶体积相变热力学：聚合物凝胶的溶胀和蜷缩可以用膨胀因子α＝（Ｖ／Ｖ０）１／３＝（ΦＴ／ΦΘ）１／３来表征，其中ΦΘ的ΦＴ分别是温度Θ和Ｔ下凝胶网络的体积分数。在平均场理论中，中...     

Microphase Separation within a Comb Copolymer with Attractive Side Chains: A Computer Simulation Study

Computer simulation modelling of a flexible comb copolymer with attractive interactions between the monomer units of the side chains is performed. The conditions for the coil‐globule transition, induced by the increase of attractive interaction, ε, between side chain monomer units, are analysed for different values of the number of monomer units in the backbone, N, in the side chains, n, and between successive grafting points, m. It is shown that the coil‐globule transition of such a copolymer corresponds to a first‐order phase transition. The energy of attraction (ε) required for the realisation of the coil‐globule transition decreases with increasing n and decreasing m. The coil‐globule transition is accompanied by significant aggregation of side chain units. The resulting globule has a complex structure. In the case of a relatively short backbone (small value of N), the globule consists of a spherical core formed by side chains and an enveloping shell formed by the monomer units of the backbone. In the case of long copolymers (large value of N), the side chains form several spherical micelles while the backbone is wrapped on the surfaces of these micelles and between them.     

Effect of correlations in the interaction along polymer chain on the globule structure

A special potential for interaction between polymer chain units, whose energy decreases with increasing distance s between the units as s ^–1, was introduced for the first time. According to Monte Carlo simulation, interactions of this type result in the formation of a globule with an equilibrium packing of domains in space. The radius of gyration of a chain segment in these globules varies with segment length according to the scaling law typical of crumpled globules.     

Monte Carlo Computer Simulation of a Single Semi-Flexible Macromolecule at a Plane Surface

Summary: The properties of a single semiflexible mushroom chain at a plane surface with a long-ranged attracting potential are studied by means of lattice Monte Carlo computer simulation using the bond fluctuation model, configurational bias algorithm for chain re-growing and the Wang-Landau sampling technique. We present the diagram of states in variables temperature T vs. strength of the adsorption potential, ε_w, for a quite short semiflexible chain consisting of N = 64 monomer units. The diagram of states consists of the regions of a coil, liquid globule, solid isotropic globule, adsorbed coil and cylinder-like liquid-crystalline globule. At low values of the adsorption strength ε_w the coil–globule and the subsequent liquid–solid globule transitions are observed upon decreasing temperature below the adsorption transition point. At high values of ε_w these two transitions change into a single transition from an adsorbed coil to a cylinder-like liquid-crystalline solid globule. We conclude that for a semiflexible chain the presence of a plane attracting surface favors the formation of a globule with internal liquid-crystalline ordering of bonds.     

Adsorption and scattering of H2 and D2 by NiAl(110)

We present quasiclassical dynamics calculations of H2 and D2 scattering by the NiAl(110) surface using a recently proposed six-dimensional potential-energy surface (PES) obtained from density-functional theory calculations. The results for dissociative adsorption confirm several experimental predictions using (rotationally hot) D2 beams, namely, the existence of a dissociation barrier, the small isotopic effect, the importance of vibrational enhancement, and the existence of normal energy scaling. The latter conclusion shows that normal energy scaling is not necessarily associated with weak corrugated surfaces. The results for rotationally elastic and inelastic diffractions are also in reasonable agreement with experiment, but they show that many more diffractive transitions are responsible for the observed structures than previously assumed. This points to the validity of the PES recently proposed [P. Riviere, H. F. Busnengo, and F. Martin, J. Chem. Phys. 121, 751 (2004)] to describe dissociative adsorption as well as rotationally elastic and inelastic diffractions in the H2NiAl(110) system.     

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.     

Dynamics of collapsed polymers under the simultaneous influence of elongational and shear flows

Collapsed polymers in solution represent an oft-overlooked area of polymer physics, however recent studies of biopolymers in the bloodstream have suggested that the physics of polymer globules are not only relevant but could potentially lead to powerful new ways to manipulate single molecules using fluid flows. In the present article, we investigate the behavior of a collapsed polymer globule under the influence of linear combinations of shear and elongational flows. We generalize the theory of globule-stretch transitions that has been developed for the specific case of simple shear and elongational flows to account for behavior in arbitrary flow fields. In particular, we find that the behavior of a globule in flow is well represented by a two-state model wherein the critical parameters are the transition probabilities to go from a collapsed to a stretched state P(g-s) and vice versa P(s-g). The collapsed globule to stretch transition is described using a nucleation protrusion mechanism, and the reverse transition is described using either a tumbling or a relaxation mechanism. The magnitudes of P(g-s) and P(s-g) govern the state in which the polymer resides; for P(g-s) ≈ 0 and P(s-g) ≈ 1 the polymer is always collapsed, for P(g-s) ≈ 0 and P(s-g) ≈ 0 the polymer is stuck in either the collapsed or stretched state, for P(g-s) ≈ 1 and P(s-g) ≈ 0 the polymer is always stretched, and for P(g-s) ≈ 1 and P(s-g) ≈ 1 the polymer undergoes tumbling behavior. These transition probabilities are functions of the flow geometry, and we demonstrate that our theory quantitatively predicts globular polymer conformation in the case of mixed two-dimensional flows, regardless of orientation and representation, by comparing theoretical results to Brownian dynamics simulations. Generalization of the theory to arbitrary three-dimensional flows is discussed as is the incorporation of this theory into rheological equations.