<div><p><strong>Dear CP2K Users,</strong></p><p>I am running a <strong>geometry optimization for a molecule adsorption system with vacuum along the y-axis</strong>, but I am struggling to achieve <strong>electronic SCF convergence even after 1000 iterations</strong>. I would appreciate any suggestions on improving the convergence.</p><p>Below is my input script:</p></div>&GLOBAL<br />  Â PRINT_LEVEL Â MEDIUM<br />  Â PROJECT_NAME mol_absorp_OPT<br />  Â RUN_TYPE Â GEO_OPT<br />  Â EXTENDED_FFT_LENGTHS T<br />&END GLOBAL<br /><br />&FORCE_EVAL<br />  METHOD Â QS<br />  &DFT<br />  Â  UKS .FALSE. <br />  Â  BASIS_SET_FILE_NAME BASIS_MOLOPT<br />  Â  POTENTIAL_FILE_NAME POTENTIAL <br />  Â  &MGRID<br />  Â  Â  NGRIDS 4<br />  Â  Â  CUTOFF 600<br />  Â  Â  REL_CUTOFF 60<br />  Â  &END MGRID<br />  Â  &QS<br />  Â  Â  METHOD GPW<br />  Â  Â  EPS_DEFAULT 1.000E-12<br />  Â  &END QS<br />  Â  &POISSON<br />  Â  Â  PERIODIC XZ<br />  Â  Â  PSOLVER MT<br />  Â  &END POISSON<br />  Â  &SCF<br />  Â  Â  &DIAGONALIZATION ON<br />  Â  Â  &END DIAGONALIZATION<br />  Â  Â  &MIXING<br />  Â  Â  Â  ALPHA 0.5<br />  Â  Â  Â  METHOD PULAY_MIXING<br />  Â  Â  Â  NBROYDEN 8<br />  Â  Â  &END MIXING<br />  Â  Â  EPS_SCF 1.0E-5<br />  Â  Â  SCF_GUESS ATOMIC<br />  Â  Â  MAX_SCF Â 5000<br />  Â  Â  ADDED_MOS 100<br />  Â  &END SCF<br />  Â  &XC<br />  Â  Â  &XC_FUNCTIONAL<br />  Â  Â  Â  Â &PBE<br />  Â  Â  Â  Â &END PBE<br />  Â  Â  &END XC_FUNCTIONAL<br />  Â  Â  &VDW_POTENTIAL<br />  Â  Â  Â  Â POTENTIAL_TYPE PAIR_POTENTIAL <br />  Â  Â  Â  Â &PAIR_POTENTIAL<br />  Â  Â  Â  Â  Â  PARAMETER_FILE_NAME dftd3.dat<br />  Â  Â  Â  Â  Â  TYPE DFTD3<br />  Â  Â  Â  Â  Â  REFERENCE_FUNCTIONAL PBE<br />  Â  Â  Â  Â  Â  R_CUTOFF [angstrom] 16<br />  Â  Â  Â  Â &END<br />  Â  Â  &END VDW_POTENTIAL<br />  Â  &END XC<br />  Â  &PRINT<br />  Â  Â  &MULLIKEN OFF<br />  Â  Â  &END<br />  Â  Â  &HIRSHFELD OFF<br />  Â  Â  &END<br />  Â  &END PRINT<br />  &END DFT<br />  &SUBSYS<br />  Â  &CELL<br />  Â  Â  A Â  24.8975601196000014 Â  Â 0.0000000000000000 Â  Â 0.0000000000000000<br />  Â  Â  B Â  0.0000000000000000 Â  Â 55.6124382018999981 Â  Â 0.0000000000000000<br />  Â  Â  C Â  0.0000000000000000 Â  Â  0.0000000000000000 Â  24.9615383147999985<br />  Â  Â  PERIODIC Â XZ<br />  Â  &END CELL<br />  Â  &TOPOLOGY<br />  Â  Â  COORD_FILE_NAME ABC.xyz<br />  Â  Â  COORD_FILE_FORMAT XYZ<br />  Â  &END<br />  Â  &KIND Â Cs<br />  Â  Â  BASIS_SET DZVP-MOLOPT-SR-GTH <br />  Â  Â  POTENTIAL GTH-PBE-q9 <br />  Â  &END KIND<br />  Â  &KIND Pb<br />  Â  Â  BASIS_SET DZVP-MOLOPT-SR-GTH <br />  Â  Â  POTENTIAL GTH-PBE-q4<br />  Â  &END KIND<br />  Â  &KIND Br <br />  Â  Â  BASIS_SET Â DZVP-MOLOPT-SR-GTH <br />  Â  Â  POTENTIAL Â GTH-PBE-q7<br />  Â  &END KIND<br />  Â  &KIND C<br />  Â  Â  BASIS_SET DZVP-MOLOPT-GTH<br />  Â  Â  POTENTIAL GTH-PBE-q4<br />  Â  &END KIND<br />  Â  &KIND N<br />  Â  Â BASIS_SET DZVP-MOLOPT-GTH<br />  Â  Â POTENTIAL GTH-PBE-q5<br />  Â  &END KIND<br />  Â  &KIND H<br />  Â  Â BASIS_SET DZVP-MOLOPT-GTH<br />  Â  Â POTENTIAL GTH-PBE-q1<br />  Â  &END KIND<br />  &END SUBSYS<br />&END FORCE_EVAL<br /><br />&MOTION<br /> &GEO_OPT<br />  Â OPTIMIZER BFGS<br />  Â MAX_ITER 500<br />  Â MAX_DR [bohr] 0.01<br />  Â MAX_FORCE [hartree*bohr^-1] 0.01<br />  Â RMS_DR [bohr] 0.01<br />  Â RMS_FORCE [hartree*bohr^-1] 0.01<br /> &END GEO_OPT<br /><br /> &CONSTRAINT<br />  Â &FIXED_ATOMS<br />  Â  Â LIST 1 2 3 4 5 6 10 11 12 13 14 15 .........................<br />  Â &END FIXED_ATOMS<br /> &END CONSTRAINT<br /><br /> &PRINT<br />  Â &TRAJECTORY<br />  Â  Â &EACH<br />  Â  Â  Â GEO_OPT 5<br />  Â  Â &END EACH<br />  Â  Â FORMAT XYZ<br />  Â &END TRAJECTORY<br />  Â &VELOCITIES OFF  <br />  Â &END VELOCITIES<br />  Â &FORCES ON  Â  Â Â <br />  Â &END FORCES<br />  Â &RESTART<br />  Â  Â BACKUP_COPIES 1<br />  Â  Â &EACH<br />  Â  Â  Â GEO_OPT 10<br />  Â  Â &END EACH<br />  Â &END RESTART<br /> &END PRINT<br />&END MOTION

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/24c304f2-56c9-4c28-bdf8-7334e3553bc0n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/24c304f2-56c9-4c28-bdf8-7334e3553bc0n%40googlegroups.com</a>.<br />