<div dir="ltr">Dear CP2K users,<div><br></div><div>Aside for my previous post, I am also interested in using CP2K to model protein crystals. Some of the proteins I am interested in are metalloproteins. Using again my test system, beta glycine, I am simply adding the following block to an input that works to preform a CELL_OPT. Please find part of this input below:</div><div><br></div><div> &FORCE_EVAL<br></div><div>   METHOD  QS</div><div>   STRESS_TENSOR  ANALYTICAL</div><div>   &DFT</div><div>     BASIS_SET_FILE_NAME /usr/local/src/cp2k/data/GTH_BASIS_SETS</div><div>     POTENTIAL_FILE_NAME /usr/local/src/cp2k/data/GTH_POTENTIALS</div><div>     MULTIPLICITY  0</div><div>     CHARGE  0</div><div>     &SCF</div><div>       MAX_SCF  20</div><div>       EPS_SCF  1.e-7</div><div>       SCF_GUESS  ATOMIC</div><div>       &OT  T</div><div>         MINIMIZER CG</div><div>         PRECONDITIONER FULL_SINGLE_INVERSE</div><div>       &END OT</div><div>       &OUTER_SCF  T</div><div>         EPS_SCF  1.e-6</div><div>         MAX_SCF  50</div><div>       &END OUTER_SCF</div><div>     &END SCF</div><div>     &QS</div><div>       LS_SCF T</div><div>       EPS_DEFAULT     1.0000000000000000E-10</div><div>       EPS_PGF_ORB     1.0000000000000000E-08</div><div>       METHOD  GPW</div><div>     &END QS</div><div>     &MGRID</div><div>       NGRIDS  4</div><div>       CUTOFF     5.0000000000000000E+02</div><div>       REL_CUTOFF     8.0000000000000000E+01</div><div>     &END MGRID</div><div>     &XC</div><div>       DENSITY_CUTOFF     1.0000000000000000E-10</div><div>       GRADIENT_CUTOFF     1.0000000000000000E-10</div><div>       TAU_CUTOFF     1.0000000000000000E-10</div><div>       &XC_FUNCTIONAL  NO_SHORTCUT</div><div>         &PBE  T</div><div>         &END PBE</div><div>       &END XC_FUNCTIONAL</div><div>     &END XC</div><div>     &POISSON</div><div>       POISSON_SOLVER  PERIODIC</div><div>       PERIODIC  XYZ</div><div>     &END POISSON</div><div>   &END DFT</div><div><br></div><div>So I tried to add this block:</div><div><br></div><div><div>     &LS_SCF</div><div>       PURIFICATION_METHOD TRS4</div><div>       EPS_FILTER 1E-6  </div><div>       EPS_SCF    1E-7</div><div>       MAX_SCF    20</div><div>       S_PRECONDITIONER ATOMIC</div><div>     &END</div></div><div><br></div><div>and it does not work. Then, I tried:</div><div><br></div><div><div>     &LS_SCF</div><div>       PURIFICATION_METHOD TRS4</div><div>       EPS_FILTER 1E-6  </div><div>       EPS_SCF    1E-7</div><div>       MAX_SCF    30</div><div>       S_PRECONDITIONER NONE</div><div>       &CURVY_STEPS</div><div>       &END CURVY_STEPS</div><div>     &END</div></div><div><br></div><div>And again, it failed. In both cases, at first the optimization goes well, I reach a point where the system is almost optimized and then everything goes wrong. If anyone can tell me what I am doing wrong, I would be grateful. More generally, if anyone could tell me which algorithms I should use to maximize the efficiency of the calculation for large systems (up to about 100000 atoms), while having decent accuracy and stability, I would be grateful. I have checked the various options in SCF, OT, OUTER_SCF, LS_SCF and MOTION, but it is unclear which ones are the most suitable and can work well together.</div><div><br></div><div>Thank you,</div><div>Pierre</div><div><br></div></div>