计算化学公社

标题: 求助吸附能计算中化学势的理解 [打印本页]

作者
Author: 莫雨萱 时间: 2022-3-25 09:39
标题: 求助吸附能计算中化学势的理解
在表面吸附模型中，计算吸附能公式是：Eads = Eadsslab − Eslab − μads，其中吸附原子化学势μads怎么理解和计算呀？

作者
Author: jiangning198511 时间: 2022-3-25 09:46
一般就是吸附分子的DFT总能

作者
Author: 莫雨萱 时间: 2022-3-25 10:00

jiangning198511 发表于 2022-3-25 09:46
一般就是吸附分子的DFT总能

因为我现在需要计算单个Fe的化学势，我目前的做法是：直接计算一个Fe单胞的总能量，然后再除以原子个数，不知道这个是否可行呀？

作者
Author: Eudaimonia 时间: 2022-3-25 10:36

莫雨萱发表于 2022-3-25 10:00
因为我现在需要计算单个Fe的化学势，我目前的做法是：直接计算一个Fe单胞的总能量，然后再除以原子个数， ...

两种算法，你这是一种，另一种使用真空中单个Fe原子的总能量

一般会定义为不同的能量，Eformation和Ebinding

作者
Author: jiangning198511 时间: 2022-3-25 10:56

Eudaimonia 发表于 2022-3-25 10:36
两种算法，你这是一种，另一种使用真空中单个Fe原子的总能量

一般会定义为不同的能量，Eformation和Eb ...

这两种好像文章中都有涉及，一般催化中吸附分子用分子的能量，而如果涉及单原子催化中单原子吸附时一般用原子晶体的总能(归一化到单个原子)，在文章中会详细说明的

作者
Author: 莫雨萱 时间: 2022-3-25 11:40

Eudaimonia 发表于 2022-3-25 10:36
两种算法，你这是一种，另一种使用真空中单个Fe原子的总能量

一般会定义为不同的能量，Eformation和Eb ...

我之前用同样的参数设置分别用这两种方法去计算过，但是发现他们计算出来的差值挺大的，最后带入吸附能公式中发现对吸附能影响也很大，就和文献上的数值吻合得很不好，所有不知道该怎么去计算这个化学势了

作者
Author: Eudaimonia 时间: 2022-3-25 14:08

莫雨萱发表于 2022-3-25 11:40
我之前用同样的参数设置分别用这两种方法去计算过，但是发现他们计算出来的差值挺大的，最后带入吸附能公 ...

差很大是很正常的，因为气相中单个金属原子不稳定或者说不应该存在，导致用这个值算出来的E一定会是负的，印象里一般会在-5 到 -10eV左右

如果你和自己比较，只要保证使用的化学势的计算方法一致即可

和文献上数值吻合不好的原因有很多，各种小细节是否设置都统一了，使用的程序是否是同一个甚至同一版本，如果出入很大甚至趋势都不同建议问一下作者

此外，我在前面说这个化学势的计算方法用的是Eformation或者EBinding，和Eads的算法不太一致，你主楼中的uads就会显得比较奇怪。因为吸附时，一般我们使用的吸附物都是反应物/生成物/中间体，存在单个原子吸附的状况比较少

作者
Author: 莫雨萱 时间: 2022-3-25 14:37

Eudaimonia 发表于 2022-3-25 14:08
差很大是很正常的，因为气相中单个金属原子不稳定或者说不应该存在，导致用这个值算出来的E一定会是负的 ...

太感谢您的回答啦，我数值的趋势与文章是一致的，只是在数值上差的比较大，我和参考文献中使用的软件也不是同一个软件，我接下来再去看看相关文章，非常感谢呀

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