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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP加电场的问题
311/4返回列表1234下一页
查看: 6444  |  回复: 30
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

whywhy01

金虫 (小有名气)

[求助] VASP加电场的问题 已有1人参与

在金属薄膜体系加电场，有15A的真空层，电场沿Z轴方向
开电场用的下面三个参数：
  EFIELD = 0.05
  IDIPOL = 3
  LDIPOL = .TRUE.
自洽计算电荷密度时发现只有LDIPOL的开关对计算过程和结果有影响，电场的正负和大小都不能改变CHGCAR和WAVECAR，感觉这个结果很不符合常理
这是完整的INCAR：
    PREC = Accurate
   ENCUT = 500
  IBRION = -1
  NELMIN = 2
   EDIFF = 1.0e-05
VOSKOWN = 1
  NBLOCK = 1
    NELM = 60
    ALGO = Normal (blocked Davidson)
   ISPIN = 2
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
ADDGRID = .FALSE.
  ISMEAR = -5
   SIGMA = 0.2
   LREAL = A
  LORBIT = 11
  MAGMOM = 4*0.749 4*0.813 4*1.974
LMAXMIX = 4
    NPAR = 6
  EFIELD = 0.3
  LDIPOL = .TRUE.
  IDIPOL = 3

其他参数应该都没有问题。希望大家能帮我看一看问题出在哪儿。。。
回复此楼

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

DFT加电场 VASP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-11-28 15:22:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-11-30 17:14:04
你怎么判断”电场的正负和大小都不能改变CHGCAR和WAVECAR，“。
你将静电势输出到LOCPOT，然后对LOCPOT处理一下，看加电场后真空区间的势的变化，是不是斜的。
赞一下(1人) 回复此楼
2楼2014-11-29 19:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
如果你的POSCAR是通过ISIF=3优化出来的话，那你的这个POSCAR就有点问题了。
对于slab的超原胞，除非你修改过vasp源代码，用ISIF=3来优化这样体系的结构是有问题的，它会调整真空层的厚度。
赞一下(1人) 回复此楼
11楼2014-11-30 20:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)
引用回帖:
19楼: Originally posted by valenhou001 at 2014-12-04 10:24:09
加减电子数，还不是模拟加静电场。 通过偶极矩、锯齿型势、平行板电容器加静电场会改变上下两个表面上的电荷分布，但是不改变体系的总电荷。

你试过把电场强度增大一些，比如0.2, 0.3或更大一点。...

侯老师，我把金属层加厚到了6层原子，发现电场能比较明显地起作用了。现在看可能是太薄的金属层会限制电场的效果。
赞一下(1人) 回复此楼
29楼2014-12-05 17:35:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

whywhy01

金虫 (小有名气)
引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-29 19:06:10
你怎么判断”电场的正负和大小都不能改变CHGCAR和WAVECAR，“。
你将静电势输出到LOCPOT，然后对LOCPOT处理一下，看加电场后真空区间的势的变化，是不是斜的。

我简单的看了一下CHGCAR和WAVECAR的大小，发现是一样的，可能这个方法不够准确。我试一下用LOCPOT处理的方法，谢谢！
赞一下 回复此楼
3楼2014-11-30 14:33:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)
引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-29 19:06:10
你怎么判断”电场的正负和大小都不能改变CHGCAR和WAVECAR，“。
你将静电势输出到LOCPOT，然后对LOCPOT处理一下，看加电场后真空区间的势的变化，是不是斜的。

是用vtotav处理LOCPOT出来的VLINE文件看吗？加正负方向电场的曲线是完全重合的
赞一下 回复此楼
4楼2014-11-30 15:21:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
引用回帖:
4楼: Originally posted by whywhy01 at 2014-11-30 16:21:11
是用vtotav处理LOCPOT出来的VLINE文件看吗？加正负方向电场的曲线是完全重合的...

不可理解你的结果。加的电场强度大小是多少？
赞一下 回复此楼
5楼2014-11-30 15:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-11-30 17:14:16
真空层在哪个方向上？
赞一下 回复此楼
6楼2014-11-30 15:34:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)
引用回帖:
5楼: Originally posted by valenhou001 at 2014-11-30 15:32:30
不可理解你的结果。加的电场强度大小是多少？...

大小是0.1 V/A，加在z轴上。运行vtota.exe是这样的：
Which direction to keep? (1-3 --- 1=X,2=Y,3=Z)
3
          12
positions read
charge density read
  2.79017870E-04
然后拿输出的VLINE直接放origin里面画图了，请问这里有没有问题？

从OUTCAR中能看到EFIELD=0.1或-0.1，程序识别了电场，但是没有明显的作用
赞一下 回复此楼
7楼2014-11-30 16:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by whywhy01 at 2014-11-30 17:23:03
大小是0.1 V/A，加在z轴上。运行vtota.exe是这样的：
Which direction to keep? (1-3 --- 1=X,2=Y,3=Z)
3
          12
positions read
charge density read
  2.79017870E-04
然后拿输出的VLINE直接放o ...

这一步的步骤没有问题。基本正确。
贴出你的POSCAR。
赞一下 回复此楼
8楼2014-11-30 20:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)
引用回帖:
8楼: Originally posted by valenhou001 at 2014-11-30 20:15:28
这一步的步骤没有问题。基本正确。
贴出你的POSCAR。...

New structure                           
   1.00000000000000     
     5.6733630828874579    0.0000000000000000    0.0000000000000000
     0.0000000000000000    4.9119752031931467    0.0013815384033351
     0.0000000000000000   -0.0028505105105781   25.1157632295913587
   Ni   Co
    8     4
Direct
  0.7500000000000000  0.6666457350900428  0.3975318092546394
  0.2500000000000000  0.6666457350900428  0.3975318092546394
  0.5000000000000000  0.1666457490900439  0.3975318092546394
  0.0000000000000000  0.1666457490900439  0.3975318092546394
  0.2500000000000000  0.3340795336386695  0.4727695484013834
  0.7500000000000000  0.3340795336386695  0.4727695484013834
  0.5000000000000000  0.8340795626386611  0.4727695484013834
  0.0000000000000000  0.8340795626386611  0.4727695484013834
  0.2500000000000000  0.9996329894720282  0.5476429281979919
  0.7500000000000000  0.9996329894720282  0.5476429281979919
  0.5000000000000000  0.4996329894720282  0.5476429281979919
  0.0000000000000000  0.4996329894720282  0.5476429281979919

这里我把薄膜放在了中间，会有影响吗？
赞一下 回复此楼
9楼2014-11-30 20:31:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
9楼: Originally posted by whywhy01 at 2014-11-30 21:31:32
New structure                           
   1.00000000000000     
     5.6733630828874579    0.0000000000000000    0.0000000000000000
     0.0000000000000000    4.9119752031931467    0.0013815384 ...

原子层放在中间，对结果没有影响。

为啥你的基矢a23, a32不是为零？　你计算的什么类型的表面（表面方向和表面大小）？
难道你的POSCAR是通过ISIF=3优化出来的吗？
赞一下 回复此楼
10楼2014-11-30 20:53:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whywhy01 的主题更新
311/4返回列表1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:29 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭