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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 请问单层石墨烯纳米带按照非对称slab静态计算时电子步不收敛是什么原因呢？
【调剂】北京石油化工学院2024年16个专业接受调剂
31/1返回列表
查看: 780  |  回复: 2
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者fengzhigu521将赠送您 50 个金币

fengzhigu521

至尊木虫 (著名写手)

[求助] 请问单层石墨烯纳米带按照非对称slab静态计算时电子步不收敛是什么原因呢？

1、构建单层石墨烯纳米带是为了算电催化反应的热力学反应路径，如果z轴方向真空层设置了20埃还需要把模型当成非对称的slab来算吗，还是直接按照对称的slab来算？
2、无论是原始的石墨烯纳米带还是加了吸附物种的石墨烯纳米带，静态计算时电子步均不收敛，我随后查了下文献，发现文献里基本都是“a vacuum space of 20 Å was introduced in the direction normal to the nanosheets in order to avoid virtual interactions between the periodic images.”或者“a vacuum layer of 20Å is chosen along the c axis direction normal to the sheer to avoid periodic interactions.”这种表述，也就是加了真空层直接忽略不对称slab中偶极的影响，请问这么处理的依据是什么？是因为电子步不收敛吗？
3、随后又查了一些二维范德华异质结的文章，也是基本按上述来表述的。
4、如果设置了偶极校正，一般会这么写“the periodic images of the slab are separated by 20 Å of vacuum and the dipole correction has been applied to the hamiltonian to decouple the electrostatic dipole interaction between periodically repeated slabs.”（h. a. hansen, v. viswanathan, j. k. nørskov, j. phys. chem. c 2014, 118, 6706.）

附：OSZICAR文件
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.145316675670E+04    0.14532E+04   -0.47310E+04  1876   0.744E+02
DAV:   2     0.778739111498E+02   -0.13753E+04   -0.12196E+04  2544   0.173E+02
DAV:   3    -0.342736648414E+03   -0.42061E+03   -0.41524E+03  2072   0.104E+02
DAV:   4    -0.366401712307E+03   -0.23665E+02   -0.23311E+02  2040   0.332E+01
DAV:   5    -0.366925498276E+03   -0.52379E+00   -0.52113E+00  2182   0.478E+00    0.432E+01
DAV:   6    -0.273381947789E+05   -0.26971E+05   -0.49736E+05  2366   0.479E+03    0.140E+02
DAV:   7     0.112889463599E+03    0.27451E+05   -0.25388E+04  2076   0.245E+02    0.407E+01
DAV:   8    -0.248287423359E+03   -0.36118E+03   -0.32647E+03  2090   0.123E+02    0.247E+01
DAV:   9    -0.296200066017E+03   -0.47913E+02   -0.19859E+02  2080   0.335E+01    0.231E+01
DAV:  10    -0.307085736739E+03   -0.10886E+02   -0.23894E+01  1984   0.137E+01    0.166E+01
DAV:  11    -0.320982407630E+03   -0.13897E+02   -0.93755E+00  1879   0.881E+00    0.133E+01
DAV:  12    -0.326774070237E+03   -0.57917E+01   -0.24821E+00  2046   0.431E+00    0.121E+01
DAV:  13    -0.339587376175E+03   -0.12813E+02   -0.13255E+01  1882   0.626E+00    0.127E+01
DAV:  14    -0.341173716380E+03   -0.15863E+01   -0.33645E+00  1832   0.266E+00    0.952E+00
DAV:  15    -0.349935296717E+03   -0.87616E+01   -0.48464E+00  1768   0.348E+00    0.948E+00
DAV:  16    -0.353821456579E+03   -0.38862E+01   -0.14643E+00  1896   0.283E+00    0.825E+00
DAV:  17    -0.359128912396E+03   -0.53075E+01   -0.19909E+00  1724   0.265E+00    0.736E+00
DAV:  18    -0.360159908992E+03   -0.10310E+01   -0.31701E-01  1872   0.127E+00    0.659E+00
DAV:  19    -0.115629214070E+05   -0.11203E+05   -0.21230E+05  2391   0.289E+03    0.105E+02
DAV:  20    -0.122382517069E+03    0.11441E+05   -0.14103E+04  1820   0.227E+02    0.559E+01
DAV:  21    -0.344572898196E+03   -0.22219E+03   -0.21981E+03  2060   0.102E+02    0.136E+01
DAV:  22    -0.361497692288E+03   -0.16925E+02   -0.13252E+02  1992   0.234E+01    0.385E+00
DAV:  23    -0.364093470793E+03   -0.25958E+01   -0.43705E+00  2016   0.470E+00    0.357E+00
DAV:  24    -0.364509436527E+03   -0.41597E+00   -0.28523E-01  2152   0.124E+00    0.351E+00
DAV:  25    -0.365029593444E+03   -0.52016E+00   -0.98409E-02  1784   0.938E-01    0.469E+00
DAV:  26    -0.364919622772E+03    0.10997E+00   -0.22244E-02  1768   0.428E-01    0.469E+00
DAV:  27    -0.364945954229E+03   -0.26331E-01   -0.16556E-03  2056   0.129E-01    0.450E+00
DAV:  28    -0.365355796530E+03   -0.40984E+00   -0.29778E-02  1783   0.523E-01    0.348E+00
DAV:  29    -0.366138780348E+03   -0.78298E+00   -0.89095E-02  1798   0.779E-01    0.249E+00
DAV:  30    -0.366796243577E+03   -0.65746E+00   -0.63014E-02  1816   0.646E-01    0.222E+00
DAV:  31    -0.164559743444E+05   -0.16089E+05   -0.13973E+05  2281   0.121E+03    0.132E+02
DAV:  32    -0.341021635291E+02    0.16422E+05   -0.19723E+04  2068   0.192E+02    0.418E+01
DAV:  33    -0.350280977509E+03   -0.31618E+03   -0.31099E+03  2035   0.111E+02    0.109E+01
DAV:  34    -0.366433388201E+03   -0.16152E+02   -0.15804E+02  1952   0.276E+01    0.307E+00
DAV:  35    -0.367153313637E+03   -0.71993E+00   -0.53935E+00  2168   0.510E+00    0.323E+00
DAV:  36    -0.367134365078E+03    0.18949E-01   -0.18102E-01  2000   0.994E-01    0.346E+00
DAV:  37    -0.367156565738E+03   -0.22201E-01   -0.24478E-02  1926   0.392E-01    0.352E+00
DAV:  38    -0.367150239336E+03    0.63264E-02   -0.18791E-02  1808   0.334E-01    0.381E+00
DAV:  39    -0.367158964805E+03   -0.87255E-02   -0.19224E-03  2032   0.114E-01    0.376E+00
DAV:  40    -0.367159327437E+03   -0.36263E-03   -0.17712E-04  2216   0.314E-02    0.378E+00
DAV:  41    -0.367163581925E+03   -0.42545E-02   -0.77818E-04  1767   0.938E-02    0.364E+00
DAV:  42    -0.361446107506E+05   -0.35777E+05   -0.67063E+05  2431   0.548E+03    0.144E+02
DAV:  43    -0.833591577116E+02    0.36061E+05   -0.25212E+04  2060   0.224E+02    0.419E+01
DAV:  44    -0.351369528165E+03   -0.26801E+03   -0.26209E+03  2053   0.104E+02    0.970E+00
DAV:  45    -0.371241410588E+03   -0.19872E+02   -0.14484E+03  1724   0.965E+01    0.222E+01
DAV:  46    -0.364303543222E+03    0.69379E+01   -0.50867E+02  1815   0.462E+01    0.562E+00
DAV:  47    -0.393129918352E+03   -0.28826E+02   -0.10909E+02  2221   0.240E+01    0.175E+01
DAV:  48    -0.363352052376E+03    0.29778E+02   -0.26349E+02  2277   0.247E+01    0.103E+01
DAV:  49    -0.366563901206E+03   -0.32118E+01   -0.43146E+01  1956   0.965E+00    0.352E+00
DAV:  50    -0.367331818244E+03   -0.76792E+00   -0.31748E+00  1971   0.400E+00    0.909E+00
DAV:  51    -0.366961099702E+03    0.37072E+00   -0.43633E-01  2016   0.172E+00    0.561E+00
DAV:  52    -0.367080146192E+03   -0.11905E+00   -0.14228E-01  1968   0.111E+00    0.407E+00
DAV:  53    -0.367134423836E+03   -0.54278E-01   -0.53102E-02  2015   0.673E-01    0.332E+00
DAV:  54    -0.178797232201E+05   -0.17513E+05   -0.41476E+05  2230   0.455E+03    0.118E+02
DAV:  55    -0.167470258655E+03    0.17712E+05   -0.18277E+04  1868   0.225E+02    0.232E+01
DAV:  56    -0.119414134652E+04   -0.10267E+04   -0.50988E+03  2315   0.106E+02    0.879E+01
DAV:  57    -0.240902185013E+03    0.95324E+03   -0.58050E+03  2708   0.116E+02    0.142E+01
DAV:  58    -0.357504505235E+03   -0.11660E+03   -0.13839E+03  1920   0.551E+01    0.409E+00
DAV:  59    -0.367008967359E+03   -0.95045E+01   -0.79327E+01  1996   0.169E+01    0.494E+00
DAV:  60    -0.367185562906E+03   -0.17660E+00   -0.38440E+00  2008   0.435E+00
   1 F= -.36718556E+03 E0= -.36717528E+03  d E =-.205598E-01  mag=     1.4547
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼2019-01-16 16:51:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

么么么哒165

禁虫 (小有名气)
本帖内容被屏蔽
2楼2019-01-17 00:21:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengzhigu521

至尊木虫 (著名写手)
引用回帖:
2楼: Originally posted by 么么么哒165 at 2019-01-17 00:21:35
试着解释一下电子步不收敛的问题，原因可能是1，步数设置的太少，2，这种结构不合理，能量无法趋于稳定，解决方法对应的是1，增加收敛步数，一般500步是上限，到500步依然不收敛，可以考虑微调一下吸附物质的放置位 ...

谢谢你的回复，原因刚才找到了，不过比较奇怪，是因为对于VASP来说，如果把2D模型放在晶胞中间的话，偶极校正就不收敛，具体来自于“Somehow VASP can't do dipole corrections if you put your slab in the center of the unit cell so try putting it at the bottom or at least a fair distance off center. ”
参见：https://cms.mpi.univie.ac.at/vasp-forum/viewtopic.php?t=13121
赞一下 回复此楼
3楼2019-01-17 09:13:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fengzhigu521 的主题更新
31/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金开始函评了吗？ +16 wych1103 2024-04-25 16/800 2024-04-26 21:32 by 淀粉搬运工
[考研] 没学上 +6 季向阳 2024-04-26 12/600 2024-04-26 21:06 by 季向阳
[考研] 0854－0855调剂 +8 shangannum1 2024-04-21 12/600 2024-04-26 16:42 by yz仔
[论文投稿] with editor 两个月了，什么原因？ +7 yiersan9 2024-04-24 15/750 2024-04-26 16:13 by jonewore
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[教师之家] 博士论文被抄袭 +25 和尚敲小木鱼 2024-04-22 42/2100 2024-04-26 13:55 by ZHONGWU_U
[考博] 25年博士申请 +6 Changzixuan 2024-04-25 11/550 2024-04-26 13:48 by 我属驴核动力驴
[硕博家园] 考研，求职还是考编? +15 xizj 2024-04-21 24/1200 2024-04-26 11:49 by Kan客
[考研] 381求调剂 +4 小刺猬987654321 2024-04-25 6/300 2024-04-26 10:57 by czl12138
[论文投稿] Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？ 10+4 yifeng11 2024-04-20 13/650 2024-04-26 09:48 by yifeng11
[教师之家] 刚博士毕业，入职普通一本高校 +19 保持肌肤 2024-04-20 20/1000 2024-04-26 09:16 by zhaojp918
[考博] 取博导收留 5+4 zzb777888 2024-04-20 10/500 2024-04-26 08:52 by polymerfriend
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[教师之家] 期末给学生划重点都是什么话术啊 +16 luokereng 2024-04-20 18/900 2024-04-25 15:46 by BusyGer
[考博] 求博导 +6 好okjh 2024-04-21 10/500 2024-04-25 14:04 by 好okjh
[博后之家] 南京大学-广州大学联合招聘博士后 欢迎广大优秀人才！！！ +4 黑魔变身啾 2024-04-20 12/600 2024-04-25 11:18 by dodonaomi
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿，黑贝
[教师之家] 大家访学都是怎么找的啊？ +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
信息提示
关闭
请填处理意见
关闭