Ab initio wavefunction analysis of electron removal quasi-particle state of NdNiO$_2$ with fully correlated quantum chemical methods

The discovery of superconductivity in hole-doped infinite-layer NdNiO$_2$ -- a transition metal (TM) oxide that is both isostructural and isoelectronic to cuprate superconductors -- has lead to renewed enthusiasm in the hope of understanding the origin of unconventional superconductivity. Here, we investigate the electron-removal states in infinite-layered Ni$^{1+}$ oxide, NdNiO$_2$, which mimics hole-doping, with the state-of-the-art many-body multireference quantum chemistry methods. From the analysis of the many-body wavefunction, we find that the hole-doped $d^8$ ground state of NdNiO$_2$ is very different from the $d^8$ ground state in isostructural cuprate analog CaCuO$_2$, although the parent $d^9$ ground states are for the most part identical. We show that the doped hole in NdNiO$_2$ mainly localizes on the Ni $3d_{x^2-y^2}$ orbital to form a closed-shell singlet, and this singlet configuration contributes to $\sim$40% of the wavefunction. In contrast, in CaCuO$_2$ the Zhang-Rice singlet configurations contribute to $\sim$65% of the wavefunction. With the help of the quantum information concept of entanglement entropy, we quantify the different types of electronic correlations in the nickelate and cuprate compounds and find that the dynamic radial-type correlations within the Ni $d$ manifold are persistent in hole-doped NdNiO$_2$. As a result, the $d^8$ multiplet effects are stronger and the additional hole foot-print is more three-dimensional in NdNiO$_2$. Our analysis shows that the most commonly used three-band Hubbard model employed to express the doped scenario in cuprates represents $\sim$90% of the $d^8$ wavefunction for CaCuO$_2$, but such a model grossly approximates the $d^8$ wavefunction for NdNiO$_2$ as it only stands for $\sim$60% of the wavefunction.