Wednesday, July 31, 2013

phymath01 mktpmath01 第一性原理 分子动力学 在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波尔兹曼分布随机生成的,由于速度的分布符合波尔兹曼统计

第一性原理目录47


合外力 on a round top=0, 合nei力 not=0, all pointing to the core!

http://zh.wikipedia.org/wiki/%E5%88%86%E5%AD%90%E5%8A%A8%E5%8A%9B%E5%AD%A6
123


第一性原理目录;简介;其他解释:补充解释展开简介;其他解释:补充解释;展开编辑本段简介根据原子核和电子互相作用的原理及;第一性原理就是从头计算,不需要任何参数,只需要一;编辑本段分子动力学简史;1957年:基于刚球势的分子动力学法(Alder;分子动力学;1971年:刚体系への拡张(RahmanandS;1980年:恒压条件下的动力学方法(Anders;19
第一性原理目录
简介
其他解释:补充解释 展开简介
其他解释: 补充解释
展开编辑本段简介 根据原子核和电子互相作用的原理及其基本运动规律,运用量子力学原理,从具体要求出发,经过一些近似处理后直接求解薛定谔方程的算法,习惯上称为第一性原理[1]。 第一性原理通常是跟计算联系在一起的,是指在进行计算的时候除了告诉程序你所使用的原子和他们的位置外,没有其他的实验的,经验的或者半经验的参量,且具有很好的移植性。作为评价事物的依据,第一性原理和经验参数是两个极端。第一性原理是某些硬性规定或推演得出的结论,而经验参数则是通过大量实例得出的规律性的数据,这些数据可以来自第一性原理(称为理论统计数据),也可以来自实验(称为实验统计数据)。 但是就某个特定的问题,第一性原理和经验参数没有明显的界限,必须特别界定。如果某些原理或数据来源于第一性原理,但推演过程中加入了一些假设(这些假设当然是很有说服力的),那么这些原理或数据就称为“半经验的”。编辑本段其他解释: 第一性原理,英文First Principle,是一个计算物理或计算化学专业名词,广义的第一性原理计算指的是一切基于量子力学原理的计算。 我们知道物质由分子组成,分子由原子组成,原子由原子核和电子组成。量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量(或离子),然后就能计算物质的各种性质。 从头算(ab initio)是狭义的第一性原理计算,它是指不使用经验参数,只用电子质量,光速,质子中子质量等少数实验数据去做量子计算。但是这个计算很慢,所以就加入一些经验参数,可以大大加快计算速度,当然也会不可避免的牺牲计算结果精度。 那为什么使用“第一性原理”这个字眼呢?据说这是来源于“第一推动力”这个宗教词汇。第一推动力是牛顿创立的,因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。如果宇宙诞生之初万事万物应该是静止的,后来却都在运动,是怎么动起来的呢?牛顿相信这是由于上帝推了一把,并且牛顿晚年致力于神学研究。现代科学认为宇宙起源于大爆炸,那么大爆炸也是有原因的吧。所有这些说不清的东西,都归结为宇宙“第一推动力”问题。 科学不相信上帝,我们不清楚“第一推动力”问题只是因为我们科学知识不完善。第一推动一定由某种原理决定。这个可以成为“第一原理”。爱因斯坦晚年致力于“大统一场理论”研究,也是希望找到统概一切物理定律的“第一原理”,可惜,这是当时科学水平所不能及的。现在也远没有答案。 但是为什么称量子力学计算为第一性原理计算?大概是因为这种计算能够从根本上计算出来分子结构和物质的性质,这样的理论很接近于反映宇宙本质的原理,就称为第一性原理了。 广义的第一原理包括两大类,以Hartree-Fork自洽场计算为基础的ab initio从头算,和密度泛函理论(DFT)计算。也有人主张,ab initio专指从头算,而第一性原理和所谓量子化学计算特指密度泛函理论计算。 补充解释
第一性原理就是从头计算,不需要任何参数,只需要一些基本的物理常量,就可以得到体系基态的基本性质的原理。
编辑本段分子动力学简史
1957年:基于刚球势的分子动力学法(Alder and Wainwright) 1964年:质点系への拡张(Rahman)

分子动力学

1971年:刚体系への拡张(Rahman and Stillinger) 1977年:约束动力学方法(Rychaert等)
1980年:恒压条件下的动力学方法(Andersenの方法、Parrinello-Rahman法)
1983年:非平衡态动力学方法(Gillan and Dixon)
1984年:恒温条件下的动力学方法(能势‐フーバーの方法) 1985年:第一原理分子动力学法(→カー?パリネロ法) 1991年:巨正则系综的分子动力学方法(Cagin and Pettit) 编辑本段基本步骤 确定起始构型
进行分子动力学模拟的第一步是确定起始构型,一个能量较低的起始构型是进行分子模拟的基础,一般分子的起始构型主要来自实验数据或量子化学计算。

分子动力学
在确定起始构型之后要赋予构成分子的各个原子速度,这一速度是根据波

尔兹曼分布随机生成的,由于速度的分布符合波尔兹曼统计,因此在这个阶段,体系的温度是恒定的。另外,在随机生成各个原子的运动速度之后须 进 行调整,使得体系总体在各个方向上的动量之和为零,即保证体系没有平动位移。 进入平衡相
由上 一步 确定的分子组建平衡相,在构建平衡相的时候会对构型、温度等参数加以监控。 进入生产相
进入生产相之后体系中的分子和分子中的原子开始根据初始速度运动,可以想象其间会发生吸引、排斥乃至碰撞,这时就根据牛顿力学和预先给定的粒子间相互作用势来对各个粒子的运动轨迹进行 计算 ,在这个过程中,体系总能量不变,但分子内部势能和动能不断相互转化,从而 体
系的 温度也不断变化,在整个过程中,体系会遍历势能面上的各个点,计算的样本正是在这个过程中抽取的。 + 计算结果
用抽样所得体系的各个状态计算当时体系的势能,进而计算构型积分。

分子动力学
作用势与动力学计算

作用势的选择与动力学计算的关系极为密切,选择不同的作用势,体系的势能面会有不同的形状,动力学计算所得的分子运动 和 分子内部运动的轨迹也会不同,进而影响到抽样的结果和抽样结果的势能计算,在计算宏观体积和微观成分关系的时候主要采用刚球模型的二体势,计算系统能量,熵等关系时早期多采用Lennard-Jones、morse势等双体势模型,对于金属计算,主要采用morse势,但是由于通过实验拟合的对势容易导致柯西关系,与实验不符,因此在后来的模拟中有人提出采用EAM等多体势模型,或者采用第一性原理计算结果通过一定的物理方法来拟合二体势函数。但是相对于二体势模型,多体势往往缺乏明确的表达式,参量很多,模拟收敛速度很慢,给应用带来很大的困难,因此在一般应用中,通过第一性原理计算结果拟合势函数的L-J,morse等势模型的应用仍然非常广泛。 时间步长与约束动力学
分子动力学计算的基本思想是赋予分子体系初始运动状态之后利用分子的自然运动在相空间中抽取样本进行统计计算,时间步长就是抽样的间隔,因而时间步长的选取对动力学模拟非常重要。太长的时间步长会造成分子间的激烈碰撞,体系数据溢出;太短的时间步长会降低模拟过程搜索
相空间的能力,因此一般选取的时间步长为体系各个自由度中最短运动周期的十分之一。

分子动力学

但是通常情况下,体系各自由度中运动周期最短的是各个化学键的振动,而这种运动对计算某些 宏观性质 并不产生影响,因此就产生了屏蔽分子内部振动或其他无关运动的约束动力学,约束动力学可以有效地增长分子动力学模拟的时间步长,提高搜索相空间的能力。 一般性步骤
以下是做模拟的一般性步骤,具体的步骤和过程依赖于确定的系统或者是软件,但这不影响我们把它当成一个入门指南:
1)首先我们需要对我们所要模拟的系统做一个简单的评估, 三个问题是我们必须要明确的:
做什么(what to do)为什么做(why to do)怎么做(how to do) 2)选择合适的模拟工具,大前提是它能够实现你所感兴趣的目标,这需要你非常谨慎的查阅文献,看看别人用这个工具都做了些什么,有没有和你相关的,千万不要做到一半才发现原来这个工具根本就不能实现你所感兴趣的idea,切记!
考虑1:软件的选择,这通常和软件主流使用的力场有关,而软件本身就具体一定的偏向性,比如说,做蛋白体系,Gromacs,Amber,Namd均可;做DNA, RNA体系,首选肯定是Amber;做界面体系,Dl_POLY比较强大,另外做材料体系,Lammps会是一个不错的选择
考虑2:力场的选择。力场是来描述体系中最小单元间的相互作用的,是用量化等方法计算拟合后生成的经验式,有人会嫌它粗糙,但是它确确实实给我们模拟大系统提供了可能,只能说关注的切入点不同罢了。常见的有三类力场:全原子力场,联合力场,粗粒化力场;当然还有所谓第一代,第二代,第三代力场的说法,这里就不一一列举了。

  1. 能量均分定理- 维基百科,自由的百科全书 - 维基百科- Wikipedia

    1. https://zh.wikipedia.org/zh-hk/能量均分定理

  2. 应用波尔兹曼统计方法可以得到:气体处于平衡态时,分子任何一个自由度的平均能量都相等,均为kT/2,这就是能量按自由度均 ... 圖二顯示的是四種惰性氣體原子速度的麥克斯韋-波茲曼分佈。 ... 均分因此預測有N個粒子的理想氣體有平均總能量(3/2) N kBT。 ..... 沒有n=1的項是因為在平衡點上並無淨力作用,所以能量第一導數為零。
    1. 物理学史- 维基百科,自由的百科全书 - 维基百科- Wikipedia

    1. zh.wikipedia.org/zh-hk/物理学史

    1. 1.1.1 力学的历史背景; 1.1.2 伽利略的动力学; 1.1.3 牛顿三大定律和万有引力定律 .....麦克斯韦玻尔兹曼引入统计诠释之前,热力学始终是基于一组唯象学定律 ...... 狄拉克考虑到薛定谔方程只含对时间的一阶导数而不具有洛伦兹协变性,他从而 ...
    1. 玻尔兹曼 - 爱问知识人- 搜索结果

    1. ishare.iask.sina.com.cn › ... › 自然科学物理

    2. 请教各位高手:忽略掉零点能,晶格振动的能量为En=nhv,利用玻尔兹曼统计(概率f正比 ... 玻尔兹曼系统的微观状态数求最概然分布即极大值,要令它一阶导数=0,二阶导数《0,如何 ... 麦克斯韦——玻尔兹曼能量分布律中,为什么左式可以化为右式 ...
    1. 詹姆斯·克拉克·麦克斯韦_百度百科

    1. baike.baidu.com/view/15809.htm

    1. 经典电动力学的创始人,统计物理学的奠基人之一。 ... 路德维希·玻尔兹曼 ... 匀速运动的势的一阶导数的两个限制,然后由泰勒展开式又得到关于稳定运动二阶导数的 ...
    1. 热力学与统计物理学

    1. 202.207.160.42/jpkc/physics/main/zxtlwtcon.asp

    2. 热力学的方法与统计物理学的方法各有什么. 特点?二者关系 ... 试由热力学基本方程导出八个偏导数麦克斯韦关系, 麦克斯韦关系的作. 用何在? 4. ... 链式关系法” 、 “复合函数微分法” 、 “混合二阶导数法” 的应用。 6. ... 什么是斯忒藩一玻耳兹曼定律?试说明辐射 .... 试导出斯特潘一玻尔兹曼定律、辐射场的状态方程。 36.何谓声子 ...


    1. 由于有显著的量子效应也不能使用麦克斯韦-玻尔兹曼 ... - phymath999

      phymath999.blogspot.com/2013/04/blog-post_3584.html - 轉為繁體網頁

      2013年4月27日 - 麦克斯韦-玻尔兹曼分布是一个概率分布,在物理学和化学中有应用。最常见的应用是 ... 如果 X_i 的分布为 \ X \sim N(0, a^2) ,那么 ... 麦克斯韦-玻尔兹曼分布可以用统计力学来推导(参见麦克斯韦-玻尔兹曼统计)。它对应于由 ..... 阱外波函数为零, 阱内波函数不为零, 它的一阶导数在边界上也不为零。可见, 波函数的一 ...
    1. [DOC]
    2. 量子力学的基本假定

    1. jpkc.qust.edu.cn/index/file/wgx/ja/4.doc

    1. 它们主要由牛顿的经典力学,麦克斯韦的电、磁和光的电磁波理论,玻耳兹曼和吉 ... 许多物理学家试图用经典热力学和统计力学方法解释黑体辐射现象。 ... 在量子化假定基础上,使振子的各本征振动的能量服从玻尔兹曼分布,得到辐射强度与波长的关系 .....二阶导数有意义,则要求要求波函数对坐标的一阶偏导数也必须是连续的。
    1. 华东理工大学插班生考试招生简章 - 恩波教育

    1. cbs.enboedu.cn/colleges/huali.htm

    2. 会计算函数、隐函数(包括由方程组确定)的偏导数二阶导数、全微分。 ... 致冷系数;热力学第二定律、熵和熵增加原理、玻尔兹曼熵关系式;统计规律、理想气体的 ... 的内能、能量按自由度均分定理;麦克斯韦速率分布律、三种统计速率;玻耳兹曼分布; ...
    1. [PDF]
    2. 第六章分子动力学方法

    1. www.bb.ustc.edu.cn/jpkc/xiaoji/jswl/skja/chapter6.pdf

    2. 将上面两式相加、减得到含二阶和一阶导数的公式. (. ) (. )] )(2. [. 1. 2 ..... 常用的初始条件可以选择为:(1)令初始位置在差分划分网格的格子上,初始速度则从玻尔兹曼.

  • No comments:

    Post a Comment