Continuum states of the hydrogen molecule with the R-Matrix method

The interaction of low-energy electrons with H and ligh with H₂ is analyzed using the R-matrix method including both coupled state and polarization effects ab initio. Particular attention is paid to the energy region containing low-lying H*₂ resonances which converge to H A ²Σ. Resonances in both elastic scattering and the photoionization asymmetry parameter, β, are presented. At low photon eneries, satisfactory results for vibrationally resolved photoionization are obtained using the adiabatic nuclei approximation. The use of nonadiabatic techniques for higher energies is discussed.     

Ab initio studies of negative ion-molecule(s) clusters present in the atmosphere. II. OH−(H2O)n for n = 0, 2, 3, 4

Ab initio calculations are performed with 6–31G basis set to study the geometry and binding of the H₃O, H₅O, H₇O, and H₉O complexes. The H₃O complex is also investigated with the 6–31 G* basis set and MP 2 (Moller–Plesset perturbation theory of second order).     

The formation of highly excited H in the reaction H2+(v) + H2 → H + H

A quasiclassical trajectory surface hopping method has been used to study H(v) + H₂ → H + H for v = 0, 3, 7, 10, 13, and 17 with an emphasis on determining the H internal energy and angular momentum distributions for high v. For v = 13 and 17, significant cross sections are found for producing H at energies above its dissociation energy. An average metastable H lifetime of 11.5 ps for v = 13 and 4.7 ps for v = 17 is found, but there is also a much longer lived component to the lifetime distributions that is more important for v = 13 than for v = 17. Some of the longer lived metastables correspond to high angular momentum orbiting states of H, but other sources of metastability are also present.     

The proton affinity and the structure of protonated species of methylsilane

An ab initio LCAO-MO-SCF calculation was made on the proton affinity (PA ) of methylsilane (CH₃SiH₃) by using STO -3G, MIDI -1, and MIDI -1* basis sets. Three types of protonated methylsilane are taken into account, and their geometrical parameters are optimized. The calculated PA of CH₃SiH₃ is 160.5 kcal/mol, which exceeds that of SiH₄ by 11.5 kcal/mol. The protonated species (I) which refers to Si—C bond protonation is shown to be most favorable, and to be a weak σ-complex between CH₄ and SiH. Other two species are also σ-complexes between H₂ molecule and SiH₃CH or CH₃SiH, and similar to CH, SiH, GeH, and C₂H.     

SCF and MC–SCF Study of CO2 with bond angle variation

Multiconfiguration (MC ) SCF calculations are reported for CO₂ for bond angles between 60° and 180°. The ground state configuration is found to be …?5a4b∣b∣a for small bending angles and …?6a3b∣b∣a for large bending angles, the change in ground state character occurring at a bond angle of about 100°. The force constant for bending obtained from the MC –SCF function is about 8.0% lower than the corresponding SCF value, and in considerably better agreement with experiment.     

Convergence accelerator approach for the evaluation of some three‐electron integrals containing explicit rij factors

A simplified analysis is presented for the evaluation of the three‐electron one‐center integrals of the form ∫rrrrrred r ₁d r ₂d r ₃, for the cases i, j, k, ≥−2, l=−2, m≥−1, n≥−1. These integrals arise in the calculation of lower bounds for energy levels and certain relativistic corrections to the energy when Hylleraas‐type basis sets are employed. Convergence accelerator techniques are employed to obtain a reasonable number of digits of precision, without excessive CPU requirements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 93–99, 1999     

Dipole,quadrupole, octupole,and dipole–octupole polarizabilities at real and imaginary frequencies for H,HE, and H2 and the dispersion-energy coefficients for interactions between them

We have calculated certain dynamic polarizabilities (for both real and imaginary frequencies) for H, He, and H₂ and the dispersion-energy coefficients for long-range interactions between them. We have done so in a sum-over-states formalism with explicitly electron-correlated wave functions to describe the states. To be precise, we have determined the dipole (α₁), quadrupole (α₂), and octupole (α₃) polarizabilities of H and He for real frequencies (ω) in a range between zero and the first electronic-transition frequency and for imaginary frequencies (iω) on a 32-point Gauss-Legendre grid running from zero to ?ω = 20 E_h, and for H₂, we have found the dipole (α), quadrupole (C), and dipole–octupole (E) polarizability tensors for the same real and imaginary frequencies. The dispersion-energy coefficients, obtained by combining the sum-over-states for-malism for the polarizabilities with analytic integration over ω, gave values of C₆, C₈, and C₁₀ for the atom–atom systems; C, C, C, C, and C for the atom–diatom systems; and C₆, C and C for the H₂? H₂ system. Nearly all the results are considered to be more reliable than those hitherto published and some have been obtained for the first time, e.g., C(iω), E(ω), and E(iω) for H₂ and C, C, and C for the H? H₂ system. © 1993 John Wiley & Sons, Inc.     

Atomic integrals containing rrr with λ, μ, ν ≥ −2

In this paper, the efficient evaluation of the atomic integrals I =∫rrrrrre^{?αr1?βr2?γr3}dτ with one or two factors r is described. These integrals are necessary for a lower-bound calculation for Li-like systems using the method of variance minimization or Temple's formula. © 1993 John Wiley & Sons, Inc.     

The spin-density-functional formalism for quantum mechanical calculations: Test on diatomic molecules with an efficient numerical method

We have applied the spin-density-functional (SDF ) formalism with the local-spin-density (LSD ) approximation to a number of small molecules with the primary aim of testing the approximation for molecular applications. A new numerical method to solve the one-electron wave equation is developed, utilizing the special features of the SDF formalism. We have calculated energy curves, dissociation energies, and equilibrium distances for some diatomic molecules [H (²Σ, ²Σ), H₂(¹Σ, ³Σ), He (¹Σ), and He₂(¹Σ)] and the vibrational frequencies of H₂. The deviations from the experimental results are typically 1/2 eV for the energies and ≤ 0.1 Å for the distances. We discuss the LSD approximation using the concept of an exchange-correlation hole and make predictions about the applicability to other molecules. The LSD approximation is compared with the Hartree-Fock and multiple-scattering-Xα methods and some difficulties in the latter methods are pointed out. It is argued that the SDF formalism within the LSD approximation has physical advantages compared to the Hartree-Fock and Xα methods and that it should provide a simple and useful method for a broad range of applications.     

Ab initio MO calculations on the stability of the CF,PF, SF,and MoF ions with the F1s core hole

For the CF, PF, SF, and MoF ions appearing after the F1s photoionization, the possibility of dissociation has been shown by the ab initio MO LCAO method within the Z + 1 core equivalent model. According to the calculations, the decay channel AF → AF + F(1s¹2p⁶) is energetically open for the ions. So the interpretation of the gas-phase emission FKα spectra, in which the bands are assigned to the discrete transition energies, can be unacceptable for these ions. The conditions and signs of such failure are discussed.     

The role of d orbitals in tetrahedral hybrid orbitals of oxy-ions

The hybrid orbitals of tetrahedral oxy-ions containing some d character have been calculated by maximum overlap method. The d characters of hybrid orbitals increase in the order of SiO, PO, SO, ClO, and decrease in order of GeO, AsO, SeO, BrO. The bond strengths are also obtained for these ions. The hybrid Orbital of VO, CrO, and MnO are of the type d³s as the result of calculation.     

Ab initio calculation of electronic states of the ions H4+ and H5+

Electronic state calculations for the ions H₄⁺ (with symmetries D ₄ and C _2v) and H (with symmetries D ₅ and D _2d) are made using the valence-bond method. All the configurations obtained from the given set of 1s-functions of Slater type are taken into account. Space functions are used throughout the computation (“spin-free quantum chemistry”). Preliminary quasidiagonalization of the secular equation is implemented by the construction of the multiplet eigenfunctions ^2S+1Γ^(α) from the initial variational functions. The results of the calculations are as follows: the ion H is unstable, the ion H is stable with equilibrium nuclear conformation of symmetry D _2d and with the energy of dissociation into H and H₂ near 4 eV.     

Polysila analogs of aromatic hydrocarbon ions: Structures and energies of Si3H3+, Si4H,and Si5H5−

Silicon analogs of aromatic monocyclic ions, (SiH) ( 4 ), (SiH) ( 5 ), and (SiH) ( 6 ) have been studied ab initio at MP 2(full)/6-31G *. The D_3h structure of Si₃H₃⁺ is the global minimum, whereas other two ions are nonplanar. The D_2d structure of (SiH) is less folded than the carbon analog and possesses a higher stabilization energy. Stabilization energies for the monocharged ions are diminished with respect to the corresponding carbons © 1993 John Wiley & Sons, Inc.     

Internal stresses and chemical binding in the H molecule

Following an interpretive formalism presented earlier, chemical binding in the H molecule has been studied in terms of the variations, with respect to R, of electrostatic field and stress components at five selected points on the interuclear axis. At three points phenomena analogous to those recorded earlier for H₂ have been observed. In particular, the existence of extremal relationships for the difference density field and the total field, as well as for the corresponding stresses, at R ? R_eq for certain specific points on the internuclear axis has been confirmed. As in the case of H₂, chemical binding in H occurs due to local variations of electrostatic pressure from point to point in such a manner as to cause the vanishing of either the total electrostatic force density or the difference force density at certain points on the internuclear axis.     

Exclusion principle and the symmetry adapted wave-functions

Wave-functions of various spin-dependent and spin-free methods are examined from the point of view of the requirements of the exclusion principle and the spin projection. It is shown that the “two-rowed” or “two-columned” requirements of the standard Young tableaux are necessary but not sufficient to replace the “antisymmetry” requirement of a wave-function and to be regarded as exclusion principle. The symmetry adapted wave-functions which are constructed from the matric basis e may not satisfy the exclusion principle and, hence, their usages are open to question. The appropriate symmetry adapted wave-functions which satisfied the exclusion principle are given for any pure spin state. We have also shown that the structure operators NP for bond functions are spin projectors in the Waller-Hartree double-antisymmetrized space only, and should not be used in the Hartree product space. Furthermore, if the corresponding matric operators PNP are used in the Hartree product space, then the wave-functions thus constructed may not be antisymmetric with respect to the permutations of indistinguishable electrons.     

Perturbational calculation of the exchange forces in the two lowest states of the hydrogen molecule

The interaction energies in the ground (X¹Σ) and the first excited (b³Σ)states of the hydrogen molecule were calculated using the modified Jansen and Byers Brown perturbation method. Calculations were performed for nine values of the internuclear distance R in the range 5 a.u. ≦ R ≦ 9 a.u. The present results were compared with the interaction energies calculated by four other perturbation methods as the expectation values of the H – E⁰ operator with the function accurate up to the first order. The results show that the method can give satisfactory values of the interaction energies.     

Binding of an electron to a molecular dipole: BeF

Results of extended-basis SCF calculations indicate that BeF may exist as a metastable species. Comparison of results obtained from SCF calculations on neutral BeF₂ with those of BeF shows that the orbital occupied by the electron of BeF is well approximated by the lowest unoccupied molecular orbital of neutral BeF₂.     

Density functional calculations of structures and ionization energies for heavy group V cluster anions

The structure and electronic structure of heavy-group V cluster anions (Sb, Bi) are calculated with density functional methods within the local spin density approximation (LSDA ). The influence of gradient corrections of the exchange and correlation energy is investigated. The calculated vertical and adiabatic ionization energies are in very good agreement with data from photoelectron spectroscopy (PES ) for Sb, whereas the relatively large deviations for Bi can be reduced by the consideration of relativistic effects in a scalar-relativistic manner. Concerning the structures, a strong similarity to the corresponding P clusters was found. In particular, the negatively charged pentamers are planar rings (with similarities to the aromatic [C₅H₅]^? anion) with especially high ionization energies. © 1995 John Wiley & Sons, Inc.     

The role of the bond midpoint electron density in homonuclear molecular binding

Electron density theory is first employed to express the potential energy curve of the H molecular ion in terms of bond midpoint properties. For heavy homonuclear diatomics, low-ordr density gradient theory is used, but now for the chemical potential, with a similar conclusion to that for H. Homonuclear clusters of alkali atoms are then treated in some detail. Finally the dissociation of doubly charged clusters is considered, by a study of supermolecular treated in some detail. Finally the dissociation of doubly charged clusters is considered, by a study of supermolecular ions (Na)₂ and (K)₂, again using a density gradient expansion in low order. The deviation between the barrier for fission and the Coulomb barrier is linear in the bond midpoint density over a substantial range of fragment separation. © 1994 John Wiley & Sons, Inc.