Non-Born-Oppenheimer electronic and nuclear densities for a Hooke-Calogero three-particle model: non-uniqueness of density-derived molecular structure.

We consider the calculation of non-Born-Oppenheimer, nBO, one-particle densities for both electrons and nuclei. We show that the nBO one-particle densities evaluated in terms of translationally invariant coordinates are independent of the wavefunction describing the motion of center of mass of the whole system. We show that they depend, however, on an arbitrary reference point from which the positions of the vectors labeling the particles are determined. We examine the effect that this arbitrary choice has on the topology of the one-particle density by selecting the Hooke-Calogero model of a three-body system for which expressions for the one-particle densities can be readily obtained in analytic form. We extend this analysis to the one-particle densities obtained from full Coulomb interaction wavefunctions for three-body systems. We conclude, in view of the fact that there is a close link between the choice of the reference point and the topology of one-particle densities that the molecular structure inferred from the topology of these densities is not unique. We analyze the behavior of one-particle densities for the Hooke-Calogero Born-Oppenheimer, BO, wavefunction and show that topological transitions are also present in this case for a particular mass value of the light particles even though in the BO regime the nuclear masses are infinite. In this vein, we argue that the change in topology caused by variation of the mass ratio between light and heavy particles does not constitute a true indication in the nBO regime of the emergence of molecular structure.

Analytical method for the representation of atoms-in-molecules densities.

We present analytic refinements and applications of the deformed atomic densities method [Fernández Rico, J.; López, R.; Ramírez, G. J Chem Phys 1999, 110, 4213-4220]. In this method the molecular electron density is partitioned into atomic contributions, using a minimal deformation criterion for every two-center distributions, and the atomic contributions are expanded in spherical harmonics times radial factors. Recurrence relations are introduced for the partition of the two-center distributions, and the final radial factors are expressed in terms of exponential functions multiplied by polynomials. Algorithms for the practical implementation are developed and tested, showing excellent performances. The usefulness of the present approach is illustrated by examining its ability to describe the deformation of atoms in different molecular environments and the relationship between these atomic densities and some chemical properties of molecules.

A finite B-spline basis set for accurate diatomic molecule calculations.

A finite basis set particularly adapted for solving the Hartree-Fock equation for diatomic molecules in prolate spheroidal coordinates has been constructed. These basis functions have been devised as products of B-splines times associated Legendre polynomials. Due to the large number of B-splines, the resulting set of eigenfunctions is amply distributed over excited states. This gives the possibility of using these basis sets to calculate sums over excited states, appearing in various orders of perturbation theory. As an illustration, the second-order corrections to the ground-state energy of some atoms and diatomic molecules with closed electron shells have been calculated.