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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助！vasp加一个负点体系的真空静电势取值
172/2返回列表上一页12
查看: 2898  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by 罪歌白夜叉 at 2016-07-01 16:19:48
是的，vasp手册里也有提到，the total energy depends linearly on the width of the vacuum, 我是参考Matthew Neurock的方法模拟电化学反应，绝对能量没有什么意义，体系的能量差值应该是有意义的...

说起Matthew Neurock的文章，他的文章中的体系能量确实经过校正，请注意J. Phys. Chem. B, Vol. 110, No. 43, 2006的Eq. 13，该文章当中其实提了几种方法吧，不过我总感觉不及刘智攀做法有特色。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
11楼2016-07-01 17:32:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

罪歌白夜叉

铁虫 (初入文坛)
引用回帖:
11楼: Originally posted by 卡开发发 at 2016-07-01 17:32:02
说起Matthew Neurock的文章，他的文章中的体系能量确实经过校正，请注意J. Phys. Chem. B, Vol. 110, No. 43, 2006的Eq. 13，该文章当中其实提了几种方法吧，不过我总感觉不及刘智攀做法有特色。...

是的，我也是打算在总能量基础上减去背景电荷作用的能量做修正，这种校正是不是能解决能量发散的问题我不也不太确定。。刘志攀他们自己写的隐式溶剂模型，包含了对带电体系的处理，我现在是用的Hennig写的VASPsol做溶剂效应，应该和刘的类似吧
赞一下 回复此楼
12楼2016-07-02 10:26:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

罪歌白夜叉

铁虫 (初入文坛)
引用回帖:
11楼: Originally posted by 卡开发发 at 2016-07-01 17:32:02
说起Matthew Neurock的文章，他的文章中的体系能量确实经过校正，请注意J. Phys. Chem. B, Vol. 110, No. 43, 2006的Eq. 13，该文章当中其实提了几种方法吧，不过我总感觉不及刘智攀做法有特色。...

另外，我后来把真空层扩大成30埃之后，作了一下静电势图，和想象的不一样，看不出两边表面附近的静电势有明显的区别，加了偶极修正之后（IDIPOL=3），也没有变化，vasp手册里面说对于非立方体系无法做Quadrupole校正，不知道是不是这个原因？（第一个图是偶极修正之前，第二个图是偶极修正之后，体系带一个正电）可以把20附近的值作为真空静电势值吗？
求助！vasp加一个负点体系的真空静电势取值
6.png


求助！vasp加一个负点体系的真空静电势取值-1
7.png
赞一下 回复此楼
13楼2016-07-02 12:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
12楼: Originally posted by 罪歌白夜叉 at 2016-07-02 10:26:55
是的，我也是打算在总能量基础上减去背景电荷作用的能量做修正，这种校正是不是能解决能量发散的问题我不也不太确定。。刘志攀他们自己写的隐式溶剂模型，包含了对带电体系的处理，我现在是用的Hennig写的VASPsol做 ...

我觉得如果用了背景电荷、偶极修正和能量修正应该可以吧。

刘智攀他们自己在SIESTA里面加了个修正的modified Poisson-Boltzmann模型，好像过渡态搜索也是自己写的吧。我感觉应该是通过控制等效的点电荷来实现的，不过真去改SIESTA也够麻烦的。

VASPsol应该类似吧，但刘智攀他们的那个介电函数的处理好像用的是Fattebert-Gygi形式的吧，我看Hennig的文章测试结果两者好像差不多。这种形式QE的Environ插件支持http://www.quantum-environment.org/，并且这个插件支持直接使用类似于MP校正的方法计算电荷体系。所以正准备试一下。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
14楼2016-07-02 13:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
13楼: Originally posted by 罪歌白夜叉 at 2016-07-02 12:23:25
另外，我后来把真空层扩大成30埃之后，作了一下静电势图，和想象的不一样，看不出两边表面附近的静电势有明显的区别，加了偶极修正之后（IDIPOL=3），也没有变化，vasp手册里面说对于非立方体系无法做Quadrupole校 ...

看起来确实似乎左右两侧的静电势差距不是很大。否则20~25这一段不修正会是个斜线。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
15楼2016-07-02 13:54:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

罪歌白夜叉

铁虫 (初入文坛)
引用回帖:
15楼: Originally posted by 卡开发发 at 2016-07-02 13:54:32
看起来确实似乎左右两侧的静电势差距不是很大。否则20~25这一段不修正会是个斜线。...

嗯嗯，谢谢前辈的指导
赞一下 回复此楼
16楼2016-07-02 21:01:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

专心学术

新虫 (初入文坛)
引用回帖:
10楼: Originally posted by 卡开发发 at 2016-07-01 17:28:29
前辈说的在理，实际我觉得是否能这么理解，电荷的发散可以通过引入背景电荷来消除，电荷体系的能量随着真空高度发散应该可以通过偶极矩修正来进行消除。其实应该还有个”电荷浓度“对体系能量的影响”也得进行消除 ...

您好，我算的是银-二氧化硅的界面体系的静电势图，加了偶极校正，真空层也足够厚（40A），可是算出来的界面两侧的真空能级不同，这点我可以理解，可是两侧还都是倾斜的，大概是什么问题，我该如何解决，我知道我的界面体系并不是电中性的。希望您有时间回复一下。
赞一下 回复此楼
17楼2018-12-20 22:58:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 罪歌白夜叉 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中科技大学，080502，354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 286求调剂 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭