求助H2分子在氧化铝负载镍原子表面上的吸附

MADAOLay

各位老师好，本人现在在测试H2分子在氧化铝负载镍原子表面上的吸附情况，但是结构优化之后H2分子总是直接解离为两个H原子吸附在表面上，请问这是正常的吗？
图中红色为氧原子，蓝色为铝原子，灰色为镍原子，白色为氢原子

喵星大佬

程序方法设置都不说嘛？

MADAOLay

喵星大佬 发表于 2020-12-17 19:44
程序方法设置都不说嘛？

不好意思，我用的是cp2k，用的是对角化计算方法。我就是想来问下氢气分子直接解离了是不是明显的计算错误，还是有这种可能性？

wzkchem5

MADAOLay 发表于 2020-12-17 19:48
不好意思，我用的是cp2k，用的是对角化计算方法。我就是想来问下氢气分子直接解离了是不是明显的计算错误 ...

泛函是什么？镍的磁矩怎么处理的？初始猜测结构是什么样的？

archer

https://www.sciencedirect.com/science/article/pii/S0360319920343366
善用文献搜索

MADAOLay

wzkchem5 发表于 2020-12-17 20:08
泛函是什么？镍的磁矩怎么处理的？初始猜测结构是什么样的？

不好意思，我刚刚入门，说的可能不太对，请大家担待。我用的泛函为PBE-D3，初猜结构第一种为氢气分子在Ni和O中间桥型吸附，第二种为在氧原子上TOP位点垂直吸附，最终结构氢气分子都解离了。关于磁矩处理，我好像没有考虑到这个，但是做了多重度为1，因为我测试过多重度为1时收敛最快。顺便问下老师“磁矩处理”一般还要考虑哪些呢？

MADAOLay

archer 发表于 2020-12-17 20:44
https://www.sciencedirect.com/science/article/pii/S0360319920343366
善用文献搜索

谢谢您

wzkchem5

MADAOLay 发表于 2020-12-17 20:57
不好意思，我刚刚入门，说的可能不太对，请大家担待。我用的泛函为PBE-D3，初猜结构第一种为氢气分子在Ni ...

设不同磁矩，取能量最低的即可。能量最低的收敛一般确实更快，但是这个关系不是必然的，所以还是要去对比能量。
初猜结构应该再试一种Ni的top位点吸附，你连O的top位点都算了，Ni的top位点没理由不算，因为只有Ni有能接受H2的sigma电子的空d轨道。

MADAOLay

wzkchem5 发表于 2020-12-17 21:16
设不同磁矩，取能量最低的即可。能量最低的收敛一般确实更快，但是这个关系不是必然的，所以还是要去对比 ...

好的，谢谢您，我算下Ni的top试试。