CN104573297A - 一种确定表面催化反应路径的方法 - Google Patents

一种确定表面催化反应路径的方法 Download PDF

Info

Publication number
CN104573297A
CN104573297A CN201310503600.5A CN201310503600A CN104573297A CN 104573297 A CN104573297 A CN 104573297A CN 201310503600 A CN201310503600 A CN 201310503600A CN 104573297 A CN104573297 A CN 104573297A
Authority
CN
China
Prior art keywords
reaction
model
energy
elementary reaction
reactant
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310503600.5A
Other languages
English (en)
Other versions
CN104573297B (zh
Inventor
李明
王晓霖
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Fushun Research Institute of Petroleum and Petrochemicals
Original Assignee
China Petroleum and Chemical Corp
Sinopec Fushun Research Institute of Petroleum and Petrochemicals
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Fushun Research Institute of Petroleum and Petrochemicals filed Critical China Petroleum and Chemical Corp
Priority to CN201310503600.5A priority Critical patent/CN104573297B/zh
Publication of CN104573297A publication Critical patent/CN104573297A/zh
Application granted granted Critical
Publication of CN104573297B publication Critical patent/CN104573297B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

本发明提供一种确定表面催化反应路径的方法，包括如下步骤：首先构建催化材料表面结构模型和反应物分子结构模型，然后构建反应物在催化材料表面吸附结构模型，选择最稳定吸附结构模型作为下一步反应的初态结构；构建反应物涉及到所有键的断裂和生成的基元反应末态吸附结构模型，选择最稳定吸附结构模型作为该基元反应的末态；以确定的初态结构和末态结构构建所有基元反应轨迹文件，利用MS软件Castep模块中的TSSearch功能对该轨迹文件进行过渡态计算，以确定的中间体依次重复上述步骤，逐步构建整个基元反应网络，直到得到最终产物，从而得到完整的表面催化反应路径。本发明方法计算迅速，操作简单，克服了实验仪器分析确定反应路径中人为推断的弊端。

Description

一种确定表面催化反应路径的方法
技术领域
本发明涉及一种确定表面催化反应路径的方法，特别是涉及一种利用计算机模拟技术对表面催化反应进行微观动力学模拟和计算确定表面催化反应路径的方法。
背景技术
催化作用是现代世界最重要的技术之一，特别是多相催化反应过程一直以来都是化学化工过程的基础。不同的催化剂对特定的反应体系具有选择性，因此有必要确定不同催化反应体系的反应机理和反应路径，以便于更好的设计催化反应过程及催化剂。世界各地的科学家们多年来想尽各种办法，在多种分析仪器上试图原位研究多相催化反应，但由于多相催化反应过程的复杂性，人们也只是在某些方面了解到了一些多相催化反应过程的细节。多相催化反应通常基于程序升温脱附、程序升温还原等实验方法以及高分辨率电子能损失谱、质谱、俄歇、核磁共振谱等表面分析技术得到催化剂结构信息，反应中间物种的检测进而推断反应机理和反应路径，已有不少实例证明可以解决一些这方面的问题。但是实验仪器分析方法需要经过繁琐的实验过程，同时仪器分析方法是通过测定表面信息来人为推测表面物种进而得到反应路径，但该物种是表面催化反应的中间体还是分支反应产物无法确定，这样就导致推测得到的反应路径具有很大的不准确性。另外，如果某一中间体在表面的反应速度很快，那么其在表面的停留时间极短，检测仪器根本无法实施检测。因此，由于表面催化反应自身的复杂性和仪器设备的限制，单通过实验仪器分析方法要确定表面催化反应路径，工作量和难度都很大，甚至是不可能进行。此外，实验方法繁琐、费时、费力，造成人力财力的巨大消耗。
量子力学是20世纪最重要的科学发现之一。在量子力学基础上发展起来的量子化学计算，为我们开辟了通向微观世界的又一个途径。量子化学是物理化学中应用量子力学基本原理和方法在原子层面讨论化学过程的分支学科。基于量子化学计算的模拟软件已成功地应用于化学科学的研究领域。目前量子化学计算对表面催化反应通常都基于实验分析提出的路径进行微观机理研究。为此，本发明提出用量子化学方法通过基元反应动力学计算确定表面催化反应路径的方法。
发明内容
针对现有技术中研究表面催化反应路径基于实验检测分析推断的技术缺陷，本发明提供一种利用计算机模拟确定表面催化反应路径的方法。
本发明提供一种确定表面催化反应路径的方法，包括如下步骤：
（1）构建催化材料表面结构模型，基于晶体结构数据库，选择催化材料晶体的晶胞结构，或者根据已有的催化材料晶体结构数据利用MS软件的Visualizer模块构建催化材料的晶胞结构并对其能量最小化处理，然后以选择或处理后的催化材料的晶胞结构构建周期晶面结构或团簇晶体结构模型，用所构建的周期晶面结构或团簇晶体结构模型模拟催化材料表面结构，并对构建的催化材料表面结构模型进行能量最小化处理，记录相应的能量EM；
（2）构建反应物分子结构模型，根据反应物的化学分子式，采用MS软件的Visualizer模块构建反应物分子结构模型，然后利用Castep模块Calculation功能下的Geometry Optimization程序对其进行能量最小化处理，相应的能量记录为ER；
（3）构建反应物在催化材料模型表面吸附结构模型，利用MS软件Visualizer模块构建反应物在催化剂模型表面所有的吸附结构模型，采用Castep模块的Geometry Optimization程序对所有吸附结构模型进行能量最小化处理，选择出能量最小的吸附结构模型，对应的能量为EI，通过计算该吸附过程的吸附能Ea，Ea= EM + ER – EI，选择吸附能最大，即最稳定吸附结构模型，并将该吸附结构模型作为下一步反应过程计算的初态结构；
（4）构建反应物所有基元反应末态吸附结构模型，利用MS软件Visualizer模块构建反应物涉及到所有键的断裂和生成的基元反应末态吸附结构模型，并采用Castep模块的Geometry Optimization程序对每一个基元反应末态的所有吸附结构模型进行能量最小化处理，选择能量最低，即最稳定吸附结构模型，作为该基元反应的末态，相应能量记为EP；
（5）构建反应物所有基元反应过程，确定其相应的基元反应过渡态，以步骤（3）确定的初态结构和步骤（4）确定的末态结构，利用软件的Reaction Preview程序进行原子match，当所有原子匹配完成后进行反应preview，得到基元反应轨迹文件，利用MS软件Castep模块中Calculation功能下的TS Search程序对该轨迹文件进行过渡态计算，找到基元反应过渡态结构并得到相应的能量ET，由每个基元反应过渡态结构对应的能量计算基元反应能垒∆E+，选取能垒最小的基元反应产物作为下一步反应的反应物，如果有两个或两个以上基元反应能垒接近，则通过基元反应末态能量与过渡态能量计算逆向反应能垒∆E-进行比较，选择逆向反应能垒高的基元反应末态作为下一步反应的反应物；
（6）构建整个基元反应网络，确定催化反应路径，以步骤（5）确定的中间体依次重复步骤（3）、（4）和（5），逐步构建整个基元反应网络，直到得到反应最终产物，从而得到完整的表面催化反应路径。
本发明方法中，步骤（1）中所述晶体结构数据库可以为MS软件自带的晶体结构数据库或ICSD等商业晶体结构数据库。
本发明方法中，步骤（1）中所述已有催化材料晶体结构数据是通过XRD等实验手段获得的。
在所述步骤（1）中，构建周期晶面结构表面模型时，表面周期模型层数可以为3～8层，优选4～5层，每层分子数4～25个，优选9～16个。构建团簇晶体结构模型时，催化材料结构可以为任意晶体结构，分子数可以为1～300，优选50～150。
在所述步骤（2）、（3）、（4）、（5）中对所述结构文件进行能量最小化处理时，计算函数可选用LDA中的CA-PZ和GGA中的PBE、PW91及RPBE方法中的一种，优选为GGA的RPBE方法。计算最大迭代次数设定为1～500步，优选为150～350步。计算精度可以设置为Media、Fine、Ultra-Fine中的一种，优选Fine。
在所述步骤（3）中构建催化材料表面吸附模型时，在催化材料表面模型上吸附1个反应物；吸附计算时对反应物在催化剂模型表面的所有高对称吸附位都给予考虑计算。
在所述步骤（5）中寻找基元反应过渡态结构时，计算函数可以是LST Maximum、Halgren-Lipscomb、LST/Optimization、Complete LST/QST、QST/Optimization中的一种，优选Complete LST/QST；计算精度可以设置为Media、Fine、Ultra-Fine中的一种，优选Fine；计算最大QST迭代次数设定为1～30步，优选为10～20步。
本发明计算方法中，其他需要设置的计算参数均采用以上设置参数设定后计算软件系统的默认值。
本发明方法中，所述能量最小化处理采用量子化学软件进行，所述量子化学软件可以为VASP、SIESTA、Material Studio的Castep和Dmol3模块中的任一种，优选为MS软件的Castep模块。
本发明方法中，在确定表面催化反应路径时，在催化剂表面结构模型和初始反应物分子结构模型的基础上，采用量子化学方法进行计算，首先对涉及到的反应物和中间体计算吸附时放出的热量Ea，Ea= EM + ER – EI，选用放出热量最大即能量最稳定的结构作为反应初态，同时抛开化学反应思想，对反应物涉及到的所有化学键断裂和生成的基元反应都予以考虑和计算，得到基元反应动力学能垒∆E+，∆E+ = ET – EI和逆反应动力学能垒（∆E-，∆E- = ET – EP），采用计算对比平行基元反应动力学能垒的方法，即选择基元反应动力学能垒最低且逆反应动力学能垒较高者进行进一步计算，从而全面考察涉及到的基元反应过程并构建反应网络，确定出反应最低势能面并得到反应路径。
本发明方法中，所述的Material Studio软件为中文正式版计算机应用软件。
采用不同量子化学计算软件，本发明方法计算时所调用的模块和设定的参数会有少许不同，但不影响本发明方法确定表面催化反应路径的思想。
本发明确定表面催化反应路径的方法与通过实验研究得到的表面催化反应路径一致，与传统实验方法和理论研究方法相比，计算迅速，操作简单，免去了繁重的实验工作，节省了大量的实验成本，同时克服了实验仪器分析确定反应路径中人为推断的弊端，全面考察反应网络，避免了传统量子化学计算研究依赖实验结论的不完整性。本发明方法应用范围广，对固体催化剂表面催化反应均可以应用。设备要求低，在普通计算机上即可以进行计算，不需要购买大量的实验设备。计算结果准确，可以深入解释实验结果，作为实验结论的重要补充。
附图说明
图1 为本发明方法流程示意图。
具体实施方式
下面结合具体实施方式对本发明作进一步的描述，并非对其保护范围的限制：
本实施例采用商业软件Material Studio（MS）实现。该软件由美国Accelrys公司研发，是专门为材料科学领域研究者开发的一款模拟软件，支持Windows、Linux等多种操作系统。
实施例1
一种确定表面催化反应路径的方法，表面催化反应选为贵金属Au催化剂表面甲醇氧化反应。反应路径确定方法如下：
（1）利用MS软件的Visualizer模块，基于MS软件自带的晶体结构数据库，导出Au元素的晶胞，并以此构建Au(111)表面周期模型，该模型选取5层结构，每层包含9个Au原子来模拟贵金属Au催化剂表面；用MS软件Castep模块Calculation功能下的Geometry Optimization程序对其进行结构优化，计算函数选取为GGA+RPBE，最大迭代次数设置为300，计算精度选为fine。
（2）利用MS软件的Visualizer模块构建甲醇分子结构，并利用Castep模块的Geometry Optimization程序对其进行结构优化，计算参数设置同上。
（3）利用MS软件的Visualizer模块构建Au(111)表面所有高对称位甲醇分子的可能吸附结构模型，并全部采用Castep模块的Geometry Optimization程序进行能量最小化处理，计算参数设置同上。计算各吸附结构相应的吸附能，记录于表1中。由表1中数据可知，甲醇分子只能稳定吸附在Au原子的顶位，因此以顶位吸附的甲醇分子作为进一步反应过程计算的初态结构。
（4）利用MS软件的Visualizer模块构建甲醇分子的所有可能基元反应末态吸附结构模型。甲醇分子结构包含有C-H、O-H、C-O三种键型，因此对这三种键断裂反应都进行考虑并构建相应产物表面吸附模型，利用Castep模块的Geometry Optimization程序进行能量最小化处理，将得到的最稳定吸附结构作为各个基元反应末态，计算参数设置同上。
（5）基于前两步确定的甲醇吸附的初态结构及不同键断裂的末态结构，利用软件的Reaction Preview程序进行match，当所有原子匹配完成后preview，从而创建甲醇的三种基元反应过程轨迹文件，利用Castep计算模块中的TS Search程序对轨迹文件进行基元反应过渡态寻找，计算函数选择Complete LST/QST，计算精度设置为Fine，计算最大QST迭代次数设定为10。由每个基元反应过渡态结构对应的能量计算基元反应能垒及逆反应能垒，记录于表2中。
（6）由表2中数据可以看出，甲醇O-H键断裂基元反应能垒最低，因此确定甲醇在Au(111)表面进行O-H键断裂反应生成CH3O中间体。进一步计算CH3O中间体在Au(111)表面的反应。
（7）重复（2）、（3）步，找到CH3O中间体在催化剂表面的所有稳定吸附构型，并记录于表1中。由表1中数据看到，CH3O中间体在表面可以稳定吸附于空位形成两种吸附结构，其中hcp空位吸附放热量较大更为稳定，则以hcp空位吸附结构作为CH3O中间体的基元反应初态。
（8）CH3O中间体包含C-H和C-O两种键型，重复（4）、（5）步，构建两种键断裂产物吸附结构并进行能量最小化作为基元反应末态结构；创建不同基元反应轨迹文件并利用Castep模块中的TS Search程序实现基元反应过渡态寻找，由每个基元反应过渡态结构对应的能量计算基元反应能垒及逆反应能垒，记录于表2中。
（9）由表2数据可知，CH3O中间体的C-H键断裂能垒较低，生成表面甲醛。以此类推，重复上述过程并构建整个反应网络，其中，对甲醛、CHO、CO和COOH中间体计算了氧化剂的氧化反应，直到生成最后的氧化产物CO2和H2O。这样得到甲醇在Au(111)表面氧化路径为CH3OH → CH3O → CH2O → CHO → CO → COOH → CO2 + H2O。
表1 Au(111)表面甲醇分子的所有吸附结构模型的吸附能
吸附物 吸附位 吸附能量 Ea （eV）
CH3OH 顶位 0.56
CH3O 空位-hcp η 1(O) 2.54
空位-fcc η 1(O) 2.52
CH2O 顶位 η 1(O) 0.45
空位-hcp η 1(C)-η 2(O) 0.97
空位-hcp η 2(C)-η 1(O) 0.97
空位-fcc η 1(C)-η 2(O) 0.90
空位-fcc η 3(C)-η 1(O) 0.91
CHO 桥位 η 1(C)-η 1(O) 2.44
空位-hcp η 2(C)-η 1(O) 2.36
空位-fcc η 2(C)-η 1(O) 2.26
CO 顶位 η 1(C) 1.93
桥位 η 2(C) 1.66
空位-hcp η 3(C) 1.71
空位-fcc η 3(C) 1.70
H2COOH 桥位 η 1(O)-η 1(O) 2.62
空位-hcp η 2(O)-η 1(O) 2.77
空位-hcp η 1(O) 2.59
空位-fcc η 1(O)-η 1(O) 2.61
空位-fcc η 1(O) 2.57
HCOOH 顶位 η 1(O) 0.66
COOH 桥位 η1 (C)-η1 (O) 2.69
空位-hcp η1 (C)-η2 (O) 2.53
空位-fcc η1 (C)-η2 (O) 2.53
CO2 桥位 η1 (C)-η1 (O) 0.05
OH* 顶位 η 1(O) 2.84
桥位 η 2(O) 3.14
空位-hcp η 3(O) 3.19
空位-fcc η 3(O) 3.19
表2 基元反应能垒和逆反应能垒
反应物 基元反应 ∆E+（eV） ∆E-（eV）
CH3OH CH3OH* → [CH3O + H]* 0.80 1.20
CH3OH* → [CH2OH + H]* 1.04 1.14
CH3OH* → [CH3 + OH]* 1.88 2.51
CH3O CH3O* → [CH2O + H]* 1.05 1.12
CH3O* → [CH3 + O]* 1.51 2.15
CH2O CH2O* → [CHO + H]* 0.74 1.22
CH2O* → [CH2 + O]* 1.36 2.04
[CH2O + OH]* → H2COOH* 1.02 1.15
CHO CHO* → [CO + H]* 0.44 1.61
CHO* → [CH + O]* 1.70 2.16
[CHO + OH]* → HCOOH* 1.91 1.92
CO CO → C + O 3.05 2.09
[CO + OH]* → COOH* 1.74 1.23
COOH COOH* → [CO2 + H]* 1.22 1.23
[COOH + OH]* → [CO2 + H2O]* 0.06 0.13
本实施方式提出的反应路径中间体与实验检测到的中间体完全一致，并对该实验给予了很好的解释，证明了上述方法的可行性及实用性。由于催化剂晶体结构参数可以从晶体结构数据库里获得，并以此构造各类不同催化剂晶体模型，加上计算机硬件和软件就可以进行计算。得到的结论既可以从理论上解释表面催化反应机理，又可以应用于实际工作确定反应路径。本发明简单合理，容易操作，成本低，具有实用价值，并对其他各类催化反应具有借鉴意义。

Claims (17)

1.一种确定表面催化反应路径的方法，包括如下步骤：
（1）构建催化材料表面结构模型，基于晶体结构数据库，选择催化材料晶体的晶胞结构，或者根据已有的催化材料晶体结构数据利用MS软件的Visualizer模块构建催化材料的晶胞结构并对其能量最小化处理，然后以选择或处理后的催化材料的晶胞结构构建周期晶面结构或团簇晶体结构模型，用所构建的周期晶面结构或团簇晶体结构模型模拟催化材料表面结构，并对构建的催化材料表面结构模型进行能量最小化处理，记录相应的能量EM；
（2）构建反应物分子结构模型，根据反应物的化学分子式，采用MS软件的Visualizer模块构建反应物分子结构模型，然后利用Castep模块Calculation功能下的Geometry Optimization程序对其进行能量最小化处理，相应的能量记录为ER；
（3）构建反应物在催化材料模型表面吸附结构模型，利用MS软件Visualizer模块构建反应物在催化剂模型表面所有的吸附结构模型，采用Castep模块的Geometry Optimization程序对所有吸附结构模型进行能量最小化处理，选出能量最小的吸附结构模型，对应的能量为EI，通过计算该吸附过程的吸附能Ea，Ea= EM + ER – EI，选择吸附能最大，即最稳定吸附结构模型，并将该吸附结构模型作为下一步反应过程计算的初态结构；
（4）构建反应物所有基元反应末态吸附结构模型，利用MS软件Visualizer模块构建反应物涉及到所有键的断裂和生成的基元反应末态吸附结构模型，并采用Castep模块的Geometry Optimization程序对每一个基元反应末态的所有吸附结构模型进行能量最小化处理，选择能量最低，即最稳定吸附结构模型，作为该基元反应的末态，相应能量记为EP；
（5）构建反应物所有基元反应过程，确定其相应的基元反应过渡态，以步骤（3）确定的初态结构和步骤（4）确定的末态结构，利用软件的Reaction Preview程序进行原子match，当所有原子匹配完成后进行反应preview，得到基元反应轨迹文件，利用MS软件Castep模块中Calculation功能下的TS Search程序对该轨迹文件进行过渡态计算，找到基元反应过渡态结构并得到相应的能量ET，由每个基元反应过渡态结构对应的能量计算基元反应能垒∆E+，选取能垒最小的基元反应产物作为下一步反应的反应物，如果有两个或两个以上基元反应能垒接近，则通过基元反应末态能量与过渡态能量计算逆向反应能垒∆E-进行比较，选择逆向反应能垒高的基元反应末态作为下一步反应的反应物；
（6）构建整个基元反应网络，确定催化反应路径，以步骤（5）确定的中间体依次重复步骤（3）、（4）和（5），逐步构建整个基元反应网络，直到得到反应最终产物，从而得到完整的表面催化反应路径。
2.按照权利要求1所述的方法，其特征在于：步骤（1）中所述晶体结构数据库为MS软件自带的晶体结构数据库或ICSD等商业晶体结构数据库。
3.按照权利要求1所述的方法，其特征在于：步骤（1）中，构建周期晶面结构表面模型时，表面周期模型层数为3～8层，每层分子数为4～25个。
4.按照权利要求1 所述的方法，其特征在于：步骤（1）中，构建周期晶面结构表面模型时，表面周期模型层数为4～5层，每层分子数为9～16个。
5.按照权利要求1所述的方法，其特征在于：步骤（1）中，构建团簇晶体结构模型时，催化材料结构为任意晶体结构，分子数为1～300。
6.按照权利要求1所述的方法，其特征在于：步骤（1）中，构建团簇晶体结构模型时，催化材料结构为任意晶体结构，分子数为50～150。
7.按照权利要求1所述的方法，其特征在于：步骤（3）中构建催化材料表面吸附模型时，在催化材料表面模型上吸附1个反应物，吸附计算时对反应物在催化剂模型表面的所有高对称吸附位都进行计算。
8.按照权利要求1所述的方法，其特征在于：步骤（2）、（3）、（4）、（5）中所述能量最小化处理采用量子化学软件进行。
9.按照权利要求8所述的方法，其特征在于：所述量子化学软件为VASP、SIESTA、Material Studio的Castep和Dmol3模块中的任一种。
10.按照权利要求8所述的方法，其特征在于：所述量子化学软件为MS软件的Castep模块。
11.按照权利要求1所述的方法，其特征在于：步骤（2）、（3）、（4）、（5）中所述进行能量最小化处理时，计算函数为LDA中的CA-PZ、GGA中的PBE、PW91和RPBE方法中的一种。
12.按照权利要求1所述的方法，其特征在于：步骤（2）、（3）、（4）、（5）中所述进行能量最小化处理时，计算最大迭代次数设定为1～500步，计算精度设置为Media、Fine、Ultra-Fine中的一种。
13.按照权利要求1所述的方法，其特征在于：步骤（5）中寻找基元反应过渡态结构时，计算函数为LST Maximum、Halgren-Lipscomb、LST/Optimization、Complete LST/QST、QST/Optimization中的一种。
14.按照权利要求1所述的方法，其特征在于：步骤（5）中寻找基元反应过渡态结构时，计算精度设置为Media、Fine、Ultra-Fine中的一种，计算最大QST迭代次数设定为1～30步。
15.按照权利要求1所述的方法，其特征在于：步骤（5）中寻找基元反应过渡态结构时，计算最大QST迭代次数设定为1～30步。
16.按照权利要求1所述的方法，其特征在于：步骤（5）中所述基元反应能垒∆E+，∆E+ = ET – EI。
17.按照权利要求1所述的方法，其特征在于：步骤（5）中所述逆反应动力学能垒∆E-，∆E- = ET – EP。
CN201310503600.5A 2013-10-24 2013-10-24 一种确定表面催化反应路径的方法 Active CN104573297B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310503600.5A CN104573297B (zh) 2013-10-24 2013-10-24 一种确定表面催化反应路径的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310503600.5A CN104573297B (zh) 2013-10-24 2013-10-24 一种确定表面催化反应路径的方法

Publications (2)

Publication Number Publication Date
CN104573297A true CN104573297A (zh) 2015-04-29
CN104573297B CN104573297B (zh) 2017-11-24

Family

ID=53089350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310503600.5A Active CN104573297B (zh) 2013-10-24 2013-10-24 一种确定表面催化反应路径的方法

Country Status (1)

Country Link
CN (1) CN104573297B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446583A (zh) * 2016-10-19 2017-02-22 南京工程学院 一种高能离子盐高压行为的预测方法
CN109002681A (zh) * 2018-05-25 2018-12-14 胡培君 模拟/预测表面催化反应中催化剂催化活性的方法及其应用
CN111048160A (zh) * 2018-10-15 2020-04-21 中国石油化工股份有限公司 预测催化剂性能的方法
CN111063395A (zh) * 2018-10-15 2020-04-24 中国石油化工股份有限公司 筛选多金属加氢催化剂的模拟方法
CN115491709A (zh) * 2022-05-23 2022-12-20 电子科技大学 一种基于二维硼碳基材料高效还原二氧化碳的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101454461A (zh) * 2005-11-16 2009-06-10 Ambrx公司 包括非天然氨基酸的方法和组合物
WO2010044191A1 (ja) * 2008-10-15 2010-04-22 国立大学法人山口大学 化学反応遷移状態探索システムとその方法とそのプログラム
CN103020413A (zh) * 2011-09-28 2013-04-03 中国石油化工股份有限公司 计算机计算芳烃加氢反应活化能及反应速率常数的方法
CN103547366A (zh) * 2011-03-04 2014-01-29 里兰斯坦福初级大学理事会 用于将二氧化碳还原成甲醇的催化剂

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101454461A (zh) * 2005-11-16 2009-06-10 Ambrx公司 包括非天然氨基酸的方法和组合物
WO2010044191A1 (ja) * 2008-10-15 2010-04-22 国立大学法人山口大学 化学反応遷移状態探索システムとその方法とそのプログラム
CN103547366A (zh) * 2011-03-04 2014-01-29 里兰斯坦福初级大学理事会 用于将二氧化碳还原成甲醇的催化剂
CN103020413A (zh) * 2011-09-28 2013-04-03 中国石油化工股份有限公司 计算机计算芳烃加氢反应活化能及反应速率常数的方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106446583A (zh) * 2016-10-19 2017-02-22 南京工程学院 一种高能离子盐高压行为的预测方法
CN106446583B (zh) * 2016-10-19 2019-03-05 南京工程学院 一种高能离子盐高压行为的预测方法
CN109002681A (zh) * 2018-05-25 2018-12-14 胡培君 模拟/预测表面催化反应中催化剂催化活性的方法及其应用
CN109002681B (zh) * 2018-05-25 2021-08-27 胡培君 模拟/预测表面催化反应中催化剂催化活性的方法及其应用
CN111048160A (zh) * 2018-10-15 2020-04-21 中国石油化工股份有限公司 预测催化剂性能的方法
CN111063395A (zh) * 2018-10-15 2020-04-24 中国石油化工股份有限公司 筛选多金属加氢催化剂的模拟方法
CN111063395B (zh) * 2018-10-15 2023-03-14 中国石油化工股份有限公司 筛选多金属加氢催化剂的模拟方法
CN111048160B (zh) * 2018-10-15 2023-05-12 中国石油化工股份有限公司 预测催化剂性能的方法
CN115491709A (zh) * 2022-05-23 2022-12-20 电子科技大学 一种基于二维硼碳基材料高效还原二氧化碳的方法

Also Published As

Publication number Publication date
CN104573297B (zh) 2017-11-24

Similar Documents

Publication Publication Date Title
Nazarian et al. Large-scale refinement of metal− organic framework structures using density functional theory
Evans et al. Computational chemistry methods for nanoporous materials
Raimondeau et al. Recent developments on multiscale, hierarchical modeling of chemical reactors
Gonzalez-Diaz et al. General theory for multiple input-output perturbations in complex molecular systems. 1. Linear QSPR electronegativity models in physical, organic, and medicinal chemistry
Coudert et al. Computational characterization and prediction of metal–organic framework properties
Vande Wouwer et al. Simulation of distributed parameter systems using a Matlab-based method of lines toolbox: chemical engineering applications
CN104573297A (zh) 一种确定表面催化反应路径的方法
Amsler et al. Anharmonic correction to adsorption free energy from DFT-based MD using thermodynamic integration
CN103793622A (zh) 一种分析反应中间体对催化剂活性影响的方法
Mondal et al. Thermal lattice expansion effect on reactive scattering of H2 from Cu (111) at T s= 925 K
Hickson et al. Quantum tunneling enhancement of the C+ H2O and C+ D2O reactions at low temperature
Liu et al. Ring polymer molecular dynamics in gas–surface reactions: Inclusion of quantum effects made simple
Lamberts et al. Formation of acetaldehyde on CO-rich ices
Göltl et al. Developing a descriptor-based approach for CO and NO adsorption strength to transition metal sites in zeolites
Zhou et al. Hybrid data-driven and mechanistic modeling approaches for multiscale material and process design
Leonardi et al. On the Maxwell− Stefan Approach to Diffusion: A General Resolution in the Transient Regime for One-Dimensional Systems
Iwasa et al. Combined automated reaction pathway searches and sparse modeling analysis for catalytic properties of lowest energy twins of Cu13
Shakouri et al. Analysis of energy dissipation channels in a benchmark system of activated dissociation: N2 on Ru (0001)
Subramanian Balashankar et al. Process optimization-based screening of zeolites for post-combustion CO2 capture by vacuum swing adsorption
Kim et al. Accelerated discovery of zeolite structures with superior mechanical properties via active learning
Horvatits et al. Characterizing adsorption sites on Ag/SSZ-13 zeolites: experimental observations and Bayesian inference
Chen et al. Dynamic and Intermediate-Specific Local Coverage Controls the Syngas Conversion on Rh (111) Surfaces: An Operando Theoretical Analysis
Gu et al. Short-and long-time dynamics of hydrogen spillover from a single atom platinum active site to the Cu (111) host surface
Zhang et al. First-principles study about the effect of coverage on H2 adsorption and dissociation over a Rh (100) surface
Gerrits Accurate simulations of the reaction of H2 on a curved Pt crystal through machine learning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant