论文部分内容阅读
Monte Carlo方法作为认知自然规律的重要科学研究方法之一,广泛应用于高分子科学中的模拟计算,不仅可以直接模拟高分子科学中存在的大量随机性问题,而且可以模拟高分子的各种微观状态。通过模拟高分子链的静态和动态性质,在分子水平上模拟高分子构象演化过程。Monte Carlo方法可以自如地控制各种物理基本规律的近似模拟。通过Monte Carlo模拟可以计算高分子体系的有关物理量,还可能发现某些新的物理规律。因此,理论和实验方法所尚未解决的重大科学问题往往是Monte Carlo方法研究的热点。高分子的表面吸附性能对大量的工业应用起到关键作用,许多生物工程也与其相关。一方面,表面对其附近的高分子构象发生了改变,另一方面,表面性质也会受吸附层的高分子构象和形状影响,然而,高分子在表面吸附性能的很多机理还不清楚。目前,实验研究对于构象信息的采集非常有限,不能全面反映高分子在表面的吸附构象信息;理论计算往往基于一定的经验公式,所计算得到的构象信息更为有限。计算机模拟因成本低和操作简单,可以得到更全面的构象信息而受到研究者关注。人们总是希望采用有效并快速的算法,模拟高分子链的动力学过程,计算其吸附和构象性质并对其进行标度,从而发现新的科学规律。采用不同的方法和模型对高分子科学进行计算机模拟,发现高分子链随着温度(或者作用强度)的变化会在表面发生临界吸附相变。计算得到的临界吸附点和交叉指数存在较大差别,还没有一种统一的标度计算方法。然而,采用不同的表面结构,不同单元性质组成的高分子链,将会产生更加复杂的效应,临界吸附标度和构象性质统计的新算法非常值得去研究。本文在国家自然科学青年基金“高分子链的表面临界吸附现象的动力学模拟研究(11304231)”,国家自然科学基金“高分子修饰的无机纳米棒自组装的模拟研究”(21171145)”,浙江省自然科学基金“共聚高分子链在表面的临界吸附性质研究(LY14B040004)”等科研项目的支持下,基于自回避行走方法对高分子链在表面临界现象的建模和标度算法中的若干关键模拟计算技术进行研究。本论文的主要研究内容和创新点如下:第一,提出一个基于Monte Carlo方法模拟高分子链在表面临界现象的模型。采用统计力学方法分析了稀溶液高分子链体系的模拟理论,采用自回避行走方法模拟高分子链的生长模型,采用键长涨落的方法模拟自回避链的微松弛运动,实现高分子链在表面临界行为的数值模拟计算。通过抽样统计体系能量、均方回转半径、均方末端距和形状因子等物理量,解决了高分子链在表面上吸附和构象性质等科学问题的模拟计算。第二,提出均聚链/共聚链有限尺寸临界吸附的标度算法。采用动态Monte Carlo方法模拟一端接枝在表面的均质高分子链AN的吸附性质,提出基于均质自回避行走链的有限尺寸临界吸附标度算法。通过计算机模拟发现高分子链接触表面的链单元数在临界吸附点附近以及在更大温度范围内仍然满足统一标度。进而,针对增加一段弱吸附作用(或无吸附作用)的高分子链段B构建共聚高分子链模型,提出一个两嵌段链AnNm的临界吸附标度算法,采用外推和逼近方法分别从吸附链单元数、构象和形状因子等方面分析和比较,解释了两嵌段共聚链AnBm与均质链An临界吸附点的一致性。采用有限尺寸标度算法解决临界温度和交叉指数的计算问题,为实验中链长未知的情况下计算临界吸附点提供便捷算法。第三,在上述研究模型中增加了平行表面以及条状纹理图案表面,提出高分子链相图的标度算法。在两个平行表面的模型上,提出了两平行表面的临界间距标度算法,解释了平行表面间高分子链的两种不同情况——强约束态和弱约束态。在周期性纹理表面的模型上,提出了高分子链在空间区域分布规律的计算方法,并采用外推法构建了无限长链在有限条纹宽度表面上的相图,分别为脱附相、多条纹区域吸附相和单条纹区域吸附相。在一端接枝的两嵌段共聚链AnBn在表面上接枝位置不同,提出了两类不同的模拟模型XAnBn和XBnAn,该模型解释了两类共聚链构象之间存在的微妙差异。最后,为高分子链的Monte Carlo模拟提出并行计算优化。基于样本独立性的特点,提出高分子链的并行Monte Carlo模拟,解决了各样本之间可以同时计算的问题。通过在集群上并行计算发现可以达到理想的线性加速比。通过在并行计算中估计非平衡态的驰豫时间和平衡态的驰豫时间,从而可以估算在每个温度的模拟系统达到平衡态和抽样统计所需的计算时间,提高计算的精度和计算资源的利用率。在本文研究工作的基础上,下一步还需要在非格子模型、Wang-Landau算法和进一步并行优化等方面开展工作。