Contraint electron in one of the degenerate orbitals

[baranas]

Hello, everyone

I am performing DFT Delta SCF excited state calculations for crystal defect system. In the ground state I see two non-occupied degenerate orbitals. I would like to perform calculations exciting one electron to one of those degenerate states. However I am stuck in an infinite SCF loop, because electron keeps hopping from one degenerate orbital to another. I have managed to relax such system, by breaking the initial symmetry. However I want to get forces and relaxation energy when geometry is kept in the ground state. Distributing electron evenly to two of those orbitals gives me additional interaction between fractional orbitals and suppress asymmetric forces. Is there any way to make electron constrained in one of those orbitals?

[support_vasp]

Hi,

We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.

Best wishes,

VASP