| 查看: 5381 | 回复: 76 | |||||||||
iamikaruk木虫 (著名写手)
|
[交流]
有人对自己做PAW有兴趣吗？ 已有17人参与
|
||||||||
|
刚刚有个想法，其实还是从版上开展的阅读源码活动而得来的想法，大家都知道使用的赝势质量的好坏对于结果有着比较重要的影响，而vasp在这一点上以其优秀的PAW赢得了广大的用户群。但是在某些功能上一些开源计算软件，例如abinit和pwscf等，都有自己的特色，但是它们使用的还是比较老的模守恒或者超软赝势方法。 所以我的问题就是，大家觉得是否有必要做一套标准的PAW（比如说交换-关联取为LDA-PZ或GGA-PBE，包含大部分常见元素）为自己和他人的研究提供便利呢？不过我也知道这种工作工作量非常大，而且很多时候吃力不讨好。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | Photochemistry | 科研宝典 | good resources |
| 计算化学-辅助工具 | 建模 |
» 猜你喜欢
0703化学调剂，求各位老师收留
已经有8人回复
-
070303一志愿西北大学学硕310找调剂
已经有8人回复
-
085600材料与化工 求调剂
已经有14人回复
-
本人考085602 化学工程 专硕
已经有13人回复
-
材料与化工求调剂
已经有3人回复
-
283求调剂
已经有13人回复
-
312求调剂
已经有4人回复
-
285求调剂
已经有6人回复
-
材料与化工专硕调剂
已经有3人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- =.= Zr_sv的PBE-PAW赝势在计算中出现了个跳跃点……
已经有21人回复
- 固体理论里的APW（缀加平面波方法）和vasp里的PAW方法是什么关系呢？
已经有9人回复
- 求助PAW91-Er 赝势 急急急！！！
已经有5人回复
- 请问有人有稀土Ce Pr Nd Sm Eu Gd Tb包含4f价电子的PAW赝势吗
已经有3人回复
- 【求助】求助Ge和Te元素的PAW势
已经有5人回复
- 【求助】vasp采用PAW势，是不是用的APW方法？
已经有5人回复
- 【求助】一个小例子，请教uspp和paw之区别
已经有19人回复
- 【求助】求ABINIT中的PAW赝势文件
已经有6人回复
- 【求助】请教：ABINIT中的PAW赝势怎么获取？
已经有8人回复
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
★ ★ ★
uuv2010: 金币+3, 鼓励高水平学术讨论，赞实践能力 2012-05-24 20:44:42
uuv2010: 金币+3, 鼓励高水平学术讨论，赞实践能力 2012-05-24 20:44:42
|
首先声明：以下结果与讨论仅限于学习和学术目的，其它任何行为与发帖个人无关 1. PAW生成 参照VASP中Fe_sv的POTCAR设置，用atompaw造了一个Fe的PAW出来 Fe的参考电子构型为 所用的input file如下： 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朵）
huangll9935楼2012-05-24 20:09:09
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★
franch: 金币+5, 鼓励交流 2012-05-23 19:46:29
franch: 金币+5, 鼓励交流 2012-05-23 19:46:29
|
好的，上传的资料分别是atompaw，gpaw和vasp的实现 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : atompawEqns.pdf
- 附件 2 : paw_note.pdf
- 附件 3 : Kresse99.pdf
2012-05-23 16:26:04, 193.83 K
2012-05-23 16:26:13, 549.42 K
2012-05-23 16:26:19, 282.72 K
21楼2012-05-23 16:26:22
lzllzl
木虫 (小有名气)
- 应助: 10 (幼儿园)
- 金币: 2493.8
- 散金: 18
- 红花: 4
- 帖子: 102
- 在线: 393.1小时
- 虫号: 300731
- 注册: 2006-11-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
PSlibrary就是为pwscf干这的吧，但不限于PAW： http://qe-forge.org/projects/pslibrary [ Last edited by lzllzl on 2012-5-22 at 08:17 ] |
8楼2012-05-22 08:16:08
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
|
8楼给的链接不错，可以研究一下。 另外也可以参考GPAW的各种参数设置：https://wiki.fysik.dtu.dk/gpaw/ 标准化的做法我已经有个大概的思路，基于atompaw+libxc生成不同交换-关联近似下的PAW，然后用VASP和基于FPLAPW的Wien2k/ELK/Fleur等软件进行验证，工作量很大。我先试着做几个PAW看，总结一下经验，然后再看是否可行。 我记得有篇论文讨论过用遗传算法的方法寻找PAW的截断半径等参数以得到比较好的PAW。 |
9楼2012-05-22 09:33:28
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
17楼2012-05-23 11:03:17
iamikaruk
木虫 (著名写手)
- 1ST强帖: 1
- 应助: 11 (小学生)
- 金币: 4950.3
- 散金: 118
- 红花: 76
- 帖子: 1362
- 在线: 543.8小时
- 虫号: 1315985
- 注册: 2011-06-05
- 专业: 无机非金属类电介质与电解
★ ★
franch: 金币+2, 鼓励交流，， 呵呵 2012-05-24 01:24:59
franch: 金币+2, 鼓励交流，， 呵呵 2012-05-24 01:24:59
|
在Kresse PRB 1999的文章中，提到了这么几个关键地方 （1）partial waves and projectors 在截断半径 这里的 （2）pseudized core charges and the local pseudopotentail local atomic pseudopotential在截断半径 截断半径 pseudized core charges在截断半径 截断半径 （3）compensation charge compensation function 并且在半径 [ Last edited by iamikaruk on 2012-5-24 at 00:26 ] |
34楼2012-05-24 00:23:43
gongjijun
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 4161.8
- 散金: 48
- 红花: 5
- 帖子: 507
- 在线: 846小时
- 虫号: 935467
- 注册: 2009-12-30
- 专业: 凝聚态物性 II ：电子结构
2楼2012-05-22 00:11:29
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3681.7
- 散金: 5802
- 红花: 79
- 沙发: 35
- 帖子: 4560
- 在线: 916小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
3楼2012-05-22 00:30:28
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
这个想法很好，虽然估计有点难以实现，因为对于平面波赝势方法来说，最珍贵的资源就是赝势了，尤其是优化好的赝势，速度快，精度高，是能拿来卖钱的～ VASP的代码要说也不见得就比那些免费的或者CASTEP好多少，但是一整套高质量并且非常软的PAW赝势，再配合不是太吓人的价格，就征服了N多的用户。 尤其那些重元素和f轨道元素，非常难以生成，并且生成成功了，也很难有人会舍得免费共享～ 另外，看了一下，对于生成PAW赝势，PWscf和Abinit共通的软件应该首选AtomPAW，不知道有多少人用过这个哦～顺便，PWscf的官方提供了PPLibrary，里面有不少PAW赝势生成的例子，虽然使用的是PWscf自带的ld1.x进行生成，但是对于学习造赝势来说，是非常好的资料了～ |
4楼2012-05-22 01:33:49
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
5楼2012-05-22 07:18:14
gongjijun
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 4161.8
- 散金: 48
- 红花: 5
- 帖子: 507
- 在线: 846小时
- 虫号: 935467
- 注册: 2009-12-30
- 专业: 凝聚态物性 II ：电子结构
6楼2012-05-22 07:38:05
fzx2008
荣誉版主 (著名写手)
- 1ST强帖: 9
- 应助: 272 (大学生)
- 贵宾: 1.045
- 金币: 21814.1
- 散金: 710
- 红花: 33
- 沙发: 1
- 帖子: 1808
- 在线: 1601小时
- 虫号: 621341
- 注册: 2008-10-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
7楼2012-05-22 08:13:56
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
10楼2012-05-22 09:37:03
