| 查看: 1880 | 回复: 2 | ||
| 【奖励】 本帖被评价2次，作者maomao7910增加金币 1.6 个 | ||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||
[资源]
自发极化设置LCALCPOL时存在的问题
|
||
|
如题，用LBERRY = .TRUE. 的方法候博已经给出详细例子， http://valenhou.blog.edu.cn/2008/133226.html ，我不再赘述。vasp5.2后加入了LCALCPOL=.TRUE.开关，可以一次性算出电子的极化贡献。不过目前这样设置后在OUTCAR中并没有出来如手册上所讲的 Ionic dipole moment: p[ion]=( 2.25510 2.25510 2.25510 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.00000 ) electrons Angst 两项，而只是给出 CALCP: cpu time 334.26: real time 341.92 Total macroscopic dipole moment: p[mac]=( 0.00000 0.00000 0.19288 ) electrons Angst 这一项，是电子极化的贡献，那么离子极化的贡献如何解决？给两个办法，如上面候博讲的求原子移动和未移动时的berry phase，但这里只求一次，即 IGPAR = 1 or2 or3就可以了。另一种方法。如果你看关于berry phase 的手册 http://cms.mpi.univie.ac.at/vasp/Berry_phase/node4.html 你会发现离子极化项其实就是 （价电荷＊原子位移）所以，自己手动一算就出来了。先在势函数里面找到你要移动的原子的价电子数 ZVAL＊（原子移动前的坐标－原子移动后的坐标）就等于离子极化的贡献。 欢迎指正！ [ Last edited by youzhizhe on 2011-7-6 at 12:59 ] |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有4人回复
-
时间戳又变了
已经有17人回复
-
科研孤儿太难了
已经有10人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 丙烯酸酯压敏胶粘剂乳液聚合时加碱中和的问题？
已经有17人回复
- 制备纳米颗粒时前驱体中pH问题
已经有10人回复
- 关于加氢后去除催化剂的问题
已经有7人回复
- 能从理论上导出极化曲线吗？
已经有5人回复
- 关于MS中通过 Brillouin Zone Path 设置和修改的问题
已经有31人回复
- 波恩（Born）模型存在什么问题？
已经有1人回复
- 电池测试中的极化现象和安全问题
已经有6人回复
- 酶法测定葡萄糖时存在的干扰问题
已经有1人回复
- 【求助】请CAD高手指点----关于CAD设置问题
已经有4人回复
- 求助关于Wifi热点设置的问题
已经有7人回复
- origin 画图时坐标设置的问题
已经有2人回复
- Word问题，设置单双页不同之后，只有偶数页有页码。求高手解决一下。
已经有1人回复
- 【求助】用castep计算BaTiO3的自发极化
已经有8人回复
- 【求助】如何解决自旋极化计算时的难以收敛问题？
已经有16人回复
高见！ 2楼2012-01-07 19:20:43
3楼2014-05-29 09:55:18
