The impact of thermal fluctuations on the translocation dynamics of a polymer chain driven through a narrow pore has been investigated theoretically and by means of extensive Molecular-Dynamics (MD) simulation. The theoretical consideration is based
on the so-called velocity Langevin (V-Langevin) equation which determines the progress of the translocation in terms of the number of polymer segments, $s(t)$, that have passed through the pore at time $t$ due to a driving force $f$. The formalism is based only on the assumption that, due to thermal fluctuations, the translocation velocity $v=dot{s}(t)$ is a Gaussian random process as suggested by our MD data. With this in mind we have derived the corresponding Fokker-Planck equation (FPE) which has a nonlinear drift term and diffusion term with a {em time-dependent} diffusion coefficient $D(t)$. Our MD simulation reveals that the driven translocation process follows a {em super}diffusive law with a running diffusion coefficient $D(t) propto t^{gamma}$ where $gamma < 1$. This finding is then used in the numerical solution of the FPE which yields an important result: for comparatively small driving forces fluctuations facilitate the translocation dynamics. As a consequence, the exponent $alpha$ which describes the scaling of the mean translocation time $<tau>$ with the length $N$ of the polymer, $<tau> propto N^{alpha}$ is found to diminish. Thus, taking thermal fluctuations into account, one can explain the systematic discrepancy between theoretically predicted duration of a driven translocation process, considered usually as a deterministic event, and measurements in computer simulations. In the non-driven case, $f=0$, the translocation is slightly subdiffusive and can be treated within the framework of fractional Brownian motion (fBm).
The force-assisted desorption kinetics of a macromolecule from adhesive surface is studied theoretically, using the notion of tensile (Pincus) blobs, as well as by means of Monte-Carlo (MC) and Molecular Dynamics (MD) simulations. We show that the ch
ange of detached monomers with time is governed by a differential equation which is equivalent to the nonlinear porous medium equation (PME), employed widely in transport modeling of hydrogeological systems. Depending on the pulling force and the strength of adsorption, three kinetic regimes can be distinguished: (i) trumpet (weak adsorption and small pulling force), (ii) stem-trumpet (weak adsorption and moderate force), and (iii) stem (strong adsorption and large force). Interestingly, in all regimes the number of desorbed beads $M(t)$, and the height of the first monomer (which experiences a pulling force) $R(t)$ above the surface follow an universal square-root-of-time law. Consequently, the total time of detachment $<tau_d>$, scales with polymer length $N$ as $<tau_d> propto N^2$. Our main theoretical conclusions are tested and found in agreement with data from extensive MC- and MD-simulations.
We suggest a theoretical description of the force-induced translocation dynamics of a polymer chain through a nanopore. Our consideration is based on the tensile (Pincus) blob picture of a pulled chain and the notion of propagating front of tensile f
orce along the chain backbone, suggested recently by T. Sakaue. The driving force is associated with a chemical potential gradient that acts on each chain segment inside the pore. Depending on its strength, different regimes of polymer motion (named after the typical chain conformation, trumpet, stem-trumpet, etc.) occur. Assuming that the local driving and drag forces are equal (i.e., in a quasi-static approximation), we derive an equation of motion for the tensile front position $X(t)$. We show that the scaling law for the average translocation time $<tau>$ changes from $<tau> sim N^{2 u}/f^{1/ u}$ to $<tau> sim N^{1+ u}/f$ (for the free-draining case) as the dimensionless force ${widetilde f}_{R} = a N^{ u}f /T$ (where $a$, $N$, $ u$, $f$, $T$ are the Kuhn segment length, the chain length, the Flory exponent, the driving force, and the temperature, respectively) increases. These and other predictions are tested by Molecular Dynamics (MD) simulation. Data from our computer experiment indicates indeed that the translocation scaling exponent $alpha$ grows with the pulling force ${widetilde f}_{R}$) albeit the observed exponent $alpha$ stays systematically smaller than the theoretically predicted value. This might be associated with fluctuations which are neglected in the quasi-static approximation.
We analyze the dynamics of desorption of a polymer molecule which is pulled at one of its ends with force $f$, trying to desorb it. We assume a monomer to desorb when the pulling force on it exceeds a critical value $f_{c}$. We formulate an equation
for the average position of the $n^{th}$ monomer, which takes into account excluded volume interaction through the blob-picture of a polymer under external constraints. The approach leads to a diffusion equation with a $p$-Laplacian for the propagation of the stretching along the chain. This has to be solved subject to a moving boundary condition. Interestingly, within this approach, the problem can be solved exactly in the trumpet, stem-flower and stem regimes. In the trumpet regime, we get $tau=tau_{0}n_d^{2}$ where $n_d$ is the number of monomers that have desorbed at the time $tau$. $tau_{0}$ is known only numerically, but for $f$ close to $f_{c}$, it is found to be $tau_{0}sim f_c/(f^{2/3}-f_{c}^{2/3})$. If one used simple Rouse dynamics, this result changes to { ormalsize $tausim f_c n_d^2/(f-f_{c})$.} In the other regimes too, one can find exact solution, and interestingly, in all regimes $tau sim n_d^2$.
The scission kinetics of bottle-brush molecules in solution and on an adhesive substrate is modeled by means of Molecular Dynamics simulation with Langevin thermostat. Our macromolecules comprise a long flexible polymer backbone with $L$ segments, co
nsisting of breakable bonds, along with two side chains of length $N$, tethered to each segment of the backbone. In agreement with recent experiments and theoretical predictions, we find that bond cleavage is significantly enhanced on a strongly attractive substrate even though the chemical nature of the bonds remains thereby unchanged. We find that the mean bond life time $<tau>$ decreases upon adsorption by more than an order of magnitude even for brush molecules with comparatively short side chains $N=1 div 4$. The distribution of scission probability along the bonds of the backbone is found to be rather sensitive regarding the interplay between length and grafting density of side chains. The life time $<tau>$ declines with growing contour length $L$ as $<tau>propto L^{-0.17}$, and with side chain length as $<tau>propto N^{-0.53}$. The probability distribution of fragment lengths at different times agrees well with experimental observations. The variation of the mean length $L(t)$ of the fragments with elapsed time confirms the notion of the thermal degradation process as a first order reaction.
We examine the phase transition of polymer adsorption as well as the underlying kinetics of polymer binding from dilute solutions on a structureless solid surface. The emphasis is put on the properties of regular multiblock copolymers, characterized
by block size M and total length N as well as on random copolymers with quenched composition p of sticky and neutral segments. The macromolecules are modeled as coarse-grained bead-spring chains subject to a short-ranged surface adhesive potential. Phase diagrams, showing the variation of the critical threshold for single chain adsorption in terms of M and p are derived from scaling considerations in agreement with results from computer experiment. Using both scaling analysis and numerical data from solving a system of coupled Master equations, we demonstrate that the phase behavior at criticality, and the adsorption kinetics may be adequately predicted and understood, in agreement with the results of extensive Monte Carlo simulations. Derived analytic expressions for the mean fraction of adsorbed segments as well as for Probability Distribution Functions of the various structural building blocks (i.e., trains, loops, tails) at time t during the chain attachment process are in good agreement with our numeric experiments and provide insight into the mechanism of polymer adsorption.
The thermally assisted detachment of a self-avoiding polymer chain from an adhesive surface by an external force applied to one of the chain ends is investigated. We perform our study in the fixed height statistical ensemble where one measures the fl
uctuating force, exerted by the chain on the last monomer when a chain end is kept fixed at height $h$ over the solid plane at different adsorption strength $epsilon$. The phase diagram in the $h - epsilon$ plane is calculated both analytically and by Monte Carlo simulations. We demonstrate that in the vicinity of the polymer desorption transition a number of properties like fluctuations and probability distribution of various quantities behave differently, if $h$ rather than $f$ is used as an independent control parameter.
We show that the structural properties and phase behavior of a self-avoiding polymer chain on adhesive substrate, subject to pulling at the chain end, can be obtained by means of a Grand Canonical Ensemble (GCE) approach. We derive analytical express
ions for the mean length of the basic structural units of adsorbed polymer, such as loops and tails, in terms of the adhesive potential and detachment force, and determine values of the universal exponents which govern their probability distributions. Most notably, the hitherto controversial value of the critical adsorption exponent $phi$ is found to depend essentially on the interaction between different loops. The chain detachment transition turns out to be of the first order, albeit dichotomic, i.e., no coexistence of different phase states exists. These novel theoretical predictions and the suggested phase diagram of the adsorption-desorption transformation under external pulling force are verified by means of extensive Monte Carlo simulations.
We study analytically and by means of an off-lattice bead-spring dynamic Monte Carlo simulation model the adsorption kinetics of a single macromolecule on a structureless flat substrate in the regime of strong physisorption. The underlying notion of
a ``stem-flower polymer conformation, and the related mechanism of ``zipping during the adsorption process are shown to lead to a Fokker-Planck equation with reflecting boundary conditions for the time-dependent probability distribution function (PDF) of the number of adsorbed monomers. The theoretical treatment predicts that the mean fraction of adsorbed segments grows with time as a power law with a power of $(1+ u)^{-1}$ where $ uapprox 3/5$ is the Flory exponent. The instantaneous distribution of train lengths is predicted to follow an exponential relationship. The corresponding PDFs for loops and tails are also derived. The complete solution for the time-dependent PDF of the number of adsorbed monomers is obtained numerically from the set of discrete coupled differential equations and shown to be in perfect agreement with the Monte Carlo simulation results. In addition to homopolymer adsorption, we study also regular multiblock copolymers and random copolymers, and demonstrate that their adsorption kinetics may be considered within the same theoretical model.
The adsorption of a single multi-block $AB$-copolymer on a solid planar substrate is investigated by means of computer simulations and scaling analysis. It is shown that the problem can be mapped onto an effective homopolymer adsorption problem. In p
articular we discuss how the critical adsorption energy and the fraction of adsorbed monomers depend on the block length $M$ of sticking monomers $A$, and on the total length $N$ of the polymer chains. Also the adsorption of the random copolymers is considered and found to be well described within the framework of the annealed approximation. For a better test of our theoretical prediction, two different Monte Carlo (MC) simulation methods were employed: a) off-lattice dynamic bead-spring model, based on the standard Metropolis algorithm (MA), and b) coarse-grained lattice model using the Pruned-enriched Rosenbluth method (PERM) which enables tests for very long chains. The findings of both methods are fully consistent and in good agreement with theoretical predictions.
