Finite chain length effects on the coil–globule transition of stiff-chain macromolecules: A Monte Carlo simulation

State Transitions and Crystalline Structures of Single Polyethylene Rings: MD Simulations

This study investigates the structural changes of cyclic polyethylene (PE) single chains during cooling through molecular dynamics simulations. The influence of topological constraint on a ring is examined by comparing it with the results of its linear counterpart. A pseudo phase diagram of state transition for PE rings based on length and temperature is constructed, revealing a consistent chain-folding transition during cooling. The shape anisotropy of short crystallized cyclic chains exhibits oscillations with chain length, leading to a more pronounced odd-even effect in single cyclic chains compared with the linear ones. A honeycomb model is proposed to elucidate the odd-even effect of chain folding in crystalline structures of single linear and cyclic chains, and we discuss its potential to predict surface tension. Analyses of the tight folding model and the re-entry modes demonstrate that a cyclic chain possesses a shorter average crystalline stem length and a more compact folded structure than its linear counterpart. The findings highlight the impact of topological change on crystallization and the odd-even effect of chain length, providing valuable insights for understanding polymer crystallization with different topologies.

Modulation of flexo-rigid balance in photoresponsive thymine grafted copolymers towards designing smart healable coating

Efficacy and durability of the photovoltaic device mandates its protection against hot, humid weather condition, high energy of UV light and unwanted scratches. Such challenges can be mitigated by smart polymeric coating with inherent properties e.g. hydrophobicity to prevent moisture, optimal viscocity for better processibility and crack-healing. The hydrophobic polymers TP1-TP4 containing pendant photo-crosslinkable thymine moieties are designed that undergo [2 + 2] photocycloaddition upon UVB irradiation and can be dynamically reverted back upon irradiation with UVC light. A judicious control of solvent environment, chain length, functionality% and concentration of the polymers regulate the aspects of photodimerization thereby, rendering intra or inter-chain collapse to form diverse nanostructures. Photodimerization of the thymine moieties renders coil to globule transformation in dilute condition whereas irradiation performed at high macromolecular concentration regime exhibits higher order nanostructures. The photoresponsive chain collapse leads to the formation of rigid crosslinked domains within flexible polymer chains akin to the hard-soft phases of thermoplastic elastomers. Such rigidification of the crosslinked segments endows a tool to photomodulate the glass transition temperature (T g) that can dynamically revert back upon decrosslinking. Further, the structural modulation of the polymers is explored towards autonomic and nonautonomic self-healing behaviour at ambient conditions. Moreover, the self-healing efficacy can be tuned with the film thickness and it remains unaltered upon using solar simulator or direct sunlight. Overall, such hydrophobic low T g polymers display photo-regulated self-healing mechanism consisting of both autonomic and non-autonomic self-healing and may find applications in designing smart protective coatings for photovoltaic devices.

Collapse transition of a heterogeneous polymer in a crowded medium

Long chain molecules can be entropically compacted in a crowded medium. We study the compaction transition of a heterogeneous polymer with ring topology by crowding effects in a free or confined space. For this, we use molecular dynamics simulations in which the effects of crowders are taken into account through effective interactions between chain segments. Our parameter choices are inspired by the Escherichia coli chromosome. The polymer consists of small and big monomers; the big monomers dispersed along the backbone are to mimic the binding of RNA polymerases. Our results show that the compaction transition is a two-step process: initial compaction induced by the association (clustering) of big monomers followed by a gradual overall compaction. They also indicate that cylindrical confinement makes the initial transition more effective; for representative parameter choices, the initial compaction accounts for about 60% reduction in the chain size. Our simulation results support the view that crowding promotes clustering of active transcription units into transcription factories.

Diagrams of States of Single Flexible-Semiflexible Multi-Block Copolymer Chains: A Flat-Histogram Monte Carlo Study

The combination of flexibility and semiflexibility in a single molecule is a powerful design principle both in nature and in materials science. We present results on the conformational behavior of a single multiblock-copolymer chain, consisting of equal amounts of Flexible (F) and Semiflexible (S) blocks with different affinity to an implicit solvent. We consider a manifold of macrostates defined by two terms in the total energy: intermonomer interaction energy and stiffness energy. To obtain diagrams of states (pseudo-phase diagrams), we performed flat-histogram Monte Carlo simulations using the Stochastic Approximation Monte Carlo algorithm (SAMC). We have accumulated two-Dimensional Density of States (2D DoS) functions (defined on the 2D manifold of macrostates) for a SF-multiblock-copolymer chain of length N=64 with block lengths b = 4, 8, 16, and 32 in two different selective solvents. In an analysis of the canonical ensemble, we calculated the heat capacity and determined its maxima and the most probable morphologies in different regions of the state diagrams. These are rich in various, non-trivial morphologies, which are formed without any specific interactions, and depend on the block length and the type of solvent selectivity (preferring S or F blocks, respectively). We compared the diagrams with those for the non-selective solvent and reveal essential changes in some cases. Additionally, we implemented microcanonical analysis in the “conformational” microcanonical (NVU, where U is the potential energy) and the true microcanonical (NVE, where E is the total energy) ensembles with the aim to reveal and classify pseudo-phase transitions, occurring under the change of temperature.

The Persistence Length of Semiflexible Polymers in Lattice Monte Carlo Simulations

While applying computer simulations to study semiflexible polymers, it is a primary task to determine the persistence length that characterizes the chain stiffness. One frequently asked question concerns the relationship between persistence length and the bending constant of applied bending potential. In this paper, theoretical persistence lengths of polymers with two different bending potentials were analyzed and examined by using lattice Monte Carlo simulations. We found that the persistence length was consistent with theoretical predictions only in bond fluctuation model with cosine squared angle potential. The reason for this is that the theoretical persistence length is calculated according to a continuous bond angle, which is discrete in lattice simulations. In lattice simulations, the theoretical persistence length is larger than that in continuous simulations.

Tuning the collapse transition of weakly charged polymers by ion-specific screening and adsorption

The experimentally observed swelling and collapse response of weakly charged polymers to the addition of specific salts displays quite convoluted behavior that is not easy to categorize. Here we use a minimalistic implicit-solvent/explicit-salt simulation model with a focus on ion-specific interactions between ions and a single weakly charged polyelectrolyte to qualitatively explain the observed effects. In particular, we demonstrate ion-specific screening and bridging effects cause collapse at low salt concentrations whereas the same strong ion-specific direct interactions drive re-entrant swelling at high concentrations. Consistently with experiments, a distinct salt concentration at which the salting-out power of anions inverts from the reverse to direct Hofmeister series is observed. At this so called isospheric point, the ion-specific effects vanish. Furthermore, with additional simplifying assumptions, an ion-specific mean-field model is developed for the collapse transition which quantitatively agrees with the simulations. Our work demonstrates the sensitivity of the structural behavior of charged polymers to the addition of specific salt beyond simple screening and shall be useful for further guidance of experiments.

The folding pathways and thermodynamics of semiflexible polymers

Inspired by the protein folding and DNA packing, we have systematically studied the thermodynamic and kinetic behaviors of single semiflexible homopolymers by Langevin dynamics simulations. In line with experiments, a rich variety of folding products, such as rod-like bundles, hairpins, toroids, and a mixture of them, are observed in the complete diagram of states. Moreover, knotted structures with a significant population are found in a certain range of bending stiffness in thermal equilibrium. As the solvent quality becomes poorer, the population of the intermediate occurring in the folding process increases, which leads to a severe chevron rollover for the folding arm. However, the population of the intermediates in the unfolding process is very low, insufficient to induce unfolding arm rollover. The total types of folding pathways from the coil state to the toroidal state for a semiflexible polymer chain remain unchanged by varying the solvent quality or temperature, whereas the kinetic partitioning into different folding events can be tuned significantly. In the process of knotting, three types of mechanisms, namely, plugging, slipknotting, and sliding, are discovered. Along the folding evolution, a semiflexible homopolymer chain can knot at any stage of folding upon leaving the extended coil state, and the probability to find a knot increases with chain compactness. In addition, we find rich types of knotted topologies during the folding of a semiflexible homopolymer chain. This study should be helpful in gaining insight into the general principles of biopolymer folding.

Self-assembly of bis-salphen compounds: from semiflexible chains to webs of nanorings

The recently-observed self-assembly of certain salphen-based compounds into neuron-like networks of microrings interconnected with nano-thin strings may suggest a new highly-potent tool for nanoscale patterning. However, the mechanism behind such phenomena needs to be clarified before they can be applied in materials design. Here we show that, in contrast with what was initially presumed, the emergence of a "rings-and-rods" pattern is unlikely to be explained by merging, collapse and piercing of vesicles as in previously reported cases of nanorings self-assembly via non-bonding interactions. We propose an alternative explanation: the compounds under study form a 1D coordination polymer, the fibres of which are elastic enough to fold into toroidal globules upon solvent evaporation, while being able to link separate chains into extended networks. This becomes possible because the structure of the compound's scaffold is found to adopt a very different conformation from that inferred in the original work. Based on ab initio and molecular dynamics calculations we propose a step-by-step description of self-assembly process of a supramolecular structure which explains all the observed phenomena in a simple and clear way. The individual roles of the compound' s scaffold structure, coordination centres, functional groups and solvent effects are also explained, opening a route to control the morphology of self-assembled networks and to synthesize new compounds exhibiting similar behaviour.

Self-Assembly of Lecithin and Bile Salt in the Presence of Inorganic Salt in Water: Mesoscale Computer Simulation

The influence of inorganic salt on the structure of lecithin/bile salt mixtures in aqueous solution is studied by means of dissipative particle dynamics simulations. We propose a coarse-grained model of phosphatidylcholine and two types of bile salts (sodium cholate and sodium deoxycholate) and also take into account the presence of low molecular weight salt. This model allows us to study the system on rather large time and length scales (up to about ∼20 μs and 50 nm) and to reveal mechanisms of experimentally observed increasing viscosity upon increasing the low molecular weight salt concentration in this system. We show that increasing the low molecular weight salt concentration induces the growth of cylinder-like micelles formed in lecithin/bile salt mixtures in water. These wormlike micelles can entangle into transient networks displaying perceptible viscoelastic properties. Computer simulation results are in good qualitative agreement with experimental observations.

On the Pseudo Phase Diagram of Single Semi-Flexible Polymer Chains: A Flat-Histogram Monte Carlo Study

Local stiffness of polymer chains is instrumental in all structure formation processes of polymers, from crystallization of synthetic polymers to protein folding and DNA compactification. We present Stochastic Approximation Monte Carlo simulations-a type of flat-histogram Monte Carlo method-determining the density of states of a model class of single semi-flexible polymer chains, and, from this, their complete thermodynamic behavior. The chains possess a rich pseudo phase diagram as a function of stiffness and temperature, displaying non-trivial ground-state morphologies. This pseudo phase diagram also depends on chain length. Differences to existing pseudo phase diagrams of semi-flexible chains in the literature emphasize the fact that the mechanism of stiffness creation matters.

Dilute Semiflexible Polymers with Attraction: Collapse, Folding and Aggregation

We review the current state on the thermodynamic behavior and structural phases of self- and mutually-attractive dilute semiflexible polymers that undergo temperature-driven transitions. In extreme dilution, polymers may be considered isolated, and this single polymer undergoes a collapse or folding transition depending on the internal structure. This may go as far as to stable knot phases. Adding polymers results in aggregation, where structural motifs again depend on the internal structure. We discuss in detail the effect of semiflexibility on the collapse and aggregation transition and provide perspectives for interesting future investigations.

Diagram of states and morphologies of flexible-semiflexible copolymer chains: A Monte Carlo simulation

A single copolymer chain consisting of multiple flexible (F) and semiflexible (S) blocks has been studied using a continuum bead-spring model by Stochastic Approximation Monte Carlo simulations, which determine the density of states of the model. The only difference between F and S blocks is the intramolecular bending potential, all non-bonded interactions are equal. The state diagrams for this class of models display multiple nematic phases in the collapsed state, characterized through a demixing of the blocks of different stiffness and orientational ordering of the stiff blocks. We observe dumbbell-like morphologies, lamellar phases, and for the larger block lengths also Saturn-like structures with a core of flexible segments and the stiff segments forming a ring around the core.

Knots as a Topological Order Parameter for Semiflexible Polymers

Using a combination of the multicanonical Monte Carlo algorithm and the replica-exchange method, we investigate the influence of bending stiffness on the conformational phases of a bead-stick homopolymer model and present the pseudophase diagram for the complete range of semiflexible polymers, from flexible to stiff. Although it is a simple model, we observe a rich variety of conformational phases, reminiscent of conformations observed for synthetic polymers or biopolymers. Depending on the bending stiffness, the model exhibits different pseudophases like bent, hairpin, or toroidal. In particular, we find thermodynamically stable knots and unusual transitions into these knotted phases with a clear phase coexistence, but almost constant mean total energy, and hence almost no latent heat.

Thermodynamics and structure of macromolecules from flat-histogram Monte Carlo simulations

Over the last decade flat-histogram Monte Carlo simulations, especially multi-canonical and Wang-Landau simulations, have emerged as a strong tool to study the statistical mechanics of polymer chains. These investigations have focused on coarse-grained models of polymers on the lattice and in the continuum. Phase diagrams of chains in bulk as well as chains attached to surfaces were studied, for homopolymers as well as for protein-like models. Also, aggregation behavior in solution of these models has been investigated. We will present here the theoretical background for these simulations, explain the algorithms used and discuss their performance and give an overview over the systems studied with these methods in the literature, where we will limit ourselves to studies of coarse-grained model systems. Implementations of these algorithms on parallel computers will be also briefly described. In parallel to the development of these simulation methods, the power of a micro-canonical analysis of such simulations has been recognized, and we present the current state of the art in applying the micro-canonical analysis to phase transitions in nanoscopic polymer systems.

Thermodynamic description of Hofmeister effects on the LCST of thermosensitive polymers

Cosolvent effects on protein or polymer collapse transitions are typically discussed in terms of a two-state free energy change that is strictly linear in cosolute concentration. Here we investigate in detail the nonlinear thermodynamic changes of the collapse transition occurring at the lower critical solution temperature (LCST) of the role-model polymer poly(N-isopropylacrylamide) [PNIPAM] induced by Hofmeister salts. First, we establish an equation, based on the second-order expansion of the two-state free energy in concentration and temperature space, which excellently fits the experimental LCST curves and enables us to directly extract the corresponding thermodynamic parameters. Linear free energy changes, grounded on generic excluded-volume mechanisms, are indeed found for strongly hydrated kosmotropes. In contrast, for weakly hydrated chaotropes, we find significant nonlinear changes related to higher order thermodynamic derivatives of the preferential interaction parameter between salts and polymer. The observed non-monotonic behavior of the LCST can then be understood from a not yet recognized sign change of the preferential interaction parameter with salt concentration. Finally, we find that solute partitioning models can possibly predict the linear free energy changes for the kosmotropes, but fail for chaotropes. Our findings cast strong doubt on their general applicability to protein unfolding transitions induced by chaotropes.

The Fundamental Role of Flexibility on the Strength of Molecular Binding

Non-covalent molecular association underlies a diverse set of biologically and technologically relevant phenomena, including the action of drugs on their biomolecular targets and self- and supra-molecular assembly processes. Computer models employed to model binding frequently use interaction potentials with atomistic detail while neglecting the thermal molecular motions of the binding species. However, errors introduced by this simplification and, more broadly, the thermodynamic consequences of molecular flexibility on binding, are little understood. Here, we isolate the fundamental relationship of molecular flexibility to binding thermodynamics via simulations of simplified molecules with a wide range of flexibilities but the same interaction potential. Disregarding molecular motion is found to generate large errors in binding entropy, enthalpy and free energy, even for molecules that are nearly rigid. Indeed, small decreases in rigidity markedly reduce affinity for highly rigid molecules. Remarkably, precisely the opposite occurs for more flexible molecules, for which increasing flexibility leads to stronger binding affinity. We also find that differences in flexibility suffice to generate binding specificity: for example, a planar surface selectively binds rigid over flexible molecules. Intriguingly, varying molecular flexibility while keeping interaction potentials constant leads to near-linear enthalpy-entropy compensation over a wide range of flexibilities, with the unexpected twist that increasing flexibility produces opposite changes in entropy and enthalpy for molecules in the flexible versus the rigid regime. Molecular flexibility is thus a crucial determinant of binding affinity and specificity and variations in flexibility can lead to strong yet non-intuitive consequences.

A Replica Exchange Molecular Dynamics Simulation of a Single Polyethylene Chain: Temperature Dependence of Structural Properties and Chain Conformational Study at the Equilibrium Melting Temperature

The conformational properties of a finite length polyethylene chain were explored over a wide range of temperatures using a replica exchange molecular dynamics simulation providing high quality simulation data representative for the equilibrium behavior of the chain molecule. The radial distribution function (RDF) and the structure factor S(q) of the chain as a function of temperature are analyzed in detail. The different characteristic peaks in the RDF and S(q) were assigned to specific distances in the chain and structural changes occurring with the temperature. In S(q), a peak characteristic for the order in the solid state was found and used to determine the equilibrium melting temperature. A detailed scaling analysis of the structure factor covering the full q range was performed according to the work of Hammouda. In the Θ region, a quantitative analysis of the full structure factor was done using the equivalent Kuhn chain, which enabled us to assign the Θ region of our chain and to demonstrate, in our particular case, the failure of the Gaussian chain approach. The chain conformational properties at the equilibrium melting temperature are discussed using conformational distribution functions, using the largest principal component of the radius of gyration and shape parameters as order parameters. We demonstrate that for the system studied here, the Landau free energy expression based on this conformational distribution information leads to erroneous conclusions concerning the thermodynamic transition behavior. Finally, we focus on the instantaneous conformational properties at the equilibrium melting temperature and give a detailed analysis of the conformational shapes using different shape parameters and a simulation snapshot. We show that the chain does not only take the lamellar rod-like and globular conformational shapes, typical of the solid and liquid states, but can also explore many other conformational states, including the toroidal conformational state. It is the first demonstration that a flexible molecule like PE can also take a toroidal conformational state, which is normally linked to stiffer chains.

Nonideal behavior of the intramolecular structure factor of dilute polymers in a theta solvent

We study the configurational properties of single polymers in a theta solvent by Monte Carlo simulation of the bond fluctuation model. The intramolecular structure factor at the theta point is found to be distinctively different from that of the ideal chain. The structure factor shows a hump around q approximately 5/Rg and a dip around q approximately 10/Rg in the Kratky plot with Rg being the radius of gyration. This feature is apparently similar to that in a melt. The theoretical expression by the simple perturbation expansion to the first order in terms of the Mayer function can be fitted to the obtained structure factor quite well, but the second virial coefficient cannot be set to zero.

The phase behavior of polyethylene ring chains

The equilibrium properties of an isolated polyethylene ring chain are studied by using molecular dynamics (MD) simulations. The results of an 80-bond linear chain are also presented, which are in agreement with previous studies of square-well chains and Lennard-Jones (LJ) homopolymers. Mainly, we focus on the collapse of polyethylene ring chains. At high temperatures, a fully oblate structure is observed for the ring chains with different chain lengths. For such an oblate structure, a shape factor of delta(*)=0.25 and a rodlike scaling relation between the radius of gyration and chain lengths could be deduced easily in theory, and the same results are obtained by our MD simulations. Such an oblate structure can be obtained by Monte Carlo simulation only for sufficient stiff ring chains. When the temperature decreases, an internal energy barrier is observed. This induces a strong peak in the heat capacity, denoting a gas-liquid-like transition. This energy barrier comes mainly from the local monomer-monomer interactions, i.e., the bond-stretching, the bond-bending, and the torsion potentials. A low temperature peak is also observed in the same heat capacity curve, representing a liquid-solid-like transition. These numerical simulation results support a two-stage collapse of polyethylene ring chains; however, the nature should be different from the square-well and LJ ring chains.

Equation of state for macromolecules of variable flexibility in good solvents: a comparison of techniques for Monte Carlo simulations of lattice models

The osmotic equation of state for the athermal bond fluctuation model on the simple cubic lattice is obtained from extensive Monte Carlo simulations. For short macromolecules (chain length N=20 ) we study the influence of various choices for the chain stiffness on the equation of state. Three techniques are applied and compared in order to critically assess their efficiency and accuracy: the "repulsive wall" method, the thermodynamic integration method (which rests on the feasibility of simulations in the grand canonical ensemble), and the recently advocated sedimentation equilibrium method, which records the density profile in an external (e.g., gravitationlike) field and infers, via a local density approximation, the equation of state from the hydrostatic equilibrium condition. We confirm the conclusion that the latter technique is far more efficient than the repulsive wall method, but we find that the thermodynamic integration method is similarly efficient as the sedimentation equilibrium method. For very stiff chains the onset of nematic order enforces the formation of an isotropic-nematic interface in the sedimentation equilibrium method leading to strong rounding effects and deviations from the true equation of state in the transition regime.

Langevin dynamics of semiflexible polyelectrolytes: rod-toroid-globule-coil structures and counterion distribution

We have investigated the nature of counterion condensation on uniformly charged semiflexible polyelectrolyte chains and the concomitant configurations by monitoring the role of chain stiffness, chain length, counterion valency, and the strength of electrostatic interaction. The counterion condensation is seen to follow the adsorption process and the effective polymer charge increases with chain stiffness. Size and shape, as calculated through the radius of gyration, effective persistence length, and hydrodynamic radius, are studied. Stable coil-like, globular, folded-chain, toroidal, and rodlike configurations are possible at suitable combinations of values of chain stiffness, chain length, electrostatic interaction strength, and the valency of counterion. For high strengths of electrostatic interactions, sufficiently stiff polyelectrolytes form toroids in the presence of multivalent counterions, whereas flexible polyelectrolytes form disordered globules. The kinetic features of the nucleation and growth of toroids are monitored. Several metastable structures are found to frustrate the formation of toroids. The generic pathway involves the nucleation of one primary loop somewhere along the chain contour, followed by a growth process where the rest of the chain is folded continuously on top of the primary loop. The dependence of the average radii of toroids on the chain length is found to be roughly linear, in disagreement with existing scaling arguments.