Protonated acetylene revisited

The potential energy surface for protonated acetylene has been re-examined with large basis sets and highly correlated methods. The energy difference of 3.6–3.8 kcal/mol between the classical structure and non-classical (bridged) structure computed with CCSD (T)/cc-pVQZ, CCSD(T)/6-311+G(3df,2pd), BD(T)/cc- pVQZ, BD(T)/6-311+G(3df,2pd) and CBS-APNO methods is in very good agreement with the best previous calculations, 3.7–4.0 kcal/mol. In contrast, BLYP, B3LYP, PW91, PBE and TPSS density functional methods do rather poorly, yielding −0.52. 0.29, 1.81, 2.16 and 0.62 kcal/mol, respectively, with the 6-311+G(3df,2pd) basis. MP2 calculations predict the classical structure to be a transition state; however, frequency calculations at the CCSD/6-311+G(3df,2pd) level of theory show that the classical structure is a local minimum. CCSD(T), BD(T) and CBS-APNO energy calculations along the MP2/6-311+G(3df,2pd) reaction path indicate that the classical structure is a shallow local minimum separated from the non-classical structure by a very small barrier of 0.11–0.13 kcal/mol. Because the barrier for proton exchange between the non-classical isomers via the classical structure is broad and nearly flat at the top, the tunneling splitting should be reduced, possibly accounting for the 15% difference between the calculated and experimental barrier heights. Contribution to the Fernando Bernardi Memorial Issue.     

CH3F分子间相互作用的ab initio研究

在MP2/6-311+ +G(3d,3p)电子相关校正水平上,对CH₃F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证,发现3种势能面上有极小点的构型.进一步在高级电子相关校正的MP4S-DTQ、CCSD(T)/6-311+ +G(3df,3pd)方法水平上,对其中总能量最小的构型进行精确计算,得到二聚体的结合能为-9.707kJ/mol.研究结果支持了由光谱实验结果推测的构型,解释了CH₃F二聚体的谐振频率的多样性.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

Benchmark of density functional theory methods on the prediction of bond energies and bond distances of noble-gas containing molecules

We have tested three pure density functional theory (DFT) functionals, BLYP, MPWPW91, MPWB95, and ten hybrid DFT functionals, B3LYP, B3P86, B98, MPW1B95, MPW1PW91, BMK, M05-2X, M06-2X, B2GP-PLYP, and DSD-BLYP with a series of commonly used basis sets on the performance of predicting the bond energies and bond distances of 31 small neutral noble-gas containing molecules. The reference structures were obtained using the CCSD(T)∕aug-cc-pVTZ theory and the reference energies were based on the calculation at the CCSD(T)∕CBS level. While in general the hybrid functionals performed significantly better than the pure functionals, our tests showed a range of performance by these hybrid functionals. For the bond energies, the MPW1B95∕6-311+G(2df,2pd), BMK∕aug-cc-pVTZ, B2GP-PLYP∕aug-cc-pVTZ, and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 2.0-2.3 kcal∕mol per molecule. For the bond distances, the MPW1B95∕6-311+G(2df,2pd), MPW1PW91∕6-311+G(2df,2pd), and B3P86∕6-311+G(2df,2pd), DSD-BLYP∕6-311+G(2df,2pd), and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 0.008-0.013 A? per bond. The current study showed that a careful selection of DFT functionals is very important in the study of noble-gas chemistry, and the most recommended methods are MPW1B95∕6-311+G(2df,2pd) and DSD-BLYP∕aug-cc-pVTZ.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

Theoretical study on structures and stability of SiCP2 isomers

The structures, energetics, spectral parameters and stability of the singlet SiCP₂ isomers are explored at the density functional theory and ab initio levels. Eight isomers connected by ten interconversion transition states are located at the CCSD(T)/6-311G(2d)//B3LYP/6-311G(d)level. The kinetically stable isomers and their relevant interconversion transition states are further refined at CCSD(T)/6-311+G(2df)//QCISD/6-311G(d) level. At QCISD/6-311G(d) level, one four-membered ring isomer cSiPCP and two linear structures PSiCP, SiCPP possess considerable kinetic stability (more than 15 kcal/mol). The valence bond structures of three kinetically stable SiCP₂ isomers are analyzed. The similarities and discrepancies in structure, energy and stability between SiCP₂ and its analogous C₂P₂, Si₂P₂, SiCN₂ and CSiNP molecules are also discussed. The predicted structures and spectroscopic properties are expected to be informative for the identification of the SiCP₂ in the laboratory and space.     

H-Atom abstraction from CH(3)NHNH(2) by NO(2): CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations

Stationary points of paths for H atom abstraction from CH(3)NHNH(2) (monomethylhydrazine) by NO(2) were characterized via CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations. Five transition states connecting CH(3)NHNH(2)-NO(2) complexes to a manifold that includes CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes were identified. Transition states that connect CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes to each other via H atom exchange and/or hindered internal rotation were also identified. The high point in the minimum energy path from the CH(3)NHNH(2) + NO(2) reactant asymptote to the manifold of HONO-containing product states is a transition state 8.6 kcal/mol above the reactant asymptote. From a kinetics standpoint, this value is considerably higher than the 5.9 kcal/mol value that was estimated for it based on theoretical results for H atom abstraction from NH(3) by NO(2).     

Theoretical Study on the Structures and Isomerization of Germylenoid GeH_2LiCl

Using geometrical optimization and DFT method at the B3LYP/6-311++G (3df,3pd) level, four equilibrium geometries and one transition state of GeH2LiCl were identified, and the structures at the MP2/6-311++G(3df,3pd) level were calculated simultaneously. We also studied the solvent effects on the structures of Germylenoid GeH2LiCl at the B3LYP/6-311++G (3df,3pd) level. The two more stable forms are suggested to be the p-complex and three-membered ring. The vibrational frequencies and infrared intensities were computed at the B3LYP/6-311++G (3df,3pd) level.     

Theoretical study on stability and properties of NC2O isomers

The structures, energetics, spectroscopies, and stabilities of the doublet NC(2)O radical are explored at density functional theory and ab initio levels. Nine minimum isomers are located connected by 22 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is bent NCCO 1 (0.0 kcal/mol) with (2)A' state followed by bent isomer CNCO 2 (16.7). Two isomers (1 and 2) and another high-lying species CCNO 4 (99.4) with bent structure are considerably stabilized by a barrier of at least 20 kcal/mol. All of the three isomers should be experimentally or astrophysically observable. This result is consistent with their indication of neutralization-reionization mass spectrometry experiments. Also, the calculated spectroscopic properties and bond distances of known NCCO 1 are consistent with recent experimental observations and theoretical studies. The bonding natures of the isomers 1, 2, and 4 are analyzed. Their molecular properties including the heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities are calculated at the higher levels G3//B3LYP, G3(MP2)//B3LYP, QCISD, and CCSD(T) (single-point). Possible formation strategies of the isomers 1, 2, and 4 in laboratory and space are also discussed in detail.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

On the structure and chemical bonding of Si6(2-) and Si6(2-) in NaSi6(-) upon Na+ coordination

Photoelectron spectroscopy was combined with ab initio calculations to elucidate the structure and bonding in Si6 2- and NaSi6 -. Well-resolved electronic transitions were observed in the photoelectron spectra of Si6 - and NaSi6 - at three photon energies (355, 266, and 193 nm). The spectra of NaSi6 - were observed to be similar to those of Si6 - except that the electron binding energies of the former are lower, suggesting that the Si6 motif in NaSi6 - is structurally and electronically similar to that in Si6 -. The electron affinities of Si6 and NaSi6 were measured fairly accurately to be 2.23+/-0.03 eV and 1.80+/-0.05 eV, respectively. Global minimum structure searches for Si6 2- and NaSi6 - were performed using gradient embedded genetic algorithm followed by B3LYP, MP2, and CCSDT calculations. Vertical electron detachment energies were calculated for the lowest Si6 - and NaSi6 - structures at the CCSD(T)/6-311+G(2df), ROVGF/6-311+G(2df), UOVGF/6-311+G(2d), and time-dependent B3LYP/6-311+G(2df) levels of theory. Experimental vertical detachment energies were used to verify the global minimum structure for NaSi6 -. Though the octahedral Si6 2-, analogous to the closo form of borane B6H6 2-, is the most stable form for the bare hexasilicon dianion, it is not the kernel for the NaSi6 - global minimum. The most stable isomer of NaSi6 - is based on a Si6 2- motif, which is distorted into C2v symmetry similar to the ground state structure of Si6 -. The octahedral Si6 2- coordinated by a Na+ is a low-lying isomer and was also observed experimentally. The chemical bonding in Si6 2- and NaSi6 - was understood using natural bond orbital, molecular orbital, and electron localization function analyses.     

B4O分子异构化稳定性的理论研究

采用CCSD(T)/6-311+G(2df)//B3LYP/6-311+G(d)方法, 系统研究了B₄O体系各个异构体的结构和能量, 以及重要异构体的解离和异构化稳定性. 结果表明, 单态平面三角形含-BO单元的异构体cB₃-BO ¹1能量最低, 其次是带状的异构体B₄O ¹2(10.9 kJ/mol), 并且¹1和¹2结构都具有良好的动力学稳定性. 因此对于B₄O体系, ¹1和¹2都是有可能存在的. 而文献报道的三态直线型结构BBBBO(146.0 kJ/mol)的能量比异构体¹1和¹2高得多.     

FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

Infrared spectroscopy of ozone-water complex in a neon matrix

Matrix isolation infrared spectroscopy has been applied to study an ozone-water complex of atmospheric interest. The complex was identified in the spectral region of three normal modes of ozone and water. Ab initio calculation at MP4(SDQ), QCISD, and CCSD(T) levels indicates the existence of only one stable conformer, which accords with the present experimental result. This conformer belongs to the Cs symmetry group where two molecular planes of ozone and water are perpendicular to the Cs symmetry plane. The binding energy was calculated to be 1.89 kcal/mol at the CCSD(T)/6-311++G(3df,3pd)//CCSD(T)/6-311++G(d,p) level of theory. The formation constant and atmospheric abundance of the ozone-water complex are estimated using the thermodynamic and spectroscopic data obtained.     

B2F2分子异构体结构的量子化学计算研究

采用abinitio方法对B₂F₂分子异构体的结构进行了计算研究,并与Al₂F₂分子进行了比较.结果表明,B₂F₂具有D_∞h对称性,³Σ_g电子态的直线型结构FBBF是B₂F₂分子的最稳定异构体,对文献的结果进行了修正.在UCCSD(T,full)/6－311+G(2df)水平下,F-B和B-B键长分别为0.12942和0.14820nm,振动频率出现在1860.00和1320.62cm^-1处.在UQCISD(T,full)/6－311+G(2df)//UMP2(full)/6－311+G(d)+ZPVE水平下,³Σ_g态的线性FBBF分子的垂直电离势为848.58kJ/mol,而由³Σ_g电子态的BF二聚为³Σ_g态的线性FBBF分子的焓变为59.86kJ/mol,此二聚化反应是放热反应,说明二聚化过程在能量上是有利的.     

An ab initio correlated study of the potential energy surface for the HOBr.H2O complex

The potential energy surface (PES) for the HOBr.H(2)O complex has been investigated using second- and fourth-order M?ller-Plesset perturbation theory (MP2, MP4) and coupled cluster theory with single and doubles excitations (CCSD), and a perturbative approximation of triple excitations (CCSD-T), correlated ab initio levels of theory employing basis sets of triple zeta quality with polarization and diffuse functions up to the 6-311++G(3dp,3df ) standard Pople's basis set. Six stationary points being three minima, two first-order transition state (TS) structures and one second-order TS were located on the PES. The global minimum syn and the anti equilibrium structure are virtually degenerated [DeltaE(ele-nuc) approximately 0.3 kcal mol(-1), CCSD-T/6-311++G(3df,3pd) value], with the third minima being approximately 4 kcal mol(-1) away. IRC analysis was performed to confirm the correct connectivity of the two first-order TS structures. The CCSD-T/6-311++G(3df,3pd)//MP2/6-311G(d,p) barrier for the syn<-->anti interconversion is 0.3 kcal mol(-1), indicating that a mixture of the syn and anti forms of the HOBr.H(2)O complex is likely to exist.