本發(fā)明涉及電數(shù)字?jǐn)?shù)據(jù)處理的,特別涉及基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法。
背景技術(shù):
1、單原子催化劑是一種將金屬原子高度分散在載體上的多相催化劑,該材料因其均一的催化活性位點(diǎn)、原子的最大利用率等優(yōu)點(diǎn),日益引起人們的重視。此外,由于其均一的催化活性,使得研究其性質(zhì)機(jī)理更為方便。
2、石墨二炔是一種二維層狀材料,是由sp雜化和sp2雜化的碳原子構(gòu)成的晶體,具有高度共軛的多孔結(jié)構(gòu),良好的導(dǎo)電性和可控電子結(jié)構(gòu)。由于其均一的多孔結(jié)構(gòu),使得其成為單原子催化劑最佳載體之一。由其構(gòu)成的單原子催化劑在各個(gè)領(lǐng)域都有涉及,例如co的氧化、no2的還原、氫化作用、氮的還原反應(yīng)、氧還原反應(yīng)等等。近年來,由石墨二炔負(fù)載單金屬原子催化劑在no2還原方面取得重大進(jìn)展。
技術(shù)實(shí)現(xiàn)思路
1、本發(fā)明的目的是,提供基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,對催化劑的結(jié)構(gòu)優(yōu)化,測試其穩(wěn)定性,分析其電子結(jié)構(gòu)以及吸附no2的性能進(jìn)行評(píng)估,為no2轉(zhuǎn)化催化劑設(shè)計(jì)提供理論基礎(chǔ)。
2、本發(fā)明設(shè)計(jì)了基于第一原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,基于第一性原理計(jì)算得出石墨二炔負(fù)載四種過渡金屬(ni,mo)的幾何構(gòu)型,計(jì)算石墨二炔負(fù)載四種金屬與其相近的碳原子形成的化學(xué)鍵強(qiáng)度來檢驗(yàn)催化劑結(jié)構(gòu)穩(wěn)定性,構(gòu)建吸附no2的模型,對模型進(jìn)行幾何優(yōu)化,確定吸附no2穩(wěn)定結(jié)構(gòu)。同時(shí)產(chǎn)生吸附no2時(shí)的波函數(shù)文件,使用multiwfn軟件對波函數(shù)進(jìn)行分析,計(jì)算出金屬元素與no2結(jié)合形成化學(xué)鍵強(qiáng)度,用來評(píng)估no2與催化劑結(jié)合方式,最后計(jì)算吸附no2的吸附能,結(jié)合化學(xué)鍵強(qiáng)度評(píng)估no2的吸附情況。本發(fā)明為適用于吸附催化no2材料設(shè)計(jì),活化no2,降低吸附催化no2所需能量,以用來節(jié)約成本,同時(shí)可以指導(dǎo)實(shí)驗(yàn)方向。
3、所述方法包含以下步驟:
4、(1)獲取石墨二炔原胞結(jié)構(gòu)文件,文件格式為.cif,原胞文件在materialsproject晶體數(shù)據(jù)庫獲得。將.cif文件轉(zhuǎn)化為.gjf文件,將石墨二炔擴(kuò)胞為2*2*1的超胞,命名為supercell.gjf;
5、(2)利用supercell.gjf生成test.inp文件,對石墨二炔進(jìn)行截?cái)嗄芎蚹點(diǎn)收斂性測試;
6、(3)利用supercell.gjf生成opt.inp文件,使用cp2k軟件對石墨二炔2*2*1的超胞進(jìn)行結(jié)構(gòu)優(yōu)化,通過得到的.restart文件生成優(yōu)化后的.gjf文件,此文件為優(yōu)化后的穩(wěn)定結(jié)構(gòu)構(gòu)型;
7、(4)分別構(gòu)建ni,fe,mo,ru負(fù)載到石墨二炔2*2*1的超胞上的結(jié)構(gòu)模型,生成.inp文件,對二種結(jié)構(gòu)進(jìn)行優(yōu)化,并對優(yōu)化好的結(jié)構(gòu)進(jìn)行能量計(jì)算,在輸出文件中得到負(fù)載二種不同金屬的石墨二炔2*2*1的超胞的能量并得到波函數(shù).molden文件;
8、(5)對步驟4所得波函數(shù).molden文件進(jìn)行分析,計(jì)算mayer鍵級(jí),并用繪圖軟件將其表現(xiàn)出來;
9、(6)構(gòu)建負(fù)載二種不同金屬的石墨二炔吸附no2的結(jié)構(gòu)模型,生成.inp文件,對結(jié)構(gòu)進(jìn)行優(yōu)化,通過輸出文件得到負(fù)載二種不同金屬的石墨二炔吸附no2的穩(wěn)定構(gòu)型,并對優(yōu)化好的結(jié)構(gòu)進(jìn)行能量計(jì)算,在輸出文件中得到二種結(jié)構(gòu)的最終能量并得到波函數(shù).molden文件;計(jì)算一個(gè)no2分子的能量,然后計(jì)算no2的吸附能;對步驟5所得波函數(shù).molden文件進(jìn)行分析,計(jì)算mayer鍵級(jí),并用繪圖軟件將其表現(xiàn)出來;
10、(7)步驟2中,截?cái)嗄軠y試結(jié)果為400ry,k點(diǎn)測試后選用221的k點(diǎn);
11、(8)步驟4中,構(gòu)建模型時(shí)在模型上方留下15埃的真空區(qū),負(fù)載金屬為單原子,構(gòu)成單原子催化劑;
12、(9)步驟5中,使用multiwfn軟件對波函數(shù)進(jìn)行分析,使用origin進(jìn)行繪圖;
13、(10)步驟6中,構(gòu)建模型時(shí)在模型上方留下15埃的真空區(qū);
14、(11)步驟7中,吸附能的計(jì)算公式為eads=etot-eadsorbate-ecatalysts,其中eads為吸附能,etot為吸附no2后的總能量,eadsorbate為no2的能量,ecatalysts為石墨二炔負(fù)載單金屬原子能量;
15、(12)步驟8中,使用multiwfn軟件對波函數(shù)進(jìn)行分析,使用vmd進(jìn)行繪圖。
16、本發(fā)明的有益效果在于,基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,從理論計(jì)算的角度設(shè)計(jì)單原子催化劑,活化no2,降低no2轉(zhuǎn)化為可利用的能源產(chǎn)品所需能量,為實(shí)驗(yàn)提供理論基礎(chǔ)。
1.一種基于第一原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于,具體按照以下步驟實(shí)施:
2.根據(jù)權(quán)利要求2所述的基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于:步驟(2)中,截?cái)嗄軠y試結(jié)果為400ry,k點(diǎn)測試后選用221的k點(diǎn)。
3.根據(jù)權(quán)利要求2所述的基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于:步驟(4)中,構(gòu)建模型時(shí)在模型上方留下15埃的真空區(qū),負(fù)載金屬為單原子,構(gòu)成單原子催化劑。
4.根據(jù)權(quán)利要求2所述的基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于:步驟(5)中,使用multiwfn軟件對波函數(shù)進(jìn)行分析,使用origin進(jìn)行繪圖。
5.根據(jù)權(quán)利要求2所述的基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于:步驟(6)中,構(gòu)建模型時(shí)在模型上方留下15埃的真空區(qū)。
6.根據(jù)權(quán)利要求2所述的基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于:步驟(7)中,吸附能的計(jì)算公式為eads=etot-eadsorbate-ecatalysts,其中eads為吸附能,etot為吸附no2后的總能量,eadsorbate為no2的能量,ecatalysts為石墨二炔負(fù)載單金屬原子能量。
7.根據(jù)權(quán)利要求2所述的基于第一性原理計(jì)算石墨二炔負(fù)載過渡金屬吸附no2的方法,其特征在于:步驟(8)中,使用multiwfn軟件對波函數(shù)進(jìn)行分析。