Semiflexible macromolecules in quasi-one-dimensional confinement: Discrete versus continuous bond angles.

The conformations of semiflexible polymers in two dimensions confined in a strip of width D are studied by computer simulations, investigating two different models for the mechanism by which chain stiffness is realized. One model (studied by molecular dynamics) is a bead-spring model in the continuum, where stiffness is controlled by a bond angle potential allowing for arbitrary bond angles. The other model (studied by Monte Carlo) is a self-avoiding walk chain on the square lattice, where only discrete bond angles (0° and ±90°) are possible, and the bond angle potential then controls the density of kinks along the chain contour. The first model is a crude description of DNA-like biopolymers, while the second model (roughly) describes synthetic polymers like alkane chains. It is first demonstrated that in the bulk the crossover from rods to self-avoiding walks for both models is very similar, when one studies average chain linear dimensions, transverse fluctuations, etc., despite their differences in local conformations. However, in quasi-one-dimensional confinement two significant differences between both models occur: (i) The persistence length (extracted from the average cosine of the bond angle) gets renormalized for the lattice model when D gets less than the bulk persistence length, while in the continuum model it stays unchanged. (ii) The monomer density near the repulsive walls for semiflexible polymers is compatible with a power law predicted for the Kratky-Porod model in the case of the bead-spring model, while for the lattice case it tends to a nonzero constant across the strip. However, for the density of chain ends, such a constant behavior seems to occur for both models, unlike the power law observed for flexible polymers. In the regime where the bulk persistence length ℓp is comparable to D, hairpin conformations are detected, and the chain linear dimensions are discussed in terms of a crossover from the Daoud/De Gennes "string of blobs"-picture to the flexible rod picture when D decreases and/or the chain stiffness increases. Introducing a suitable further coarse-graining of the chain contours of the continuum model, direct estimates for the deflection length and its distribution could be obtained.

Coil-bridge transition in a single polymer chain as an unconventional phase transition: theory and simulation.

The coil-bridge transition in a self-avoiding lattice chain with one end fixed at height H above the attractive planar surface is investigated by theory and Monte Carlo simulation. We focus on the details of the first-order phase transition between the coil state at large height H ⩾ Htr and a bridge state at H ⩽ Htr, where Htr corresponds to the coil-bridge transition point. The equilibrium properties of the chain were calculated using the Monte Carlo pruned-enriched Rosenbluth method in the moderate adsorption regime at (H/Na)tr ⩽ 0.27 where N is the number of monomer units of linear size a. An analytical theory of the coil-bridge transition for lattice chains with excluded volume interactions is presented in this regime. The theory provides an excellent quantitative description of numerical results at all heights, 10 ⩽ H/a ⩽ 320 and all chain lengths 40 < N < 2560 without free fitting parameters. A simple theory taking into account the effect of finite extensibility of the lattice chain in the strong adsorption regime at (H/Na)tr ⩾ 0.5 is presented. We discuss some unconventional properties of the coil-bridge transition: the absence of phase coexistence, two micro-phases involved in the bridge state, and abnormal behavior in the microcanonical ensemble.

Free Standing Dry and Stable Nanoporous Polymer Films Made through Mechanical Deformation.

A new straight forward approach to create nanoporous polymer membranes with well defined average pore diameters is presented. The method is based on fast mechanical deformation of highly entangled polymer films at high temperatures and a subsequent quench far below the glass transition temperature Tg . The process is first designed generally by simulation and then verified for the example of polystyrene films. The methodology does not need any chemical processing, supporting substrate, or self assembly process and is solely based on polymer inherent entanglement effects. Pore diameters are of the order of ten polymer reptation tube diameters. The resulting membranes are stable over months at ambient conditions and display remarkable elastic properties.

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.

How to define variation of physical properties normal to an undulating one-dimensional object.

One-dimensional flexible objects are abundant in physics, from polymers to vortex lines to defect lines and many more. These objects structure their environment and it is natural to assume that the influence these objects exert on their environment depends on the distance from the line object. But how should this be defined? We argue here that there is an intrinsic length scale along the undulating line that is a measure of its stiffness (i.e., orientational persistence), which yields a natural way of defining the variation of physical properties normal to the undulating line. We exemplify how this normal variation can be determined from a computer simulation for the case of a so-called bottle-brush polymer, where side chains are grafted onto a flexible backbone.

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 /=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.