Pulling self-interacting polymers in two dimensions

Solvent-dependent critical properties of polymer adsorption

Advanced chain-growth computer simulation methodologies have been employed for a systematic statistical analysis of the critical behavior of a polymer adsorbing at a substrate. We use finite-size scaling techniques to investigate the solvent-quality dependence of critical exponents, critical temperature, and the structure of the phase diagram. Our study covers all solvent effects from the limit of super-self-avoiding walks, characterized by effective monomer-monomer repulsion, to poor solvent conditions that enable the formation of compact polymer structures. The results significantly benefit from taking into account corrections to scaling.

Stretching tethered polymer chains: density functional approach

We propose application of density functional theory to calculate the force acting on a selected segment of a tethered polymer chain that leads to stretching the chain. The density functional allows one to determine the effects due to the presence of other chains and solvent molecules. For high and moderate solvent densities the plot of the force versus the distance of the segment from the surface exhibits oscillatory behavior that has not been predicted by other approaches.

Adsorption of a single polymer chain on a surface: effects of the potential range

We investigate the effects of the range of adsorption potential on the equilibrium behavior of a single polymer chain end-attached to a solid surface. The exact analytical theory for ideal lattice chains interacting with a planar surface via a box potential of depth U and width W is presented and compared to continuum model results and to Monte Carlo (MC) simulations using the pruned-enriched Rosenbluth method for self-avoiding chains on a simple cubic lattice. We show that the critical value U(c) corresponding to the adsorption transition scales as W(-1/ν), where the exponent ν=1/2 for ideal chains and ν≈3/5 for self-avoiding walks. Lattice corrections for finite W are incorporated in the analytical prediction of the ideal chain theory U(c)≈(π(2)/24)(W+1/2)(-2) and in the best-fit equation for the MC simulation data U(c)=0.585(W+1/2)(-5/3). Tail, loop, and train distributions at the critical point are evaluated by MC simulations for 1≤W≤10 and compared to analytical results for ideal chains and with scaling theory predictions. The behavior of a self-avoiding chain is remarkably close to that of an ideal chain in several aspects. We demonstrate that the bound fraction θ and the related properties of finite ideal and self-avoiding chains can be presented in a universal reduced form: θ(N,U,W)=θ(NU(c),U/U(c)). By utilizing precise estimations of the critical points we investigate the chain length dependence of the ratio of the normal and lateral components of the gyration radius. Contrary to common expectations this ratio attains a limiting universal value <R(g[perpendicular])(2)>/<R(g[parallel])(2)>=0.320±0.003 only at N~5000. Finite-N corrections for this ratio turn out to be of the opposite sign for W=1 and for W≥2. We also study the N dependence of the apparent crossover exponent φ(eff)(N). Strong corrections to scaling of order N(-0.5) are observed, and the extrapolated value φ=0.483±0.003 is found for all values of W. The strong correction to scaling effects found here explain why for smaller values of N, as used in most previous work, misleadingly large values of φ(eff)(N) were identified as the asymptotic value for the crossover exponent.

Stretching semiflexible polymer chains: evidence for the importance of excluded volume effects from Monte Carlo simulation

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 ε(b) for chain bending. In the absence of excluded volume interactions, the persistence length l(p) of the polymers would then simply be l(p) = l(b)(2d - 2)(-1)q(b) (-1) with q(b) = exp(-ε(b)/k(B)T), the bond length l(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 ≤ q(b) ≤ 1, N ≤ 50,000), 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 l(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) ≪ 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.

Theory of Biopolymer Stretching at High Forces

We provide a unified theory for the high force entropic elasticity of biopolymers solely in terms of the persistence length, ξ_p , and the monomer spacing, a. When the force f>ℱ _h ~ k_BTξ_p /a² the biopolymers behave as freely jointed chains (FJCs) while in the range ℱ _l ~ k_BT/ξ_p <f<ℱ _h the worm-like chain (WLC) is a better model. We show that ξ_p can be estimated from the force extension curve (FEC) at the extension x ≈ ¹/₂ (normalized by the contour length of the biopolymer). After validating the theory using simulations, we provide a quantitative analysis of the FECs for a diverse set of biopolymers (dsDNA, ssRNA, ssDNA, polysaccharides, and unstructured PEVK domain of titin) for x ≥ ¹/₂. The success of a specific polymer model (FJC or WLC) to describe the FEC of a given biopolymer is naturally explained by the theory. Only by probing the response of biopolymers over a wide range of forces can the f-dependent elasticity be fully described.

Adsorption of externally stretched two-dimensional flexible and semiflexible polymers near an attractive wall

We study analytically a model of a two-dimensional partially directed flexible or semiflexible polymer, attached to an attractive wall which is perpendicular to the preferred direction. In addition, the polymer is stretched by an externally applied force. We find that the wall has a dramatic effect on the polymer. For wall attraction epsilon1 smaller than the nonsequential nearest-neighbor attraction epsilon, the fraction of monomers at the wall is zero and the model is the same as that of a polymer without a wall. However, for epsilon1 greater than epsilon, the fraction of monomers at the wall undergoes a first-order transition from unity at low temperature and small force, to zero at higher temperatures and forces. We present phase diagram for this transition. Our results are confirmed by Monte Carlo simulations.