锂硫电池台阶图自由能校正中锂离子以及小分子该如何处理？

yuzj41d

文献中Zhang X, Chen M, Zhang Q, et al. Applied Surface Science, 2021, 563: 150172.在催化剂表面 Li2S2 → Li2S转化的步骤如图所示，
我在进行类似的工作，请问图中这种自由能校正要如何进行呢 ？
可以看到中间有两步分别引入了一个锂离子，以及最后一步Li4S2裂解为两个Li2S，
这三步涉及原子数增减，对于锂离子，首先我是将 锂离子 进行了近似处理 ：Li+ + e- = Li ；
然后有三种处理思路：1 ，两个锂原子构成的单胞，能量为-3.99 eV 一个就是 -1.99 eV ；
2，将吸附构型的其余原子删去，只留一个锂进行优化，基本是一个15 x 15 x 15 A的一个原胞，能量是 -0.298 eV；
3， 优化出同样Slab 单独吸附一个锂原子 的结构，然后在能量变化时用 如：delat_G=G(Li-Li2S2@NMC)-G(Li2S2@NMC)-G(Li@NMC)+G(NMC)
也就是重复算一次Slab然后减去重复的能量。

我觉得这三种都有一定的道理但处理之后的数据截然不同，在此请教各位老师哪一个是更合理的呢？

第二个问题是：想请问这种吸附体系下，G(Li)与G(Li2S)如何校正呢？
Slab平板上的吸附构型我均采用刘锦成博士 先计算频率 后使用VASPKIT 501 校正，这个方法对上述小分子同样适用吗?
烦请各位老师解惑！

jiangning198511

我觉得是：
1.delat_G=G(Li-Li2S2@NMC)-G(Li2S2@NMC)-G(Li-BULK)
Li-BULK为Li-BCC结构归一化得能量
2.可以用

yuzj41d

jiangning198511 发表于 2022-7-5 10:21
我觉得是：
1.delat_G=G(Li-Li2S2@NMC)-G(Li2S2@NMC)-G(Li-BULK)
Li-BULK为Li-BCC结构归一化得能量

1，你好，请问一下 Li-BULK 为什么要是 BCC结构的呢 ？虽然我用的就是，但是Materials Project中还有很多其他结构的 ，为什么用BCC的呢 ？2，最后一步的能量也是用delat_G=G（Li2S@La2O3）+G(Li2S)-G(Li4S2)吗？其中G(Li2S)也是用单胞能量吗？ 我前面是将 Li2S在 15 x 15 x15 的格子里算的吸附能，这里是不是错了 。

jiangning198511

yuzj41d 发表于 2022-7-5 11:39
1，你好，请问一下 Li-BULK 为什么要是 BCC结构的呢 ？虽然我用的就是，但是Materials Project中还有很多 ...

1. 我记得是BCC为实验结构，所以选择BCC
2. 你可以重复一下文献的计算结果，看看是否能对的上，我这么也没什么经验

yinlina

我最近刚学习Li-S电池，关于吉布斯自由能的计算，想请教一下，请问您的G(Li)和G(Li2S)如何矫正的呢？