Fe—PNP螯合物催化酮氢化反应机理的DFT研究
DFT研究C-H活化及Huisgen环加成反应机理
DFT研究C-H活化及Huisgen环加成反应机理DFT研究C-H活化及Huisgen环加成反应机理摘要：C-H活化反应是有机合成中的一种重要策略，能够直接将碳氢键转化为碳键，为有机合成提供了一种高效且环境友好的方法。
Huisgen环加成反应能够在无需催化剂的情况下实现合成1,2,3-三唑化合物，因此成为了有机合成中的重要工具。
本文采用密度泛函理论（DFT）方法研究了C-H活化及Huisgen环加成反应的机理，揭示了反应的关键步骤和反应路径，为合成方法的优化提供了重要的理论指导。
1. 引言C-H活化是一种将碳氢键转化为碳键的反应，具有高效、环境友好等优点，成为有机合成领域的热点研究方向。
Huisgen环加成反应是一种无需催化剂即可生成1,2,3-三唑化合物的重要反应，具有广泛的应用价值。
研究这两种反应的机理，可以深入了解反应的速率控制步骤及反应路径，从而指导实验研究，提高反应的效率。
2. 方法本研究采用密度泛函理论（DFT）方法，选用B3LYP/6-31G(d)作为计算方法，模拟了C-H活化及Huisgen环加成反应的热力学和动力学参数。
结合实验数据，揭示了反应机理中的关键步骤和反应路径。
3. 结果与讨论3.1 C-H活化机理C-H活化反应的机理可以分为氧化还原步骤、C-H活化步骤和反应产物生成步骤。
本研究以苯为例，探究了C-H活化反应的机理。
首先，在氧化还原步骤中，苯与氧分子发生反应生成苯氧自由基；其次，在C-H活化步骤中，苯氧自由基与金属催化剂反应，将苯环上的氢原子转化为碳氧键，得到苯酚产物；最后，在反应产物生成步骤中，苯酚与反应溶剂或其他反应物发生进一步反应生成目标产物。
3.2 Huisgen环加成机理Huisgen环加成反应可以分为亲核进攻步骤、环合步骤和负离子迁移步骤。
本研究以苯乙烯为例，探究了Huisgen环加成反应的机理。
首先，在亲核进攻步骤中，苯乙烯与亲核试剂进行反应，生成激发态的中间体；其次，在环合步骤中，中间体发生环合反应，生成1,2,3-三唑化合物；最后，在负离子迁移步骤中，负离子发生迁移，得到最终产物。
fe-n-c催化剂机理及研究
Fe-N-C催化剂机理及研究
Fe-N-C催化剂是一类非贵金属催化剂，由铁、氮和碳组成，具有重要的电催化性能，可用于燃料电池、水电解和电化学CO₂还原等领域。
其机理和研究主要包括以下几个方面：
1. 海绵铁的制备：通常通过加热有机铁盐（如铁酸铵）与碳源（如多孔炭材料或有机聚合物）反应，生成海绵状的Fe-N-C前驱体。
其中，碳源起到模板的作用，可决定催化剂的孔隙结构和表面活性。
2. 铁原子和氮原子的锚定和配位：在高温和气氛下，铁原子被嵌入到碳基质中，并经过一系列的反应与氨基衍生物发生配位，形成Fe-N-C结构。
氮原子的锚定不仅有助于稳定催化剂的活性位点，还可以调节催化剂的电子结构，增强其催化活性。
3. 催化活性位点：Fe-N-C催化剂的活性位点通常为铁原子与氮原子之间的键合结构，其中Fe-Nx键（x表示配位数）具有较高的电子密度和可调节的电子结构。
这些活性位点对氧还原反应、氢析出反应和CO₂还原反应等具有良好的催化活性。
4. 机理研究：为了深入理解Fe-N-C催化剂的电催化性能，研究者通常采用密度泛函理论（DFT）等计算方法，模拟催化剂表面的电子结构和反应机理。
此外，通过实验手段，如电化学原位技术、X射线吸收和散射等，还可以研究催化剂的结构演化、活性位点的形成和失活机制等。
目前，对Fe-N-C催化剂的研究主要集中于改进催化剂的活性位点和稳定性，提高其催化活性和耐久性，以实现其在能源转换和储存领域的应用。
ESI-MSn及DFT理论研究氧化脱肟、异喹啉酮及异吲哚啉酮合成的反应机理的开题报告
ESI-MSn及DFT理论研究氧化脱肟、异喹啉酮及异吲哚啉酮合成的反应机理的开题报告一、研究背景氧化脱肟、异喹啉酮及异吲哚啉酮合成是有机合成中常用的重要反应。
随着质谱、量子化学计算等技术的发展，研究反应机理越来越受到关注。
本文将采用质谱技术和DFT理论方法探究这些反应的反应机理，以期为有机合成的研究提供理论支持。
二、研究内容1. 介绍反应的背景和研究现状，概述本文的研究内容。
2. 介绍ESI-MSn技术及其在有机反应研究中的应用。
将采用该技术研究氧化脱肟、异喹啉酮及异吲哚啉酮合成的反应中生成的中间体以及转化过程，揭示反应机理。
3. 利用DFT计算方法对上述反应中的中间体和过渡态进行计算，进一步分析反应机理。
4. 根据实验结果和理论分析，总结反应机理及其特点，探讨反应条件、反应性质及反应特异性等问题。
5. 在总结研究成果的基础上，提出可行的优化方案，改进反应方法，提高反应效率和选择性。
三、研究意义1. 揭示氧化脱肟、异喹啉酮及异吲哚啉酮合成反应的反应机理，为合成具有特定结构和活性的有机化合物提供理论支持。
2. 探究反应条件、反应性质及反应特异性等问题，为优化反应条件提供可行方案，提高反应效率和选择性。
3. 对质谱技术和DFT计算方法在有机合成研究中的应用进行探索和应用，为有机合成研究提供新的思路和方法。
四、研究方法1. 采用高分辨质谱仪进行反应中间体和转化过程的分析，揭示反应机理。
2. 利用DFT方法计算反应中间体和过渡态的结构和能量，分析反应机理。
3. 结合实验和理论分析，总结反应机理及其特点，提出优化方案，改进反应方法。
五、预期结果1. 揭示氧化脱肟、异喹啉酮及异吲哚啉酮合成反应的反应机理。
2. 提出可行的反应优化方案，提高反应效率和选择性。
3. 探索质谱技术和DFT计算方法在有机合成研究中的应用。
无机化学催化反应的研究进展
无机化学催化反应的研究进展无机化学催化反应是一种高效的化学反应方式，具有广泛的应用前景。
在过去的几十年中，无机化学催化反应已经成为了化学领域中的重要研究方向。
本文将介绍无机化学催化反应的研究进展，包括反应机理、催化剂设计、反应体系等方面。
一、反应机理催化反应机理是无机化学催化反应研究的核心问题。
反应机理的解析对于催化剂的设计和反应条件的优化具有重要意义。
在反应机理的研究中，理论计算成为了一种重要工具。
目前，密度泛函理论（DFT）是计算无机化学催化反应机理的主要方法。
DFT 可以计算分子在特定区域的电子结构和分子间相互作用等相关性质，从而提供反应物和产物之间的能量差和反应中间体的产生和消失情况等反应机理。
二、催化剂设计催化剂是无机化学催化反应的关键因素。
随着对催化机理的理解逐渐加深，人们开始开发更加高效、选择性更好的催化剂。
常用的催化剂有金属催化剂、非金属催化剂、纳米催化剂等。
其中，金属催化剂具有丰富的反应底物适应性和反应选择性，在无机化学催化反应中起到了重要作用。
非金属催化剂具有较好的稳定性和环保性，已经成为了无机化学催化反应研究的热点之一。
三、反应体系反应体系包括反应底物、溶剂、温度、压力、ph值等因素。
在无机化学催化反应研究中，反应体系的优化对于反应的高效进行有着至关重要的作用。
例如，在氢化反应中，溶剂的选择和反应温度可以显著影响反应的选择性和活性。
此外，反应的酸碱度和ph值对于催化剂的稳定性和反应的选择性也有着很大的影响。
四、应用前景无机化学催化反应是一种高效的化学反应方式，在有机合成、环境治理、能源转化等领域具有广泛的应用前景。
例如，在有机合成中，无机催化剂可以用于各种酯、酰胺和醛烷化等反应中，可以获得高产率和高选择性。
在环境治理领域，氧化反应已经成为了一种有前途的污染物处理方式。
在能源转化领域，无机催化剂可以用于燃料电池、氢化制氢等过程中，可以达到高效转化和能源利用的目的。
综上所述，无机化学催化反应研究在理论和应用方面均具有很高的价值。
一种铁配合物催化的烯烃氢化反应机理的研究
一种铁配合物催化的烯烃氢化反应机理的研究作者：王方方赵英杰张洪铨来源：《山东工业技术》2016年第24期摘要：采用密度泛函理论（DFT）研究了铁配合物（MesCNC）Fe（N2）2催化的烯烃氢化反应的机理。
利用TPSSTPSS方法对反应体系中势能面上的各驻点的几何构型进行了全参数优化。
计算结果表明催化过程有两条可能的反应路径。
反应路径一：催化剂Cat先与乙烯配位再与H2分子配位生成中间体INT1，INT1中的两个氢原子发生均裂，生成含有乙基的化合物C。
反应路径二：Cat先与H2分子配位得到中间体INT2，INT2中的两个氢原子发生均裂分别与Fe成键生成化合物B，B与乙烯配位的同时一个氢原子转移到乙烯上生成化合物C。
C中的另一个氢原子转移到乙基上得到乙烷与Fe的配合物D，最后D脱去一分子的乙烷生成最初的Cat。
整个催化过程的决速步骤是化合物D的生成。
关键词：密度泛函理论；Fe的配合物；反应机理；烯烃氢化DOI：10.16640/ki.37-1222/t.2016.24.0290 引言烯烃加氢还原生成烷烃的反应在有机合成化学特别是精细化工，香料和医药合成中起着重要的作用。
催化氢化反应的催化剂主要以钌、铑、铱这类贵金属为活性中心[1-2]。
高成本，低含量以及由这些稀有金属所引起的环境污染促使人们寻找高效低毒及地球含量丰富的金属催化剂，而铁配合物成为很有应用前景的替代品。
虽然文献上很早就有Fe（CO）5及其衍生物催化烯烃氢化反应的相关报道，但是能在温和条件下表现出高活性和广泛底物范围的铁催化剂很少[3-6]。
近来，Chirik课题组合成了富电子的二咪唑吡啶铁双氮配合物（MeCNC）Fe（N2）2（MeCNC=2，6-（2，6-Me2-C6H3-imidazol-2-ylidene）2−C5H3N）及（MesCNC）Fe（N2）2（MesCNC=2，6-（2，4，6-Me3-C6H2-imidazol）2−C5H3N），这两种化合物在温和条件下（4 atm H2，23℃）对非官能团化的烯烃的氢化反应展现出良好的催化性能。
配合物的反应动力学及反应机理
SN1速度
减小
增加
(zēngjiā)
不影响
不影响
减小
增加 增加 增加
SN2速度
增加
增加
增加 减小 减小
减小 减小 减小
* 近似地适用于半满或全满d壳层情况，不适用于具有共价键或键配合物
共八十四页
2. 配体的取代(qǔdài)反应
2.3 八面体配合(pèihé)物的取代反应的速率方程
1. 离解(D)机理
共八十四页
1. 基本概念
1.6 势能(shìnéng)曲线 活化（自由）能
过渡态：瞬间存在的不稳定物种；
中间体：相对稳定的过渡(guòdù)化合物，在慢反应中能检测其存在；
共八十四页
1. 基本概念
1.7 活化(huóhuà)参数
通常反应机理(jī lǐ)不能够完全由速率方程决定，还需从实验获得的活化参数予以证实。
例：[Fe(H2O)6]3+ + Y 比 [Fe(H2O)5(OH)]2+ + Y 反应速率大2个数量级
反应速率只与水合金属离子溶度有关，与外来配体L无关。
反应机理：
r k[M(H2O)x ]
[M(H2O)x ]n slow[M(H2O)x1]n H2O [M(H2O)x1]n Lm fast[M(H2O)x1L]nm
键的断裂是反应速率的决定步骤
共八十四页
2. 配体的取代(qǔdài)反应
（2）缔合（置换(zhìhuàn)）机理（A机理，SN2机理）
1. 中心原子的电子结构
d0-d2电子：有空的d轨道，易接受外来配体的电子，生成 配位数为7 的活化配合物。
Cr3 (d3 ) :
d3-d6电子：惰性。（需要较高活化能）
Fe—PNP螯合物催化酮氢化反应机理的DFT研究
主要 有三个步骤 ： （ １ ） 负氢转移过程 ； （ ２ ） 一分子 氢气异裂加成过程 ； （ ３ ） 负氢转移 同时 Ｆ ｅ － Ｏ断裂生成醇并使催化 剂再生过程。整个催化过程 决
速步骤 是氢异 裂加 氢过程 ． 氢异裂加氢反应过程的活化能 垒为 ５ ２ ． ６ ｋ Ｊ ｍｏ ｌ － Ｉ 放 出热量 ７ ６ ． ２ １ ｄ ｏ ｔ ｏ ｌ ～ 。
【 摘 要】 本文采 用密度泛函理 论（ Ｄ 兀） ， 对Ｆ ｅ — Ｐ Ｎ Ｐ螯合物催化 酮氢化反应机理作 了深入的理论研究。 利用密度涟函理论 中的 Ｂ ３ Ｌ Ｙ Ｐ方法 对反应物 、 产物 、 过渡 态（ １ 、 ｓ ） 的几何构 型进行全参数优化 ， 并通过振动频 率分析确认各稳定点和过渡 态的真 实性。计 算结果表明 ， 整个催化循环
（ ２ ）一分子氢气异裂加到配合物 Ｃ上加经过渡态 Ｔ Ｓ Ｃ — Ｄ生成配合物
Ｄ ： （ ３ ） 随后配合物 Ｄ碳上 的氢进攻羰基氧同时 Ｆ ｅ — Ｏ断开生成醇水平 上优化得到 的各物质 的几何构型如
图３ 所示 ． 相对能量曲线见图 ４
２ ． １ 催 化剂和底物 的结构 首先 对模型催 化剂 Ｆ ｅ — Ｐ Ｎ Ｐ 螫 合物和底 物苯 乙酮在 Ｂ ３ Ｌ Ｙ Ｐ ／ Ｂ Ｓ Ｉ
化过程主要有 三个步骤 ： （ １ ） Ｆ ｅ — Ｐ Ｎ Ｐ 螫合物脱溴与苯乙酮配位生成 Ｂ 后Ｆ ｅ 上 的 Ｈ原子经 过渡态 Ｔ Ｓ Ｂ — Ｃ转移 到羰基碳 上生成 配合物 ｃ ：
dft计算有机反应机理
dft计算有机反应机理全文共四篇示例，供读者参考第一篇示例：近年来，DFT（密度泛函理论）计算技术在有机化学领域中得到了越来越广泛的应用，特别是在研究有机反应的机理方面。
DFT是一种基于量子力学的计算方法，能够准确地描述分子的结构、反应活性和能量。
在有机反应机理研究中，通过DFT计算可以揭示反应过程中的中间体、过渡态和活化能等关键参数，从而帮助化学家更好地理解反应的机理和动力学过程。
DFT计算在有机反应机理研究中的应用主要包括以下几个方面：1. 结构优化：DFT可以通过优化分子的几何构型，确定最稳定的分子结构和构型。
这对于研究反应物、产物和中间体的结构及构型变化十分重要，可以帮助研究人员预测反应的可能路径和产物。
2. 过渡态搜索：DFT计算在寻找反应过渡态的过程中也发挥着重要作用。
过渡态是反应活化能的关键参数，通过计算过渡态的结构和活化能可以预测反应的速率常数和反应速率。
DFT计算可以帮助研究人员快速、准确地找到反应过渡态，从而揭示反应机理。
3. 活化能计算：DFT计算可以准确地计算反应物转变为过渡态的能量，即活化能。
活化能是一个反应进行的临界能量，通过DFT计算可以帮助研究人员评估反应的难易程度和速率常数。
这对于设计新的催化剂和优化反应条件具有重要意义。
4. 分子轨道分析：DFT计算还可以对分子轨道进行分析，分析分子的电子结构和反应机制。
通过分析分子轨道的相互作用，可以揭示反应的电子转移过程，帮助揭示反应的原子层面机理。
DFT计算在有机反应机理研究中的应用为研究人员提供了一种准确、高效和可靠的手段，可以帮助他们揭示反应机理、预测反应路径和优化反应条件。
随着计算技术的不断发展和改进，DFT计算在有机反应机理研究中的应用前景将更加广阔，为有机化学领域的发展提供更多可能性和机遇。
第二篇示例：DFT计算在有机反应机理研究中的应用一、DFT计算简介DFT是一种用于计算电子结构和化学性质的理论方法。
它基于量子力学原理，通过求解电子的薛定谔方程来描述分子的结构和性质。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Fe—PNP螯合物催化酮氢化反应机理的DFT研究
【摘要】本文采用密度泛函理论（DFT），对Fe-PNP螯合物催化酮氢化反应机理作了深入的理论研究。
利用密度泛函理论中的B3LYP方法对反应物、产物、过渡态（TS）的几何构型进行全参数优化，并通过振动频率分析确认各稳定点和过渡态的真实性。
计算结果表明，整个催化循环主要有三个步骤：（1）负氢转移过程；（2）一分子氢气异裂加成过程；（3）负氢转移同时Fe-O断裂生成醇并使催化剂再生过程。
整个催化过程决速步骤是氢异裂加氢过程，氢异裂加氢反应过程的活化能垒为52.6kJ mol-1，放出热量76.2kJ mol-1。
【关键词】构型优化；密度泛函理论；Fe-PNP螯合物；反应机理
0 引言
酮加氢还原生成醇的反应在有机合成化学特别是精细化工，香料和医药合成中起着重要的作用。
目前，催化氢化反应的高效双官能团催化剂主要以钌、铑、铱这类贵金属为活性中心，稀有，价格昂贵，并且具有一定的毒性。
有相似活性的铁配合物成为近年人们研究的热点，铁是地球上含量最丰富的过渡金属，其毒性低，价格低，铁配合物对环境友好等，其应用前景宽广[1，2]。
以铁为活性组分的催化剂在最近取得了长足的进展，但在铁催化还原酮领域的研究报道还比较少[3，4]。
最近Milstein 等合成出一种新的催化剂[Fe（Br）（CO）（H）（PNP-iPr）]，他们通过实验方法对催化反应的机理进行了初步探索，推出了反应的中间体，但催化过程过渡态不明确[5]。
本文拟采用理论方法对该反应机理进行分析，以期为以后研究铁配合物的催化提供较好的参考和理论指导。
1 研究方法
利用密度泛函理论（DFT）中的B3LYP[6-8]方法对反应物、产物、过渡态（TS）的几何构型进行全参数优化，其中对Fe原子采用LANL2DZ[9]基组，其它原子用6-31+G（d，p）基组进行全面优化，这两个基组的组合记做BS1。
然后通过振动频率分析确认各稳定点（没有虚频）和过渡态（有且仅有一个虚频）的真实性，得到各驻点的零点能（ZPE）。
除非另有说明，本文中报道的能量值都是在B3LYP/BS1水平上获得的电子能量加上零点振动能（ZPE），以kJ mol-1为单位。
以上所有的计算均由Gaussian 03程序[10]完成。
2 结果与讨论
2.1 催化剂和底物的结构
首先对模型催化剂Fe-PNP螯合物和底物苯乙酮在B3LYP/BSI水平上进行结构优化和频率分析，优化得到的Fe-PNP螯合物和底物苯乙酮的几何构型和结构参数如图1所示。
将优化得到的Fe-PNP螯合物构型与实验上通过X-射线分析得到的Fe-iPr PNP晶体结构参数进行了对比，发现把异丙基简化为甲基计算得到的
几何参数与原物质的X-ray分析结果吻和的很好，可见对催化剂的简化和采用的计算方法是可行和可靠的。
底物的几何构型（键长单位为？魡）
2.2 催化循环过程
根据实验研究和计算的结果得到的催化反应机理如图2。
整个催化过程主要有三个步骤：（1）Fe-PNP螯合物脱溴与苯乙酮配位生成B后Fe上的H原子经过渡态TSB-C转移到羰基碳上生成配合物C；（2）一分子氢气异裂加到配合物C 上加经过渡态TSC-D生成配合物D；（3）随后配合物D碳上的氢进攻羰基氧同时Fe-O断开生成醇并使催化剂再生。
在B3LYP/BS1水平上优化得到的各物质的几何构型如图3所示，相对能量曲线见图4。
2.2.1 催化剂与酮配位后负氢转移过程
配合物A与苯乙酮配位生成B后铁上的氢经过渡态TSB-C转移到羰基碳上生成配合物C。
在TSB-C中Fe-O1和H-C1的键长分别为2.091？魡和2.025？魡（图3）。
频率分析显示TSB-C仅有一个193.22i cm-1的虚频，主要是H1由Fe 向C1的振动。
由B经TSB-C生成C反应的活化能垒为31.9kJ mol-1，放出热量65.8 kJ mol-1（见图4）。
几何构型（键长单位为？魡）
2.2.2 氢异裂过程
一分子氢异裂后经由过渡态TSC-D分别加到C1和Fe上生成配合物D。
其中过渡态TSC-D中H1-H2和C1-H2的键长分别是0.908？魡和1.652？魡（图3），频率分析显示TSC-D中仅有一个虚频756.62i cm-1，主要是H1和H2分别向Fe 和C1的振动。
氢异裂加成反应过程的活化能垒为52.6kJ mol-1，放出热量76.2kJ mol-1。
2.2.3 催化剂再生过程
配合物D经过渡态TSD-A H1转移到O1上同时Fe-O1断开生成醇同时催化剂再生。
其中过渡态TSD-A中H1-O1和H1-C1的键长分别是1.205？魡和1.462？魡（图3），频率分析显示TSD-A中仅有一个虚频1051.80i cm-1，主要是H1由C1向O1的振动。
该反应过程的活化能垒为47.3kJ mol-1，吸收热量39.9kJ mol-1。
根据计算结果可知整个反应过程的绝速步骤是氢异裂加成过程根据下面的方程（基于过渡态理论）可以估算出反应速率常数：
K=■exp（-ΔG/RT）
其中h是普朗克常数其值为 6.626×10-34 Js，kB是波尔兹曼常数其值为1.381×10-23 J K-1，R是气体常数其值为R=8.314J mol-1 K-1，根据计算结果可以算出反应速率常数为2.5h-1。
3 结论
采用密度泛函方法研究了Fe-PNP螯合物催化酮氢化反应机理。
根据计算结果，整个催化过程由Fe-PNP配合物失去Br后与苯乙酮配合开始，主要包括如下步骤：（1）Fe-PNP螯合物脱溴与苯乙酮配位后生成的B上Fe的H原子经过渡态TSB-C转移到羰基碳上生成配合物C；（2）一分子氢气异裂加到配合物C 上加经过渡态TSC-D生成配合物D；（3）随后配合物D碳上的氢进攻羰基氧同时Fe-O断开生成醇并使催化剂再生。
对整个催化循环来说，氢异裂加成过程是整个催化过程的决速步骤，活化能为52.6 kJ mol-1。
催化反应的速率常数为2.5h-1。
【参考文献】
[1]S. Enthaler，K. Junge，M. Beller. Sustainable Metal Catalysis with Iron：From Rust to a Rising Star？[J]. Angew. Chem. Int. Ed.，2008，47（18）：3317-3321.
[2]C.Bolm，J. Legros，J. Le Paih，L. Zani. Iron-Catalyzed Reactions in Organic Synthesis[J]. Chem. Rev. 2004，104（12）：6217-6254.
[3]C. P. Casey and H. Guan. An Efficient and Chemoselective Iron Catalyst for the Hydrogenation of Ketones[J]. J. Am. Chem. Soc.，2007，129（18）：5816-5817.
[4]C. Sui-Seng，F. N. Haque，A. Hadzovic，A. Ptz，V. Reuss，N. Meyer，
A. J. Lough，M. Zimmer-De Iuliis，R. H. Morris. Synthesis and Characterization of Iron（II）Complexes with Tetradentate Diiminodiphosphine or Diaminodiphosphine Ligands as Precatalysts for the Hydrogenation of Acetophenone[J]. Inorg. Chem.，2009，48（2）：735-743.
[5]Robert Langer，Gregory Leitus，Yehoshoa Ben-David，and David Milstein. Efficient Hydrogenation of Ketones Catalyzed by an Iron Pincer Complex[J]. Angew. Chem.，2011，123（9）：2168-2172.
[6]R. G. Parr，W. Yang. Density-Functional Theory of Atoms and Molecules[M]. Oxford University press：New York，1989：47-66.
[7]A. D. Becke. Density-functional thermochemistry. III. The role of exact exchange[J]. J. Chem. Phys.，1993，98（7）：5648-5652.
[8]A. D. Becke. Density-functional exchange-energy approximation with correct asymptotic behavior[J]. Phys. Rev. A，1988，38（8）：3098-3100.
[9]C. Lee，W. Yang and R. G. Parr. Development of the Colle-Salvetti
correlation-energy formula into a functional of the electron density[J]. Phys. Rev. B，1988，37（2）：785-789.
[10]M.J.Frisch et al.，Gaussian 03，Revision B05，Gaussian，Inc.，Pittsburgh PA，2003.。