Solutions Of Semiflexible Polymers Research Articles

Cluster Formation in Solutions of Polyelectrolyte Rings.

We use molecular dynamics simulations to explore concentrated solutions of semiflexible polyelectrolyte ring polymers, akin to the DNA mini-circles, with counterions of different valences. We find that the assembly of rings into nanoscopic cylindrical stacks is a generic feature of the systems, but the morphology and dynamics of such a cluster can be steered by the counterion conditions. In general, a small addition of trivalent ions can stabilize the emergence of clusters due to the counterion condensation, which mitigates the repulsion between the like-charged rings. Stoichiometric addition of trivalent ions can even lead to phase separation of the polyelectrolyte ring phase due to the ion-bridging effects promoting otherwise entropically driven clustering. On the other hand, monovalent counterions cause the formation of stacks to be re-entrant with density. The clusters are stable within a certain window of concentration, while above the window the polyelectrolytes undergo an osmotic collapse, disfavoring ordering. The cluster phase exhibits characteristic cluster glass dynamics with arrest of collective degrees of freedom but not the self-ones. On the other hand, the collapsed phase shows arrest on both the collective and single level, suggesting an incipient glass-to-glass transition, from a cluster glass of ring clusters to a simple glass of rings.

Orientational Wetting and Topological Transitions in Confined Solutions of Semiflexible Polymers

Despite their considerable practical and biological applications, the link between molecular properties, assembly conditions and self-organized structure in confined polymer solutions remains elusive. Here, we explore the lyotropic nematic ordering of semi-flexible chains in spherical confinement for multiple contour lengths across a wide regime of concentrations. We uncover an original crossover from two distinct quadrupolar states, both characterized by regular tetrahedral patterns of surface topological defects, to either longitudinal, latitudinal or spontaneously-twisted bipolar structures with increasing densities. These transitions, along with the intermediary arrangements that they involve, are attributed to the combination of orientational wetting with subtle variations in their liquid-crystal (LC) elastic anisotropies. Our molecular simulations are corroborated by density functional calculations, and are quantified through the introduction of several order parameters as well as an unsupervised learning scheme for the localization of topological defects. Our results agree quantitatively with the predictions of continuum nematic elasticity theories, and evidence the extent to which the folding of macromolecules and the self-assembly of low-molecular-weight LCs may be guided by the same, universal principles.

Semiflexible polymer solutions. II. Fluctuations and Frank elastic constants.

We study the collective elastic behavior of semiflexible polymer solutions in a nematic liquid-crystalline state using polymer field theory. Our polymer field-theoretic model of semiflexible polymer solutions is extended to include second-order fluctuation corrections to the free energy, permitting the evaluation of the Frank elastic constants based on orientational order fluctuations in the nematic state. Our exact treatment of wormlike chain statistics permits the evaluation of behavior from the nematic state, thus accurately capturing the impact of single-chain behavior on collective elastic response. Results for the Frank elastic constants are presented as a function of aligning field strength and chain length, and we explore the impact of conformation fluctuations and hairpin defects on the twist, splay, and bend moduli. Our results indicate that the twist elastic constant Ktwist is smaller than both bend and splay constants (Kbend and Ksplay, respectively) for the entire range of polymer rigidity. Splay and bend elastic constants exhibit regimes of dominance over the range of chain stiffness, where Ksplay > Kbend for flexible polymers (large-N limit) while the opposite is true for rigid polymers. Theoretical analysis also suggests the splay modulus tracks exactly to that of the end-to-end distance in the transverse direction for semiflexible polymers at intermediate to large-N. These results provide insight into the role of conformation fluctuations and hairpin defects on the collective response of polymer solutions.

Wetting transitions of polymer solutions: Effects of chain length and chain stiffness.

Wetting and drying phenomena are studied for flexible and semiflexible polymer solutions via coarse-grained molecular dynamics simulations and density functional theory calculations. This study is based on the use of Young's equation for the contact angle, determining all relevant surface tensions from the anisotropy of the pressure tensor. The solvent quality (or effective temperature, equivalently) is varied systematically, while all other interactions remain unaltered. For flexible polymers, the wetting transition temperature Tw increases monotonically with chain length N, while the contact angle at temperatures far below Tw is independent of N. For semiflexible polymer solutions, Tw varies non-monotonically with the persistence length: Initially, Tw increases with increasing chain stiffness and reaches a maximum, but then a sudden drop of Tw is observed, which is associated with the isotropic-nematic transition of the system.

The depletion thickness in solutions of semi-flexible polymers near colloidal surfaces: analytical approximations.

We derive a simple, yet accurate approximate mean-field expression for the depletion thickness δsf of a solution of dilute semi-flexible polymers next to a hard surface. In the case of a hard wall this equation has the simple form δsf = δ0[1 − tanh(psf/δ0)], where psf accounts for the degree of flexibility and δ0 is the depletion thickness in the case of fully flexible polymers. For fixed polymer coil size, increasing the chain stiffness leads to a decrease in the depletion thickness. The approach is also extended to include higher polymer concentrations in the semidilute regime. The analytical expressions are in quantitative agreement with numerical self-consistent field computations. A remarkable finding is that there is a maximum in the depletion thickness as a function of the chain stiffness in the semidilute concentration regime. This also means that depletion attractions between colloidal particles reach a maximum for a certain chain stiffness, which may have important implications for the phase stability of colloid–polymer mixtures. The derived equations could be useful for the description of interactions in- and phase stability of mixtures of colloids and semi-flexible polymers.

Cylindrical confinement of solutions containing semiflexible macromolecules: surface-induced nematic order versus phase separation.

Solutions of semiflexible polymers confined in cylindrical pores with repulsive walls are studied by Molecular Dynamics simulations for a wide range of polymer concentrations. Both the case where both lengths are of the same order and the case when the persistence length by far exceeds the contour length are considered, and the enhancement of nematic order along the cylinder axis is characterized. With increasing density the character of the surface effect changes from depletion to the formation of a layered structure. For binary 50 : 50 mixtures of the two types of polymers an interplay between surface enrichment of the stiffer component and the isotropic-nematic transition is found, and a phase separated structure with cylindrical symmetry occurs, with the isotropic phase located around the cylinder axis. For melt densities the mixed nematic phase forms at the wall a layer with a screw-like structure of a tilted smectic phase. The observed behavior is tentatively interpreted in terms of the competition of the chain orientational entropy with entropy of mixing and excluded volume due to the wall.

Correction to Simulation Study of Entanglement in Semiflexible Polymer Melts and Solutions

Correction to Simulation Study of Entanglement in Semiflexible Polymer Melts and Solutions

Correction to Unified Entanglement Scaling for Flexible, Semiflexible, and Stiff Polymer Melts and Solutions

ADVERTISEMENT RETURN TO ISSUEPREVAddition/CorrectionNEXTORIGINAL ARTICLEThis notice is a correctionCorrection to Unified Entanglement Scaling for Flexible, Semiflexible, and Stiff Polymer Melts and SolutionsScott T. Milner*Scott T. MilnerMore by Scott T. MilnerCite this: Macromolecules 2020, 53, 20, 9152Publication Date (Web):October 15, 2020Publication History Published online15 October 2020Published inissue 27 October 2020https://doi.org/10.1021/acs.macromol.0c02231Copyright © 2020 American Chemical SocietyRIGHTS & PERMISSIONSArticle Views534Altmetric-Citations-LEARN ABOUT THESE METRICSArticle Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated. Share Add toView InAdd Full Text with ReferenceAdd Description ExportRISCitationCitation and abstractCitation and referencesMore Options Share onFacebookTwitterWechatLinked InReddit PDF (420 KB) Get e-Alerts Get e-Alerts

Simulation Study of Entanglement in Semiflexible Polymer Melts and Solutions

Multiple scaling arguments have been proposed to describe how the entanglement molecular weight depends on polymer architecture and concentration. The Lin–Noolandi (LN) scaling argument, well suppo...

Unified Entanglement Scaling for Flexible, Semiflexible, and Stiff Polymer Melts and Solutions

The entanglement length Ne is a key parameter for all entangled polymer fluids for which no comprehensive scaling theory yet exists. We have pieces of a theory; the Lin-Noolandi (LN) argument predicts Ne scaling for flexible chains that agrees with data on polymer melts. There are arguments for how Ne should depend on polymer concentration, but which are not obviously consistent with LN. Morse scaling describes entanglement for solutions of stiff chains, consistent with data. Everaers proposed an ansatz that Ne depends only on the arclength concentration, as if chains were uncrossable threads of vanishing thickness. This ansatz is consistent with simulations of bead-spring chains, but not with LN, as it has no role for packing length, the central parameter in LN scaling. We propose a comprehensive scaling theory that includes LN in one limit, thread ansatz in another, and reduces to Morse scaling for stiff chains. One new ingredient is that the typical distance of closest approach between two chains is governed by the packing length or chain diameter, whichever is larger. If a chain is sufficiently flexible and bulky, the packing length is relevant; but for stiffened bead-spring chains without side groups, the packing length is smaller than the chain diameter, so thread scaling applies. Our approach presents a consistent physical picture of entanglements in all regimes as close encounters between two chains. For solutions, we determine the entanglement probability between chain segments, and consistently describe the crossover between the Edwards and semidilute regimes.

Phase behavior of flexible and semiflexible polymers in solvents of varying quality.

The interplay of nematic order and phase separation in solutions of semiflexible polymers in solvents of variable quality is investigated by density functional theory (DFT) and molecular dynamics (MD) simulations. We studied coarse-grained models, with a bond-angle potential to control chain stiffness, for chain lengths comparable to the persistence length of the chains. We varied both the density of the monomeric units and the effective temperature that controls the quality of the implicit solvent. For very stiff chains, only a single transition from an isotropic fluid to a nematic is found, with a phase diagram of "swan-neck" topology. For less stiff chains, however, also unmixing between isotropic fluids of different concentration, ending in a critical point, occurs for temperatures above a triple point. The associated critical behavior is examined in the MD simulations and found compatible with Ising universality. Apart from this critical behavior, DFT calculations agree qualitatively with the MD simulations.

Viscoelastic Properties of Tightly Entangled Semiflexible Polymer Solutions

The dynamic viscoelasticity and birefringence of well-characterized semiflexible polymer solutions—cellulose tris(phenyl carbamate)/tricresyl phosphate solutions—were investigated over a wide concentration region ranging from a dilute to a tightly entangled regime where the entanglement spacing is smaller than the Kuhn segment size. The stress optical rule did not hold true in the tightly entangled regime, indicating that the molecular origin of stress is not simply attributed to the orientation of segments but to the bending of chains, which does not contribute to birefringence significantly. At high frequencies, the power law behavior with the exponent 3/4 due to the tension stress was observed for the complex shear moduli at all the concentrations.

Nematic order in solutions of semiflexible polymers: Hairpins, elastic constants, and the nematic-smectic transition.

Coarse-grained models of lyotropic solutions of semiflexible polymers are studied by both molecular dynamics simulations and density functional theory calculations, using an implicit solvent bead-spring model with a bond-angle potential. We systematically vary the monomer density, persistence length, and contour length over a wide range and explore the full range from the isotropic-nematic transition to the nematic-smectic transition. In the nematic regime, we span the entire regime from rigid-rod like polymers to thin wormlike chains, confined in effective straight tubes caused by the collective nematic effective ordering field. We show that the distribution of bond angles relative to the director is well described by a Gaussian, irrespective of whether the chains are rod-like or rather flexible. However, the related concept of "deflection length" is shown to make sense only in the latter case for rather dilute solutions since otherwise the deflection length is of the order of about two bond lengths only. When the solution is semi-dilute, a substantial renormalization of the persistence length occurs, while this effect is absent in the isotropic phase even at rather high monomer densities. The effective radii of the "tubes" confining the chains in the related description of orientational ordering are significantly larger than the distances between neighboring chains, providing evidence for a pronounced collective character of orientational fluctuations. Hairpins can be identified close to the isotropic-nematic transition, and their probability of occurrence agrees qualitatively with the Vroege-Odijk theory. The corresponding theoretical predictions for the elastic constants, however, are not in good agreement with the simulations. We attribute the shortcomings of the theories to their neglect of the coupling between local density and orientational fluctuations. Finally, we detected for this model a transition to a smectic phase for reduced monomer densities near 0.7.

Incorporating particle flexibility in a density functional description of nematics and cholesterics

ABSTRACTWe describe a general implementation of the Fynewever–Yethiraj density functional theory (DFT) for the investigation of nematic and cholesteric self-assembly in arbitrary solutions of semi-flexible polymers. The basic assumptions of the theory are discussed in the context of other generalised Onsager descriptions for flexible polyatomic systems. The location of the isotropic-to-nematic phase transition is found to be in good agreement with molecular simulations for elongated chains up to relatively high polymer flexibilities, although the predictions of the theory in the nematic regime lead to gradual underestimations of order parameters with decreasing particle stiffness. This shortcoming is attributed to increasing overestimations of the molecular conformational entropy in higher-density phases, which may not be easily addressed in the formalism of DFT for realistic particle models. Practical consequences of these limitations are illustrated through the application of DFT to systems of near-persistence-length DNA duplexes, whose cholesteric behaviour is found to be strongly contingent on their detailed accessible conformational space in concentrated solutions.

Adsorption and structure formation of semiflexible polymers on spherical surfaces

Rigid spheres with a short-range attractive potential are taken as a coarse-grained model of vesicles, which contain a solution of semiflexible polymers in their interior. Assuming good solvent conditions with an implicit description of the solvent, effective monomers experience bond-length and bond-angle potentials as well as excluded-volume interaction. Due to the attractive vesicle surface, phase separation occurs between a thin shell of adsorbed monomers at the surface and a rather dilute, and therefore, disordered polymer solution in the sphere interior. While at a planar attractive surface the wormlike chains would exhibit liquid crystalline (nematic and smectic) order, the curvature of the sphere frustrates this order, causing the existence of defects in the orientational order. In the dense monolayer of surface-attached polymers, the stiff chains acquire a curvature corresponding to banana-like shapes for rather small contour length, causing a local nematic order of biaxial type. These orderings are compared to related phenomena in spheres densely filled with semiflexible polymers. We complement the results obtained by Molecuar Dynamics simulations by corresponding Density Functional Theory calculations.

Fiber Network Formation in Semi-Flexible Polymer Solutions: An Exploratory Computational Study.

The formation of gels through the bundling of semi-flexible polymer chains into fiber networks is ubiquitous in diverse manufactured and natural materials, and, accordingly, we perform exploratory molecular dynamics simulations of a coarse-grained model of semi-flexible polymers in a solution with attractive lateral interchain interactions to understand essential features of this type of gel formation. After showing that our model gives rise to fibrous gels resembling real gels of this kind, we investigate how the extent of fiber bundling influences the “melting” temperature, , and the emergent rigidification of model bundled fibers having a fixed number of chains, N, within them. Based on our preliminary observations, we suggest the fiber size is kinetically selected by a reduced thermodynamic driving force and a slowing of the dynamics within the fibers associated with their progressive rigidification with the inclusion of an increasing number of chains in the bundle.

Disentangling entanglements in biopolymer solutions

Reptation theory has been highly successful in explaining the unusual material properties of entangled polymer solutions. It reduces the complex many-body dynamics to a single-polymer description, where each polymer is envisaged to be confined to a tube through which it moves in a snake-like fashion. For flexible polymers, reptation theory has been amply confirmed by both experiments and simulations. In contrast, for semiflexible polymers, experimental and numerical tests are either limited to the onset of reptation, or were performed for tracer polymers in a fixed, static matrix. Here, we report Brownian dynamics simulations of entangled solutions of semiflexible polymers, which show that curvilinear motion along a tube (reptation) is no longer the dominant mode of dynamics. Instead, we find that polymers disentangle due to correlated constraint release, which leads to equilibration of internal bending modes before polymers diffuse the full tube length. The physical mechanism underlying terminal stress relaxation is rotational diffusion mediated by disentanglement rather than curvilinear motion along a tube.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics.

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

Equilibrium Dynamics and Shear Rheology of Semiflexible Polymers in Solution

We study the structure and dynamics of semidilute solutions of semiflexible polymers at rest and under shear using hybrid molecular dynamics simulations that take hydrodynamic interactions into account. We show that the polymer center-of-mass diffusion coefficient significantly decreases with increasing chain stiffness at fixed monomer density. The zero-shear viscosity shows a corresponding increase due to the intermolecular interactions of stiffer chains. We apply steady shear flow to the polymer solutions and show that at high shear rates the flow properties become almost independent of polymer stiffness. We characterize the polymer conformations under shear and find that in this regime polymers are elongated and aligned along the flow direction, but semiflexible polymers surprisingly form an increased number of hairpin-like conformations due to increased tumbling.

Dynamics of Semiflexible Polymer Solutions in theTightly EntangledConcentration Regime

I corroborate an important experimental evidence reported by Schuldt et al. [Phys. Rev. Lett. 2016, 117, 197801] revealing the incapability of the current theoretical framework to fully describe the dynamics of semiflexible polymer solutions in the tightly entangled concentration regime. These results have been endorsed here by means of previously published, but overlooked, data by Tassieri et al. [Phys. Rev. Lett. 2008, 101, 198301; Biophys. J. 2008, 94, 2170]. The ensemble of information provides a strong evidence that the scaling law of the plateau modulus as a function of polymer’s concentration and persistence length, i.e., G0 ∝ cαLpβ, should have both the exponents positive, in contrast with all the existing predictions converging on similar values of α > 0 and β < 0—hence the need of new theoretical models able to better interpret the contribution of the polymer bending rigidity to the viscoelastic properties of the polymer network.