应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 有人对自己做PAW有兴趣吗？
774/8返回列表首页上一页234567下一页末页
查看: 5083  |  回复: 76
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tangosnow

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我现在在做GIPAW，好像和GPAW还是两个概念。我其实不知道PAW是啥意思，汗，我是用pwscf的ld1.x生成的。
问个问题，同一个元素的赝势在某种结构中算出来和文献对应较好，在另外结构中，就不一定能对应得好了？
赞一下 回复此楼
31楼2012-05-23 21:57:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
31楼: Originally posted by tangosnow at 2012-05-23 21:57:54:
我现在在做GIPAW，好像和GPAW还是两个概念。我其实不知道PAW是啥意思，汗，我是用pwscf的ld1.x生成的。
问个问题，同一个元素的赝势在某种结构中算出来和文献对应较好，在另外结构中，就不一定能对应得好了？

GIPAW是？
你的问题就是说PAW方法的可移植性问题了？这的确需要仔细考虑，尤其是3d过渡族金属
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
32楼2012-05-23 23:11:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
26楼: Originally posted by 信任开花 at 2012-05-23 17:41:20:
对于POTCAR中有关PAW的信息

不知道可有什么提示的没有？

我写信问一个老外，结果他回复说所有的PAW信息都在POTCAR中

某有给我更详细的说明。。。

有的，比较重要就是那几个cutoff了。
正在看Kresse PRB 1999那篇文章，不过我觉得奇怪的是它的PAW Fe_sv构型是，而不是文章中所说的
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
33楼2012-05-23 23:21:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★
franch: 金币+2, 鼓励交流，， 呵呵 2012-05-24 01:24:59
在Kresse PRB 1999的文章中，提到了这么几个关键地方
（1）partial waves and projectors
在截断半径内用RRKJ scheme将PS partial wave写成

这里的对应着POTCAR的description部分的关键字RCUT
（2）pseudized core charges and the local pseudopotentail
local atomic pseudopotential在截断半径内取成以下形式

截断半径对应着POTCAR中关键字RCLOC
pseudized core charges在截断半径内取成以下形式

截断半径对应着POTCAR中关键字RPACOR
（3）compensation charge
compensation function 定义为

并且在半径内满足

与radial integration sphere 有关，在POTCAR中关键字为RDEP，而由RDEP/RAUG决定。

[ Last edited by iamikaruk on 2012-5-24 at 00:26 ]
赞一下(1人) 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
34楼2012-05-24 00:23:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★ ★
uuv2010: 金币+3, 鼓励高水平学术讨论，赞实践能力 2012-05-24 20:44:42
首先声明：以下结果与讨论仅限于学习和学术目的，其它任何行为与发帖个人无关
1. PAW生成
参照VASP中Fe_sv的POTCAR设置，用atompaw造了一个Fe的PAW出来
Fe的参考电子构型为
所用的input file如下：
CODE:
Fe 26
LDA-PW scalarrelativistic loggrid 2001
4 3 3 0 0 0 ! Fe:[1s2 2s2 2p6] 3s2 3p6 3d6 4s2
3 2 7
4 0 1
0 0 0
c ! 1s
c ! 2s
v ! 3s
v ! 4s
c ! 2p
v ! 3p
v ! 3d
2
1.909 1.4685 0.801 1.600
n
y
1.0
n
y
0.0
n
custom rrkj vanderbiltortho besselshape
bessel
1.9
1.9
1.7
1.7
1.9
1.9
1
3 2 6
4 0 2
0 0 0
2
prtcorewf nooptim usexcnhat  ! generation of abinit PP file with Kresse's formalism
0

projector、pseudo wavefunction和logderiv就不放在这里了，有兴趣的可以自己生成看看。

2. 测试结果
测试中用了三个软件，分别是FPLAPW的Fleur，VASP和abinit，计算了无磁BCC Fe的能量-晶格常数以及实验晶格常数下铁磁BCC Fe的spin
FPLAPW的muffin tin radius取为PAW的augmentation半径1.909 Bohr，没有core-leakage，平面波截止取为4.5 a.u.，k-points为8x8x8
VASP和abinit的平面波截断能量分别取为550 eV和20 Hartree，abinit中PAW的能量cutoff取为40 Hartree，k-points均为8x8x8，空带数取为16，smear采用Methfessel–Paxton方法，所加smearing为0.01 eV。
无磁BCC Fe的能量-晶格常数对比结果请见附件，可以看到VASP的PAW质量好的吓人，跟FPLAPW的结果几乎重合了
实验晶格常数下BCC Fe的spin分别为FPLAPW:2.14645 Bohr， VASP:2.166 Bohr，abinit:2.174 Bohr
对于这个结果，各位有什么看法？
赞一下(3人) 回复此楼

» 本帖附件资源列表

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

» 本帖已获得的红花（最新10朵）

huangll99
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
35楼2012-05-24 20:09:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
uuv2010: 金币+1, 鼓励讨论 2012-05-24 20:45:06
引用回帖:
35楼: Originally posted by iamikaruk at 2012-05-24 20:09:09:
首先声明：以下结果与讨论仅限于学习和学术目的，其它任何行为与发帖个人无关
1. PAW生成
参照VASP中Fe_sv的POTCAR设置，用atompaw造了一个Fe的PAW出来
Fe的参考电子构型为 3s^2 3p^6 3d^7 4s^1
所用的input fi

强烈要求上各种弹性数据～或者上声子谱对比～～～～～
一般生成的质量尚可的赝势在总能方面的差别都不会非常大，但是如果考虑总能的梯度变化，则可能会显示出优劣。所以对于弹性数据和声子色散，都是对总能的微分形式，对检验赝势质量非常重要。尤其是一整套稳定态单质的弹性数据，一般是比较标准的赝势测试结果～
俺自己生成的超软赝势，只做E-V曲线时，看起来比ELK的结果好，但是一做弹性计算，就比ELK差点事了～不过还好差别不大，并且GGA对于纯金属的计算很多时候就不如LDA～
至于磁性数据，没有机会接触磁性计算，木有发言权啊……
赞一下 回复此楼
思想重于技巧，内涵重于表象
36楼2012-05-24 20:30:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送鲜花一朵
引用回帖:
35楼: Originally posted by iamikaruk at 2012-05-24 20:09:09:
首先声明：以下结果与讨论仅限于学习和学术目的，其它任何行为与发帖个人无关
1. PAW生成
参照VASP中Fe_sv的POTCAR设置，用atompaw造了一个Fe的PAW出来
Fe的参考电子构型为 3s^2 3p^6 3d^7 4s^1
所用的input fi

膜拜一下大牛，虽然都没看明白怎么回事
赞一下 回复此楼
37楼2012-05-24 20:31:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
36楼: Originally posted by souledge at 2012-05-24 20:30:21:
强烈要求上各种弹性数据～或者上声子谱对比～～～～～
一般生成的质量尚可的赝势在总能方面的差别都不会非常大，但是如果考虑总能的梯度变化，则可能会显示出优劣。所以对于弹性数据和声子色散，都是对总能的微

算弹性数据和声子谱之类的实在是太累了，先简单和其它对一下能量再说。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
38楼2012-05-24 21:00:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

小木虫灌水大队队长
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
内容已删除
回复此楼
人生是一场旅行，在欣赏风景的同时，别忘了你也是别人的风景。
39楼2012-05-24 21:51:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

小木虫灌水大队队长
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
33楼: Originally posted by iamikaruk at 2012-05-23 23:21:50:
有的，比较重要就是那几个cutoff了。
正在看Kresse PRB 1999那篇文章，不过我觉得奇怪的是它的PAW Fe_sv构型是3p^63d^74s^1，而不是文章中所说的3p^63d^{6.5}4s^{1.5}

我没有看这篇文章，不过，获取赝势的时候，可以取 4s1.5 这样的构型。
赞一下 回复此楼
人生是一场旅行，在欣赏风景的同时，别忘了你也是别人的风景。
40楼2012-05-24 21:55:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 iamikaruk 的主题更新
774/8返回列表首页上一页234567下一页末页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭