<div dir="ltr">Hello,<div><br></div><div>Is there any specific parameters or settings I should pay more attention for a charged SIC calculation?  I am asking because I struggled a lot for the SCF convergence problem. Any ideas?  </div><div>Btw, I have tried different mixing schemes (pulay and broyden) and mixing parameters alpha. Following is my input:</div><div>----------------------------------------------------------------------------------------------</div><div><div>&GLOBAL</div><div>  PROJECT_NAME TEST</div><div>  RUN_TYPE ENERGY</div><div>  PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div>  METHOD QS</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div>    UKS</div><div>    CHARGE -1</div><div>    &SIC</div><div>     ORBITAL_SET UNPAIRED</div><div>     SIC_METHOD MAURI_US</div><div>    &END SIC</div><div>    &MGRID</div><div>      CUTOFF 250</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>    &SCF</div><div>      ADDED_MOS 40</div><div>      MAX_SCF 100</div><div>      EPS_SCF 1.0e-6</div><div>      CHOLESKY INVERSE</div><div>      SCF_GUESS ATOMIC</div><div>      &DIAGONALIZATION</div><div>       ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &MIXING</div><div>       METHOD Broyden_MIXING</div><div>       ALPHA 0.1</div><div>       NBUFFER 5</div><div>      &END MIXING</div><div>      &SMEAR</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE [K] 300</div><div>      &END SMEAR</div><div>    &END SCF</div><div>    &PRINT</div><div>    &MO_CUBES</div><div>      WRITE_CUBE F</div><div>      NHOMO 2</div><div>      NLUMO 2</div><div>    &END MO_CUBES</div><div>    &END PRINT</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC [angstrom] 10.5345 10.5345 10.5345</div><div>      ALPHA_BETA_GAMMA [deg] 90 90 90</div><div>      PERIODIC XYZ</div><div>      SYMMETRY CUBIC</div><div>    &END CELL</div><div>    &COORD</div><div>.....................</div><div>      SCALED T</div><div>    &END COORD</div><div>    &KIND Cs</div><div>      ELEMENT Cs</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q9</div><div>    &END KIND</div><div>    &KIND In</div><div>      ELEMENT In</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q13</div><div>    &END KIND</div><div>    &KIND Ag</div><div>      ELEMENT Ag</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q11</div><div>    &END KIND</div><div>    &KIND Cl</div><div>      ELEMENT Cl</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q7</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div></div><div>---------------------------------------------------------------------------------</div><div><br></div><div>Best,</div><div>Xiaoming Wang</div></div>