Thank you for the kind reply.
I could gain Va matrix successfully via upgrading the psi4 version as you told.
Actually, I'm trying to reproduce the DFT energies with matrices.
For the RHF, I could gain HF energy as
Etot=(H+F)*D+E_nuc (where Etot = total energy, H is core-Hamiltonian, T+Vne and E_nuc is nuclear-nuclear repulsion one.)
But in the case of DFT, I couldn't get the same value with that.
I saw the tutorial you've mentioned but it is kind of hard to me since I'm lack of information to know what the psi4 module is doing such as Vpot.properties(), Vpot.functional(), superfunctional, etc.
My goal is calculating DFT like the psi4numpy tutorial 3a.(https://github.com/psi4/psi4numpy/blob/master/Tutorials/03_Hartree-Fock/3a_restricted-hartree-fock.ipynb)
And, therefore, I'm trying to modify that tutorial 3a to do the DFT.
I've found ks_helper.py in the Tutorial 04 but couldn't use it for the same reason above...
So, my main question is if it is possible to create Exc directly with DFT Fock matrix (=H+2J+Vxc).
Thank you very much.