DFT电子结构计算李震宇

《DFT电子结构计算李震宇》由会员分享，可在线阅读，更多相关《DFT电子结构计算李震宇（51页珍藏版）》请在金锄头文库上搜索。

1、 电子结构理论与计算 李震宇 中国科学技术大学合肥微尺度物质科学国家实验室 Version: 20110812 1主要内容 主要内容 第一章第一章 电子结构理论基础电子结构理论基础 BO 近似， HF 理论， 电子关联与后 HF 理论， 密度泛函理论 （DFT） ，交换关联泛函，赝势，半经验电子结构理论 第二章第二章 电子结构计算方法电子结构计算方法 原子基组，平面波基组与周期性边界条件，自洽场（SCF） 第三章第三章 电子结构分析电子结构分析 电荷密度，分子轨道，能量 第四章第四章 电子结构理论与计算的应用电子结构理论与计算的应用 动力学与统计力学性质 附录附录 量子力学基本原理，常见数值算

2、法 2前言 前言 理念：1. 实用：介绍计算参数的意义和选择策略；2. 简洁：注重物理概念，不过分纠缠于细节。 参考资料： 1 Richard M Martin, Electronic Structure: Basic Theory and Practical Methods. Cambridge University Press, Cambridge, 2004. 2 Christopher J Cramer, Essentials of Computational Chemistry: Theories and Models. 2nd Ed. Wiley, England, 2004. 感

3、谢所有帮助我将本讲义电子化的同学，以及所有其他选修本课程的同学。讲义中错误和疏漏之处在所难免，如有任何意见与建议，请发送到，十分感谢！ 3第一章 电子结构理论基础 第一章 电子结构理论基础 由于计算能力较弱， 早期的科学研究以理论和实验的互动为主， 其中理论研究通常使用 唯象理论。 随着科技的发展， 现在已经有可能完全从量子力学基本原理出发来计算体系的各 种性质。 这就使得计算成为一种重要的独立的科学研究手段， 与理论和实验形成三足鼎立之 势。 图示：新的科学研究范式 三者之间有紧密的联系：计算可以帮助理解实验结果，例如扫描隧道显微镜(STM)观测到的 是局域态密度，要从中获得几何结构的信息，

4、需要计算模拟来配合；计算也可以验证理论猜 想， 例如通过计算机实验来检验材料设计原理。 我们希望通过科学计算获得理论设计与计算 表征的能力。 体系的大多数性质都决定于它的电子结构，即电子的空间与能量分布。 1.1 基于波函数的量子化学方法基于波函数的量子化学方法 1.1.1 波恩波恩-奥本海默奥本海默(BO)近似近似 在原子中，核的质量比电子大很多，运动也比电子慢得多。1 amu (12C/12)是电子静止 质量的 1822.83 倍。因此，通常在研究某一瞬间电子结构时可以忽略原子核的速度，假定 它和原子核长期固定在某位置时的电子结构一样。这就是 BO 近似。其核心思想是把核的运 动和电子的运

5、动分开处理：处理电子运动时，认为核是固定不动的；处理核运动时，认为快 速运动的电子建立一个平均化了的负电荷分布，核在电子的负电荷平均场中运动。 对于包含许多电子和原子核的体系，其定态薛定谔方程可写为： ,HE r Rr R 其中 22222 22 22IJI Ii IiiIijIJIiiIijIJe Z Ze ZeHmmrrr- 上述哈密顿量中第一项核动能项为小项。引入 BO 近似，对原子核和电子进行分离变量 ,;Nelr RRr R elelelHE RrRr 其中 NR是描写原子核状态的波函数，它只和所有核的位置有关；而 elH R和4;elr R是描写电子状态的哈密顿量和波函数，以原子核

6、位置为参数【输入参数：POSCAR】 。 E R和原子核的位置有关，通常被称为势能面(PES)。知道势能面以后可以从其极小点得到几何构型，从其势阱深度得到反应能，从其鞍点得到过渡态，等等。 BO 近似又称为“绝热近似” ，它认为电子总是不断地适应核位置的变化，或者说跟随 着核且始终围绕着核来运动， 运动过程中不会发生从一个势能面 E1至另一个势能面 E2的跃 迁，即吸收或放出能量发生电子能级跃迁。这是因为电子态涉及到的能量量级为电子伏特， 远大于室温对应的能量，所以电子通常处于基态。 图示：石墨烯能量随 C-C 键长变化的势能面 1.1.2 单电子近似与平均场近似单电子近似与平均场近似 由于电

7、子之间存在相互作用，所以当电子数较多时，即使在 BO 近似下薛定谔方程仍然 很难求解，需要作进一步的简化。最粗糙的近似是假设电子之间没有任何相互作用。这样， 体系的哈密顿算符简化为： i iHh 其中 222I ii IiiIe Zhmr 于是，我们可以先只考虑体系中的某一个电子，其定态薛定谔方程为iiiih ，本征函数i为单电子波函数，通常也称为分子轨道(MO)，本征值i为对应的分子能级。下面的问题是怎么从单电子波函数写出总的多体波函数呢？ 5对无相互作用电子，考虑泡利不相容原理，体系的波函数及其对应的能量本征值为 12.HPn i iE 这种由分子轨道直接相乘得到的多体波函数称为 Hart

8、ree 积(HP)。 单电子近似并不意味着需要完全忽略电子间的相互作用。 对严格的 Hamiltonian 采用 HP多体波函数作为试探波函数能量(elH在多体波函数下的期望值)变分可以得到新的单电子哈密顿算符 22 2 2I iii IiiIe Zhvmr 其中 j ij j iijvdrr 2jje 显然， 该单电子哈密顿算符中的第三项可以看成体系中其他电子的库仑相互作用的一个平均 场近似，该平均场不随时间而改变。因为在求平均场时其他电子的坐标都是积分变量，所以 ih仅与第 i 个电子的坐标有关。对这种采用 HP 多体波函数的平均场近似，体系的总能量是 222| |1 2ijHP iij

9、iijijeEd drr r 其中第二项为重复计数修正。 我们看到，在平均场近似中，需要求解单电子本征方程才能得到体系的解，但要求解单 电子方程又必须知道各个轨道的电荷分布。 这看起来自相矛盾。 在实际操作中我们可以从 N个初始单粒子波函数 0 i出发，建立相应的(0) ih，然后通过解单电子本征方程，得到一组新的分子轨道 1 i，再用这组 1 i来建立新的(1) ih，重复这个过程，直到最后得到的 n i与上一步的1n i相等或近似相等。这样一个过程称为自洽场(SCF)方法【输入参数：EDIFF】 。 1.1.3 Slater 行列式与行列式与 Hatree-Fock 自洽场自洽场 显然，H

10、P 多体波函数不满足电子（费米子）全同性所要求的波函数的反对称性。为简 单起见，先讨论无相互作用两电子体系 12Hhh 其中 2 2 1112hvm r6 2 2 2222hvm r假设 ,12ababE 是两电子体系本征方程的解，则 12122112abababP 也是体系本征方程的解。这样，对应于同一个本征值，有两个线性独立的本征函数，但这两 个本征函数都不满足反对称条件，为此作一线性组合，得： 112212Aabab为体系满足波函数反对称条件的解，用行列式的形式表示为 111 222ab SD ab该行列式被称为 Slater 行列式(SD)。对于由 N 个电子组成的全同粒子体系，可以将

11、上述结 果进行直接推广，相应的定态波函数为 1212 121212111 2221,!NN SDNNNN NNN r rr 上式的反对称性是很明显的， 因为任意两个粒子坐标的互换， 相当于行列式中相应的两 列元素的互换，将改变行列式的符号。同时，若有两个或两个以上单粒子波函数相等，则0SD，意味着这样的状态是不存在的，这就是 Pauli 不相容原理。 Fock 建议将 Hartree-SCF 过程推广到采用 Slater 行列式作为多体波函数的情形。对两 电子体系的电子库仑相互作用算符，我们有【作业】 12 121SDSDababd dJKrr r 其中 22 1212 12121112abH

12、PHPabJd dd drrr rr r 12 1211122abababKd drr r 分别代表电子间的库仑排斥能和交换能。交换能 Kab是采用 Slater 行列式以后比 HP 多出来 的同自旋电子间的一种相互作用。 Roothann 提出可以基组表象来建立了 HFSCF 的矩阵代数。 如果将分子轨道写成一些基函数的线性组合， 则单电子哈密顿量在这组基组下可以表示为一 个矩阵。常用的基组有原子轨道和平面波。对闭壳层体系，最后可以得到久期方程det(F-ES)=0，其中单电子哈密顿量f 通常被称为 Fock 算符，其矩阵元为72 211 22I IeIFZPmr 其中 212 121122

13、ed dr r r 2occii iPa a由于 Fock 矩阵的表达式中包含密度矩阵，我们需要自洽求解。 至此，我们通过如下步骤得到了作为量子化学基础的 HartreeFork 理论。 1.1.4 电子关联与电子关联与 post-HF 方法方法 HF 理论采用单行列式多体波函数，它假定一个电子在由原子核和其它电子形成的平均 势场中独立运动，没有考虑电子间的瞬时相关作用，也就是电子关联。但是，HF 理论中采 用的 Slater 行列式波函数满足 Pauli 不相容原理，所以它在一定程度上包含了同自旋电子的 关联作用，这部分电子关联被单独命名为交换相互作用。所以，通常所说电子关联能是指精 确的基

14、态能量与 HF 能量之差：Ecorr=E-EHF。 关联能一般只占体系电子总能量的 0.32%，但是它可以与化学反应的反应热或者活化 能相当甚至更大。一般说来，化学反应途径上电子关联程度不同。所以，在考虑化学反应和 电子激发等问题时电子关联十分重要。 由电子瞬时作用产生的电子关联称为动态关联。 此外， 还有非动态关联，它指体系几个对称性相同的电子组态函数（单行列式）接近简并时导致 HF 方法给出很大的误差。常用的电子关联计算方法【输入参数：计算方法】有 组态相互作用（CI）：Configuration Interaction 耦合簇方法（CC）：Coupledcluster method 多体微扰理论（MPx）：MollerPlesset perturbation theory 密度泛函理论（DFT）：Density function theory 前面三种都是以 HF 的为基础，经常被称为 postHF 方法。而 DFT 基于完全不同的思路，将 在下一节中专门介绍。 1. CI 理论 可以通过引入激发组态来考虑电子关联 800|.vir occvir occ aaabab iiijij aia b ijCCC 一般用 S, D, T, Q 等字母来表示从低到高的激发组态，其中最简单的情况是只引入单激发组 态(CIS)。CIS 可用来计算激发能，但对基态能量没有修正【作业】。最复