第一性原理

第一性原理的理解及其应用

第一性原理，英文First Principle，是一个计算物理或计算化学专业名词，广义的第一性原理计算指的是一切基于量子力学原理的计算。

我们知道物质由分子组成，分子由原子组成，原子由原子核和电子组成。量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量（或离子），然后就能计算物质的各种性质。

从头算（ab initio）是狭义的第一性原理计算，它是指不使用经验参数，只用电子质量，光速，质子中子质量等少数实验数据去做量子计算。但是这个计算很慢，所以就加入一些经验参数，可以大大加快计算速度，当然也会不可避免的牺牲计算结果精度。

根据原子核和电子互相作用的原理及其基本运动规律，运用量子力学原理，从具体要求出发，经过一些近似处理后直接求解薛定谔方程的算法，习惯上称为第一性原理。

广义的第一原理包括两大类，以Hartree-Fork自洽场计算为基础的ab initio从头算，和密度泛函理论（DFT）计算。也有人主张，ab initio专指从头算，而第一性原理和所谓量子化学计算特指密度泛函理论计算。

第一性原理通常是跟计算联系在一起的，是指在进行计算的时候除了告诉程序你所使用的原子和他们的位置外，没有其他的实验的，经验的或者半经验的参量，且具有很好的移植性。作为评价事物的依据，第一性原理和经验参数是两个极端。第一性原理是某些硬性规定或推演得出的结论，而经验参数则是通过大量实例得出的规律性的数据，这些数据可以来自第一性原理(称为理论统计数据)，也可以来自实验(称为实验统计数据)。

但是就某个特定的问题，第一性原理和经验参数没有明显的界限，必须特别界定。如果某些原理或数据来源于第一性原理，但推演过程中加入了一些假设(这些假设当然是很有说服力的)，那么这些原理或数据就称为“半经验的”。

那为什么使用“第一性原理”这个字眼呢？据说这是来源于“第一推动力”这个宗教词汇。第一推动力是牛顿创立的，因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。如果宇宙诞生之初万事万物应该是静止的，后来却都在运动，是怎么动起来的呢？牛顿相信这是由于上帝推了一把，并且牛顿晚年致力于神学研究。现代科学认为宇宙起源于大爆炸，那么大爆炸也是有原因的吧。所有这些说不清的东西，都归结为宇宙“第一推动力”问题。

科学不相信上帝，我们不清楚“第一推动力”问题只是因为我们科学知识不完善。第一推动一定由某种原理决定。这个可以成为“第一原理”。爱因斯坦晚年致力与“大统一场理论”研究，也是希望找到统概一切物理定律的“第一原理”，可惜，这是当时科学水平所不能及的。现在也远没有答案。

但是为什么称量子力学计算为第一性原理计算？大概是因为这种计算能够从根本上计算出来分子结构和物质的性质，这样的理论很接近于反映宇宙本质的原理，就称为第一原理了。

[1]

第一性原理计算方法的应用

核物理七班张志昌 学号 32100706

1，体系的能量。

进行第一性原理计算前�ナ紫刃枰�确定体系模型�ゼ茨Ｐ偷木О�和晶胞中原子的坐标。对于晶体具有周期对称性�ゾ哂腥�个基矢方向和基矢大小�ゾЦ癯Ｊ��ァＳ捎诶砺奂扑闳范ǖ钠胶饩Ц癯Ｊ�和实验值有所差别�ソ�立模型前需要确定平衡晶格常数。晶格常数的确定采用如下步骤�� 通过改变三个基矢的大小�ジ谋涞グ�的体积��81-119%�ァＭü�第一性原理计算可以得到具有不同体积的模型的能量。通过拟合Murnaghan方程�サ玫骄�体的晶格常数以及单胞的能量�� '0'0 0 00 0 0' '0 0( ) ( ) 1 1( 1)BB V V VE V E V BB B V V�� �オ� �� �� �オ� �オ� �� �� �� �オ� �� �� �オオ� �� �� �オ� �オ� �� 其中��0V为基态平衡体积��0 0( )E V为基态下体系的结合能�ハ喽杂诙杂�孤立原子能量�ァ�V为原胞体积��0B为体模量��'0B为体模量对压强的导数。如课件中图形所示�タ梢匀范ㄔ谝欢ㄌ寤�下体系的能量达到极小值�ゼ刺逑档幕�态能量�ニ�对应的体积为体系的平衡体积�ソ�而可以得到模型三个基矢的大小确定晶体的平衡晶格常数。这里需要指出的是不同的第一性原理计算方法给出的能量�ゴ�表的物理意义不同�サ�是本质上都可以反应体系的稳定性。如总能指构成体系的原子孤立时的能量减去原子成键放出的能量�ソ岷夏苁且怨铝⒃�子的能量为零点�ヌ逑稻哂械淖苣塥ゼ丛�子构成晶体时放出的能量。 在上面求得的晶格常数的基础上�ジ�据要研究的物理问题�ト范ㄌ逑抵邪�含原子数目的多少�ソ�立第一性原理计算模型。第一性原理计算的模型通常选取一个或几个单胞�コ�单胞�プ魑�模型�パ∪〉哪Ｐ途哂腥�个基矢方向�ビΡＶぱ刈湃�个基矢方向平移可以构成无限大的晶体。

第一性原理计算输入的原子坐标有两种坐标形式�ヒ恢质堑芽ǘ�坐标��Cartesian coordinates�オヒ恢质欠质�坐标��fractional coordinates�ァＨ缍杂贜i3Al高温合金�ゾ哂腥缤妓�示的晶体结构�ヂ猎�子位于立方体的顶点�ツ�原子位于立方体的面心位置。如果取一个单胞作为研究模型�ピ蛉�个基矢1a�オ�2a�オ�3a�シ直鹞���a, 0 , 0�オ�0, a, 0�オ�0, 0, a�オテ渲衋为体系的晶格常数。单胞中包含四个不等价原子�ト�个Ni和一个Al。采用笛卡尔坐标四个原子的坐标可表示为��0��0��0�オオ�a/2, a/2, 0�オオ�a/2, 0 , a/2�オオ�0, a/2, a/2�ァＨ绻�采用分数坐标表示为��0��0��0�オオ�1/2, 1/2, 0�オオ�1/2, 0 , 1/2�オオ�0, 1/2, 1/2�ァ５峡ǘ�坐标(x,y,z)和分数坐标(a,b,c)之间关系为1a�ァ�a+1a�ァ�b+1a�ァ�c=(x, y, z)�テ渲衋,b,c为一个原子的三个分数坐标��x,y,z为该原子的笛卡尔坐标。图中所示各点表示将晶格常数的大小a取不同值时得到的单胞体积作为横轴�ザ�纵轴表示对应体积下将原子坐标输入进行第一性原理计算求得的体系能量。拟合后得到Ni3Al的平衡基矢大小以及体系能量。 对于研究合金中的掺杂问题�ビ捎诓粼釉�素的量很少�ニ�以建立的模型需要取多个单胞�コ�晶胞�ァＫ孀拍Ｐ椭性�子数目的增加�サ谝恍栽�理计算方法的计算量指数增加�ザ杂诓粼恿亢艿偷那榭霆ト�0.1%�バ枰�模型中至少取1000个原子来和实际相符合�フ獬�出了第一性原理计算在目前的计算机上的计算能力��100原子左右�オニ�以建立模型时需要考虑能够反映要研究的实际问题就可以。假设一个超单胞中只存在一个掺杂原子�フ庋�相邻两个超晶胞中掺杂原子的间距为超晶胞的基矢大小。一般两个原子之间相隔三到四个原子层�ピ�子之间的相互作用就可以认为非常小了。所以选取八个单胞构成的超单胞就可以基本反应掺杂量很低的掺杂问题了。 将建立好的模型�ソ�行第一性原理计算可以得到体系的总量�ザ宰苣芙�行变换可以定义体系的内聚能、形成能以及择优占位能�ソ�而可以对掺杂是否有利于形成�バ纬刹粼雍蠖蕴逑滴榷ǖ挠跋於�进行分析�ト缈渭�中所列。内聚能�ヌ逑档淖苣芗跞ニ�有原子孤立时的能量�ゼ从捎谠�子之间的相互作用而放出能量�ゴ佣�内聚能为负值�ピ叫”硎拘纬傻奶逑翟轿榷āＭü�和没有掺杂体系的内聚能相比较可以看出掺杂元素对体系稳定性的影响。形成能�ヌ逑档淖苣芗跞ヌ逑抵懈髯栽�素对应的晶体中原子的能量。形成能可表示各种金属组成合金的能力。另外通过比较掺杂原子替代合金中不同元素原

核物理七班张志昌 学号 32100706

子时体系的形成能可以得到掺杂原子倾向位于合金的什么位置�フ飧霾钪稻涂梢远ㄒ逦�择优占据能。这里需要强调�ジ髦帜芰渴歉�据要研究的具体问题来定义的�ケ热缥颐且�就掺杂原子倾向于位于合金的什么位置�ナ褂米苣苁遣荒艿玫姜ヒ蛭�超单胞模型中各种原子的数

目不相同�ザ�每种原子的能量是不一样的�ッ挥锌杀刃元ニ�以定义了择优占位能。

2、电荷密度

电荷密度就是晶体中电子密度的分布。通过电荷密度可以知道晶体中原子间的成键状态�ト缃鹗艏�、共价键、离子健、van der Waals键和氢键。为了更好的表示原子形成晶体后原子间的电荷转移和成键情况�ヒ�入差分电荷密度�ゼ戳礁鎏逑抵械绾擅芏鹊牟钪怠Ｕ饬礁鎏逑涤Ω镁哂邢嗤�的超单胞�コ�单胞中原子类型可以不一样�ザ�原子位置要基本一致�ト缈渭�中所示�ビ肗i3Al中一个Ni被掺杂元素替代时体系的电荷密度减去没有掺杂的Ni3Al的电荷密度�ザ�得到差分电荷密度�ネü�图形可以清楚地看出由于掺杂元素的存在导致的电

子分布状态的改变。再如�ノ颐强梢越玁i3Al晶体的电荷密度减去Ni和Al原子放在超单胞相同位置时孤立Ni和Al原子的电荷密度�タ梢缘玫絅i3Al中Ni和Al原子成键过程中电子密度分布的变化�ゴ佣�更好的观察原子之间的成键情况。所以差分电荷密度相当于一个体系的电荷密度相对于另一个体系的电荷密度改变�ツ康氖俏�了更好的研究体系中的成键状态。

3、能带

能带理论是目前研究固体中电子运动的一个主要理论基础。是在用量子力学研究金属电导理论的过程中开始发展起来的。最初的成就在于定性地阐明了晶体中电子运动的普遍性的特点�ダ�如固体为什么会有导体、非导体的区别�ゾ�体中电子的平均自由程为什么会远大于原子的间距等。在半导体技术上�ツ艽�论提供了分析半导体理论问题的基础�ビ欣�的推动了半导体技术的发展。 能带理论是一个近似的理论。在固体中存在大量的电子�ニ�们的运动是相互关联的�ッ扛龅缱拥脑硕�都要受其他电子运动的牵连。价电子是人们最关心的。在原子结合成固体的过程中价电子的运动状态发生了很大的变化�ザ�内层电子的变化比较小�タ梢园言�子核和内层电子近似看成是一个离子实�ゼ鄣缱涌梢钥闯稍谝桓龅刃�势场中运动。能带理论的出发点是固体中的电子不再束缚于个别原子�ザ�是在整个固体内运动�コ莆�共有化电子。电子的能量状态从处于一个电子能级变到在一个能量范围内都会存在。我们下面从对自由电子的能量讨论�サ贸瞿艽�的表示方法。 自由电子的能量�ザ�能��2/ 2kE E p m= =�ジ�据德布罗依波长与动量关系��/p h kl= = h�テ渲衕为普朗克常量��h为/ 2hp��2 /kp l=成为波数或波矢。所以电子能量可表示为2 2/ 2E k m= h�ゼ吹缱幽芰课�波矢k的函数�ト缈渭�中图形所示。我们讨论固体中电子的能量通常是在k空间�サ箍占洫ツ诮�行的。 对于晶体具有周期对称性�テ涠杂Φ牡箍占湟簿哂兄芷诙猿菩元ザ杂谝晃�的情况�ナ悼占涞闹芷谖猘�ピ虻箍占涞闹芷谖�2 /ap�ザㄒ宓谝桓鲋芷讵�/ap-��/ap�ノ�第一布里渊区。则电子能量在实空间分布的随晶体的周期对称性变化�プ�化到k空间�サ缱幽芰克娴箍占涞闹芷诙猿菩苑⑸�变化�ゼ吹缱幽芰吭诘谝徊祭镌ㄇ�随k的变化在整个k空间中周期性重复。因此可以得到课件中能量随k变化的图形�ハ嗟庇谧杂傻缱拥哪艽�经过周期性势场调制后的结果。所以我们讨论能带只需要考虑电子能量在第一布里渊区随k的变化关系就可以。由于计算所取晶体超单胞形状的不同会导致第一布里渊形状的变化�ヒ约熬�体中原子化合状态的不同导致电子所受到的周期性势场不同�ニ�以电子能量随k的变化关系在不同的晶体中是

不同的�ゼ床煌�的晶体具有不同的能带结构�ゴ佣�反映出不同的物理性质。

4、电子状态密度

由前面的讨论可知晶体中电子能量状态可以取一定的能量范围�ピ诖四芰糠段�内在不

核物理七班张志昌 学号 32100706

同的能量区间电子能量状态的多少或填充这些能量状态的电子数目是不一样的。电子状态密度反映了这个不同�ゼ丛谝荒芰壳�间内电子状态�ナ�目�サ亩嗌侏�/d dEh��dh为在dE能量区间内电子能量状态的数目。通过对体系电子状态密度的分析可以得到晶体中原子间的电子杂化情况。

最后感谢老师这学期的课程，的确开阔了我的视野，让我对量子力学有了初步的了解。我会继续学习计算物理的相关知识，巩固这学期学到的知识。

核物理七班张志昌

学号 32100706