Thank you for your suggestion. However, what I am referring to is not the charge defect calculation. Basically, I believe charge-point defects result from adding/removing some electrons to the structure. I knew this feature QuantumATK works very well. But I would like to manually occupy the kohn-sham levels. For example, I would like to excite an electron from band#100 to band#101 (assume the ground state is a closed shell singlet). How can I implement it?