1-甲基胸腺嘧啶A-带结构动力学

获取了1-甲基胸腺嘧啶(MT)涵盖紫外光谱中A带和B带吸收的共5 个激发波长的共振拉曼光谱, 并结合密度泛函理论方法研究了MT的电子激发和Franck-Condon 区域结构动力学. 在TD-B3LYP/6-311++G(d,p)计算水平下, A带和B带吸收被分别指认为π_H→π_L^*/π_H-2→π_L+2^*和π_H→π_L+2^→/π_H-2→π_L^*跃迁. 甲基参与嘧啶环的共轭使MT的A带最大吸收波长λ_max相对于胸腺嘧啶(T)发生明显红移, 并对Franck-Condon区域的动态结构产生一定影响. A带和B带共振拉曼光谱分别被指认为14 个振动模式和11 个振动模式的基频、泛频和组合频. C5＝C6伸缩+C6H12面内弯曲振动v₉, 环变形振动v₁₆和N3C2N1反对称伸缩+C4C5C10反对称伸缩振动v₁₈占据了A带共振拉曼光谱强度的绝大部分. 这表明¹π_Hπ_L^*激发态结构动力学主要沿这些反应坐标展开. 考察了溶剂对共振拉曼光谱的影响, 结果表明, C4＝O9伸缩+N3H11面内弯曲振动v₈的活性与溶剂性质有关, 其激发态位移量随溶剂性质的变化规律与胸腺嘧啶一致.     

6-N,N-二甲基腺嘌呤激发态结构动力学的共振拉曼光谱

采用共振拉曼光谱技术和密度泛函理论方法研究了6-N,N-二甲基腺嘌呤(DMA)的A带和B带电子激发和Franck-Condon 区域结构动力学. π_H→π_L^*跃迁是A带吸收的主体, 其振子强度约占整个A带吸收的79%.由弥散轨道参与的n→Ryd 和π_H→Ryd 跃迁在B带跃迁中扮演重要角色, 其振子强度约占B带吸收的62%,而在A带吸收中占主导的π_H→π_L^*跃迁的振子强度在B带吸收中仅占33%. 嘌呤环变形伸缩+C8H/N9H面内弯曲振动ν₂₃和五元环变形伸缩+C8H弯曲振动ν₁₃的基频、泛频和合频占据了A带共振拉曼光谱强度的绝大部分, 说明¹π_Hπ_L^*激发态结构动力学主要沿嘌呤环的变形伸缩振动, N9H/C8H/C2H弯曲振动等反应坐标展开, 而ν₁₀, ν₂₉, ν₂₁, ν₂₆和ν₄₀的基频、泛频和合频占据了B带共振拉曼光谱强度的主体部分, 它们决定了B带激发态的结构动力学. A带共振拉曼光谱中ν₂₆和ν₁₂被认为与1nπ*/1ππ*势能面锥型交叉有关. B带共振拉曼光谱中ν₂₁的激活与¹ππ^*/¹πσ_N9H^*势能面锥型交叉相关.     

二聚体C2H4-nFn···LiH (n=0, 1, 2)中π锂键的性质和弱方向性

采用二级微扰理论(MP2)量子化学研究方法, 对C₂H_4-nF_n···LiH (n=0, 1, 2)二聚体的结构和π锂键性质进行了分析. 结果表明氟原子的取代改变了乙烯分子的π电子云形状, 从而使二聚体体系中的π锂键发生偏移、伸长和弯曲. 通过与类似的π氢键体系C₂H_4-nF_n···HF (n=0, 1, 2)比较, 发现π锂键在二级弱相互作用的影响下, 发生了明显的弯曲, 表现出弱的方向性. 在CCSD(T)/6-311++G(3df, 3pd)理论水平下, 二聚体的相互作用能强弱顺序为: 33.85 kJ·mol^-1 (C₂H₄-LiH)>27.32 kJ·mol^-1 (C₂H₃F-LiH)>21.34 kJ·mol^-1 (cis-C₂H₂F₂-LiH)>20.25 kJ·mol^-1 (g-C₂H₂F₂-LiH), 说明氟取代削减了乙烯分子与LiH之间的相互作用.     

AuCl3催化2-(1-炔基)-2-烯基酮合成多取代呋喃的机理

采用密度泛函理论的B3LYP泛函对AuCl₃催化的2-(1-炔基)-2-烯基酮与亲核试剂反应的机理进行了研究, 得到了反应的最优路径. 结果表明, 整个反应的决速步骤是羟基H转移到AuCl₃的配体Cl上, 其活化能为49.3 kJ·mol^-1. 通过计算发现, 催化剂AuCl₃的配体Cl原子在反应中有重要的作用, 它不仅稳定配合物, 而且直接参与反应, 协助质子的转移, 显著降低质子转移的活化能(由71.5 kJ·mol^-1降低到49.3 kJ·mol^-1). 另外还讨论了HBF₄不能催化此反应的可能原因, 计算结果与实验结果一致.     

统计方法校正O3LYP方法计算的生成热

由于引入各种内在近似, 密度泛函理论存在固有误差. 本文采用O3LYP/6-311+G(3df, 2p)//O3LYP/6-31G(d)计算了220个中小型有机分子的生成热(Δ_fH_calc^Θ), 随后应用神经网络(ANN)和多元线性回归(MLR)方法对Δ_fH_calc^Θ进行校正. 采用计算得到的生成热、零点能、分子中原子总数、氢原子个数、双中心成键电子数、双中心反键电子数、单中心价层孤对电子数、单中心内层电子数作为ANN和MLR的描述符. 以180个分子作为训练集构造ANN或MLR模型, 并对40 个独立测试集分子的Δ_fH_calc^Θ进行了预测. 结果表明: 经过ANN和MLR校正后,训练集分子生成热的理论计算值和实验值间的均方根偏差(RMSD)从24.7 kJ·mol^-1分别降低到11.8、13.0 kJ·mol^-1; 独立测试集分子的RMSD从21.3 kJ·mol^-1分别降低到10.4、12.1 kJ·mol^-1. 因此ANN模型的拟合和预测能力要明显优于MLR模型.     

金(I)配合物催化乙烯加氢的反应机理

采用密度泛函理论B3LYP方法, 对两类金(I)配合物AuX (X=F, Cl, Br, I)和AuPR₃⁺(R=F, Cl, Br, I, H, Me,Ph)催化C₂H₄加氢反应的机理进行了理论研究. 计算显示Au(I)配合物对C₂H₄氢化具有较好的催化效果, 其作用下的加氢反应存在“活化H―H键后再与C₂H₄反应”和“活化C＝C键后再与H₂反应”两种途径, 前者的活化能较后者低90-120 kJ·mol^-1, 因而具有明显的能量优势. 研究表明AuPR₃⁺ 的催化能力明显强于AuX. 此外, X/PR3基团供、吸电子能力的变化对配合物的催化能力也具有较为显著的影响. 电子结构分析显示Au(I)配合物在C₂H₄ 加氢反应中不仅能够削弱H―H、C＝C 键的强度, 还使H₂ σ_H―H^*、C₂H₄ π_C＝C^* 轨道能级下降, 从而缩小了π_C＝C-σ_H―H^*或σ_H―H-π_C＝C^*轨道间的能级差, 促进了C₂H₄-H₂反应中的电子离域, 从而降低禁阻反应发生的难度.σ_H―H^*、π_C＝C^*轨道能级改变量与加氢反应活化能E_a的降低值之间存在较好的一致性关系, 因此使上述轨道能级下降幅度越大的Au(I)配合物可以获得较好的催化效果.     

KI/NH3协同催化CO2与环氧丙烷合成碳酸丙烯酯的反应机理

采用密度泛函理论(DFT)研究了NH₃/KI、KI及无催化剂条件下, CO₂与环氧丙烷(PO)合成碳酸丙烯酯(PC)的反应机理. 在B3LYP/6-311++G^**基组水平上(I采用MIDIX基组)优化了反应过程中的反应物、中间体、过渡态和产物, 通过振动分析及内禀反应坐标(IRC)确定中间体和过渡态的真实性. 同时, 在相同基组水平应用自然键轨道(NBO)理论和分子中的原子(AIM)理论分析了这些化合物的轨道间相互作用和成键特征. 研究结果表明: 在无催化剂条件下非协同反应通道PO+CO₂→M0a→TS0c→M0c→TS0c′→PC为最有利通道, 其活化能为200.65 kJ·mol^-1; KI催化下活化能降低至187.40 kJ·mol^-1, 反应速率较小; 在KI/NH₃协同催化下, 除KI的催化作用外, NH₃中的氢原子还能与CO₂或PO中的氧原子形成氢键, 活化反应, 活化能降低至154.64 kJ·mol^-1, 大幅度提高了PO与CO₂环合生成产物PC的反应速率, 理论计算与实验结果一致.     

糠醛在Pt(111)表面的吸附和脱碳反应

采用广义梯度近似的密度泛函理论并结合平板模型的方法, 优化了糠醛分子在Pt(111)面的吸附模型,并探究了糠醛脱碳反应形成呋喃的机理. 结果表明: 吸附后糠醛分子环上的C―H(O)键及支链―CHO相对于金属表面倾斜上翘, 分子平面被扭曲, 易于呋喃的形成; 同时, 糠醛分子向Pt表面转移电子0.765e, 环中的大π键与Pt(111)表面的d轨道发生较强的相互作用, 使得糠醛的芳香性被破坏, 环上的碳原子呈现准sp³杂化. 此外, 对糠醛脱碳反应中的各反应步骤进行过渡态搜索, 通过比较各步骤的活化能, 得出糠醛更易先失去支链上的H形成酰基中间体(C₄H₃O)CO, 中间体继续脱碳加氢形成产物呋喃. 该过程的控速步骤为(C₄H₃O)CO^*+^*→C₄H₃O^*+CO^* (^*为吸附位),活化能为127.65 kJ·mol^-1.     

8-羟基喹啉银(铂)金属配合物电子光谱与非线性光学性质

运用密度泛函理论B3LYP方法对8-羟基喹啉(银、铂)(AgQ、PtQ₂)金属配合物及其衍生物的非线性光学性质进行理论计算研究. 结果表明: 引入取代基使铂配合物的最强吸收波长产生较大红移. 最低能量跃迁吸收来自最高占据分子轨道(HOMO)到最低空分子轨道(LUMO)的d→π^*和π→π^*跃迁, 属于金属配体电荷转移(MLCT)与配体配体电荷转移(LLCT). 金属银和铂掺杂8-羟基喹啉使其三阶非线性光学系数γ值明显增大, 并且在配合物上引入―Ph, ―PhOCH₃, ―PhF₂, ―PhF₅基团将进一步增大γ值. 引入基团的供电子性越强, γ值增大的幅度越大, 引入基团的吸电子性越强, γ值增大的幅度越小.     

C2H3与CH3F氢抽提反应机理与动力学性质

采用双水平直接动力学方法对C₂H₃与CH₃F氢抽提反应进行了研究. 在QCISD(T)/6-311++G(d, p)//B3LYP/6-311G(d, p)水平上, 计算的三个反应通道R1、R2和R3的能垒(ΔE^≠)分别为43.2、43.9和44.1 kJ·mol^-1, 反应热为-38.2 kJ·mol^-1. 此外, 利用传统过渡态理论(TST)、正则变分过渡态理论(CVT)和包含小曲率隧道效应(SCT)的CVT, 分别计算了200-3000 K温度范围内反应的速率常数k^TST、k^CVT和k^CVT/SCT. 结果表明: (1) 三个氢抽提反应通道的速率常数随温度的增加而增大, 其中变分效应的影响可以忽略, 隧道效应则在低温段影响显著; (2) R1反应是主反应通道, 但随着温度的升高, R2反应的竞争力增大, 而R3反应对总速率常数的影响很小.     

2-乙酰基-1-甲基吡咯激发态结构动力学及溶剂效应的共振拉曼光谱和密度泛函理论计算

获取了覆盖紫外光谱中A带和B带吸收的共7个不同激发波长的共振拉曼光谱, 并结合密度泛函理论方法研究了2-乙酰基-1-甲基吡咯(2-Ac-NMP)的A带和B带电子激发和Franck-Condon区域结构动力学. 在TD-B3LYP/6-311++G(d,p)计算水平上, A带和B带吸收的跃迁主体为π→π^* . A带和B带共振拉曼光谱分别指认为13个振动模式和8个振动模式的基频、泛频和组合频, 其中C＝O伸缩振动(ν₈)、C3-C4-C5不对称伸缩振动+C2-C6伸缩振动(ν₁₄)及环上CH面内摇摆(ν₁₈)对拉曼光谱强度贡献最大, 表明2-Ac-NMP的S_π激发态结构动力学主要沿反应坐标展开. 考察了溶剂对共振拉曼光谱强度模式的影响, 结果表明, 在同一溶剂中, 随激发波长由长变短, C＝O伸缩振动模(ν₈)的强度呈现出由强变弱再变强的现象. 这种变化规律与Franck-Condon区域S_n/S_π态混合或势能面交叉相关, 并受溶剂的有效调控.     

N-甲基吡咯-2-甲醛激发态结构动力学及其溶剂效应的共振拉曼光谱和密度泛函理论研究

获取了覆盖N-甲基吡咯-2-甲醛(NMPCA)A-带和B-带电子吸收共7个激发波长的共振拉曼光谱,并结合含时密度泛函理论(TD-DFT)方法研究了的A-带和B-带电子激发和Franck-Condon区域结构动力学.TD-B3LYP/6-311++G(d,p)计算表明:A-带和B-带电子吸收的跃迁主体为π→π*.共振拉曼光谱可以指认为,11-13振动模式(A-带激发)或者7-11振动模式(B-带激发)的基频、倍频和组合频,其中C=O伸缩振动(ν7)、环的变形振动+N1-C6伸缩振动(ν17)、环的变形振动(ν21)和C6-N1-C2/C2-C3-C4不对称伸缩振动(ν14)占据了绝大部分.这表明NMPCA的Sπ激发态结构动力学主要沿C=O伸缩振动、环的变形振动和环上N1-C6伸缩振动等反应坐标展开.在同一溶剂的共振拉曼光谱中随激发波长由长变短,ν7与ν14的强度比呈现出由强变弱再变强的现象,这种变化规律被认为与Franck-Condon区域Sn/Sπ态混合或势能面交叉有关.溶剂对Sn/Sπ态混合或势能面交叉具有调控作用.     

HOCl…HCOCl复合物的结构和电子性质

在DFT-B3LYP/6-311++G**水平上求得HOCl+HCOCl复合物势能面上的四种稳定构型(S1,S2,S3和S4).其中,在复合物S1和S3中,HOCl单体的5H原子作为质子供体,与HCOCl单体中作为质子受体的10原子相互作用,形成红移氢键复合物;在复合物S4中,HOCl单体的7Cl原子作为质子供体,与HCOCl单体中作为质子受体的IO原子相互作用,形成红移卤键复合物;而在复合物S2中,同时存在2C-3H…6O蓝移氢键和4Cl…5O相互作用.在MP2/6-311++G**水平上计算的单体间的相互作用能考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正,其值在-5.05与-14.76 kJ·mol-1之间.采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

Circular dichroism, optical rotation and absolute configuration of 2-cyclohexenone-cis-diol type phenol metabolites: redefining the role of substituents and 2-cyclohexenone conformation in electronic circular dichroism spectra

The absolute configurations of 2-cyclohexenone cis-diol metabolites resulting from the biotransformation of the corresponding phenols have been determined by comparison of their experimental and calculated circular dichroism spectra (TDDFT at the PCM/B2LYP/Aug-cc-pVTZ level), optical rotations (calculated at the PCM/B3LYP/Aug-cc-pVTZ level) and by stereochemical correlation. It is found that circular dichroism spectra and optical rotations of 2-cyclohexenone derivatives are strongly dependent on the ring conformation (M or P sofa S(5) or half-chair), enone non-planarity and the nature and positions of the hydroxy and alkyl substituents. The effect of non-planarity of the enone chromophore, including the distortion of the C=C bond, is determined for the model structures by TDDFT calculations at the PCM/B2LYP/6-311++G(2d,2p) level. Non-planarity of the C=C bond in the enone chromophore is commonly encountered in 2-cyclohexenone derivatives and it is a source of significant rotatory strength contribution to the electronic circular dichroism spectra. It is shown that the two lowest-energy transitions in acrolein and 2-cyclohexenone and its derivatives are n(C=O)-π(C=O)* and π(C=C)-π(C=O)*, as expected, while the shorter-wavelength (below 200 nm) transitions are of more complex nature. In 2-cyclohexenone and its alkyl derivatives it is predominantly a mixture of π(C=C)-π(C=C)* and π(C=C)-σ* transitions, whereas the presence of hydroxy substituent results in a dominant contribution due to the n(OH)-π(C=O)* transition. A generalized model for correlation of the CD spectra of 2-cyclohexenones with their structures is presented.     

Theoretical Studies on the Intermolecular Interactions of Aza-calix[2]arene[2]-triazines with RDX

Six fully optimized structures of the aza-calix[2]arene[2]-triazines/RDX supramo-lecular complexes have been obtained at the DFT-B3LYP/6-311++G** level,and the corresponding intermolecular interactions have been investigated using the B3LYP,mPWPW91 and MP2 methods at the 6-311++G** level,respectively.The natural bond orbital(NBO) and atoms in molecules(AIM) analyses have been performed to reveal the origin of interactions.To our interest,the result indicates that the strongest interaction is up to -22.34 kJ/mol after basis set superposition error(BSSE) and zero point energy(ZPE) correction at the MP2/6-311++G** level.Furthermore,the intermolecular interactions between aza-calix[2]arene[2]-triazines with the substituted amidos and RDX are stronger than those of other complexes.Thus,the complexes with amidos can be used as the candidates to increase the stability of explosive and eliminate the explosive wastewater.     

Structure and vibrational assignment of the enol form of 1,1,1-trifluoro-2,4-pentanedione

Molecular structure of 1,1,1-trifluoro-pentane-2,4-dione, known as trifluoro-acetylacetone (TFAA), has been investigated by means of Density Functional Theory (DFT) calculations and the results were compared with those of acetylacetone (AA) and hexafluoro-acetylacetone (HFAA). The harmonic vibrational frequencies of both stable cis-enol forms were calculated at B3LYP level of theory using 6-31G** and 6-311++G** basis sets. We also calculated the anharmonic frequencies at B3LYP/6-31G** level of theory for both stable cis-enol isomers. The calculated frequencies, Raman and IR intensities, and depolarization ratios were compared with the experimental results. The energy difference between the two stable cis-enol forms, calculated at B3LYP/6-311++G**, is only 5.89 kJ/mol. The observed vibrational frequencies and Raman and IR intensities are in excellent agreement with the corresponding values calculated for the most stable conformation, 2TFAA. According to the theoretical calculations, the hydrogen bond strength for the most stable conformer is 57 kJ/mol, about 9.5kJ/mol less than that of AA and about 14.5 kJ/mol more than that of HFAA. These hydrogen bond strengths are consistent with the frequency shifts for OH/OD stretching and OH/OD out-of-plane bending modes upon substitution of CH(3) groups with CF(3) groups. By comparing the vibrational spectra of both theoretical and experimental data, it was concluded that 2TFAA is the dominant isomer.     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

(CH3)2S与HOCl分子间的卤键和氢键相互作用

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…ClOH卤键复合物和(CH3)2S…HOCl氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOCl相比, 在两种复合物中, 10Cl—11O和12H—11O键伸缩振动频率发生显著的红移. 经MP2/6-311++G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为-11.69和-24.16 kJ·mol-1. 自然键轨道理论(NBO)分析表明, 在(CH3)2S…ClOH卤键复合物中, 引起10Cl—11O键变长的因素包括两种电荷转移: (i) 孤对电子LP(1S)1→σ*(10Cl—11O); (ii) 孤对电子LP(1S)2→σ*(10Cl—11O), 其中孤对电子LP(1S)2→σ*(10Cl—11O)转移占主要作用, 总的结果是使σ*(10Cl—11O)的自然布居数增加0.14035e, 同时11O原子的再杂化使其与10Cl成键时s成分增加, 即具有与电荷转移作用同样的“拉长效应”; 在(CH3)2S…HOCl氢键复合物中也存在类似的电荷转移, 但是11O原子的再杂化不同于前者. 自然键共振理论(NRT)进行键序分析表明, 在卤键复合物和氢键复合物中, 10Cl—11O和12H—11O键的键序都减小. 通过分子中原子理论(AIM)分析了复合物中卤键和氢键的电子密度拓扑性质.     

Cl+F2→ClF+F和F+ClF′→ClF+F′反应机理的密度泛函理论研究

用密度泛函理论(DFT)B3LYP方法, 在6-311G基组下,计算研究了反应Cl+F2→ClF+F和对称反应F+ClF′→ClF+F′的机理。求得前者的过渡态为三角形,活化能为15.57 kJ*mol-1;后者的过渡态为线形和三角形,活化能分别为11.52和196.25 kJ*mol-1。结果均经过振动分析和IRC计算验证。