计算物理方法.pptVIP

对于相空间中点集的一切互补的对偶集，存在着和，使得。这是相空间区域的连通性和遍历性的陈述。由于概率正定性的要求，对于一切，。概率归一化要求，对所有的，。由于要求极限分布为平衡分布，所以对所有的x，可以证明，只要游动所选的”过渡几率”满足如下的细致平衡条件．就可以达到平衡时分布为这样的目的：第30页，共63页，星期日，2025年，2月5日第五项中的细致平衡条件实际只是一个充分条件，并不是一个必要条件。该条件并不能唯一确定过渡几率。所以，过渡几率的选择有很大的自由度。不同的选取即是不同的方法；在Metropolis方法中一般采用一个简单的选择过渡几率的方法：第31页，共63页，星期日，2025年，2月5日设游动到了点，以下的具体操作步骤为：（1）首先选取一个试探位置，假定该点位置为：，其中为在间隔内均匀分布的随机数。（2）计算的数值。（3）如果不等式满足，那么进行这一游动，并取，返回（1）开始对游动到点的试探。（4）如果，产生一个区间内的随机数。第32页，共63页，星期日，2025年，2月5日（5）如果此时，那么还接受这步游动，并取这步游动所到达的点为，然后返回到步骤（1），开始下一步到达点的游动。（6）如果此时，就拒绝游动，仍留在的位置不变。（7）返回到步骤（1），重新开始对游动到点的又一次试探。第33页，共63页，星期日，2025年，2月5日必须指出：采用这样的游动过程，只有在产生了大量的点后．才能得到收敛到满足分布的集。这里有一个明显的重要问题，就是如何选择的大小，才能提高游动的效率?如果选得太大，那么绝大部分试探的步子都将会被舍弃，就很难达到平衡分布：反之，如果取得太小，那么绝大部分试探步子都会被接受，这同样难以达到所要求的平衡分布。根据实际应用的经验，选取的一个粗略标准应当是：选择适当大小的原则是要在游动的试探过程中，有1／3到1／2的试探步子将被接受。按照这样的标准选择得到的，就可以大大提高游动的效率。第34页，共63页，星期日，2025年，2月5日另一个在Metropolis方法中的问题是：进行这样的随机游动，从哪一点出发才可以比较快地达到平衡分布呢?原则上讲，从任何一个初始位置出发均可达到平衡分布，但是为了尽快地达到平衡分布，我们最好是要选择一个合适的初始位置．这个初始位置应当是在游动范围内所要求的几率分布密度最大的区域。第35页，共63页，星期日，2025年，2月5日四、例子---一维简谐振子假定有一个质量为m的粒子，在一维简单简谐势中运动。取为长度单位，为时间中的单位。有第36页，共63页，星期日，2025年，2月5日首先，选择任意的、连接个时间间隔、且的一条路径，计算上式的能量，然后再选一系列路径，每条路径与前条路径最多只有在一个时刻（例如），有不相同的空间点（见图）。采用Metropolis方法来确定满足上面要求的新径迹。其中将随机定下的坐标改变到的过渡几率为，为分别包括在时刻坐标为和的两条径迹的能量差，可由上式算出。这样的随机游动抽样得到的径迹也许会与前一个径迹相同。每当新径迹选出后，就利用前式计算被积函数