Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Correction: Folding kinetics of a polymer

Bart Vorselaars *, Štěpán Růžička , David Quigley and Michael P. Allen *
Department of Physics, University of Warwick, Coventry CV4 7AL, UK. E-mail: m.p.allen@warwick.ac.uk

Received 5th March 2013 , Accepted 5th March 2013

First published on 6th February 2017


Abstract

Correction for ‘Folding kinetics of a polymer’ by Štěpán Růžička et al., Phys. Chem. Chem. Phys., 2012, 14, 6044–6053.


This correction article replaces the “Addition” for the above article (DOI: 10.1039/C2CP00051B) that was originally published on the 9th May 2013.

In our original article, a convergence problem resulted in an averaging error in computing the entropy from a set of Wang–Landau1 (WL) Monte-Carlo simulations. Here we report corrected results for the freezing temperature of the homopolymer chain as a function of the range of the non-bonded interaction. We find that the previously reported forward-flux sampling2 (FFS) and brute-force (BF) simulation results are in agreement with the revised Wang–Landau calculations. This confirms the utility of FFS for computing crystallisation rates in systems of this kind.

We previously3 found a small but significant discrepancy in the folding–unfolding transition temperature Tf of a polymer chain: WL results were shifted with respect to both FFS and BF results. Here Tf for the WL results follow from an equal-area Maxwell construction of the derivative of the entropy with respect to energy, ∂S(E)/∂E. The entropy was determined by averaging the results of 20 independent WL simulations i, S[S with combining macron](E) = 〈Si(E)〉, where image file: c7cp90027a-t1.tif and Nsim = 20. The derivative of the entropy for each simulation and for the average are displayed in Fig. 1 as solid blue lines and a dash-dotted red line, respectively. The entropy varies substantially between the runs, and hence the averaging method influences the final result. For example, one could average the density of states (DOS) W(E) = exp(S(E)) instead. Here the DOS is normalized, so that its integral over all energies equals 1. The entropy of the DOS averaged over the simulations equals S[W with combining macron](E) = log〈exp(Si(E))〉 (its derivative is the black dashed line in Fig. 1). The derivative of S[W with combining macron](E) differs substantially from the one of S[S with combining macron](E) with respect to the position of the upwards trend in 1/T(E). This also has consequences on Tf: the difference between the two different averages is about 1% for χ = 1.07 using the Maxwell construction (horizontal lines in Fig. 1).


image file: c7cp90027a-f1.tif
Fig. 1 (Fig. 2 revised) The inverse temperature T(E)−1 = ∂S(E)/∂E for 20 independent WL simulations with m = 20 iterations (blue), obtained from entropy averaging (dash-dotted red), from density-of-states averaging (dashed black) and for a simulation with m = 28 iterations (solid yellow), all for χ = 1.07. Accompanying horizontal equal-area lines from a Maxwell construction have the same colour.

The noise in the WL results and the resulting influence of the average procedure has been reduced as follows. We increased the number of WL iterations from 20 to at least 28, following Taylor et al.,4 and we tightened up the flatness criterion: the next iteration is started when each entry in the ‘visit’-histogram is within 10% of the average as opposed to allowing that the minimum deviates at most 20% from the average. There are also a few other minor changes to the algorithm, to be published elsewhere. The result for T(E)−1 from one of these improved simulations is shown in Fig. 1, yellow line. Two other independent runs gave almost indistinguishable results. Observe that in the transition region, near E = −225, the derivative of S[W with combining macron](E) converges faster to the extended simulations in comparison to S[S with combining macron](E), but that in the flat region the fluctuations in the latter are smaller. Nevertheless, the transition temperature from the DOS-average is much closer to the improved results than the entropy-average: the difference is respectively 0.1% and 1% for χ = 1.07. Results for other values of χ show the same tendency. The reason that averaging the DOS gives a better result might be that in our WL simulations the entropy of the lower energies is usually underestimated before convergence. Weighting with the DOS effectively punishes these low-entropy results.

The phase diagram with the revised WL transition temperatures is given in Fig. 2. The maximum discrepancy in Tf between the FFS/BF and WL data has reduced from approximately 4% to 0.4%, and is now within the statistical error of the simulations. The original differences may seem small, but note that a difference in the transition temperature of approximately 2% results in an apparent shift of one of the rate constants by more than 4 orders of magnitude, as was previously the case for χ = 1.06. This is caused by the relatively large energy difference between the folded A and unfolded B state, EBEA ≈ 270 for χ = 1.07 (Fig. 3). As Tf = 0.4964, a change in the transition temperature by 0.3% corresponds to a shift in the relative occupancy of the state A to B by a factor (PT1(EB)/PT1(EA))/(PT2(EB)/PT2(EA)) ≈ 5 (where PT(E) ∝ W(E)exp(−E/T)) and hence a similar change in one of the rate constants, if the reverse rate is assumed constant. This sensitivity might contribute to the observed discrepancy in the rates between simulations and experiments.5


image file: c7cp90027a-f2.tif
Fig. 2 (Fig. 4 revised) Phase diagram of the N = 128 bead polymer chain. Here χ is the ratio of the maximum attraction cut-off to the hard-sphere diameter, and Tf the freezing temperature determined from three different methods.

image file: c7cp90027a-f3.tif
Fig. 3 (Fig. 3 revised) Free energy relative to its minimum (black) at T = TWLf = 0.4964, and energy distribution functions at T = 0.498 for WL (purple), for the various collision dynamics runs in the folded (blue) and unfolded state (red), all for χ = 1.07. The respective averaged distributions are coloured darker. The probability axis is not shown, but each distribution function has equal area. The vertical lines show the subsequent FFS interfaces for A → B (blue) and B → A (red).

The updated chevron plot is given in Fig. 4. For χ = 1.06 the apparent shift in one of the rate constants, assuming the reverse rate is exact, is indicated. It is now approximately half an order of magnitude and within the noise of the WL and FFS results (e.g. discarding the first red data point halves the difference). Upon decreasing the statistical noise other small inaccuracies might be revealed. For example, there is a small discrepancy between WL and FFS in the energy distribution of the folded states (purple vs. dark blue line in Fig. 3). As opposed to a pure 1D potential with two minima, there are multiple folded states which can not be reached from each other within the time scale of our simulations. Therefore one has to switch to ensemble-averaging. The ensemble of configurations that we used as starting points for the crystallized A states is based on the configuration with the lowest energy that we obtained at the last FFS interface in the folding direction EBm. This still gives a slightly too high average energy. An improved scheme might require the first and last FFS interface to be located closer to their respective free energy minimum, and hence knowledge of the latter should be at hand. In this respect the WL and FFS methods complement each other. Nevertheless, as can be judged from the phase diagram, this is only a minor effect.


image file: c7cp90027a-f4.tif
Fig. 4 (Fig. 5 revised) Chevron plots for various χ. The schematic for χ = 1.06 indicates the unfolding rate for it to bring TFFSf into agreement with TWLf.

We conclude that within statistical error the three methods give consistent results and confirm that forward-flux sampling (FFS) is also a viable and accurate method for use in crystallisation studies of this kind.

The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and readers.

References

  1. F. Wang and D. P. Landau, Phys. Rev. Lett., 2001, 86, 2050–2053 CrossRef CAS PubMed.
  2. R. J. Allen, D. Frenkel and P. R. ten Wolde, J. Chem. Phys., 2006, 124, 024102 CrossRef PubMed.
  3. S. Ruzicka, D. Quigley and M. P. Allen, Phys. Chem. Chem. Phys., 2012, 14, 6044–6053 RSC.
  4. M. P. Taylor, W. Paul and K. Binder, J. Chem. Phys., 2009, 131, 114907 CrossRef PubMed.
  5. S. Auer and D. Frenkel, Nature, 2001, 409, 1020–1023 CrossRef CAS PubMed.

Footnote

Current address: School of Mathematics and Physics, University of Lincoln, Lincolnshire LN6 7TS, United Kingdom. E-mail: bvorselaars@lincoln.ac.uk

This journal is © the Owner Societies 2017
Click here to see how this site uses Cookies. View our privacy policy here.