Force spectroscopy of polymer desorption: Theory and Molecular Dynamics simulation
Abstract
Forced detachment of a single polymer chain, stronglyadsorbed on a solid substrate, is investigated by two complementary methods: a coarsegrained analytical dynamical model, based on the Onsager stochastic equation, and Molecular Dynamics (MD) simulations with Langevin thermostat. The suggested approach makes it possible to go beyond the limitations of the conventional BellEvans model. We observe a series of characteristic force spikes when the pulling force is measured against the cantilever displacement during detachment at constant velocity (displacement control mode) and find that the average magnitude of this force increases as grows. The probability distributions of the pulling force and the endmonomer distance from the surface at the moment of final detachment are investigated for different adsorption energy and pulling velocity . Our extensive MDsimulations validate and support the main theoretical findings. Moreover, the simulation reveals a novel behavior: for a strongfriction and massive cantilever the force spikes pattern is smeared out at large . As a challenging task for experimental biopolymers sequencing in future we suggest the fabrication of stiff, superlight, nanometersized AFM probe.
I Introduction
In recent years singlemolecule pulling techniques based on the use of laser optical tweezers (LOT) or atomic force microscope (AFM) have gained prominence as a versatile tool in the studies of noncovalent bonds and selfassociating biomolecular systems Ritort (); Franco (); Hugel (); Butt (); Gao (); Janshof (); Carrion (); Fisher (); Gaub (). The latter could be exemplified by the basepair binding in DNA as well as by ligandreceptor interactions in proteins and has been studied recently by means of Brownian dynamic simulations and the master equation approach Alexander_1 (); Alexander_2 (). The LOT and AFM methods are commonly used to manipulate and exert mechanical forces on individual molecules. In LOT experiments, a micronsized polystyrene or silica bead is trapped in the focus of the laser beam by exerting forces in the range . Typically, AFM (which covers forces interval in range) is ideal for investigations of relatively strong inter or intramolecular interactions which are involved in pulling experiments in biopolymers such as polysaccharides, proteins and nucleic acids. On the other hand, due to the relatively small signaltonoise ratio, the AFM experiments have limitations with regard to the mechanochemistry of weak interactions in the lower piconewton regime.
The method of dynamic force spectroscopy (DFS) is used to probe the forceextension relationship, rupture force distribution, and the force vs loading rate dependence for singlemolecule bonds or for more complicated multiplybonded attachments. Historically, the first theoretical interpretation of DFS has been suggested in the context of single cell adhesion by Bell Bell () and developed by Evans Evans_1 (); Evans_2 (); Evans_3 (). The consideration has been based on the semiphenomenological Arrhenius relation which describes surface detachment under timedependent pulling force, , with being the loading rate. It was also assumed that the effective activation energy, , may be approximated by a linear function of the force, i.e., . Here is the distance between the bonded state and the transition state where the activation barrier is located. The resulting BellEvans (BE) equation then gives the mean detachment force as a function of temperature and loading rate , i.e., , where is the desorption rate in the absence of applied pulling force.
As one can see from this BEequation, the simple surmounting of BEactivation barrier results in a linear dependence of detachment force on the logarithm of loading rate, provided one uses the applied force as a governing parameter in the detachment process (i.e., working in an isotensional ensemble when is controlled and the distance from the substrate to the clamped endmonomer of the polymer chain fluctuates). For multiplybonded attachments the interpretation problem based on this equation becomes more complicated since a nonlinear relationship is observed Merkel (). In this case chain detachment involves passages over a cascade of activation barriers. For example, Merkel et al. Merkel () suggested that the net rate of detachments can be approximated by a reciprocal sum of characteristic times, corresponding to jumps over the single barriers. In particular, regarding the detachment of biotinstreptavidin single bonds, it was suggested that two consecutive barriers might be responsible for the desorption process.
A simple example of multiplybonded bioassembly is presented by a singestranded DNA (ssDNA) macromolecule, strongly adsorbed on graphite substrate. The forcedinduced desorption (or peeling) of this biopolymer has been studied analytically and by means of Brownian dynamics (BD) simulation by Jagota et al. Jagota_1 (); Jagota_2 (); Jagota_3 (); Jagota_4 (). In ref. Jagota_1 () the equilibrium statistical thermodynamics of ssDNA forcedinduced desorption under force control (FC) and displacement control (DC) has been investigated. In the latter case one works in an isometric ensemble where is controlled and fluctuates. It has been demonstrated that the force response under DC exhibits a series of spikes which carry information about the underlying base sequence of ssDNA. The Brownian dynamics (BD) simulations Jagota_2 () confirmed the existence of such force spikes in the forcedisplacement curves under DC.
The nonequilibrium theory of forced desorption has been developed by Kreuzer et al. Kreuzer_1 (); Kreuzer_2 (); Kreuzer_3 () on the basis of Master Equation approach for the cases of constant velocity and forceramp modes in an AFMexperiment. The authors assumed that individual monomers detachments represent a fast process as compared to the removal of all monomers. This justifies a twostate model where all monomers either remain on the substrate or leave it abruptly. The corresponding transition rates (which constitute a necessary input in the Master Equation approach) must satisfy detailed balance. As a result of the Master Equation solution, the authors obtained a probability distribution of detachment heights (i.e., distances between the cantilever tip and the substrate) as well as an average detachment height as a function of the pulling velocity.
Irrespective of all these efforts, a detailed theoretical interpretation of the dynamic force spectroscopy experiments is still missing. For example, in terms of Kramers reactionrate theory Hanggi () the Arrheniuslike BE  model holds only when the effective activation energy . On the other hand, it is clear that for large forces (which we experience in AFM), the case when occurs fairly often. In this common case the general approach, based on the BEmodel, becomes questionable. Besides, it can be shown Hanke (), that the activation energy vs force dependence, , is itself a nonlinear function, so that the conventional BE  model, based on the linear approximation, , should be limited to small forces. Moreover, the Arrhenius  like relationship for the detachment rate, which was used in the BE  model, is a consequence of a saddlepoint approximation for the stationary solution of Fokker Planck equation Hanggi (). This contradicts the typical loading regimes, used in experiments, where applied force or distance grow linearly with time.
The present paper is devoted to the theoretical investigation of a single molecule desorption dynamics and aimed at interpretation of AFM  or LOT  based dynamic force spectroscopy in the DC constantvelocity mode. The organization of the paper is as follows: in Sec. II we give the equilibrium theory of detachment for the case of strong polymer adsorption. The mean force (measured at the cantilever tip) versus displacement diagram is discussed in detail. In particular, the characteristic force“spikes” structure (which was first discussed in Ref. Jagota_1 (); Jagota_2 ()) can be clearly seen. In Sec. III we give a dynamical version of the detachment process. Our approach rests on construction of general free energy functions, depending on coarsegrained variables, which govern the nonlinear response and structural bonding changes in presence of external forces. The corresponding freeenergybased stochastic equations (known as Onsager equations Puri ()) are derived and solved numerically. This solution makes it possible to provide not only forcedisplacement diagrams and the ensuing dependence on cantilever displacement velocity but also the detachment force probability distribution function (PDF). In Sec. IV the main theoretical results are then checked against extensive Molecular Dynamics (MD) simulation. A brief discussion of results is offered in Sec. V.
Ii Equilibrium theory at the strong adsorption case
Recently we suggested a theory of the forceinduced polymer desorption (for relatively weak adsorption energy) in the isotensional Bhattacharya_1 (); Bhattacharya_2 () and isometric Bhattacharya_3 () equilibrium ensembles supported by extensive Monte Carlo (MC) simulations. In the former case, the fraction of adsorbed monomers changes abruptly (undergoes a jump) when one varies the adsorption energy or the external pulling force. In the second case, the order parameter varies steadily with changing height of the AFMtip, even though the phase transition is still of first order. The total phase diagram in terms of adsorption energy  pulling force, or, adsorption energy  endmonomer height, has been discussed theoretically and in terms of MCsimulations.
On the other hand, the AFM experiments deal with relatively strong forces ( Ritort ()) so that in the case of a single molecule desorption experiment only a really strong adsorption energy is essential. This limit has been discussed in the recent papers by Jagota et al. Jagota_1 (); Jagota_2 (); Jagota_3 (); Jagota_4 () and Kreuzer ey al. Kreuzer_1 (); Kreuzer_2 (); Kreuzer_3 (). Here we consider this problem in a slightly more general form. In so doing we distinguish between two different models: with frictionless and strongfriction substrates, as indicated in Fig. 1.
ii.1 Frictionless substrate
This case has been considered in Refs. Jagota_1 (); Jagota_2 (); Kreuzer_1 (); Kreuzer_2 (); Kreuzer_3 () and is based on the assumption that the force resisting sliding is sufficiently small, i.e., the cantilever tip and the contact point are both placed along the same axis (see Fig. 1 (left panel)). The total partition function for a fixed cantilever distance , i.e., , is a product of partition functions of the adsorbed part , , of the desorbed portion (a stretched polymer portion), , and of the cantilever itself, , where is the number of desorbed polymer segments, and denotes the distance between the clamped end of this desorbed portion and the substrate. As a result,
(1) 
where the integration interval, , and the stepfunction, , imply that restrictions, and , should be applied simultaneously. In this representation is the control variable (which is monitored by the corresponding AFM operating mode) whereas and are coarsegrained dynamic variables which should be integrated (in our case, an integral over , and summation over ) out. Moreover, if we introduce the function
(2) 
then Eq. (1) can be rewritten as
(3) 
In the strong adsorption regime, attains a simple form
(4) 
where the dimensionless adsorption energy . The cantilever manifests itself as a harmonic spring with a spring constant , i.e., the corresponding partition function reads
(5) 
Finally, we derive the partition function of the desorbed part of the polymer as function of the dynamic variables and , based of the Freely Jointed Bond Vector (FJBV) model Lai (); Schurr (). The corresponding Gibbs free energy (i.e., the free energy in the isotensionalensemble) is
(6) 
where the dimensionless force . The corresponding distance .
(7) 
where the so called Langevin function has been used. In the isometricensemble, the proper thermodynamic potential is the Helmholtz free energy, , which is related to by Legendre transformation,
(8) 
where . Taking the Gibbs free energy, Eq. (6), into account and the relation Eq. (7) for the Helmholtz free energy, we have
(9)  
where the function . As a result, Eq. (9) along with Eq. (7) parametrically define , and the corresponding partition function
(10) 
as function of and .
By making use of Eqs. (4), (5), (10), the total partition function given by Eq. (1) reads
(11) 
In Eq. (11) the force should be expressed in terms of as follows: , where denotes the inverse Langevin function. The corresponding effective free energy function in terms of and reads
(12) 
In the limit of a very stiff cantilever, , the cantilever partition function approaches a function Kreuzer_1 ():
(13) 
and Eq. (11) takes the form
(14) 
where and the stepfunction ensures that the condition holds. It is this very stiff cantilever limit that was considered in ref. Jagota_1 (); Jagota_2 ().
For the isometric ensemble, i.e., in the ensemble, the average force , measured by AFMexperiment, is given by
(15)  
where is given by Eq. (11).
The numerical results, which follow from Eq. (15), are shown in Fig. 2. One can immediately see the ”sawtooth“, or forcespikes structure on the forcedisplacement diagram as it was also found by Jagota et al. Jagota_1 () in the limit of very stiff cantilever. Physically, spikes correspond to the reversible transitions , during which the release of polymer stretching energy is balanced by the adsorption energy. The corresponding thermodynamic condition reads . This condition also leads to the spikes amplitude law Jagota_1 (), i.e. the spikes amplitude gradually decreases in the process of chain detachment (i.e., with growing ).
This structure is more pronounced at larger adsorption energy and cantilever spring constant . Thus, while the force oscillates, its mean value remains nearly constant in a broad interval of distances , exhibiting a kind of plateau. Complementary information (for fixed at different values of ) is given on Fig. 3. One can verify that the plateau height is mainly determined by whereas the spikes amplitude is dictated by the cantilever spring constant .
ii.2 Strong polymersubstrate friction
In this limit one has to take into account the specific geometry of an AFM experiment, shown in Fig. 1 (right panel). For simplicity, an infinite friction of the polymer at the surface is assumed. The adsorbed polymer portion may be considered as a twodimensional selfavoiding chain comprising segments. The last contact point (marked as in Fig. 1) can move due to adsorption or desorption elementary events. In Ref. Netz () this was classified as the sticky case. In Fig. 1, is the distance from the cantilever base to to the substrate, is the height of the cantilever tip above the substrate, and is the distance between the cantilever tip and the contact point . Eventually, is the lateral distance between cantilever base and the contact point . One may assume that initially the desorbed portion of segments has occupied a distance of which, due to selfavoiding configurations of an adsorbed chain, equals (where ).
The specific geometry of the AFM experiment in the case of strong polymersubstrate friction (shown in Fig. 1) brings about changes only in the cantilever partition function, i.e., instead of Eq. (5), one has
(16)  
As a result, the total partition function in this case is given by
(17) 
where again the variable should be excluded in favor of by means of the relation . In Eq. (17) the following constraints
(18) 
have been taken into account.
The corresponding free energy functional in terms of dynamical variables and has the following form
(19) 
The average force, which is measured in AFMexperiments, is given by
(20)  
Iii Dynamics of desorption
In our recent paper Paturej () we have studied a single polymer forceinduced desorption kinetics by making use of the notion of tensile blobs as well as by means of Monte Carlo and Molecular Dynamics simulations. It was clearly demonstrated that the total desorption time scales with polymer length as .
In order to treat a realistic AFM experiment in which the cantileversubstrate distance changes with constant velocity , i.e., , one has to consider the AFM tip dynamics. With this in mind, we will develop a coarsegrained stochastic model based on the freeenergy functional Eq. (12). Before proceeding any further, we need to define the adsorptiondesorption potential profile . This plays the role of the potential of mean force (PMF) which depends on .
iii.1 Stochastic Model
In the Helmholtz freeenergy functional , given by Eq. (12) and Eq. (19), the free energy of the adsorbed portion is given by a simple contact potential, , where is an integer number in the range . Considering desorption dynamics (see below), we will treat as a continuous variable with a corresponding adsorptiondesorption energy profile satisfying the following conditions:

For integer values the energy profile has minima whereby we use the contact potential .

For halfinteger values of the adsorption potential goes over maxima.

The activation barrier for monomer desorption, , and the corresponding adsorption activation barrier, , are proportional to the adsorption strength of the substrate whereby .

The adsorptiondesorption energy profile satisfies the boundary conditions: (a fully adsorbed chain), and (an entirely detached chain).
One may show that the following energy profile, given as
(21) 
meets the conditions (1)  (4).
The minima and maxima of Eq. (21) are located in the points defined by with denoting the continuous index of a monomer.
As a result,
(22) 
In Eq. (22) the first term, , is very small and could be neglected. Thus, the minima and maxima are located at the integer and half integer points respectively (see Fig. 4)
In order to calculate the activation barriers, we determine first at the halfinteger points, i.e.,
(23) 
as well as at the integer points
(24) 
Therefore, the activation barriers for the detachment, , and adsorption, , are given by
(25) 
i.e., . Finally, one may readily see that and which is in line with condition (iv).
The total Helmholtz free energy for the frictionless substrate model is given by
(26) 
whereas for the strong polymersubstrate friction model we have
(27) 
These Helmholtz free energy functions govern the dissipative process which is described by the stochastic (Langevin) differential equations
(28) 
where and are the Onsager coefficients. The random forces and describe Gaussian noise with means and correlators given by
(29) 
Equations (28) are usually referred to as the Onsager equations Puri ().
The set of stochastic differential equations, Eq.(28) can be treated by a time integration scheme. Each realization of the solution provides a time evolution of and . In order to get mean values of the observables, these trajectories should be averaged over many independent runs . For example, in order to obtain the average force, Eq. (20), one should average over the runs
(30) 
iii.2 Thermodynamic forces
The thermodynamic forces which arise in Eq. (28), i.e., and , could be calculated explicitly. For example, for the free energy function, Eq. (26), one has
(31)  
where we have used (recall that )
On the other hand, a direct calculation shows that
(32) 
so that
(33) 
Thus, for the force , given by Eq. (31), one has
(34) 
In the strong friction case, Eq. (27) leads to a more complicated expression for the thermodynamic force:
(35) 
For the model, given by Eq. (27), the corresponding force reads
(38) 
Finally, the variable should be expressed in terms of by making use of the relationship , where is the inverse Langevin function. A very good approximation for the inverse Langevin function, published in Ref. Cohen (), is given by
(39) 
iii.3 Quasistationary approximation
It could be shown that for a strongly stretched desorbed portion of the polymer chain, the variable rapidly relaxes to its quasistationary value (see Appendix A). In other words, can quickly adjust to the slow evolution of (governed by the Kramers process). In this quasistationary approximation , and from Eq. (37) one has , so that the following nonlinear equation for emerges
(40) 
This could be represented as
(41) 
i.e., the height is instantaneously coupled to the number of desorbed beads, . Inserting Eq. (39) into Eq. (40), one obtains
(42) 
Eventually, we get a system of socalled semiexplicit differentialalgebraic equations (DAE) Ascher ()
(44) 
In this particular form of DAE one can distinguish between the differential variable and the algebraic variable . Eq. (44) can be solved numerically by making use of an appropriate Runge  Kutta (RK) algorithm, as shown in the Appendix B.
iii.4 Results
We have solved numerically our stochastic model, given by Eq. (43) and Eq. (41), for the case of frictionless substrate. To this end we used the second order Runge  Kutta (RK) algorithm for stochastic differentialalgebraic equations (see Appendix B for more details). The advantage of the stochastic differential equations approach as compared to the Master Equation method Kreuzer_3 () is that the former one gives a more detailed (not averaged) dynamic information corresponding to each individual forcedisplacement trajectory (as is often in an experiment). The result of averaging over runs is shown in Fig. 5 (left).
Fig. 5 (right panel) shows the resulting force  displacement diagram for and different detachment velocities. It it worth noting that the ”sawtooth” pattern can be seen for all investigated detachment velocities ranging between and . For larger velocities the plateau height of the force grows substantially. In other words, the mean detachment force increases as the AFMtip velocity gets higher and the bonds stretching between successive monomers becomes stronger.
We have also studied the detachment force behavior as well as that of the cantilever tip distance from the substrate at the moment of a full detachment, (i.e. when ), by repeating the detachment procedure times and plotting the probability distribution functions (PDF) for different adsorption energies and detachment velocities  Fig. 6. As one can see from Fig. 6a, b, both the average and the dispersion of detachment force grow with which agrees with findings for reversible (i.e., when a broken bond can rebind) bondbreaking dynamics Diezemann (). In contrast, the mean cantilever tip distance variance decreases and its average value increases with growing (cf. Fig. 6c, d).
The average detachment force dependence on cantilever velocity is a widely covered subject in the literature in the context of biopolymers unfolding Kreuzer_4 (); Eom_1 (); Eom_2 () or forced separation of two adhesive surfaces Leckband (); Seifert (); Diezemann (). Figure 7a, which shows the result of our calculations, has the characteristic features discussed also in ref. Leckband (). One observes a well expressed crossover from a shallowslope for relatively small detachment rates to a steepslope region as detachment speed increases. One remarkable feature is that this crossover practically does not depend on the adsorption energy : the curve is merely shifted upwards upon increasing of . Therefore, the crossover is not related to a competition between the Kramers rate and the cantilever velocity but rather accounts for the highly nonlinear chain stretching as the velocity increases. The corresponding detachment distance of the cantilever tip (detachment height), Fig. 7b, reveals a specific sigmoidal shape in agreement with the results based on the Master Equation Kreuzer_3 (). At low velocities of pulling, , when the chain still largely succeeds in relaxing back to equilibrium during detachment, an interesting entropy effect is manifested in Fig. 7b: the (effectively) stiffer coil at leaves the substrate at lower values of than in the case of the colder system, . As the pulling velocity grows, however, this entropic effect vanishes and the departure from the substrate is largely governed by the stretching of the bonds rather than of the coil itself whereby the difference in behavior between and disappears.
Eventually, as it can be seen from Fig. 8, the total detachment (peel) time vs. velocity relationship has a welldefined powerlaw behavior, , with the power , in line with previous theoretical findings Seifert ().
Iv MD simulations
iv.1 The model
In our MDsimulations we use a coarsegrained model of a polymer chain of beads connected by finitely extendable elastic bonds. The bonded interactions in the chain is described by the frequently used KremerGrest potential, . The FENE (finitely extensible nonlinear elastic) potential is given by
(45) 
with and .
In order to allow properly for excluded volume interactions between bonded monomers, the repulsion term is taken as WeeksChandlerAnderson (WCA) potential (i.e., the shifted and truncated repulsive branch of the LennardJones potential,) given by
(46) 
with or for , or , and , . The overall potential has a minimum at bond length . The nonbonded interaction between monomers are taken into account by means of the WCA potential, Eq. 46. Thus, the interactions in our model correspond to good solvent conditions.
The substrate in the present investigation is considered simply as a structureless adsorbing plane, with a LennardJones potential acting with strength in the perpendicular –direction, . In our simulations we consider as a rule the case of strong adsorption , where is a temperature of Langevin thermal bath described below.
The dynamics of the chain is obtain by solving the Langevin equations of motion for the position of each bead in the chain,
(47) 
which describes the Brownian motion of a set of bonded particles.
The influence of solvent is split into slowly evolving viscous force and rapidly fluctuating stochastic force. The random Gaussian force is related to friction coefficient by the fluctuationdissipation theorem. The integration step is and time in measured in units of , where denotes the mass of the polymer beads, . In all our simulations the velocityVerlet algorithm was used to integrate equations of motion (47).
The molecule is pulled by a cantilever at constant velocity . The cantilever is imitated by two beads connected by harmonic spring and attached to one of the ends of the chain. ^{1}^{1}1This setup is different than the one used by S. Iliafar et al. Jagota_3 ()); in their study a harmonic spring was connected to a ”big” monomer (with large friction coefficient) on one side, and to a mobile wall on the other side. In our case the harmonic spring spans two beads.
The mass of beads , forming the cantilever, was set either to or to . The equilibrium size of this harmonic spring was set to and the spring constant was varied in the range . The hydrodynamics radius of beads composing the cantilever was varied by changing the friction coefficient , taking into account the Stokes’ law, , where is the solvent viscosity.
Taking the value of the thermal energy J at K, the typical Kuhn length of nm and the mass of coarsegrained monomer as kg setups the unit of time in our simulations which is given in ps. The velocities used in simulations are in units of . Spring constants of our cantilever in real units are: .
Two typical snapshots of a polymer chain during slow detachment from an adsorbing substrate with different strengths of adsorption, and are shown in Fig. 9. Evidently, the chain is much more stretched for the stronglyattractive substrate where all adsorbed monomers stick firmly to the surface.
iv.2 MDresults
As we have already seen in Sec. III.4, the averaging of the force profile over many runs reveals the inherent sawtoothstructure of the force vs distance dependence (see Fig. 5) which is otherwise overshaded by thermal noise. Our MDsimulation result, depicted in Fig. 10, show the same tendency against the noisy background of a single detachment event. Therefore, for better clarity and physical insight, all our graphic results that are given below result from such averaging procedure.
Figure 11a shows how adsorption energy affects the force vs distance relationship. Apparently, with increasing the mean force (plateau height) is found to grow in agreement with our equilibrium theory results, given in Fig. 3. As suggested by our recent theory Bhattacharya_1 (); Bhattacharya_3 (), the plateau height goes up as , or as , for relatively small or large values, respectively. The amplitude of spikes increases with growing too, in line with the equilibrium findings (see Fig. 3). Moreover, as found by Jagota et al. Jagota_1 (), the amplitude of spikes follows an exponential law, , where is the number of desorbed polymer segments. On the other hand, the comparison of Fig. 11b and Fig.2 suggests that the stiffness of the cantilever spring constant affects mainly the spike amplitude especially at large .
Eventually, we demonstrate the impact of cantilever velocity, , as well as of its mass, , and friction coefficient, , on the forcedistance profile. Apparently, these parameters affect differently strong the observed force  distance relationship. Similar to the results, obtained for our coarsegrained model in Sec. II, in the MDsimulation data the plateau height grows less than twice upon velocity increase of three orders of magnitude (see Fig. 12)! Only for a very massive, (), and strongfriction, (), cantilever, the plateau height grows significantly and gains a slight positive slope (see Fig.12d) whereby oscillations vanish. This occurs for the fastest detachment . Evidently, this effect is related to the combined role of the friction force in the case of rapid detachment along with the much larger inertial force () whereby the substrateinduced oscillations are overshadowed by the increased effort of pulling. In contrast, neither Fig. 12b, nor Fig. 12c indicate any major qualitative changes in the vsbehavior when mediumfriction, or mass cantilever alone are drastically changed.
The PDF of the detachment force and its velocity dependence are shown in Fig. 13. Similarly as in Sec. II, the average value and dispersion grow with increasing speed of pulling and this is weakly sensitive with regard to the adsorption strength of the substrate . Remarkably, the mean detachment force shows a similar nonlinear dependence on (cf. Fig. 7a). The crossover position does not change practically as the adhesion strength is varied, and the variation of the other parameters () towards a massive and strongfriction cantilever render this crossover considerably more pronounced.
The complementary PDF for the detachment height is given in Fig. 14a together with the corresponding average vs relationship. As predicted by our analytic model, cf. Section II, the height of final detachment of the chain from the substrate becomes larger for faster peeling and stronger adhesion , which is consistent with the MD data. One can see again the typical sigmoidalshape in the vs dependence.
The two panels for different temperature, shown in Fig. 14b, indicate a smaller increase in at the higher temperature, provided the pulling velocity is sufficiently small too. This can be readily understood in terms entropic (rubber) elasticity of polymers and represents a case of delicate interplay between entropy and energydominated behavior. It is well known that a polymer coil becomes less elastic (i.e., it contracts) upon a temperature increase, cf. the lowest (grey) curve in Fig. 14b, (left panel) at , so that is smaller than in the corresponding lowest curve for in the right panel of Fig. 14b. This occurs at low values of . On the other hand, the softer chain (at ) stretches more easily and, therefore, goes up to for the highest speed instead of for . This entropic effect is well expressed at weak attraction to the surface, , which does not induce strong stretching of the bonds along the chain backbone. In contrast, at high , the bonds extend so strongly that the chain turns almost into a string and entropy effects become negligible. The energy cost of stretching then dominates and leads to higher values of at the higher temperature (cf. upper most green symbols in Fig. 14b) since it is now the elasticity of the bonds between neighboring segments which governs the physics of detachment. In this case the elastic constant of the bonds effectively decreases with an increase of so that the distance of detachment in the left panel of Fig. 14b for is higher than that for in the right panel.
V Discussion
We have demonstrated in this paper that a simple theory, based on the Onsager stochastic equations, yields an adequate description of a typical AFM experiment within the displacementcontrol mode. This approach makes it possible to relax most of the restrictions inherent in the BEmodel. For example, this approach also holds for small desorption activation barriers (i.e., for ), and also for nonlinear barrier vs. force dependence. It naturally takes into account the reversible desorptionadsorption events Diezemann () which are neglected in BEmodel. Moreover, it does not rest on the stationary approximation (which is customary in the standard Kramers rate calculation Hanggi ()) and is, therefore, ideally suited for description of driven force (FC) or displacementcontrol (DC) regimes. One of the principal results in this analytic treatment is the predicted existence of characteristic spikes the mean force vs distance profile, observed in the DCregime. These spikes depend on the adsorption energy , cantilever spring constant as well as on the cantilever velocity . In equilibrium, this has been found earlier by Jagota and coworkers Jagota_1 (). The PDF of detachment forces and detachment distances are been thoroughly investigated. The relevant mean detachment force is found to be a strongly nonlinear function of which is mainly governed by the nonlinear chain stretching upon increasing . The average full detachment (peeling) time scales which is supported by early theoretical findings Seifert ().
Some of these predictions were checked by means of MDsimulation and found in a qualitative agreement with the results, gained by the analytic method. Most notably, this applies to properties like the characteristic force oscillations pattern and the mean force vs cantilever velocity dependence. On the other hand, our MDsimulation reveals a very strong increase in the magnitude of the force plateau for a strongfriction () and massive () cantilever. Interestingly, in this case the spikes pattern is almost totally smeared out. This might be the reason why the force spikes pattern is not seen in laboratory detachment experiments. We recall that in a recent Brownian dynamic simulation (which totally ignores inertia forces) Jagota_3 (), the friction coefficient of the cantilever was times larger than the friction coefficient of the chain segments. It was shown that for this highfriction cantilever and large velocity of pulling, the force spikes pattern was significantly attenuated Jagota_3 () so that information on the base sequence was hardly assessable. Therefore, fabrication of a stiff and superlight, nanometersized AFM probe would be a challenging task for future developments of biopolymer sequencing.
As an outlook, our coarsegrained Onsager stochastic model could be generalized to encompass investigations of forced unfolding of a multidomain, selfassociating biopolymers Kreuzer_4 ().
Acknowledgement
This work was supported by grant SFB 625 from German Research Foundation (DFG). Computational time on PLGrid Infrastructure is acknowledged. A.M. thanks the Max Planck Institute for Polymer Research in Mainz, Germany, and CECAM  MPI for hospitality during his visit at the Institute.
Appendix A The separation as an instantaneously adjustable variable
Due to strong adsorption, the desorbed portion of polymer chain is expected to be strongly stretched. One could simplify the force,