Department of Chemistry, Kyung Hee University, 1 Seochun-Dong, Giheung-Gu, Yongin-Si, Gyeonggi-Do, 446-701, Korea.
J Chem Theory Comput. 2011 Mar 8;7(3):642-57. doi: 10.1021/ct100647b. Epub 2011 Jan 19.
A systematic study of the excited-state tautomerization of 7-azaindole-(H2O)n (n = 1 and 2) complexes in both gas and solution phases were investigated theoretically. Electronic structures and energies for the reactant, transition state (TS), and product were computed using the time-dependent density functional theory (TDDFT) and complete active space self-consistent field (CASSCF) levels with 6-31G (d,p), 6-311G(d,p), and 6-311+G(d,p) basis sets. Barrier heights and tautomerization energies were corrected by the second-order multireference perturbation theory (MRPT2) to consider the dynamic electron correlation. The solvent effect decreased the tautomerization barrier height in the 7-azaindole-H2O complex. In the 7-azaindole-(H2O)2 complex, two transition states were found for two asynchronous but concerted paths: in the first, the pyrole ring proton moved first to water; in the second, the water proton moved first to the pyridine ring. The CASSCF level with the MRPT2 correction clearly showed that the former path was much preferable to the latter. The preferable barrier height was only 1.6 kcal/mol with a zero-point energy correction, which would make the excited-state tautomerization possible. At all TDDFT levels, the TS structures and barrier heights depended on both the basis set used and the solvent effect. Most TDDFT methods failed to reproduce the CASSCF structures and MRPT2 energies. Only two methods, WB97XD/6-31G(d,p) and M062X/6-311+G(d,p), predicted two TSs for the two asynchronous paths in the 7AI-(H2O)2 complex but failed to reproduce the energetics. Further systematic study is necessary to test whether current TDDFT methods, including solvent effects, can be used to understand excited-state proton transfer reactions.
采用含时密度泛函理论(TDDFT)和完全活性空间自洽场(CASSCF)方法,并在 6-31G(d,p)、6-311G(d,p)和 6-311+G(d,p)基组上,对气相和溶液相中 7-氮杂吲哚-(H2O)n(n=1 和 2)复合物的激发态互变异构的反应历程进行了理论研究。用二级多参考微扰理论(MRPT2)对反应物、过渡态(TS)和产物的电子结构和能量进行了修正,以考虑动态电子相关。溶剂效应对 7-氮杂吲哚-H2O 复合物的互变异构势垒高度有降低作用。在 7-氮杂吲哚-(H2O)2 复合物中,发现了两条异步但协同的反应途径存在两个过渡态:在第一个途径中,吡咯环质子首先转移到水中;在第二个途径中,水合质子首先转移到吡啶环。带有 MRPT2 校正的 CASSCF 能级清楚地表明,前一条途径比后一条途径更可取。带有零点能校正的有利势垒高度仅为 1.6 kcal/mol,这使得激发态互变异构成为可能。在所有 TDDFT 水平上,TS 结构和势垒高度取决于所使用的基组和溶剂效应。大多数 TDDFT 方法都无法重现 CASSCF 结构和 MRPT2 能量。只有两种方法,WB97XD/6-31G(d,p)和 M062X/6-311+G(d,p),可以预测 7AI-(H2O)2 复合物中两条异步途径的两个 TS,但无法重现能量。需要进一步的系统研究来检验当前的 TDDFT 方法(包括溶剂效应)是否可以用于理解激发态质子转移反应。