Ionization potential ofBe9calculated including nuclear motion and relativistic corrections

Abstract

Variational calculations employing explicitly correlated Gaussian functions have been performed for the ground states of $^{9}\mathrm{Be}$ and $^{9}\mathrm{Be}^{+}$ including the nuclear motion [i.e., without assuming the Born-Oppenheimer (BO) approximation]. An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed, leading to energies of the two systems noticeably improved over those found in the recent paper of Pachucki and Komasa [Phys. Rev. A 73, 052502 (2006)]. The non-BO wave functions were used to calculate the ${\ensuremath{\alpha}}^{2}$ relativistic corrections $(\ensuremath{\alpha}={e}^{2}∕\ensuremath{\hbar}c)$. With those corrections and the ${\ensuremath{\alpha}}^{3}$ and ${\ensuremath{\alpha}}^{4}$ corrections taken from Pachucki and Komasa, a new value of the ionization potential (IP) of $^{9}\mathrm{Be}$ was determined. It agrees very well with the most recent experimental IP.