常温分子动力学中的多尺度建模与几何数值方法

来源 :中国科学院数学与系统科学研究院 | 被引量 : 0次 | 上传用户:liongliong533
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文作者的研究兴趣是将分子动力学中的应用建模与数值积分有机地结合到一个统一的理论框架之下。数学建模技术不仅要呈现自然之美,也应该有助于科学计算的有效性和可行性。另一方面,科学计算中所面临的困难与挑战也推动了相应建模技术的发展。二者间的交互作用已经在分子动力学多尺度建模和几何(保结构)数值积分的发展中得到了清楚的体现。本文着重研究常温(正则系综)分子动力学,并对其进一步的发展提出了一些基本思路。   在第一章,我们针对常温分子动力学的多尺度建模问题,发展了一个有几种变形的动力学框架。所阐述的方法可以适用于多种应用,包括轻重粒子的混合体和刚性势能模型。在遍历性假设的前提下,正则系综属性理论上得到了证明。我们也讨论了这种建模方法对相关数值方法发展的影响和建议,并以一些简单模型问题(包括一个非调和的双原子气体模型)的数值实验来验证了所研究的技术。   在第二章,我们讨论了固体物理中多尺度建模的声子反射问题。多空间和时间尺度的材料现象的研究往往需要对问题做粗粒度展示。多尺度在常温条件下带来问题一不同尺度交界处的声子反射会导致最小尺度(原子)区域的非物理的局部过热。我们建立了一个物理模型,其有效地控制了粗粒度区域的温度并保持了原子区域的正确动力学行为。一个多尺度的加热器系统被有效实施来达成目标。模型的合理性在一个粗粒度化的一维Lennard-Jones链中得到验证。   在第三章,我们回顾和扩展了当前在多尺度平均积分子方面的研究,其主要应用于物理多体问题,如分子动力学、材料建模和天体力学。很多方法,如平缓冲量和可逆平均方法,已经被提出来解决特定结构的多尺度问题的直接数值积分。常温粗粒度化分子动力学,作为我们研究兴趣之所在,其特点要求对当前体系框架做进一步的扩展。同时,平均值的计算在一些应用中至关紧要,但是现存的方法在这方面都具有缺陷。我们展示了一个新的方法,通过引入影子变量(作为真实物理变量的映象)来提高、简化平均值的精度和计算,从而增强多尺度方法的适用性。影子变量可以从一个辅助方程计算获得。尽管构造扩展相空间的几何数值积分子是可能的,但是在实践中我们发现方法的一个变形一使用耗散控制影子变量的偏移,获得了长期的能量保持。这导致了扩展相空间中的形式几何不变量一时间可逆和相体积的丧失,却稳定了物理变量的相关属性。数值方法应用于一个重力域的三体问题和一个部分温控的稀疏双原子气体模型。   在第四章,我们讨论了分子动力学正则系综下的高级温度控制技术。我们首先概述了动力学温控方法的现状,方法包括Nosé-Poincaré和推广的加热器技术(引入一个更复杂的扩展模型以获得更好的遍历性)。阐述了一个通常的温度控制模型一投影温控分子动力学(PTMD),它灵活地包含了现存的其它温控方法如Nosé-Poincaré,Nosé-Hoover(有链或无链),Bulgac-Kusnezov,递归Nosé-Poincaré等等。(PTMD)为热力学量的计算提供了可能的优势,有助于多尺度建模和仿真技术的发展。此外,PTMD通过避免感兴趣变量受到温控器的虚假影响,为实现选定自由度的真实非平衡态行为,提供了一个有意义的技术思路。   在第五章,我们主要讨论了以前章节中相关建模和数值技术的进一步研究发展。研究主题包含了Nosé动力学的扩展,稀有事件的增强采样,投影常温常压(Gibbs)分子动力学,随机耗散粒子动力学和Langevin动力学等等。
其他文献
在本文中,我们研究了一个输运方程的反问题。它与脉冲中子伽马辐射捕获测井的问题有关。在第一个问题中,我们将考虑来自脉冲中子源的快中子的分布,其初始能量为0.5MeV。我们还研
本文主要以带有稀疏结构的集合约束优化问题为研究对象,借助切锥和Fréchet法锥来刻画稀疏优化问题的四种必要条件:S-稳定点、L-稳定点、B-稳定点、BF点,并且分析了四者之间的
“删繁就简三秋树,领异标新二月花”让我们简简单单教语文吧!所谓的简单,不是方法单一,不是手段呆板;而是简而精,就是把繁杂的的教学方法变得简单、清晰、流畅,这些都需要教
博物馆藏品具有深厚的社会以及人文历史价值,这也是博物馆投入大量资源和精力进行藏品收藏的原因.为了有效发挥出博物馆藏品的内在价值,必须做好博物馆藏品的管理保护以及合
本文主要做了以下三方面的工作。 第一部分即论文第二章,我们考虑了在修改的Chebyshev节点阵下,|x|a的Lagrange插值多项式在零点的收敛速度可以小于等于c/n2α,其中c是仅与α
本文主要分为四个章节。首先本文简单地介绍了算子代数、算子理论以及Kadison-Singer代数的发展情况,给出一些基本概念,并对本文所要重点研究的问题作了描述。  然后本文构造
积分微分方程在生物科学、应用数学、物理学等诸多领域中有着广泛的应用,一直被许多国内外研究学者所关注.如何高效准确的求出积分微分方程数值解既是数值分析研究领域的重点
顶点代数理论是数学中一个非常活跃和重要的领域.类比特征零时的研究方法,本文对素特征域上的顶点代数的张量积的性质进行了研究,证明了素特征域上的顶点代数的张量积仍然具有
本文首先简单介绍了一些代数的基本概念以及投影算子的一些基本性质,然后介绍了保零积映射的概念以及基本性质,并且应用它们以得到我们所要的主要结果。设φ为从有单位元的代数
随着医疗水平和计算机运行速度的发展,利用基因序列中的信息完成对癌症患者进行亚型分类的研究成为生物统计的一大热点。在基因水平上对癌症患者进行亚型分类,可以为患者的医疗决策进行辅助,从而达到精准医疗,提高治愈概率,降低医疗成本,因此这无论对于诊断还是医疗都至关重要。然而,基因数据同时还有小样本,高维度,分布不均衡,大量噪音基因等特点,甚至单个相关基因位点对于癌症的相关性也比较微弱以至于许多模型无法准确
学位