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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP加电场的问题
北京石油化工学院2026年研究生招生接收调剂公告
311/4返回列表1234下一页
查看: 5966  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600，材料与化工321分，求调剂 +11 大馋小子 2026-03-27 11/550 2026-04-01 16:10 by cymywx
[考研] 调剂 +5 好好读书。 2026-03-28 7/350 2026-04-01 15:32 by 王亮_大连医科大
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 085601材料工程找调剂 +18 oatmealR 2026-03-29 19/950 2026-04-01 14:36 by chenqifeng666
[考研] 311求调剂一志愿合肥工业大学 +14 秋二十二 2026-03-30 14/700 2026-04-01 11:45 by chemdavid
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
[考研] 271求调剂 +11 勒布朗@ 2026-03-31 13/650 2026-04-01 11:18 by 逆水乘风
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 一志愿南师大0703化学 275求调剂 +5 Ripcord上岸 2026-03-27 5/250 2026-03-31 19:52 by mg1014
[考研] 一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂! +9 真名有冰 2026-03-29 10/500 2026-03-31 18:01 by xhai2011
[考研] 329求调剂，一志愿西北工业大学，材料工程（085601） +6 小小机灵虫 2026-03-29 12/600 2026-03-31 16:58 by 记事本2026
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 085601 材料工程 313分 求调剂 +6 Ong3 2026-03-27 6/300 2026-03-31 15:18 by yanflower7133
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +13 辰沐5211314 2026-03-26 13/650 2026-03-31 14:37 by 记事本2026
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 16:05 by shuang5186
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 070300化学求调剂 +4 起个名咋这么难 2026-03-27 4/200 2026-03-27 21:39 by 83503孙老师
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
信息提示
关闭
请填处理意见
关闭