The ZPL energy of internal optical transition in defect is estimated to be the difference between the excited state (calculated by constraining the excited state electronic configuration during relaxation) total energy and the ground state total energy. However, the constrain at the excited state can be done with LDIAG=.FALSE. or without it. LDIAG=.FALSE. do not allow rotation of orbitals during electronic optimization process and vice versa. Which is more accurate for calculating ZPL?