计算化学公社

标题: 求助使用VASP计算CO2、H2在Ag（110）上的吸附时，结构优化完分子距离slab表面很远 [打印本页]

作者
Author: calamity 时间: 2020-6-3 18:45
标题: 求助使用VASP计算CO2、H2在Ag（110）上的吸附时，结构优化完分子距离slab表面很远
如题，使用如图1参数，1 1 1的k点优化附件所示POSCAR结构时，初始C原子距离表面1.95A左右，然后优化完成后达到了3A左右已经超过了化学吸附的范畴..这个INCAR是在筛选可能的吸附结构，但是我总共弄了6个结构，都优化后到了3A左右的距离，包括H2在Ag（110）的吸附也出现了这样的问题。请问各位这是什么原因引起的？是泛函选择问题？或是其他输入参数问题？还是初始结构存在问题？

作者
Author: hakuna 时间: 2020-6-3 23:10
1.k-mesh太差
2.收敛精度太低
3.考虑色散矫正

作者
Author: calamity 时间: 2020-6-4 13:39

hakuna 发表于 2020-6-3 23:10
1.k-mesh太差
2.收敛精度太低
3.考虑色散矫正

谢谢！我马上试试

作者
Author: calamity 时间: 2020-6-4 21:17

hakuna 发表于 2020-6-3 23:10
1.k-mesh太差
2.收敛精度太低
3.考虑色散矫正

你好，我之前尝试了修改k点密度、收敛精度并加入了DFT-D3修正，但是还是出现吸附分子远离slab表面，难道这是我slab模型的问题？有没有可能就是这样的模型不能使得CO2和H2产生化学吸附

作者
Author: hakuna 时间: 2020-6-4 23:24
CO2、H2本身很稳定，吸附本来很弱，Ag上会更弱。CO2吸附时，本身结构也会变化，某些情况下呈V字形，垂直或平行吸附，也会存在mock state，一种距离上判断是吸附的，但吸附能看则不利的吸附结构。CO2研究很多，多看看文献。

作者
Author: calamity 时间: 2020-6-5 19:45

hakuna 发表于 2020-6-4 23:24
CO2、H2本身很稳定，吸附本来很弱，Ag上会更弱。CO2吸附时，本身结构也会变化，某些情况下呈V字形，垂直或 ...

受教了，不过我现在也是在复现一篇CO2RR文献里的工作，他里面也提到了CO2吸附在表面键角发生变化，但是目前由于模型连吸附都优化不出来，无法向下继续计算，请问我应该怎么办？

作者
Author: hakuna 时间: 2020-6-8 16:10
看看别人怎么摆放CO2的，研究一下不同构型CO2的轨道空间趋向，怎么摆放才和表面原子轨道匹配.......

作者
Author: aaronzjw 时间: 2020-6-8 16:21
carbon dioxide 本来就是若吸附，远离表面正常

作者
Author: reductioncat. 时间: 2020-6-11 03:14
同学建议你先算算吸附能从你的POSCAR看 Ag-O键长2.3 已经不像是物理吸附了。二氧化碳线性分子吸附结构不是非常多的。你可以试试End-On 的构型。

作者
Author: calamity 时间: 2020-6-26 16:27

reductioncat. 发表于 2020-6-11 03:14
同学建议你先算算吸附能从你的POSCAR看 Ag-O键长2.3 已经不像是物理吸附了。二氧化碳线性分子吸附结构 ...

好的，我试一下，谢谢

作者
Author: calamity 时间: 2020-6-26 16:28

aaronzjw 发表于 2020-6-8 16:21
carbon dioxide 本来就是若吸附，远离表面正常

之前联系了一个文献作者，认为Ag表面CO2物理吸附很正常，化学吸附很难很难，看来应该也是正常情况。不过也有看到有人优化得到了化学吸附

作者
Author: calamity 时间: 2020-6-26 16:28

hakuna 发表于 2020-6-8 16:10
看看别人怎么摆放CO2的，研究一下不同构型CO2的轨道空间趋向，怎么摆放才和表面原子轨道匹配.......

好的，我后来找到了一篇相关文献，正在尝试，谢谢！

作者
Author: 将军金甲夜不脱 时间: 2020-6-26 17:52
楼主，想问下你计算时，前边已经用1 1 1 K点优化了结构，后来看你又修改了K点，想问下你是怎么修改的？是重新算了吗还是在111基础上变为更大的k点？

作者
Author: calamity 时间: 2020-6-26 21:22

将军金甲夜不脱发表于 2020-6-26 17:52
楼主，想问下你计算时，前边已经用1 1 1 K点优化了结构，后来看你又修改了K点，想问下你是怎么修改的？是重 ...

在111k点优化后的结构上进行更高k点的优化的，就相当于111k点进行低精度的快速优化然后再用更密的k点进行高精度优化

作者
Author: archer 时间: 2020-7-1 20:00
H2通常只能在表面发生弱吸附，你发现距离远很正常。对于CO2，应该有一个物理吸附和化学吸附构型。现在你只找到了物理吸附点，对于化学吸附应当是，CO2呈现V形，O朝下接近金属表面。

作者
Author: aaronzjw 时间: 2020-7-2 11:32

calamity 发表于 2020-6-26 16:28
之前联系了一个文献作者，认为Ag表面CO2物理吸附很正常，化学吸附很难很难，看来应该也是正常情况。不过 ...

那说明人家的体系就是化学吸附，你试试给体系加电子，说不定也能实现

作者
Author: calamity 时间: 2020-7-8 10:08

archer 发表于 2020-7-1 20:00
H2通常只能在表面发生弱吸附，你发现距离远很正常。对于CO2，应该有一个物理吸附和化学吸附构型。现在你只 ...

谢谢！

作者
Author: Xu98H 时间: 2020-8-17 08:45
您好，关于CO2在Ag上的吸附，我也遇到了相同的问题，方便请教一下吗？

作者
Author: 崔振莹 时间: 2020-10-29 19:20
你好，请问问题解决了吗，我现在在优化H2在Fe表面吸附的时候也遇到了同样的问题，优化之后都是远离表面的，请问方便交流一下吗，谢谢

作者
Author: calamity 时间: 2020-11-20 21:57

Xu98H 发表于 2020-8-17 08:45
您好，关于CO2在Ag上的吸附，我也遇到了相同的问题，方便请教一下吗？

确实不容易化学吸附，如果要化学吸附的话可以找文献看看怎么放

作者
Author: calamity 时间: 2020-11-20 21:57

崔振莹发表于 2020-10-29 19:20
你好，请问问题解决了吗，我现在在优化H2在Fe表面吸附的时候也遇到了同样的问题，优化之后都是远离表面的， ...

还是那样，就物理吸附的算，你的情况可能差不多，H2很难吸附在Fe表面，不一定是错误的

作者
Author: Yaohaohao 时间: 2024-9-25 20:49

aaronzjw 发表于 2020-7-2 11:32
那说明人家的体系就是化学吸附，你试试给体系加电子，说不定也能实现

请问使用NELECT命令来给体系加电子吗

欢迎光临计算化学公社 (http://bbs.keinsci.com/)