Stretching Semiflexible Polymer Chains: Evidence for the Importance of Excluded Volume Effects from Monte Carlo Simulation

Abstract

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice ($d=3$ dimensions) and square lattice ($d=2$ dimensions), varying chain stiffness by an energy penalty ${\epsilon }_b$ for chain bending. In the absence of excluded volume interactions, the persistence length ${\ell }_p$ of the polymers would then simply be ${\ell }_p={\ell }_b(2d-2{)}^{-1}{q}_b^{-1}$ with $q_b=\exp (-{\epsilon }_b/k_{B}T)$, the bond length ${\ell }_b$ being the lattice spacing, and $k_{B}T$ is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both $q_b$ and the chain length $N$ are varied over a wide range $(0.005\le q_b\le 1,N\le 50000$), and also a stretching force $f$ is applied to one chain end (fixing the other end at the origin). In the absence of this force, in $d=2$ a single crossover from rod-like behavior (for contour lengths less than ${\ell }_p$) to swollen coils occurs, invalidating the Kratky-Porod model, while in $d=3$ a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in $d=2$, while theories based on the Kratky-Porod model are found to work in $d=3$ for stiff chains in an intermediate regime of chain extensions. While for $q_b\ll 1$ in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g. bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.

23 pages, 17 figures, 2 tables, to be published in J. Chem. Phys. (2011)

References in corpus (7)

• Polymer chain stiffness versus excluded volume: A Monte Carlo study of the crossover towards the wormlike chain model
• A review of Monte Carlo simulations of polymers with PERM
• Breakdown of the Kratky-Porod Wormlike Chain Model for Semiflexible Polymers in Two Dimensions
• Stretching Homopolymers
• Forced-induced desorption of a polymer chain adsorbed on an attractive surface - Theory and Computer Experiment
• Force induced stretched state: Effects of temperature
• Pulling self-interacting polymers in two-dimensions

Cited by in corpus (13)

• Conformational properties of active semiflexible polymers
• Semiflexible Polymers Under Good Solvent Conditions Interacting With Repulsive Walls
• Scattering function of semiflexible polymer chains under good solvent conditions
• Conformations, Transverse Fluctuations and Crossover Dynamics of a Semi-Flexible Chain in Two Dimensions
• Monte Carlo Simulations of Lattice Models for Single Polymer Systems
• Rigidity-induced scale invariance in polymer ejection from capsid
• Dynamic scaling theory of the forced translocation of a semi-flexible polymer through a nanopore
• Elasticity of a DNA chain dotted with bubbles under force
• Transitions of tethered chain molecules under tension
• Off-Lattice Markov Chain Monte Carlo Simulations of Mechanically Driven Polymers
• Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy
• Lattice models of directed and semiflexible polymers in anisotropic environment
• The extraordinary importance of self-avoiding behavior in two-dimensional polymers: Insights from large-deviation theory