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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助！vasp加一个负点体系的真空静电势取值
172/2返回列表上一页12
查看: 2806  |  回复: 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
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭