论文部分内容阅读
分子动力学(molecular dynamics,MD)模拟已成为研究微米或纳米尺度动力学、流体及固体统计特性的重要工具。目前,分子动力学主要分为两大类:平衡分子动力学(equilibrium molecular dynamics,EMD)和非平衡分子动力学(nonequilibrium molecular dynamics,NEMD)。本文针对MD模拟提出一种新的基于能量转换的速度修正法(energy conversion based correction,ECBC)。不同于传统的依赖于宏观参数的系统性速度修正方法,ECBC方法基于分子或原子对间的相互作用物理过程中的能量守恒与转换特性发展而来。采用ECBC修正方法进行MD模拟时,由于不用施加额外势能或引入温度等宏观参数,因此毋需针对NEMD采用特殊的处理方法即可直接用于EMD和NEMD模拟,特别是非保守(开放)系统的MD模拟。该方法突破了现有MD计算方法基于保守平衡系统思想引起的许多应用限制,可以极大的拓展MD应用范围。 微纳颗粒在流体中的布朗运动主要取决于流体分子对颗粒的作用力。由于颗粒尺寸极小,并且流体分子与其碰撞频率极高,因此难以通过实验手段测得颗粒所受布朗力。本文采用MD方法研究作用于颗粒的布朗力,发现一些很有意义的新现象:(1)对于不同的1/Kn数,布朗力的概率分布都符合高斯分布且具有自相似性;(2)与常用于郎之万模型不同,布朗力的频谱图并非白噪声,而是具有优先频率的类似对数正态分布的噪声;(3)随着微纳颗粒的尺寸变化,作用于颗粒的布朗力呈现出非单调性。当1/Kn数由小变大时,作用于颗粒的布朗力首先增加并且在1/Kn≈250达到最大值,之后随着1/Kn数继续增大,布朗力逐渐减小。本文发现的布朗力的变化规律对于颗粒的布朗运动的动力学特性研究将起到很好的促进作用。