ERRORL2 in pw_grids:pw_grid_setup_internal processor
Svante
svante.... at gmail.com
Fri Feb 26 09:06:54 UTC 2010
Hi!
I'm trying to do a geometry and cell optimization of a periodic NiO-
cluster with GPW. I'm quite new to cp2k, so I've been fiddling around
a bit with the input files, and searching (this group mostly) for
answers when I got stuck. I know NiO is quite hard to model, partly
since it's antiferromagnetic, but since cp2k has support for UKS and
you can manually change the occupation of alpha and beta electrons it
should be possible.
I've pasted my input file below, and I realize many of the parameters
are not the ideal for my situation, but as I said, I'm quite new to
QM-calculations in general and cp2k in particular. I get an error
during the first optimization step, when the multigrid is printing,
and the error seems to have to do with the plan wave grid somehow.
This is what the error looks like:
--------------------------
OPTIMIZATION STEP: 1
--------------------------
-------------------------------------------------------------------------------
---- MULTIGRID
INFO ----
-------------------------------------------------------------------------------
count for grid 1: 24048 cutoff
[a.u.] 140.00
count for grid 2: 10524 cutoff
[a.u.] 46.67
count for grid 3: 5828 cutoff
[a.u.] 15.56
count for grid 4: 1616 cutoff
[a.u.] 5.19
total gridlevel count : 42016
*
*** 19:45:53 ERRORL2 in pw_grids:pw_grid_setup_internal
processor 0 ***
*** err=-300 condition FAILED at line
451 ***
*
===== Routine Calling Stack =====
10 pw_grid_setup
9 pw_env_rebuild
8 qs_env_rebuild_pw_env
7 cp_eval_at
6 linmin_2pnt
5 cg_linmin
4 cp_cg_main
3 geoopt_cg
2 cp_cell_opt
1 CP2K
CP2K| condition FAILED at line 451
CP2K| Abnormal program termination, stopped by process number 0
And here is the input file:
&FORCE_EVAL
METHOD QS
STRESS_TENSOR NUMERICAL
&DFT
BASIS_SET_FILE_NAME /home/tore/cp2k/tests/QS/BASIS_SET
POTENTIAL_FILE_NAME /home/tore/cp2k/tests/QS/POTENTIAL
# &POISSON
# PERIODIC XYZ
# &END POISSON
&MGRID
&END MGRID
&QS
METHOD GPW
&END QS
&SCF
# ADDED_MOS 20
EPS_SCF 1.0E-5
MAX_SCF 5
CHOLESKY INVERSE
&DIAGONALIZATION
ALGORITHM STANDARD
&END DIAGONALIZATION
# &PRINT
# &PROGRAM_RUN_INFO MEDIUM
# &END
# &TOTAL_DENSITIES MEDIUM
# &END
# &MO_MAGNITUDE MEDIUM
# &END
# &END PRINT
&END SCF
&XC
&XC_FUNCTIONAL PADE
&END XC_FUNCTIONAL
&END XC
UKS
&END DFT
&SUBSYS
&CELL
ABC 4.178 4.178 4.178
PERIODIC XYZ
&END CELL
&COORD
Ni1 0 0 0
Ni2 2.089 2.089 0
Ni1 0 2.089 2.089
Ni2 2.089 0 2.089
O 2.089 0 0
O 0 2.089 0
O 0 0 2.089
O 2.089 2.089 2.089
&END COORD
&KIND Ni1
ELEMENT Ni
BASIS_SET DZV-GTH-PADE
POTENTIAL GTH-PADE-q18
&BS
&ALPHA
#tar bort en spinn upp-elektron från 4s och 1 från 3d
NEL -1 -1
L 0 2
N 4 3
&END ALPHA
&BETA
#tar bort en spinn ner-elektron från 4s och adderar en 3d
NEL -1 1
L 0 2
N 4 3
&END BETA
&END BS
&END KIND
&KIND Ni2
ELEMENT Ni
BASIS_SET DZV-GTH-PADE
POTENTIAL GTH-PADE-q18
&BS
&ALPHA
#tar bort en spinn upp-elektron från 4s och adderar en 3d
NEL -1 1
L 0 2
N 4 3
&END ALPHA
&BETA
#tar bort en spinn ner-elektron från 4s och 1 från 3d
NEL -1 -1
L 0 2
N 4 3
&END BETA
&END BS
&END KIND
&KIND O
BASIS_SET DZVP-GTH-PADE
POTENTIAL GTH-PADE-q6
&BS
&ALPHA
NEL +1
L 1
N 2
&END ALPHA
&BETA
NEL +1
L 1
N 1
&END BETA
&END BS
&END KIND
# &PRINT
# &CELL MEDIUM
# &END CELL
# &KINDS MEDIUM
# &END KINDS
# &END PRINT
&END SUBSYS
# &PRINT
# &DISTRIBUTION MEDIUM
# &END DISTRIBUTION
# &TOTAL_NUMBERS MEDIUM
# &END TOTAL_NUMBERS
# &END PRINT
&END FORCE_EVAL
&GLOBAL
PRINT_LEVEL MEDIUM
PROJECT_NAME NiOtest10
RUN_TYPE CELL_OPT
EXTENDED_FFT_LENGTHS
&END GLOBAL
&MOTION
&CELL_OPT
TYPE GEO_OPT
OPTIMIZER CG
# KEEP_ANGLES TRUE
&CG
&LINE_SEARCH
TYPE 2PNT
&2PNT
&END 2PNT
&END LINE_SEARCH
&END CG
# &PRINT
# &CELL MEDIUM
# &END CELL
# &PROGRAM_RUN_INFO MEDIUM
# &END PROGRAM_RUN_INFO
# &END PRINT
&END CELL_OPT
&GEO_OPT
MAX_ITER 300
OPTIMIZER LBFGS
# &PRINT
# &PROGRAM_RUN_INFO MEDIUM
# &END PROGRAM_RUN_INFO
# &END PRINT
&END GEO_OPT
&END MOTION
What am I doing wrong here, and what exactly does the errormessage
mean? Does it have to do with the EXTENDED_FFT_LENGHTS if cp2k wasn't
properly linked to that FFT-library (I think it is, but I'm not
totally sure)?
Thanks in advance
/Svante
More information about the CP2K-user
mailing list