| 查看: 5060 | 回复: 24 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【原创】密度泛函理论与能带理论(一)
|
|||
|
1，密度泛函理论(DFT)在CASTEP中的应用可以这样简单地描述：用平面波赝势方法求解一组Schrodinger (Kohn-Sham)方程→波函数→电子密度→total energy and other properties 2，CASTEP最基本的计算结果就是total energy，许多其他的物理性质是以这一结果为基础得到的。例如，计算total energy对atomic positions的导数可以得到forces，计算total energy对lattice parameters的导数可以得到stresses 上面两点都是我们必须了解的基础知识。当然，要想更好地运用第一性原理，还必须进一步深入：“了解CASTEP是如何把能带理论融入密度泛函理论的”，要弄清楚能带结构(band structure)和态密度(DOS)在CASTEP中是如何计算得到的，它们是在DFT这一理论框架的什么位置切入的，是通过波函数的计算结果得到的还是通过电子密度的计算结果得到的（或者其他途径）？ 搞清楚这些问题是很有必要的，对于计算参数的选择以及计算结果的分析都至关重要！ 在此希望发挥集体的力量，希望各位虫友建言献策、提供文献，小弟愿效犬马之劳：为大家做一个清楚的总结！ |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有5人回复
-
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):各位专家都很厉害 12-7 15:04
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):各位专家都很厉害 12-7 15:04
|
根据波恩-奥本海默绝热近似，把离子实和价电子的运动分离开来；电子的哈密顿算符中存在电子之间的交换关联效应，无法分离变量求解，因此采用一种近似方法，考虑电子之间相互作用的平均效应，经过一些处理，得到K-S方程。由于处理电子相互作用时采用了一些近似，因此引入修正项Vex，而针对修正项则有很多近似，LDA和GGA等，LDA则有参数化的CA-PZ等方案，GGA包括PBE、RPBE、PW等。求解K-S方程主要是采用自洽场（SCF）方法，先初始一组波函数（例如原子轨道的线性组合），然后带入K-S方程组中，求出Vk-s，然后再求出一组新的波函数，再带入求出新的Vk-s，反复迭代，直到体系收敛（自洽），这叫做哈特利-福克自洽场方法。 在求解时，把原子看成是离子实和价电子，波函数用平面波基组展开，如果考虑原子的所有电子，则是全电子方法，而多数时候，体系的成键等性质只和外层电子有关，因此采用赝势法，即把原子势能的内层平滑化，截断半径以外则和全电子结构一致，这种简化减少了展开采用的平面波数目，从而简化计算，同时全电子和赝势得到的能量本征值相同，可以模拟体系的性质。赝势一般为模守恒势，但其较“硬”，从而耗时，于是又产生了超软势。 得到波函数后，用算符得到力学量的数值，然后求的体系的总能和其他性质。 个人理解，仅供参考！ |
14楼2009-12-07 12:47:29
2楼2009-12-07 01:09:24
3楼2009-12-07 08:19:07
4楼2009-12-07 08:48:49
20