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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » DMol3 计算离子在金属表面的吸附能的问题
71/1返回列表
查看: 3217  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

bay__gulf

金虫 (著名写手)

刘苏州

[求助] DMol3 计算离子在金属表面的吸附能的问题 已有2人参与

离子放在周期性的Ag金属表面上, 有20 Ang 左右的真空层. 吸附能用总的能量分别减掉离子的能量和金属的能量
E(ion-Ag) - E(ion) - E(Ag)
金属Ag原子个数 5*5*4, 泛函PBE, 全原子和DSPP都试过, 限制自旋, 电荷在离子和离子-金属两个体系中+1, 只有金属的体系0.
做过多个正负离子的吸附能, 按照这样的方法都有问题:
遇到问题是, 在离表面很远的地方(10+ Ang) 的吸附能不为 0.
过程中其他参数包括金属原子的坐标都是不变的. 想问的是为什么会有这样大的偏差, 这种情况下吸附能如何算?
扫描的曲线大概是下面的样子
DMol3 计算离子在金属表面的吸附能的问题

PS 看了离子在10 Ang 左右位置的电荷(Mulliken/Hirshfeld), 发现对+1 价离子带电量~ +0.5 e, 对-1 价离子带电量~ -0.5 e.
而在计算吸附能E(ion-Ag) - E(ion) - E(Ag) 设置的ion 的电荷数+/- 1.0 e, Ag 金属表面不带电荷. 这个问题是否会导致这样的状况, 是否有必要设置在远离表面的位置上离子的带电量就是它的形式电量+- 1.0 e.
回复此楼

» 收录本帖的淘帖专辑推荐

MS计算模拟

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-09-06 13:46:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
不清楚原因，只是想问一下，会不会你那个原子并没有离表面10埃，不知道考虑边界条件没。
赞一下 回复此楼
2楼2015-09-06 15:14:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州
引用回帖:
2楼: Originally posted by xxxxxxxfc at 2015-09-06 15:14:40
不清楚原因，只是想问一下，会不会你那个原子并没有离表面10埃，不知道考虑边界条件没。

你好, 试过更大的真空层50 Ang, 原子离表面25 Ang, 这个位置依然有~100 kcal/mol 吸附能. 三个方向都是周期边界, 每个体系的边界大小都是一样的.
赞一下 回复此楼
3楼2015-09-06 16:44:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)
★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。。 2015-09-15 11:04:25
引用回帖:
3楼: Originally posted by bay__gulf at 2015-09-06 16:44:19
你好, 试过更大的真空层50 Ang, 原子离表面25 Ang, 这个位置依然有~100 kcal/mol 吸附能. 三个方向都是周期边界, 每个体系的边界大小都是一样的....

我意思是说虽然你的离子离界面有10埃，但该离子与其相邻的重复单元里面的原子距离有没有10埃？如果也满足条件的话，我就不太清楚了，可能是其它原因吧。
赞一下 回复此楼
4楼2015-09-06 16:52:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
bay__gulf: 金币+20, ★★★很有帮助 2015-09-07 09:06:36
franch: 金币+2, 谢谢回帖交流。。。。 2015-09-15 11:04:47
离子吸附？就是说体系带电是不是？
按道理讲这种计算是有问题的，详细请看
http://cms.mpi.univie.ac.at/vasp/vasp/FAQ.html
看最后几段
赞一下 回复此楼
5楼2015-09-06 20:42:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州
引用回帖:
5楼: Originally posted by luzihen at 2015-09-06 20:42:24
离子吸附？就是说体系带电是不是？
按道理讲这种计算是有问题的，详细请看
http://cms.mpi.univie.ac.at/vasp/vasp/FAQ.html
看最后几段

谢谢, 我去看看
回复此楼
6楼2015-09-07 09:07:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai1826369

银虫 (正式写手)
引用回帖:
6楼: Originally posted by bay__gulf at 2015-09-07 03:07:25
谢谢, 我去看看...

你好，打扰了，你的那个Hirshfeld电荷用dmol3怎么计算的啊？
赞一下(1人) 回复此楼
奋斗，奋斗，为了理想而奋斗
7楼2018-04-19 08:57:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bay__gulf 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭