Ab initio MO study on the global minimum structure of C3H3 anion: propynl-1-yl versus allenyl anions

The geometrical structures of the C₃H⁻₃ anion are surveyed at the coupled-cluster doubles (CCD) level of theory with the aug-cc-pVDZ basis set. To clarify the CCD geometries, the stable two isomers -- propynl-l-yl 1 and allenyl 2 anions -- are further optimized at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory both with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. The final energies are calculated at the CCSD(T) and the complete active space self-consistent field (CASSCF) multi-reference internally contracted CI (MRCI) levels of theory with the aug-cc-pVTZ basis set. At the MRCI level of theory including both the corrections due to the cluster energies (MRCI+Q) and the zero-point vibrational energies, the allenyl anion 2 is about 1.3 kcal mol⁻¹ lower in energy than the propynl-l-yl anion 1. These results contrast with the previous theoretical estimates, where the propynl-l-yl anion 1 is 2-3 kcal mol⁻¹ lower in energy than the allenyl anion 2. The activation energies of the intramolecular hydrogen transfer in the 1 → 2 conversion reactions are 63.5 kcal mol⁻¹ at the MRCI+Q level of theory with the aug-cc-pVTZ basis set including the zero-point energy corrections. The adiabatic electron affinity of the planer propargyl (H₂CCCH) radical, which is the global minimum of the C₃H₃ radical, is calculated to be 0.976 eV (after correction for the zero-point energy changes) at the CCSD(T) level of theory with the aug-cc-pVTZ basis set. The present electron affinity is in fairly good agreement with the experimental one (0.893 eV) observed by Oakes and Ellison.     

The use of MIDI basis set at the correlated level

The H₂O + H₂O, the NH₃ + NH₃, the BH₃ + H₂O and the Ne + Ne systems have been studied in the supermolecule approach, using several medium sized basis sets (especially the so-called MIDI basis set). The calculations have been carried out by the use of localized molecular orbitals (LMOs).

The dispersion interaction energies have been computed by a new method (Kozmutza and Kapuy; Int. J. Quantum Chem., 38 (1990) 665), whose essence lies in the use of LMO contributions at the correlated level.

The method proposed seems to be useful for the H₂O + H₂O, the NH₃ + NH₃, and the Ne + Ne systems at different intermolecular distances, using the MIDI basis, but fails in describing correctly the correlation energy for the BH₃ + H₂O system.     

Ab initio computation of the potential energy surfaces of the water·hydrocarbon complexes H2O·C2H2, H2O·C2H4 and H2O·CH4: minimum energy structures, vibrational frequencies and hydrogen bond energies

On the basis of ab initio MP2/6–31 + + G(2d,2p) calculations, we examined the potential energy surfaces of the water·hydrocarbon complexes H₂O·CH₄, H₂O·C₂H₂ and H₂O·C₂H₂ to locate all the minimum energy structures and estimate the hydrogen bond energies and vibrational frequencies associated with the C(spⁿ)---H·O and the O---H·C(spⁿ) bonds (n = 1−3). Our calculations show that H₂O·C₂H₂, H₂O·C₂H₄ and H₂O·CH₄ have two minimum energy structures (i.e., the C---H·O and O---H·C hydrogen bond forms), but H₂O·C₂H₄ has only one when the vibrational motion is taken into account, the O---H·C hydrogen bond form. We have also computed the barrier for the interconversion from one minimum to the other. The fully optimized geometries of H₂O·CH₄, H₂O·C₂H₄ and H₂O·C₂H₂ as well as the vibrational shifts of the C---H stretching frequencies in their C---H·O hydrogen-bonded forms are in good agreement with the available experimental data. The calculated hydrogen bond energies show that the C(spⁿ---H·O bond strengths decrease in the order C(sp)---H·O>C(sp²)---H·O>C(sp³)---O>C(sp³---H·O, which is also consistent with the available experimental data.     

加拿大油砂沥青减压渣油在CO/H2-H2O作用下的热改质特性研究

以加拿大油砂沥青大于420℃的减压渣油（BVR）为原料,对比研究其在CO/H₂-H₂O和N₂体系中的热改质特性,通过系统分析BVR在H₂-H₂O、CO-H₂O、N₂-H₂O等不同氢源下的热改质特性以揭示CO/H₂-H₂O对渣油热改质的作用机制,最后探讨合成气压力、含水量以及温度对BVR临CO/H₂-H₂O改质生焦倾向的影响。结果表明,与临氮改质相比,相同反应条件下,合成气和水可使BVR热改质的生焦诱导期延长3.5-6.5 min;相同生焦率（约0.1%）时,合成气和水可显著提升BVR热改质降黏率,410℃时相对临氮改质的降黏率为29.1%,而420℃时可达54.6%。比较不同氢源下BVR热改质的生焦诱导期、改质油黏度和安定性、渣油转化率发现,H₂-H₂O、CO-H₂O、N₂-H₂O等均对BVR热改质表现出与CO/H₂-H₂O相同的促进效果,各氢源作用活性的大小顺序为H₂-H₂O > CO/H₂-H₂O > CO-H₂O > N₂-H₂O。由此可知,CO/H₂-H₂O对渣油热改质的促进作用可归因于氢气、CO水热变换新生氢和水热裂解的综合效应,且其中氢气的作用仍最显著。合成气压力、含水量和反应温度可通过影响不同氢源的贡献而调控BVR临CO/H₂-H₂O改质生焦倾向。低成本易获取的合成气可以提供BVR热降黏改质所需氢源,水能够通过CO水热变换反应供出新生活泼氢而协同合成气实现BVR高效改质。     

Long-range dispersion coefficients for like centrosymmetric linear molecules and an application to H2-H2

The real spherical tensor theory of long-range intermolecular coefficients developed in previous papers is applied to derive explicit formulae for the first three dispersion coefficients for like centrosymmetric linear molecules. The expansion of angle-dependent coefficients in associated Legendre polynomials allows one to identify the isotropic and anisotropic components of the dispersion interaction in terms of London dispersion constants, the treatment of higher coefficients being simplified by the coupling of the elementary (l, l′)-polarizations to resultant angular momenta L_A and L_B onto each molecule. The contributions from all coupling schemes are given explicitly for C₆, C₈, C₁₀, and numerical results are presented for H₂-H₂ using two-term reduced spectra values from the Kaiserlautern group.     

四元体系H3BO3-MgSO4-MgCl2-H2O在273.15K的溶解度等温线

作者考察了H₃BO₃-MgSO₄-MgCl₂-H₂O四元体系在273.15±0.05K温度下的溶解度等温线。四元体系共饱和点的液相组成为33.44%(重量)MgCl₃,1.77%MgSO₄和0.65%H₃B0₃。     

Ab initio study of the two competitive channels HO2 + H → H2 + O2 and HO2 + H → H2O + O

Saddle point geometries and barrier heights have been calculated for the H abstraction reaction HO₂(²A″)+H(²S) → H₂(¹Σ⁺_g)+O₂(³Σ⁻_g) and the concerted H approach-O removing reaction HO₂ (²A″)+H(²S) → H₂O(¹A₁)+O(³P) by using SDCI wavefunctions with a valence double-zeta plus polarization basis set. The saddle points are found to be of C_s symmetry and the barrier heights are respectively 5.3 and 19.8 kcal by including size consistent correction. Moreoever kinetic parameters have been evaluated within the framework of the TST theory. So activation energies and the rate constants are estimated to be respectively 2.3 kcal and 0.4×10⁹ ℓ mol⁻¹ s⁻¹ for the first reaction, 20.0 kcal and 5.4.10⁻⁵ ℓ mol⁻¹ s⁻¹ for the second. Comparison of these results with experimental determinations shows that hydrogen abstraction on HO₂ is an efficient mechanism for the formation of H₂ + O₂, while the concerted mechanism envisaged for the formation of H₂O + O is highly unlikely.     

Distance dependence analysis of the electron transfer reactivity for the metastable FeOH2/FeOH2 system: an ab initio investigation

A novel theoretical scheme and ab initio application in discussing the electron transfer (ET) reactivity are presented in this paper, and are also calibrated in terms of the mono-hydrated iron ion system, Fe²⁺–OH₂/Fe³⁺–OH₂. The detailed geometry optimizations have been made at UMP2(full)/6-311+G^* level, and the activation geometrical configuration and the energy have been obtained at this level of theory using the activation model and the ab initio potential energy surface fitted from MP2(full)/6-311+G^* single point energies. The corresponding energy quantities (such as the activation energy, and dissociation energy) have also been obtained at different levels of theory (HF, MP2, MP3, MP4, QCISD and PUHF, PMP2 and PMP3 with the spin-projection) and a same basis set (6-311+G^*). The electron correlation calculations include the all electron correlation and the valence electron correlation. The electronic transmission coefficient is calculated using the ab initio potential energy surface slopes and the coupling matrix element determined from the two-state model and the Slater-type d-electron wave functions. The pair distribution function is calculated using two different schemes. Taking the pair distribution function and the local ET rate into account, a statistically averaged overall observed ET rate scheme and a spherically averaged local ET rate scheme are proposed. The relevant kinetic parameters are obtained in terms of these new schemes at different ab initio calculational levels. The contact distance dependence of these parameters and the applicability of the presented models and ab initio calculational method are also discussed.     

高能体系分子间相互作用研究: 含NNO~2和NH~2混合物

以abinitioHF/6-31G^*计算求得NH~3+NH~2NO~2的两种优化构型,经MP4电子相关能校正和Boys-Bernardi方案校正基组叠加误差求得精确的分子间相互作用能。还用PM3方法计算研究TATB(均三氨基三硝基苯)分别与HMX(奥克托金)和RDX(黑索金)的混合体系,经色散能校正电子相关近似地求得分子间相互作用能。结果表明,NH~3与NH~2NO~2之间的最大结合能为-38.32kJ/mol;分子间相互作用增强了N-NO~2键强度;TATB与HMX,RDX的结合能远大于石墨与HMX或RDX的结合能,表明TATB对HMX和RDX的润湿和钝感作用较石墨更强。     

Theoretical investigations of structure and enzymatic mechanisms of aspartyl proteinases : Part I. Ab-initio calculations on an active site model: hydrogen diformiate with H2O and H3O

The active site of aspartyl proteinases (Asp) was modelled as two formiates connected with a proton and set in geometry corresponding to Asp 32 and Asp 215 side chain carboxylate groups of endothiapepsin. The shared solvent molecule was alternatively H₂O and H₃O⁺. Their positions and those of hydrogen-bonded protons were optimized using the STO-3G basis set. Full geometry optimizations were made of the hydrogen diformiate complexes with H₂O and H₃O⁺. Asymmetric hydrogen-bonded structures resulted from these calculations, except for the fully optimized complex with H₂O. In the complexes with H₃O⁺, one proton moved consistently to the proximate carboxylic oxygen yielding a neutral, hydrated formic acid dimer. Interaction energies and proton potential energy curves were calculated using the 4-31G basis set. The interaction energy with H₂O was found to be 20.49 kcal mol⁻¹ and 202.75 kcal mol⁻¹ with H₃O⁺.     

Valence and core ionization energies for SiH4, H3SiCl and H3CCl from SCF and PNO-CEPA calculations

Valence and core ionization energies for the molecules SiH₄, H₃SiCl and H₃CCl have been obtained from RHF-SCF and PNO-CEPA wavefunctions. The calculated vertical ionization energies IE_n^v(n = 1, 2, 3) agree within about 0.2 eV with the PE spectroscopically determined values for SiH₄, H₃SiCl and H₃CCl. The increase in IE₁^v for H₃SiCl relative to H₃CCl, which has been previously ascribed to a different d-orbital participation in the two molecules, is already reproduced by calculations with (s, p) basis sets. The ΔSCF method has been used to predict several core ionization energies. In addition, relative intensities calculated for the ionizations into the ²A₁ states of SiH₄⁶ are compared with available experimental data, and some spectroscopic constants for SiH₄ and its ionic states are predicted.     

3-硝基-1,2,4-三唑-5-酮与NH3及H2O分子间相互作用的理论研究

在DFT-B3LYP/6-311＋＋G**水平上, 求得3-硝基-1,2,4-三唑-5-酮(NTO)/NH₃和NTO/H₂O两种超分子体系势能面上5种全优化构型. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得NTO与NH₃和H₂O的分子间最大相互作用能依次为－37.58和－30.14 kJ/mol, 表明NTO与NH₃的分子间相互作用强于与H₂O的作用. 超分子体系中电子均由NH₃或H₂O向NTO转移, 相互作用能主要由强氢键所贡献, 由自然键轨道分析揭示了相互作用的本质. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围从单体形成超分子的热力学性质变化. 发现由NTO和NH₃形成超分子II和III在常温下可自发进行; 而NTO和H₂O只在低温下才能自发形成IV, V和VI超分子.     

Basis set effects on the stability of the Cl2O3 isomers using B3P86 and B3LYP methods of density functional theory

Four Cl₂O₃ isomers have been calculated using density functional theory with B3P86 and B3LYP functionals and various basis sets. The energy hypersurfaces of Cl₂O₃ are very flat and the relative energies of the isomers which have hypervalent characters such as ClOCl(O)O and ClClO₃ are strongly dependent on the basis sets. The stability for Cl₂O₃ isomers are in the order of (1)ClOOOCl(C_s), (2)ClOOOCl(C₂), (3)ClClO₃ and (4)ClOCl(O)O with ClOCl(O)O being most stable. We suggest that at least the cc-pV6Z(-ghi) basis set for Cl and the cc-pVTZ basis set for O are required to obtain reliable relative stabilities of Cl₂O₃ isomers with hypervalent characters.     

AB INITIO STUDIES ON THE INTERACTION BETWEEN CH3OH AND H2CO

The interaction between CH₃OH and H₂CO has been studied by abinitio method (6-311G basis set).It has been found that there are twopossible complexes:(Ⅰ)H-bond Br?nsted complex CH₃OH...OCH₂(Ⅱ)Donorscceptor Lewis complex CH₃OH...CH₂P.HF/6-311G stabilization energiesfor complex(Ⅰ)and complex(Ⅱ)are 25.1 and 17.1 kJ/mol respectively.     

Ab initio studies on photochemical reactions of Al atoms with H2 molecules

Potential energy surfaces are computed for the five lowest electronic states of the Al + H₂ system in its symmetric nuclear arrangement. Mechanisms of photochemical reactions of Al atoms with H₂ molecules are proposed, based on the calculated potential energy surfaces. The insertion reaction of the ground-state Al atom into the H₂ molecule is difficult under normal conditions. However, photoexcited Al atoms are capable of reacting with H₂ molecules along different pathways. The results obtained are consistent with experimental findings. The potential energy profiles of the dissociation reaction, AlH₂ → AlH + H, are traced by employing the UMP2 energy gradient method. Photocexcited Al atoms react with H₂ molecules along the 2 ²A₁ state pathway, and the AlH₂(²Σ_g⁺) formed dissociates easily into AlH(¹Σ) and H(²S). The dissociation reaction of ground-state AlH₂ is difficult.     

The interaction of H2, CO, CO2, H2O and NH3 on ZnO surfaces: an Oniom Study

We have used the Oniom method with three layers in order to study the interaction of CO, H₂, H₂O, NH₃ and CO₂ molecules with the ZnO surfaces using a (ZnO)₃₄₈ cluster model. The layers are divided into the high layer at the CCSD level, the medium layer at the RHF level and the low level layer using the UFF force field method. The orbital and binding energies of the adsorbed molecules, Mulliken and ChelpG charges as well as geometrical parameters were analyzed and compared with the available experimental and theoretical data.     

Ab initio study of H2SO4 rotamers

This work uses ab initio calculations to obtain harmonic frequencies and anharmonic constants for the O–H symmetric and asymmetric stretches of H₂SO₄ in its C₂, C_s, C_1a, and C_1b configurations. In addition, a high-resolution potential energy surface is calculated as a function of both O=S–O–H dihedral angles in order to accurately obtain minimum and saddle point energies. The resulting peak positions and Boltzmann populations are compared to experimental frequencies and intensities and provide evidence for the assignment of rotamers in H₂SO₄ as suggested in recent work.     

Ab initio prediction of MnOH2/MnOH2 electron transfer reactivity at electron correlation level

An ab initio calculation of electron transfer reactivity of Mn²⁺OH₂/Mn³⁺OH₂ redox system, in gas and in solution, has been made in this paper. The detailed geometry optimization and the scanning of the potential energy surfaces have been carried out at UMP2 (full)/6-311+G^* level. The relevant energy quantities (such as the activation energy and the binding energy) are calculated at different levels of theory (HF, MP2, MP3, MP4 and QCISD and corresponding spin-projection PUHF, PMP2 and PMP3) with the same basis set (6-311+G^*). Both all-electrons and valence electrons have been correlated in energy calculations. The electronic transmission coefficient is calculated using the ab initio potential energy surface slopes and the coupling matrix element determined from the two-state model and the Slater-type d-electron wave functions. The pair distribution function is obtained using both a classical sphere approximation scheme and a novel scheme. The relevant kinetic parameters are obtained at different ab initio levels in terms of the new model. The contact-distance dependence of these parameters and the applicability of the presented models are also discussed.     

A theoretical study of the reaction of SiH2 with C2H2 and C2D2

Potential energy surfaces of the reaction of SiH₂ and C₂H₂ (and C₂D₂) have been calculated by means of ab initio molecular orbital theory at the QCISD/6-311G++(2df, 2p)//MP2/6-31G(d, p) level with corrections for the triple excitations to the QCISD energies. The barrier heights for the two reaction channels of the adduct, thus calculated, were further utilized for the dynamical calculation of the rate constants in the framework of quantum statistical Rice-Ramsperger-Kassel theory. Contributions of the rate constants of the various pathways to the total rate constant (K_T) for the disappearance of the reactants are critically examined and compared with experiment. The pressure dependence of K_T(C₂H₂) is primarily due to the formation of silirene. K_T(C₂D₂) is consistently higher than K_T(C₂H₂). The standard heat of formation of silirene is predicted to be 72.1 ± 3 kcal/mol. Rearrangement of silirene to vinylsilylene requires an activation energy smaller than that to silylacetylene.     

叠氮化氢二聚体的分子间相互作用

通过ａｂ ｉｎｉｔｉｏＨＦ／６－３１Ｇ＊计算求得叠氮化氢二聚体（ＨＮ３）２势能面上四种优化结构型,经ＭＰ４ＳＤＴＱ电子相关校正和基组叠加误差２校正求得这些构型下的分子间相互作用能。结果表明,在四种优化构型中Ｎ３Ｈ．．．ＮＨＮ２最稳定,其分子间相互作用能为－１６．０７ＫＪ．ｍｏｌ＾－１。