<div dir="ltr"><div>Dear All,</div><div><br></div><div>I'm trying to optimize a solid state structure in a periodic model using PBE level of DFT including dispersion; this system has about 180 atoms.</div><div>when I try to optimize the system, it crashes immediately and I come up with an error saying "On entry to DLASCL parameter number  4 had an illegal value". </div><div><br></div><div>I would appreciate if someone could have a look and comment on it.</div><div><br></div><div>Thank you in advance.</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>Here is the input file:</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT TEST</div><div>  PRINT_LEVEL MEDIUM</div><div>  RUN_TYPE GEO_OPT</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME /home/programs/cp2k/cp2k-2.5.<wbr>1/tests/QS/BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME /home/programs/cp2k/cp2k-2.5.<wbr>1/tests/QS/GTH_POTENTIALS</div><div>    CHARGE 0</div><div>    MULTIPLICITY 1</div><div>  &MGRID</div><div>      NGRIDS 4</div><div>      CUTOFF 500</div><div>      REL_CUTOFF 40</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-10</div><div>      WF_INTERPOLATION PS</div><div>      EXTRAPOLATION_ORDER 3</div><div>      METHOD GPW</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 1.0E-7</div><div>       &OT ON</div><div>         MINIMIZER DIIS</div><div>       &END OT</div><div>      MAX_SCF 30</div><div>      &OUTER_SCF</div><div>        MAX_SCF 30</div><div>        EPS_SCF 1.0E-7</div><div>      &END OUTER_SCF</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div></div><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>            PARAMETER_FILE_NAME /home/programs/cp2k/cp2k-2.5.<wbr>1/tests/QS/dftd3.dat</div><div>            REFERENCE_FUNCTIONAL PBE</div><div>         &END PAIR_POTENTIAL</div><div>      &END vdW_POTENTIAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 11.4081 11.6597 16.1620</div><div>      ALPHA_BETA_GAMMA 80.452 72.505 79.100</div><div>      PERIODIC XYZ</div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>      COORDINATE XYZ</div><div>      COORD_FILE_NAME <a href="http://input.xyz" target="_blank">input.xyz</a></div><div>      CONNECTIVITY OFF</div><div>    &END TOPOLOGY</div><div><br></div><div>    &KIND Cu</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q11</div><div>    &END KIND</div><div>    &KIND C</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND</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 B</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q3</div><div>    &END KIND</div><div>    &KIND F</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q7</div><div>    &END KIND</div><div>    &KIND H</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q1</div><div>    &END</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&MOTION</div><div>  &GEO_OPT</div><div>    OPTIMIZER BFGS</div><div>    MAX_FORCE 1.0E-4</div><div>  &END GEO_OPT</div><div>  &PRINT</div><div>     &TRAJECTORY</div><div>      &EACH</div><div>         GEO_OPT 1</div><div>      &END EACH</div><div>    &END TRAJECTORY</div><div>  &END</div><div>&END MOTION</div></div><div><br></div></div>