第1性原理计算简述.docx

第一性原理计算简述 第一性原理，英文First Principle，是一个计算物理或计算化学专业名词，广义的第一性原理计算指的是一切基于量子力学原理的计算。我们知道物质由分子组成，分子由原子组成，原子由原子核和电子组成。量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量（或离子），然后就能计算物质的各种性质。从头算（ab initio）是狭义的第一性原理计算，它是指不使用经验参数，只用电子质量，光速，质子中子质量等少数实验数据去做量子计算。但是这个计算很慢，所以就加入一些经验参数，可以大大加快计算速度，当然也会不可避免的牺牲计算结果精度。那为什么使用“第一性原理”这个字眼呢？据说这是来源于“第一推动力”这个宗教词汇。第一推动力是牛顿创立的，因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。如果宇宙诞生之初万事万物应该是静止的，后来却都在运动，是怎么动起来的呢？牛顿相信这是由于上帝推了一把，并且牛顿晚年致力于神学研究。现代科学认为宇宙起源于大爆炸，那么大爆炸也是有原因的吧。所有这些说不清的东西，都归结为宇宙“第一推动力”问题。科学不相信上帝，我们不清楚“第一推动力”问题只是因为我们科学知识不完善。第一推动一定由某种原理决定。这个可以成为“第一原理”。爱因斯坦晚年致力与“大统一场理论”研究，也是希望找到统概一切物理定律的“第一原理”，可惜，这是当时科学水平所不能及的。现在也远没有答案。但是为什么称量子力学计算为第一性原理计算？大概是因为这种计算能够从根本上计算出来分子结构和物质的性质，这样的理论很接近于反映宇宙本质的原理，就称为第一原理了。 　　广义的第一原理包括两大类，以Hartree-Fork自洽场计算为基础的ab initio从头算，和密度泛函理论（DFT）计算。也有人主张，ab initio专指从头算，而第一性原理和所谓量子化学计算特指密度泛函理论计算。 根据原子核和电子互相作用的原理及其基本运动规律，运用量子力学原理，从具体要求出发，经过一些近似处理后直接求解薛定谔方程的算法，习惯上称为第一???理。第一性原理通常是跟计算联系在一起的，是指在进行计算的时候除了告诉程序你所使用的原子和他们的位置外，没有其他的实验的，经验的或者半经验的参量，且具有很好的移植性。作为评价事物的依据，第一性原理和经验参数是两个极端。第一性原理是某些硬性规定或推演得出的结论，而经验参数则是通过大量实例得出的规律性的数据，这些数据可以来自第一性原理(称为理论统计数据)，也可以来自实验(称为实验统计数据)。但是就某个特定的问题，第一性原理和经验参数没有明显的界限，必须特别界定。如果某些原理或数据来源于第一性原理，但推演过程中加入了一些假设(这些假设当然是很有说服力的)，那么这些原理或数据就称为“半经验的”。 最早期的从头算主要是一些量子化学的专家在研究。其实这还要从头说起，薛定谔搞出了薛定谔方程。方程写起来虽然简单，而且直观。（动能算符+势能算符）*波函数=能量的本征值*波函数。这是定态波函数，这个方程主要是用来求基态。方程虽然容易写，但是求解起来实在是太难了。有点像经典力学里的混沌现象。虽然每个粒子都服从牛顿三大定律，但是求解是不可能的。首先，求解的难题是电子和核的相互作用项难以分开。还好奥本海默提出了绝热近似，即核的质量和电子的质量不是一个数量级的，所以速度也不是一个数量级的。认为核缓慢的能跟上电子的运动。这样，我们即可将核的方程和电子的方程分开。而在真实的材料中，电子的作用是很大。例如成键主要是指电子和电子的相互作用。因此我们只需要仔细求解电子的薛定谔方程即可。对于N个电子的系统，方程是3N维的。其求解仍然很难。所以，科学家们就想办法。最后的办法是，利用自洽的方法求解。自洽其实是一个很有意思的办法，说白了，就是指，A和B有数学的依赖关系，但A和B的方程是解不出来的。给出一个A我可以求B，有B又反求出一个A，这样循环下去直到|A_new - A_old|无穷小，或者|B_new -B_old|无穷小。方程自洽，运算结束。这只是简单的举例。首先，这个由hartree和fock提出了hartree-fock近似，他们将多电子的波函数写成了单个电子波函数的积的形式，考虑了全同粒子的交换性，而哈密顿量写成了密度的泛函，进而求得基态能量，可以说这是一个很大的进步。但确定是没有考虑电子间的关联作用。接着，由kohn等人进一步改善了薛定谔方程，首先他们认为基态可以写成电子密度的泛函，真实的材料的所有性质都可以有电子密度求出。这就是密度泛函理论，详细的推导还是很繁琐的。k-s方程中，动能是用的是无关联的自由电子动能，势能用的是无交换的电子相互作用能和离子对电子的势能，剩下的就是交换关联能。虽然方程是相当完美的，也很准确，但