def METHOD_NAME(self, cell=None, dm=None, dm_last=0, vhf_last=0, hermi=1,
             kpt=None, kpts_band=None):
    vj, vk = self.get_jk(cell, dm, hermi, kpt, kpts_band, True, True)
    vhf = vj - vk
    return vhf