基于随机漫步模型模拟二维及三维含边界布朗运动-东南大学.PDF

基于随机漫步模型模拟二维及三维含边界 布朗运动 余浩斌 （东南大学物理系，南京 211189） 摘要：布朗运动是分子物理的重要内容, 它对人们认识物质的结构, 热 运动的本质及其规律都有重要意义。本文利用随机漫步模型,通过蒙特卡 罗法在计算机中模拟二维及三维含边界布朗运动，并与二维及三维自由 边界布朗运动对应比较，探究布朗粒子位移平方的平均值与时间的关系。 关键词：随机漫步模型，布朗运动，蒙特卡罗法 1 引 言 1827年英国植物学家布朗在显微镜下观察到悬浮在水中的花粉都在不停的作无规则运 动, 最初以为是外界影响所致, 后来精确实验表明, 在排除外界干扰的情况下, 只要悬浮在 液体或气体中的微粒足够小, 总能观察到这种无规则运动, 这称之为布朗运动。布朗的发现 是一个新奇的现象，它的原因是什么,人们是迷惑不解的。 在布朗之后，这一问题一再被提出，为此有许多学者进行过长期的研究。一些早期的 研究者简单地把它归结为热或电等外界因素引起的。最早隐约指向合理解释的是维纳 (1826——1896)，1863 年他提出布朗运动起源于分子的振动，他还公布了首次对微粒速度 与粒度关系的观察结果。不过他的分子模型还不是现代的模型，他看到的实际上是微粒的位 移，并不是振动。1905 年，爱因斯坦依据分子运动论的原理提出了布朗运动的理论。就在 差不多同时，斯莫卢霍夫斯基也作出了同样的成果。他们的理论圆满地回答了布朗运动的本 质问题。1908 到1913 年期间，贝兰通过测量布朗运动的实验验证了爱因斯坦理论的正确性。 就这样，布朗运动自发现之后，经过多半个世纪的研究，人们逐渐接近对它的正确认识。布 1 朗运动也成为分子运动论和统计力学发展的基础。布朗运动代表了一种随机涨落现象，它的 理论在其他领域也有重要应用，如对测量仪器的精度限度的研究、高倍放大电讯电路中的背 景噪声的研究、金融股票的研究等。 在实际应用中，布朗运动是有边界的。本文利用随机漫步模型，通过蒙特卡洛法，模 拟含矩形边界的二维及三维布朗运动。在分析数据时，发现含边界布朗粒子位移平方的平均 值与时间的关系与自由边界布朗粒子存在差异。 2 基本原理及方法 2.1 随机漫步模型 关于随机漫步，可以分为三类，第一类是简单的随机行走(random walk)，假如有一个 灰尘粒子，在空气中做随机行走，不考虑重力的影响，灰尘向各个方向运动的几率都是一样 的，如在二维空间中，粒子有四种选择，即上下左右，每个方向的概率都是0.25。当粒子变 得更聪明了，即有了简单的记忆之后，就可以变成更加复杂一点的随机行走，例如人走路， 每一步会时而左，时而右，总体方向朝前，几乎不会原地连身体都不转一下就朝后行走（除 非有目的的朝后走），又例如，一款经典的小游戏—贪吃蛇中，蛇可以向前，向上，向下行 走，但是，它不能直接向后，即此时随机行走有了一定的记忆性，它记住前一步的位置并避 免朝回走(non-reversal random walk)。那么进一步，当粒子不仅仅能够记住之前的一步的位 置，而是此时，粒子的记忆力增强了，它记住了自己走过的所有地方，由于一种特殊的爱好， 这个粒子只喜欢新鲜的没有走过的位置，即每次粒子都不能够再走那些它走过的地方。此时 称之为自回避随机行走（self-avoiding random walk ）。 我们采用第一类随机漫步模型, 用计算机来模拟二维布朗粒子的运动。对于布朗粒子所 受到的各种力, 根据分子的统计规律, 我们首先采取一个简单的模型, 即假设每次碰撞后 粒子移动的距离相同, 布朗粒子的运动方向则是随机的, 由计算机产生的均匀分布的随机 2 数来决定, 由此来确定布朗粒子每走一步所在的位置。如图 1所示。 设某个时刻 t, 粒子位于O( x , y ) 位置, 当粒子受到随机力后, 其位移为 r k , 我 i i 们取k