<div dir="ltr">Dear CP2K experts:<div>When I use CP2K to perform an AIMD calculation, I want to print the HOMO-LUMO GAP during MD calculation. However,the program tends to stuck when output hte EIGENVALUE, and I don't kown how th fix it ;<br></div><div><div><br></div><div><br></div><div>&GLOBAL</div><div>  PROJECT ${PROJECT_NAME}</div><div>  RUN_TYPE MD</div><div>  PRINT_LEVEL LOW </div><div>&END GLOBAL</div><div><br></div><div>&MOTION</div><div>  &MD</div><div>    ENSEMBLE NVT</div><div>    STEPS 10000</div><div>    TIMESTEP 1   </div><div>    TEMPERATURE 300</div><div>    &THERMOSTAT</div><div>      REGION MASSIVE</div><div>      &NOSE </div><div>        TIMECON 1000 </div><div>      &END NOSE</div><div>    &END</div><div>  &END MD</div><div>  &GEO_OPT</div><div>    OPTIMIZER BFGS </div><div>    MAX_FORCE 2E-4</div><div>    MAX_ITER 500</div><div>  &END GEO_OPT</div><div>  &PRINT</div><div>    &TRAJECTORY</div><div>      &EACH</div><div>        MD 1</div><div>      &END EACH</div><div>    &END TRAJECTORY</div><div>    &VELOCITIES OFF</div><div>    &END VELOCITIES</div><div>    &FORCES OFF</div><div>    &END FORCES</div><div>    &RESTART_HISTORY</div><div>      &EACH</div><div>        MD 500</div><div>      &END EACH                                                                                                                    </div><div>    &END RESTART_HISTORY                              </div><div>    &RESTART</div><div>    BACKUP_COPIES 3                </div><div>      &EACH     </div><div>        MD 1                   </div><div>      &END EACH   </div><div>    &END RESTART      </div><div>  &END PRINT     </div><div>&END MOTION</div><div>  #STRESS_TENSOR  ANALYTICAL #CELL_OPT should be opened</div><div>&FORCE_EVAL  </div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME ${data}/BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME ${data}/GTH_POTENTIALS</div><div>    CHARGE 0</div><div>    MULTIPLICITY 1</div><div>    &MGRID</div><div>      CUTOFF 500</div><div>      NGRIDS 5</div><div>      REL_CUTOFF 50</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-6</div><div>      METHOD GPW</div><div>      EXTRAPOLATION ASPC</div><div>    &END QS</div><div>    &SCF</div><div>      MAX_SCF 30</div><div>      #SCF_GUESS RESTART</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 1.0E-5</div><div>      CHOLESKY INVERSE_DBCSR ##</div><div>      &OUTER_SCF</div><div>        MAX_SCF 40</div><div>        EPS_SCF 1.0E-5</div><div>      &END</div><div>      &OT </div><div>        MINIMIZER DIIS</div><div>        PRECONDITIONER  FULL_SINGLE_INVERSE</div><div>        STEPSIZE 0.01</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>          POTENTIAL_TYPE PAIR_POTENTIAL</div><div>          &PAIR_POTENTIAL</div><div>            PARAMETER_FILE_NAME ${data}/dftd3.dat    </div><div>            TYPE DFTD3</div><div>            REFERENCE_FUNCTIONAL PBE</div><div>            R_CUTOFF [angstrom] 16   </div><div>          &END PAIR_POTENTIAL </div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div>    &PRINT</div><div>      &E_DENSITY_CUBE OFF</div><div>      &END E_DENSITY_CUBE</div><div>      &MO_CUBES</div><div>        NLUMO 2</div><div>        NHOMO 2</div><div>        WRITE_CUBE .FALSE.</div><div>        &EACH</div><div>          MD 10</div><div>        &END EACH</div><div>      &END MO_CUBES</div><div>    &END PRINT  </div><div>  &END DFT</div><div>  # &MM</div><div>  #  &FORCEFIELD</div><div>  #    &SPLINE</div><div>  #      EMAX_SPLINE 5E-0.01</div><div>  #    &END SPLINE</div><div>  #  &END FORCEFIELD</div><div>  #&END MM</div><div><br></div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 17.37820 17.45740 50.00000 </div><div>      PERIODIC XYZ</div><div>    &END CELL</div><div><br></div><div>    &TOPOLOGY</div><div>      COORD_FILE_FORMAT XYZ</div><div>      COORD_FILE_NAME configure.xyz</div><div>    &END</div><div><br></div><div>   &KIND C </div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>   &END KIND</div><div><br></div><div>   &KIND N</div><div>     BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>     POTENTIAL GTH-PBE-q5</div><div>   &END KIND </div><div>     </div><div>   &KIND H</div><div>     BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>     POTENTIAL GTH-PBE-q1</div><div>   &END KIND</div><div>     </div><div>   &KIND I      </div><div>     BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>     POTENTIAL GTH-PBE-q7</div><div>   &END KIND</div><div><br></div><div>   &KIND Pb</div><div>     BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>     POTENTIAL GTH-PBE-q4</div><div>   &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div></div></div>