求助吸附能计算中化学势的理解

莫雨萱 · 发表于 Post on 2022-3-25 09:39:19

在表面吸附模型中，计算吸附能公式是：Eads = Eadsslab − Eslab − μads，其中吸附原子化学势μads怎么理解和计算呀？

jiangning198511 · 发表于 Post on 2022-3-25 09:46:13

一般就是吸附分子的DFT总能

莫雨萱 · 发表于 Post on 2022-3-25 10:00:15

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

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

Eudaimonia · 发表于 Post on 2022-3-25 10:36:21

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

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

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

jiangning198511 · 发表于 Post on 2022-3-25 10:56:59

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

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

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

莫雨萱 · 发表于 Post on 2022-3-25 11:40:02

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

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

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

Eudaimonia · 发表于 Post on 2022-3-25 14:08:39

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

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

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

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

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

莫雨萱 · 发表于 Post on 2022-3-25 14:37:25

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合讨论] 求助吸附能计算中化学势的理解