24小时热门版块排行榜    

查看: 5121  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 重要来源:本周末出结果 +3 瞬息宇宙 2026-08-12 3/150 2026-08-13 06:40 by 玲珑133
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +8 majunge000 2026-08-10 9/450 2026-08-12 22:17 by 虫友是什么虫
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见