应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 急救：用VASP算电偶极矩出问题
81/1返回列表
查看: 3560  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

May.Be

新虫 (初入文坛)

[求助] 急救：用VASP算电偶极矩出问题

我用VASP算一个分子HF的电偶极矩，算出来以后比文献的大几百倍，是不是我的参数设置有问题，特别是DIPOL这个参数
INCAR如下：
SYSTEM = HF
PREC=Medium
EDIFF  = 1E-6
ISMEAR = 0   
DIPOL=1.0  1.0  1.0
LCALCPOL= .TRUE
POSCAR如下：
CIF file                                
   1.00000000000000     
    10.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000   10.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000   10.0000000000000000
   H    F
     1     1
Direct
  0.4688379198180487  0.4880179885076180  0.4823762497841139
  0.5457720731819519  0.5262220074923755  0.5200837532158821

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
KPOINTS如下：
Automatic generation
0
Monhkorst-Pack
1            1        1
0.0        0.0      0.0
算出的电偶极矩项如下：
Ionic dipole moment: p[ion]=(    27.10758    28.28428    28.77037 ) electrons Angst

Total electronic dipole moment: p[elc]=(     3.21352     1.87427     1.38611 ) electrons Angst
太大了，还有这个单位 electrons Angst是不是1.6*10（-29）Cm,哪位大神能指点一下，只要有帮助金币可以加。
回复此楼

» 收录本帖的淘帖专辑推荐

第一性

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-06-17 08:48:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)
现在都不会计算偶极矩的飘过……
赞一下 回复此楼
长枪一横花飘零，松风追月伴我行。无双人间世如梦，倾奇万世永留名。
2楼2014-06-19 20:45:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaobin521

木虫 (小有名气)
这个还是挺简单的，单位是没错的，但在计算的时候，需要减掉晶格常数，比如对27.10758 需要减掉30  得到的结果是p[ion]=(-2.9   -1.7   -1.3)
                                                  p[elc]=(3.2   1.9   1.4)
这样的话就属于正常了，在算总的极矩时，还要注意除以他们的体积，得到单位为C/m2，这样的话就圆满了。希望我的回复对你有所帮助。。。
赞一下(1人) 回复此楼
好好生活，好好工作。
3楼2015-04-26 10:53:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)
引用回帖:
3楼: Originally posted by gaobin521 at 2015-04-26 10:53:14
这个还是挺简单的，单位是没错的，但在计算的时候，需要减掉晶格常数，比如对27.10758 需要减掉30  得到的结果是p=(-2.9   -1.7   -1.3)
                                                  p=(3.2   1.9   1.4)
...

兄弟，不仅回答的好！人长的也帅！
赞一下 回复此楼
长枪一横花飘零，松风追月伴我行。无双人间世如梦，倾奇万世永留名。
4楼2015-04-26 14:26:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gswylq

木虫 (正式写手)
引用回帖:
3楼: Originally posted by gaobin521 at 2015-04-26 10:53:14
这个还是挺简单的，单位是没错的，但在计算的时候，需要减掉晶格常数，比如对27.10758 需要减掉30  得到的结果是p=(-2.9   -1.7   -1.3)
                                                  p=(3.2   1.9   1.4)
...

请问，为什么要减去晶格常数？
赞一下 回复此楼
5楼2015-04-28 16:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houxiami

银虫 (初入文坛)
引用回帖:
5楼: Originally posted by gswylq at 2015-04-28 16:43:00
请问，为什么要减去晶格常数？...

同问这个问题。。。还有减掉的晶格常数怎么是30 ？
这种情况是不是和DIPOL选取的数值有一定的定量关系呢？
赞一下(1人) 回复此楼
6楼2016-03-22 20:40:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿克飞顿

金虫 (正式写手)
引用回帖:
2楼: Originally posted by 前田庆次 at 2014-06-19 20:45:18
现在都不会计算偶极矩的飘过……

已经半年了
回复此楼
7楼2019-12-30 17:14:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿克飞顿

金虫 (正式写手)
引用回帖:
3楼: Originally posted by gaobin521 at 2015-04-26 10:53:14
这个还是挺简单的，单位是没错的，但在计算的时候，需要减掉晶格常数，比如对27.10758 需要减掉30  得到的结果是p=(-2.9   -1.7   -1.3)
                                                  p=(3.2   1.9   1.4)
...

请问下是将各方向的p[ion]和p[elc]分别相加吗，再除以体积，得到各个方向的总偶极矩？还有就是平时我们说的电偶极矩也就是偶极矩和这里的哪个数值对应呢，总的还是p[elc]?谢谢
赞一下 回复此楼
8楼2020-01-03 10:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 May.Be 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭