Novel theory of the HD dipole moment. II. Computations

Charge asymmetry in the rovibrationally excited HD molecule

The recently developed method for performing all-particle non-Born-Oppenheimer variational calculations on diatomic molecular systems excited to the first excited rotational state and simultaneously vibrationally excited is employed to study the charge asymmetry and the level lifetimes of the HD molecule. The method uses all-particle explicitly correlated Gaussian functions. The nonlinear parameters of the Gaussians are optimized with the aid of the analytical energy gradient determined with respect to these parameters.

Classical calculations of radiative association in absence of electronic transitions

A formula for the cross section of radiative association where no electronic transitions take place is derived and tested for diatomic molecules. The approach is based on classical mechanics and therefore it is valid for direct, i.e., non-resonant, radiative association. For the formation of carbon monoxide (CO) and the cyano radical (CN), in the X1Σ+ and A1Π states, respectively, the treatment reproduces the baselines of the cross sections obtained using quantum mechanical perturbation theory. The method overestimates the formation cross section of potassium sodide (NaK) by about 8%. For the lower mass diatoms hydrogen fluoride (HF) and deuterium hydride (HD), the formula overestimates the cross sections by 12% and 60%, respectively. The formula can be used alone for estimates of radiative association rate constants, or in combination with Breit-Wigner theory to include resonance contributions.

Charge asymmetry in HD+

Expanding the wave functions of the ground and excited states of HD(+) (or pde) in terms of spherically symmetric explicitly correlated Gaussian functions with preexponential multipliers consisting of powers of the internuclear distance, and using the variational method, we performed very accurate nonadiabatic calculations of all bound states of this system corresponding to the zero total angular momentum quantum number (vibrational states; v=0-22). The total and the transition energies obtained agree with the best available calculations. For each state we computed the expectation values of the d-p, d-e, and p-e interparticle distances. This is the first time these quantities were computed for HD(+) using rigorous nonadiabatic wave functions. While up to the v=20 state some asymmetry is showing in the d-e and p-e distances, for v=21 and v=22 we observe a complete breakdown of the Born-Oppenheimer approximation and localization of the electron almost entirely at the deuteron.

A new algorithm to handle finite nuclear mass effects in electronic calculations: The ISOTOPE program

We report the development of a simple algorithm to modify quantum chemistry codes based on the LCAO procedure, to account for the isotope problem in electronic structure calculations. No extra computations are required compared to standard Born-Oppenheimer calculations. An upgrade of the Gamess package called ISOTOPE is presented, and its applicability is demonstrated in some examples.

Non-Born-Oppenheimer isotope effects on the polarizabilities of H2

We present non-Born-Oppenheimer quantum-mechanical calculations of the behavior of isolated molecules of the H2 isotopomer series in static electric fields. Some conceptual aspects of such calculations are discussed. The values for polarizabilities and hyperpolarizabilities of the H2 isotopomers which we present are the first ever fully nonadiabatic calculated values of these properties.