In a previous work Sims and Hagstrom [J Chem Phys 140,224312(2014)] reported
Hylleraas-configuration interaction (Hy-CI) method variational calculations for the 1S
ground states of the beryllium isoelectronic sequence with an estimated accuracy of 10
to 20 nanohartrees (nHa). In this work the calculations have been extended to the five
higher states of the neutral beryllium atom, 3 1S, 4 1S, 5 1S, 6 1S, and 7 1S. The best
non-relativistic energies obtained for these states are -14.4182 4034 6, -14.3700 8789
0, -14.3515 1167 6, -14.3424 0357 8, and -14.3372 6649 96 Ha, respectively. The 6 1S
result is superior to the known reference energy for that state, while for the 7 1S
state there is no other comparable calculation.