We introduce a reweighting scheme for the path ensembles in the transition interface sampling framework. The reweighting allows for the analysis of free energy landscapes and committor projections in any collective variable space. We illustrate the reweighting scheme on a two dimensional potential with a nonlinear reaction coordinate and on a more realistic simulation of the Trp-cage folding process. We suggest that the reweighted path ensemble can be used to optimize possible nonlinear reaction coordinates.

1.
G. M.
Torrie
and
J. P.
Valleau
,
Chem. Phys. Lett.
28
,
578
(
1974
).
2.
H.
Grubmüller
,
Phys. Rev. E
52
,
2893
(
1995
).
3.
A.
Laio
and
M.
Parrinello
,
Proc. Natl. Acad. Sci. U.S.A.
99
,
12562
(
2002
).
4.
D.
Frenkel
and
B.
Smit
,
Understanding Molecular Simulation
, 2nd ed. (
Academic
,
San Diego, CA
,
2002
).
5.
C. H.
Bennett
, in
Algorithms for Chemical Computations
,
ACS Symposium Series
Vol.
46
, edited by
R.
Christofferson
(
American Chemical Society
,
Washington, D.C.
,
1977
).
6.
D.
Chandler
,
J. Chem. Phys.
68
,
2959
(
1978
).
7.
A. F.
Voter
,
J. Chem. Phys.
106
,
4665
(
1997
).
8.
A. F.
Voter
,
Phys. Rev. Lett.
78
,
3908
(
1997
).
9.
M. R.
Sørensen
and
A. F.
Voter
,
J. Chem. Phys.
112
,
9599
(
2000
).
10.
A. F.
Voter
,
Phys. Rev. B
57
,
R13985
(
1998
).
11.
C.
Dellago
,
P. G.
Bolhuis
, and
P. L.
Geissler
,
Adv. Chem. Phys.
123
,
1
(
2002
).
12.
C.
Dellago
,
P. G.
Bolhuis
,
F. S.
Csajka
, and
D.
Chandler
,
J. Chem. Phys.
108
,
1964
(
1998
).
13.
T. S.
van Erp
,
D.
Moroni
, and
P. G.
Bolhuis
,
J. Chem. Phys.
118
,
7762
(
2003
).
14.
R. J.
Allen
,
D.
Frenkel
, and
P. R.
ten Wolde
,
J. Chem. Phys.
124
,
024102
(
2006
).
15.
D.
Moroni
,
T. S.
van Erp
, and
P. G.
Bolhuis
,
Phys. Rev. E
71
,
056709
(
2005
).
16.
A. M.
Ferrenberg
and
R. H.
Swendsen
,
Phys. Rev. Lett.
63
,
1195
(
1989
).
17.
T. S.
van Erp
,
Phys. Rev. Lett.
98
,
268301
(
2007
).
18.
P. G.
Bolhuis
,
J. Chem. Phys.
129
,
114108
(
2008
).
19.
D.
Minh
and
J.
Chodera
,
J. Chem. Phys.
131
,
134110
(
2009
).
20.
W.
Lechner
,
J.
Rogal
,
J.
Juraszek
,
B.
Ensing
, and
P.
Bolhuis
,
J. Chem. Phys.
133
,
174110
(
2010
).
21.
J.
Juraszek
and
P. G.
Bolhuis
,
Proc. Natl. Acad. Sci. U.S.A.
103
,
15859
(
2006
).
22.
J.
Juraszek
and
P. G.
Bolhuis
,
Biophys. J.
95
,
4246
(
2008
).
23.
These Markovian probability distributions depend on the dynamics. For instance, they are delta functions for Hamiltonian dynamics. For the Langevin dynamics employed in this paper they are given by bivariate Gaussian expressions (Ref. 12).
24.
P. G.
Bolhuis
,
C.
Dellago
, and
D.
Chandler
,
Proc. Natl. Acad. Sci. U.S.A.
97
,
5877
(
2000
).
25.
B.
Peters
and
B. L.
Trout
,
J. Chem. Phys.
125
,
054108
(
2006
).
26.
B.
Peters
,
G. T.
Beckham
, and
B. L.
Trout
,
J. Chem. Phys.
127
,
034109
(
2007
).
27.
W.
E
,
W.
Ren
, and
E.
Vanden-Eijnden
,
Phys. Rev. B
66
,
052301
(
2002
).
28.
E.
Vanden-Eijnden
and
M.
Venturoli
,
J. Chem. Phys.
130
,
194103
(
2009
).
30.
A.
Berezhkovskii
and
A.
Szabo
,
J. Chem. Phys.
122
,
014503
(
2005
).
31.
Note that λ0 usually defines the boundary of stable state A, likewise λn defines B. Often λ0=λ1 and λn1=λn, but this is not required, and in fact not the case in this work. See Ref. 32.
32.
T. S.
van Erp
and
P. G.
Bolhuis
,
J. Comput. Phys.
205
,
157
(
2005
).
You do not currently have access to this content.