[CP2K-user] [CP2K:14313] Optimisation of Si(100) slab does not converge
fa...@gmail.com
fabia... at gmail.com
Mon Nov 30 09:36:28 UTC 2020
Dear Dmitrii,
Your cell is not correct, all angles should be 90 degree.
The total energy you reported is positive. This indicates that your
structure is far from stable. Where did you get the coordinates from? Did
you optimize bulk Si with cp2k before creating the slab?
Cheers,
Fabian
On Monday, 30 November 2020 at 04:52:41 UTC+1 dre... at gmail.com wrote:
> Thank you for your reply!
> I need to compare adsorption energy for organic molecules on different
> surfaces (Na, Au, Si, and C). I successfully used diagonalization for Na
> and Au as they are metals. However, I am not sure if I use OT for Si and
> C, it would be accurate to compare Eads with those of Na and Au. What do
> you think?
> Do you have any suggestion how to expand Si surface to allow
> reconstructing?
> Do you think fo graphite I also cannot use unitary cell?
>
> Regards,
> Dmitrii
>
>
>
> On Monday, November 30, 2020 at 2:21:43 PM UTC+11 Lucas Lodeiro wrote:
>
>> Hi DMITRII,
>>
>> I am not an expert on silicon systems, but I read plenty times that (001)
>> silicon surfaces presents reconstruction when the surface is expanded (the
>> reconstruction is not possible for the unitary surface cell), as your case.
>> This could explain why do you have this problem with silicon and not with
>> Au and Na. The first time that I read about this was "density functional
>> theory: a practical introduction" of Janice Steckel and David Scholl, and
>> in some papers, like this:
>> https://www.researchgate.net/figure/Pictorial-view-of-the-100-surface-of-silicon-a-unrelaxed-surface-b-relaxed_fig5_282648906
>>
>> Another thing is, Why do you use diagonalization instead of OT method? I
>> expect that silicon system (even surface systems) were not metals... then
>> the OT method will be faster.
>>
>> Regards
>> Lucas Lodeiro
>>
>> El dom, 29 nov 2020 a las 22:24, DMITRII Drugov (<dre... at gmail.com>)
>> escribió:
>>
>>> Good day dear CP2K users, could you please have a look at my .in file to
>>> see did I specify ABC cell parameters and ALPHA_BETA_GAMMA for Si(100) with
>>> 16.29210*16.29210*5.76 A^3. I am using TZVP-MOLOPT-GTH and I allowed two
>>> top layer to relax while other bulk part is frozen. In total it only 42
>>> atoms our of 85 is under relaxation. My system didn't reach any convergence
>>> yet, despite I run it for 30 hour on supercomputer with 120 GB memory and
>>> 36 CPU. However, I used same .in file for gold and Na slabs and with around
>>> 70-60 atoms with DZVP, and it converged within 1-2 hours. Do you think it's
>>> normal for TZVP, should I wait longer? Or smth wrong with my .in setting?
>>>
>>> Regards,
>>> Dmitrii
>>>
>>> -------- Informations at step = 3 ------------
>>> Optimization Method = LBFGS
>>> Total Energy = 2.9038335460 <(903)%20833-5460>
>>> Real energy change = -14.0555995686
>>> Decrease in energy = YES
>>> Used time = 15200.012
>>>
>>> Convergence check :
>>> Max. step size = 0.3482436711
>>> Conv. limit for step size = 0.0030000000
>>> Convergence in step size = NO
>>> RMS step size = 0.0816981844
>>> Conv. limit for RMS step = 0.0015000000
>>> Convergence in RMS step = NO
>>> Max. gradient = 1.6813091629 <(681)%20309-1629>
>>> Conv. limit for gradients = 0.0004500000
>>> Conv. for gradients = NO
>>> RMS gradient = 0.4396446909
>>> Conv. limit for RMS grad. = 0.0003000000
>>> Conv. for gradients = NO
>>>
>>> my setting are:
>>> &GLOBAL
>>> PROJECT Si100_optimisation
>>> RUN_TYPE GEO_OPT
>>> PRINT_LEVEL MEDIUM
>>> &END GLOBAL
>>> &FORCE_EVAL
>>> METHOD QS
>>> &DFT
>>> BASIS_SET_FILE_NAME BASIS_MOLOPT
>>> POTENTIAL_FILE_NAME GTH_POTENTIALS
>>> CHARGE 0
>>> MULTIPLICITY 1
>>> &MGRID
>>> CUTOFF 800
>>> NGRIDS 5
>>> REL_CUTOFF 70
>>> &END MGRID
>>> &QS
>>> EPS_DEFAULT 1.0E-12
>>> WF_INTERPOLATION ASPC
>>> &END QS
>>> &SCF
>>> SCF_GUESS ATOMIC
>>> EPS_SCF 1.0E-7
>>> MAX_SCF 1000
>>> CHOLESKY INVERSE
>>> ADDED_MOS 20
>>> &SMEAR ON
>>> METHOD FERMI_DIRAC
>>> ELECTRONIC_TEMPERATURE [K] 300
>>> &END SMEAR
>>> &DIAGONALIZATION
>>> ALGORITHM STANDARD
>>> &END DIAGONALIZATION
>>> &MIXING
>>> METHOD BROYDEN_MIXING
>>> ALPHA 0.4
>>> NBROYDEN 8
>>> &END MIXING
>>> &END SCF
>>> &XC
>>> &XC_FUNCTIONAL
>>> &PBE
>>> &END PBE
>>> &END XC_FUNCTIONAL
>>> &vdW_POTENTIAL
>>> DISPERSION_FUNCTIONAL PAIR_POTENTIAL
>>> &PAIR_POTENTIAL
>>> PARAMETER_FILE_NAME dftd3.dat
>>> TYPE DFTD3
>>> REFERENCE_FUNCTIONAL PBE
>>> R_CUTOFF 15.0
>>> &END PAIR_POTENTIAL
>>> &END vdW_POTENTIAL
>>> &END XC
>>> &POISSON
>>> PERIODIC xy
>>> POISSON_SOLVER ANALYTIC
>>> &END POISSON
>>> &END DFT
>>> &SUBSYS
>>> &CELL
>>> ABC 16.404 16.404 50.0
>>> ALPHA_BETA_GAMMA 90.0 90.0 120.0
>>> PERIODIC xy
>>> &END CELL
>>> &COORD
>>> Si 17.64977 -17.64977 -4.07302
>>> Si 17.64977 -12.21907 -4.07302
>>> Si 17.64977 -6.78837 -4.07302
>>> Si 17.64977 -1.35767 -4.07302
>>> Si 12.21907 -17.64977 -4.07302
>>> Si 16.29210 -16.29210 -5.43070
>>> Si 13.57675 -13.57675 -5.43070
>>> Si 14.93442 -14.93442 -4.07302
>>> Si 12.21907 -12.21907 -4.07302
>>> Si 16.29210 -10.86140 -5.43070
>>> Si 13.57675 -8.14605 -5.43070
>>> Si 14.93442 -9.50372 -4.07302
>>> Si 12.21907 -6.78837 -4.07302
>>> Si 16.29210 -5.43070 -5.43070
>>> Si 13.57675 -2.71535 -5.43070
>>> Si 14.93442 -4.07302 -4.07302
>>> Si 12.21907 -1.35767 -4.07302
>>> Si 6.78837 -17.64977 -4.07302
>>> Si 10.86140 -16.29210 -5.43070
>>> Si 8.14605 -13.57675 -5.43070
>>> Si 9.50372 -14.93442 -4.07302
>>> Si 6.78837 -12.21907 -4.07302
>>> Si 10.86140 -10.86140 -5.43070
>>> Si 8.14605 -8.14605 -5.43070
>>> Si 9.50372 -9.50372 -4.07302
>>> Si 6.78837 -6.78837 -4.07302
>>> Si 10.86140 -5.43070 -5.43070
>>> Si 8.14605 -2.71535 -5.43070
>>> Si 9.50372 -4.07302 -4.07302
>>> Si 6.78837 -1.35767 -4.07302
>>> Si 1.35767 -17.64977 -4.07302
>>> Si 5.43070 -16.29210 -5.43070
>>> Si 2.71535 -13.57675 -5.43070
>>> Si 4.07302 -14.93442 -4.07302
>>> Si 1.35767 -12.21907 -4.07302
>>> Si 5.43070 -10.86140 -5.43070
>>> Si 2.71535 -8.14605 -5.43070
>>> Si 4.07302 -9.50372 -4.07302
>>> Si 1.35767 -6.78837 -4.07302
>>> Si 5.43070 -5.43070 -5.43070
>>> Si 2.71535 -2.71535 -5.43070
>>> Si 4.07302 -4.07302 -4.07302
>>> Si 1.35767 -1.35767 -4.07302
>>> Si 17.64977 -14.93442 -1.35767
>>> Si 17.64977 -9.50372 -1.35767
>>> Si 17.64977 -4.07302 -1.35767
>>> Si 14.93442 -17.64977 -1.35767
>>> Si 13.57675 -16.29210 -2.71535
>>> Si 16.29210 -13.57675 -2.71535
>>> Si 14.93442 -12.21907 -1.35767
>>> Si 12.21907 -14.93442 -1.35767
>>> Si 13.57675 -10.86140 -2.71535
>>> Si 16.29210 -8.14605 -2.71535
>>> Si 14.93442 -6.78837 -1.35767
>>> Si 12.21907 -9.50372 -1.35767
>>> Si 13.57675 -5.43070 -2.71535
>>> Si 16.29210 -2.71535 -2.71535
>>> Si 14.93442 -1.35767 -1.35767
>>> Si 12.21907 -4.07302 -1.35767
>>> Si 9.50372 -17.64977 -1.35767
>>> Si 8.14605 -16.29210 -2.71535
>>> Si 10.86140 -13.57675 -2.71535
>>> Si 9.50372 -12.21907 -1.35767
>>> Si 6.78837 -14.93442 -1.35767
>>> Si 8.14605 -10.86140 -2.71535
>>> Si 10.86140 -8.14605 -2.71535
>>> Si 9.50372 -6.78837 -1.35767
>>> Si 6.78837 -9.50372 -1.35767
>>> Si 8.14605 -5.43070 -2.71535
>>> Si 10.86140 -2.71535 -2.71535
>>> Si 9.50372 -1.35767 -1.35767
>>> Si 6.78837 -4.07302 -1.35767
>>> Si 4.07302 -17.64977 -1.35767
>>> Si 2.71535 -16.29210 -2.71535
>>> Si 5.43070 -13.57675 -2.71535
>>> Si 4.07302 -12.21907 -1.35767
>>> Si 1.35767 -14.93442 -1.35767
>>> Si 2.71535 -10.86140 -2.71535
>>> Si 5.43070 -8.14605 -2.71535
>>> Si 4.07302 -6.78837 -1.35767
>>> Si 1.35767 -9.50372 -1.35767
>>> Si 2.71535 -5.43070 -2.71535
>>> Si 5.43070 -2.71535 -2.71535
>>> Si 4.07302 -1.35767 -1.35767
>>> Si 1.35767 -4.07302 -1.35767
>>> &END COORD
>>> &KIND Au
>>> BASIS_SET DZVP-MOLOPT-SR-GTH
>>> POTENTIAL GTH-PBE-q11
>>> &END KIND
>>> &KIND F
>>> BASIS_SET TZVP-MOLOPT-GTH
>>> POTENTIAL GTH-PBE-q7
>>> &END KIND
>>> &KIND O
>>> BASIS_SET TZVP-MOLOPT-GTH
>>> POTENTIAL GTH-PBE-q6
>>> &END KIND
>>> &KIND C
>>> BASIS_SET TZVP-MOLOPT-GTH
>>> POTENTIAL GTH-PBE-q4
>>> &END KIND
>>> &KIND Si
>>> BASIS_SET TZVP-MOLOPT-GTH
>>> POTENTIAL GTH-PBE-q4
>>> &END KIND
>>> &KIND S
>>> BASIS_SET TZVP-MOLOPT-GTH
>>> POTENTIAL GTH-PBE-q6
>>> &END KIND
>>> &KIND N
>>> BASIS_SET TZVP-MOLOPT-GTH
>>> POTENTIAL GTH-PBE-q5
>>> &END KIND
>>> &KIND Na
>>> BASIS_SET DZVP-MOLOPT-SR-GTH
>>> POTENTIAL GTH-PBE-q9
>>> &END KIND
>>> &END SUBSYS
>>> &END FORCE_EVAL
>>> &MOTION
>>> &CONSTRAINT
>>> &FIXED_ATOMS
>>> COMPONENTS_TO_FIX XYZ
>>> LIST 1..43
>>> &END FIXED_ATOMS
>>> &END CONSTRAINT
>>> &GEO_OPT
>>> OPTIMIZER LBFGS
>>> MAX_ITER 300
>>> &END GEO_OPT
>>> &END MOTION
>>> [image: Screen Shot 2020-11-30 at 12.12.28 pm.png]
>>>
>>> --
>>> 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 cp... at googlegroups.com.
>>> To view this discussion on the web visit
>>> https://groups.google.com/d/msgid/cp2k/162639f5-a816-408c-b798-55a572c13fffn%40googlegroups.com
>>> <https://groups.google.com/d/msgid/cp2k/162639f5-a816-408c-b798-55a572c13fffn%40googlegroups.com?utm_medium=email&utm_source=footer>
>>> .
>>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20201130/007fa197/attachment.htm>
More information about the CP2K-user
mailing list