| 查看: 1566 | 回复: 4 | ||
jk射手座新虫 (著名写手)
|
[求助]
求教LCALCPOL计算极化的结果怎么处理的，多出来一行
|
|
我先贴下我的INCAR： SYSTEM=PTOLAO ENCUT=450 ISMEAR=-5 SIGMA=0.02 EDIFF=1E-3 ISTART=0 ICHARG=2 IBRION=2 PREC=Accurate LCALCPOL=.TRUE. DIPOL=0.3 0.3 0.3 IDIPOL=3 LDAU=.TRUE. LDAUTYPE=2 LDAUL=-1 -1 3 -1 -1 LDAUU=0 0 6 0 0 LDAUJ=0 0 0 0 0 LDAUPRINT=2 K网格是自动 6 6 4 的 POTCAR用的是PAW的 POSCAR有两个， 其一： data_PTOLAO 1.00000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 8.1400000000000000 Pb Ti La Al O 1 1 1 1 6 Selective dynamics Direct 0.00 0.00 0.00 F F T 0.50 0.50 0.25 F F T 0.00 0.00 0.50 F F T 0.50 0.50 0.75 F F T 0.50 0.00 0.25 F F T 0.00 0.50 0.25 F F T 0.50 0.50 0.00 F F T 0.50 0.00 0.75 F F T 0.00 0.50 0.75 F F T 0.50 0.50 0.50 F F T 其二： data_PTOLAO 1.00000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.8100000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 7.9772000000000000 Pb Ti La Al O 1 1 1 1 6 Selective dynamics Direct 0.0000000000000000 0.0000000000000000 -0.0444396041984977 F F T 0.5000000000000000 0.5000000000000000 0.2327315837951796 F F T 0.0000000000000000 0.0000000000000000 0.4971976926301170 F F T 0.5000000000000000 0.5000000000000000 0.7423875204210301 F F T 0.5000000000000000 0.0000000000000000 0.2778793310366013 F F T 0.0000000000000000 0.5000000000000000 0.2778793310366013 F F T 0.5000000000000000 0.5000000000000000 0.0092994307667890 F F T 0.5000000000000000 0.0000000000000000 0.7482258333968913 F F T 0.0000000000000000 0.5000000000000000 0.7482258333968913 F F T 0.5000000000000000 0.5000000000000000 0.5106130477183970 F F T 然后计算出的结果我提取的 'moment'.提取出来是酱婶儿的： angular momentum for each species LDAUL = -1 -1 3 -1 -1 dipolmoment 0.000000 0.000000 -21.359368 electrons x Angstroem Ionic dipole moment: p[ion]=( 12.57300 12.57300 -20.04233 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.28231 ) electrons Angst [root@node16 -1.5calc]# grep 'moment' OUTCAR1 angular momentum for each species LDAUL = -1 -1 3 -1 -1 dipolmoment 0.000000 0.000000 -24.451270 electrons x Angstroem Ionic dipole moment: p[ion]=( 12.57300 12.57300 -24.01300 ) electrons Angst Total electronic dipole moment: p[elc]=( 0.00000 0.00000 0.15750 ) electrons Angst [root@node16 -1.5calc]# grep 'volume' OUTCAR2 volume/ion in A,a.u. = 11.64 78.54 volume of cell : 116.39 volume of cell : 116.39 请问 electrons x Angstroem那行的dipolemoment 和下边两行电子与离子的偶极矩是啥关系啊？多设置了一个IDIPOLE=3之后就多了 electrons x Angstroem的那一行； 搜了些相关的说明与讨论，如果我理解为： electrons x Angstroem那行dipolemoment是 总P（也就是p[ion]+p[elc]）在Z方向上的分量，这样理解对吗？我计算了下，分别是3.1027和4.20832，那么其与z轴的夹角为arccos(3.1027/4.20832),大概比45°小一点；我算的是两个四方相组合的超晶格，总P的方向应该离体对角线不远，不知这样理解对不？ 请各路大神帮忙指点下！！！！ （另外还想问下，LBERRY大家都知道是BERRY PHASE方法，那么设置LCALCPOL算极化，这个方法有名字吗？？） [ Last edited by jk射手座 on 2014-7-2 at 10:12 ] |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
明天放榜？
已经有4人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 氧还原计算塔菲尔（Tafel）斜率，极限扩散电流是从LSV极化曲线上读出来吗？
已经有13人回复
- 用DPF CAM-B3LYP方法计算极化率的输入文件要怎么写
已经有5人回复
- 求教极化电阻是怎么个概念
已经有4人回复
- 关于电池极化和电极极化的问题？求教
已经有7人回复
- 半导体的极化率怎么计算
已经有6人回复
- 【求助】VASP计算出的电极化强度数据，怎么做后处理
已经有15人回复
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
2楼2014-07-02 10:12:43
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
3楼2014-07-04 15:41:47
youqin5912
金虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 1183.4
- 散金: 785
- 红花: 15
- 帖子: 397
- 在线: 594.4小时
- 虫号: 1141030
- 注册: 2010-11-07
- 性别: GG
- 专业: 同步辐射技术及其应用
楼主得到答案没？我顶一下 4楼2014-07-11 21:49:16
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
5楼2014-07-16 11:29:08
