Mechanical desorption of a single chain: unusual aspects of phase coexistence at a first-order transition

Elasticity of a Grafted Rod-like Filament with Fluctuating Bending Stiffness

Quite often polymers exhibit different elastic behavior depending on the statistical ensemble (Gibbs vs. Helmholtz). This is an effect of strong fluctuations. In particular, two-state polymers, which locally or globally fluctuate between two classes of microstates, can exhibit strong ensemble inequivalence with negative elastic moduli (extensibility or compressibility) in the Helmholtz ensemble. Two-state polymers consisting of flexible beads and springs have been studied extensively. Recently, similar behavior was predicted in a strongly stretched wormlike chain consisting of a sequence of reversible blocks, fluctuating between two values of the bending stiffness (the so called reversible wormlike chain, rWLC). In this article, we theoretically analyse the elasticity of a grafted rod-like semiflexible filament which fluctuates between two states of bending stiffness. We consider the response to a point force at the fluctuating tip in both the Gibbs and the Helmholtz ensemble. We also calculate the entropic force exerted by the filament on a confining wall. This is done in the Helmholtz ensemble and, under certain conditions, it yields negative compressibility. We consider a two-state homopolymer and a two-block copolymer with two-state blocks. Possible physical realizations of such a system would be grafted DNA or carbon nanorods undergoing hybridization, or grafted F-actin bundles undergoing collective reversible unbinding.

Tensile elasticity of a freely jointed chain with reversible hinges

Many biopolymers exhibit reversible conformational transitions within the chain, which affect their bending stiffness and their response to a stretching force. For example, double stranded DNA may have denatured "bubbles" of unzipped single strands which open and close randomly. In other polymers, the transitions may be due to the reversible attachment and detachment of ligands on ligand-receptor complexes along the backbone. Semiflexible bundles under tension formed by the reversible attachment of cross-linkers, on a coarse-grained level, exhibit similar behaviour. The simplest theoretical model which captures what the above mentioned systems have in common is a freely jointed chain (FJC) with reversible hinges. Each hinge can be open, as in the usual FJC, or closed forcing the adjacent segments to align (stretch). In this article, we analyse it in the Gibbs ensemble. Remarkably, even though the usual FJC in the thermodynamic limit exhibits ensemble equivalence, the reversible FJC exhibits ensemble inequivalence. Even though a mean field treatment suggests a continuous phase transition to a fully hinged state at a certain force, the generating function method ("necklace model") shows that there is no phase transition. However, there is a crossover between the two states with clearly different responses. In the low force (linear response) regime, the reversible FJC has higher tensile compliance than its usual counterpart. In contrast, in the strong force regime, the tensile compliance of the reversible FJC is much lower than that of the usual FJC.

Statistical ensemble inequivalence for flexible polymers under confinement in various geometries

The problem of statistical ensemble inequivalence for single polymers has been the subject of intense research. In a recent publication, we show that even though the force-extension relation of a free Gaussian chain exhibits ensemble equivalence, confinement to half-space due to tethering to a planar substrate induces significant inequivalence [S. Dutta and P. Benetatos, Soft Matter, 2018, 14, 6857-6866]. In the present article, we extend that work to the conformational response to confining forces distributed over surfaces. We analyze in both the Helmholtz and the Gibbs ensemble the pressure-volume equation of state of a chain in rectangular, spherical, and cylindrical confinement. We especially consider the case of a directed polymer in a cylinder. We also analyze the case of a tethered chain inside a rectangular box, a sphere, and outside a sphere. In general, confinement causes significant ensemble inequivalence. Remarkably, we recover ensemble equivalence in the limit of squashing confinement. We trace the ensemble inequivalence to the persistence of strong fluctuations. Our work may be useful in the interpretation of single molecule experiments and caging phenomena.

Inequivalence of fixed-force and fixed-extension statistical ensembles for a flexible polymer tethered to a planar substrate

Recent advances in single macromolecule experiments have sparked interest in the ensemble dependence of force-extension relations. The thermodynamic limit may not be attainable for such systems, which leads to inequivalence of the fixed-force and the fixed-extension ensembles. We consider an ideal Gaussian chain described by the Edwards Hamiltonian with one end tethered to a rigid planar substrate. We analytically calculate the force-extension relation in the two ensembles and we show their inequivalence, which is caused by the confinement of the polymer to half space. The inequivalence is quite remarkable for strong compressional forces. We also perform Monte-Carlo simulations of a tethered wormlike chain with contour length 20 times its persistence length, which corresponds to experiments measuring the conformations of DNA tethered to a wall. The simulations confirm the ensemble inequivalence and qualitatively agree with the analytical predictions of the Gaussian model. Our analysis shows that confinement due to tethering causes ensemble inequivalence, irrespective of the polymer model.

Phase transitions in single macromolecules: Loop-stretch transition versus loop adsorption transition in end-grafted polymer chains

We use Brownian dynamics simulations and analytical theory to compare two prominent types of single molecule transitions. One is the adsorption transition of a loop (a chain with two ends bound to an attractive substrate) driven by an attraction parameter ε and the other is the loop-stretch transition in a chain with one end attached to a repulsive substrate, driven by an external end-force F applied to the free end. Specifically, we compare the behavior of the respective order parameters of the transitions, i.e., the mean number of surface contacts in the case of the adsorption transition and the mean position of the chain end in the case of the loop-stretch transition. Close to the transition points, both the static behavior and the dynamic behavior of chains with different length N are very well described by a scaling ansatz with the scaling parameters (ε - ε*)N^ϕ (adsorption transition) and (F - F^*)N^ν (loop-stretch transition), respectively, where ϕ is the crossover exponent of the adsorption transition and ν is the Flory exponent. We show that both the loop-stretch and the loop adsorption transitions provide an exceptional opportunity to construct explicit analytical expressions for the crossover functions which perfectly describe all simulation results on static properties in the finite-size scaling regime. Explicit crossover functions are based on the ansatz for the analytical form of the order parameter distributions at the respective transition points. In contrast to the close similarity in equilibrium static behavior, the dynamic relaxation at the two transitions shows qualitative differences, especially in the strongly ordered regimes. This is attributed to the fact that the surface contact dynamics in a strongly adsorbed chain is governed by local processes, whereas the end height relaxation of a strongly stretched chain involves the full spectrum of Rouse modes.

Attractive and repulsive polymer-mediated forces between scale-free surfaces

We consider forces acting on objects immersed in, or attached to, long fluctuating polymers. The confinement of the polymer by the obstacles results in polymer-mediated forces that can be repulsive (due to loss of entropy) or attractive (if some or all surfaces are covered by adsorbing layers). The strength and sign of the force in general depends on the detailed shape and adsorption properties of the obstacles but assumes simple universal forms if characteristic length scales associated with the objects are large. This occurs for scale-free shapes (such as a flat plate, straight wire, or cone) when the polymer is repelled by the obstacles or is marginally attracted to it (close to the depinning transition where the absorption length is infinite). In such cases, the separation h between obstacles is the only relevant macroscopic length scale, and the polymer-mediated force equals Ak_{B}T/h, where T is temperature. The amplitude A is akin to a critical exponent, depending only on geometry and universality of the polymer system. The value of A, which we compute for simple geometries and ideal polymers, can be positive or negative. Remarkably, we find A=0 for ideal polymers at the adsorption transition point, irrespective of shapes of the obstacles, i.e., at this special point there is no polymer-mediated force between obstacles (scale free or not).

Anomalous critical slowdown at a first order phase transition in single polymer chains

Using Brownian dynamics, we study the dynamical behavior of a polymer grafted onto an adhesive surface close to the mechanically induced adsorption-stretching transition. Even though the transition is first order (in the infinite chain length limit, the stretching degree of the chain jumps discontinuously), the characteristic relaxation time is found to grow according to a power law as the transition point is approached. We present a dynamic effective interface model which reproduces these observations and provides an excellent quantitative description of the simulation data. The generic nature of the theoretical model suggests that the unconventional mixing of features that are characteristic for first-order transitions (a jump in an order parameter) and features that are characteristic of critical points (an anomalous slowdown) may be a common phenomenon in force-driven phase transitions of macromolecules.

Polymers undergoing inhomogeneous adsorption: Order parameters for a partially directed walk model

We consider partially directed walk models of polymers undergoing inhomogeneous adsorption. The inhomogeneity can be in the polymer, in the surface, or in both. For the cases where the polymer is either a homopolymer or a strictly alternating copolymer and where the surface is either homogeneous or has stripes of width 1, we calculate detailed order parameters and show that these provide important information about the ways in which the polymer adsorbs.

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.

Pulling alternating copolymers adsorbed on a striped surface

We consider a partially directed walk model of a strictly alternating copolymer adsorbing on a striped surface where the energy is associated with the numbers of the two types of monomers adsorbed on the two types of surface sites. A force is applied to the last monomer and the polymer responds to this force, sometimes by desorbing. The force can be applied at various angles, with the surface component parallel or perpendicular (or at some other angle) to the stripe direction. The desorption behavior is strongly dependent on the force direction and the response gives information about the shape and direction of the polymer adsorbed on the surface, especially at low temperatures. In some cases the ground state is degenerate and this also has an important effect on the temperature dependence of the critical force needed for desorption. We give a complete solution of the problem using generating function techniques and an approximate treatment that is especially useful at low temperatures and helps in our physical understanding of the situation.