第一性原理计算
• Lagrange乘法因子被写为-ij
E ij Sij 0
i j
• 多电子系统和单电子系统的最大区别在于存在电子之间的 作用,包括库仑和交换积分 • 找到多电子系统最好的波函数,保留轨道图像 • 找到一个解,当一个电子的轨道变化时由于电子耦合作用 会影响其他自旋轨道上的电子的运动,在上述情况下,这 个解能保证所考虑的电子的运动
4.3 闭核系统的能量
•在分子模型中,关心分子的基态能量
•他们大部分具有闭壳系统
•N/2个轨道,N个电子的闭壳系统 •每一个空间轨道i具有两个自旋轨道:i,i
• Hiicore:电子在裸露原子核作用场中运动的能量。 对于在轨道i运动的电子:如果一个轨道具有两个电子,那 么对N/2个电子而言单个电子能量为2 Hiicore.
•每一个自旋轨道写为单个电子轨道的线性组合
i c vi v
v 1 k
•一个电子的波函数v称为基函数,常对应原子轨道. •将用希腊符号来标记.
•共有K个基函数。我们希望得到K个分子轨道
•根据变分思想,我们得到对应最低能量的系数
•Cvi通过改变系数得到波函数
•对于给定的基函数以及波函数的形式,最好的系数应满足:能量最小原理。
1 1
K
K
闭壳系统Fock矩阵中每一项元素为
F H P [2(uv ) ( )]
core
K
K
1 1
5.4 解Roothann-Hall方程
•Fock矩阵微K乘k阶矩阵,如果运用实函数,则是对称的 •Roothaan-Hall方程可以写为矩阵形式
对给定的轨道,屏蔽常数有下列贡献之和确定 • 如果到原子核的距离比到上述各组远,取0; • 同一组的其他电子贡献为0.35;但若其他轨道为1s,则贡献为0.3。 • 电子的主量子数比现有轨道小2或更多,则贡献为1。 • 如果每一个电子具有的主量子数比现有轨道小1:现有轨道为d或f,取 1.0;现有轨道为s,p取0.85。
core E 2 H ii (4 J ij 2 K ij ) J ij i 1 i 1 j i 1 i 1 N /2 N /2 N /2 N /2
• jii=Kii,
E 2 H
i 1
N /2
core ii
(2 J ij K ij )
i 1 j 1
Fi i ij j
j 1
N
Fi (1) H
core
(1) J j (1) K j (1)
j 1
Fock算符 :多电子系统有效的哈密顿算符
对于闭核系统
Fi (1) H
core
(1) 2 J j (1) K j (1)
j 1
N /2
• 除非i等于j,否则Lagrange乘法因子为0 .
v 1 v 1
K
K
每一侧都乘以u(1),(也是一基函数),积分
c d (1)F (1) (1) c d (1) (1)
v 1 vi 1 i v i v 1 vi 1 v
K
K
d 1 (1)v (1)
Fuv d 1 (1) Fi (1)v (1)
•两个基函数之间的交叠积分,写为Suv •两个基函数的交叠项并不一定为0
Fock矩阵
闭壳系统的Fock矩阵元素
Fuv d 1 (1) H
core
(1)v (1)
j 1
N /2
d (1)[2 J
1
j
(1) K j (1)](1)v (1)
Fock矩阵地元素可以写为心部,库仑作用,和交换能项的综合
d (1)H
1
core
ZA 1 2 M (1)v (1) d 1 (1)[ ](1)v (1) 2 A1 r1 RA
Huvcore每一个元素可以通过对应单个电子哈密顿项的动能积分和势能积分得到 库仑以及交换作用项
N /2 j 1
d (1)[2J
N /2N /2
5 hartree-Fock方程
• 变分方法 • 真实函数近似计算的能量总是高于真实函数得到的能量。 • 波函数越好,则能量越低。 • 等能量最低时,得到最好的波函数。 • 在最小点时,能量的一次微分,E=0.
• 这一条件作用于能量表达式
,得到Hartree-Fock方程
• 正交归一化条件可以写为交叠矩阵Sij
5.1 原子的H-F计算以及Slater规则
• 如果假设电子的分布实球形对称的 • 近似分析函数
Rnl (r)Ylm ( , )
•氢原子得到的径向函数不能直接用于多电子原子。
•因为内壳层电子对原子核电荷具有屏蔽作用。
•如果考虑屏蔽效应把轨道指数作适当地修正,仍可以采用氢原子的波函 数形式。
•交叠矩阵中非对角线元素可能不为0
•为了用标准方法解R-H方程,必须对之进行变换。 这相当于变换基函数,使之正交归一化。 •寻找一个矩阵X,使XTSX=I. •对称正交归一化条件下交叠矩阵是对角化的
U T SU D diag(1,k )
E 0 cvi
5.3 闭壳系统和Roothann-Hall方程
• 具有N/2轨道N个电子的闭壳系统
• 系统的H-F方程最早由Roothaan和Hall独立提出所谓的 Roothaan或Roothann-Hall方程 • 矩阵形式,它可以用标准的方法求解
E 2 H
i 1
N /2
core ii
• Slater
Rnl (r) (2 )n1/ 2 [(2n)!]1/ 2 r n1 er
R1s (r) (2 )3 / 2 e r
R2s (r) R2 p (r) (4 5 / 3)3/ 2 re r
R3s (r) R3 p (r) R3d (r) (8 7 / 45)1/ 2 r 2er
FC=SCE E为由轨道能量组成的对角矩阵
c1,1 c 2 ,1 C ck ,1
c1, 2 c1,k c 2 ,k c k ,k
1 0 0 0 0 0 2 E 0 0 k
•标准本征方法需要矩阵形式FC=CE •R-H方程中只有交叠矩阵S取单位矩阵才为标准形式
1
j
(1) K j (1)](1)v (1)
库仑算符Jj(1)由自旋轨道j的交互作用引起
J j (1) d 2 j (2)
1 j ( 2) r12
利用轨道函数的基函数线性组合的方法
1 k J j (1) d 2 cj (2) cj (2) r12 1 1
Z n*
n* :
Z为原子数,为屏蔽常数,n*为有效的主量子数
n=1,2,3时取与n相同的值;n=4,5,6s时,取3.7,4.0,4.2 屏蔽常数按下列规则确定 : 把轨道划分如下组: (1s);(2s,2p);(3s,3p);(3d);(4s,4p);(4f);(5s,5p);(5d)
• 假设一个电子位于原子核作用场中一个自旋轨道I中, 其他电子位于轨道j中
1 2 M ZA 1 1 ( 1 ) [ d ( 2 ) ( 2 ) ] ( 1 ) [ d ( 2 ) ( 2 ) ] i (1) ij j (1) i 1 2 j j i 2 j i 2 r r r A1 iA j i j i j 12 12
交换算符
H core (1) i (1) J j (1) i (1) K j (1) i (1) ij j (1)
j i j i i
N N core H (1) J j (1) K j (1) i (1) ij j (1) j 1 j 1
N /2 i 1 core 2 H ii
• 电子与电子的作用项
• i和j中存在4个电子。 • 在一个轨道上的2个电子以库仑作用的方式与另一个轨道 上的两个电子发生作用。记为4Jij。 • 但是在这些组合中,包含了成对电子的两种作用方式,给 出了总的交换能-Kij。 • 最后,在同一轨道上的一对电子库仑作用也要考虑在内, 但是这些电子具有成对的自旋因此不存在交换能。
• 完整的轨道,我们要用角度部分乘以径向函数部分
1s ( r ) 3 / exp( r )
2 s ( r ) 5 / 3 r exp( r )
2 p ( r ) 5 / exp( r ) cos
z
• Slater给出了一系列的轨道指数的确定原则
S ij i j d ij
• 能量最小的限制条件可以通过Lagrange因子解决 • 最小化的函数的微分加上限制条件乘以Lagrange乘法因子 的微分。其和设置为0 • 在每一个正交归一条件下的Lagrange乘法因子记为ij
E ij Sij 0
i j
cj cj [ 2(uv ) ( )]
j 1 1 1
N /2 K
通过引入密度矩阵
P 2 ci ci
11 N /2
计算某一点的电子密度
P 2 ci ci
11
N /2
( r ) P ( r ) ( r )
• 硅的价电子的屏蔽常数计算 电子结构为(1s2)(2s22p6)(3s23p2) • 规则b得到30.35; • 根据规则(c),得2.0; • 根据规则(d),得80.85;
• 我们得到原子核数为14时,Z-为4.15.
5.2 H-F方程中的原子轨道线性组合的方法(LCAO)
第一性原理计算在金属材料研究中的应用
第一性原理计算在金属材料研究中的应用在过去,金属材料的研究主要依靠实验来进行。
而如今,第一性原理计算已经成为了一种新的技术,可用于模拟金属材料的结构、性质和反应。
因此,第一性原理计算已经成为金属材料研究的一种重要工具。
本文将介绍第一性原理计算在金属材料研究中的应用及其优点与局限性。
一、第一性原理计算简介第一性原理计算是使用量子化学理论以及密度泛函理论来计算材料的性质。
据此,材料的电子结构和固有性质可以直接从基本定律得到。
这直接破除了传统材料科学中需要大量实验和经验来建立新材料的做法。
二、第一性原理计算在金属材料研究中的应用第一性原理计算玩家能够提供一个完整的金属材料体系,其中包含各种金属结构以及它们特有的热力学、电子、机械和磁学特性。
第一性原理计算还可以提供材料之间相互作用的关键细节,以及元素和合金的更好理解。
这样一个完整的体系,可以用来预测结构、构确性质和跟踪反应。
以下介绍了第一性原理计算在金属材料研究中所扮演的具体角色:1.预测材料性质:第一性原理计算可以预测材料的结构和电子性质,包括能带、电荷密度分布、电子能级结构和振动特性等。
这种预测使得研究者可以更好地了解材料的性质和反应。
2.设计新材料:第一性原理计算可以预测新合金或材料的性质,并提供一些重要信息,例如新材料的制造条件和可能发生的反应等。
3.优化现有材料:第一性原理计算也可以用于优化现有材料的物理和化学性质,以提高制造效率和性能。
三、第一性原理计算的优点与局限性1.优点(1)准确性高:第一性原理计算可以从基本原理出发精确计算材料的性质。
(2)可重复性强:第一性原理计算的结果可以通过重复实验得到。
(3)节省经费和时间:相比于传统实验,第一性原理计算更加经济高效,减少了材料研究的成本和时间。
2.局限性(1)计算复杂度:第一性原理计算需要处理极其复杂的数学和物理理论,并且需要耗费大量时间来进行计算。
(2)计算结果与实验结果相差较大:由于第一性原理计算过于理论化,因此在与实验结果进行对比时,会有一些误差。
第一性原理计算
第一性原理计算引言第一性原理计算是一种基于量子力学原理的计算方法,用于研究材料的性质和行为。
它通过解析薛定谔方程,从头开始计算材料的性质,而不依赖于经验参数或已知的实验数据。
这使得第一性原理计算成为研究材料性质的重要工具,也为材料设计和开发提供了新的途径。
原理和方法第一性原理计算的核心是薛定谔方程的求解。
薛定谔方程描述了量子力学系统的行为,通过求解薛定谔方程可以得到体系的能量、电子结构、晶体结构、力学性能等信息。
然而,薛定谔方程的精确求解是不可行的,因此需要使用一些近似方法来简化计算过程。
其中最常用的方法是密度泛函理论(DFT)。
密度泛函理论的基本思想是将体系中的电子密度视为基本变量,通过最小化体系的总能量来确定电子密度。
这可以通过Kohn-Sham方程来实现,其中包括了交换-相关能的近似处理。
通过求解Kohn-Sham方程,可以得到体系的电子结构和能量。
此外,还有一些其他的方法被用于提高计算精度,如GW近似、自洽Poisson方程、多体微扰理论等。
这些方法的选择取决于研究问题的特点和需要。
应用领域第一性原理计算在材料科学、物理学和化学等领域有着广泛的应用。
1.材料设计:第一性原理计算可以用于预测新材料的性质,从而加速材料的设计和开发过程。
它可以通过计算和优化材料的能带结构、晶体结构等来寻找具有特定性能的材料。
2.反应动力学:第一性原理计算还可以用于研究化学反应的动力学过程。
通过计算反应的势能面和反应路径,可以预测反应速率和产物选择性。
3.催化剂设计:催化剂是许多化学反应中的关键组分。
第一性原理计算可以帮助设计和优化催化剂的表面结构和活性位点,从而提高催化剂的效率和选择性。
4.电子器件:第一性原理计算在电子器件领域的应用也日益重要。
它可以用于模拟和优化半导体器件的性能,如晶体管、太阳能电池等。
5.生物物理学:第一性原理计算在生物物理学研究中也发挥着重要作用。
它可以用于预测蛋白质的结构和稳定性,研究生物分子的相互作用以及药物分子的设计等。
第一性原理计算
5.1 原子的H-F计算以及Slater规则
• 如果假设电子的分布实球形对称的 • 近似分析函数
Rn(lr)Ylm (,)
•氢原子得到的径向函数不能直接用于多电子原子。 •因为内壳层电子对原子核电荷具有屏蔽作用。 •如果考虑屏蔽效应把轨道指数作适当地修正,仍可以采用氢原子的波函 数形式。
• Slater
第一性原理计算
• Hiicore:电子在裸露原子核作用场中运动的能量。
• 对于在轨道 i运动的电子:如果一个轨道具有两个电子, 则对N/2个电子而言单个电子能量为2 Hiicore.
•
N /2
2
H
c ii
ore
i 1
• 电子与电子的作用项
• i和j中存在4个电子。
• 在一个轨道上的2个电子以库仑作用的方式与另一个轨道 上的两个电子发生作用。记为4Jij。
•需要猜想密度矩阵P。最简单的方法是采用空矩阵 •H-F计算的结果是一系列K原子轨道,K是计算的基函数的数量 •N个电子被填充到这些轨道上,从最低能量的轨道开始,一个轨道2个电子进行填充。
H-F方程给出了一系列的轨道轨道能量,i,
N/2
i Hiciore (2JijKij)
j1
整个基态的电子能量
N/2
d1(1)2 [Jj(1)Kj(1)]1)(v(1)
j1
库仑算符Jj(1)由自旋轨道j的交互作用引起
Jj(1) d2j(2)r112j(2)
利用轨道函数的基函数线性组合的方法
k
Jj(1)d2 cj
1
1k (2)r12 1cj
(2)
交换项可以写为
K j(1 )i(1 ) [d2 k 1 c j (2 )r 1 12j(2 )k ] 1 c j (2 )
化学反应动力学的第一性原理计算方法
化学反应动力学的第一性原理计算方法化学反应是物质变化的一种形式,通常是指原子或分子之间的化学键被打破或形成,从而形成新的化合物。
化学反应动力学研究的是化学反应速率的研究,也就是反应物转变为产物的速率。
动力学的研究对于理解化学反应机理和制定化学反应工艺有着重要的意义。
在现代化学研究中,化学反应动力学的第一性原理计算方法已经成为重要的工具。
化学反应动力学的第一性原理计算方法指的是运用量子力学原理和分子动力学模拟技术对化学反应动力学过程进行精确的计算和模拟。
这种方法无需依靠实验数据,而是直接从微观层面分析分子之间的相互作用。
通过对分子结构和动力学过程的分析,可以计算得到反应动力学的速率常数、反应机理、反应能垒等其它重要参数,从而能够深入理解化学反应的本质。
化学反应动力学的第一性原理计算方法主要应用于分子动力学模拟和量子化学计算两个方面。
其中,分子动力学模拟方法主要是基于原子力场,通过数值积分求解牛顿方程,模拟反应过程。
它可以计算物质的结构、能量以及动力学过程。
量子化学计算方法则是基于量子力学理论,通过求解薛定谔方程,计算分子间的相互作用和反应机理。
这种方法可以计算各种化学反应的能垒、活化能、反应速率以及反应机理。
对于化学反应动力学的第一性原理计算方法,其中一个比较重要的问题就是如何评估理论计算的准确性。
实际上,在计算化学的过程中,化学反应动力学的第一性原理计算方法也不能完全避免计算误差。
因此,如何评估计算误差以及如何优化理论计算模型是这个领域研究者一直在关注的问题。
面对以上问题,化学反应动力学的第一性原理计算方法的研究者们借鉴了机器学习的思想,开发出了一种基于数据库和机器学习的化学反应动力学数据驱动模型。
该模型基于已有的关于反应动力学的实验数据和理论计算数据,通过机器学习方法对反应动力学模型进行训练、验证和优化。
这种模型可以有效地降低计算误差,提高计算准确性,并能够提高计算速度。
总之,化学反应动力学的第一性原理计算方法是指在量子力学和分子动力学的基础上,通过计算和模拟分子间的相互作用和反应过程来研究化学反应动力学的的方法。
第一性原理计算是什么意思
第一性原理计算是什么意思简介第一性原理计算(First Principles Calculation)是一种基于量子力学原理的计算方法,用于研究材料和分子的性质及其相互作用。
通过求解薛定谔方程,第一性原理计算可以预测和解释材料的结构、能量、电子结构、磁性、光学性质等。
这种计算方法是建立在非经验的基础上,仅依赖于原子核和电子之间的相互作用,因此被称为“第一性原理”。
原理第一性原理计算的基础是量子力学中的薛定谔方程。
该方程描述了粒子的行为,并可以用于计算材料的性质。
在第一性原理计算中,薛定谔方程被用来描述系统的电子结构,通过求解薛定谔方程,可以得到材料的电子能级、原子间的相互作用等信息。
第一性原理计算基于密度泛函理论(Density Functional Theory,DFT),该理论通过体系的电子密度来描述材料的电子结构。
根据克斯特兰–库尔(Hohenberg-Kohn)定理和克斯特兰–库尔–夏姆(Kohn-Sham)方程,DFT可以将多体问题简化为一个单体问题,使得计算变得可行。
薛定谔方程的求解需要进行数值计算,常用的方法包括平面波基组法(Plane Wave Basis Set)和赝势法(Pseudo-potential Method)。
平面波方法将波函数展开为平面波的线性组合,可以较好地描述材料的周期性结构。
赝势方法则通过引入有效势能的概念,去除了原子核与内层电子的相互作用,从而大大简化了计算。
应用第一性原理计算可以应用于许多领域,尤其在材料科学和化学领域中发挥着重要作用。
1.新材料的设计与发现:通过第一性原理计算,可以预测新材料的结构稳定性、电子结构、能量等性质,从而指导新材料的设计与合成。
例如,通过计算优选的材料组合,可以设计出具有特定电子结构和物理化学性质的材料,如催化剂、光电材料等。
2.催化剂的研究与设计:第一性原理计算可以揭示催化反应中的活性位点和反应机理,从而指导催化剂的设计和优化。
第一性原理计算简述
第一性原理计算简述第一性原理,英文Firs t Principle,是一个计算物理或计算化学专业名词,广义的第一性原理计算指的是一切基于量子力学原理的计算。
我们知道物质由分子组成,分子由原子组成,原子由原子核和电子组成。
量子力学计算就是根据原子核和电子的相互作用原理去计算分子结构和分子能量(或离子),然后就能计算物质的各种性质。
从头算(ab initio)是狭义的第一性原理计算,它是指不使用经验参数,只用电子质量,光速,质子中子质量等少数实验数据去做量子计算。
但是这个计算很慢,所以就加入一些经验参数,可以大大加快计算速度,当然也会不可避免的牺牲计算结果精度。
那为什么使用“第一性原理”这个字眼呢?据说这是来源于“第一推动力”这个宗教词汇。
第一推动力是牛顿创立的,因为牛顿第一定律说明了物质在不受外力的作用下保持静止或匀速直线运动。
如果宇宙诞生之初万事万物应该是静止的,后来却都在运动,是怎么动起来的呢?牛顿相信这是由于上帝推了一把,并且牛顿晚年致力于神学研究。
现代科学认为宇宙起源于大爆炸,那么大爆炸也是有原因的吧。
所有这些说不清的东西,都归结为宇宙“第一推动力”问题。
科学不相信上帝,我们不清楚“第一推动力”问题只是因为我们科学知识不完善。
第一推动一定由某种原理决定。
这个可以成为“第一原理”。
爱因斯坦晚年致力与“大统一场理论”研究,也是希望找到统概一切物理定律的“第一原理”,可惜,这是当时科学水平所不能及的。
现在也远没有答案。
但是为什么称量子力学计算为第一性原理计算?大概是因为这种计算能够从根本上计算出来分子结构和物质的性质,这样的理论很接近于反映宇宙本质的原理,就称为第一原理了。
广义的第一原理包括两大类,以Hartr ee-Fork自洽场计算为基础的abinitio从头算,和密度泛函理论(DFT)计算。
第一性原理计算方法讲义
第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。
而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。
第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。
量子力学是反映微观粒子运动规律的理论。
量子力学的出现,使得人们对于物质微观结构的认识日益深入。
原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。
量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。
以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。
目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。
但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。
绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。
Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。
但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。
1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。
它建立在非均匀作为基本变量。
1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。
第一性原理计算在物理实验中的应用案例
第一性原理计算在物理实验中的应用案例引言:第一性原理计算是一种通过量子力学和电子结构理论对物质进行全面建模和计算的方法。
它基于物质的基本原子结构和电子行为进行计算,不依赖于任何经验参数或已知实验数据。
在物理实验中,第一性原理计算已经成为一种强大的工具,可以帮助研究人员预测材料的性质、理解实验现象的本质,甚至设计更好的实验方法。
本文将通过几个案例介绍第一性原理计算在物理实验中的应用。
1. 电子结构计算在新材料研究中的应用近年来,第一性原理计算在新材料研究领域得到了广泛应用。
例如,在半导体材料研究中,通过计算电子结构和能带结构,可以预测材料的导电性、光学性质和磁性行为。
研究人员可以利用这些计算结果指导实验设计,合成具有特定性能的新材料。
此外,第一性原理计算还可以帮助解释实验中观察到的奇特现象,如高温超导现象。
通过计算电子结构和电子-声子相互作用,可以理解超导材料中电子配对的机制,从而为实验提供更深入的理论基础。
2. 力学性质计算在纳米材料研究中的应用随着纳米技术的发展,纳米材料的研究成为热门领域。
第一性原理计算在研究纳米材料的力学性质方面发挥了重要作用。
通过计算纳米材料的力学性质,如应力-应变关系、杨氏模量等,可以预测纳米材料在不同应力下的力学行为。
这对于合成更强韧的纳米材料、设计新型传感器和制备高性能纳米器件具有重要价值。
此外,第一性原理计算还可以帮助研究人员解决纳米材料在实验中的制备和加工过程中的问题,指导实验进行更精确的参数调控。
3. 表面反应计算在催化剂研究中的应用催化剂在化学工业生产中起着重要作用,催化剂表面的原子结构和反应机制对于催化过程的理解至关重要。
第一性原理计算可以模拟催化剂表面的原子结构和可能的反应路径,并预测催化反应的活性和选择性。
这为催化剂的设计和优化提供了理论基础。
例如,在有机合成领域,通过计算催化剂表面和底物的相互作用能量,可以确定最佳的反应条件,提高催化剂的效率和产率。
第一性原理计算
材料是由大量的原子组成的多体体系,而原子又是由中子和质子所组成的原子核和核外电子所组成的。
材料的性能主要由核外电子之间的相互作用所决定。
原则上,如果可以写出构成材料的多体薛定愕方程,并求出该方程的解,就可得到材料的许多基本性质,如电导率、磁有序、振动谱、光学介电函数等。
但是,可以解析求解的系统仅限于氢原子,而由两个氢原子的氢分子和两个电子加两个质子组成的氦原子就己经无法求解了。
Hohenberg和Sham在1964年提出了一个重要的计算思想,证明了电子能量由电子密度决定。
所以就可以通过电子密度得到所有电子结构的信息而无需再处理复杂的多体电子波函数,只用三个空间变量就可描述电子结构,这种方法称为电子密度泛函理论。
按照该理论,粒子的哈密顿量由局域的电子密度决定,由此得到局域密度近似方法,基于该方法的自洽计算被称为第一性原理方法。
基于局域密度泛函的第一性原理方法对于电子基态的计算是非常准确的,与基态相关的电子能带结构、声子谱、结合能等都能用此种方法进行定量的计算。
第一性原理计算方法,例如密度泛函理论(DFT)计算,它将问题归结为对电子密度函数的描述,只需要将各类原子位置和个数作为参数输入计算即可。
它是一种预先定义的方法,它适用于周期表上的所有元素,而且大量的文献证明了它的准确可靠性。
和其他量子力学方法相同,第一性原理计算结果包含所有原子的位置,力场,电子结构(即“电子云”的描述),和体系的能量。
从第一性原理计算得到的基本结果,以及它们随时间演化的规律,我们能推出几乎材料所有的性质。
所以第一原理计算方法己经成为研究固体性质的一种重要的理论方法[20]。
费米能级是指费米子系统在趋于绝对零度时的化学位;但是在半导体物理和电子学领域中,费米能级则经常被当做电子或空穴化学势的代名词。
费米子可以是电子、质子、中子(自旋为半整数的粒子)。
晶体中电子所能具有的能量范围,在物理学中往往形象化地用一条条水平横线表示电子的各个能量值。
第一性原理计算方法讲义
第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中,求解的过程中需要知道一些物理参量,如温度场方程中的热传导系数和浓度场方程中的扩散系数等,这些参量随着材料的不同而改变,需要通过实验或经验来确定,所以这些方法也叫做经验或者半经验方法。
而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等,而不需要知道那些经验或半经验的参数。
第一性原理计算方法的理论基础是量子力学,即对体系薛定额方程的求解。
量子力学是反映微观粒子运动规律的理论。
量子力学的出现,使得人们对于物质微观结构的认识日益深入。
原则上,量子力学完全可以解释原子之间是如何相互作用从而构成固体的。
量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。
以量子力学为基础而发展起来的固体物理学,使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题,引发了通讯技术和计算机技术的重大变革。
目前,结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。
但是固体是具有~1023数量级粒子的多粒子系统,具体应用量子理论时会导致物理方程过于复杂以至于无法求解,所以将量子理论应用于固体系统必须采用一些近似和简化。
绝热近似(Born-Oppenheimei近似)将电子的运动和原子核的运动分开,从而将多粒子系统简化为多电子系统。
Hartree-Fock近似将多电子问题简化为仅与以单电子波函数(分子轨道)为基本变量的单粒子问题。
但是其中波函数的行列式表示使得求解需要非常大的计算量;对于研究分子体系,他可以作为一个很好的出发点,但是不适于研究固态体系。
1964年,Hohenberg和Kohn 提出了严格的密度泛函理论(Density Functional Theory, DFT)。
它建立在非均匀作为基本变量。
1965年,Kohn和Sham 电子气理论基础之上,以粒子数密度()r提出Kohn-Sham方程将复杂的多电子问题及其对应的薛定谔方程转化为相对简单的单电子问题及单电子Kohn-Sham方程。
