24小时热门版块排行榜    

查看: 5205  |  回复: 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.
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +5 13108017953 2026-09-21 5/250 2026-09-27 16:46 by Leogzhya
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +3 zju2000 2026-09-26 3/150 2026-09-27 16:26 by lala1008
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +5 oxidpl 2026-09-20 7/350 2026-09-27 10:43 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[教师之家] 课题组招2027级博士 上海工程技术大学 激光智能制造方向 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:37 by oxidpl
信息提示
请填处理意见