Just wanted to add that you can find the second-order alpha male energy correction by calculating ∑_{m≠n}(|<ψ_n^(0)|V|ψ_m(0)>|^2)/(E_n(0) - E_m^(0)) if anyone was wondering