<div dir="ltr"> <br>Dear all<br><br>HI. I am trying to calculate Band-gap of ZnO by means of Hybrid-DFT (PBE and HF) with reference to this example file (https://www.cp2k.org/exercises:2014_ethz_mmm:tio2_gap)<br>But, an Error message appeared immediately after starting this calculation. the message is as follows. What do I have to do to solve this problem? I guess this question is very stupid but I have no idea. <br><br><br>**********************<br> *** WARNING in hfx_energy_potential.F:626 :: The Kohn Sham matrix is not  ***<br> *** 100% occupied. This may result in uncorrect Hartree-Fock results. Try ***<br> *** to decrease EPS_PGF_ORB and EPS_FILTER_MATRIX in the QS section.      ***<br>***********************<br><br>****input****<br><br>&GLOBAL<br>   PROJECT_NAME zno<br>   RUN_TYPE  CELL_OPT<br> &END GLOBAL<br> &MOTION<br>   &GEO_OPT<br>     OPTIMIZER BFGS<br>     MAX_ITER  500<br>     MAX_DR     1E-6<br>     MAX_FORCE     1E-06<br>   &END GEO_OPT<br>   &CELL_OPT<br>     OPTIMIZER BFGS<br>     MAX_ITER  500<br>     MAX_DR     1E-6<br>     MAX_FORCE     1E-06<br>   &END CELL_OPT<br> &END MOTION<br> &FORCE_EVAL<br>   STRESS_TENSOR ANALYTICAL<br>   METHOD  QS<br>   &DFT<br>    BASIS_SET_FILE_NAME /cp2k-3.0/data/BASIS_MOLOPT<br>    POTENTIAL_FILE_NAME /cp2k-3.0/data/POTENTIAL<br>    &SCF<br>      SCF_GUESS  RESTART<br>      EPS_SCF 1.0E-8<br>      MAX_SCF 10<br>      &OUTER_SCF<br>        EPS_SCF 1.0E-8<br>        MAX_SCF 1000<br>      &END OUTER_SCF<br>      &OT ON<br>        PRECONDITIONER FULL_SINGLE_INVERSE<br>        MINIMIZER DIIS<br>      &END OT<br>    &END SCF<br>    &QS<br>      METHOD  GPW<br>    &END QS<br>    &MGRID<br>      CUTOFF    500<br>    &END MGRID<br><br>  &PRINT              ! section required to obtain the HOMO-LUMO gap.<br>     &MO_CUBES<br>       WRITE_CUBE .FALSE. ! no cube file is generated<br>       NHOMO 1            ! but we require 1 HOMO and 1 LUMO  in the output<br>       NLUMO 1            ! so that we get the band gap<br>     &END<br>  &END<br>    &XC                      ! this is the section to define the electronic exchange<br>      &XC_FUNCTIONAL          ! our functional is hybrid<br><br>         &PBE                 ! it has 75% of PBE<br>           SCALE_X 0.75<br>           SCALE_C 1.0<br>         &END<br>         &PBE_HOLE_T_C_LR<br>           SCALE_X 0.25       ! + 25% of truncated PBE0 functional - that includes exact hfx<br>           CUTOFF_RADIUS 3.5  ! that has interaction truncated at 3.5 A from the atomic core<br>         &END<br>      &END XC_FUNCTIONAL<br>      &HF<br>        FRACTION 0.25         ! this is the hfx section. The amount of hfx must be consistent with above<br>        &SCREENING            ! Screening of the electronic repulsion up to the given threshold. This section is needed<br>          EPS_SCHWARZ 1.0E-6<br>          SCREEN_ON_INITIAL_P TRUE  ! having an external wave-function, a preliminary screening<br>        &END                        ! can be performed to speed up calculations<br>        &INTERACTION_POTENTIAL      ! Sets up interaction potential between the two regions<br>          POTENTIAL_TYPE TRUNCATED  ! the potential is  truncated<br>          CUTOFF_RADIUS 3.5         ! at 3.5 A (see above)<br>#          T_C_G_DATA ./t_c_g.dat    ! external file with parameters needed to truncate the potential<br>        &END<br>      &END<br>     &END<br> &END DFT<br>   &SUBSYS<br><br>&CELL<br>  PERIODIC XYZ<br>  A   2.814062947057155  -1.624700000000000   0.000000000000000<br>  B   0.000000000000000   3.249400000000000   0.000000000000000<br>  C   0.000000000000000   0.000000000000000   5.203800000000000<br>       MULTIPLE_UNIT_CELL  3 3 3    ! a 2x2x1 system is required to get realistic results<br>&END CELL<br><br>&COORD<br>Zn  0.938020982352385   1.624700000000000   0.000000000000000<br>Zn  1.876041964704770   0.000000000000000   2.601900000000000<br>O  0.938020982352385   1.624700000000000   1.988371980000000<br>O  1.876041964704770   0.000000000000000   4.590271980000001<br>&END COORD<br><br>&TOPOLOGY<br>       MULTIPLE_UNIT_CELL  3 3 3<br>&END TOPOLOGY<br><br>&TOPOLOGY<br>       MULTIPLE_UNIT_CELL  3 3 3<br>&END TOPOLOGY<br><br><br>@set FUNCTIONAL PBE<br>     &KIND Zn<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       POTENTIAL GTH-${FUNCTIONAL}-q12<br>     &END KIND<br>     &KIND O<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       POTENTIAL GTH-${FUNCTIONAL}-q6<br>     &END KIND<br>   &END SUBSYS<br> &END FORCE_EVAL<br>                               <br>****************************************************<br></div>