密度泛函理论如何在企业落地

密度泛函理论如何在企业落地
密度泛函理论及其应用一、密度泛函理论（Density Functional Theory ：DFT ）
VASP 的理论基础是电荷密度泛函理论在局域电荷密度近似（LDA ）或是广义梯度近似（GGA ）的版本。DFT 所描述的电子气体交互作用被认为是对大部分的状况都是够精确的，并且它是唯一能实际有效分析周期性系统的理论方法。
11 单电子薛定谔方程式
一个稳定态（与时间无关）的单一粒子薛定谔方程式可表示为一个本征值问题（暂略动能项的 ）： /2m ()()H r E r ψψ=（1）
2[]()()V r E r ψψ-∇+=（2）多体量子系统 （如双电子的薛定谔方程式）： 2212121212[(,)](,)(,)V r r r r E r r ψψ-∇-∇+=（3）在普遍的状况下，里的是无法分离变量的，因此，即便简单如12(,)V r r 12,r r 双电子的薛定谔方程式就己经没有解析解了。而任何的计算材料的量子力学问题，都需要处理大量数目的电子。 12 Hohenberg-Kohn 定理量子力学作为20世纪最伟大的发现之一，是整个现代物理学的基石。量子力学最流行的表述形式是薛定谔的波动力学形式，它的核心是波函数及其运动方程薛定谔方程。对一个给定的系统，我们可能得到的所有信息都包含在系统
的波函数当中。对一个外势场v (r)中的N 电子体系，量子力学的波动力学范式可
以表示成：
v (r) Ψ (r1; r2; …; r N ) 可观测量
⇒⇒（4）
即，对给定的外势，将其代入薛定谔方程可以得到电子波函数，进一步通
过波函数计算力学量算符的期望值可以得到所有可观测量的值。电荷密度是这
些可观测量中的一个：
333232()(,)N N n r N d r d r d r r r r =ψ⎰⎰⎰ 2(,)N r r r ψ （5）
如前所述，任何的计算材料的量子力学问题，都需要处理大量数目的电子。
而，对于超过两个电子以上的体系，薛定谔方程就已经难以严格求解了。对于
实际物质的这样一种每立方米中有数量级的原子核和电子的多粒子系统，
2910我们是更不可能由薛定谔方程来严格求解其体系的电子结构的。但，建立于
Hohenberg-Kohn 定理上的密度泛函理论不但给出了将多电子问题简化为单电子
问题的理论基础，同时也成为分子和固体的电子结构和总能量计算的有力工具。
因此，密度泛函理论是多粒子系统理论基态研究的重要方法。
密度泛函理论的基本想法是原子、分子和固体的基态物理性质可以用粒子密度函数来描述，这源于HThomas 和E·费米1927年的工作。密度泛函理论基础
是建立在PHohenberg 和WKohn 的关于非均匀电子气理论基础上的，它可归结
为两个基本定理：
定理一：不计自旋的全同费米子系统的基态能量是粒子数密度函数的()n r 唯一泛函。
它的推论是，任何一个多电子体系的基态总能量都是电荷密度的唯一()n r
泛函，唯一确定了体系的（非简并）基态性质。
()n r
由于电荷密度与电子数直接联系:，这样决定多电子薛定谔
N ()n r dr N =⎰
方程解的电子数和外势场都由电荷密度唯一确定，因此基态波函数
N ()n r
以及其它的电子结构性质都由电荷密度唯一确定。
[]F n 由于决定了哈密顿量，多电子体系的基态是的唯一泛函，自然
()V r ψ()n r
动能和库仑能也是的泛函，那么体系的所有性质也将是基态密度的泛函。
()n r
于是定义一个普适泛函，有：
[]F n 2,,2
2,()1
(1)()2()l ps l ps l l ps d r l l V r E r dr r ⎡⎤Φ+⎡⎤
=
+-⎢⎥⎢⎥Φ⎣
⎦⎢⎥⎣⎦ []ˆˆ()F n r T
U ≡<ψ+ψ> （6）
适用于任何外场下的具有任意电子数的体系。所以系统基态的能量可表示
为泛函的形式：
[]()()()()()()ext ee HK E n T n V n V n n r V r dr F n =++=+⎰
（7）
这里:
（8）
[][][]HK ee F n T n V n =+非经典项
[][]ee V n J n =+（9）
[]()()121212
112J n n r n r dr dr r =
⎰
（10）
其中是经典电子排斥能。非经典项是一个难以理解而又非常重要的量，
[]J n 它是交换-关联能的主要来源。
定理二：能量泛函在粒子数不变条件下对正确的粒子数密度函数
[]E n 取极小值，并等于基态能量。
()n r
在电子数恒定的约束条件下：，按照Hohenberg-Kohn 第二定理，
()n r dr N =⎰
基态能量满足如下条件：，即
[]{
}
()0E n n r dr N δμ⎡⎤--=⎣⎦
⎰
[][][]{}()()
ee ext
T n V n E n V n r n r δδμδδ+==+ （11）
因而只要知道和的泛函形式，就可以通过上式来求解电子
[]()T n r []()ee V n r
结构。
利用上述性质，我们会想利用各种方法猜测并代入 E GS[] 求值，
()n r ()n r
只要一直试到产生最低的能量，则该能量保证是基态的总能，且该电荷分布保
证是基态的电荷密度分布。
Hohenberg-Kohn 的密度泛函理论（DFT ）只有对基态才是严格成立的。但，即使只是获得基态，都已经足以预测很多性质了。例如，分子的键长，振动频
率，固体的晶胞边长、d性系数张量，甚至是化学键的断裂或是生成，对电子
而言都是基态的性质。因此，能预测系统的基态是非常有用的。
13 Kohn-Sham 方程
有了上述两个定理，剩下的问题就是能量泛函的具体表述形式。Kohn 等人
引进了一个与相互作用多电子体系有相同电子密度的假想的非相互作用多电子
体系。因为电子密度一般可以表示成轨道形式，这个假想的非相互作用体系的
动能算符期望值可以非常简单的写成各电子动能的和。
[]232
()()2N S i i i
T n d r r r m φφ=-∇∑⎰ （
1
2）
其中是密度函数对应的Khon-Sham(KS)轨道。将的主要部分写成：
()i r φ
ee V
[]233()()
2ee eeH q n r n r V V n d r d r r r '≈='
-⎰⎰
（13）
至此，我们得到一个很自然的关于能量泛函中未知项（交换相关泛函）的
定义： ()()xc tot s ext eeH S ee ee E E T V V T T V V =---=-+-（14）
将能量泛函对KS 轨道进行变分可以得到著名的KS 方程：
21(()()())2
ext H xc i i i v r v r v r ε-∇+++Φ=Φ
（15）
其中、、分别是外势、Hartree 势和交换相关势。在KS
()ext v r ()H v r ()xc v r
方程中，有效势由电子密度决定，而电子密度又由方程的本征
eff ext H xc v v v v =++函数—KS 轨道求得所以我们需要自洽求解KS 方程。这种自洽求解过程通常被
称为自洽场（SCF ）方法
当我们得到一个自洽收敛的电荷密度后，我们就可以得到系统的总能
0n
[]233
300000()()()()2N
i xc xc i
n r n r q E d r d r d rv r n r E n r r ε''=--+'-∑⎰⎰⎰
（16）
其中是KS 方程的本征值。
i εKohn-Sam 能量泛函使我们有可能通过近似方法来描述与电子密度有关的交
换关联能。而密度泛函理论的发展就是以寻找合适的交换相关泛函近似形式为
主线的。人们发展了局域密度近似（LDA ）、局域自旋密度近似（LSDA ）、
广义梯度近似、X3LYP 等等，并取得了很好的结果。
14 赝势
赝势是一个用来模拟离子实对价电子作用的有效势。其物理本质在于价态
芯态正交条件对价态的贡献等效于一个排斥势，它与芯区的势对价电子的强烈
吸引相互抵消，使得构造一个相对平缓的有效赝势成为可能。赝势方法的发展
经历了从经验赝势、模守恒赝势到超软赝势的几个阶段。
VASP 中使用的是超软赝势与缀加平面波势。
使用赝势可以帮助我们方便的处理电子-离子间的
交互作用。那么，赝势是如何产生的呢？ 如图10示，实线分别是真实位势与 All-
/Z r electron 价电子波函数，我们要取距原子中心
V ψ处为划分点，以上波函数完全一样保留，而以
c r c r c r 内则对波函数加以改造。主要是要把振荡剧烈的波函
数改造为一种变化缓慢的波函数，而它须要是没有节点

的，如虚线的所示。少了剧烈振荡不但允许只以相对很少的平面波来展
pseuto ψ开波函数，没有节点的（径向）波函数也意味着没有比它本征值更低的量子态
来与它正交。求解内层电子的需要就自动消失了。
￥
5
百度文库VIP限时优惠现在开通,立享6亿+VIP内容
立即获取
密度泛函理论及其应用
密度泛函理论及其应用一、密度泛函理论（Density Functional Theory ：DFT ）
VASP 的理论基础是电荷密度泛函理论在局域电荷密度近似（LDA ）或是广义梯度近似（GGA ）的版本。DFT 所描述的电子气体交互作用被认为是对大部分的状况都是够精确的，并且它是唯一能实际有效分析周期性系统的理论方法。
11 单电子薛定谔方程式
一个稳定态（与时间无关）的单一粒子薛定谔方程式可表示为一个本征值问题（暂略动能项的 ）： /2m ()()H r E r ψψ=（1）
第 1 页
2[]()()V r E r ψψ-∇+=（2）多体量子系统 （如双电子的薛定谔方程式）： 2212121212[(,)](,)(,)V r r r r E r r ψψ-∇-∇+=（3）在普遍的状况下，里的是无法分离变量的，因此，即便简单如12(,)V r r 12,r r 双电子的薛定谔方程式就己经没有解析解了。而任何的计算材料的量子力学问题，都需要处理大量数目的电子。 12 Hohenberg-Kohn 定理量子力学作为20世纪最伟大的发现之一，是整个现代物理学的基石。量子力学最流行的表述形式是薛定谔的波动力学形式，它的核心是波函数及其运动方程薛定谔方程。对一个给定的系统，我们可能得到的所有信息都包含在系统

dft密度泛函理论用于计算原子和分子的电子结构。

密度泛函理论(Density functional theory ，缩写DFT)是一种研究多电子体系电子结构的方法。密度泛函理论在物理和化学上都有广泛的应用，特别是用来研究分子和凝聚态的性质，是凝聚态物理计算材料学和计算化学领域最常用的方法之一。

虽然密度泛函理论的概念起源于Thomas-Fermi模型，但直到Hohenberg-Kohn定理提出之后才有了坚实的理论依据。Hohenberg-Kohn第一定理指出体系的基态能量仅仅是电子密度的泛函。

理论应用：

自1970年以来，密度泛函理论在固体物理学的计算中得到广泛的应用。在多数情况下，与其他解决量子力学多体问题的方法相比，采用局域密度近似的密度泛函理论给出了非常令人满意的结果。直到二十世纪九十年代，理论中所采用的近似被重新提炼成更好的交换相关作用模型。密度泛函理论是目前多种领域中电子结构计算的领先方法。

尽管密度泛函理论得到了改进，但是用它来恰当的描述分子间相互作用，特别是范德瓦尔斯力，或者计算半导体的能隙还是有一定困难的。

我都没听过什么叫“拉曼光谱”“密度泛函基组”

网上流传的答案是没有“直接了当”的解答，还是牛头不对马嘴？

题主是从哪里得到这个计算任务的？有没有例子？有没有老师？有没有同学？

---