Jarosław
Paturej
ab,
Johan L. A.
Dubbeldam
c,
Vakhtang G.
Rostiashvili
d,
Andrey
Milchev
de and
Thomas A.
Vilgis
d
aDepartment of Chemistry, University of North Carolina, Chapel Hill, NC 27599, USA. E-mail: paturej@live.unc.edu
bInstitute of Physics, University of Szczecin, Wielkopolska 15, 70451 Szczecin, Poland
cDelft University of Technology, 2628CD Delft, The Netherlands
dMax Planck Institute for Polymer Research, 10 Ackermannweg, 55128 Mainz, Germany
eInstitute of Physical Chemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
First published on 18th February 2014
Forced detachment of a single polymer chain, strongly adsorbed on a solid substrate, is investigated by two complementary methods: a coarse-grained analytical dynamical model, based on the Onsager stochastic equation, and Molecular Dynamics (MD) simulations with a Langevin thermostat. The suggested approach makes it possible to go beyond the limitations of the conventional Bell–Evans model. We observe a series of characteristic force spikes when the pulling force is measured against the cantilever displacement during detachment at constant velocity vc (displacement control mode) and find that the average magnitude of this force increases as vc increases. The probability distributions of the pulling force and the end-monomer distance from the surface at the moment of the final detachment are investigated for different adsorption energies ε and pulling velocities vc. Our extensive MD simulations validate and support the main theoretical findings. Moreover, the simulations reveal a novel behavior: for a strong-friction and massive cantilever the force spike pattern is smeared out at large vc. As a challenging task for experimental bio-polymer sequencing in future we suggest the fabrication of a stiff, super-light, nanometer-sized AFM probe.
The method of dynamic force spectroscopy (DFS) is used to probe the force–extension relationship, rupture force distribution, and the force vs. loading rate dependence for single-molecule bonds or for more complicated multiply bonded attachments. Historically, the first theoretical interpretation of DFS has been suggested in the context of single cell adhesion by Bell12 and developed by Evans.13–15 The consideration has been based on the semi-phenomenological Arrhenius relationship which describes surface detachment under time-dependent pulling force, f = rlt, with rl being the loading rate. It was also assumed that the effective activation energy, Eb(f), may be approximated by a linear function of the force, i.e., Eb(f) = E(0)b − xβf. Here xβ is the distance between the bonded state and the transition state where the activation barrier is located. The resulting Bell–Evans (BE) equation then gives the mean detachment force as a function of temperature T and loading rate rl, i.e., where κ0 is the desorption rate in the absence of applied pulling force.
As one can see from this BE equation, the simple surmounting of the BE-activation barrier results in a linear dependence of the detachment force on the logarithm of the loading rate, provided one uses the applied force as a governing parameter in the detachment process (i.e., working in an isotensional ensemble when f is controlled and the distance D from the substrate to the clamped end-monomer of the polymer chain fluctuates). For multiply bonded attachments the interpretation problem based on this equation becomes more complicated since a non-linear f–lnrl relationship is observed.16 In this case chain detachment involves passage over a cascade of activation barriers. For example, Merkel et al.16 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 biotin–streptavidin single bonds, it was suggested that two consecutive barriers might be responsible for the desorption process.
A simple example of multiply bonded bio-assembly is presented by a singe-stranded DNA (ssDNA) macromolecule, strongly adsorbed on the graphite substrate. The force-induced desorption (or peeling) of this biopolymer has been studied analytically and by means of Brownian dynamics (BD) simulations by Jagota et al.17–20 In ref. 17 the equilibrium statistical thermodynamics of ssDNA force-induced desorption under force control (FC) and displacement control (DC) have been investigated. In the latter case one works in an isometric ensemble where D is controlled and f 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) simulations18 confirmed the existence of such force spikes in the force–displacement curves under DC.
The nonequilibrium theory of forced desorption has been developed by Kreuzer et al.21–23 on the basis of the master equation approach for the cases of constant velocity and force-ramp modes in an AFM experiment. The authors assumed that the detachment of individual monomers represents a fast process as compared to the removal of all monomers. This justifies a two-state 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 a 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 reaction-rate theory24 the Arrhenius-like BE model holds only when the effective activation energy Eb(f) ≫ kBT. On the other hand, it is clear that for large forces (which we experience in AFM), the case when Eb(f) ≈ kBT occurs fairly often. In this common case the general approach, based on the BE-model, becomes questionable. Besides, it can be shown,25 that the activation energy vs. force dependence, Eb(f), is itself a nonlinear function, so that the conventional BE-model, based on the linear approximation, Eb(f) ≈ E(0)b − xβf, 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 saddle-point approximation for the stationary solution of the Fokker–Planck equation.24 This contradicts the typical loading regimes, used in experiments, where the applied force or distance grows 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 constant-velocity mode. The organization of the paper is as follows: in Section 2 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. 17 and 18) can be clearly seen. In Section 3 we give a dynamical version of the detachment process. Our approach rests on construction of general free energy functions, depending on coarse-grained variables, which govern the non-linear response and structural bonding changes in the presence of external forces. The corresponding free-energy-based stochastic equations (known as Onsager equations26) are derived and solved numerically. This solution makes it possible to provide not only force–displacement diagrams and the ensuing dependence on the cantilever displacement velocity vc but also the detachment force probability distribution function (PDF). In Section 4 the main theoretical results are then checked against extensive Molecular Dynamics (MD) simulations. A brief discussion of results is offered in Section 5.
On the other hand, the AFM experiments deal with relatively strong forces (20 pN–10 nN (ref. 1)) 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.17–20 and Kreuzer et al.21–23 Here we consider this problem in a slightly more general form. In doing so we distinguish between two different models: with frictionless- and strong-friction substrates, as indicated in Fig. 1.
(1) |
(2) |
(3) |
In the strong adsorption regime, Ξads(n) attains a simple form
Ξads(n) = exp[ε(N − n)], | (4) |
(5) |
Finally, we derive the partition function of the desorbed part of the polymer as a function of the dynamic variables n and R, based on the Freely Jointed Bond Vector (FJBV) model.30,31 The corresponding calculations using the Legendre transformation are relegated to Appendix A. Thus, the polymer partition function takes the form
Ξpol(n, R) = exp[n()], | (6) |
By making use of eqn (4)–(6), the total partition function given by eqn (1) reads
(7) |
The corresponding effective free energy function in terms of n and R reads
(8) |
In the limit of a very stiff cantilever, kcb2/kBT ≫ 1, the cantilever partition function approaches a δ-function:21
(9) |
(10) |
For the isometric ensemble, i.e., in the D-ensemble, the average force 〈fz〉, measured by the AFM experiment, is given by
(11) |
The numerical results, which follow from eqn (11), are shown in Fig. 2. One can immediately see the “sawtooth”-, or force-spikes structure in the force–displacement diagram as it was also found by Jagota et al.17 in the limit of a very stiff cantilever. Physically, spikes correspond to the reversible transitions n ⇄ n + 1, during which the release of the polymer stretching energy is balanced by the adsorption energy. The corresponding thermodynamic condition reads (n, R) = (n + 1, R). This condition also leads to the spike amplitude law famp ∝ exp(ε/n),17i.e. the spike amplitude gradually decreases in the process of chain detachment (i.e., with growing n – cf.Fig. 2).
Fig. 2 The equilibrium force–displacement diagrams calculated according to eqn (11). The sawtooth structure becomes more pronounced with increasing adsorption energy ε and spring constant kc. |
This structure is more pronounced at larger adsorption energy ε and cantilever spring constant kc. Thus, while the force oscillates, its mean value remains nearly constant in a broad interval of distances D, exhibiting a kind of plateau. Complementary information (for fixed kc at different values of ε) is given in Fig. 3. One can verify that the plateau height is mainly determined by ε whereas the spike amplitude is dictated by the cantilever spring constant kc.
Fig. 3 The equilibrium force–displacement diagrams calculated according to eqn (11). The same as in Fig. 2 but for fixed kc and different ε. |
It should be noted that similar spikes have been indeed observed in experiments on single-molecule stretching of proteins due to unfolding of the multidomain biopolymer structure.21 As far as the size of each such domain is much larger than that of a single segment, resolving this “sawtooth” behavior is significantly easier than in the present case of monomer peeling.
The specific geometry of the AFM experiment in the case of strong polymer–substrate friction (shown in Fig. 1) brings about changes only in the cantilever partition function, i.e., instead of eqn (5), one has
(12) |
As a result, the total partition function in this case is given by
(13) |
(14) |
The corresponding free energy functional in terms of dynamical variables n and R has the following form
(15) |
The average force, which is measured in AFM experiments, is given by
(16) |
In order to treat a realistic AFM experiment in which the cantilever–substrate distance changes with constant velocity vc, i.e., D(t) = D0 + vct, one has to consider the AFM tip dynamics. With this in mind, we will develop a coarse-grained stochastic model based on the free-energy functional eqn (8). Before proceeding any further, we need to define the adsorption–desorption potential profile Fads(n). This plays the role of the potential of mean force (PMF) which depends on n.
1. For integer n-values the energy profile has minima whereby we use the contact potential Fads(n) = −kBTε(N − n).
2. For half-integer values of n the adsorption potential goes over maxima.
3. The activation barrier for monomer desorption, ΔE+ = Fads(n + 1/2) − Fads(n), and the corresponding adsorption activation barrier, ΔE− = Fads(n + 1/2) − Fads(n + 1), are proportional to the adsorption strength ε of the substrate whereby ΔE+ > ΔE−.
4. The adsorption–desorption energy profile satisfies the boundary conditions: Fads(0) = −kBTεN (a fully adsorbed chain) and Fads(N) = 0 (an entirely detached chain).
One can show (see Appendix B for more details) that the following energy profile, given as (cf.Fig. 4)
Fads(n) = Tε{1 + cos[(2n + 1)π] + n} − kBTεN, | (17) |
The total Helmholtz free energy for the frictionless substrate model is given by
(18) |
(19) |
These Helmholtz free energy functions govern the dissipative process, which is described by the stochastic (Langevin) differential equations
(20) |
(21) |
Eqn (20) are usually referred to as the Onsager equations.26
The set of stochastic differential equations, eqn (20), can be treated by a time integration scheme. Each realization (l) of the solution provides a time evolution of n(l)(t) and R(l)(t). In order to obtain mean values of the observables, these trajectories should be averaged over many independent runs l = 1, 2, ߪ, . For example, in order to obtain the average force, eqn (16), one should average over the runs
(22) |
(23) |
For the strong friction case, eqn (19) leads to a more complicated expression for the thermodynamic force:
(24) |
The force fR for the frictionless substrate model is given by (see Appendix C for more details):
(25) |
For the model, given by eqn (19), the corresponding force reads
(26) |
Finally, the variable should be expressed in terms of R/bn by making use of the relationship = −1(R/bn), where −1(x) is the inverse Langevin function. A very good approximation for the inverse Langevin function, published in ref. 34, is given by
(27) |
(28) |
This could be represented as
(29) |
i.e., the height R is instantaneously coupled to the number of desorbed beads, n. Inserting eqn (27) into eqn (28), one obtains
(30) |
The Onsager equation for the slow variable n is given as
(31) |
Eventually, we get a system of so-called semi-explicit differential-algebraic equations (DAE)35
(32) |
In this particular form of DAE one can distinguish between the differential variable n(t) and the algebraic variable R(t). Eqn (32) can be solved numerically by making use of an appropriate Runge–Kutta (RK) algorithm, as shown in Appendix E.
Fig. 5 (right panel) shows the resulting force–displacement diagram for ε = 5 and different detachment velocities. It is worth noting that the “sawtooth” pattern can be seen for all investigated detachment velocities ranging between vc = 5 × 10−4 and vc = 10−2. For larger velocities the plateau height of the force increases substantially. In other words, the mean detachment force increases as the AFM-tip velocity increases and the bond stretching between successive monomers becomes stronger.
We have also studied the detachment force behavior as well as the cantilever tip distance from the substrate at the moment of a full detachment (i.e. when n = N), by repeating the detachment procedure 104 times and plotting the probability distribution functions (PDF) for different adsorption energies ε and detachment velocities vc – Fig. 6. As one can see from Fig. 6a and b, both the average and the dispersion of detachment force grow with vc which agrees with findings for reversible (i.e., when a broken bond can rebind) bond-breaking dynamics.36 In contrast, the mean cantilever tip distance R variance decreases and its average value increases with growing vc (cf.Fig. 6c and d).
The average detachment force dependence on cantilever velocity vc is a widely covered subject in the literature in the context of biopolymer unfolding37–39 or forced separation of two adhesive surfaces.36,40,41Fig. 7a, which shows the results of our calculations, shows the characteristic features discussed also in ref. 40. One observes a well expressed crossover from a shallow-slope for relatively small detachment rates to a steep-slope region as the 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 ε. 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 vc increases. The corresponding detachment distance of the cantilever tip R (detachment height), Fig. 7b, reveals a specific sigmoidal shape in agreement with the results based on the master equation.23 At low velocities of pulling, vc, 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 T = 1.0 leaves the substrate at lower values of R than in the case of the colder system, T = 0.1. As the pulling velocity increases, 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 T = 1.0 and T = 0.1 disappears.
Eventually, as can be seen from Fig. 8, the total detachment (peel) time τdetvs. velocity vc relationship has a well-defined power-law behavior, τdet ∼ 1/vcα, with the power α ≈ 1, in line with previous theoretical findings.41
Fig. 8 The average detachment (peeling) time τdetversus detachment velocity vc for two different adsorption energies ε = 10 and ε = 16. The inversely proportional dependence, τdet ∼ 1/vc, agrees well with previous theoretical findings.41 |
(33) |
In order to allow properly for excluded volume interactions between bonded monomers, the repulsion term is taken as the Weeks–Chandler–Anderson (WCA) potential (i.e., the shifted and truncated repulsive branch of the Lennard-Jones potential) given by
VWCA(r) = 4εLJ[(σ/r)12 − (σ/r)6 + 1/4]θ(21/6σ − r) | (34) |
The substrate in the present investigation is considered simply as a structureless adsorbing plane, with a Lennard-Jones potential acting with strength εs in the perpendicular z-direction, VLJ(z) = 4εs[(σ/z)12 − (σ/z)6]. In our simulations we consider as a rule the case of strong adsorption εs/kBT = 5–20, where kBT is a temperature of the Langevin thermal bath described below.
The dynamics of the chain is obtained by solving the Langevin equations of motion for the position rn = [xn, yn, zn] of each bead in the chain,
mn = Fnj + FnWCA − γṙn + Rn(t)(1, …, N) | (35) |
The influence of solvent is split into slowly evolving viscous force and rapidly fluctuating stochastic force. The random Gaussian force Rn is related to friction coefficient γ = 0.25 τ−1 by the fluctuation–dissipation theorem. The integration step is τ = 0.005 and time is measured in units of where m denotes the mass of the polymer beads, m = 1. In all our simulations the velocity-Verlet algorithm was used to integrate equations of motion (35).
The molecule is pulled by a cantilever at constant velocity V = [0, 0, vc]. The cantilever is imitated by two beads connected by a harmonic spring and attached to one of the ends of the chain.†
The mass of beads mc, forming the cantilever, was set either to mc = 1 or to 25. The equilibrium size of this harmonic spring was set to 0 and the spring constant was varied in the range kc = 50–400εLJ/σ2. The hydrodynamic radius a of beads composing the cantilever was varied by changing the friction coefficient γc = 0.25–25 τ−1, taking into account the Stokes' law, γc = 6πηa, where η is the solvent viscosity.
Taking the value of the thermal energy kBT ≈ 4.11 × 10−21 J at kBT = 300 K, the typical Kuhn length of σ = 1 nm and the mass of the coarse-grained monomer as m ≈ 10−25 kg sets the unit of time in our simulations which is given in 10−12 s = 1 ps. The velocities used in simulations are in units of 10−4 ÷ 10−1 nm ps−1 ≈ 10−1 ÷ 102 m s−1. Spring constants of our cantilever in real units are: kc = 50 ÷ 400 kBT nm−2 = 0.2 ÷ 1.6 N m−1.
Two typical snapshots of a polymer chain during slow detachment from an adsorbing substrate with different strengths of adsorption, εs = 2.5 and εs = 20 are shown in Fig. 9. Evidently, the chain is much more stretched for the strongly attractive substrate where all adsorbed monomers stick firmly to the surface.
Fig. 11a shows how adsorption energy εs affects the force f vs. distance D relationship. Apparently, with increasing εs the mean force (plateau height) is found to increase in agreement with our equilibrium theory results, given in Fig. 3. As suggested by our recent theory,27,29 the plateau height goes up as fp ∝ εs1/2, or as fp ∝ εs, for relatively small or large εs values, respectively. The amplitude of spikes increases with growing εs too, in line with the equilibrium findings (see Fig. 3). Moreover, as found by Jagota et al.,17 the amplitude of spikes follows an exponential law, famp ∝ exp(εs/n), where n is the number of desorbed polymer segments. On the other hand, the comparison of Fig. 11b and 2 suggests that the stiffness of the cantilever spring constant kc affects mainly the spike amplitude especially at large εs.
Eventually, we demonstrate the impact of the cantilever velocity, vc, as well as of its mass, mc, and friction coefficient, γc, on the force–distance profile. Apparently, these parameters affect to different extents the observed force–distance relationship. Similar to the results, obtained for our coarse-grained model in Section 2, in the MD simulation data the plateau height increases less than twice upon the velocity increase of three orders of magnitude (see Fig. 12). Only for a very massive (mc = 25) and strong friction (γc = 25), cantilever, the plateau height increases significantly and gains a slight positive slope (see Fig. 12d) whereby oscillations vanish. This occurs for the fastest detachment vc = 0.1σ/τ. 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 (mc = 25) whereby the substrate-induced oscillations are overshadowed by the increased effort of pulling. In contrast, neither Fig. 12b nor Fig. 12c indicates any major qualitative changes in the f- vs. D-behavior when medium-friction or mass cantilever alone is drastically changed.
The PDF of the detachment force and its velocity vc dependence are shown in Fig. 13. Similar to that in Section 2, the average value and dispersion increase with increasing speed of pulling and this is weakly sensitive with regard to the adsorption strength of the substrate εs. Remarkably, the mean detachment force 〈fd〉 shows a similar nonlinear dependence on ln vc (cf.Fig. 7a). The crossover position does not change practically as the adhesion strength is varied, and the variation of the other parameters (mc = 1 → 25, γc = 0.25 → 25) towards a massive and strong-friction cantilever renders this crossover considerably more pronounced.
The complementary PDF for the detachment height R is given in Fig. 14a together with the corresponding average 〈R〉 vs. vc relationship. As predicted by our analytical model, cf. Section 2, the height of the final detachment of the chain from the substrate becomes larger for faster peeling vc and stronger adhesion ε, which is consistent with the MD data. One can see again the typical sigmoidal-shape in the 〈R〉 vs. vc dependence.
The two panels for different temperatures, shown in Fig. 14b, indicate a smaller increase in 〈R〉 at the higher temperature, provided the pulling velocity vc is sufficiently small too. This can be readily understood in terms of entropic (rubber) elasticity of polymers and represents a case of delicate interplay between entropy and energy-dominated 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 T = 1.0, so that R is smaller than in the corresponding lowest curve for T = 0.1 in the right panel of Fig. 14b. This occurs at low values of vc. On the other hand, the softer chain (at T = 0.1) stretches more easily and, therefore, R goes up to ≈90 for the highest speed vc = 10−1 instead of R ≈ 85 for T = 1.0, vc = 10−1. This entropic effect is well expressed at weak attraction to the surface, εs/kBT = 2.25, which does not induce strong stretching of the bonds along the chain backbone. In contrast, at high εs/kBT = 20, 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 R 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 T so that the distance of detachment R in the left panel of Fig. 14b for T = 1.0 is higher than that for T = 0.1 in the right panel.
Some of these predictions were checked by means of MD simulations and found to be in qualitative agreement with the results, gained by the analytic method. Most notably, this applies to properties like the characteristic force oscillation pattern and the mean force vs. cantilever velocity vc dependence. On the other hand, our MD simulations reveal a very strong increase in the magnitude of the force plateau for a strong-friction (γc = 25) and massive (mc = 25) cantilever. Interestingly, in this case the spike pattern is almost totally smeared out. This might be the reason why the force spike pattern is not seen in laboratory detachment experiments. We recall that in recent Brownian dynamic simulations (which totally ignore inertial forces),19 the friction coefficient of the cantilever was 70 times larger than the friction coefficient of the chain segments. It was shown that for this high-friction cantilever and large velocity of pulling, the force spike pattern was significantly attenuated19 so that information on the base sequence was hardly assessable. Therefore, fabrication of a stiff and super-light, nanometer-sized AFM probe would be a challenging task for future developments of biopolymer sequencing.
As an outlook, our coarse-grained Onsager stochastic model could be generalized to encompass investigations of forced unfolding of multi-domain, self-associating biopolymers.37 In doing so one should go beyond the FJBV chain model and take into account the stiffness, which is typical for most of the biopolymers.
(36) |
(37) |
Fpol(n, R) = Gpol(n, ) + fR, | (38) |
(39) |
(40) |
In eqn (40) the first term, (1/2π) arcsin (1/2π) ≈ 0.025, 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 Fads(s) at the half-integer points, i.e.,
Fads(n + 1/2) = kBTε(n + 5/2) − kBTεN, | (41) |
(42) |
Therefore, the activation barriers for the detachment, ΔE+, and adsorption, ΔE−, are given by
(43) |
(44) |
On the other hand, a direct calculation shows that
(45) |
(46) |
(47) |
Taking into account eqn (46), one finally derives eqn (25).
(48) |
This equation can be easily solved and the corresponding solution has the form
(49) |
(50) |
This relaxation time should be compared to the characteristic time, τkram, of the slow variable n(t) which is governed by the Kramers process. According to the semi-phenomenological Bell model,12 the characteristic time of unbonding (that is, desorption in our case) is given by τKram = τ0exp[(ΔE − r0fp)/kBT] where τ0 = ξ0b2/kBT is the segmental time, ΔE = F1 − F2 is the activation energy for single monomer desorption, r0 stands for the width of adsorption potential, and fp is the plateau height. The free energies in the desorbed, F1, and in the adsorbed, F2, states are given by F1 = −kBT lnμ2 and F2 = −kBTε − kBT lnμ3 where μ2 and μ3 are the so-called connective constants in two- and three dimensions respectively.43
As mentioned in Section 4.2, for large adsorption energies ε the dimensionless plateau height . Taking this into account, one could represent τkram in the following form:
(51) |
(52) |
This condition holds for all typical values of the relevant parameters.
(53) |
(54) |
The integral over the deterministic force in eqn (53) within this 2nd order approximation reads
(55) |
This is so-called trapezoidal rule for approximation of the integral. In order to calculate fn(ni, Ri), one should first take the initial value ni, and find Ri through the solution of the nonlinear equation G(ni, Ri, ti) = 0. For the calculation of fn(ni+1, Ri+1), one can use the forward Euler method of order 1, i.e., nEi+1 = ni + hλnfn(ni, Ri) + h1/2(2λnkBT)1/2Zn and REi+1 are obtained as solution of the equation G(nEi+1, REi+1, ti+1) = 0. Here the random variable Zn is Gaussian with the zero mean value and with variance
〈Zn2〉 = 1 | (56) |
As a result, the recursive procedure which relates the i-th and (i + 1)-th grid points can be defined as:
1. For a given initial value of ni, go to eqn (29) or eqn (30) and solve this nonlinear equation (e.g. G(ni, Ri, ti) = 0) with respect to Ri.
2. Compute g1 = fn(ni, Ri).
3. Compute ni+1 and Ri+1 within the Euler approximation, i.e., calculate first nEi+1 = ni + hλnfn(ni, Ri) + h1/2(2λnkBT)1/2Zn and then solve G(nEi+1, REi+1, ti+1) = 0 with respect to REi+1.
4. Compute g2 = fn(nEi+1, REi+1).
5. Compute the corrected ni+1, i.e.
6. Finally, with the value of ni+1, go to item 1 and solve the nonlinear equation G(ni+1, Ri+1, ti+1) = 0 with respect to Ri+1.
Footnote |
† This setup is different from the one used by S. Iliafar et al.19 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. |
This journal is © The Royal Society of Chemistry 2014 |