[CP2K-user] First NEB run
Hasan Al-Mahayni
hasanal... at gmail.com
Mon Jun 22 21:05:40 UTC 2020
Hello,
I am trying to test NEB for the first time. I have two slabs (inital.xyz
and final.xyz) and I have included the input file as well (Cp2kNEB). When I
run the job, it creates 239 input and output files called
"slabNEB-r-239.inp/out". The job quickly runs out of memory. What is the
problem?
Thank you,
Hasan.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200622/b9a55bd5/attachment.htm>
-------------- next part --------------
!&EXT_RESTART
!RESTART_FILE_NAME slanNEB-1.restart
!&END
&GLOBAL
PRINT_LEVEL Medium
PROJECT slabNEB
RUN_TYPE BAND # Nudged elastic band calculation
&END GLOBAL
&MOTION
&BAND
NUMBER_OF_REPLICA 5 # Number of "replica" geometries along the path
K_SPRING 0.05
&OPTIMIZE_BAND
OPT_TYPE DIIS
&DIIS
MAX_STEPS 1000
&END
&END
&CONVERGENCE_CONTROL
MAX_FORCE 0.0045
RMS_FORCE 0.003
MAX_DR 0.002
RMS_DR 0.001
&END
BAND_TYPE CI-NEB # Climbing-image NEB
&CI_NEB
NSTEPS_IT 3 # First take 5 normal steps, then start CI
&END
&REPLICA
COORD_FILE_NAME slab1.xyz
&END
&REPLICA
COORD_FILE_NAME slab2.xyz
&END
&PROGRAM_RUN_INFO
INITIAL_CONFIGURATION_INFO
&END
&CONVERGENCE_INFO
&EACH
BAND 1
REPLICA_EVAL 1
&END
&END
&END BAND
&CONSTRAINT
&FIXED_ATOMS
COMPONENTS_TO_FIX XYZ
LIST 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
&END
&END CONSTRAINT
&END MOTION
&FORCE_EVAL
METHOD Quickstep
&DFT
UKS .TRUE. !Requests a spin-polarized calculation using alpha and beta orbitals, i.e. no spin restriction is applied
CHARGE 0 !Total charge
BASIS_SET_FILE_NAME BASIS_MOLOPT
POTENTIAL_FILE_NAME GTH_POTENTIALS
WFN_RESTART_FILE_NAME slabNEB-RESTART.kp
&MGRID
NGRIDS 5
RELATIVE_CUTOFF 50
&END MGRID
&QS
EPS_DEFAULT 1.0E-12
METHOD GPW
EXTRAPOLATION USE_GUESS
&END QS
&SCF
EPS_SCF 1.0E-7
SCF_GUESS RESTART !Use the RESTART file as an initial guess (and ATOMIC if not present).
MAX_SCF 30
ADDED_MOS 200
CHOLESKY OFF
&OUTER_SCF F
EPS_SCF 1.0E-5
MAX_SCF 30
&END OUTER_SCF
&SMEAR T
METHOD FERMI_DIRAC
ELECTRONIC_TEMPERATURE 3.0000000000000000E+02
&END SMEAR
&MIXING T
METHOD BROYDEN_MIXING
ALPHA 4.0000000000000002E-01
BETA 1.5000000000000000E+00
NMIXING 5
NBUFFER 8
&END MIXING
&END SCF
&XC
FUNCTIONAL_ROUTINE NEW
DENSITY_CUTOFF 1.0e-12
GRADIENT_CUTOFF 1.0e-12
TAU_CUTOFF 1.0e-12
&XC_FUNCTIONAL
&END XC_FUNCTIONAL
&XC_GRID
USE_FINER_GRID .TRUE.
&END XC_GRID
&END XC
&POISSON
POISSON_SOLVER PERIODIC
PERIODIC XYZ
&END POISSON
&KPOINTS
SCHEME MONKHORST-PACK 3 3 1
FULL_GRID .TRUE.
&END KPOINTS
&END DFT
&SUBSYS
&CELL
ABC 10.22390 10.22390 30.17390
ALPHA_BETA_GAMMA 90.000000 90.000000 120.000000
PERIODIC XYZ
&END CELL
&KIND Cu
BASIS_SET DZVP-MOLOPT-SR-GTH-q11
POTENTIAL GTH-PBE-q11
&END KIND
&KIND O
BASIS_SET DZVP-MOLOPT-SR-GTH-q6
POTENTIAL GTH-PBE-q6
&END KIND
&KIND C
BASIS_SET DZVP-MOLOPT-SR-GTH-q4
POTENTIAL GTH-PBE-q4
&END KIND
&TOPOLOGY
COORD_FILE_FORMAT XYZ
COORD_FILE_NAME slab1.xyz
CONN_FILE_FORMAT off
&END TOPOLOGY
&END SUBSYS
&END FORCE_EVAL
-------------- next part --------------
A non-text attachment was scrubbed...
Name: slab1.xyz
Type: chemical/x-xyz
Size: 5687 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200622/b9a55bd5/attachment.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: slab2.xyz
Type: chemical/x-xyz
Size: 5685 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200622/b9a55bd5/attachment-0001.xyz>
More information about the CP2K-user
mailing list