I'm trying to run a Hartree-Fock calculation with an embedding potential. This potential could be applied by modifying the core Hamiltonian, and then doing an scf calculation as normal.
So is there any way that I can run an scf calculation with a modified core Hamiltonian?(preferably a python-side method).
Any help would be appreciated (including just telling this it isn’t reasonably possible).
Thanks in advance