Here is the input and structure file.<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年8月25日星期四 UTC+8 15:45:28<zhouk...@gmail.com> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear all,<div><br></div><div>Now I am simulation the MXene (Ti3C2O2)+water system.  But SCF can not converge. </div><div>The initial structure is from the MD simulation using VASP (PBE+D3) and thus the problem is not from bad positions of atoms.</div><div><br></div><div>Can anyone give me some suggestions?<br></div><div><br></div><div><div>Best,</div><div>Justin</div></div><div>-----------------------------------------------------------------</div><div>&FORCE_EVAL<br>  ! the electronic structure part of CP2K is named Quickstep<br>  METHOD Quickstep<br>  STRESS_TENSOR ANALYTICAL<br>  &DFT<br>    BASIS_SET_FILE_NAME   ./BASIS_MOLOPT  <br>    POTENTIAL_FILE_NAME   ./POTENTIAL <br>    WFN_RESTART_FILE_NAME ./MX-RESTART.wfn<br>    &MGRID<br>       CUTOFF [Ry] 400 <br>    &END<br>    &QS<br>       ! use the GPW method (i.e. pseudopotential based calculations with the Gaussian and Plane Waves scheme).<br>       METHOD GPW <br>       ! default threshold for numerics ~ roughly numerical accuracy of the total energy per electron,<br>       ! sets reasonable values for all other thresholds.<br>       EPS_DEFAULT 1.0E-10 <br>       ! used for MD, the method used to generate the initial guess.<br>       EXTRAPOLATION ASPC <br>       EXTRAPOLATION_ORDER 4<br>    &END<br><br>    &POISSON<br>       PERIODIC XYZ ! the default, gas phase systems should have 'NONE' and a wavelet solver<br>    &END<br>    &PRINT<br>     ! output atomic charge<br>     &MULLIKEN <br>       ! FILENAME ${SYSTEM}<br>       &EACH<br>            MD 10<br>       &END EACH<br>     &END MULLIKEN<br>     ! output others <br>    &END<br>    ! use the OT METHOD for robust and efficient SCF, suitable for all non-metallic systems.<br>    &SCF                              <br>      SCF_GUESS RESTART !!!#####ATOMIC ! can be used to RESTART an interrupted calculation<br>      MAX_SCF 1000<br>      EPS_SCF 1.0E-7 <br>      &OUTER_SCF ! repeat the inner SCF cycle 10 times<br>        MAX_SCF 10<br>        EPS_SCF 1.0E-7<br>      &END<br><br>      &OT ON<br>        MINIMIZER DIIS<br>      &END OT<br><br><br><br>    &END SCF<br><br>    ! specify the exchange and correlation treatment<br>    &XC<br>      ! use a PBE functional <br>      &XC_FUNCTIONAL PBE<br>      &END XC_FUNCTIONAL<br>      ! adding Grimme's D3 correction (by default without C9 terms) <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>         &END<br>      &END VDW_POTENTIAL<br><br>      &XC_GRID<br>         XC_SMOOTH_RHO NN50<br>         XC_DERIV NN50_SMOOTH<br>      &END<br><br>    &END XC<br>  &END DFT<br><br> <br>  ! description of the system<br>  &SUBSYS<br>    &CELL <br>      ! unit cells that are orthorhombic are more efficient with CP2K<br>      ABC [angstrom] 12.0944890976000003 10.4741353988999997 16.2000007629000002<br>    &END CELL<br><br>    ! atom coordinates can be in the &COORD section,<br>    ! or provided as an external file.<br>    &TOPOLOGY<br>      COORD_FILE_NAME MX-wt.xyz <br>      COORD_FILE_FORMAT XYZ<br>    &END<br><br>    ! MOLOPT basis sets are fairly costly,<br>    ! but in the 'DZVP-MOLOPT-SR-GTH' available for all elements<br>    ! their contracted nature makes them suitable<br>    ! for condensed and gas phase systems alike.<br><br>    &KIND C<br>      BASIS_SET DZVP-MOLOPT-SR-GTH <br>      POTENTIAL GTH-PBE<br>    &END KIND<br>    &KIND Ti<br>      BASIS_SET DZVP-MOLOPT-SR-GTH <br>      POTENTIAL GTH-PBE<br>    &END KIND<br>    &KIND H<br>      BASIS_SET DZVP-MOLOPT-SR-GTH <br>      POTENTIAL GTH-PBE<br>    &END KIND<br>    &KIND O<br>      BASIS_SET DZVP-MOLOPT-SR-GTH <br>      POTENTIAL GTH-PBE<br>    &END KIND<br>    &KIND Li<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      POTENTIAL GTH-PBE<br>    &END KIND<br>  &END SUBSYS<br>&END FORCE_EVAL<br></div><div><br></div><div>-----------------------</div><div><br></div><div><br></div><div><br></div></blockquote></div>

<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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/23f3d72c-da77-410f-bc99-b37f8434fd03n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/23f3d72c-da77-410f-bc99-b37f8434fd03n%40googlegroups.com</a>.<br />