带电粒子形成胶体晶体的布朗动力学模拟.pdfVIP

带电粒子形成胶体晶体的布朗动力学模拟 2) 顾 凌云，欧阳 文泽 ,，徐 升华， 孙祉伟 （中国科学院力学研究所微重力重点实验室，北京 100190) 摘要：本论文由布朗动力学模拟方法，观察了带电的胶体粒子有序化形成晶体的过程。根据DLVO 理论，带电胶 体粒子间的相互作用势采用Yukawa 排斥势模型。模拟中，我们通过计算体系的RDF 曲线和均方位移来分析体系 结晶过程和晶体结构随时间的演变。当带电粒子间具有长程排斥相互作用时，体系在较低的粒子数密度下，就能 够自发形成有序的体心立方 （BCC ）结构。当带电粒子间具有短程相互作用时，虽然相图中的稳定结构为面心立 方 （FCC ）结构，但是在直接的动力学模拟中很难实现，体系会陷于无序的类玻璃态结构，从而阻碍了其向有序 化结构的转变。进一步的模拟结果显示，当玻璃态结构被越过后，体系可以形成FCC 晶体结构，并且比该条件下 的BCC 结构更加稳定。 关键词：胶体晶体；带电粒子；Yukawa 作用势；布朗动力学模拟 引 言 复杂流体(ComplexFluids)，即物理学中的软凝聚态物质，是介于理想固体和流体之间具有复杂本构 【1】 关系的物质，在自然界和工业生产中广泛存在，包括表面活性剂、液晶、高分子溶液等 。胶体悬浮粒子 体系作为一种相对简单的复杂流体体系，能够涵盖很多典型性质，是一个很好的研究复杂流体的模型体系。 胶体晶体(colloidalcrystals)是胶体粒子在溶液中自组织形成的三维或者二维有序结构。近年来，胶体 【2】 晶体的相关研究已成为热门课题，它不但可以作为原子、分子晶体研究的模型体系 ，而且具有重要的应 【3】 【4】 用价值，例如可以作为模板，用来制备多孔材料 ；制备具有三维结构的光子晶体 ，以用于纳米技术、 计算机技术、芯片等领域。 在带电粒子形成胶体晶体的研究方面，由于带电粒子之间相互作用力的强度和作用范围都可以变化， 其相变行为比硬球体系更加复杂。在实验中，相互作用力的强度和作用范围并不是独立的参量，而是受到 胶体粒子尺寸大小、表面性质和表面电荷密度、溶剂离子浓度等多种因素的影响，存在一定的函数关系， 【5】 其具体数值需要通过实验测定来确定 。在计算机模拟中，可以通过改变相互作用势能的参数，方便地改 【6-7】 变粒子间相互作用力，从而进行比较研究。文献中已经采用蒙特卡洛(Monte Carlo, MC) 和分子动力 【8-9】 学等方法(MolecularDynamicsmethod,MD) 来研究带电胶体粒子形成结晶，并取得了一定的进展。在 胶体粒子体系中，由于溶剂分子远远小于胶体粒子，并且溶剂分子运动的特征时间尺度比胶体粒子也小几 个数量级，因而要在溶剂分子的时间尺度上计算体系的物理参量和观察体系的动态演化过程，从所需的计 1）国家自然科学基金和中国科学院创新工程 （KJCX2-YW-L08 ）资助项目 2 ）E-mail ：oywz@imech.ac.cn 算时间和计算量角度看不但很难实行，而且程序会非常复杂。布朗动力学方法作为一种介观的模拟方法， 关注胶体粒子的运动行为，忽略溶剂分子的具体运动信息，更好的解决了胶体粒子和溶剂分子的时间尺度 差异问题，从而提高了计算效率。这一处理方法更加恰当地模拟了胶体粒子的行为，给出了一种简单而有 效的胶体粒子运动模型。 本论文主要采用布朗动力学模拟方法，研究了带电粒子形成胶体晶体的过程以及其中的相关问题。带 电粒子有