Computational investigation of the weakly bound dimers HOX...SO(3) (X = F, Cl, Br)

The electronic structure and thermochemical stability of the HOX-SO(3) (X = F, Cl, Br) complexes is studied using second-order M?ller-Plesset perturbation theory (MP2). The calculated dissociation energies of the HOF-SO(3), HOCl-SO(3), and HOBr-SO(3) complexes are 5.43, 6.02, and 5.98 kcal mol(-1) at MP2/6-311++G(3df,3pd) level, respectively. Anharmonic OH stretching frequencies of the HOX (X = F, Cl, Br) moieties along with the frequency shifts upon complex formation are calculated at the MP2/6-311++G(2df,2p) level. AIM and NBO analyses were also performed. Theoretical data strongly encourage performing of matrix-isolation studies of the title complexes and their spectroscopic identification.     

SiH_4与Na、Mg和Be等金属氢化物分子间反向氢键相互作用

在B3LYP/6-311++g**、MP2/6-311++g(3df,3pd)及MP2/aug-cc-pvtz水平上分别求得H3SiH…MeHn(Me=Na,Mg,Be;n=1或2)复合物势能面上的3个稳定构型,探讨了以Si-H为电子供体的红移反向氢键相互作用(IHB).经MP2/6-311++g(3df,3pd)水平的计算,在3个复合物中,含基组重叠误差(BSSE)校正的单体间相互作用能分别为-5.98、-8.65和-3.96kJ.mol-1,与MP2/aug-cc-pvtz水平下计算得到的-6.18、-9.12和-4.28kJ·mol-1接近,可见3个反向氢键复合物的相对稳定性顺序为:SiH4...MgH2SiH4...NaHSiH4...BeH2.NBO分析及对相关原子化学位移的计算表明,在复合物中,电子流向总体表现为SiH4→MeHn(n=1或2),且直接参与反向氢键形成的H3的化学位移向低场移动.与传统氢键相比,这里Si1-H3既是氢键供体,又是电子供体,从而形成反向氢键相互作用.另外,采用分子中原子理论(AIM)分别对各复合物中相关键鞍点处的电子密度拓扑性质进行了分析,结果表明3个复合物中均存在以静电性质为主的分子间反向氢键弱相互作用.     

Properties of the halogen-hydride interaction: an ab initio and "atoms in molecules" analysis

X-Cl...H-Y interactions are analyzed by applying ab initio methods as well as the Bader theory. All calculations were performed using Pople's basis sets (6-311++G(2df,2pd) and 6-311++G(3df,3pd)) as well as the Dunning-type bases (aug-cc-pVDZ and aug-cc-pVTZ) within the MP2 method. For the complexes analyzed here, X-Cl and H-Y may be treated as a Lewis acid and a Lewis base, respectively. The Cl...H interactions are rather weak or at most moderate since, for the strongest interaction of the F3...HLi complex, the binding energy calculated at the MP2/6-311++G(3df,3pd) level of approximation amounts to -3.4 kcal/mol, and the H...Cl distance is equal to 2.65 A, less than the corresponding sum of van der Waals radii. These interactions may be classified as halogen-hydride interactions. However, some of the complexes analyzed, especially F3SiCl...HBeF and F3SiCl...HBeF, are very weakly bound, probably by typical van der Waals interactions.     

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.     

The molecular structure of fluoromalononitrile as determined by gas-phase electron diffraction and ab initio calculations

The molecular structure of fluoromalononitrile was studied by means of gas-phase electron diffraction and quantum mechanical methods using HF/6-31G(d), MP2/6-311++G(2df,2pd) and DFT/B3LYP/6-31G(d), B3PW91/6-31G(d), B3LYP/6-311++G(2df,2pd) and B3PW91/6-311++G(2df,2pd). The r(g) and angle(alpha) structural parameters we obtained from the present analysis are: CC=1.487(5) A, CN=1.157(3) A, CF=1.386(5) A, CH=1.096 A (ass.), angleCCC=106.7(1.0) degrees , angleCCF=108.0(0.7) degrees , angleCCN=177.6(2.0) degrees . Uncertainties in parenthesis are 3sigma.     

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

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

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

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.     

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.     

Microwave spectrum of 3-butyne-1-thiol: evidence for intramolecular S-H...pi hydrogen bonding

The microwave spectrum of 3-butyne-1-thiol has been studied by means of Stark-modulation microwave spectroscopy and quantum-chemical calculations employing the B3LYP/6-311++G(3df,2pd), MP2/aug-cc-pVTZ, MP2/6-311++G(3df,2pd), and G3 methods. Rotational transitions attributable to two conformers of this molecule were assigned. One of these conformers possesses an antiperiplanar arrangement of the atoms S-C1-C2-C3, while the other is synclinal and stabilized by the formation of an intramolecular hydrogen bond between the H-atom of the thiol group and the pi-electrons of the C[triple bond]C triple bond. The energy difference between these conformers was estimated to be 1.7(4) kJ mol(-1) by relative intensity measurements, with the hydrogen-bonded conformer being lower in energy. The spectra of five vibrationally excited states of the synclinal conformer were observed, and an assignment of these states to particular vibrational modes was made with the aid of a density functional theory (DFT) calculation of the vibrational frequencies at the B3LYP/6-311++G(3df,2pd) level of theory.     

An Ab Initio Study of the Complexes Formed by the Positive Acetylene Ion with the Hydrogen and Nitrogen Molecules and the Argon Atom

The complexes formed by the positive acetylene ion with the hydrogen molecule, the nitrogen molecule, and the argon atom are investigated with ab initio calculations using the 6-311G** and the 6-31+G(2df,2pd) basis sets. MP2/6-311G** energies and optimum geometries are obtained, as well as single-point MP3, MP4, and QCISD(T) energies with the MP2/6-311G** optimized geometries. Single-point calculations are performed with the 6-31+G(2df,2pd) basis set at MP2/6-311G** optimized geometries.     

Theoretical study and rate constant calculation for the reactions of SH (SD) with Cl2, Br2, and BrCl

The mechanisms of the SH (SD) radicals with Cl2 (R1), Br2 (R2), and BrCl (R3) are investigated theoretically, and the rate constants are calculated using a dual-level direct dynamics method. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6-311G(d,p) and MPW1K/6-311G(d,p) levels. Higher-level energies are obtained at the approximate QCISD(T)/6-311++G(3df, 2pd) level using the MP2 geometries as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MPW1K geometries. Complexes with energies less than those of the reactants or products are located at the entrance or the exit channels of these reactions, which indicate that the reactions may proceed via an indirect mechanism. The enthalpies of formation for the species XSH/XSD (X = Cl and Br) are evaluated using hydrogenation working reactions method. By canonical variational transition-state theory (CVT), the rate constants of SH and SD radicals with Cl2, Br2, and BrCl are calculated over a wide temperature range of 200-2000 K at the a-QCISD(T)/6-311++G(3df, 2pd)//MP2/6-311G(d, p) level. Good agreement between the calculated and experimental rate constants is obtained in the measured temperature range. Our calculations show that for SH (SD) + BrCl reaction bromine abstraction (R3a or R3a') leading to the formation of BrSH (BrSD) + Cl in a barrierless process dominants the reaction with the branching ratios for channels 3a and 3a' of 99% at 298 K, which is quite different from the experimental result of k3a'/k3' = 54 +/- 10%. Negative activation energies are found at the higher level for the SH + Br2 and SH + BrCl (Br-abstraction) reactions; as a result, the rate constants show a slightly negative temperature dependence, which is consistent with the determination in the literature. The kinetic isotope effects for the three reactions are "inverse". The values of kH/kD are 0.88, 0.91, and 0.69 at room temperature, respectively, and they increase as the temperature increases.     

PH~3...H~2O体系分子间相互作用的从头计算MO研究

在MP2/6-311++G(3d,3p)水平,对PH~3...H~2O体系可能存在的氢键复合物进行了全自由度能量梯度优化,发现PH~3...H~2O体系存在三个能量极小结构A,B和C。其中结构A和B以H~2O作为质子授体,结构C以PH~3作为质子授体,结构A较结构B和C分别稳定6.52kJ/mol和8.18kJ/mol。结构A具有C~s对称性,其结构中P原子和O原子间距离为354.78nm,键角OHP为171.35ⅲ,属于接近于直线的传统型氢键结构。进一步在高级电子相关校正的MP4SDTQ下,用6-311++G(3df,3pd)基加上键函数{3s3p2d1f},通过BSSE校正,精确计算了结构A的结合能为-10.84kJ/mol。     

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.     

Formic acid catalyzed hydrolysis of SO3 in the gas phase: a barrierless mechanism for sulfuric acid production of potential atmospheric importance

Computational studies at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels are performed to explore the changes in reaction barrier height for the gas phase hydrolysis of SO(3) to form H(2)SO(4) in the presence of a single formic acid (FA) molecule. For comparison, we have also performed calculations for the reference reaction involving water assisted hydrolysis of SO(3) at the same level. Our results show that the FA assisted hydrolysis of SO(3) to form H(2)SO(4) is effectively a barrierless process. The barrier heights for the isomerization of the SO(3)···H(2)O···FA prereactive collision complex, which is the rate limiting step in the FA assisted hydrolysis, are found to be respectively 0.59 and 0.08 kcal/mol at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels. This is substantially lower than the ~7 kcal/mol barrier for the corresponding step in the hydrolysis of SO(3) by two water molecules--which is currently the accepted mechanism for atmospheric sulfuric acid production. Simple kinetic analysis of the relative rates suggests that the reduction in barrier height facilitated by FA, combined with the greater stability of the prereactive SO(3)···H(2)O···FA collision complex compared to SO(3)···H(2)O···H(2)O and the rather plentiful atmospheric abundance of FA, makes the formic acid mediated hydrolysis reaction a potentially important pathway for atmospheric sulfuric acid production.     

The change in the vibrational anharmonicity for FHRg, FArHRg, FArHN2 and FArHP2 (Rg=Ne, Ar, Kr)

The change in the vibrational anharmonic constant for the linear hydrogen-bonded complexes FHRg, FArHRg, FArHN₂, FArHP₂ (Rg=Ne, Ar, Kr) was determined at the MP2/6-311++G(3df,3pd) level of theory. It was found that in the vibrationally red-shifted complexes, hydrogen bond formation increases the anharmonic constant while in the vibrationally blue-shifted species (except FArHNe) the anharmonic constant decreases on complexation.     

[Li…X]e-[1](X=FH,OH2, NH3)的光电性质从头算

在MP2理论水平上采用6-311G基组系列计算了一价阴离子van der Waals复合物[Li…X]e-[1](X=FH, OH2, NH3)的偶极矩(μ)、平均极化率(α)以及平均一阶超极化率(β), 讨论了基组效应和电子相关效应对计算结果的影响, 比较了价电子对复合物一阶超极化率的贡献. 在MP4(SDQ)/6-311++G(2df, 2pd)水平上计算得到[Li…FH]e-[1]的μ=2.5633 a.u., α=1.0476×10³ a.u., β=1.0948×10⁵ a.u.；[Li…OH2]e-[1] 的μ=2.3204 a.u., α=1.2201×10³ a.u., β=2.1410×10⁵ a.u.；[Li…NH3]e-[1]的μ=2.4687 a.u., α=1.4817×10³ a.u., β=3.4040×10⁵ a.u.. 计算结果表明, 三种一价阴离子复合物分子均具有非常大的一阶超极化率, 而一个价电子对复合物的一阶超极化率的贡献超过1.0×10⁵ a.u..     

Ab initio studies of the conformations of ethynyl formate, cyano formate and related ethynyl and cyano compounds

The optimized geometries, relative free energies and related thermodynamic properties, harmonic frequencies, and dipole moments have been calculated at the HF and MP2 levels for ethynyl formate (1a), ethynyl acetate (1b), cyano formate, HCO₂CN (1c), cyano acetate (1d), S-ethynyl thioformate (2a), S-ethynyl thioacetate (2b), S-cyano thioformate (2c), S-cyano thioacetate (2d), N-ethynylformamide (3a), N-ethynylacetamide (3b), N-cyanoformamide (3c), and N-cyanoacetamide (3d) with the gaussian 98 program. For ethynyl formate, the calculation for 25 °C at the MP2/6-311++G(df,pd) level predicts that the Z isomer is more stable by 1.23 kcal/mol. For S-ethynyl thioformate, calculations at the MP2/6-311++G(2d,2p) level predict that the E isomer is favored by 0.71 kcal/mol at 25 °C. The E isomers of N-ethynylformamide and N-ethynylacetamide were found at all levels to be more stable than the Z isomers at 25 °C. For cyano formate and cyano acetate, calculations at the MP2/6-311++G(df,pd) level predict that the Z isomers are more stable at 25 °C by 1.50 and 2.72 kcal/mol, respectively. At this level and temperature, the Z isomers of 2c, 2d, 3c, and 3d are predicted to have free energies of 0.46, −0.07, 1.22, and 2.28 kcal/mol, respectively, relative to the E conformations. Z to E free-energy barriers at 25 °C of 8.63, 10.64, 17.63, 7.39, and 14.03 kcal/mol were calculated for 1a, 2a, 3a, 1c, and 3c at the HF/6-311G(d,p) level, and at the HF/6-311+G(d,p) level, the free-energy barrier for 2c was 7.08 kcal/mol.     

Dihydrogen trioxide clusters, (HOOOH)n (n = 2-4), and the hydrogen-bonded complexes of HOOOH with acetone and dimethyl ether: implications for the decomposition of HOOOH

Hydrogen-bonded gas-phase molecular clusters of dihydrogen trioxide (HOOOH) have been investigated using DFT (B3LYP/6-311++G(3df,3pd)) and MP2/6-311++G(3df,3pd) methods. The binding energies, vibrational frequencies, and dipole moments for the various dimer, trimer, and tetramer structures, in which HOOOH acts as a proton donor as well as an acceptor, are reported. The stronger binding interaction in the HOOOH dimer, as compared to that in the analogous cyclic structure of the HOOH dimer, indicates that dihydrogen trioxide is a stronger acid than hydrogen peroxide. A new decomposition pathway for HOOOH was explored. Decomposition occurs via an eight-membered ring transition state for the intermolecular (slightly asynchronous) transfer of two protons between the HOOOH molecules, which form a cyclic dimer, to produce water and singlet oxygen (Delta (1)O 2). This autocatalytic decomposition appears to explain a relatively fast decomposition (Delta H a(298K) = 19.9 kcal/mol, B3LYP/6-311+G(d,p)) of HOOOH in nonpolar (inert) solvents, which might even compete with the water-assisted decomposition of this simplest of polyoxides (Delta H a(298K) = 18.8 kcal/mol for (H 2O) 2-assisted decomposition) in more polar solvents. The formation of relatively strongly hydrogen-bonded complexes between HOOOH and organic oxygen bases, HOOOH-B (B = acetone and dimethyl ether), strongly retards the decomposition in these bases as solvents, most likely by preventing such a proton transfer.     

HeLi^n^+(n=0, 1)簇合物的稳定性和垂直激发态光谱

本文用ab initio研究了簇合物HeLi^n^+(n=0, 1)的几何构型和成键性质。在MP2(FULL)/6-31G**, 水平优化所得LeLi^+的平衡键长为0.2062nm, 与实验值0.205nm十分吻合。比较了HeLi^+(X^1∑^+和a^3∑^+), HeLi(X^2∑^+和a^4II)以及HLi(X^1∑^+)的稳定性, 计算了HeLi^+基态的相关能, 势能曲线和垂直激发态光谱。计算采用了6-31G**, 6-311G**,6-311G(2df, 2pd), 6-311G(3df, 2pd)和6-311+G(3df, 2pd)基组; 采用的方法包括MP2(FULL), MP4, MCSCF, MRSDCI, CCD和ST4CCD。计算表明, 同价HeLi^n^+中激发态的离解能均远比基态的大, 其中HeLi^+(a^3∑^+)的离解能最高(60.49kj/mol),说明激发态是稳定束缚态。HeLi^+基态比等电子体HLi分子基态的稳定性小得多。HeLi^+由A^1∑^+到B^1II的垂直跃迁(3σ→1π)振子强度较大而垂直跃迁能较小。