<div dir="ltr">Hello Folks,<div><br></div><div>I'm trying to simulate an Al-C interaction, 4 Al and 1 C using the andreussi solvation model. I'm also trying to use AIMD (DFT) as well to obtain overall bond lengths of the Al-C via radial distribution function post-MD trajectory analysis. During my system setup and start of CP2K, I find that SCF convergences take a long time and convergence exceeds 1, rather than minimizing towards 0. Could someone take a look at my input file and suggest what I might be able to do to achieve convergence of SCF/SCCS?</div><div><br></div><div>Overall, I have done this exact experiment using an explicit experiment with 60+ Al atoms and 1 C. I'd like to see if I can achieve the same results using an implicit model, because this would allow me to reduce simulation time from 2+ months to something less. </div><div><br></div><div>Input also below:</div><div><div>&GLOBAL</div><div> PROJECT AlC1</div><div> RUN_TYPE MD</div><div> PRINT_LEVEL low</div><div>&END GLOBAL</div><div><br></div><div><br></div><div>&FORCE_EVAL</div><div> METHOD Quickstep</div><div> STRESS_TENSOR analytical</div><div> &DFT</div><div>    BASIS_SET_FILE_NAME /shared/centos7/cp2k/cp2k-6.1.0/data/BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME /shared/centos7/cp2k/cp2k-6.1.0/data/POTENTIAL</div><div>  </div><div>  &QS</div><div>    METHOD GPW</div><div>    EPS_DEFAULT 1.0E-12</div><div>    EXTRAPOLATION ASPC</div><div>  &END QS</div><div>  </div><div>  &MGRID</div><div>    NGRIDS 5</div><div>    CUTOFF [Ry] 280</div><div>    REL_CUTOFF 60</div><div>  &END MGRID</div><div><br></div><div>  &SCCS</div><div>   ALPHA [N*m^-1] 0.0</div><div>   BETA [GPa] 0.0</div><div>   DELTA_RHO 2.0E-5</div><div>   DERIVATIVE_METHOD cd5</div><div>   DIELECTRIC_CONSTANT 80</div><div>   GAMMA 0.0</div><div>   EPS_SCCS 1.0E-8</div><div>   EPS_SCF 0.3</div><div>   MAX_ITER 200</div><div>   METHOD Andreussi</div><div>!   METHOD Fattebert-Gygi</div><div>   MIXING 0.6</div><div>   &ANDREUSSI</div><div>    RHO_MAX 0.0001</div><div>    RHO_MIN 0.00001</div><div>   &END ANDREUSSI</div><div>!   &FATTEBERT-GYGI</div><div>!    BETA 1.3</div><div>!    RHO_ZERO 0.00078</div><div>!   &END FATTEBERT-GYGI</div><div>  &END SCCS</div><div>  </div><div>  &SCF</div><div>    SCF_GUESS ATOMIC</div><div>    MAX_SCF 300</div><div>    EPS_SCF 1.0E-7</div><div>    ADDED_MOS 200</div><div>    &DIAGONALIZATION</div><div>      ALGORITHM STANDARD</div><div>    &END DIAGONALIZATION</div><div>    &SMEAR TRUE</div><div><span style="white-space:pre">             </span>    METHOD FERMI_DIRAC</div><div><span style="white-space:pre">              </span>    ELECTRONIC_TEMPERATURE 1000</div><div>    &END SMEAR</div><div>    &MIXING</div><div>      METHOD BROYDEN_MIXING</div><div>      ALPHA 0.1</div><div>      BETA 1.0</div><div>      NBROYDEN 8</div><div>    &END MIXING</div><div>    &PRINT</div><div>      &RESTART OFF</div><div>      &END</div><div>    &END</div><div>  &END SCF</div><div>  </div><div>  &XC</div><div>   &XC_FUNCTIONAL PBE</div><div>   &END XC_FUNCTIONAL</div><div>  &END XC</div><div>  </div><div>  &PRINT</div><div>    &PDOS  SILENT</div><div>      APPEND  T</div><div>      &EACH</div><div>        MD  100</div><div>      &END EACH</div><div>    &END PDOS</div><div>    &SCCS</div><div>      &DIELECTRIC_FUNCTION</div><div>        &EACH</div><div>          MD 100</div><div>        &END EACH</div><div>      &END DIELECTRIC_FUNCTIOn</div><div>    &END SCCS</div><div>  &END PRINT</div><div>  </div><div> &END DFT</div><div> </div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC [angstrom] 15 15 15</div><div>!      PERIODIC XYZ</div><div>!      MULTIPLE_UNIT_CELL 1 1 1</div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME AlC1.xyz</div><div>      COORD_FILE_FORMAT XYZ</div><div>!      MULTIPLE_UNIT_CELL 1 1 1</div><div>    &END</div><div>    &KIND C</div><div>      ELEMENT C</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q4</div><div>    &END KIND</div><div>    &KIND Al</div><div>      ELEMENT Al</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q3</div><div>    &END KIND</div><div>  &END SUBSYS</div><div><br></div><div>&END FORCE_EVAL</div><div><br></div><div>&MOTION</div><div>  &MD</div><div>    ENSEMBLE NPT_F</div><div>    TEMPERATURE [K] 1000</div><div>    TIMESTEP [fs] 0.5</div><div>    STEPS 25000</div><div>    &THERMOSTAT</div><div>      REGION GLOBAL</div><div>      TYPE CSVR</div><div>      &CSVR</div><div>        TIMECON 50</div><div>      &END CSVR</div><div>    &END THERMOSTAT</div><div>    &BAROSTAT</div><div>      PRESSURE [bar] 1.01</div><div>      TIMECON 2000</div><div>    &END BAROSTAT</div><div>  &END MD</div><div>  &PRINT</div><div>    &TRAJECTORY</div><div>      &EACH</div><div>        MD 1</div><div>      &END EACH</div><div>    &END TRAJECTORY</div><div>    &VELOCITIES</div><div>      &EACH</div><div>        MD 1</div><div>      &END EACH</div><div>    &END VELOCITIES</div><div>    &FORCES</div><div>      &EACH</div><div>        MD 1</div><div>      &END EACH</div><div>    &END FORCES</div><div>    &CELL</div><div>      &EACH</div><div>        MD 1</div><div>      &END EACH</div><div>    &END CELL</div><div>    &RESTART_HISTORY</div><div>      &EACH</div><div>        MD 10</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><div><br></div><div><br></div><div>Thanks,</div><div>Dev</div></div>