Dear Colleagues:
I am trying to perform a geometry optimization within the two-component spin-orbit DFT using the LRC-ωPBEh functional in the 7.7 version. It gets the first energy right but then calculates unreasonable gradients. I do not observe that behavior in the 7.4 version with the same functional or in 7.7 with the regular PBE0 one. It makes me think that something is wrong with the gradient module in the 7.7 version. Could you please help me with that? I'll be happy to provide the necessary inputs.
Thanks,
Alex
