The lowest order relativistic corrections for the hydrogen molecule

We have calculated the lowest order relativistic corrections for the X ¹Σ, B¹Σ, a³Σ, b³Σ, I¹Π_g, C¹Π_u, i³Π_g, c³Π_u, J¹Δ_g, and j³Δ_g states of the hydrogen molecule using variational Monte Carlo methods and compact, explicitly correlated trial wavefunctions. Our values are in good agreement with earlier calculations on the X¹Σ and B¹Σ states. For the other states, our work provides the first estimate of these properties. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Comparison of theory and experiment for excited singlet states of the H2 molecule

This paper presents a survey of published and unpublished ab initio calculations of the vibrational structures of the ten lowest electronic singlet states of the hydrogen molecule up to the H(n = 1) + H(n = 2) dissociation limit. The data are based on adiabatic potential functions (clamped-nuclei electronic energies and nuclear-mass-dependent diagonal corrections). Nonadiabatic coupling has been treated ab initio within the five states. of ¹Λ symmetry (X,EF, GK, HH?) and ¹Σ I.¹Π_g. The accuracies of the theoretical energies are determined by comparisons with experimental data for H₂, HD, and D₂. The level shifts and predissociation probabilities of the excited ¹Σ states, generated by nonadiabatic coupling with the discrete and continuous vibrational structure of the ground state, and radiative properties have also been calculated.     

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 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.     

Theoretical spectroscopic studies of InI and InI+

Relativistic configuration interaction calculations are carried out to study the electronic structure and spectroscopic properties of InI and InI⁺. Potential energy curves of the ground and a number of low‐lying states are constructed. Spectroscopic parameters of the bound states of both species are computed and compared with the experimental and other theoretical data. Effects of spin‐orbit coupling on the spectroscopic properties are studied. Because of the presence of the heavy atoms the effect is large. The spin‐orbit splitting of the ground state (X²Π) of InI⁺ is more than 8350 cm^?1. As a result of the strong spin‐orbit interaction between X²Π and A²Σ⁺ of InI⁺, the potential energy curve of A²Σ becomes repulsive. Radiative lifetimes for the spin‐forbidden transitions such as A³Π?X¹Σ and B³Π₁ ?X¹Σ of InI and spin‐allowed transitions such as B²Σ⁺?A²Σ⁺, C²Π?A²Σ⁺, and B²Σ⁺?X²Π are calculated. Vertical and adiabatic ionization energies of InI and the electric dipole moments of both the neutral and ionic species are estimated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

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.     

Evaluation of two‐center,three‐ and four‐electron integrals over Slater‐type orbitals in elliptical coordinates

An algorithm for evaluation of two‐center, three‐electron integrals with the correlation factors of the type rr and rrr as well as four‐electron integrals with the correlation factors rrr and rrr in the Slater basis is presented. This problem has been solved here in elliptical coordinates, using the generalized and modified form of the Neumann expansion of the interelectronic distance function r for k ≥ ?1. Some numerical results are also included. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

D∞h B2(BS)2−/2− and Td B(BS)4−: Boron sulfide clusters containing BB multiple bonds and B− tetrahedral centers

A density functional theory investigation on the geometrical and electronic properties of B₄S (B₂(BS)) and B₅S (B(BS)) clusters has been performed in this work. Both the doublet B₂(BS) ([S?B? BB? B?S]^?) (D_∞h, ²Π_u) and the singlet B₂(BS) ([S?B? B?B? B?S]^2?) (D_∞h, ¹Σ) proved to have perfect linear ground‐state structures containing a multiply bonded BB core (BB or B?B) terminated with two BS groups, while T_d B(BS) turned out to possess a perfect B^? tetrahedral center directly corrected to four BS groups, similar to the corresponding boron hydride molecules of D_∞h B₂H, D_∞h B₂H, and T_d BH, respectively. B₄S₂ and B₅S₄ neutrals, however, appeared to be much different: they favor a planar fan‐shaped C_2v B₄S₂ (a di‐S‐bridged B₄ rhombus) and a planar kite‐like C_2v B₅S₄ (a di‐S‐bridged B₃ triangle bonded to two BS groups), respectively. One‐electron detachment energies and symmetrical stretching vibrational frequencies are calculated for D_∞h B₂(BS) and T_d B(BS) monoanions to facilitate their future characterizations. Neutral salts of B₂(BS)₂Li₂ with an elusive B?B triple bond and B(BS)₄Li containing a tetrahedral B^? center are predicted possible to be targeted in experiments. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

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.     

An insight into optical and EPR properties of AgCl and AgF complexes through MS–Xα and SCCEH calculations

MS-Xα and SCCEH calculations on the Ag²⁺ complexes AgF and AgCl (displaying an elongated D_4h symmetry) have been carried out for a better understanding of their experimental optical and EPR properties. As salient features, the present work supports that the unpaired electron in AgCl spends a little more time on ligands than on Ag²⁺, in agreement with the previous analysis of EPR and optical data for KCl:Ag²⁺. Furthermore, the five experimental optical transitions observed in that case are reasonably assigned. The first transition (observed at 12,500 cm^?1) is assigned to a jump involving the 5a_1g orbital built mainly (∽70%) from 3p orbitals of axial ligands, a fact that reflects the distinct level scheme for AgCl when compared to that for more ionic complexes. Calculations on AgF and AgF performed as a function of the equatorial Ag²⁺ –F^? distance led to a reasonable understanding of experimental gyromagnetic and superhyperfine tensors displayed by Ag²⁺ in fluorides. The different relative decrease undergone by g‖– go (8%) and g ? – go (28%) on passing from CsCdF₃:Ag²⁺ to RbCdF₃:Ag²⁺ is shown to be consistent with the formation of AgF and AgF complexes, respectively, related to the different substitutional position of Ag²⁺ in such lattices. The decrement of about 8.5% experienced by both g‖ – go and g? – go values on going from CsCdF₃:Ag²⁺ to NaF:Ag²⁺ is pointed out to reflect the different electrostatic potential (exerted by the rest of the lattice upon the complex) seen by AgF embedded in NaCl or perovskite-type lattices. © 1994 John Wiley & Sons, Inc.     

Structure and stability of B6, B,and B clusters

We investigated various isomers of B₆, B, and B clusters with ab initio [Hartree–Fock (HF), MP2)] and density functional theory (DFT) methods. Ten B₆ isomers, 6 B isomers, and 6 B isomers are determined to be local minima on their potential energy hypersurfaces by the HF, B3LYP, B3PW91, and MP2 methods. Fourteen of these structures are first reported. The most stable neutral B₆ cluster is the capped pentagonal pyramid (C_5v), in agreement with the results reported previously. Hexagon B (C_2h) isomer and fan‐shaped B (C_2v) isomer are found to be the most stable on the cationic and anionic energy hypersurfaces, respectively. Natural bond orbital analysis suggests that there are three‐centered bonds in the most stable B₆ neutral and ionic clusters. The multicentered bonds are responsible for the special stability of the lowest‐energy isomer. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 269–278, 2003     

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     

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.     

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.     

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.     

Ab initio MO calculations of hyperfine coupling constants. A basis set study for 19F, 35Cl, 19F,and 35Cl

Hyperfine coupling constants (HFCC ) of the ¹⁹F and ³⁵Cl atoms and the ¹⁹F and ³⁵Cl radical anions have been calculated by the unrestricted Hartree–Fock (UHF ) method using polarization and diffuse functions with contracted double-zeta as well as uncontracted basis sets. The A_dip values are fairly insensitive to changes in the basis set and show good accordance with experimental and other theoretical studies. The isotropic HFCCS a_N of ¹⁹F, ¹⁹F, and ³⁵Cl show strong dependence on d functions and the state of contraction of the s, p set. Spin-projected UHF wave functions lead to better agreement with experiment.     

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.     

Explicit expression of the franck–condon factors in terms of the potentials of the two states

The calculus of the overlap integral for two states represented by the vibrational wave functions ψ and ψ is reduced to that of the Franck–Condon integral ?(0, x) = ∫ ψψ (t) dt. It is proved that for “numerical potentials” (as well as for a Dunham potential), this integral is given on each interval by a simple analytic expression in terms of the two potentials. The Franck–Condon factors are well determined by “coupling constants” related uniquely to the coordinates of the turning points of the potentials. An application to the band system BII? XΣ of Nα₂ is compared with the usual numerical methods.