共查询到20条相似文献,搜索用时 62 毫秒
1.
We present a sub-0.3 K accuracy, ground-state one-dimensional potential energy curve of the metastable linear configuration of the (H(2))(2) cluster calculated exclusively with explicitly correlated Gaussian functions with shifted centers. The H(2) internuclear distance is kept at the isolated H(2) vibrational ground-state average value of 1.448 736 bohr and the intermonomer separation is varied between 2 and 100 bohrs. The analytical gradient of the energy with respect to the nonlinear parameters of the Gaussians (i.e., the exponents and the coordinates of the shifts) has been employed in the variational optimization of the wave function. Procedures for enlarging the basis set and for adjusting the centers of the Gaussians to the varying intermonomer separation have been developed and used in the calculations. 相似文献
2.
应用量子化学从头计算和密度泛函理论(DFT)对HO2+C2H2反应体系的反应机理进行了研究.在B3LYP/6-311G**和CCSD(T)/6-311G**水平上计算了HO2+ C2H2反应的二重态反应势能面.计算结果表明,主要反应方式为自由基HO2的H原子和C2H2分子中的C原子结合,经过一系列异构化,最后分解得到主要产物P1 (CH2O+ HCO).此反应是放热反应,化学反应热为-321.99 kJ·mol-1.次要产物为P2 (CO2 +CH3),也是放热反应. 相似文献
3.
The six-dimensional (6D) potential energy surface (PES) for the H(2) molecule interacting with a clean Ru(0001) surface has been computed accurately for the first time. Density functional theory (DFT) and a pseudopotential based periodic plane-wave approach have been used to calculate the electronic interactions between the molecule and the surface. Two different generalized gradient approximation (GGA) exchange-correlation functionals, PW91 and RPBE, have been adopted. Based on the DFT/GGA calculated potential energies, an analytical 6D PES has been constructed using the corrugation reducing procedure. A very accurate representation of the DFT/GGA data has been achieved, with an average error in the interpolation of about 3 meV and a maximum error not larger than about 30 meV. The top site is found to be the most reactive site for both functionals used, but PW91 predicts a higher reactivity than RPBE, with lower-energy and earlier-located dissociation barriers. The energetic corrugation displayed by the RPBE PES is larger than the PW91 PES while the geometric corrugation is smaller. The differences between the two PESs increase as the distance of the molecular center of mass to the surface decreases. A direct comparison with experimental investigations on H(2)/Ru(0001) could shed light on the suitability of these XC potentials often used in DFT calculations. 相似文献
4.
An improved three-dimensional potential energy surface for the H(2)-Kr system is determined from a direct fit of new infrared spectroscopic data for H(2)-Kr and D(2)-Kr to a potential energy function form based on the exchange-Coulomb model for the intermolecular interaction energy. These fits require repetitive, highly accurate simulations of the observed spectra, and both the strength of the potential energy anisotropy and the accuracy of the new data make the "secular equation perturbation theory" method used in previous analyses of H(2)-(rare gas) spectra inadequate for the present work. To address this problem, an extended version of the "iterative secular equation" method was developed which implements direct Hellmann-Feynman theorem calculation of the partial derivatives of eigenvalues with respect to parameters of the Hamiltonian which are required for the fits. 相似文献
5.
Hinde RJ 《The Journal of chemical physics》2008,128(15):154308
We present a six-dimensional potential energy surface for the (H(2))(2) dimer based on coupled-cluster electronic structure calculations employing large atom-centered Gaussian basis sets and a small set of midbond functions at the dimer's center of mass. The surface is intended to describe accurately the bound and quasibound states of the dimers (H(2))(2), (D(2))(2), and H(2)-D(2) that correlate with H(2) or D(2) monomers in the rovibrational levels (v,j)=(0,0), (0,2), (1,0), and (1,2). We employ a close-coupled approach to compute the energies of these bound and quasibound dimer states using our potential energy surface, and compare the computed energies for infrared and Raman transitions involving these states with experimentally measured transition energies. We use four of the experimentally measured dimer transition energies to make two empirical adjustments to the ab initio potential energy surface; the adjusted surface gives computed transition energies for 56 experimentally observed transitions that agree with experiment to within 0.036 cm(-1). For 26 of the 56 transitions, the agreement between the computed and measured transition energies is within the quoted experimental uncertainty. Finally, we use our potential energy surface to predict the energies of another 34 not-yet-observed infrared and Raman transitions for the three dimers. 相似文献
6.
We report on a global potential energy hypersurface for the upper sheet of the lowest triplet state of H3+. The analytic representation is based on the double many-body expansion theory. The ab initio data points, calculated with a large cc-pV5Z basis, are represented with a root mean square deviation of only 5.54 cm(-1) in the energy region below the H(+)+2H(2S) dissociation threshold. The quasi-bound vibronic states supported by this surface have also been calculated. 相似文献
7.
A double many-body expansion potential energy surface is reported for the electronic ground state of HS(2) by fitting accurate multireference configuration interaction energies calculated with aug-cc-pVTdZ and aug-cc-pVQdZ basis sets upon separate extrapolation of the complete-active-space self-consistent field and dynamical correlation components of the total energy to the complete basis set limit. The major topographical features of the potential energy surface are examined in detail, and the model function is used for thermalized calculations of the rate constants for the S + SH → H + S(2) reaction at 298 and 400 K. A value of (1.44 ± 0.06) × 10(-11) cm(3) s(-1) is obtained at 298 K, providing perhaps the most reliable estimate of the rate constant known thus far for such a reaction. 相似文献
8.
Frank Jensen 《Theoretical chemistry accounts》1998,99(5):295-300
A total of 36 stationary points have been located on the H2CO potential energy surface by means of gradient extremal following. These 36 points are believed to represent all the important
stationary points on this surface. There is no indication that the structure of the surface becomes less complicated as the
size of the basis set is enlarged at the Hartree-Fock level of theory, but many of the second- and third-order saddle points
disappear when electron correlation is introduced. Of the ten first-order saddle points (transition structures) located, the
majority have reaction paths entering the associated minima in a side-on approach, i.e. these cannot be located by uphill
walking from the minimum.
Received: 5 February 1998 / Accepted: 21 May 1998 / Published online: 29 July 1998 相似文献
9.
An eight-dimensional time-dependent quantum dynamics wave packet approach is performed for the study of the H2+C2H-->H+C2H2 reaction system on a new modified potential energy surface (PES) [L.-P. Ju et al., Chem. Phys. Lett. 409, 249 (2005)]. This new potential energy surface is obtained by modifying Wang and Bowman's old PES [J. Chem. Phys. 101, 8646 (1994)] based on the new ab initio calculation. This new modified PES has a much lower transition state barrier height at 2.29 kcal/mol than Wang and Bowman's old PES at 4.3 kcal/mol. This study shows that the reactivity for this diatom-triatom reaction system is enhanced by vibrational excitations of H2, whereas the vibrational excitations of C2H only have a small effect on the reactivity. Furthermore, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. The comparison of the rate constant between this calculation and experimental results agrees with each other very well. This comparison indicates that the new modified PES corrects the large barrier height problem in Wang and Bowman's old PES. 相似文献
10.
Landera A Krishtal SP Kislov VV Mebel AM Kaiser RI 《The Journal of chemical physics》2008,128(21):214301
Ab initio CCSD(T)cc-pVTZ//B3LYP6-311G(**) and CCSD(T)/complete basis set (CBS) calculations of stationary points on the C(6)H(3) potential energy surface have been performed to investigate the reaction mechanism of C(2)H with diacetylene and C(4)H with acetylene. Totally, 25 different C(6)H(3) isomers and 40 transition states are located and all possible bimolecular decomposition products are also characterized. 1,2,3- and 1,2,4-tridehydrobenzene and H(2)CCCCCCH isomers are found to be the most stable thermodynamically residing 77.2, 75.1, and 75.7 kcal/mol lower in energy than C(2)H + C(4)H(2), respectively, at the CCSD(T)/CBS level of theory. The results show that the most favorable C(2)H + C(4)H(2) entrance channel is C(2)H addition to a terminal carbon of C(4)H(2) producing HCCCHCCCH, 70.2 kcal/mol below the reactants. This adduct loses a hydrogen atom from the nonterminal position to give the HCCCCCCH (triacetylene) product exothermic by 29.7 kcal/mol via an exit barrier of 5.3 kcal/mol. Based on Rice-Ramsperger-Kassel-Marcus calculations under single-collision conditions, triacetylene+H are concluded to be the only reaction products, with more than 98% of them formed directly from HCCCHCCCH. The C(2)H + C(4)H(2) reaction rate constants calculated by employing canonical variational transition state theory are found to be similar to those for the related C(2)H + C(2)H(2) reaction in the order of magnitude of 10(-10) cm(3) molecule(-1) s(-1) for T = 298-63 K, and to show a negative temperature dependence at low T. A general mechanism for the growth of polyyne chains involving C(2)H + H(C[triple bond]C)(n)H --> H(C[triple bond]C)(n+1)H + H reactions has been suggested based on a comparison of the reactions of ethynyl radical with acetylene and diacetylene. The C(4)H + C(2)H(2) reaction is also predicted to readily produce triacetylene + H via barrierless C(4)H addition to acetylene, followed by H elimination. 相似文献
11.
Tianlv Xu James Farrell Yuning Xu Roya Momen Steven R. Kirk Samantha Jenkins David J. Wales 《Journal of computational chemistry》2016,37(31):2712-2721
Using the quantum theory of atoms in molecules a near complete combined directed spanning quantum topology phase diagram (QTPD) was constructed from the nine (H2O)5 reaction‐pathways and five unique Poincaré–Hopf solutions that were found after an extensive search of the MP2 potential energy surface. Two new energy minima that were predicted from earlier work are found and include the first (H2O)5 conformer with a 3‐DQT quantum topology. The stress tensor Poincaré–Hopf relation indicated a preference for 2‐DQT (H2O)5 topologies as well as the presence of coupling between shared‐shell O? H BCPs to the hydrogen‐bond BCPs that share an H NCP. The complexity of the near complete combined QTPD was explained in terms of the O…O bonding interactions that were found in six of the nine (H2O)5 reaction‐pathways and for all points of the combined QTPD. The stabilizing role of the O…O bonding interactions from the values of the total local energy density was explored. © 2016 Wiley Periodicals, Inc. 相似文献
12.
13.
Calculations for the cumulative reaction probability N(E) (for J=0) and the thermal rate constant k(T) of the H+CH(4)-->H(2)+CH(3) reaction are presented. Accurate electronic structure calculations and a converged Shepard-interpolation approach are used to construct a potential energy surface which is specifically designed to allow the precise calculation of k(T) and N(E). Accurate quantum dynamics calculations employing flux correlation functions and multiconfigurational time-dependent Hartree wave packet propagation compute N(E) and k(T) based on this potential energy surface. The present work describes in detail the various convergence test performed to investigate the accuracy of the calculations at each step. These tests demonstrate the predictive power of the present calculations. In addition, approximate approaches for reaction rate calculations are discussed. A quite accurate approximation can be obtained from a potential energy surface which includes only interpolation points on the minimum energy path. 相似文献
14.
Hernández-Rojas J Bretón J Gomez Llorente JM Wales DJ 《The journal of physical chemistry. B》2006,110(27):13357-13362
Likely candidates for the global potential energy minima of C60(H2O)n clusters with n < or = 21 are found using basin-hopping global optimization. The potential energy surfaces are constructed using the TIP4P intermolecular potential for the water molecules, a Lennard-Jones water-fullerene potential, and a water-fullerene polarization potential, which depends on the first few nonvanishing C60 multipole polarizabilities. This combination produces a rather hydrophobic water-fullerene interaction. As a consequence, the water component of the lowest C60(H2O)n minima is quite closely related to low-lying minima of the corresponding TIP4P (H2O)n clusters. In most cases, the geometrical substructure of the water molecules in the C60(H2O)n global minimum coincides with that of the corresponding free water cluster. Exceptions occur when the interaction with C60 induces a change in geometry. This qualitative picture does not change significantly if we use the TIP3P model for the water-water interaction. Structures such as C60@(H2O)60, in which the water molecules surround the C60 fullerene, correspond to local minima with much higher potential energies. For such a structure to become the global minimum, the magnitude of the water-fullerene interaction must be increased to an unphysical value. 相似文献
15.
Fifteen structures of the (H2)2 dimer have been investigated at the MP2/[4s3p] level. The SCF and MP2 (2nd order Møller-Plesser treatment) interaction energies have been corrected for the basis set superposition error. Only the T-shaped structure has been established as a minimum on the potential energy surface. Two equivalent T-shaped structures are connected by a saddle point with a rhomboid structure.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday 相似文献
16.
Huang X Braams BJ Bowman JM Kelly RE Tennyson J Groenenboom GC van der Avoird A 《The Journal of chemical physics》2008,128(3):034312
We report a new full-dimensional potential energy surface (PES) for the water dimer, based on fitting energies at roughly 30,000 configurations obtained with the coupled-cluster single and double, and perturbative treatment of triple excitations method using an augmented, correlation consistent, polarized triple zeta basis set. A global dipole moment surface based on Moller-Plesset perturbation theory results at these configurations is also reported. The PES is used in rigorous quantum calculations of intermolecular vibrational frequencies, tunneling splittings, and rotational constants for (H2O)2 and (D2O)2, using the rigid monomer approximation. Agreement with experiment is excellent and is at the highest level reported to date. The validity of this approximation is examined by comparing tunneling barriers within that model with those from fully relaxed calculations. 相似文献
17.
Potential energy surface (PES) intersection seams of two or more electronic states from the 1 1A', 2 1A', 3 1A', 1 1A", and 2 1A" states in the C(1D)H2 reactive system are investigated using the internally contracted multireference configuration interaction method and the aug-cc-pVQZ basis set. Intersection seams with energies less than 20 kcal/mol relative to the C(1D) + H2 asymptote are searched systematically, and finally several seam lines (at the linear H-C-H, linear C-H-H, and C(2v), geometries, respectively) and a seam surface (at Cs geometries) are discovered and determined. The minimum energy crossing points on these seams are reported and the influences of the PES intersections, in particular, conical intersections, on the CH2 spectroscopy and the C(1D) + H2 reaction dynamics are discussed. In addition, geometries and energies of the 1 1A2 and 1 1B2 states of methylene biradical CH2 are reported in detail for the first time. 相似文献
18.
《Chemical physics》2005,308(3):259-266
A new ground state potential energy surface has been developed for the F+H2 reaction. Using the UCCSD(T) method, ab initio calculations were performed for 786 geometries located mainly in the exit channel of the reaction. The new data was used to correct exit channel errors that have become apparent in the potential energy surface of Stark and Werner [J. Chem. Phys. 104 (1996) 6515]. While the entrance channel and saddlepoint properties of the Stark–Werner surface are unchanged on the new potential, the exit channel behavior is more satisfactory. The exothermicity on the new surface is much closer to the experimental value. The new surface also greatly diminishes the exit channel van der Waals well that was too pronounced on the Stark–Werner surface. Several preliminary dynamical scattering calculations were carried out using the new surface for total angular momentum equal to zero for F+H2 and F+HD. It is found that gross features of the reaction dynamics are quite similar to those predicted by the Stark–Werner surface, in particular the reactive resonance for F+HD and F+H2 survive. However, the most of the exit channel van der Waals resonances disappear on the new surface. It is predicted that the differential cross-sections at low collision energy for the F+H2 reaction may be drastically modified from the predictions based on the Stark–Werner surface. 相似文献
19.
In the present paper, kinetic isotope effects of the title reaction are studied with canonical variational transition state theory on the modified Wang Bowman (MWB) potential energy surface (PES) (Chem Phys Lett 2005, 409, 249) and the ab initio calculations at the quadratic configuration interaction (QCISD (T, full))/aug‐cc‐pVTZ//QCISD (full)/cc‐pVTZ level. The calculated rate constants for the isotopic variants of this title reaction on the MWB PES have good agreement with those of the present ab initio calculations over the temperature range of 20–5000 K for the forward reactions and 800–5000 K for the reverse reactions, respectively. In particular, the forward rate constants for the title reaction and its isotopically substituted reactions have negative temperature dependences at about 40 K. Rate expressions are presented for all the studied reactions. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 289–298, 2010 相似文献
20.
Ab initio modified Gaussian-2 G2M(RCC,MP2) calculations have been performed for various isomers and transition states on the singlet C4H4 potential energy surface. The computed relative energies and molecular parameters have then been used to calculate energy-dependent rate constants for different isomerization and dissociation processes in the C4H4 system employing Rice-Ramsperger-Kassel-Marcus theory and to predict branching ratios of possible products of the C2(1Sigmag+)+C2H4, C(1D)+H2CCCH2, and C(1D)+H3CCCH reactions under single-collision conditions. The results show that C2 adds to the double C=C bond of ethylene without a barrier to form carbenecyclopropane, which then isomerizes to butatriene by a formal C2 "insertion" into the C-C bond of the C2H4 fragment. Butatriene can rearrange to the other isomers of C4H4, including allenylcarbene, methylenecyclopropene, vinylacetylene, methylpropargylene, cyclobutadiene, tetrahedrane, methylcyclopropenylidene, and bicyclobutene. The major decomposition products of the chemically activated C4H4 molecule formed in the C2(1Sigmag+)+C2H4 reaction are calculated to be acetylene+vinylidene (48.6% at Ecol = 0) and 1-buten-3-yne-2-yl radical [i-C4H3(X2A'), H2C=C=C=CH*]+H (41.3%). As the collision energy increases from 0 to 10 kcal/mol, the relative yield of i-C4H3+H grows to 52.6% and that of C2H2+CCH2 decreases to 35.5%. For the C(1D)+allene reaction, the most important products are also i-C4H3+H (55.2%) and C2H2+CCH2 (30.1%), but for C(1D)+methylacetylene, which accesses a different region of the C4H4 singlet potential energy surface, the calculated product branching ratios differ significantly: 65%-69% for i-C4H3+H, 18%-14% for C2H2+CCH2, and approximately 8% for diacetylene+H2. 相似文献