Ab initio valence bond calculations. I. Methylene

An ab initio valence bond (VB ) method for calculating energies and wave functions for simple electronic systems is described and the results of its application to methylene are given together with extensive comparisons with previous theoretical results.     

Ab initio valence bond calculations and the spin-paired diradical character of S 4 2+

The results of some ab initio valence bond calculations, with STO-6G basis sets for the s and p orbitals, are reported for the ground state of cyclic S ₄ ²⁺ . The sum of the weights for two long-bond (or spin-paired diradical) structures is approximately 50% of the total.     

Ab initio valence bond calculations. X. Vertical valence ionization potentials of allene and butatriene

The ground and vertical valence ionized states of allene and butatriene have been studied in the ab initio valence bond framework using the 6–31G basis set after contraction and introducing the core–valence shell separation. The final wave functions have been analyzed in terms of VB structures by means of population analysis.     

Ab initio valence bond calculations. XI. Pyridine: Ground and ionized states

The ab initio valence bond method has been used to study the ground and the lowest vertical valence ionized states of pyridine. On the basis of our calculations the first two ionization potentials are assigned to π and n electron removals, respectively. The final wave functions have been interpreted in terms of valence bond structures by means of an appropriate population analysis.     

Ab initio calculations on H transfer in the HF trimer

Systematic ab initio LCAO MO calculations with basis sets of increasing size have been performed to analyze the changes in relative stabilization energies of the two cyclic (D_3h and C_3h symmetry) and the noncyclic HF trimers. The results form the basis for a discussion of cyclic H-transfer reactions.     

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.     

Ab initio MO calculations on 1,2-dithietes and valence isomers

The valence isomerization of the 1,2-dithiete parent compound to the open-chain dithial was studied by CASSCF multiconfiguration methods including the CASPT2 perturbational treatment. The isomerization energy remains small at the highest level of theory. In agreement with Jonas and Frenking, the cyclic structure is only then preferred over the acyclic ones if f-functions on the sulfur atoms are considered. If they are included, the 1,2-dithiete is more stable by 3.8 kcal/mol and the barrier amounts to 24.9 kcal/mol at the CASPT2(8,8)/6-31G(2df)//MP2/6-31G* level of theory. According to MP2/6-31G* geometry optimizations, substitution of H by NH₂ and CH₃ reduces the stability of the 1,2-dithiete ring structures relative to the open-chain dithiocarbonyl structures, whereas the inverse holds for acceptor substitution by CN and CF₃. A higher stability of benzodithiete relative to ortho-dithiobenzoquinone is predicted at all employed levels of theory. This is in good agreement with conclusions drawn from experimental results. Whereas the experimental microwave geometry of 1,2-dithiete is well reproduced theoretically, the experimental electron diffraction geometry of 3,4-bis-(trifluoromethyl)-1,2-dithiete differs from the calculated one. © 1996 John Wiley & Sons, Inc.     

Benzooxirene. Ab initio calculations

Ab initio calculations have been performed on benzooxirene, the corresponding oxo carbene (“ketocarbene”), and the transition state linking the two. At the highest level used, QCISD(T)/6-31G^*//MP2(FULL)/6-1G^* with MP2(FULL)/ 6-31G^* zero point energy corrections, the relative energies of the oxirene, the transition state and the carbene are 0, 24.6, and −17.8 kJ mol⁻¹. Correlation energy effects are very important in this system: at the QCISD(T) level the oxirene lies above the carbene, as at the MP4 and HF levels, but at the MP2 level the ordering is reversed. Benzooxirene is probably slightly nonplanar: the HF/6-31G^* geometry is C_2v but the MP2(Fermi contact)/6-31G^* geometry is C_s with a 6-/3-ring coplanarity deviation of about 6.9 °, although in the MP2(FULL)/6-31G^* geometry this is reduced to about 3.1 °.     

Ab initio investigation of the potential energy surfaces of C2H2F2 and C2H2F2+

Ab initio MO calculations have been performed for neutral and cationic C₂H₂F₂ structures. Olefinic and carbene structures are investigated for the neutral isomers, while olefinic, carbene, and fluoronium-type cations are found. Stability orders and rotational barriers are discussed in terms of orbital and Coulomb interaction. Contrary to previous studies, the higher stability of the geminal isomers is interpreted to be caused by Coulomb attraction.     

Ab initio HF and DFT calculations on an organic non-linear optical material

The molecular geometric optimization, vibrational frequencies, and gauge-including atomic orbital (GIAO) ¹H and ¹³C chemical shift values of 3-[(1E)-N-ethylethanimidoyl]-4-hydroxy-6-methyl-2H-pyran-2-one have been investigated by using ab initio Hartree–Fock (HF) and density functional method (B3LYP: Becke-3-Lee–Yang–Parr) with 6–31G(d) and 6–31++G(d,p) basis sets. Also, the first hyperpolarizabilities have been calculated at the HF and B3LYP levels employing the corresponding basis sets. To understand this phenomenon in the context of molecular orbital picture, we examined the molecular HOMOs and molecular LUMOs generated via HF and B3LYP levels. The computed vibrational frequencies are used to determine the types of molecular motions associated with each of the experimental bands observed. Data of 3-[(1E)-N-ethylethanimidoyl]-4-hydroxy-6-methyl-2H-pyran-2-one display significant second-order molecular nonlinearity and provide the basis for design of efficient nonlinear optical materials.     

水分子垂直电离能的从头计算

分子的芯电子和价电子电离能可用XPS和UPS测得,并可作理论上的近似计算,其中ab initio Hartree-Fock SCF LCAO-MO法是目前广泛采用的一种方法.它所作的三个近似(非相对论近似、Born-Oppenheimer近似和单电子近似)很明确.除了几个基本物理常数外,不再引进任何参数.但是计算中选用的基函数对所得结果的精确度有较大影响.因此,常常需要以扩大的Slater型函数(STF)或Gaussian型函数(GTF)为基函数.由于所需计算的双电子积分的数目正比于基函数数目的四次方,这样就会使计算化费的计算机时间大大增     

Ab initio cluster calculations for the absorption energies of F and F+ centers in bulk ZnO

We present the results of a series of ab initio calculations on the ground states and the low lying excited states of the F and F+ centers in bulk ZnO. Both types of F centers are oxygen vacancies, causing rather strong distortions of the local geometries. The calculations were performed by means of wave function based methods, mostly at the CASSCF level. Dynamic correlation was included for the first two coordination shells of the F centers. The calculated absorption energy for the F+ center (3.19 eV) is in excellent agreement with the experimental value of 3.03 eV. For the emission from the 3T2 state of the F center to the 1A1 ground state we obtained a transition energy of 2.73 eV. Experimentally, a green photoluminescence is observed at 2.38-2.45 eV. We estimated that the errors in our calculation should be even smaller in the latter case than for the F+ state, where the calculated transition energy differs by less than 0.2 eV from the experimental value. Therefore, we assume that the 3T2 to 1A1 transition is not the origin of the green luminescence.     

Ab initio evaluation of bond orders and valence numbers in some substituted benzenes

Mayer's definition of bond order and valence number of an atom in a molecule in the ab initio SCF theory have been applied to calculate bond orders and valence numbers using different ab initio methods for some fluorobenzenes. The results have been discussed in the light of Mulliken population analysis. The traditional view that the valency of an atom is directly related to the atomic charge is found to be invalid in the present ab initio calculations.     

Ab initio investigations of the bound rovibrational levels of NeH 2 + , NeHD+, and NeD 2 +

Summary Bound rovibrational levels have been calculated for NeH ₂ ⁺ , NeHD⁺, and NeD ₂ ⁺ using three recent fits to an accurateab initio PES. The NeH ₂ ⁺ molecule behaves essentially as a linear molecule, the predicted rotational constant is 2.57 cm^–1. The fundamental frequencies are 811, 1189, and 1748 cm^–1 for the Ne-H ₂ ⁺ stretch, the Ne-H ₂ ⁺ bend and H ₂ ⁺ stretching modes, respectively.Dedicated to the 60th birthday of Prof. W. Kutzelnigg, Bochum     

Ab initio calculations on sulfur-containing compounds. II. One-electron properties of H2S

Closed-shell RHF one-electron properties are calculated for H₂S using a total of 41 different s, p basis sets and two polarized basis sets (6–31G* and 6–31G**). Total energies and geometries alone are not a comprehensive criteria for selecting the best basis sets. It is shown here that the comparison of a number of one-electron properties can serve as an excellent criteria for testing basis sets. The quality or reliability of a basis set is taken as being its agreement with a large uncontracted s, p basis set (s, p limit).     

Ab initio calculations of doublet states of NH+

The eight low-lying doublet states of the NH⁺ ion are investigated with an ab initio configuration interaction method including all single and double excitations from a multi-reference configuration space (MRSD CI). The spectroscopic constants for the X²Π, A²Σ⁻,B²Δ and C²Σ⁺ states and the transition moments for X²Π-A ²Σ⁻¹ and X²Π-B²Δ are calculated. The results are compared with experiments and other calculations.     

Ab initio direct dynamics trajectory simulation of C2H5F-->C2H4 + HF product energy partitioning

Direct dynamics classical trajectory simulations were performed to study product energy partitioning in C(2)H(5)F-->C(2)H(4)+HF dissociation. The intrinsic reaction coordinate potential energy curve, reaction energetics, and transition state (TS) properties were calculated for this reaction at different levels of electronic structure theory, and MP2/6-31G( *) was chosen as a meaningful and practical method for performing the direct dynamics. The trajectories show that the HF bond, uncoupled from the other degrees of freedom, is formed within the first 10 fs as the system moves from the TS towards products. The populations of the HF vibration states, determined from the simulations, decrease monotonically as found from experiments. However, the simulation's populations for the low and high energy vibration states are larger and smaller, respectively, than the experimental results. The HF rotational temperature found from the simulations is in agreement with experiment. Increasing the TS's excess energy gives higher rotational temperatures for both C(2)H(4) and HF. Energy is partitioned to the products from both the excess energy in the TS and the potential energy release in the exit channel. Partitioning from these two energy sources is distinguished by varying the TS's excess energy. An analysis of the simulation's energy disposal shows that the fractions of the excess energy partitioned to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation, are 0.17, 0.64, 0.076, 0.067, and 0.046, respectively, and are in good agreement with previous simulations on empirical potentials and experiments. The partitioning found for the potential energy release is 81%, <0.05%, 5%, 11%, and 3% to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation. This result is substantially different than the deduction from experiments, which summarizes the partitioning as 20%, 45%, 24%, and <12% to relative translation, C(2)H(4) vibration+rotation, HF vibration, and HF rotation. Possible origins of the difference between the simulations and experiments in the release of the potential energy is discussed.