Hamilton力学的辛算法.ppt

Hamilton力学的辛算法和分子动力学模拟 陈敏伯 中国科学院上海有机化学有机所 计算化学课题组 2006年10月 内容 冯康对世界科学的重大贡献 Euclid空间 辛空间 Hamilton力学的辛结构 正则变换的辛结构 辛算法应用实例 Schr?dinger：“Hamilton原理已经成为现代物理学的基石。” Hamilton原理将不同的物理规律纳入了统一的数学形式。 现在问题就归结到：怎样才能对Hamilton力学的运动方程作正确的数值计算。 一切Hamilton体系的动力学演化都使辛度量保持不变，即都是辛（正则）变换。 一切解Hamilton方程 “正确”的离散算法都应当是辛变换的。 （冯康，1997年国家自然科学一等奖“哈密尔顿系统辛几何算法” ） Lax：“他的声望是国际性的。” 丘成桐：“中国…在数学历史上很出名的有三个：一个是陈省身教授在示性类方面的工作，一个是华罗庚在多复变函数方面的工作，一个是冯康在有限元计算方面的工作。” (1998年3月11日《中国科学报》) “冯氏大定理” 同一物理定律的不同的数学表述，尽管在物理上是等价的；但在计算上是不等价的。 冯康：如果在算法中能够保持辛几何的对称性，将可避免人为耗散性这类算法的缺陷，成为具有高保真性的算法。 在天体力学的轨道计算，粒子加速器中的轨道计算和分子动力学计算中得到广泛的应用。 冯康（1920－1993）的学术成就 1965年发表论文“基于变分原理的差分格式”。国际学术界承认冯康独立发展了有限元方法。 （仅获1982年国家自然科学二等奖。冯康得悉非常难过，曾打算将申请撤回。） 前国际数学会理事长J. –L. Lions教授1981年说：“中国学者在对外隔绝的环境下独立创造了有限元，在世界上是最早之列。今天这一贡献已为全人类所共享。” 1984年以后创建的“哈密尔顿系统的辛几何算法”。 （1991年评为国家自然科学奖二等奖。冯康获悉后撤回申请。直到1997年底，在冯康去世四年之后，终于授予了国家自然科学一等奖。 ） 石钟慈：“国际上最早系统地研究并建立辛几何算法的。” 数学地位 外微分 辛几何 辛几何的基础是外微分形式。 外微分形式是如下概念推广到高维的产物： 1、作功—在场中沿某一路径所作的功； 2、流量—单位时间内流体穿过某曲面的量 3、面积或体积—平行四边形面积 或平行六面体体积。 外微分形式中有“1-形式”、“2-形式”等 辛构造就是非简并的闭2-形式。 Euclid空间 对称性： 线性： （ k 为任意实数） （ c是V中的任意向量） 非简并性： ，当且仅当 时才 辛空间（Simplectic Space） 反对称性： 双线性： 非简并性：若向量a对于W中的任意向量b均有 ，则 辛空间 度量：作功、面积（或体积）、流量等 辛内积： 2维： a、b平行四边形面积 2n维： 单位辛矩阵 的性质 若 A 为对称阵，且 ，则 Euclid空间和辛空间的对应关系 Hamilton力学的辛结构 正则变换的辛结构 正则变换M的性质 无穷小辛阵 定义：若 ， 则该2n阶矩阵 称为“无穷小辛阵” 设 为对称阵，当且仅当 时 ， 为无穷小辛阵。 （证明略） 若 为无穷小辛阵， 则 为辛阵 。 若 为无穷小辛阵，又若 非奇异， 则 为辛阵 辛阵 线性Hamilton体系的辛差分格式 线性Hamilton体系——Hamilton函数是 的二次型 中点Euler法的辛格式 可分、线性Hamilton体系的中点Euler公式 —— “可分、线性Hamilton体系” 演绎细节： 基于Padé逼近的辛格式 线性Hamilton体系相流 有理Padé逼近： 用以下 构造的差分格式都是辛格式 可分线性Hamilton体系的交叉显式辛格式