Assume I have the electron density represented as the expansion coefficients under an DF auxiliary basis, which come from machine learning prediction. Although the kinetic energy cannot be evaluated without an AO basis, theoretically the Kohn-Sham effective potential has no kinetic term and should be able to be constructed from DF density (assume I am using PBE). Then I might restart the DFT from this effective potential.

Is this possible in psi4? And if so, may I have some hints on how to do this?