| 查看: 5122 | 回复: 16 | ||||||
| 【奖励】 本帖被评价11次，作者luzihen增加金币 8.8 个 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
[资源]
何谓Pseudopotential,APW,PAW方法
|
||||||
|
心烦意乱.感觉这个博士读的真是....动心忍性... 闲话,下面正题,记录一下看资料后整理出的思路.如果有错请指出. ================================================================== DFT计算的时候需要求解K-S方程. 而在求解K-S方程的时候,我们并不是直接求解. 我们把方程投影到不同的基组上进行求解,比如VASP里面常用的平面波基组,再比如其他软件的gaussian基组. 在用平面波求解问题的时候,一个显著的问题是,近核区域的波函数起伏很大(对应1s,2s等轨道),需要很多组平面波才能表示物理的波函数(近核电子波函数).如此,求解方程会非常慢,需要非常大的内存. 这个问题目前主要有两种解决方法 1. Pseudopotential:在K-S方程中势能V上动手脚 2. Augmented Plane Wave:在基组上动手脚 3. Projector Augmented Plane Wave 第一种Pseudopotential不说了,就是近核区域的-q/r换掉,换成Pseudopotential 第二种是把Augmented Plane Wave, 平面波对近核区域表现不好,一些局域基组对近核区域表现比较好,于是两者组合一下,近核区域用局域基组,之外的区域用平面波基组.如此达到全电子计算. 在处理两个基组耦合问题时,有了Linearized Augmented Plane Wave method,LAPW有了和Pseudopotential差不多的效率 第三种Projector Augmented Plane Wave,文献中说和第二种差不多,我的理解是:如上所述,求解K-S方程时候进行了一次投影,即将波函数投影到平面波基组上求解。这里为了减少平面波数量，再这之间加一层投影。这里中间夹的这一层投影算子可以把物理的波函数（近核区域起伏很大，有零区域）投影成一个比较平滑，适合平面波基组求解的波函数（减少平面波数量）。利用这个夹在中间的投影算子可以来回投影平滑的波函数和物理的波函数，所以需要的物理量可以得出。 注：上面说的投影都是展开之意，有时候基组并不完全，可以利用最小二程确定系数。 [1] PRB上PAW方法的原始文献不容易看懂，附上一篇BLÖCHL的比较容易看懂的文章。其中nodal areas意思为areas where wave-function equal 0. Blöchl, Peter E., Clemens J. Först, and Johannes Schimpl. "Projector augmented wave method: ab initio molecular dynamics with full wave functions." Bulletin of Materials Science 26.1 (2003): 33-41. [2] PRB的原始文献： Blöchl, Peter E. "Projector augmented-wave method." Physical Review B 50.24 (1994): 17953. |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : PAW.pdf
2015-08-07 22:13:08, 169.75 K
» 收录本帖的淘帖专辑推荐
 VASP | 第一性原理计算-tg | 计算 | soft study |
» 猜你喜欢
重要来源：本周末出结果
已经有3人回复
-
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
7楼2015-08-10 06:47:23
2楼2015-08-07 22:24:35
3楼2015-08-08 17:21:22
4楼2015-08-08 19:35:11
简单回复
Quan.6楼
2015-08-09 06:06
回复
五星好评 顶一下，感谢分享！
seer9楼
2016-09-12 21:07
回复
五星好评 顶一下，感谢分享！
awmc20085楼
2015-08-08 22:00
回复
五星好评 顶一下，感谢分享！
墨兰清12楼
2016-11-28 15:42
回复
五星好评 顶一下，感谢分享！
wurenglai13楼
2016-11-30 09:53
回复
五星好评 顶一下，感谢分享！
wurenglai14楼
2016-11-30 09:54
回复
顶一下，感谢分享！
sixiangzhe1411楼
2016-10-26 20:03
回复
五星好评 顶一下，感谢分享！
q96524001010楼
2016-09-24 19:17
回复
五星好评 顶一下，感谢分享！
剑雨孤星8楼
2016-03-24 17:49
回复
五星好评 顶一下，感谢分享！
