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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 急救：用VASP算电偶极矩出问题
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 3017  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +5 Amber00 2026-03-31 5/250 2026-03-31 22:43 by zhengzhi8909
[考研] 352分-085602-一志愿985 +6 海纳百川Ly 2026-03-29 6/300 2026-03-31 21:06 by yuq
[考研] 285求调剂 +3 FZAC123 2026-03-30 3/150 2026-03-31 17:49 by 热情沙漠
[考研] 307分求调剂 +6 (o~o) 2026-03-31 6/300 2026-03-31 17:22 by 唐沐儿
[考研] 材料科学与工程调剂 +13 深V宿舍吧 2026-03-30 14/700 2026-03-31 16:28 by hypershenger
[考研] 354求调剂 +3 lxb598 2026-03-31 4/200 2026-03-31 13:42 by sophie2180
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-31 4/200 2026-03-31 11:56 by jp9609
[考研] 一志愿大连理工大学，机械工程学硕，341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 求调剂 +4 图鉴212 2026-03-30 4/200 2026-03-31 10:20 by cal0306
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 262求调剂 +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
[考研] 085600材料与化工调剂 +16 kikiki7 2026-03-30 16/800 2026-03-31 10:03 by 氯化亚硝酰
[考研] 083000学硕274求调剂 +12 Li李鱼 2026-03-26 12/600 2026-03-31 10:01 by cal0306
[考研] 一志愿大连理工大学材料求调剂 +6 Gymno 2026-03-30 6/300 2026-03-31 07:26 by 无际的草原
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
[考研] 337求调剂 +6 《树》 2026-03-29 6/300 2026-03-30 10:15 by herarysara
[考研] 295求调剂 +5 wei-5 2026-03-26 5/250 2026-03-30 08:34 by 探123
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +7 AZMK 2026-03-28 11/550 2026-03-28 21:01 by xxxsssccc
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 一志愿厦门大学化学学硕307求调剂 +10 y7czhao 2026-03-26 10/500 2026-03-28 14:23 by 唐沐儿
信息提示
关闭
请填处理意见
关闭