析氧反应（OER）是电解水中的重要反应，而开发低成本、高活性和稳定性的OER催化剂仍极具挑战性。近日，中国科学院孙公权、王素力等人 使用密度泛函理论（DFT）方法计算了具有不同结构（MN4 C8 、MN4 C10 、MN4 C12 ）的金属氮碳（MNC，M=Co，Ru，Rh，Pd，Ir）基电催化剂的OER活性和稳定性。 作者使用 Dmol3 程序包 进行DFT计算，并通过Grimme或TS方法进行自旋极化和长距离色散校正。 在所有模拟中，作者使用类导体屏蔽模型（COSMO）来模拟水溶液，并将介电常数设置为78.54。在AIMD计算中，作者采用NVE系综，并使用1 fs、2 fs和5 fs的时间步长，以及将温度设置为300 K、350 K、400 K、450 K和500 K。此外，作者应用局部电场来模拟双电层（EDL）。
结果与讨论
ΔG*OH 是筛选OER电催化剂的重要参数，即ΔG*OH 可以直接评估*OH在可逆电极电位下的吸附/解吸能力。而MNC催化剂中M−N4 的局部配体和几何结构对其电子和几何性质以及OER中间体（*OH，*O，*OOH）的吸附能都很重要。 因此，作者计算了MN4 C8 、MN4 C10 和MN4 C12 （M=Co、Ru、Rh、Pd、Ir）电催化剂的ΔG*OH ，具体如图1 所示。 如图2a所示，作者将ΔG*OH 分为三组，即ΔG*OH >1.53 eV（PdN4 C8 ，PdN4 C10 ，PdN3 C12 ），ΔG*OH <1.23eV（RuN4 C8 、RuN4 C10 、RuN3 C12 、CoN4 C8 和CoN4 C10 ）和1.23eV<ΔG*OH <1.53eV（RhN4 C8 ，RhN4 C10 ，RhN3 C12 ，IrN4 C8 ），对于ΔG*OH >1.53eV，ΔG1 通常是RDS。对于ΔG*OH <1.23eV，在OER工作环境下，产物的脱附变得困难，即过渡金属原子将被氧化。 过渡金属原子和氮掺杂碳基底之间的结合能是耐久性的常用描述符。如图2b所示，具有不同局部结构的MNC的结合能顺序为：MN4 C8 <MN4 C10 <MN4 C12 。M原子的第二近邻碳原子的数量对MNC催化剂的稳定性起决定性作用。 除了微观结构，金属原子是影响MNC稳定性的另一个重要因素。 对于MN 4 C 8 和MN 4 C 10 ，结合能按顺序为： Pd<Ir<Rh<Ru<Co。 然而，在MN 4 C 12 体系中，RuNC比CoNC更稳定。 总之，RuN 4 C 12 具有最高的结合能，而PdN 4 C 8 在真空环境下是最不稳定的。 作者进一步分析了RhN4 C8 、RhN4 C10 、RhN4 C12 、IrN4 C8 ，IrN4 C10 ，IrN4 C12 和CoN4 C12 （1.23eV<ΔG<1.53eV）上*OH、*O和*OOH的吸附结构以及OER自由能图（图3a），发现MNC电催化剂的RDS不同，并且RhN4 C8 具有最高的OER活性，相应MNC电催化剂的过电位如图3b所示。 对于电化学反应，三相界面中的局部电场是影响电催化剂耐久性的重要参数（图4），作者详细研究了电场强度对电催化剂的影响。对于纯MNC，ΔEb 在各种电场下几乎保持不变，然而，当MNC吸附*OOH、*OH和*O时，情况变得不同。在正电场（阴极条件）下，ΔEb 将随着电场强度的增加而减小，而在负电场（阳极条件）下观察到相反的趋势。 此外，作者还计算了温度的影响，即在300K下用*O吸附MN4 C12 进一步与电场耦合，即瞬时最小ΔEb 变得更低。作者还考虑了在各种环境下用2*OH吸附的MN4 C12 的稳定性。 最小瞬时ΔEb 可能低至−2.68eV，而稳定性趋势分别为CoN4 C12 （-3.44eV）<RhN4 C12 （-3.61eV）<IrN4 C12 (-3.84 eV) 和CoN4 C12 (-1.81 eV) < IrN4 C12 (-2.51 eV) < RhN4 C12 (-2.68 eV)。 如图5a所示。对于不同的RuNC电催化剂，结构演变经历了不同的亚稳态/中间态结构，包括*OH-RuNC、2*OH-RuNC，*O&*OH-RuNC、2*O-RuNC和最终的*O-RuNC。在该过程中， OH*+2O*→ OOH*+O*步骤是RuN 4 C 8 （2.05eV）、RuN 4 C 10 （1.92eV）和RuN 4 C 12 （2.09eV）的RDS。 作者进一步计算了 RuNC和*O-RuNC的OER活性（图5b）。 然而，在结构演化之前和之后，过电势显著不同。 因此，工作条件效应导致的结构演变对OER活性有显著影响。 与RhNC电催化剂类似，RuON4 C12 的ΔEb 在正电场下随着电场强度的增加而减小，而在负电场下随着*OOH、*OH和*O的吸附而增加。然而，RuO NC的稳定性变化比RhNC更显著。 当被*O原子吸附时，RhN4 C12 的ΔΔEb （ΔEb的变化量）从−0.40 V/Å到0.40 V/Å仅为0.70 eV，远低于RuO N4 C12 的7.74 eV。将最小瞬时ΔEb （300 K动态，0.40 V/Å）与*O吸附进行比较，稳定性趋势为RuON4 C12 （-1.55 eV）<CoN4 C12 （-3.44 eV） < RhN4 C12 (-3.61 eV) < IrN4 C12 (-3.84 eV) (图6)。 RuON4 C12 在工作条件下的ΔEb 值（ΔEb =-1.55eV）远低于其固有稳定性（ΔEb =-8.84eV），在水电解过程中其结构演变对其有显著影响。 对于所有PdNC的OER活性，它们的RDS相同，因为它们对*OH的吸附能低（图7a）。PdNC稳定性在电场环境中的变化趋势与其他MNC相似。然而，PdN4 C12 的2*OH的ΔEb 变化率比O*更显著。 当*O原子吸附时，PdN4 C12 的结合能从-0.40V/Å到0.40V/Å的变化仅为0.77eV，这远低于2*OH-PdN4 C12 （图7b）的6.53eV。PdN4 C12 在工作条件下的ΔEb 值（ΔEb =1.99eV）低于其固有稳定性（ΔEb =-6.09eV），这是因为其受到水电解过程中中间产物吸附的影响。 因此，稳定性趋势为RuO N4 C12 （-1.55 eV）<PdN4 C12 （-1.99 eV）<CoN4 C12 (-3.44 eV) < RhN4 C12 (-3.61 eV) < IrN4 C12 (-3.84 eV)（图7c）。 如图8所示，ΔG*OH >1.53eV或ΔG*OH <1.23eV的MNC电催化剂在工作条件下表现出较低的稳定性，这是由于其在工作条件中的低稳定性或结构演化。
根据ΔG*OH 的值，作者将这些电催化剂分为三组，即ΔG*OH >1.53eV（PdN4 C8 、PdN4 C10 、PdN4 C12 ）、ΔG*OH <1.23eV（RuN4 C8 ，RuN4 C10 ，RuN4 C12 ，CoN4 C8 ）和1.23eV<ΔG*OH <1.53eV（RhN4 C8 、RhN4 C10 和RhN4 C12 ，IrN4 C8 和IrN4 C10 ），并且ΔG*OH 决定了其在OER过程中是否会出现结构演变。结果证明，由于反应中间体与MNC之间的中等结合能，1.23eV<ΔG*OH <1.53eV的MNC（M=Rh，Ir）具有较高的OER活性。 此外，这些催化剂可以在工作条件（高温、动态条 件、局部电场和强吸附）下保持M-N-C结构，而无需进一步氧化和结构演变，因此具有优异的稳定性。 然而，ΔG *OH >1.53eV或ΔG *OH <1.23eV的MNC电催化剂在工作条件下表现出较低的稳定性，这是由于其在工作条件中的低稳定性或结构演变。 总之， 作者提出了一种MNC电催化剂的综合评价方法，即将ΔG *OH 作为OER活性和稳定性的筛选标准，并将工作条件下的ΔE b 作为稳定性的描述符。 这对在工作条件下设计和筛选ORR、OER和HER电催化剂具有重要意义。
Xiaoming Zhang etr.al Structure evolution and durability of Metal-Nitrogen-Carbon (M = Co, Ru, Rh, Pd, Ir) based oxygen evolution reaction electrocatalyst: A theoretical study Journal of Colloid and Interface Science 2023 https://doi.org/10.1016/j.jcis.2023.02.103
原创文章，作者：计算搬砖工程师，如若转载，请注明来源华算科技，注明出处：https://www.v-suan.com/index.php/2023/11/17/77176290bc/
