Reduced first order density matrix for the Be ground state

A reduced first order density matrix for the Be ground state is computed from an extensive configuration interaction (CI) wave function. A sequence of increasingly accurate CI wave functions Φq converging towards the exact Ψ is used to assess the quality of the results which include approximate bounds for the overlaps 〈Φq|Ψ〉, electron–nuclear coalescence cusp data, Weinhold's overlap between density matrices, virial ratios, occupation number spectra, and some expectation values. The nuclear magnetic shielding constant and the molar diamagnetic susceptibility are determined with 2.0 and 1.5% of uncertainty, respectively.