一种金属-氮-碳固氮电催化剂的设计方法与流程

一种金属
‑
氮
‑
碳固氮电催化剂的设计方法
技术领域
1.本发明涉及化学和材料领域，尤其涉及一种基于第一性原理的具有催化氮还原反应活性的石墨烯基电催化剂的设计方法。


背景技术：

2.氨是一种重要的化工原料及能量载体，广泛应用于农业、化工、能源、医药等领域。目前工业合成氨仍采用传统的的haber
‑
bosch方法，该方法需要在铁基催化剂的作用下，在高温高压条件下进行，能耗大，还会产生大量的二氧化碳，排放高。利用电化学方法在环境条件下催化氮还原反应合成氨，具有低耗能、清洁、绿色等优点，成为替代haber
‑
bosch法的工业合成氨的潜在方法之一。然而，现有电催化剂在催化氮还原反应时，受自身性能的限制以及竞争性析氢反应的影响，存在产率低、法拉第效率差等问题，尚无法满足工业合成氨的需要。因此，开发高活性、高稳定性、高选择性的氮还原反应电催化剂，是近年来的研究热点。
3.过渡金属具有丰富的d电子结构，价格低廉，储量丰富，展现出较高的电催化固氮潜能。但部分过渡金属导电率低。氮掺杂碳具有较高的导电性、较强的稳定性、较大的比表面积，在电催化领域有着重要应用。通过氮原子的锚定作用将过渡金属原子均匀的分散在氮掺杂碳表面，可形成单原子、双原子或单团簇等多种类型的催化剂。这些催化剂具有暴露的活性位点，可调的电子结构，是理想的电催化固氮催化剂之一。然而，这些催化剂在设计时缺乏理论指导，针对性不强，在制备方面存在一定困难，需要大量的试错实验，导致成本高，效率低。因此，亟需开发具有高效电化学固氮活性的过渡金属嵌入的氮掺杂碳材料。


技术实现要素：

4.受到实验条件和成本的约束，目前过渡金属氮掺杂碳材料的催化性能研究多以先设计制备、再表征性能的方式，使得催化剂的设计缺乏指导、制备缺乏定向、性能无法预判。
5.本发明基于第一性原理模拟方法，开发出一种筛选具有氮还原性能金属
‑
氮掺杂碳材料的方法，为后期实验制备提供必要的理论指导，加快新型催化剂的研究进度，进而揭示催化反应的机理。
6.本发明中一种金属
‑
氮
‑
碳固氮电催化剂的设计方法是通过下述步骤实现的：
7.步骤1、通过金属负载类型和不同的氮配位构建金属
‑
氮掺杂碳材料的模型(如图1所示)；
8.步骤2、对步骤1构建的模型进行优化；
9.步骤3、计算金属
‑
氮掺杂碳材料的结合能；
10.步骤4、根据步骤3计算的结合能判断金属
‑
氮掺杂碳材料的稳定性，如果稳定执行步骤5，如果不稳定执行步骤13；
11.步骤5、计算氮气吸附金属
‑
氮掺杂碳材料表面的不同活性位点时的吉布斯自由能；
12.步骤6、根据步骤5计算的吉布斯自由能判断氮气吸附能力，当吉布斯自由能小于0时执行步骤6，当吉布斯自由能大于0，执行步骤13；
13.步骤7、计算氢原子在金属
‑
氮掺杂碳材料表面的吸附吉布斯自由能；
14.步骤8、如果步骤5的吉布斯自由能小于步骤7的吉布斯自由能，执行步骤9，否则执行步骤13；
15.步骤9、计算金属
‑
氮掺杂碳材料的第一个氢化步骤(*n2 h

e
‑
→
*n2h)和产生最后一个nh3分子步骤(*nh2 h

e
‑
→
*nh3)的吉布斯自由能变化值；
16.步骤10、如果第一个氢化步骤的吉布斯自由能变化值和产生最后一个nh3分子步骤的吉布斯自由能变化值同时满足小于1.0，执行步骤11，否则执行步骤13；
17.步骤11、根据固氮反应机制构建及优化各个中间体吸附态，并计算固氮反应机制的每一步的吉布斯自由能变化值，在所有吉布斯自由能变化值为正值的步骤中，吉布斯自由能变化值最大的步骤即为整个固氮反应的决速步骤；
18.步骤12、如果金属
‑
氮
‑
碳材料决速步骤的吉布斯自由能变化值小于0.5ev，则为性能优异的氮还原催化剂，如果决速步骤的吉布斯自由能变化值大于等于0.5ev执行步骤13；
19.步骤13、放弃。
20.进一步地限定，步骤1中通过金属负载类型和不同的氮配位利用materials studio软件的结构库构建金属
‑
氮掺杂碳材料模型，并得到相应的cif文件，然后结合vesta软件将cif文件转化为vasp软件计算所需的poscar文件。
21.进一步地限定，步骤2中利用步骤1的poscar文件并结合设置计算所需的incar和kpoints文件，优化金属
‑
氮掺杂碳材料的几何结构，得到相应的contcar文件，利用vesta软件得到结构的关键键长、键角和能量的信息。可通过键长键角信息有助于了解材料的结构。
22.进一步地限定，步骤3中根据步骤2获得的能量信息，利用公式(1)计算金属
‑
氮掺杂碳材料的结合能e
b
：所有能量的单位为ev：
23.e
b
＝e(mn4‑
mn4/c)
–
e(nc)
–
e(m
总
) (1)
24.其中，e(mn4‑
mn4/c)为金属
‑
氮掺杂碳材料整体的能量，e(nc)为不包含金属原子的氮掺杂碳材料的能量，e(m
总
)为所有金属单原子的能量总和，m
总
代表所有金属单原子。
25.进一步地限定，步骤5中根据公式(2)计算氮气吸附金属
‑
氮掺杂碳材料表面的不同活性位点时的吉布斯自由能
[0026][0027]
其中，*代表金属
‑
氮掺杂碳材料，*n2代表n2吸附在金属
‑
氮掺杂碳材料表面，每种结构的吉布斯自由能为g＝e e
zpe
–
ts，e为结构的电子能量，e
zpe
为零点能，s为熵，温度t的取值为298.15k。
[0028]
氮气吸附在金属
‑
氮
‑
碳材料主要有两种构型，即竖直和水平方式，参见图3；氮气竖直吸附在表面时，一个n原子与表面连接；氮气水平吸附在表面时，两个n原子与表面连接。
[0029]
进一步地限定，步骤7中根据公式(3)计算氢原子在金属
‑
氮掺杂碳材料表面的吸附吉布斯自由能δg
*h
：
[0030]
[0031]
其中，*代表金属
‑
氮掺杂碳材料，*h代表h原子吸附在金属
‑
氮掺杂碳材料表面，质子和电子的能量和等于1/2氢气的能量，每种结构的吉布斯自由能为g＝e e
zpe
–
ts，e为结构的电子能量，e
zpe
为零点能，s为熵，温度t的取值为298.15k。
[0032]
进一步地限定，步骤9中利用公式(4)计算第一个氢化步骤(*n2 h

e
‑
→
*n2h)的吉布斯自由能变化值δg(*n2 h

e
‑
→
*n2h)：
[0033][0034]
利用公式(5)计算产生最后一个nh3分子步骤(*nh2 h

e
‑
→
*nh3)的吉布斯自由能变化值δg(*nh2 h

e
‑
→
*nh3)：
[0035][0036]
其中，*代表金属
‑
氮掺杂碳材料，*n2h和*nh3分别代表n2h和nh3吸附在金属
‑
氮掺杂碳材料表面，质子和电子的能量和等于1/2氢气的能量，每种结构的吉布斯自由能为g＝e e
zpe
–
ts，e为结构的电子能量，e
zpe
为零点能，s为熵，温度t的取值为298.15k。
[0037]
进一步地限定，步骤11中固氮反应机制为远端、交替、连续和酶促机制。
[0038]
高稳定性是一切性能的基础。过渡金属和氮掺杂碳材料种类繁多，结构多样，且金属与氮原子的配位模式直接影响材料的稳定性。设计不同组成、结构、配位模式的金属
‑
氮掺杂碳材料，结合第一性原理计算，以金属
‑
氮掺杂碳材料的结合能作为判断依据，筛选出稳定性高的金属
‑
氮掺杂碳材料。
[0039]
利用第一性原理计算方法研究稳定性高的金属
‑
氮掺杂碳材料的电催化固氮性能，筛选出竞争性析氢能力弱、固氮性能强的催化剂。
[0040]
本发明提供了一种基于第一性原理筛选固氮金属
‑
氮
‑
碳催化剂的方法和流程，为实验研究提供直接的理论指导，避免大量试错实验造成的时间和成本损失。
附图说明
[0041]
图1是本发明固氮金属
‑
氮
‑
碳催化剂的流程图；
[0042]
图2是金属
‑
氮掺杂碳材料模型；
[0043]
图3是氮气在不同活性位点的吸附构型；(a)、(c)——竖直吸附构型，(b)、(d)——水平吸附构型；
[0044]
图4是固氮反应的四种机制(远端、交替、连续和酶促机制)
[0045]
图5是crn4‑
mnn4/c材料远端机制和交替机制的结构图。
具体实施方式
[0046]
实施例1:本实施中一种金属
‑
氮
‑
碳固氮电催化剂的设计方法是通过下述步骤实现的：
[0047]
步骤1、通过金属负载类型和不同的氮配位利用materials studio软件的结构库构建金属
‑
氮掺杂碳材料模型，并得到相应的cif文件，然后结合vesta软件将cif文件转化为vasp软件计算所需的poscar文件；
[0048]
步骤2、对步骤1构建的模型进行优化：利用步骤1的poscar文件并结合设置计算所需的incar和kpoints文件，优化金属
‑
氮掺杂碳材料的几何结构，得到相应的contcar文件，
利用vesta软件得到结构的关键键长、键角和能量的信息；参见表1。
[0049]
表1mn4‑
m’n4/c材料的关键键长和结合能(ev)
[0050][0051]
步骤3、计算金属
‑
氮掺杂碳材料的结合能：根据步骤2获得的能量信息，利用公式(1)计算金属
‑
氮掺杂碳材料的结合能e
b
：
[0052]
e
b
＝e(mn4‑
mn4/c)
–
e(nc)
–
e(m)
–
e(m’) (1)
[0053]
其中，e(mn4‑
mn4/c)为金属
‑
氮掺杂碳材料整体的能量，e(nc)为不包含金属原子的氮掺杂碳材料的能量，e(m)为金属单原子m的能量，e(m’)为金属单原子m’的能量，m和m’表示为相同金属单原子，也可以表示不同的金属原子；
[0054]
步骤4、根据步骤3计算的结合能判断的稳定性，如果稳定执行步骤5，如果不稳定执行步骤13；
[0055]
步骤5、根据公式(2)计算氮气吸附金属
‑
氮掺杂碳材料表面的不同活性位点时的吉布斯自由能
[0056][0057]
其中，*代表金属
‑
氮掺杂碳材料，*n2代表n2吸附在金属
‑
氮掺杂碳材料表面，每种结构的吉布斯自由能为g＝e e
zpe
–
ts，e为结构的电子能量，e
zpe
为零点能，s为熵，温度t的取值为298.15k。
[0058]
步骤6、根据步骤5计算的吉布斯自由能判断氮气吸附能力，当吉布斯自由能小于0时执行步骤6，当吉布斯自由能大于0，执行步骤13；
[0059]
步骤7、根据公式(3)计算氢原子在金属
‑
氮掺杂碳材料表面的吸附吉布斯自由能δg
*h
：
[0060][0061]
其中，*代表金属
‑
氮掺杂碳材料，*h代表h原子吸附在金属
‑
氮掺杂碳材料表面，质子和电子的能量和等于1/2氢气的能量，每种结构的吉布斯自由能为g＝e e
zpe
–
ts，e为结构的电子能量，e
zpe
为零点能，s为熵，温度t的取值为298.15k；
[0062]
步骤8、如果δg
*n2
<δg
*h
，执行步骤9，否则执行步骤13；
[0063]
步骤9、利用公式(4)计算第一个氢化步骤(*n2 h

e
‑
→
*n2h)的吉布斯自由能变化值δg(*n2 h

e
‑
→
*n2h)：
[0064][0065]
利用公式(5)计算产生最后一个nh3分子步骤(*nh2 h

e
‑
→
*nh3)的吉布斯自由能变化值δg(*nh2 h

e
‑
→
*nh3)：
[0066][0067]
其中，*代表金属
‑
氮掺杂碳材料，*n2h和*nh3分别代表n2h和nh3吸附在金属
‑
氮掺杂碳材料表面，质子和电子的能量和等于1/2氢气的能量，每种结构的吉布斯自由能为g＝e e
zpe
–
ts，e为结构的电子能量，e
zpe
为零点能，s为熵，温度t的取值为298.15k；
[0068]
表2mn4‑
m’n4/c材料的氮气吸附位点和构型、氮气和氢原子吸附自由能(ev)
[0069][0070]
步骤10、如果δg(*n2→
*n2h)<1.0同时满足δg(*nh2→
*nh3)<1.0，执行步骤11，否则执行步骤13；
[0071]
步骤11、根据固氮反应机制构建及优化各个中间体吸附态，并计算固氮反应机制的每一步的吉布斯自由能变化值，在所有吉布斯自由能变化值为正值的步骤中，吉布斯自由能变化值最大的步骤即为整个固氮反应的决速步骤δg
max
；
[0072]
其中，吉布斯自由能变化值计算根据反应机制构每一步进行的具体反应选择公式(2)、公式(3)、公式(4)、公式(5)进行计算。
[0073]
表3催化剂在不同机制中每一步的吉布斯自由能数值(δg/ev)
[0074][0075][0076]
步骤12、如果金属
‑
氮
‑
碳材料决速步骤的吉布斯自由能变化值δg
max
小于0.5ev，则为性能优异的氮还原催化剂，如果决速步骤的吉布斯自由能变化值δg
max
大于等于0.5ev执行步骤13；
[0077]
其中，表3中决速步骤δg
max
小于0.5ev，则为性能优异的氮还原催化剂；
[0078]
步骤13、放弃。