Kohn-Sham orbitals and potentials from quantum Monte Carlo molecular densities

In this work we show the possibility to extract Kohn-Sham orbitals, orbital energies, and exchange correlation potentials from accurate Quantum Monte Carlo (QMC) densities for atoms (He, Be, Ne) and molecules (H2, Be2, H2O, and C2H4). The Variational Monte Carlo (VMC) densities based on accurate Jastrow Antisymmetrised Geminal Power wave functions are calculated through different estimators. Using these reference densities, we extract the Kohn-Sham quantities with the method developed by Zhao, Morrison, and Parr (ZMP) [Phys. Rev. A 50, 2138 (1994)]. We compare these extracted quantities with those obtained form CISD densities and with other data reported in the literature, finding a good agreement between VMC and other high-level quantum chemistry methods. Our results demonstrate the applicability of the ZMP procedure to QMC molecular densities, that can be used for the testing and development of improved functionals and for the implementation of embedding schemes based on QMC and Density Functional Theory.