Slow relaxation and solvent effects in the collapse of a polymer

Investigation of Polymer Aging Mechanisms Using Molecular Simulations: A Review

Aging has a serious impact on the properties of functional polymers. Therefore, it is necessary to study the aging mechanism to prolong the service and storage life of polymer-based devices and materials. Due to the limitations of traditional experimental methods, more and more studies have adopted molecular simulations to analyze the intrinsic mechanisms of aging. In this paper, recent advances in molecular simulations of the aging of polymers and their composites are reviewed. The characteristics and applications of commonly used simulation methods in the study of the aging mechanisms (traditional molecular dynamics simulation, quantum mechanics, and reactive molecular dynamics simulation) are outlined. The current simulation research progress of physical aging, aging under mechanical stress, thermal aging, hydrothermal aging, thermo-oxidative aging, electric aging, aging under high-energy particle impact, and radiation aging is introduced in detail. Finally, the current research status of the aging simulations of polymers and their composites is summarized, and the future development trend has been prospected.

Study of the conformation of polyelectrolyte aggregates using coarse-grained molecular dynamics simulations

The conformation of polyelectrolyte aggregates as a function of the backbone rigidity is investigated by coarse-grained molecular dynamics simulation. The polyelectrolyte is represented by a bead-spring chain with charged side chains. The simulations start from the uniform distributions of the polyelectrolytes, and the resultant polyelectrolyte conformation after a few microseconds exhibits spherical self-aggregates, clusters, or bending bundle-like aggregates, depending on the backbone rigidity. The interaggregate structures on a large scale are featured by the static structure factor (SSF). The simulated SSFs of the bending bundle-like aggregates are consistent with those of the small angle X-ray scattering (SAXS) measurement so we successfully assign the microscopic structures of polyelectrolytes to the SAXS measurement. The power-law of the SSFs for the bundle conditions is steeper than that of the conventional cylinder model. The present study finds that such discrepancy in the power-law results from the bending of the bundle-like aggregates. In addition, the relaxation behavior includes slow dynamics. The present study proposes that such slow dynamics results from diffusion-limited aggregation and from gliding processes to reduce local metastable folding within the aggregates.

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.

Methods for Monte Carlo simulations of biomacromolecules

The state-of-the-art for Monte Carlo (MC) simulations of biomacromolecules is reviewed. Available methodologies for sampling conformational equilibria and associations of biomacromolecules in the canonical ensemble, given a continuum description of the solvent environment, are reviewed. Detailed sections are provided dealing with the choice of degrees of freedom, the efficiencies of MC algorithms and algorithmic peculiarities, as well as the optimization of simple movesets. The issue of introducing correlations into elementary MC moves, and the applicability of such methods to simulations of biomacromolecules is discussed. A brief discussion of multicanonical methods and an overview of recent simulation work highlighting the potential of MC methods are also provided. It is argued that MC simulations, while underutilized biomacromolecular simulation community, hold promise for simulations of complex systems and phenomena that span multiple length scales, especially when used in conjunction with implicit solvation models or other coarse graining strategies.

Relaxation of ultralarge VWF bundles in a microfluidic-AFM hybrid reactor

The crucial role of the biopolymer "Von Willebrand factor" (VWF) in blood platelet binding is tightly regulated by the shear forces to which the protein is exposed in the blood flow. Under high-shear conditions, VWFs ability to immobilize blood platelets is strongly increased due to a change in conformation which at sufficient concentration is accompanied by the formation of ultra large VWF bundles (ULVWF). However, little is known about the dynamic and mechanical properties of such bundles. Combining a surface acoustic wave (SAW) based microfluidic reactor with an atomic force microscope (AFM) we were able to study the relaxation of stretched VWF bundles formed by hydrodynamic stress. We found that the dynamical response of the network is well characterized by stretched exponentials, indicating that the relaxation process proceeds through hopping events between a multitude of minima. This finding is in accordance with current ideas of VWF self-association. The longest relaxation time does not show a clear dependence on the length of the bundle, and is dominated by the internal conformations and effective friction within the bundle.

Kinetically driven helix formation during the homopolymer collapse process

Using Langevin simulations, we find that simple "generic" bead-and-spring homopolymer chains in a sufficiently bad solvent spontaneously develop helical order during the process of collapsing from an initially stretched conformation. The helix formation is initiated by the unstable modes of the straight chain, which drive the system towards a long-lived metastable transient state. The effect is most pronounced if hydrodynamic interactions are screened.

Stepwise unfolding of collapsed polymers

Motivated by recent experimental data on DNA stretching in presence of polyvalent counterions, we study the force-induced unfolding of a homopolymer on and off lattice. In the fixed force ensemble the globule unravels via a series of steps due to surface effects which play an important role for finite-size chains. This holds both for flexible and stiff polymers. We discuss in a qualitative way how this result may impact on the interpretation of DNA stretching experiments showing peaks in the characteristic curves, by extracting from the raw data the corresponding elongation- versus-force characteristic curves. Furthermore, approximate analytical and numerical calculations, valid in a quasi-equilibrium fixed stretch ensemble, and if the initial low-temperature state is ordered in a spool, show that the average force versus elongation displays peaks related to the geometry of the initial configuration. We finally argue how the proposed mechanisms identified for the arising of peaks may couple in the experiments, and comment on the role of dynamic effects.

Stretching an adsorbed polymer globule

Using molecular dynamic simulation, we study the stretching of an adsorbed homopolymer in a poor solvent with one end held at a distance ze from the substrate. We measure the vertical force f on the end of the chain as a function of the extension ze and the substrate interaction energy w. The force reaches a plateau value at large extensions. In the strong adsorption limit, we show that the plateau value increases linearly in w in good agreement with a theoretical model. In the weak adsorption limit, a polymer globule with a layered structure is formed and elastically deformed when stretched. In both cases a simple theoretical model permits us to predict the relation between the necessary force to fully detach the polymer and its critical extension.