Hohenberg-Kohn 定理:Hohenberg-Kohn定理指出,多粒子系统的所有基态性质都是密度泛函的唯一泛函,即总能可以通过求解最有利的基态电子密度分布函数得到。因此系统的总能量可表达为:T[n]代表无相互作用的粒子的动能,U[n]代表一般的库仑能,Exc[n]代表多体效应的交换-关联能之和。Hohenberg-Kohn定理还未能解决,即如何确定粒子数密度函数ρ(r),动能泛函T[ρ]和交换关联能泛函[ρ]xcE 。
Kohn-Sham 方程:利用变分原理,将能量密度针对电子密度进行变分,可以得到Kohn-Sham方程。Kohn-Sham 近似的核心思想:
(1) 动能的大部分通过相同电子密度的无相互作用体系来计算
(2) 电子相互作用中库仑作用占据了主要部分,而交换相关是相对次要的
(3) 非经典的交换和相关作用,动能校正项,自相互作用折入交换相关泛中
交换关联泛函:在求解Kohn-Sham 方程时必须要知道Vxc(r) 的准确解和便于表达的解析形式才能求解该方程。由于交换相关项包含许多非经典项,至今仍没有准确的函数描述。一般把交换相关项分为两个部分,即交换部分Ex 和相关部分Ec。 此外,对于动能的近似也被归并到交换相关项中。交换项和相关项的比重分别为90%和10%,即交换项起着更重要的作用。交换相关泛函包括局域密度近似泛函,广义梯度近似泛函和杂化泛函。
局域密度近似(LDA)泛函:即泛函在空间上的能量密度只和每一点上的电子密度有关。该类泛函常用的有SVMN和SVWN5等。对于电子密度ρ(r)空间变化不大的体系而言,采用均匀电子气模型的局域密度近似(LDA)是一种很好的近似。但是对于电子
密度相对空间变化较大的电子气体系,应用局域密度理论去完全消除电子与其自身的相互作用时的失效而引起的这些误差被证明是比较大的。
广义梯度近似(GGA)泛函:泛函在空间上的能量密度不仅与每一点上的电子密度有关,而且还与电子密度的梯度有关。常用的GGA泛函有BLYP、PBE等。除了密度梯度以外,在泛函中包含其他自变量(如动能密度或是局域自旋密度的
二阶导数)而得到的泛函称为基于高阶广义梯度近似的meta-GGA泛函,如M06-L、TPSS 等。而泛函中包含的信息越多,就越能正确地描述客观体系。
杂化泛函: 用HF方法中交换能和DFT方法中交换能做线性组合得到计算体系的交换关联泛函。杂化泛函分为单杂化泛函和双杂化泛函。
单杂化泛函:单杂化泛函大致分为全局杂化泛函和局域杂化泛函两类。
全局杂化泛函:以B3LYP为主要代表,其核心思想是在交换泛函部分引入了Hertree-Fock形式的交换能,在分子体系中最为常用的泛函有B3LYP、PBE0、X3LYP、TPSSh、M06等。
局域杂化泛函:局域杂化泛函是利用分区方案,将电子-电子间的相互作用部分分为短程与长程部分分别处理。分子体系中著名的有wB97X-D泛函,可以均衡的处理价键和能垒等。固体中最成功的为HSE泛函,可以更好地计算固体的带隙。
双杂化泛函:这一类泛函不仅要包含占据轨道的信息,而且要能恰当地考虑未占据轨道的信息。其核心思想是在交换部分引入部分HF形式交换能的同时,在相干泛函部分引入部分二级微扰形式MP2相关能,使得此类方法描述体系性质的误差有所降低。这一类的代表泛函有B2PLYP,mPW2PLYP和MC3BB等。
基于高阶广义梯度近似的meta-GGA泛函:MO6-L
全局杂化泛函:B3LYP
笔记
有CP2K使用常识的人都知道,在做杂化泛函计算之前应当先做纯泛函计算,以其收敛的波函数作为杂化泛函计算的初猜波函数。额外做的纯泛函计算耗时只占杂化泛函计算的一个零头,却能令杂化泛函容易收敛得多得多得多,并因此巨幅节约耗时。为实现这个目的,先用Multiwfn产生一个PBE泛函的输入文件PBE.inp并进行计算(基组和PBE0计算时相同),算完后会得到PBE-RESTART.wfn文件,然后在PBE0的输入文件里把# WFN_RESTART_FILE_NAME to_be_specified这一行改为WFN_RESTART_FILE_NAME PBE-RESTART.wfn,然后把# SCF_GUESS RESTART这一行的#去掉。这样PBE0计算开始时就会把PBE-RESTART.wfn里的PBE波函数当初猜了。
每天读点密度泛函理论 (David S. Sholl著,李健译)
微量杂质引发的金属脆化
相比于纯金属所具有的晶体结构,晶粒之间的接触区域存在有大量复杂结构。这些接触区域存在于所有多晶体材料中,称为晶界。
早在100年前,人们就已经知道:在铜中加入少量的某种杂质,可以将其从塑性状态变成脆性断裂方式(断裂之前没有塑性变形)。例如,当在铜(Cu)中铋(Bi)的含量低于100ug/g时,铜就会发生脆化。对于铅(Pb)和汞(Hg)杂质,也观察到类似的效果。
当杂质元素Bi的含量很低时就会发生脆化现象,这并不是什么难以置信的事情,其原因在于Bi几乎不溶于体相Cu中。也就是说,Bi原子更有可能偏聚在晶界处,而不是存在于晶粒内部,因此晶界上Bi的浓度要比该材料整体的平均浓度高得多。
研究者们已经给出了几种不同的解释。其中:第一类解释把发生这一现象的原因归结于电荷作用。例如,Bi原子与邻近Cu原子形成键合作用,使这些Cu原子要比纯Cu更硬,从而降低了Cu晶格平滑变形的能力。第二类解释是紧邻晶界的杂质原子可能改变了晶界处原子的电子结构,从而弱化晶界处原本的成键。第三种解释则归因于尺寸效应,认为Bi原子比Cu原子大得多,如果Bi原子存在于晶界处,则会分离界面两侧的Cu原子,使二者之间的自然间距增大。被拉大的键长将会弱化这些原子键,并使晶界处更容易断裂。第二和第三种解释都涉及到了晶界附近原子键的弱化,但对弱化给出了不同的根本原因。直接采用实验方法区分和辨别这些可能的原因机理是非常困难的。
范德华半径:若原子之间是通过范德华力而产生相互作用,且没有形成化学键,则相邻两个原子核之间距离的1/2为范德华半径。