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
First published on 6th February 2017
Correction for ‘Folding kinetics of a polymer’ by Štěpán Růžička et al., Phys. Chem. Chem. Phys., 2012, 14, 6044–6053.
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(E) = 〈Si(E)〉, where
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
(E) = log〈exp(Si(E))〉 (its derivative is the black dashed line in Fig. 1). The derivative of S
(E) differs substantially from the one of S
(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).
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(E) converges faster to the extended simulations in comparison to S
(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, EB − EA ≈ 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
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.
![]() | ||
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.
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 |