基于Monte Carlo方法模拟自回避链的计算方法及并行优化

来源 :山东大学 | 被引量 : 0次 | 上传用户:or4108432566
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
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算法和进一步并行优化等方面开展工作。
其他文献
文章运用第五次、第六次全国人口普查上海市居村委数据,建立上海人口地理信息库,发现上海城市扩张主要的空间范围及前沿地带为外环线附近工业集聚带以内地区,中心城空心化与
以TiOSO4·2H2O为原料,采用沉淀法制备了TiO2催化剂,考察焙烧温度对TiO2催化剂催化燃烧二氯甲烧(DCM)性能的影响.采用X射线衍射(XRD)、BET比表面积测试、扫描电子显微镜(SEM)、
近年来,我国电网工程任务饱满,施工建设如火如荼,电力物资的采购供应工作越来越受到重视,而做好采购工作的首要任务就是签订买卖合同。虽然在生活中存在很多现成的买卖合同样
作为现代企业三大核心财务活动之一,股利分配决定着企业的利润分配水平、分配方式、分配时间等,它与剩余价值的分配有关,对于企业而言有着不可或缺的重要意义。本文对上市公
现代产权理论的主要内容现代产权理论产生于 20 世纪 30年代末期,60 年代开始在美国兴起,随后得到了大力发展.由于该理论弥补了传统经济学的某些缺陷和不足,重新确立了经济理
期刊
当前,多核处理器已经广泛应用于从巨型机、桌面处理器到移动终端等各个计算领域,计算技术已经全面迈入多核和众核时代。本文回顾了高性能微处理器的发展趋势:工艺的影响、体
引言SDR(special drawing right)特别提款权,是国际货币基金组织于1969年设立的一种储蓄资产与记账单位.人民币加入SDR是中国经济进一步融入全球金融体系的重要里程碑,这不仅
期刊
景观生态与高等教育在系统上具有耦合性,景观生态学的理论与方法可以为高等教育学的研究提供全新视角与技术支持,高等教育系统的整体性和异质性研究应是高等教育研究的重点,
文章研究了一种基于最小二乘支持向量机(LSSVM)和高斯混合模型(GMM)的传动系统故障程度量化方法。在传动系统故障机理的基础上,利用风电机组的风速、环境温度、有功功率和4个