<div dir="ltr">For insulators there is no definition of the Fermi energy. <div>Your system shows a clear gap, indeed it is an insulator.</div><div>In your plot, the zero is most probably set at the energy of the highest occupied state.</div><div>But this choice is arbitrary.</div><div>Regards</div><div>Marcella<br><br>On Saturday, June 9, 2018 at 5:38:57 PM UTC+2, Jannen wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Dear CP2K Users,</div><div><br></div><div>I am calculate the DOS for a crystal structure. It is a insulator. The density of state was found at Fermi energy level as shown in attached picture. I was really confused by the results. The DOS should be 0 at Fermi energy level. Even though the plots was broadening using deskinsgroup code, but that will just cause a slightly shifting above the Fermi energy level. Could anyone help me on the results? Any problem with the input? </div><div><br></div><div>Please see the attached settings for DFT section.</div><div><br></div><div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div>    CHARGE 0</div><div>    MULTIPLICITY 1</div><div>    UKS TRUE</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-14</div><div>    &END QS</div><div>    &MGRID</div><div>      CUTOFF 400</div><div>      NGRIDS 4</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &SCF</div><div>      SCF_GUESS RESTART</div><div>      EPS_SCF 1.0E-07</div><div>      MAX_SCF 1000</div><div>      &OT TRUE</div><div>       PRECONDITIONER FULL_ALL</div><div>       MINIMIZER CG</div><div>       ENERGY_GAP 0.001</div><div>     &END OT</div><div>      &PRINT</div><div>        &RESTART OFF</div><div>        &END RESTART</div><div>      &END PRINT</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>         &VDW_POTENTIAL</div><div>           DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>             &PAIR_POTENTIAL</div><div>              TYPE DFTD3</div><div>              REFERENCE_FUNCTIONAL PBE</div><div>              PARAMETER_FILE_NAME dftd3.dat</div><div>              R_CUTOFF 10.0</div><div>             &END PAIR_POTENTIAL</div><div>          &END VDW_POTENTIAL</div><div>    &END XC</div><div>    &PRINT</div><div>         &PDOS</div><div>            NLUMO -1</div><div>            COMPONENTS TRUE</div><div>            &LDOS</div><div>               COMPONENTS TRUE</div><div>               LIST  1..486    </div><div>            &END LDOS</div><div>         &END</div><div>     &END PRINT<br></div><div>  &END DFT</div></div><div><br></div><div>Best regards,</div></div></blockquote></div></div>