<div dir="ltr"><font size="4">Dear Sir/Madam,</font><div><font size="4"><br></font></div><div><font size="4">I was running a large AIMD simulation (about 5200 basis functions) for Au surface supported CeO</font><font size="1">2</font><font size="4"> cluster ( (Au(111) surface has</font><span style="font-size: large;">180 Au atoms</span><font size="4">, and Ce</font><font size="1">10</font><font size="4">O</font><font size="1">20</font><font size="4"> cluster). </font></div><div><font size="4"><br></font></div><div><font size="4">At the beginning, simulation goes well, then it has a very serious convergence problems at around step number of 350. it can't get converged no matter what I do to the SCF control. </font></div><div><font size="4"><br></font></div><div><font size="4">Even if I make it converged coincidentally, but it will not converge after considering Kinetic energy and temperature.</font></div><div><font size="4"><br></font></div><div><font size="4">I tried every method I can find， but I can‘t figure it out.</font></div><div><font size="4"><br></font></div><div><font size="4">I am depressed and frustrated now, hope someone could help me out with it.  Any suggestions are appreciated!!</font></div><div><font size="4"><br></font></div><div><font size="4">The input output file were enclosed in attachment, part of input file is as follows</font></div><div><font size="4"><br></font></div><div><font size="4"><div>     </div><div><br></div><div>     MULTIPLICITY  1</div><div>     PLUS_U_METHOD  MULLIKEN_CHARGES</div><div>     &SCF</div><div>       MAX_SCF  100</div><div>       EPS_SCF     7.0E-05</div><div>       CHOLESKY  INVERSE</div><div>       SCF_GUESS  ATOMIC</div><div>       ADDED_MOS  400</div><div>       &DIAGONALIZATION  T</div><div>         ALGORITHM  STANDARD</div><div>       &END DIAGONALIZATION</div><div>       &OUTER_SCF  T</div><div>         EPS_SCF     8.0000000000000001E-05</div><div>         MAX_SCF  5</div><div>       &END OUTER_SCF</div><div>       &SMEAR  T</div><div>    #     METHOD  ENERGY_WINDOW</div><div>         METHOD  FERMI_DIRAC</div><div>         ELECTRONIC_TEMPERATURE     3.0000000000000000E+02</div><div>       &END SMEAR</div><div>       &MIXING  T</div><div>         METHOD  BROYDEN_MIXING</div><div>         ALPHA     0.1</div><div>         BETA      1.5</div><div>         NBUFFER  8</div><div>       &END MIXING</div><div>       &PRINT</div><div>         &RESTART  SILENT</div><div>           FILENAME =RESTART.wfn</div><div>           &EACH</div><div>             MD  5</div><div>           &END EACH</div><div>         &END RESTART</div><div>       &END PRINT</div><div>     &END SCF</div><div>     &QS</div><div>       EPS_DEFAULT     1.0000000000000000E-11</div><div>       EPS_PGF_ORB     1.0000000000000001E-05</div><div>       EXTRAPOLATION  ASPC</div><div>       EXTRAPOLATION_ORDER  3</div><div>     &END QS</div><div>     &MGRID</div><div>       CUTOFF     5.0000000000000000E+02</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_GRID</div><div>         XC_SMOOTH_RHO  NN10</div><div>         XC_DERIV  SPLINE2_SMOOTH</div><div>       &END XC_GRID</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>     &PRINT</div><div>       &MULLIKEN  SILENT</div><div>         FILENAME =md.mulliken</div><div>         &EACH</div><div>           MD  5</div><div>         &END EACH</div><div>       &END MULLIKEN</div><div>     &END PRINT</div><div>   &END DFT</div><div>   &SUBSYS</div></font></div><div><font size="4"><br></font></div><div><font size="4"><br></font></div><div><font size="4"><br></font></div><div><font size="4"><br></font></div></div>