论文部分内容阅读
在利用分子动力学模拟研究聚合物的相分离、多级自组装、以及其它一些性质的领域，一个最突出的问题是模拟的时间和空间尺度都很难达到一些现象发生的实际范围。例如具有缠结的高分子链的松弛和一些蛋白质的折叠等现象经常发生在微秒的时间尺度和微米的空间尺度下。由于全原子分子动力学模拟考虑了太多的分子细节，目前很难对一些发生在介观尺度的现象进行模拟。因此增大分子模拟的空间和时间尺度是很有必要的，对此我们从如下两方面做出了一些努力：一方面通过发展高效率的分子动力学程序包GALAMOST，充分利用GPU（图形处理器，Graphical Processing Unit）的高性能计算能力，来提高分子模拟的运算能力；另一方面同时发展粗粒化模拟方法，通过体现分子的主要特征来研究我们所关心的问题，以达到增大分子模拟的空间和时间尺度的目的。本文的主要内容包括：(1)首先介绍我们自主开发的高效率的程序包GALAMOST。它是一个充分利用GPU运算能力的分子动力学模拟程序包。除了一般的分子动力学模拟的功能外，它通过利用一些最新的模拟技术，可以用来研究发生在介观尺度的自组装、相转变、以及一些其它关于高分子体系的性质。为了加速模拟过程和程序运行速度，一种杂化的粒子密度场分子动力学方法包含在程序包内。这种方法结合了自洽场和分子动力学的思想，用粒子和密度场的相互作用取代了粒子和粒子之间的成对相互作用，而这种成对相互作用是分子动力学中计算最耗时的部分。因此这种方法可以大大加速体系的相形貌的演化过程。除了一些解析势能形式外，一种通过自下而上的粗粒化方法得到的数值势能也可以借助GALAMOST应用在模拟之中。数值势能可以通过迭代波尔兹曼或一些其它基于结构的粗粒化方法从全原子模拟或者实验数据中得到。利用得到的数值势能我们可以在同样的热力学条件下模拟更大尺度的体系。此外，GALAMOST还包含两种特定的模型，软的各向异性粒子模型和链增长聚合反应模型。应用这两种模型可以分别研究软的各向异性粒子的多级自组装和一些与聚合反应相关的问题。由于不断优化在GPU上执行的算法，GALAMOST包含的每种方法的运算效率都得到了很大的提升。(2)除了在程序方面的努力外，我们还发展了一种超高粗粒化模型来研究高分子熔体。在这个模型中一个粒子代表一条高分子链，缠结的产生和消失被看作是在单位时间间隔内的几率事件。我们建立了单位时间间隔内任意两条邻近分子链之间发生缠结的几率和缠结变化速率之间的联系。通过调节缠结消失的几率来保持总的缠结数量在目标数量附近。这个模型不仅可以用来反映缠结的一些基本特征以及高分子熔体的一些宏观动力学性质，还可以研究微观分子特征和宏观性质之间的关系。作为一个示例，我们应用这个模型模拟了聚乙烯C1000H2002在450K下的动力学性质，并且把模拟得到的结果和其它实验数据和模拟数据做了比较。(3)对于一些体积更大的胶束、树枝状分子以及胶体粒子，我们通过简谐势能来描述它们之间的相互作用。我们发现这种简单的势能形式可以用来描述球形可互相穿透的粒子的热力学性质和动力学行为。我们结合耗散粒子动力学和热力学积分方法，通过计算其液相的化学势以及固相的各种晶体结构的化学势，来绘制简谐势能粒子的相图。在相图中，我们除了发现了经常出现的面心立方和体心立方等晶格结构外，还发现了四方、六方、斜方和金刚石等稳定结构。在相图中我们也找出了多个凹形熔化区域和多形态的晶体相转变。