Reaction paths of some open-shell reactive intermediates

The intermolecular interaction energy for reacting systems in singlet, triplet and doublet states was partitioned by the perturbation expansion method into the chemically meaningful five interaction terms: the Coulomb, exchange-repulsion, induction, dispersion, and charge-transfer energies. In the local ZDO approximation, these energy terms were evaluated for the dimerization of methylenes (^1,3CH₂), the additions of carbenes (^1,3CH₂ and^1,3CF₂) as well as amino radicals (²NH₂ and²NF₂) toward ethylene, and the hydrogen abstractions by methylenes (^1,3CH₂), nitrene (³NH), and hydroxyl radical (²OH) from methane. It has been found that the reaction path is much influenced by the spinmultiplicity, and that the charge-transfer and exchange-repulsion terms play a dominant role in determining the course of reactions.     

A theoretical study of electrophilic reactions

The electronic structures of protonated formyl and acetylium cations and their deprotonation paths leading to HCO⁺, COH⁺ and CH₃CO⁺have been studied by means of ab initio calculations. The results support Olah's theory that dipositive species can be the de facto reagents in electrophilic reactions.     

Ab initio SCF and CI study of the electronic spectrum of pyridine N-oxide

Configuration interaction (CI) studies of ground, n *, * * electronically excited states are reported for pyridine N-oxide. The transition energy to the lowest * excited ¹ B ₂ state is calculated at 4.35 eV, compared to the experimental spectrum range of 3.67–4.0 eV. This state lies below the lowest n * excited ¹ A ₂ state calculated at 4.81 eV above the ground state. The only experimentally reported triplet state at 2.92 eV above the ground state is predicted to be the ³ A ₁ (*) state. The calculated energy lies at 3.27 eV. Numerous other high-lying singlet states as well as the triplet states have also been calculated. The intramolecular charge transfer character of the ground and the excited states have been studied in terms of the calculated dipole moment and other physical properties.     

Theoretical study of the electronic structure of diazomethane

The least-energy dissociation path of the ground state of CH₂N₂ was determined fromab initio calculations using in a complementary way basis sets of minimal size (STO-3G) and double-zeta (DZ) quality. The results indicate that the least-energy point of attack of the N₂ molecule on CH₂ (¹ A ₁) is roughly perpendicular to the molecular plane (93 °), the C and N atoms being almost co-linear (angle C-N-N203 ° with outermost N atom pointing away from CH₂). The potential barrier of 1.2 eV found previously on theC _2v dissociation path, disappears completely along the least-energy dissociation path (point groupC _s (out-of-plane)). These findings corroborate the Woodward-Hoffman rules for this process since the outermost orbitals of the two intersecting states found in point groupC _2v (...2b ₁ and ...8a ₁) both correlate to the same irreducible representation (10á) in point groupC _s (out-of-plane).Larger basis set calculations (DZ + polarization functions on all centers, 3d _c and 3d _N developed here), were also carried out on CH₂N₂ (¹ A ₁,³ A ₂ and¹ A ₂) at the¹ A ₁ equilibrium geometry and on CH₂ (³ B ₁) and N₂ (¹ _g ⁺ ) at their respective equilibrium geometries. These calculations, together with consideration of correlation energy differences, yieldD ₀ ⁰ (CH₂N₂,¹ A ₁) = 19 kcal/mole and vertical excitation energies of 67 and 73 kcal/mole for the³ A ₂ and¹ A ₂ states respectively. The latter value is in good agreement with the measured experimental value: 72.4 kcal/mole corresponding to the maximum of intensity in the¹ A ₂¹ A ₁ absorption band.     

Band structure of even and ions of odd polyenes

The analysis of experimental data for singlet transitions (E _n) of even polyenes (I), cations (II) and anions (III) of odd polyenes show that for infinite chains E (I)/E (II)=E (I)/E (III) = 2:1. It is shown that the energy gap is equal for the three systems. In cases (II) and (III) there is a level (NBMO) in the gap which is vacant in (II) and occupied in (III). That is why the first optical transition in (II) and (III) depends on the semiwidth of the gap.     

Electronic states of the [2 n ]cyclophanes

Low energy singlet and triplet states for a series of [2 _n ]cyclophanes are discussed in terms of the results of a simple model calculation. Experimental trends can be explained under the assumption of significant - interaction involving the saturated bridges. This interaction destabilizes low energy excimer states, in contrast to the usual red shift observed for alkylbenzenes. The observed near-constancy of the onset of the absorption spectra can be explained by near-cancelation of through-bond and through-space contributions.From 01.01.1984, Royal Danish School of Educational Studies, Department of Chemistry, DK-2400 Copenhagen NV, Denmark     

Ab initio Calculations of the potential energy surface of the reaction of singlet methylene with the hydrogen molecule

The potential energy surface for the insertion of singlet methylene into H₂ has been computed on theab initio SCF level as well as with inclusion of electron correlation by means of the CEPA method. The results are compared with those of previous semiempirical,ab initio SCF and CI calculations. The system is a prototype of a reaction where an allowed and a symmetry-forbidden path can compete. The electron correlation energy was found to be very different for different regions of the surface, but did not have much influence on the optimum reaction path. From the computed heat of the reaction, the heat of formation of singlet methylene was estimated to be 101.5 kcal/mol. According to the calculations the reaction does not need any activation energy.     

Theoretical studies on the conformation of saccharides

Conformations of 2-methoxytetrahydropyran as a model for the six-membered ring in aldopyranosides have been calculated by the PCILO method using the algorithm of the conjugated gradient to optimize the geometry. The calculated geometry of the fourteen basic forms of 2-methoxytetrahydropyran was found to be in agreement with the available data obtained by X-ray diffraction of pyranosides. The results indicate differences in the geometry of 2-methoxytetrahydropyran resulting from the change of the axial vs. equatorial position of the methoxyl group. These changes are particularly meaningful in the values of bond angles and they are in agreement with the anomeric and exoanomeric effects. The experimentally found differences in the energies of an axial (⁴ C ₁) and equatorial (¹ C ₄) conformer, G = 2.9–3.7 kJ/mol, and the dipole moment, = 1.20 ± 0.05 D (1D = 3.33 10^–30mAs) agree well with the calculated values E = 3.18 kJ/mol and <> = 1.18 D which, in turn, suggest that the axial conformer is preferred over the equatorial one by a ratio a:e = 78:22.     

Multipole expansion of diatomic overlap

The systematic extension of Ruedenberg's expansion formula proposed in Part I [1] is applied to a series of diatomic and polyatomic molecules (BH, NH, HF, Be₂, C₂, F₂, CO, BH₃, CH₄, NH₃, H₂O, HCN and H₂CO). In general, good agreement with the results of full SCF calculations with the same minimum STO basis set is achieved. Thus, the errors due to this integral approximation scheme called MEDO (Multipole Expansion of Diatomic Overlap) are almost negligible compared to those introduced by basis set truncation.     

Triplet states via intermediate neglect of differential overlap: Benzene,pyridine and the diazines

The intermediate neglect of differential overlap technique is modified and applied to the calculation of excited triplet states. The resulting method generally reproduces the transition energies of the better-classified observations within a rms error of 1000 cm^–1. Trends are well reproduced, and the calculated orders ofn-* and -* triplet states are in good accord with the experimental information to date.The method is applied to benzene and the azines. The lowest four triplet states of benzene are calculated in good accord with experiment. Pyridine is calculated to have an-* triplet nearly degenerate with the lowest lying -* triplet, corroborating suggestions of Japar and Ramsay based on experimental information. A detailed analysis is made of the diazines, and assignments are suggested for the higher lying triplet states not yet classified or not yet observed.     

Numerical integration of overlap and ligand contributions to the electric field gradient

The total electric field gradient (EFG) tensor V _pq is calculated by numerical integration of threedimensional integrals. Each of them is solved a) by integrating over one dimension analytically and b) by integrating over the remaining two dimensions on the basis of a Gauss-type integration rule. The use of 100 abscissas in the twodimensional numerical integration scheme yields satisfactory accuracy which was checked by evaluating overlap integrals; an increase to 400 abscissas does not increase the result drastically. Calculating quadrupole splittings E _Q from numerically integrated electric field gradient tensors V _pq we observe that depending a) on the amount of covalency and b) on the amount of deviation from octahedral or tetrahedral symmetry, involved in a molecular system, overlap and ligand contributions to V _pq play an important role. Especially for the sandwich compound ferrocene, Fe(C₅H₅)₂, we find a significant difference between E _Q ^{num. int.} which follows from the numerical integration method, and E _Q ^conventional which is derived from effective charges.Supported by Deutsche Forschungsgemeinschaft     

Clusters of hydrogen molecules: Ab initio SCF calculations corrected semiempirically for correlation energies

Stabilization energy of the (H₂) _n clusters (n = 2–8) was calculated as a sum of the SCF interaction energy and the semiempirical interaction correlation energy estimated according to Sinanolu and Pamuk. Optimum successive attachment of hydrogen molecules leads to the formation of a gas-phase solvation shell consisting of seven hydrogen molecules. Basis set effect has been found to be important with all clusters under study. The non-additivity effect was investigated with the (H₂)₄ cluster. Vertical ionization potentials of the clusters considered are predicted to be 0.4–0.6 eV lower than the ionization potential of the parent H₂ molecule.     

On the affinity of cytosine towards electrophiles

Ab initio SCF computations on the intrinsic preferences of the H⁺, CH ₃ ⁺ and C₂H ₅ ⁺ cations towards the two principal sites of protonation or alkylation on cytosine, N₃ or O₂, show that this preference undergoes a continuous modification with the increase in size and complexity of the cation. N₃ is the preferred site of fixation of H⁺, O₂ the preferred site of C₂H ₅ ⁺ , while CH ₃ ⁺ has no marked preference. The exchange repulsion term of the binding energy appears responsible for the preference of C₂H ₅ ⁺ for O₂.This work was supported by the Ligue Francaise contre le Cancer and the National Foundation for Cancer Research (USA)     

Ab initio Studies on polymers

Ab initio crystal orbital calculations have been performed on the infinite all-trans polyene. A structure with r _c=c= 1.346 Å, r _c-c = 1.446 Å, r _c-h = 1.08 Å, and CCC = 125.3 ° was found to be most stable. The most important force constants, the band structure and the density of states were determined as well.     

Non-additivity in water-ion-water interactions

The three-body system Li⁺(H₂O)₂ was analyzed to study that non-additive part of the interaction potential which can be obtained by the Hartree-Fock approximation.For long and intermediate distances the three-body correction was found to be well represented by the induction energy, where bond dipoles are induced on each water molecule by point charges located on the (unpolarizable) lithium ion and on the other molecule respectively: for shorter distances this approximation was corrected by means of an exponential repulsive term. Such a potential model for non-additive interactions was extended to the more general situation Li⁺(H₂O)_n, and Monte-Carlo calculations were carried out on clusters containing up to six water molecules; comparison with other simulation results and with available data showed a significantly improved agreement with experiment. Tentative values for H are presented for n =7, 8,..., 20, where experimental data are not available.     

Calculation of harmonic force constants,vibrational frequencies and integrated intensities of some organic molecules using the MINDO—Forces method

The harmonic force constants, vibrational frequencies and integrated intensity ratios of CH₂, H₂O, CH₂O, C₂H₂, CO₂, HCN, CH₃, CH₄, and C₂H₄ have been calculated using the MINDO—FORCES program and the Pulay method for the calculation of the molecular force constants. The results obtained are in general quite satisfactory when compared with available literature values. The results are, however, not as satisfactory in case of molecules containing heteroatoms, due to the neglect of some dipolar repulsion integrals for the heteroatoms by the MINDO/3 method. Calculated integrated intensities for CH₃ and C₂H₄ agree well with experimental results. The calculated integrated intensities for other molecules are obtained for the first time and no comparison with published data is therefore possible.Part of the M.Sc. Thesis of K. H. A. 1978.     

Excited states of phenyl carbonyl compounds

In an attempt to clarify the origin of the dual phosphorescence in phenyl alkyl ketones, we have made some calculation (within the C.I.P.S.I. method in an excitonic scheme) to elucidate the conformation of both ground states and excited states of propiophenone. Our calculations have shown the presence of two stable isomers in the ground state, first n ^* state, and first ^* singlet and triplet states. So our work suggests that the origin of the dual phosphorescence of propiophenone could be related to the conformational change of the molecule in the n ^* state, because the most stable conformations in the n ^* state and in the ground state are different.     

Screened potential in π-electron systems and its application to the estimation of excitation energies

A diagrammatic technique was developed for the estimation of the screened potential of -electron systems. The screened potential was expanded in terms of the polarization propagators which were constructed from either the singlet, , or triplet vertex part, . These vertex parts correspond to the singlet or triplet excitations, respectively, in the Random Phase Approximation (RPA) containing exchange diagrams. The excitation energies were calculated by using the screened potential in the framework of RPA with exchange. The excitation energies of several conjugated molecules with or without a hetero atom are in agreement with the experimental data.     

Conformational analysis of dimethylphosphate with quantum-mechanical and classical methods

In order to shed light on the conformational behavior of polynucleotide chains, and in particular to clarify the origins of the barriers to internal rotation in the phosphodiester linkage, we computed, with a quantum-mechanical ab initio procedure, the energies associated to 86 combinations of the two torsion angles in the dimethylphosphate anion (CH₃O)₂PO₂ ^–, and then we sought for an analytical expression apt to reproduce these energies with the highest possible accuracy. An excellent agreement (standard deviation of the fitted energies from the ab initio energies 0.28 kcal/mole) with the quantum-mechanical calculations was reached with a potential consisting of four terms: 1) a 6–12 Lennard-Jones contribution, in which different parameters are used to describe the interactions of methyls with the ester oxygens and with the anionic oxygens; 2) a contribution with twofold periodicity, accounting for the anomeric effects connected to the interactions between the lone pair electrons and the polar bonds of phosphorus with the anionic oxygens; 3) a contribution with threefold periodicity, representing the usual bond-staggering term; and 4) a Coulombic contribution, arising from electrostatic interactions between partially charged atoms.