[CP2K-user] [CP2K:16422] Cell Optimization of Platinum (Pt) with fcc (111) surface
GENG YUAN
yuangeng0821 at gmail.com
Mon Jan 3 20:34:56 UTC 2022
Dear CP2K Users,
Hope you are going great and happy new year.
I am running the cell optimization for a single bulk Pt (134 atoms) with
fcc 111 surface. By checking my output file, I notice that the cell
optimization takes 446 steps to converge (which is almost close to the
maximum number of optimization steps:500). Besides, a certain degree of
structure deformation is observed during optimization (please see the
attached files and screenshots).
I am wondering if there are any other ways to let the cell optimization
converge faster (in other words, converge in fewer steps) and whether the
above-mentioned deformation looks reasonable?
Below is the input: (attached are input file, output file, original
coordinate, trajectory, and the screenshots of the trajectory)
------------------------------------------------------------
&GLOBAL
PROJECT Pt_cellopt
RUN_TYPE CELL_OPT
PRINT_LEVEL LOW
&END GLOBAL
&FORCE_EVAL
METHOD QS
STRESS_TENSOR ANALYTICAL
&SUBSYS
&CELL
ABC 13.7 13.7 13.7
SYMMETRY CUBIC
&CELL_REF
ABC 13.7*1.5 13.7*1.5 13.7*1.5
&END
&END CELL
&TOPOLOGY
COORD_FILE_NAME ./Pt134.xyz
COORDINATE xyz
&END
&KIND Pt
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q18
&END KIND
&END SUBSYS
&DFT
BASIS_SET_FILE_NAME ./BASIS_MOLOPT
POTENTIAL_FILE_NAME ./GTH_POTENTIALS
&QS
EPS_DEFAULT 1.0E-12
&END QS
&MGRID
CUTOFF 400
NGRIDS 5
REL_CUTOFF 40
&END MGRID
&SCF
SCF_GUESS ATOMIC
EPS_SCF 1.0E-06
MAX_SCF 200
ADDED_MOS 100
&OUTER_SCF
MAX_SCF 50
EPS_SCF 1.0E-6
&END
&DIAGONALIZATION T
ALGORITHM STANDARD
&END DIAGONALIZATION
&MIXING T
METHOD BROYDEN_MIXING
ALPHA 0.1
NBROYDEN 8
&END MIXING
&SMEAR ON
METHOD FERMI_DIRAC
ELECTRONIC_TEMPERATURE [K] 500
&END SMEAR
&PRINT
&RESTART ON
&END RESTART
&END PRINT
&END SCF
&XC
&XC_FUNCTIONAL
&LIBXC
FUNCTIONAL XC_GGA_X_PBE_R
&END
&LIBXC
FUNCTIONAL XC_GGA_C_PBE
&END
&END XC_FUNCTIONAL
&VDW_POTENTIAL
POTENTIAL_TYPE PAIR_POTENTIAL
&PAIR_POTENTIAL
LONG_RANGE_CORRECTION
PARAMETER_FILE_NAME dftd3.dat
TYPE DFTD3
REFERENCE_FUNCTIONAL revPBE
&END PAIR_POTENTIAL
&END VDW_POTENTIAL
&XC_GRID
XC_DERIV NN50_SMOOTH
&END XC_GRID
&END XC
&END DFT
&END FORCE_EVAL
&MOTION
&CELL_OPT
EXTERNAL_PRESSURE 0
TYPE DIRECT_CELL_OPT
KEEP_SYMMETRY TRUE
MAX_DR 1.0E-05
MAX_FORCE 1.0E-05
RMS_DR 1.0E-05
RMS_FORCE 1.0E-05
MAX_ITER 500
OPTIMIZER BFGS
&END CELL_OPT
&END MOTION
Many thanks in advance,
Sincerely,
Geng
--
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/c8bce8c7-3186-41b6-a64d-1c4296879204n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220103/b5705b7c/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: traj.xyz
Type: chemical/x-xyz
Size: 3915273 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220103/b5705b7c/attachment.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: out
Type: application/octet-stream
Size: 2334591 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220103/b5705b7c/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cell.inp
Type: chemical/x-gamess-input
Size: 2007 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220103/b5705b7c/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Pt134.xyz
Type: chemical/x-xyz
Size: 3276 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220103/b5705b7c/attachment-0001.xyz>
More information about the CP2K-user
mailing list