Nonequilibrium fluctuation-dissipation theorem of ... - LY Chen

Nonequilibrium fluctuation-dissipation theorem of ... - LY Chen Nonequilibrium fluctuation-dissipation theorem of ... - LY Chen

jilll.utsa.edu
from jilll.utsa.edu More from this publisher
22.01.2015 Views

144113-4 L. Y. Chen J. Chem. Phys. 129, 144113 2008 Wx 0 ,x 1 , ...,x N + Wx N ,x N−1 , ...,x 0 N−1 =− x n+1 − x n F n+1 − F n . 21 n=0 The assumption of microscopic reversibility means that the right hand side of Eq. 21 is equal to zero. Now one might be attempted to conclude that the right hand side of Eq. 21 vanishes as dt→0 because x n+1 −x n dt and F n+1 −F n dt. But that conclusion is invalid. In fact, it does not vanish as long as the external force Fx,t is not independent of the system’s coordinates. In general, we have Wx 0 ,x 1 , ...,x N + Wx N ,x N−1 , ...,x 0 N−1 =− x n+1 − x n 2 F n , 22 n=0 where F n Fxt n ,t n /x i . The right hand side of Eq. 22 would be zero if x n+1 −x n 2 were proportional to dt 2 , but the stochastic nature of the system makes x n+1 −x n 2 proportional to dt instead of dt 2 . Therefore the right hand side of Eq. 22 does not vanish in the dt→0 limit. It is equal to twice the dissipative work instead. This means that the CFT and JE are invalid unless the dissipative work is vanishingly small. Now one more question needs to be asked: What about the case when the external force is independent of the system’s coordinates In such a case, the right hand side of Eq. 22 does vanish. However, with such kind of external forces, the system may not transition from state A to state B because of the stochastic nature of its dynamics. In the experiments of mechanically unfolding biological molecules at constant speed, the forces acting on the system are indeed functions of the coordinates of the terminus atoms attached to the tips of the atomic force microscope or the optical tweezers. CFT and JE are inapplicable to those experiments unless the pulling speeds are slow enough that the system is always near equilibrium. In summary, a new FDT has been derived from the Brownian dynamics without invoking any assumption. This new FDT is valid in the near-equilibrium or far nonequilibrium regime, wherever the Brownian dynamics is valid. It is also shown that the CFT and the JE are valid only in the near-equilibrium regime where the dissipative work is nearly zero. The author thanks Robert Buckley for a critical reading of the manuscript. He acknowledges support from a NIH grant Grant No. SC3 GM084834 and the Texas Advanced Computing Center. 1 See, for example, J. Liphardt, B. Onoa, S. B. Smith, I. Tinoco, Jr., and C. Bustamante, Science 292, 733 2001; G. Hummer and A. Szabo, Proc. Natl. Acad. Sci. U.S.A. 98, 3658 2001; S. Park, F. Khalili-Araghi, E. Tajkhorshid, and K. Schulten, J. Chem. Phys. 119, 3559 2003; A.Imparato and L. Peliti, Phys. Rev. E 72, 046114 2005; D. Collin, F. Ritort, C. Jarzynski, S. B. Smith, I. Tinoco, and C. Bustamante, Nature London 437, 231 2005. 2 F. M. Ytreberg, R. H. Swendsen, and D. M. Zuckerman, J. Chem. Phys. 125, 184114 2006. 3 C. Jarzynski, Phys. Rev. Lett. 78, 2690 1997. 4 G. E. Crooks, Phys. Rev. E 61, 2361 2000. 5 E. G. D. Cohen and D. Mauzerall, Mol. Phys. 103, 2923 2005. 6 I. Kosztin, B. Barz, and L. Janosi, J. Chem. Phys. 124, 064106 2006. 7 J. M. G. Vilar and J. M. F. Rubi, Phys. Rev. Lett. 100, 020601 2008. 8 L. Y. Chen, J. Chem. Phys. 129, 091101 2008. 9 H. Risken, The Fokker-Planck Equation Springer-Verlag, New York, 1989. 10 J. C. Phillips, R. Braun, W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kale, and K. Schulten, J. Comput. Chem. 26, 1781 2005. Author complimentary copy. Redistribution subject to AIP license or copyright, see http://jcp.aip.org/jcp/copyright.jsp

144113-4 L. Y. <strong>Chen</strong> J. Chem. Phys. 129, 144113 2008<br />

Wx 0 ,x 1 , ...,x N + Wx N ,x N−1 , ...,x 0 <br />

N−1<br />

=− x n+1 − x n F n+1 − F n .<br />

21<br />

n=0<br />

The assumption <strong>of</strong> microscopic reversibility means that the<br />

right hand side <strong>of</strong> Eq. 21 is equal to zero. Now one might<br />

be attempted to conclude that the right hand side <strong>of</strong> Eq. 21<br />

vanishes as dt→0 because x n+1 −x n dt and F n+1 −F n <br />

dt. But that conclusion is invalid. In fact, it does not vanish<br />

as long as the external force Fx,t is not independent <strong>of</strong> the<br />

system’s coordinates. In general, we have<br />

Wx 0 ,x 1 , ...,x N + Wx N ,x N−1 , ...,x 0 <br />

N−1<br />

=− x n+1 − x n 2 F n ,<br />

22<br />

n=0<br />

where F n Fxt n ,t n /x i . The right hand side <strong>of</strong> Eq. 22<br />

would be zero if x n+1 −x n 2 were proportional to dt 2 , but the<br />

stochastic nature <strong>of</strong> the system makes x n+1 −x n 2 proportional<br />

to dt instead <strong>of</strong> dt 2 . Therefore the right hand side <strong>of</strong><br />

Eq. 22 does not vanish in the dt→0 limit. It is equal to<br />

twice the dissipative work instead. This means that the CFT<br />

and JE are invalid unless the dissipative work is vanishingly<br />

small.<br />

Now one more question needs to be asked: What about<br />

the case when the external force is independent <strong>of</strong> the system’s<br />

coordinates In such a case, the right hand side <strong>of</strong> Eq.<br />

22 does vanish. However, with such kind <strong>of</strong> external<br />

forces, the system may not transition from state A to state B<br />

because <strong>of</strong> the stochastic nature <strong>of</strong> its dynamics. In the experiments<br />

<strong>of</strong> mechanically unfolding biological molecules at<br />

constant speed, the forces acting on the system are indeed<br />

functions <strong>of</strong> the coordinates <strong>of</strong> the terminus atoms attached<br />

to the tips <strong>of</strong> the atomic force microscope or the optical<br />

tweezers. CFT and JE are inapplicable to those experiments<br />

unless the pulling speeds are slow enough that the system is<br />

always near equilibrium.<br />

In summary, a new FDT has been derived from the<br />

Brownian dynamics without invoking any assumption. This<br />

new FDT is valid in the near-equilibrium or far nonequilibrium<br />

regime, wherever the Brownian dynamics is valid. It is<br />

also shown that the CFT and the JE are valid only in the<br />

near-equilibrium regime where the dissipative work is nearly<br />

zero.<br />

The author thanks Robert Buckley for a critical reading<br />

<strong>of</strong> the manuscript. He acknowledges support from a NIH<br />

grant Grant No. SC3 GM084834 and the Texas Advanced<br />

Computing Center.<br />

1 See, for example, J. Liphardt, B. Onoa, S. B. Smith, I. Tinoco, Jr., and C.<br />

Bustamante, Science 292, 733 2001; G. Hummer and A. Szabo, Proc.<br />

Natl. Acad. Sci. U.S.A. 98, 3658 2001; S. Park, F. Khalili-Araghi, E.<br />

Tajkhorshid, and K. Schulten, J. Chem. Phys. 119, 3559 2003; A.Imparato<br />

and L. Peliti, Phys. Rev. E 72, 046114 2005; D. Collin, F. Ritort,<br />

C. Jarzynski, S. B. Smith, I. Tinoco, and C. Bustamante, Nature London<br />

437, 231 2005.<br />

2 F. M. Ytreberg, R. H. Swendsen, and D. M. Zuckerman, J. Chem. Phys.<br />

125, 184114 2006.<br />

3 C. Jarzynski, Phys. Rev. Lett. 78, 2690 1997.<br />

4 G. E. Crooks, Phys. Rev. E 61, 2361 2000.<br />

5 E. G. D. Cohen and D. Mauzerall, Mol. Phys. 103, 2923 2005.<br />

6 I. Kosztin, B. Barz, and L. Janosi, J. Chem. Phys. 124, 064106 2006.<br />

7 J. M. G. Vilar and J. M. F. Rubi, Phys. Rev. Lett. 100, 020601 2008.<br />

8 L. Y. <strong>Chen</strong>, J. Chem. Phys. 129, 091101 2008.<br />

9 H. Risken, The Fokker-Planck Equation Springer-Verlag, New York,<br />

1989.<br />

10 J. C. Phillips, R. Braun, W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa,<br />

C. Chipot, R. D. Skeel, L. Kale, and K. Schulten, J. Comput. Chem. 26,<br />

1781 2005.<br />

Author complimentary copy. Redistribution subject to AIP license or copyright, see http://jcp.aip.org/jcp/copyright.jsp

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!