<div dir="ltr">Dear CP2K users,<div><br></div><div>I would like to perform a DFT+U calculation with LaFeO3 perovskite. However, I have met the following error for KPOINTS setting and I am not sure where is the problem from. Could you kindly help me with this? It would be significantly appreciated for your help.</div><div><br></div><div>Another question is that is it ok to set two BASIS_SET_FILE_NAME within one input file, as I checked that there is no basis set for La in the BASIS_MOLOPT file.</div><div><br></div><div>For your reference, I run this calculation on my laptop with cp2k.popt version 5.1. I also attach the input file here for your information. Thanks.</div><div><img src="cid:f09f3286-fdcb-4a02-90cd-69e752a46e8c" alt="WeChat Screenshot_20191121194541.png" width="320" height="236" style="text-align: center; margin-left: 1em; margin-right: 1em;"></div><div><br></div><div><div>!*************************GLOBAL SETTINGS**********************************</div><div>&GLOBAL</div><div>  PRINT_LEVEL LOW</div><div>  PROJECT_NAME lfo_rlx</div><div>  RUN_TYPE CELL_OPT</div><div>  SAVE_MEM T</div><div>&END GLOBAL</div><div><br></div><div>!*************************MOTION SETTINGS**********************************</div><div>&MOTION</div><div>  &CELL_OPT</div><div>    OPTIMIZER LBFGS</div><div>    TYPE DIRECT_CELL_OPT</div><div>  &END CELL_OPT</div><div>&END MOTION</div><div><br></div><div>!************************FORCE_EVAL SETTINGS********************************</div><div>&FORCE_EVAL</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  METHOD QS</div><div>!********************FORCE_EVAL -> DFT SETTINGS*****************************</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT_UCL</div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div>    UKS T</div><div>    CHARGE 0</div><div>    PLUS_U_METHOD MULLIKEN</div><div>!****************FORCE_EVAL -> DFT -> SCF SETTINGS**************************</div><div>    &SCF</div><div>      ADDED_MOS 100</div><div>      MAX_SCF 200</div><div>      EPS_SCF 1.0E-6</div><div>      SCF_GUESS RESTART</div><div>      &DIAGONALIZATION T</div><div>        ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &OUTER_SCF</div><div>        EPS_SCF 1.0E-5</div><div>        MAX_SCF 5</div><div>      &END OUTER_SCF</div><div>      &SMEAR T</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE [K] 300</div><div>      &END SMEAR</div><div>      &MIXING T</div><div>        METHOD PULAY_MIXING</div><div>        ALPHA 0.1</div><div>        BETA [bohr^-1] 0.001</div><div>        NMIXING 5</div><div>        NBUFFER 8</div><div>      &END MIXING</div><div>    &END SCF</div><div>!****************FORCE_EVAL -> DFT -> QS SETTINGS***************************</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-12</div><div>      EXTRAPOLATION USE_GUESS</div><div>      METHOD GAPW</div><div>    &END QS</div><div>!***************FORCE_EVAL -> DFT -> GRID SETTINGS**************************</div><div>    &MGRID</div><div>      NGRIDS 5</div><div>      CUTOFF [Ry] 600</div><div>      REL_CUTOFF [Ry] 60</div><div>    &END MGRID</div><div>!*****************FORCE_EVAL -> DFT -> XC SETTINGS**************************</div><div>    &XC</div><div>      FUNCTIONAL_ROUTINE NEW</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>!***************FORCE_EVAL -> DFT -> POISSON SETTINGS************************</div><div>    &POISSON</div><div>      POISSON_SOLVER PERIODIC</div><div>      PERIODIC XYZ</div><div>    &END POISSON</div><div>!***************FORCE_EVAL -> DFT -> KPOINTS SETTINGS************************</div><div>    &KPOINTS</div><div>      SCHEME  MONKHORST-PACK  2  2  2</div><div>      FULL_GRID  .TRUE.</div><div>    &END KPOINTS</div><div>  &END DFT</div><div>!********************FORCE_EVAL -> SUBSYS SETTINGS***************************</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC [angstrom] 7.91435 7.91435 7.91435</div><div>      ALPHA_BETA_GAMMA [deg] 90 90 90</div><div>      PERIODIC XYZ</div><div>    &END CELL</div><div>    &KIND La</div><div>      ELEMENT La</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q11</div><div>    &END KIND</div><div>    &KIND Fe</div><div>      ELEMENT Fe</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q16</div><div>      MAGNETIZATION 2.0</div><div>      &DFT_PLUS_U T</div><div>        L 2</div><div>        U_MINUS_J [eV] 5.3</div><div>      &END DFT_PLUS_U</div><div>    &END KIND</div><div>    &KIND O</div><div>      ELEMENT O</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME lfo_222.xyz</div><div>      COORD_FILE_FORMAT XYZ</div><div>    &END TOPOLOGY</div><div>  &END SUBSYS</div><div><br></div><div>!********************FORCE_EVAL -> PRINT SETTINGS****************************</div><div>  &PRINT</div><div>    &FORCES ON</div><div>    &END FORCES</div><div>  &END PRINT</div><div>&END FORCE_EVAL</div></div><div><br></div><div><br></div></div>