24小时热门版块排行榜    

查看: 6083  |  回复: 76
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

iamikaruk

木虫 (著名写手)

[交流] 有人对自己做PAW有兴趣吗? 已有17人参与

刚刚有个想法,其实还是从版上开展的阅读源码活动而得来的想法,大家都知道使用的赝势质量的好坏对于结果有着比较重要的影响,而vasp在这一点上以其优秀的PAW赢得了广大的用户群。但是在某些功能上一些开源计算软件,例如abinit和pwscf等,都有自己的特色,但是它们使用的还是比较老的模守恒或者超软赝势方法。
所以我的问题就是,大家觉得是否有必要做一套标准的PAW(比如说交换-关联取为LDA-PZ或GGA-PBE,包含大部分常见元素)为自己和他人的研究提供便利呢?不过我也知道这种工作工作量非常大,而且很多时候吃力不讨好。
回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
已阅   回复此楼   关注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
对于这个结果,各位有什么看法?

» 本帖附件资源列表

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

» 本帖已获得的红花(最新10朵)

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

gongjijun

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这个想法非常好,就是不知道如何下手?
2楼2012-05-22 00:11:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主


小木虫: 金币+0.5, 给个红包,谢谢回帖
不错,这是一个好想法。
需要提出一个方案出来,让大家做起来
不求速度,但要有行动
3楼2012-05-22 00:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这个想法很好,虽然估计有点难以实现,因为对于平面波赝势方法来说,最珍贵的资源就是赝势了,尤其是优化好的赝势,速度快,精度高,是能拿来卖钱的~
VASP的代码要说也不见得就比那些免费的或者CASTEP好多少,但是一整套高质量并且非常软的PAW赝势,再配合不是太吓人的价格,就征服了N多的用户。
尤其那些重元素和f轨道元素,非常难以生成,并且生成成功了,也很难有人会舍得免费共享~
另外,看了一下,对于生成PAW赝势,PWscf和Abinit共通的软件应该首选AtomPAW,不知道有多少人用过这个哦~顺便,PWscf的官方提供了PPLibrary,里面有不少PAW赝势生成的例子,虽然使用的是PWscf自带的ld1.x进行生成,但是对于学习造赝势来说,是非常好的资料了~
思想重于技巧,内涵重于表象
4楼2012-05-22 01:33:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见