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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 急救：用VASP算电偶极矩出问题
81/1返回列表
查看: 3425  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +47 AI8RaGaPaSCR 2026-08-07 50/2500 2026-08-14 02:15 by QhTuz5gvDira
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
关闭
请填处理意见
关闭