Some low-lying states of the nine-valence-electron systems HBeP À and HPBe À anions have been studied for the¯rst time using three methods CASSCF, CASPT2 and B3LYP with the contracted atomic natural orbital (ANO) and cc-pVTZ basis sets. The geometries of all stationary points along the potential energy surfaces were optimized at the CASSCF/ANO, CASPT2/ ANO and B3LYP/cc-pVTZ levels. The potential energy curves of isomerization reactions between HBeP À and HPBe À were calculated as a function of HBeP bond angle. The ground and the¯rst excited states of HBeP À are predicted to be X 2 Å and A 2 AE þ states, respectively. The X 2 AE þ and A 2 Å states of the linear HPBe À are both¯rst-order saddle points because they have unique imaginary frequency. Two bent minima M1 and M2 were found along the 1 2 A 0 and 1 2 A 00 potential energy surfaces, respectively. The calculated results indicated that the ground-state HBeP À is linear, while the ground-state HPBe À is bent.