密度泛函理论(DFT)(共4页)
4页1、精选优质文档-倾情为你奉上一、 计算方法密度泛函理论(DFT)、含时密度泛函理论(TDDFT)二、 计算方法原理1 计算方法出处及原理本计算方法设计来源于量子化学理论中的BornOppenheimer 近似,给近似下认为原子核不动, 这样电子就相当于在一个由核产生的外部的静态势场 V 中运动。那么一个固定的电子态可以用波函数 (, ,), 并且满足多 N 电子体系薛定谔方程: (2-3)其中,l , 哈密顿算符;l E, 体系总能量;l , 动能项;l , 由带正电的原子核引起的外场势能项;l , 电子电子相互作用能。通常把 和 叫做通用算符, 因为对于任何一个 N 电子体系, 表达式都相同.而势能函数 与体系密切相关。由于电子相互作用项 的存在, 复杂的多体系的薛定谔方程公式 2-3并不能拆分为简单的单电子体系的薛定谔方程。根据 DFT 的核心理念, 对于一个归一化的波函数 , 电子的密度 n() 可以定义为: (2-4)更重要的是, DFT 的核心理念告诉我们, 对于一个给定的基态, 如果基态的电子密度是知道的话, 那么基态的波函数就唯一确定。也就是说, 基态的波函数是基态电子密度
2、的泛函11, 表达为: (2-5)既然有以上的假定, 那么对于基态的任何一个观测量, 它的数学期望就应该是的泛函: (2-6)特别的, 基态的能量也是的泛函: (2-7)这里外部势能的贡献可以通过基态的电子密度来精确表达: (2-8)或者外部势能可以用电子密度 n 来表达: (2-9)泛函 T n 和 U n 被称作通用泛函, 而势能泛函 V n 被称做非通用泛函, 因为它与当前研究的系统息息相关。对于一个给定的体系, 就存在一个对应的,相应的, 该体系的能量可以表达为: (2-10)假定, 已经得到了T n 和 U n 的表达式, 那么对于公式 2-10, 以 为自变量, 求解 E n 的最小值, 就可以得到基态的对应的能量 E0 , 同样也能得到其他的基态的客观测量。求解能量最小值的变分问题可以通过 Lagrangian 乘数待定法 32 来轻松解决12。首先, 假定, 不考虑电子电子相互作用的体系, 能量可以表达为: (2-11)其中, 是不包含电子电子相互作用的体系动能项, 是不包含电子电子相互作用情况下的电子所处的外部有效势能。很明显, 如果我们将表达为: (2-12)那么可
3、以把不考虑电子相互作用情况下的电子密度定义为: (2-13)这样我们就得到一个不含电子电子相互作用体系的所谓的 KohnSham 方程: (2-14)通过该式公式 2-14可以得到分子轨道, 得到分子轨道之后, 当然可以得到原来的包含电子电子相互作用体系的电子密度, 如下: (2-15)这时, 可以把有效单粒子的势能精确地表达为: (2-16)上式的第二项通常被称作 Hartree 项, 描述的电子与电子之间的库仑斥力作用。最后一项,描述的是电子交换相关势能 (exchangecorrelation potential)。在公式 2-16中, 包含多体体系中的所有的相互作用。由于 Hartree 项, 项都是的函数; 而电子密度又是波函数的函数, 同时波函数反过来又是的函数。这样, 求解 KohnSham 方程公式 2-14就成了一个自洽的过程。落实到量子化学中的具体计算中, 就是先猜测一个初始的电子密度, 然后计算对应的并求解 KohnSham 方程公式 2-14得到波函数。既然有了波函数, 反过来就有了此波函数对应的电子密度, 可以用这个新得到的电子密度, 然后再去求解新的波函数, 以及电子密度。什么时候达到所谓的收敛呢? 就是你当前循环猜测的和基于此猜测值通过 KohnSham 方程公式 2-14 求解出来的波函数所确定的电子密度一致, 就是所谓的收敛。2 计算方法应用领域此方法多用于材料合成领域前期材料性能预测,以及后期材料性能分析。专心-专注-专业
《密度泛函理论(DFT)(共4页)》由会员M****1分享,可在线阅读,更多相关《密度泛函理论(DFT)(共4页)》请在金锄头文库上搜索。
2022年公司财务年终工作总结范文
民用建筑电气毕业设计论文
达成遂成遂渝线隧道应急抢险施工组织
2020年小学学年第一学期德育工作总结..doc
2023出纳辞职报告范文15篇
山东省高院审理人身损害纠纷
2014-2015学年度第一学期高一数学考试卷6
传染病防治应急预案标准版本(2篇)
物资材料采购计划表
职业学校工会半年工作总结(三篇).doc
笑猫日记读后感15篇7
钢结构课程设计报告.doc
古典诗词写景七法
2023医院工作总结样本(7篇).doc
南阳收获机械项目投资计划书_范文模板
洁净厂房管理规程
物料性质及安全
北京关于成立驱动系统产品公司可行性报告
金融统计与分析相关资料
2015年高中英语 VOA习惯用语 第259讲 pinhead素材
2023-07-22 5页
2023-03-10 5页
2023-07-25 9页
2023-04-18 6页
2022-11-07 31页
2022-10-26 12页
2023-05-11 14页
2023-01-02 58页
2023-05-27 4页
2023-05-12 7页