C2F51分子C—I键解离的自旋-轨道从头算研究

采用考虑相对论效应的6—311G^＊＊全电子基组与多参考微扰理论,计算了该分子的包含自旋-轨道耦合效应的垂直激发能和基态、激发态C—I键解离势能曲线．理论计算发现,势能曲线3^3A＂与1^1A＂,2^1A＇出现交叉,交叉区域在C—I键长为0．241nm附近;基态1^A＇到激发态3^3A＂（^3Q0）的垂直激发能为4．658eV,与实验值4．662eV非常吻合．讨论了C2F5I分子作为碘激光介质的可行性．     

正、负和中性TiP10团簇结构与电子性质的密度泛函研究

采用密度泛函理论的B3LYP方法研究了正、负和中性TiP10团簇的几何构型和电子结构. 计算结果表明, 中性TiP10团簇的基态构型为金属夹心结构, 正、负离子团簇同样具有该基态稳定结构. 通过对基态稳定结构的分子轨道分析表明, δ键对形成夹心结构起到重要作用. 理论计算得到的中性TiP10团簇的垂直和绝热电离能分别为7.84和7.68 eV, 垂直和绝热电子亲和势分别为3.18和3.35 eV.     

S3分子的结构与最可几解离过程理论研究

用量子化学从头算方法,计算得到了S3分子优化的几何结构(基态和几个低激发态);研究了其在C2v与Cs群对称性下的解离过程,给出两种群对称性下基态和激发态的势能曲面及等高投影图;在两种群对称性下均得到在相应势能曲面上的优化的解离路径.     

NaC分子基态及其低电子激发态

运用包含Davidson修正的多参考组态相互作用(MRCI+Q)方法结合6-311++G(3df,3pd)基组计算了NaC分子基态(X4∑)以及三个低电子激发态(a2Π, b2∑, A4Π)的势能曲线(PECs), 确定出相应态的平衡键长Re和垂直激发能Te. 然后将PECs拟合到Murrel-Sorbie(MS)解析势能函数形式, 继而获得各态的光谱数据: 谐振频率ωe、离解能De、非谐性常数ωeΧe、转动常数Be、Drot和振转耦合常数αe. 计算结果表明: X4∑、a2Π、b2∑是三个束缚电子态. 基态X4∑的平衡键长为0.2259 nm, 谐振频率为431 cm-1, 离解能为1.92 eV, 目前计算值与实验结果和其它理论值一致. a2Π和b2∑激发态的核间距、谐振频率分别为0.2447、0.2369 nm 和329、335 cm-1, Te分别为1.58 和1.75eV, De则为0.71和0.42 eV. A4Π态为排斥态, 其相对基态的垂直激发能为2.48 eV. 通过求解核运动的薛定谔方程找到了转动量子数J=0时NaC分子三个低电子态(X4∑, a2Π, b2∑)的全部振动能级和转动惯量.     

SSCH3光解离的多参考理论计算

摘要运用多参考微扰理论(Multireference MΦller-Plesset theory)计算了SSCH3垂直激发能及其S-C与S-S两种断键方式的绝热(Adiabatic)与非绝热(Diabatic)的基态和激发态势能曲线, 研究了在193nm激光作用下SSCH3的光解离过程, 理论计算值与实验值相符.     

CH3CH2S自由基H迁移异构化及裂解反应的理论研究

采用密度泛函方法(B3LYP)在6-311+G(d,p)基组水平上研究了CH3CH2S自由基H迁移异构化以及裂解反应的微观动力学机理. 在QCISD(T)/6-311++G(d,p)//B3LYP/6-311+G(d,p)+ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了在200~2000 K温度区间内的速率常数kTST和kCVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数kCVT/SCT. 研究结果表明, CH3CH2S自由基1,2-H迁移、1,3-H迁移、C—C键断裂和β-C—H键断裂反应的势垒ΔE≠分别为149.74, 144.34, 168.79和198.29 kJ/mol. 当温度低于800 K时, 主要发生1,2-H迁移反应, 高于1800 K时, 主要表现为C—C键断裂反应, 在1300—1800 K范围内, 1,3-H迁移反应是优势通道, 在计算的整个温度段内, β-C—H键断裂反应可以忽略.     

高分辨光解碎片平动能谱研究C2H5I的紫外光解

利用新研制的小型光解碎片平动能谱仪(飞行距离仅5 cm), 研究了在279.71, 281.73, 304.02和304.67 nm等四个波长下C₂H₅I的光解动力学. 在I*(²P_1/2)、I(²P_3/2)的不同光解通道, 获得了分辨出一些振动峰的平动能谱. 这些振动峰归属于光解碎片C₂H₅I的(v₂ = 0, 1, 2,…)各伞形振动激发态(v₂, 540 cm-1), 得到了这些振动激发态的布居. 并且, 获得了C₂H₅I中C—I键的解离能D₀(C—I) = 2.314 ± 0.03 eV. 还得到了C₂H₅I在上述波长下光解后的能量分配: 对于I*通道, 在281.73, 304.02 nm下, 内能占可资用能的比例分别为22.1%, 22.4%;对于I通道, 在279.71, 304.67 nm下, 分别为25.2%, 25.9%. 这些比例值更为可靠, 因将一些振动峰分辨.     

CH2SH与NO2双自由基反应机理的理论研究

在B3LYP/6-311＋＋G(2df,p)水平上优化了标题反应驻点物种的几何构型, 并在相同水平上通过频率计算和内禀反应坐标(IRC)分析对过渡态结构及连接性进行了验证. 采用双水平计算方法HL//B3LYP/6-311＋＋G(2df,p)对所有驻点及部分选择点进行了单点能校正, 构建了CH2SH＋NO2反应体系的单重态反应势能剖面. 研究结果表明, CH2SH与NO2反应体系存在4条主要反应通道, 两个自由基中的C与N首先进行单重态耦合, 形成稳定的中间体HSCH2NO2 (a). 中间体a经过C—N键断裂和H(1)—O(2)形成过程生成主要产物P1 (CH2S＋trans-HONO), 此过程需克服124.1 kJ•mol－1的能垒. 中间体a也可以经过C—N键断裂及C—O键形成转化为中间体HSCH2ONO (b), 此过程的能垒高达238.34 kJ•mol－1. b再经过一系列的重排异构转化得到产物P2 (CH2S＋cis-HONO), P3 (CH2S＋HNO2)和P4 (SCH2OH＋NO). 所有通道均为放热反应, 反应能分别为－150.37, －148.53, －114.42和－131.56 kJ•mol－1. 标题反应主通道R→a→TSa/P1→P1的表观活化能为－91.82 kJ•mol－1, 此通道在200～3000 K温度区间内表观反应速率常数三参数表达式为kCVT/SCT＝8.3×10－40T4.4 exp(12789.3/T) cm3•molecule－1•s－1.     

CH3SH与NO2的反应机理及动力学

在G3B3, CCSD(T)/6-311＋＋G(d,p)//B3LYP/6-311＋＋G(d,p)水平上详细研究了CH3SH与基态NO2的微观反应机理. 在B3LYP/6-311＋＋G(d,p)水平得到了反应势能面上所有反应物、过渡态和产物的优化构型, 通过振动频率分析和内禀反应坐标(IRC)跟踪验证了过渡态与反应物和产物的连接关系. 在CCSD(T)/6-311＋＋G(d,p)和G3B3水平计算了各物种的能量, 得到了反应势能面. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)并结合小曲率隧道效应模型(SCT), 分别计算了在200～3000 K温度范围内的速率常数kTST, kCVT和kCVT/SCT. 研究结果表明, 该反应体系共存在5个反应通道, 其中N进攻巯基上H原子生成CH3S＋HNO2的通道活化势垒较低, 为主要反应通道. 动力学数据也表明, 该通道在200～3000 K计算温度范围内占绝对优势, 拟合得到的速率常数表达式为k1CVT/SCT＝1.93×10－16T0.21exp(－558.2/T) cm3&;#8226;molecule－1&;#8226;s－1.     

N_2H体系的势能面及振动激发态的理论研究

本文采用含Davidson校正的多参考组态相互作用MRCI+Q方法以及aug-cc-pVQZ基组计算了N2H体系的基态与第一激发态近24000个从头算能量点,使用三次样条插值法构建了高精度的全域绝热势能面.基态势能面具有一个极小点以及高于极小点0.44 eV的鞍点.基于该势能面上发生的N(~4S)+NH(X~3∑~-)→H(~2S)+N_2(X~1∑_g~g)反应是一个无能垒的放热反应,放热量为6.172 eV,比过去的势能面更接近实验值(5.956 eV).基态与第一激发态势能面间存在明显的锥形交叉,说明两个态之间存在非绝热耦合.此外,采用Lanczos算法预测了基态N_2H分子的振动能级.     

HO2自由基电子激发态的理论研究

采用CASPT2/CASSCF方法对HO2自由基进行统计算, 优化了三个电子态的稳定点几何构型, 得到详细的频率数据. 利用垂直激发计算确定了3个里德堡态、11个价电子态的电子结构以及在三种理论水平上(CASSCF, SS-CASPT2和MS-CASPT2)的能量信息. 计算中使用了ANO-L和ANO-L+基组, 验证了已知实验数据的同时, 通过与其它理论计算结果的对比, 揭示了应用弥散轨道系数对于该体系激发态研究的重要性.     

Electronic spectrum and photodissociation of ClONO in comparison to BrONO

A theoretical study of the low-lying singlet and triplet states of ClONO is presented. Calculations of excitation energies and oscillator strengths are reported using multireference configuration interaction, MRD-CI, methods with the cc-pVDZ + sp basis set. The calculations predict the dominant transition, 4(1)A' <-- 1(1)A', at 5.70 eV. The transition 2(1)A' <-- 1(1)A', at 4.44 eV, with much lower intensity nicely matches the experimental absorption maximum observed around 290 nm (4.27 eV). The potential energy curves for both states are found to be highly repulsive along the Cl-O coordinate implying that direct and fast dissociation to the Cl + NO2 products will occur. Photodissociation along the N-O coordinate is less likely because of barriers on the order of 0.3 eV for low-lying excited states. A comparison between the calculated electronic energies related to the two dominant excited states of ClONO and BrONO indicates that the transitions lie about 0.6 eV higher if bromine is replaced by chlorine. The stratospheric chemistry implications of ClONO and BrONO are discussed.     

Theoretical studies on low-lying electronic states of cyanocarbene HCCN and its ionic states

Geometries of 10, 7, and 6 low-lying states of the HCCN neutral radical, its anion and cation, were optimized by using the complete active space self-consistent field (CASSCF) method in conjunction with the aug-cc-pVTZ basis set, respectively. Taking the further correlation effects into account, the second-order perturbations (CASPT2) were carried out for the energetic correction. Vertical excitation energies (T(v)) at the ground state geometry of the HCCN neutral radical were calculated for 11 states. The results of our calculations suggest that the spin-allowed transitions of HCCN at 4.179, 4.395, 4.579, 4.727 and 5.506 eV can be attributed to X(3)A' --> 2(3)A', X(3)A' --> (3)A', X(3)A --> 3(3)A', X(3)A' --> 2(3)A', and X(3)A' --> 3(3)A', respectively. The singlet-triplet splitting gap of HCCN is calculated to be 0.738 eV. The vertical and adiabatic ionization energies were obtained to compare with the PES data. The results we obtained were consistent with the available experiment results.     

On the electronically excited states of uracil

Vertical excitation energies in uracil in the gas phase and in water solution are investigated by the equation-of-motion coupled-cluster and multireference configuration interaction methods. Basis set effects are found to be important for converged results. The analysis of electronic wave functions reveals that the lowest singlet states are predominantly of a singly excited character and are therefore well described by single-reference equation-of-motion methods augmented by a perturbative triples correction to account for dynamical correlation.Our best estimates for the vertical excitation energies for the lowest singlet n --> pi* and pi --> pi* are 5.0 +/- 0.1 eV and 5.3 +/- 0.1 eV, respectively. The solvent effects for these states are estimated to be +0.5 eV and +/- 0.1 eV, respectively. We attribute the difference between the computed vertical excitations and the maximum of the experimental absorption to strong vibronic interaction between the lowest A" and A' states leading to intensity borrowing by the forbidden transition.     

Quantum dynamical study of the O(1D)+HCl reaction employing three electronic state potential energy surfaces

Quantum dynamical calculations are reported for the title reaction, for both product arrangement channels and using potential energy surfaces corresponding to the three electronic states, 1 1A', 2 1A', and 1 1A", which correlate with both reactants and products. The calculations have been performed for J=0 using the time-dependent real wavepacket approach by Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. Reaction probabilities for both product arrangement channels on all three potential energy surfaces are presented for total energies between 0.1 and 1.1 eV. Product vibrational state distributions at two total energies, 0.522 and 0.722 eV, are also presented for both channels and all three electronic states. Product rotational quantum state distributions are presented for both product arrangement channels and all three electronic states for the first six product vibrational states.     

The electronic spectrum of chloroformic acid in comparison to formic acid

The electronic spectra of chloroformic acid ClCOOH and formic acid HCOOH are computed in large-scale multireference configuration interaction (MRD-CI) calculations. The computed spectrum of formic acid is in reasonable agreement with prior calculations and experimental data. The first electronic transition of ClCOOH is computed at 6.41 eV (193.4 nm), about 0.5 eV higher than in HCOOH. Together with five strong transitions calculated at 7.66 eV (161.9 nm; 2(1)A' <-- X(1)A'), 8.36 eV (148.3 nm; 3(1)A' <-- X(1)A'), 8.49 eV (146.0 nm; 4(1)A' <-- X(1)A'), 9.00 eV (137.8 nm; 5(1)A' <-- X(1)A'), and 9.44 eV (131.3 nm; 7(1)A' <-- X(1)A'), this can serve as a guideline for experimental search of ClCOOH.     

Bromine-loss and hydrogen-loss dissociations in low-lying electronic states of the CH3Br+ ion studied using multiconfiguration second-order perturbation theory

Complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with an ANO-RCC basis were performed for the 1(2)A', 1(2)A", 2(2)A', and 2(2)A" states of the CH3Br+ ion. The 1(2)A' state is predicted to be the ground state. The 2(2)A' state is predicted to be a bound state. The adiabatic and vertical excitation energies and the relative energies at the molecular geometry were calculated, and the energetic results for 2(2)A' and 2(2)A" are in reasonable agreement with the experimental data. Potential energy curves (PECs) for Br-loss and H-loss dissociations from the four C(s) states were calculated at the CASPT2//CASSCF level and the electronic states of the CH3(+) and CH2Br(+) ions as the dissociation products were determined by checking the relative energies and geometries of the asymptote products along the PECs. In the Br-loss dissociation, the 1(2)A', 1(2)A", and 2(2)A' states correlate with CH3(+) (X1A1') and the 2(2)A" state correlates with CH3(+) (1(3)A"). The energy increases monotonically with the R(C-Br) value along the four Br-loss PECs. In the H-loss dissociation the 1(2)A', 1(2)A", 2(2)A', and 2(2)A" states correlate with the X(1)A(1), 1(3)A", 1(3)A', and 1(1)A" states (1(3)A' lying above 1(1)A") of CH2Br(+), respectively. Along the 2(2)A" H-loss PEC there is an energy barrier and the CASSCF wave functions at large R(C-H) values have shake-up ionization character. Along the 2(2)A' H-loss PEC there are an energy barrier and a minimum. At the end of the present paper we present a comprehensive review on the electronic states and the X-loss and H-loss dissociations of the CH(3)X(+) (X = F, Cl, and Br) ions on the basis of our previous studies and the present study.     

Photodissociation dynamics of the ethoxy radical investigated by photofragment coincidence imaging

The photodissociation dynamics of the ethoxy radical (CH3CH2O) have been studied at energies from 5.17 to 5.96 eV using photofragment coincidence imaging. The upper state of the electronic transition excited at these energies is assigned to the C2A'state on the basis of electronic structure calculations. Fragment mass distributions show two photodissociation channels, OH + C2H4 and CH3 + CH2O. The presence of an additional photodissociation channel, identified as D + C2D4O, is revealed in time-of-flight distributions from the photodissociation of CD3CD2O. The product branching ratios and fragment translational energy distributions for all of the observed mass channels are nonstatistical. Moreover, the significant yield of OH + C2H4 product suggests that the mechanism for this channel involves isomerization on the excited-state surface. Photodissociation at a much lower yield is seen following excitation at 3.91 eV, corresponding to a vibronic band of the B2A' <-- X2A' transition.     

PdO0,±1, PdH0,±1及PdOH分子结构和电子性质的从头计算

利用多种从头计算方法对PdO0,±1, PdH0,±1以及PdOH的结构和电子性质进行了理论研究. 优化了PdO分子两个能量相近的态(3∏和3∑-), 结果表明PdO 的3∏态较3∑-态更为稳定, 在单双激发耦合簇理论(CCSD)水平下3∑-态的总能量比3∏态的高出0.286 eV. 在结构优化的基础上计算了PdO以及PdH的绝热电离能(AIE)和电子亲和能(AEA), 计算结果与实验值符合得很好. PdOH的基态为Cs对称性的角型结构Pd—OH(2A'态), 另外还优化得到两个2A'态的亚稳结构, 分别对应于Pd—OH和O—Pd—H, CCSD水平下两个亚稳态的总能量较基态分别高0.405和2.284 eV, 优化得到了连接这两个2A'态的过渡态, 并计算了相应的反应能垒.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.