分子动力学模拟第五章统计力学及分析方法.pdf

Molecular Dynamics Modelling and Simulations

r v 1 A ens = N òò A (p , r )p (p , r )dpdr

, ens(ensemble), p(p,r) ( ) A 1 2 2 2 2 p (A) = exp[-(A - A ) / 2s ] s = (A - A ) s 2p NVE NVT é - H (p , r ) ù p (p , r ) =d (H (p , r ) - E ) p (p , r ) = expê ú k T ë B û é - H (p , r ) + PV ù é - H (p , r ) +µ N ù p NPT (p , r ) = expê ú pµ VT (p , r ) = expê ú k T k T ë B û ë B û

(ergodic hypothesis) A = A ens time 1 1 A time = lim t ò A(p ,r )dt = t åA(p ,r ) obs obs v µ 2 2 (µ ) E = N p i = (3N - N c )kB T E = 1 n N (p i ) K å 2m 2 K n - n å å 2m i =1 i 0 µ >n0 i =1 i v 2 N E T = 1 å p i T = d K (3N - N c )kB i =1 mi Nk B 2 ü d d 3 v 1 n U = å åV(r (µ ) ) ij n - n 0 µ >n0 i <j v r c v g (r) ¥ 2 U / N = 2pr ò V(r )g (r )r dr 0 v ! W = åx p i xi v -3NkB T N N dU( r ) W = -3 PV + å å r ij