24小时热门版块排行榜    

查看: 3957  |  回复: 16
【奖励】 本帖被评价11次,作者luzihen增加金币 8.8

luzihen

金虫 (正式写手)


[资源] 何谓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.
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:libolin3@tal.com
  • 附件 1 : PAW.pdf
  • 2015-08-07 22:13:08, 169.75 K

» 收录本帖的淘贴专辑推荐

VASP 第一性原理计算-tg 计算 soft study

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)


对了,忘了最初的问题。

那么如此说来PAW方法不能称为Pseudopotential,就是文章里这样写不好?
"PAW-PBE pseudopotentials were employed."
2楼2015-08-07 22:24:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)


我的理解PAW是pseudopotential的一种方法,属于plane wave basis set
PBE是一种GGA exchange correlation泛函
3楼2015-08-08 17:21:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)


引用回帖:
3楼: Originally posted by awdawdawdd at 2015-08-08 17:21:22
我的理解PAW是pseudopotential的一种方法,属于plane wave basis set
PBE是一种GGA exchange correlation泛函

问题就是paw似乎并不是淹势。
不过很多文献里面就这么写了,包括一些prl,所以应该不是很要要紧。

[ 发自手机版 http://muchong.com/3g ]
4楼2015-08-08 19:35:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)


引用回帖:
4楼: Originally posted by luzihen at 2015-08-08 19:35:11
问题就是paw似乎并不是淹势。
不过很多文献里面就这么写了,包括一些prl,所以应该不是很要要紧。
...

PAW确实就是生成pseudopotential的一种方法啊。除此之外用的多的还有NCPP 和 USPP。。。你去看R.M. Martin的electronic structure 11.11 写得很清楚啊。。。。GGA是交换关联泛函Exc[n]的一种近似,由Exc[n]可以求出Vxc,然后Vxc是属于potential的,因此在生成pseudopotential的时候考虑在里面了,所以可以直接写PAW-PBE pseudopotential  反正我是这么理解的。。
7楼2015-08-10 06:47:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
awmc20085楼
2015-08-08 22:00   回复  
五星好评  顶一下,感谢分享!
Quan.6楼
2015-08-09 06:06   回复  
五星好评  顶一下,感谢分享!
2016-03-24 17:49   回复  
五星好评  顶一下,感谢分享!
seer9楼
2016-09-12 21:07   回复  
五星好评  顶一下,感谢分享!
q96524001010楼
2016-09-24 19:17   回复  
五星好评  顶一下,感谢分享!
2016-10-26 20:03   回复  
五星好评  顶一下,感谢分享!
墨兰清12楼
2016-11-28 15:42   回复  
五星好评  顶一下,感谢分享!
wurenglai13楼
2016-11-30 09:53   回复  
五星好评  顶一下,感谢分享!
wurenglai14楼
2016-11-30 09:54   回复  
顶一下,感谢分享!
Yolonuts15楼
2017-02-27 14:31   回复  
五星好评  顶一下,感谢分享!
邹辉16楼
2018-08-23 12:18   回复  
五星好评  顶一下,感谢分享!
env17楼
2019-11-15 09:33   回复  
五星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 luzihen 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+3 BruceChum 2024-06-15 20/1000 2024-06-18 01:22 by BruceChum
[考博] 2025考博 +7 自强不息a?a 2024-06-15 12/600 2024-06-17 22:51 by 自强不息a?a
[找工作] 药学硕士找不到工作,打算去做科研助理了 +5 pom戴墨镜 2024-06-14 12/600 2024-06-17 21:26 by 雨后春笋!
[基金申请] 化学会年会改成一年开一次了? +7 babu2015 2024-06-12 15/750 2024-06-17 20:28 by 落笔画枫
[论文投稿] 审稿问题:为什么荧光激发波长和紫外吸收波长差的大? 10+5 sdawege 2024-06-14 10/500 2024-06-17 18:54 by HH-探针
[基金申请] 面青地会评时间 +8 tanjydd 2024-06-15 8/400 2024-06-17 17:08 by 小龙虾2008
[有机交流] 车间生产,真空度很高,温度很高,但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 38/1900 2024-06-17 16:46 by 科研虫子小助手
[基金申请] 博后面上今天有bug可以看到是否资助? +22 lyfbangong 2024-06-12 33/1650 2024-06-17 16:17 by carolloo
[有机交流] 苯酚和炔醇反应,mitsunobu反应 5+3 TONGMEIMEI 2024-06-12 9/450 2024-06-17 15:00 by dschong
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +10 zju2000 2024-06-12 16/800 2024-06-17 14:54 by 大天尊
[催化] 镍负载氧化铝的保存问题 8+3 lwn0130 2024-06-15 4/200 2024-06-17 10:48 by adaihao
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[基金申请] 博后基金,博管会会提前知道消息吗? +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[硕博家园] 考博应该从哪里开始准备 +6 亮姐2019 2024-06-11 8/400 2024-06-13 12:53 by 帕秋莉GO
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
信息提示
请填处理意见