本发明涉及催化剂分子的开发与利用领域，具体涉及一种催化剂表面活性相与其脱硫催化活性关系的解析方法，其特征是利用密度泛函理论(dft)识别活性相并建立其与脱硫活性构效关系的一种方法。
背景技术：
国家出台政策严格限制石油中硫元素的含量，发展高效脱硫催化剂成为一项紧急的任务。目前钼基硫化物由于其高效的脱硫加氢选择性广泛应用于石油工业脱硫中，但是催化剂在加氢脱硫条件下展现出不同的形貌，过渡金属原子的修饰同样影响原油脱硫速率以及其产物选择性。目前实验手段下很难观察认识多种不同处理手段下的钼基硫化物与其展现出的优异脱硫活性的构效关系，借助dft理论计算模拟，可以从电子层次和原子层次上认识不同处理手段下的催化剂的活性相及其对应的脱硫关系，从而建立活性相与脱硫活性之间的构效关系，为工业生产高效选择性的脱硫催化剂提供理论指导。
技术实现要素：
本发明旨在建立催化剂表面活性相与其脱硫催化活性构效关系的一种解析方法。利用dft研究催化剂的不同活性相特性及模型分子在表面其催化脱硫的反应网络。通过研究不同活性相的催化剂的表面性质，模型分子在其表面的吸附以及脱硫反应网络，结合表面电子结构特性分析、反应势垒与反应速率及热力学与动力学分析，认识催化剂活性相与脱硫活性构效关系。
本发明采用的计算模拟软件为美国accelrys公司研发的materialsstudio(ms)软件，其中的dmol3可用来模拟表面、气相及固体等过程而应用于本发明中。
本发明的内容主要包括催化剂的不同活性相中模型分子的吸附及通过脱硫反应分析认识活性相与催化活性关系两部分，具体内容如下：
(1)催化剂的不同活性相中的模型分子的吸附
系统研究催化剂的不同活性相的结构与性质，通过ms构建需要研究的不同活性相的催化剂，通过前线轨道分析与静电势面的研究，确定其各自的活性相及其本质特性；研究含硫化合物在其各自的活性位上的吸附，分析其之间的电荷转移情况，并分析各种吸附构型的成键与键的活化情况，根据键的活化情况初步预测其活性相所对应的脱硫活化机理；根据态密度、前线轨道、静电势以及密里根电荷布居分析的结果，相互验证吸附能以及成键的本质因素，明确吸附过程的电荷转移情况，初步认识活性相与成键活化的关系。
(2)认识活性相与催化活性关系
选择不同活性相上分别对应的最稳定吸附结构，研究各自结构的脱硫基元反应，认识各基元反应中键的活化与电荷转移特点，明确各反应基元步骤的最佳反应位点，构建各自的脱硫反应网络；结合吸附结构与各自的不同的反应网络，分析不同基元步的发生以及同一基元步不同反应能量与活性相的关系；最终对比分析活性相对脱硫反应的选择性及其各自的最优反应，认识催化剂活性相与脱硫活性之间的本质联系。
本发明的优点在于：提供了一种系统的方法去认识加氢脱硫催化剂活性相与脱硫活性的关系，本计算方法简介快速，结果准确，通过理论计算构建不同的催化剂活性相，认识各自的活性相的特性，明确各自活性相对应的脱硫反应选择性，从而构建催化剂活性相与脱硫反应活性的构效关系，为工业生产脱硫催化剂提供理论指导。
附图说明
本说明书包括如下附图：
图1技术流程图
具体实施方式
本发明所涉及到的相关技术步骤如附图1所示，结合附图1与以下的具体一个实例，说明本发明的具体应用步骤。
实例基于密度泛函理论计算，系统研究了过渡金属co原子两种不同位置修饰下的mos2三角形团簇的表面性质及相关脱硫反应势能面，通过分析比较掌握边缘与角位两种不同活性相所对应的不同脱硫活性，进一步说明本发明的技术方案和优点。具体实施方式包括以下4步：
1.模型建立
1)根据相关实验与理论模拟文献，建立边缘修饰与角位修饰下的两种mos2纳米团簇模型，以噻吩作为汽油模型分子，并对模型分子与催化剂进行试算，与相关文献进行对比验证分析；
2)查阅相关计算文献，选取合适的计算方法，如本案例采取的dft方法与gga-pw91基组作为交换相关泛函计算。
2.计算模拟
1)利用ms软件materialsvisualizer模块，建立h2和噻吩分子模型。根据相关文献建立的三角形mos2团簇的边缘均为mo边缘，且均为50％-s覆盖，利用dmol3模块对催化剂及其它模型分子进行结构优化得到最稳定结构，将边缘一个mo原子进行co原子替代，根据替换的mo原子的位置的不同，分为边缘修饰和角位修饰两种模型，并产生相应的不同空缺位，对其进行优化计算；
2)对催化剂模型进行静电势与轨道计算，另分别在边缘与角位修饰的mos2上进行噻吩分子的吸附以及h原子共吸附时的结构进行优化计算，并进行相关电子结构分析，对反应中间体进行吸附计算，寻找最稳定中间体进行反应计算；
3)以不同修饰位分别的最稳定噻吩吸附态作为初态，各断键加氢中间体为末态，采用ms中的lst/qst方法进行过渡态的搜寻，得到各自的反应势能面，并对相关的级元反应的反应速率进行计算以及相关的热学性质；
4)绘制各自的hds势能面网络，明确各自的脱硫活性以及反应速率。
3.结果分析
1)认识催化剂活性相：明确催化剂表面的静电势面，比较分析不同催化剂对应的lumo轨道，认识路易酸位的存在，比较其吸附空间位阻；
2)吸附与反应基元步分析：比较模型分子在催化剂上的吸附类型、吸附能、模型分子的活化情况，比较相同基元步的不同活化能垒，不同催化剂的不同基元步，并对比反应速率以及各自的速率控制步骤，明确不同修饰位对应的不同反应机理以及脱硫活性与选择性；
3)电子结构分析：结合静电势、前线分析轨道、电荷布局、差分电荷、态密度、能带结构分析等分析活性位及吸附过程电子结构特性；计算不同修饰位催化剂相应的熵变、焓变及自由能等热力学性质。
4.成果提炼
1)揭示co原子修饰位对催化剂活性位的影响；
2)阐释不同位置co修饰催化剂的脱硫活性及各自反应机理；
3)深入认识co修饰位与催化脱硫活性的关系。
以上所述，仅为本发明的实施例而已。并非对本发明作任何形式上的限制；凡熟悉本专业的普通技术人员均可按说明书附图和以上所述而顺畅地实施本发明；但是，凡熟悉本专业的技术人员在不脱离本发明技术方案范围内，可利用以上所揭示的技术内容而作出的些许更改、修饰与演变的等同变化，均为本发明的等效实施例；同时，凡依据本发明的实施技术对以上实施例所作的任何等同变化的更动、修饰与演变等，均属于本发明的技术方案的保护范围之内。