Wavefunction and energy of the 1s²2sns configuration in a beryllium atom

A new set of trial functions for 1s^22sns configurations in a beryllium atom is suggested. A Mathematica program based on the variational method is developed to calculate the wavefunctions and energies of 1s^22sns (n=3-6) configurations in a beryllium atom. Non-relativistic energy, polarization correction and relativistic correction which include mass correction, one- and two-body Darwin corrections, spin-spin contact interaction and orbit-orbit interaction, are calculated respectively. The results are in good agreement with experimental data.