[CP2K-user] [CP2K:17544] Bad convergence performance for the MXene system

[zhouk...@gmail.com]

Dear all,

Now I am simulation the MXene (Ti3C2O2)+water system.  But SCF can not 
converge. 
The initial structure is from the MD simulation using VASP (PBE+D3) and 
thus the problem is not from bad positions of atoms.

Can anyone give me some suggestions?

Best,
Justin
-----------------------------------------------------------------
&FORCE_EVAL
  ! the electronic structure part of CP2K is named Quickstep
  METHOD Quickstep
  STRESS_TENSOR ANALYTICAL
  &DFT
    BASIS_SET_FILE_NAME   ./BASIS_MOLOPT  
    POTENTIAL_FILE_NAME   ./POTENTIAL 
    WFN_RESTART_FILE_NAME ./MX-RESTART.wfn
    &MGRID
       CUTOFF [Ry] 400 
    &END
    &QS
       ! use the GPW method (i.e. pseudopotential based calculations with 
the Gaussian and Plane Waves scheme).
       METHOD GPW 
       ! default threshold for numerics ~ roughly numerical accuracy of the 
total energy per electron,
       ! sets reasonable values for all other thresholds.
       EPS_DEFAULT 1.0E-10 
       ! used for MD, the method used to generate the initial guess.
       EXTRAPOLATION ASPC 
       EXTRAPOLATION_ORDER 4
    &END

    &POISSON
       PERIODIC XYZ ! the default, gas phase systems should have 'NONE' and 
a wavelet solver
    &END
    &PRINT
     ! output atomic charge
     &MULLIKEN 
       ! FILENAME ${SYSTEM}
       &EACH
            MD 10
       &END EACH
     &END MULLIKEN
     ! output others 
    &END
    ! use the OT METHOD for robust and efficient SCF, suitable for all 
non-metallic systems.
    &SCF                              
      SCF_GUESS RESTART !!!#####ATOMIC ! can be used to RESTART an 
interrupted calculation
      MAX_SCF 1000
      EPS_SCF 1.0E-7 
      &OUTER_SCF ! repeat the inner SCF cycle 10 times
        MAX_SCF 10
        EPS_SCF 1.0E-7
      &END

      &OT ON
        MINIMIZER DIIS
      &END OT



    &END SCF

    ! specify the exchange and correlation treatment
    &XC
      ! use a PBE functional 
      &XC_FUNCTIONAL PBE
      &END XC_FUNCTIONAL
      ! adding Grimme's D3 correction (by default without C9 terms) 
      &VDW_POTENTIAL
         POTENTIAL_TYPE PAIR_POTENTIAL 
         &PAIR_POTENTIAL
            PARAMETER_FILE_NAME dftd3.dat
            TYPE DFTD3
            REFERENCE_FUNCTIONAL PBE
         &END
      &END VDW_POTENTIAL

      &XC_GRID
         XC_SMOOTH_RHO NN50
         XC_DERIV NN50_SMOOTH
      &END

    &END XC
  &END DFT

 
  ! description of the system
  &SUBSYS
    &CELL 
      ! unit cells that are orthorhombic are more efficient with CP2K
      ABC [angstrom] 12.0944890976000003 10.4741353988999997 
16.2000007629000002
    &END CELL

    ! atom coordinates can be in the &COORD section,
    ! or provided as an external file.
    &TOPOLOGY
      COORD_FILE_NAME MX-wt.xyz 
      COORD_FILE_FORMAT XYZ
    &END

    ! MOLOPT basis sets are fairly costly,
    ! but in the 'DZVP-MOLOPT-SR-GTH' available for all elements
    ! their contracted nature makes them suitable
    ! for condensed and gas phase systems alike.

    &KIND C
      BASIS_SET DZVP-MOLOPT-SR-GTH 
      POTENTIAL GTH-PBE
    &END KIND
    &KIND Ti
      BASIS_SET DZVP-MOLOPT-SR-GTH 
      POTENTIAL GTH-PBE
    &END KIND
    &KIND H
      BASIS_SET DZVP-MOLOPT-SR-GTH 
      POTENTIAL GTH-PBE
    &END KIND
    &KIND O
      BASIS_SET DZVP-MOLOPT-SR-GTH 
      POTENTIAL GTH-PBE
    &END KIND
    &KIND Li
      BASIS_SET DZVP-MOLOPT-SR-GTH
      POTENTIAL GTH-PBE
    &END KIND
  &END SUBSYS
&END FORCE_EVAL

-----------------------



-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/5794c0f1-7f6b-405e-b32c-4a88f1d3afa4n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220825/c3c77c21/attachment-0001.htm>