Ab initio molecular calculations with pseudopotentials: Higher quality calculations on Li2, LiH,and BeH2

New psudopotentials of higher quality have been used in ab initio molecular calculations. Simple molecular systems like Li₂, LiH, and BeH₂ have been investigated with varying basis sets. Energy, geometrical parameters, and dipole moments have been determined. The convergence properties of the various quantities have been investigated. They show characteristic behavior, e.g., stationary value for the energy where such behavior may be expected.     

Ab initio vibrational state calculations with a quartic force field: applications to H2CO, C2H4, CH3OH, CH3CCH, and C6H6

For polyatomic molecules, n-mode coupling representations of the quartic force field (nMR-QFF) are presented, which include terms up to n normal coordinate couplings in a fourth-order polynomial potential energy function. The computational scheme to evaluate third-and fourth-order derivatives by finite differentiations of the energy is fully described. The code to generate the nMR-QFF has been implemented into GAMESS program package and interfaced with the vibrational self-consistent field (VSCF) and correlation corrected VSCF (cc-VSCF) methods. As a demonstration, fundamental frequencies have been calculated by the cc-VSCF method based on 2MR-QFF for formaldehyde, ethylene, methanol, propyne, and benzene. The applications show that 2MR-QFF is a highly accurate potential energy function, with errors of 1.0-1.9% relative to the experimental value in fundamental frequencies. This approach will help quantitative evaluations of vibrational energies of a general molecule with a reasonable computational cost.     

Ab initio molecular fragment calculations with pseudopotentials: Model peptide studies

The energetics of rotation about the N? C′(ω); N? C^α(φ), and C^α? C′(?) bonds of the peptide unit have been investigated in the pseudo-FSGO fragment scheme on model compounds formamide and N-methylacetamide. The results indicated that the position of the minimum in ω is in the near vicinity of 0°, i.e., the planar arrangement of the peptide unit. The minimum in φ (C′? N? C^α? H) has been found to be 180° and in ψ(H? C^α? C′? N) to be 60°, in good agreement with PCILO and Gaussian-70 results.     

Ab initio molecular fragment calculations with pseudopotentials. Some fragments containing nitrogen and oxygen

Parameters for the OH (sp²) and NH₃ (planar, sp²) pseudopotential–FSGO molecular fragments have been established in this paper. Calculations for the molecules of formic acid, formaldehyde, and formamide show good agreement with experiment.     

Ab initio calculations ofp-Cl2C6H4 molecule and its35Cl NQR parameters

Ab initio calculations ofp-dichlorobenzene molecule were carried out using the Hartree-Fock method in the 6–31 G* valence-split basis set. The molecule was also calculated by the MNDO method in the valence sp-basis set for comparison. The populations of the valent p-orbitals of the C and CI atoms were analyzed. The optimized geometry of the molecule as well as its³⁵Cl NQR frequency and the asymmetry parameter of the electric field gradient at the³⁵CI nuclei calculated using the populations of the less diffuse components of the valent p-orbitals of the Cl atoms are in agreement with the corresponding experimental values for the -modification of 1,4-Cl₂C₆H₄.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 2177–2179, September, 1996.     

Ab initio calculations ofo- andm-Cl2C6H4 molecules and their35CI NQR parameters

The MNDO andab initio (Hartree—Fock method, split 6-31G^* basis set, full optimization of geometry) calculations for theo- andm-dichlorobenzene molecules have been performed. The populations of the valent p-orbitals of the carbon and chlorine atoms have been analyzed. Estimations of the³⁵Cl NQR frequencies and the asymmetry parameters of the electric field gradient at the³⁵CI nuclei have been made. It was confirmed that these values are not determined by the contribution of the whole valent p-shell, but only by its less diffuse component.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 2180–2183, September, 1996.     

Ab initio quantum-chemical calculations of interactions of ions and hydrides of alkali metals with C3H6 and C4H8 molecules

Calculations of the C₃H₆ · LiH, C₄H₈ · M⁺, and C₄H₈ · MH systems and of C₂H₂ · MH complexes (M = Li or Na) were carried out by the unrestricted Hartree-Fock-Roothaan (UHF) method with partial optimization of the geometry using fixed geometric parameters of the C₃H₆ and C₄H₈ molecules. The standard 3-21G and 6-31G* basis sets were used. Unlike the C₃H₆ · LiH structure, the C₄H₈ · M⁺ and C₄H₈ · MH systems are typical complexes. It was found that the C₄H₈ · M⁺, C₄H₈ · MH, and C₂H₂ · MH complexes are similar in coordination of M⁺ ions and MH molecules by carbon atoms in spite of considerable differences in the interatomic distances (–1 A) between these atoms in the C₄H₈ and C₂H₂ molecules. The heats of formation (Q), which were calculated in the UHF/6-31G* approximation and using second- and fourth-order Möller-Plesset perturbation theory taking into account the electron correlation energy in the MP2/6-31G*. MP4(SDQ)/6-31G*, and MP4(SDTQ)/6-31G* approximations, satisfy the following relationships: Q(C₂H₃ · MH) < Q(C₄H₈ · MH) < Q(C₄H₈ · M⁺). It was observed that in going from Li to Na the corresponding values of Q tend to decrease.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya. No. 7, pp. 1636–1640, July, 1996.     

Ab initio rate constants from hyperspherical quantum scattering: application to H+C2H6 and H+CH3OH

The dynamics and kinetics of the abstraction reactions of H atoms with ethane and methanol have been studied using a quantum mechanical procedure. Bonds being broken and formed are treated with explicit hyperspherical quantum dynamics. The ab initio potential energy surfaces for these reactions have been developed from a minimal number of grid points (average of 48 points) and are given by analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimized using the second order perturbation theory method with a correlation consistent polarized valence triple zeta basis set. Single point energies are calculated on the optimized geometries with the coupled cluster theory and the same basis set. The reaction of H with C2H6 is endothermic by 1.5 kcal/mol and has a vibrationally adiabatic barrier of 12 kcal/mol. The reaction of H with CH3OH presents two reactive channels: the methoxy and the hydroxymethyl channels. The former is endothermic by 0.24 kcal/mol and has a vibrationally adiabatic barrier of 13.29 kcal/mol, the latter reaction is exothermic by 7.87 kcal/mol and has a vibrationally adiabatic barrier of 8.56 kcal/mol. We report state-to-state and state-selected cross sections together with state-to-state rate constants for the title reactions. Thermal rate constants for these reactions exhibit large quantum tunneling effects when compared to conventional transition state theory results. For H+CH3OH, it is found that the CH2OH product is the dominant channel, and that the CH3O channel contributes just 2% at 500 K. For both reactions, rate constants are in good agreement with some measurements.     

Marron-weare variational computation of energy differences for motions of C2H6, B2H6, CH3CH inf2 sup+ , and CH3BH2

The passage of D_3dC₂H₆ and B2H6 toward a D₂ h bridged structure, and the motion of a methyl proton maintaining C symmetry in C₂H _inf5 ^sup+ and CH₃BH₂ are described by integral Hellmann-Feynman computations in a Frost floating spherical Gaussian basis. Marron and Weare's variational corrections to the integral Hellmann-Feynman formula forAE between statesA andB are evaluated with variational functions of the form η(ψ_A/S_AB-ψ_B)) used to refine the stateB. An analogous function ξ(ψ_B/S_AB-ψ_A) refines state A. Both η and ξ are chosen variationally to minimize Marron and Weare's functional. No obvious advantage of the variational method became apparent in this simple application.     

Ab initio calculations and molecular dynamics simulation of H2 adsorption on CN3Be3+ cluster

Hydrogen adsorption properties of the CN₃Be₃⁺ cluster have been studied using density functional theory and MP2 method with a 6–31++G** basis set. Five hydrogen molecules get adsorbed on the CN₃Be₃⁺ cluster with a hydrogen storage capacity of 10.98 wt%. Adsorption of three H₂ molecules on one of the three Be atoms in a cluster is reported for the first time. It is due to the more positive charge on this Be atom than the remaining two. The average value for H₂ adsorption energy in CN₃Be₃⁺ (5H₂) complexes is 0.41 (0.43) eV/H₂ at MP2 (wB97XD) level, which fits well within the ideal range. Adsorption energy from electronic structure calculations plays an important role in retaining the number of H₂ molecules on a cluster during atom-centered density matrix propagation (ADMP) molecular dynamics (MD) simulations. According to ADMP-MD simulations, out of five H₂ adsorbed molecules on CN₃Be₃⁺, four and two H₂ molecules remain absorbed on CN₃Be₃⁺ cluster at 275 K and 350 K, respectively, during the simulation.

     

Ab initio SCF MO calculations on the CH3Br molecule

Self-consistent fieldab initio calculations have been performed for CH₃Br. The calculated equilibrium conformation is in good agreement with experiment. Valence and core level ionization potentials, Mulliken population analysis and electronic properties are presented. The ionization potentials are in good agreement with the experimental values, except for one case in which the experimental value may be wrongly assigned. The calculated dipole moment, 2.43 Debye, is 34% or 0.6 Debye larger than the experimental value.Two of us (C.W.B. and G. del C.) wish to acknowledge the hospitality of the IBM San José Research Laboratory and joint study agreements between the IBM and the Lawrence Berkeley Laboratory and the National University of Mexico respectively.     

Ab initio valence-bond calculations of H2O

The first and second bond dissociation energies for H₂O have been calculated in anab initio manner using a multistructure valence-bond scheme. The basis set consisted of a minimal number of non-orthogonal atomic orbitals expressed in terms of gaussian-lobe functions. The valence-bond structures considered properly described the change in the molecular system as the hydrogen atoms were individually removed to infinity. The calculated equilibrium geometry for the H₂O molecule has an O-H bond length of 1.83 Bohrs and an HOH bond angle of 106.5°. With 49 valence-bond structures the energy of H₂O at this geometry was ?76.0202 Hartrees. The calculated equilibrium bond length for the OH radical was 1.86 Bohrs and the energy, using the same basis set, was ?75.3875 Hartrees. After correction for zero point energies the calculated bond dissociation energies are: H₂O → OH + H, D₁=75.38 kcal/mole and OH → O+H, D₂=54.79 kcal/mole.     

Accurate ab initio predictions of ionization energies of hydrocarbon radicals: CH2, CH3, C2H, C2H3, C2H5, C3H3, and C3H5

The ionization energies for methylene (CH2), methyl (CH3), ethynyl (C2H), vinyl (C2H3), ethyl (C2H5), propargyl (C3H3), and allyl (C3H5) radicals have been calculated by the wave-function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled-cluster level with single and double excitations plus a quasiperturbative triple excitation [CCSD(T)]. When it is appropriate, the zero-point vibrational energy correction, the core-valence electronic correction, the scalar relativistic effect correction, the diagonal Born-Oppenheimer correction, and the high-order correlation correction have also been made in these calculations. The comparison between the computed ionization energy (IE) values and the highly precise experimental IE values determined in previous pulsed field ionization-photoelectron (PFI-PE) studies indicates that the CCSD(T)/CBS method is capable of providing accurate IE predictions for these hydrocarbon radicals achieving error limits well within +/-10 meV. The benchmarking of the CCSD(T)/CBS IE predictions by the PFI-PE experimental results also lends strong support for the conclusion that the CCSD(T)/CBS approach with high-level energy corrections can serve as a valuable alternative for reliable IE determination of radicals, particularly for those radicals with very unfavorable Franck-Condon factors for photoionization transitions near their ionization thresholds.     

Ab initio direct dynamics studies on the reaction of H atom with CH3CH2Cl

The multiple channel reaction H + CH(3)CH(2)Cl --> products has been studied by the ab initio direct dynamics method. The potential energy surface information is calculated at the MP2/6-311G(d,p) level of theory. The energies along the minimum energy path are further improved by single-point energy calculations at the PMP4(SDTQ)/6-311+G(3df,2p) level of theory. For the reaction, four reaction channels (one chlorine abstraction, one alpha-hydrogen abstraction, and two beta-hydrogen abstractions) have been identified. The rate constants for each reaction channel are calculated by using canonical variational transition state theory incorporating the small-curvature tunneling correction in the temperature range 298-5000 K. The total rate constants, which are calculated from the sum of the individual rate constants, are in good agreement with the experimental data. The calculated temperature dependence of the branching fractions indicates that for the title reaction, H-abstraction reaction is the major reaction channel in the whole temperature range 298-5000 K.     

Ab initio VB Studies of the Ground and Low-lying Excited States of BeH and BH

A scheme has been proposed to classify valence bond(VB) wave functions for the calculations of ground and excited states,according to the symmetry properties of one-electron orbitals which are involved in the construction of VB wave functions.This scheme is illustrated by the examples of BeH and BH.Ab initio VB computations of these two test molecules in combination with the present classification scheme give reliable results.For example,calculation results show that the state C2∑ of BeH is stable,with the bonding energy 0.87 eV and bond length 0.238nm,which are in good agreement with those obtained by Gerratt et al.The bonding features of ground and low-lying excited states of BeH and BH are discussed.     

Ab initio studies on reaction H2C=C(OH)Li+CH3 + (CH3 -)

The possible structures and isomerizations of H₂C=C(OH)Li are studied theoretically by the gradient analytical method at RHF/6-31+G level. According to these results, reactions of H₂C=C(OH)Li with CH₃ ⁺ and CH ₃ ^- are investigated thoroughly. When H₂C=C(OH)Li reacts with CH ₃ ⁺ , H_zC=C(OH)Li firstly changes from structure1 to structure4, and then combines with CH₃ ⁺. In this reaction, the configuration of central carbon is retained. When H₂C=C(OH)Li reacts with CH ₃ ^- , structure1 firstly breaks its C-O bond to give contact ion-pair. Then through transition state16 which is similar to structure2, the attack of CH ₃ ^- from the opposite side of^-OH replaces^-OH group and inverts the configuration of carbenoid carbon atom. All the results show that the ambident reactivity of carbenoid has close relationship with the stability of special structures. Project supported by the National Natural Science Foundation of China (Grant No. 29773025).     

Ab initio calculations for ground states of CH2Li− and CH2Be

The ground state of CH₂Li^? and CH₂Be molecules has been investigated by an SCF calculation using a contracted Gaussian basis set. Only for the second system a bound state with respect to the ground states of the molecular fragments has been found.