Hylleraas-Correlated Complex-Scaling Variational Calculations of Ground-State Energies in Four-Electron Beryllium-Like Systems

We present a variational study of the ground-state energies of the four-electron members of the beryllium isoelectronic sequence (Be, B⁺, C²⁺) using explicitly correlated Hylleraas-type trial wave functions formulated within a complex-scaling (complex-rotation) framework. The trial function incorporates all six interelectronic distances of the four-electron problem together with a restricted set of nonlinear variational parameters, and the nonrelativistic Hamiltonian is evaluated after the coordinate transformation $r\to r{\text{e}}^{\left\{\text{i}\theta \right\}}$ . Analytic matrix elements were derived and evaluated numerically with the computer-algebra system Maxima. For the (1s²2s²) ¹S ground state, the resulting energies are compared with Hartree-Fock and variational Monte Carlo values, and with high-accuracy Hylleraas-configuration-interaction (Hy-CI) benchmark energies available for the Be-like sequence. The present energy for neutral beryllium agrees with the Hy-CI reference to within 8 × 10^？5 a.u., and the B⁺ energy agrees to within 1.1 × 10^？3 a.u. For C²⁺, however, a deviation of about 0.41 a.u. (≈11 eV) is found, revealing a clear limitation of the present implementation for the more highly charged member of the sequence. The evolution of the optimized nonlinear parameters along the sequence is examined and is found to coincide with this breakdown, without our data allowing any causal conclusion to be drawn. The complex-scaling formalism used to construct the trial function is also discussed; the nonzero imaginary parts obtained for the nominally bound ground states are not physically defensible as resonance widths and are therefore not analyzed here as a physical result. The results show close numerical agreement with benchmark values for the lower members of the sequence, while also revealing implementation limitations that prevent the present calculations from being regarded as fully validated variational ground-state energies. These findings identify the numerical and formal improvements required before the present framework can be extended reliably to genuine four-electron autoionizing resonances.