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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP加电场的问题
北京石油化工学院2026年研究生招生接收调剂公告
311/4返回列表1234下一页
查看: 5975  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +4 末未mm 2026-04-02 5/250 2026-04-03 14:20 by 王亮_大连医科大
[考研] 333求调剂 +9 wfh030413@ 2026-04-03 9/450 2026-04-03 14:20 by 啊俊！
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 071000生物学调剂 +8 知昭蔓 2026-04-02 8/400 2026-04-03 10:36 by macy2011
[考研] 294求调剂 +3 Grey_Ey 2026-04-02 5/250 2026-04-02 21:34 by swb0218
[考研] 0832食品科学与工程学硕282调剂 +4 鱼在水中游a 2026-04-02 7/350 2026-04-02 14:12 by baoball
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 08开头看过来！！！ +4 wwwwffffff 2026-03-31 6/300 2026-04-02 11:42 by 均值回归
[考研] 311求调剂 +9 勇敢的小吴 2026-04-02 9/450 2026-04-02 11:37 by Sammy2
[考研] 285求调剂 +11 AZMK 2026-04-01 11/550 2026-04-01 22:40 by peike
[考研] 一志愿南昌大学324求调剂 +12 hanamiko 2026-03-27 12/600 2026-04-01 22:34 by 运气yunqi
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 284求调剂 +12 小熊～～ 2026-03-31 12/600 2026-04-01 20:23 by 花??
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 调剂考研 +3 王杰一 2026-03-29 3/150 2026-03-29 08:09 by fmesaito
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
信息提示
关闭
请填处理意见
关闭