A procedure is described for the precise nonrelativistic evaluation of the dipole polarizabilities of H2 + and D2 + that avoids any approximation based on the size of the electron mass relative to the nucleus mass. The procedure is constructed so that sum rules may be used to assess the accuracy of the calculation. The resulting polarizabilities are consistent with experiment within the error bars of the measurements and are far more precise than values obtained by other theoretical methods.PACS numbers: 33.15. Kr, 31.15.Ar The separation of nuclear and electronic motion is the underlying principle of the theory of molecular structure. The theory is challenged by recent measurements of Jacobson et al.[1] of the electric dipole polarizabilities of H 2 + and D 2 + which have a precision beyond that obtained in the Born-Oppenheimer approximation. The measurements stimulated the introduction of methods [2][3][4][5] that take into account the diabatic coupling omitted in the earlier calculations and they led to polarizabilities that agree with the measured values within the combined experimental and theoretical uncertainties. We present here new theoretical predictions of much greater accuracy which in turn pose a significant challenge to experiment. The accuracy of our method can be assessed by the use of sum rules and we predict nonrelativistically the polarizabilities of H 2 + and D 2 + to a precision well beyond that achieved by the experiments. The method is general and it should be possible to apply it to manyelectron diatomic molecules.Separating out the center of mass motion we may write for the Hamiltonian of H 2 + or D 2 + in an electric field F = Fn lying along the Z-axis of the space-fixed framewhere R is the vector joining the nuclei, r is the position vector of the electron measured from the midpoint of R, M is the mass of the proton or deuteron, V (r, R) is the electrostatic interaction potential and (1 + ǫ) = [1 + (1 + 2M ) −1 ]. We use atomic units throughout. The change in energy of the system for small values of the applied field is given by ∆E = − 1 2 α d F 2 , where α d is the polarizability. Thus if Ψ (0) (r, R) is the eigenfunction of the unperturbed system with Hamiltonian H 0 and E 0 is the eigenvalue, the polarizability can be writtenwhereAlternatively Ψ (1) can be determined from the stationary value of the functionalIf we write Ψ (1) (r, R) as an expansion over some chosen basis set ψ n (r, R),assumed to diagonalize the unperturbed Hamiltonian H 0 so that ψ n |H 0 |ψ n ′ = E n δ nn ′ , the polarizability may be writtenThis expression for the polarizability is stationary with respect to first order errors in Ψ (1) and is bounded from below.The completeness of the set ψ n (r, R) can be assessed by inspecting other sum rules. Introduce the oscillator strengthand define the sumso thatThen provided the ψ n form a complete set,and S(0) = 1.The eigenfunctions Ψ (0) (r, R) and ψ n (r, R) can be written as sums of products of nuclear and electronic wave functions of the form 1