| 查看: 2033 | 回复: 11 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
gemucai木虫 (正式写手)
|
[求助]
用自己生成的赝势算能带效果很好，算声子却一团糟
|
|||
| 我用的是PWSCF,能带计算和声子计算都比较熟悉了。但最近我自己生成赝势之后，算出的能带与别人文献比较几乎无差别，但gamma点的声子频率却都是负的几千1/cm。我做过了优化，力和压力都几乎为0了，但还是有如此大的虚频。为何电子和声子计算会差别如此之大？有可能是我的赝势哪里没做好才导致了DFPT过程的异常？ |
回复此楼» 猜你喜欢
求调剂
已经有3人回复
-
化工270求调剂
已经有7人回复
-
276求调剂
已经有5人回复
-
295求调剂
已经有8人回复
-
材料学硕318求调剂
已经有16人回复
-
268求调剂
已经有4人回复
-
此成果不能导入原因：元数据必填信息不完整，可 进行补充。
已经有5人回复
-
材料工程专硕283求调剂
已经有6人回复
-
281求调剂
已经有8人回复
-
求调剂
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 采用vasp ，phononpy或frophon计算声子谱的具体方法是什么？
已经有5人回复
- 今天用CASTEP计算声子谱，出现下列提示，有能看懂的大侠帮忙解释下
已经有11人回复
- 为何用晶格常数的实验值计算声子？
已经有18人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- Abinit计算声子问题
已经有10人回复
- ZnO的声子计算
已经有18人回复
- castep计算Fe2O3声子谱失败，原因求教~~
已经有10人回复
- 自旋体系的声子计算应该用什么软件呢？
已经有5人回复
- 【求助】使用vasp进行带心冻结声子计算的几个问题~~~
已经有9人回复
- 【求助】有人用有限差分方法计算过G点声子谱吗？
已经有10人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
- 【求助】请问大家在做声子计算时，这个 qpt是按照什么设置的？？？？
已经有9人回复
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
|
那两个赝势好像有问题。我测试4p5 5s1组态的时候，ld1.x报错 WARNING! Expected number of nodes: 0= 3- 2- 1, number of nodes found: 3. Setting wfc to zero for this iteration. (This warning will only be printed once per wavefunction) WARNING! Expected number of nodes: 0= 1- 0- 1, number of nodes found: 3. Setting wfc to zero for this iteration. (This warning will only be printed once per wavefunction) Message from routine run_pseudo: Errors in PS-KS equations WARNING! Expected number of nodes: 0= 2- 1- 1, number of nodes found: 3. Setting wfc to zero for this iteration. (This warning will only be printed once per wavefunction) 相反，对于Mo的阳离子组态却误差非常小： n l nl e AE (Ry) e PS (Ry) De AE-PS (Ry) 3 2 4D 1( 5.00) -0.78364 -0.78177 -0.00187 1 0 5S 1( 0.00) -0.74605 -0.74501 -0.00105 2 1 5P 1( 0.00) -0.44623 -0.44648 0.00024 ----------------------------------------------------------------------------------------- n l nl e AE (Ry) e PS (Ry) De AE-PS (Ry) 3 2 4D 1( 4.00) -1.54772 -1.55039 0.00267 1 0 5S 1( 0.00) -1.32339 -1.31887 -0.00452 2 1 5P 1( 0.00) -0.94261 -0.94129 -0.00132 至于放在晶体中的情形还没有测试。 |
7楼2012-04-17 11:10:34
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
2楼2012-04-07 10:05:19
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
3楼2012-04-07 14:02:42
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+2, 谢谢分享 2012-04-07 16:13:37
gemucai: 金币+7, ★★★很有帮助, 我又做了一个稍微好点的，会跟这两个比较。先谢谢了！ 2012-04-07 16:16:17
fzx2008: 金币+2, 谢谢分享 2012-04-07 16:13:37
gemucai: 金币+7, ★★★很有帮助, 我又做了一个稍微好点的，会跟这两个比较。先谢谢了！ 2012-04-07 16:16:17
|
Mo的赝势很不好生成，尤其考虑到半芯态时。 这里有两个从CPMD转来的LDA模守恒赝势，06的是没有半芯态，14的是有半芯态的。我近来在做其他的计算，还没时间测试这两个赝势的效果和差别，如果有兴趣，可以测试一下～ 顺便，如果测试了，还希望早点发上测试结果图来，也有利于其他使用QE的人～免得赝势不好让太多人浪费了时间，或者免得大家都得自己跑测试～ |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Mo.LDA06.UPF
- 附件 2 : Mo.LDA14.UPF
2012-04-07 16:08:30, 295.5 K
2012-04-07 16:08:34, 295.5 K
4楼2012-04-07 16:08:37
10