烷烃与多种自由基（如Cl、O(3P)、OH、NH和CN）在低温区发生反应时，速率常数与温度呈现出特殊的负相关性，即随着温度的升高反应速率常数反而降低。
为了解释这一反常现象，美国能源部Sandia国家实验室的Georgievskii等人借助CASPT2、QCISD(T)等高精度量化方法对氰自由基（CN）与乙烷分子的反应机理进行研究，并分析了这类反应的动力学性质。研究证实，在低温区（<200k）这类自由基反应体系的决速步骤是通过弱相互作用（范德华力）形成复合物，而在高温区的决速步骤则涉及到化学键的断裂和生成。通过仔细考察反应过程中涉及的两种过渡态机制，Georgievskii等人提出一种双过渡态模型，其理论计算的结果准确地重现了实验中观察到的特殊温度效应曲线。
图1、CASPT2方法优化得到的氰自由基与
乙烷分子反应的过渡态结构图2、双过渡态模型与其它过渡态模型的结果对比
Ref: Y. Georgievskii; S. J. Klippenstein, J. Phys. Chem. A, 2007, 111, 3802-3811.
使用软件：Molpro
近年来叠氮自由基（N3）与氟原子（F）的反应机理得到了广泛的关注，该反应体系在新型激光器的研制方面有着潜在的用途。
中科院化学所边文生等人用CASSCF和MRCI方法对N3+F反应体系进行研究，通过扫描单重态和三重态的势能曲线确定了两态之间的最低能量交叉点，从而澄清了叠氮自由基与氟原子的反应机制。计算结果表明，该反应体系的单重态和三重态之间存在着绝热和非绝热两种转化方式，其中绝热反应方式（图1, 式i）将占主导优势，这与实验观测到的现象相吻合。
图1、氟原子与叠氮自由基反应的两种产物
图2、产物NF分子的单重态和三重态势能曲线
图3、叠氮自由基N3的单重态和三重态势能曲线
Ref: H. T. Ma; X. J. Liu; W. S. Bian; L. P. Meng; S. J. Zheng, Chemphyschem, 2006, 7, 1786-1794.
使用软件：Molpro
QM/MM方法混合使用量子力学和分子力学，能够显著降低计算量，近年来已成为模拟酶催化反应的一种重要手段。由于酶催化反应涉及的体系通常比较大，即使采用了QM/MM方法，需要用量子力学处理的原子数仍然比较多，因而目前多采用精度相对较低的计算方法，如半经验方法或密度泛函理论（DFT）。在处理较大分子体系的时候半经验方法是一个不错的选择，但是其精度相对较低，计算误差有可能达到十几个kcal/mol-1甚至更高。相对而言，DFT方法能够在一定程度上提高准确性，但是对于色散力等关键的弱相互作用仍然难以处理，因此DFT方法通常会低估反应能垒约几个千卡。
图1、对羟苯甲酸羟化酶PHBH的过渡态结构 图2、不同计算方法扫描得到的CM体系的势能曲线对比
最近，英国布里斯托大学的Mulholland、德国马普协会煤化所的Thiel等人采用QM/MM方式对分支酸变位酶（chorismate mutase，CM）和对羟苯甲酸羟化酶（para-hydroxybenzoate hydroxylase，PHBH）在溶液环境中的反应过程进行模拟，计算过程中这两种酶催化体系分别含有7千和2万3千个原子。他们用Molpro软件的B3LYP、LMP2和LCCSD(T0)等高精度计算方法处理QM部分的原子，成功地预测了这两种酶在催化反应过程中的活化焓和自由能。计算结果表明，LCCSD(T0)的计算结果最为准确，其预测值与实验观测值几乎完全吻合，而B3LYP和LMP2方法在不同程度上都低估了酶催化反应的活化能垒。
Ref: F. Claeyssens; J. N. Harvey; F. R. Manby; R. A. Mata; A. J. Mulholland; K. E. Ranaghan; M. Schutz; S. Thiel; W. Thiel; H. J. Werner, Angew. Chem. Int. Ed., 2006, 45, 6856-6859.
使用软件：Molpro
树枝状聚苯乙炔（Phenylacetylene Dendrimer，PAD）是一种带有多个发色团的支化大分子，在光照激发之后各个聚合单元之间能够迅速而有效地传递能量。传统的理论认为间位连接的苯乙炔单元会阻塞电子态的离域，因此PAD的能量传递将被局限在树枝状分子的分枝上。
美国伊利诺斯大学的Martinez和Bardeen等人认为这种理论仅仅考虑了基态分子的电子结构，而对于激发态的PAD来说，情况可能会有很大差异。为此他们采用CASSCF和CASPT2等方法对PAD分子进行研究，结果表明，尽管PAD各单元之间在处于基态时耦合作用非常微弱，但是在处于激发态时由于构型的变化耦合作用将会大幅度增强。在此基础上，Martinez和Bardeen等人对传统的Forster能量传递模型提出了修正建议。
图1、三乙炔基苯处于激发态时的最优构型 图2、PAD单元的激发、驰豫和发射示意图
图3、苯乙炔单元在吸收、发射状态的键长和跃迁偶极矩
Ref: K. M. Gaab; A. L. Thompson; J. J. Xu; T. J. Martinez; C. J. Bardeen, J. Am. Chem. Soc., 2003, 125, 9288-9289.
使用软件：Molpro
众所周知，由于氮气的氢化反应具有非常高的活化能垒，在没有催化剂存在时，这个反应只有在极端条件下才可能发生。目前人们尚未完全掌握氮气加氢的详细反应机制，仅知道氮气分子与第一个氢分子的反应是整个氢化过程中的决速步骤，因此，葡萄牙科英布拉大学Varandas等人利用MCSCF/aug-cc-pVQZ方法对偶氮烯（N2H2）分子进行高精度计算，详细地分析了N2H2分子的cis、trans和异构模型的最优构象，并计算了这几种构型之间相互转变的过渡态能垒，以求获得精确而全面的基态势能曲线，从而有助于了解氮气氢化过程的反应机理。
 |
 |
| 图1.氮气氢化过程的两种反应方式 | 图2.偶氮烯三种构型的分解反应路径 |
 |
| 图3.偶氮烯的cis、trans和iso三种构型的最优构象，以及三者之间相互转变的过渡态结构 |
Ref: M. Biczysko; L. A. Poveda; AJC Varandas, Chem. Phys. Lett., 2006, 424, 46-53.
使用软件：Molpro