[CP2K-user] [CP2K:18073] SCF slow down after hundreds of steps
Dragoon
longdragoonae at gmail.com
Thu Nov 24 05:29:01 UTC 2022
Hi,
I'm new to CP2K and trying to setup an AIMD simulation for 1 Mn2+ and 69
water molecules. The job runs on 144 cpus with 5GB/cpu memory.
I noticed that after hundreds of MD steps, the calculation time for each
step increases from 4-8 s/step to tens of seconds. I even got 2453 seconds
and 5114 seconds for two steps.
The 2453s long step took 2 SCF steps:
1 OT DIIS 0.15E+00 1.7 0.00001472 -1292.4330648654
-1.29E+03
2 OT DIIS 0.15E+00 667.8 0.00000884 -1292.4330876796
-2.28E-05
The 5114s long step took 3 SCF steps:
1 OT DIIS 0.15E+00 2766.3 0.00002727 -1292.4373004578
-1.29E+03
2 OT DIIS 0.15E+00 154.7 0.00001826 -1292.4373848508
-8.44E-05
3 OT DIIS 0.15E+00 1.7 0.00000377 -1292.4374551945
-7.03E-05
I also remember that previous runs with 3 GB/cpu or less stopped with
out-of-memory error after hundreds of steps.
I pasted my input file contents below. Any help is greatly appreciated.
Also, any suggestions on cutoff, basis and functional (as well as any other
aspects) are greatly appreciated.
Best,
Zhuoran
&GLOBAL
PROJECT mn_water_nvt
RUN_TYPE MD
PRINT_LEVEL MEDIUM
FFTW_PLAN_TYPE EXHAUSTIVE
&END GLOBAL
&FORCE_EVAL
METHOD Quickstep
&SUBSYS
&KIND Mn
ELEMENT Mn
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q15
&END KIND
&KIND H
ELEMENT H
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q1
&END KIND
&KIND O
ELEMENT O
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q6
&END KIND
&CELL
SYMMETRY CUBIC
ABC 12.777326 12.777326 12.777326
&END CELL
&COORD
@INCLUDE ../mn_water_aimd_init.xyz
&END COORD
&END SUBSYS
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
POTENTIAL_FILE_NAME GTH_POTENTIALS
CHARGE +2
MULTIPLICITY 6
UKS
&QS
EPS_DEFAULT 1.0E-12
EPS_PGF_ORB 1.0E-07
&END QS
&MGRID
NGRIDS 4
CUTOFF 450
REL_CUTOFF 60
&END MGRID
&XC
&XC_FUNCTIONAL PBE
&END XC_FUNCTIONAL
&XC_GRID
XC_DERIV NN10_SMOOTH
XC_SMOOTH_RHO NN10
&END XC_GRID
&vdW_POTENTIAL
DISPERSION_FUNCTIONAL PAIR_POTENTIAL
&PAIR_POTENTIAL
TYPE DFTD3
PARAMETER_FILE_NAME dftd3.dat
REFERENCE_FUNCTIONAL PBE
&END PAIR_POTENTIAL
&END vdW_POTENTIAL
&END XC
&SCF
SCF_GUESS ATOMIC
EPS_SCF 1.0E-5
MAX_SCF 100
&OT
PRECONDITIONER FULL_ALL
MINIMIZER DIIS
ENERGY_GAP 0.001
&END OT
&OUTER_SCF
MAX_SCF 50
EPS_SCF 1.0E-5
&END
&PRINT
&RESTART
&EACH
MD 1000
&END EACH
&END
&END
&END SCF
&END DFT
&PRINT
&GRID_INFORMATION ON
&END GRID_INFORMATION
&END PRINT
&END FORCE_EVAL
&MOTION
&MD
ENSEMBLE NVT
STEPS 10000
TIMESTEP 0.5
TEMPERATURE 300
&THERMOSTAT
TYPE NOSE
REGION MASSIVE
&NOSE
LENGTH 6
MTS 4
TIMECON 10.00
YOSHIDA 9
&END NOSE
&END THERMOSTAT
&END MD
&PRINT
&RESTART
BACKUP_COPIES 0
&EACH
MD 1000
&END EACH
&END RESTART
&RESTART_HISTORY
&EACH
MD 0
&END EACH
&END RESTART_HISTORY
&TRAJECTORY
FORMAT XYZ
&EACH
MD 1
&END EACH
ADD_LAST NUMERIC
&END TRAJECTORY
&STRESS
&EACH
MD 1
&END EACH
ADD_LAST NUMERIC
&END STRESS
&END PRINT
&END MOTION
&EXT_RESTART
RESTART_FILE_NAME mn_water_nvt.restart
&END
--
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/aff5ccfe-f185-4674-8e7a-3cf4c70878aen%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20221123/e39198eb/attachment.htm>
More information about the CP2K-user
mailing list