约束动力学中的修正RATTLE方法及其在正则模态分析中的应用

来源 :中国科学院数学与系统科学研究院 | 被引量 : 0次 | 上传用户:momoww
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在分子动力学模拟中,时间步长受限于被模拟分子中的键长伸缩和键角张合这类高频运动的周期。这使得分子动力学模拟的时间步长非常小,通常为1飞秒。约束动力学通过约束键长或键角来达到延长时间步长的目的。由于当约束所有的键角时,约束方程的个数急剧增加,现有的约束动力学算法不适用于约束键角的动力学模拟,并且其中被最广泛应用的约束动力学算法SHAKE和RATTLE对求解键角约束存在收敛性问题。   在论文中,提出一种新的约束动力学算法:修正RATTLE方法(M-RATTLE)。在修正RATTLE方法中,引入一种新的校正所有原子位置坐标来满足约束条件的方法。这种新的校正方法理论上等价于在RATTLE方法中计算拉格朗日乘子,从而它不会引起总能量的耗散。而且这种校正方法特别适用于一般分子的刚体动力学和只有二面角未被约束的树状结构分子。在这些情形下,未被约束的自由度的运动空间的一组基能直接被写出,从而修正RATTLE方法只要求在每一步求解k维的线性方程组,其中k是未被约束的自由度的个数。这比直接求解m维约束方程的传统方法的计算量要少得多,因为k《m。本文将修正RATTLE方法应用于吲哚分子和苯分子的刚体动力学模拟及键长和键角都被约束的丙氨酸多肽。计算效果与传统方法的比较表明修正RATTLE方法具有更高的效率和更好的稳定性。并且将修正RATTLE方法应用于生物大分子的正则模态分析,提出了修正RATTLE子空间模态分析方法。对丙氨酸多肽和蛋白质:BPTI,1HHP,1VII的计算结果表明修正RATTLE方法子空间正则模态分析方法能够正确地得到生物大分子的低频模态。   论文的最后一部分是对蛋白质alpha-t-alpha(1ABZ)的粗粒化动力学模拟的工作。
其他文献
波分复用技术在提升全光网络的利用效率、充分发挥光纤巨大带宽等方面起到了非常重要的作用。波分复用技术在一条光纤上利用不同的波长建立多个通道、传输若干路信号,对于每个
学位
本文主要研究了针对双曲守恒律方程的大时间步长格式,证明了一维大时间步长格式的一些性质,并且说明了将大时间步长格式推广到高维问题的几种途径,最后给出了应用大时间步长格式
本文内容涉及Hamilton系统辛几何算法的三个方面:线性多步方法步推映射的辛性、Hamilton系统辛算法形式能量的有效计算、时域Maxwell方程的辛方法。主要成果如下:   1.基于
本文首先介绍了模糊中位数的定义,然后介绍了将模糊中位数应用于图像平滑滤波的方法,也就是模糊中值滤波。本文将模糊中值滤波的模拟结果与中值滤波进行了系统比较,并对结果做了
R.Lashof&S.Smale在1958年将超曲面的Gauss-Bonnet定理推广到一般的欧氏空间的子流形中,本文将采用同调论和示性类的方法,对该结果给出一个简单证明和一些应用。        
本论文主要研究低阶非协调有限元在一般四边形网格上的精度.   网格条件在工程计算中起着重要的作用,本文分析了一类非常实用而且在理论上也很有意义的四边形网格条件即(1+
本文将考虑下面的非线性椭圆方程。   在第一章,介绍了上述方程的背景并给出了主要定理。   在第二章,考虑上述方程正解的紧致性定理,首先采用F.Pacard的思想,建立在H1(Ω)
本论文研究了冯·诺依曼代数的生成元问题,首先给出了一些经典结论。生成元问题指的是可分希尔伯特空间H上的任何冯·诺依曼代数M是否由单个元生成,即是否存在A∈M,使得M={A,A*}"
2013-05-09国际浆纸网报道:在美国高盛银行4月底的报告中显示,近期大宗商品前景有所下调,高盛认为从中国到欧美市场需求前景表现疲软,弱于预期的宏观经济数据,增强其对全球经
本文主要讨论MMP(Mathematics Mechanization Platform)的系统结构及其高层编程语言实现与应用。MMP是由国家973项目资助的大型数学机械化平台软件,其核心功能是符号计算及其