eng Science and Education Publishing International Journal of Physics 2333-4576 2026-10-21 14 4 86 92 10.12691/ijp-14-4-2 IJP20261442 article Boubacar Sow boubacar1.sow@ucad.edu.sn 1 Babou Diop 1 Malick Sow 1 Omar Absatou Niass 1 Department of Physics, Atoms Laser Laboratory, University Cheikh Anta Diop, Dakar, Senegal Total energies of three-electron atomic systems along the lithium isoelectronic sequence are computed with the variational Monte Carlo (VMC) method using compact, fully correlated Slater-Jastrow trial wave functions. The calculation addresses the (1s2 2s) 2Se ground state, the (1s2 2p) 2Po first excited state and the metastable doubly excited (1s2s2p) 4Po quartet state, and revisits, with an independent stochastic method, the states studied previously with a hypergeometric correlation factor. The trial wave functions combine a double-zeta 1s orbital, Slater-type 2s and 2p orbitals and a Padé-Jastrow correlation factor satisfying the electron-electron cusp conditions exactly; all parameters are optimised by a mixed energy-variance criterion on correlated samples. The VMC energies reproduce the full core plus correlation (FCPC) reference values to within 0.07% or better over the whole sequence, the virial ratio remains within |¦Ç - 2| ¡Ü 0.023 for all twenty-four states, and the 2s ¡ú 2p resonance excitation energies agree with the reference differences to about 0.1 eV or better. For the quartet state, the computed excitation energy of the lithium atom, 57.54 eV, reproduces the experimental term value 57.41 eV to 0.13 eV, and a Z-expansion fit of the sequence yields a leading coefficient 0.75022, in excellent agreement with the exact hydrogenic limit 3/4. The variational upper-bound property of the VMC energies is used to examine critically some previously published values for the 2Po series, and the specific difficulties raised by autoionizing doublet states such as (1sns2) 2Se, (1s2sns) 2Se and (1s2snp) 2Po are discussed in detail. https://pubs.sciepub.com/ijp/14/4/2/ijp-14-4-2.pdf Variational Monte Carlo Correlated wave functions Lithium isoelectronic sequence Doubly excited states Jastrow factor Cusp conditions