Stiffening Transition in Semiflexible Cyclic Macromolecules

Conformational Transition of Semiflexible Ring Polyelectrolyte in Tetravalent Salt Solutions: A Simple Numerical Modeling without the Effect of Twisting

In this work, the conformational behaviors of ring polyelectrolyte in tetravalent salt solutions are discussed in detail through molecular dynamics simulation. For simplification, here we have neglected the effect of the twisting interaction, although it has been well known that both bending and twisting interactions play a deterministic in the steric conformation of a semiflexible ring polymer. The salt concentration CS and the bending energy b take a decisive role in the conformation of the ring polyelectrolyte (PE). Throughout our calculations, the b varies from b = 0 (freely joint chain) to b = 120. The salt concentration CS changes in the range of 3.56 × 10-4 M ≤ CS ≤ 2.49 × 10-1 M. Upon the addition of salt, ring PE contracts at first, subsequently re-expands. More abundant conformations are observed for a semiflexible ring PE. For b = 10, the conformation of semiflexible ring PE shifts from the loop to two-racquet-head spindle, then it condenses into toroid, finally arranges into coil with the increase of CS. As b increases further, four phase transitions are observed. The latter two phase transitions are different. The semiflexible ring PE experiences transformation from toroid to two racquet head spindle, finally to loop in the latter two phase transitions. Its conformation is determined by the competition among the bending energy, cation-bridge, and entropy. Combined, our findings indicate that the conformations of semiflexible ring PE can be controlled by changing the salt concentration and chain stiffness.

Structural Behavior of a Semiflexible Polymer Chain in an Array of Nanoposts

The structural properties of a flexible and semiflexible circular chain confined in an array of parallel nanoposts with a square lattice cross-sectional projection were studied using coarse-grained molecular dynamics simulations. To address the effect of the circular topology, a comparison with linear analogs was also carried out. In the interpretation of the chain structural properties, the geometry of the post array is considered as a combination of a channel approximating the interstitial volume with the diameter dc and a slit approximating the passage aperture with the width wp. The number of interstitial volumes occupied by a chain monotonically increases with the decreasing ratio dc/wp regardless of the way the geometry of the post array is varied. However, depending on how the array geometry is modified, the chain span along the posts displays a monotonic (constant post separation) or a non-monotonic behavior (constant passage width) when plotted as a function of the post diameter. In the case of monotonic trend, the width of interstitial spaces increases with the increasing chain occupation number, while, in the case of non-monotonic trend, the width of interstitial spaces decreases with the increasing chain occupation number. In comparison with linear topology, for circular topology, the stiffness affects more significantly the relative chain extension along the posts and less significantly the occupation number. The geometrical parameters of the post arrays are stored in the single-chain structure factors. The characteristic humps are recognized in the structure factor which ensue from the local increase in the density of segments in the circular chains presented in an interstitial volume or from the correlation of parallel chain fragments separated by a row of posts. Although the orientation correlations provide qualitative information about the chain topology and the character of confinement within a single interstitial volume, information about the array periodicity is missing.

Comparison of a stripe and slab confinement for ring and linear macromolecules in nanochannel

The combined effects of the channel asymmetry and the closed chain topology on the chain extension, structure factor, and the orientation correlations were studied using coarse-grained molecular dynamics simulations for moderate chain lengths. These effects are related to applications in linearization experiments with a DNA molecule in nanofluidic devices. According to the aspect ratio, the channels are classified as a stripe or slabs. The chain segments do not have any freedom to move in the direction of the narrowest stripe size, being approximately the same size as the segment size. The chains of both ring and linear topologies are extended more in a stripe than in a slab; this effect is strengthened for a ring. For a ring in a stripe, the extension-confinement strength dependence leads to effective Flory exponents even larger than 3/4, which is characteristic for a self-avoiding two-dimensional chain. While the chain extension-confinement strength dependence for both topologies conforms to the de Gennes regime in a stripe, a linear chain undergoes gradual transition to the pseudoideal regime as the slab height increases in the slab-like confinement. For a confined circle, the onset of the pseudoideal regime is shifted to larger slab heights. The structure factor confirms the absence of the pseudoideal and extended de Gennes regime in a stripe and the transition from the extended to the pseudoideal regime of a circular and linear chain upon increasing the slab heights.

Maximum entropy principle applied to semiflexible ring polymers

Based on the success of the maximum entropy principle (MEP) in the study of semiflexible treelike polymers [M. Dolgushev and A. Blumen, J. Chem. Phys. 131, 044905 (2009)], it is of much interest to establish MEP's potential for general semiflexible polymers which contain loops. Here, we embark on this endeavor by considering discrete semiflexible polymer rings in a Rouse-type scheme. Now, for treelike polymers a beads-and-bonds (i.e., a discrete) picture is essential for an easy inclusion of branching points. Moreover, one may envisage (similar to our former work [M. Dolgushev and A. Blumen, J. Chem. Phys. 131, 044905 (2009)]) to impose for each angle between two bonds a distinct stiffness condition. Working in this way leads already for a polymer ring to a complicated problem. Hence, we follow a reduced variational approach as applied earlier to polymer chains, in which a single Lagrange multiplier is used for each set of identical conditions imposed on topologically equivalent bonds and bonds' orientations. In this way, we obtain for the discrete ring an analytically closed form which involves Chebyshev polynomials. This expression turns out to lead to a series of solutions: Apart from the regular solution, several other solutions appear. One may be tempted to discard the other solutions, since for them the potential energy matrix is not positive definite. A more careful analysis based on topological features suggests, however, that such solutions can be assigned to rings displaying knots. Monte Carlo simulations which take excluded volume interactions into account agree with our interpretation.