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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP加电场的问题
北京石油化工学院2026年研究生招生接收调剂公告
311/4返回列表1234下一页
查看: 5972  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 265求调剂11408 +4 刘小鹿lu 2026-03-27 4/200 2026-04-02 11:27 by guyan1000
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 266求调剂 +10 阳阳哇塞 2026-04-01 10/500 2026-04-02 09:07 by vgtyfty
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 求调剂推荐 材料 304 +21 荷包蛋hyj 2026-03-26 21/1050 2026-04-01 21:09 by lijunpoly
[考研] 320分，材料与化工专业，求调剂 +14 一定上岸aaa 2026-03-27 18/900 2026-04-01 20:10 by 积极调剂的小学�
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
[考研] 070300求调剂306分 +5 26要上岸 2026-03-27 5/250 2026-04-01 11:09 by oooqiao
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 086000调剂 +5 7901117076 2026-03-26 5/250 2026-03-31 17:45 by 544594351
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
[考研] 南京大学化学调剂 +11 景随风 2026-03-29 16/800 2026-03-31 10:14 by herarysara
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 085602化工求调剂（331分） +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
信息提示
关闭
请填处理意见
关闭