Pseudo potential regex example
Geoffrey Wood
ge.f... at gmail.com
Fri Nov 21 19:18:17 UTC 2014
Hello ---
I was running the Pseudo potential input file in the regex examples (O-PBE0-q6.inp)
it runs fine except the final pseudo potential file doesn't seem to have
enough parameters in it (i..e it is missing the last projector). The file
should look something like this:
2 4
0.24455430 2 -16.66721480 2.48731132
2
0.22095592 1 18.33745811
0.21133247 0
but I only get:
2 4 0 0
0.244554300000 2 -16.667214800000 2.487311320000
1
0.220955920000 1 18.337458110000
I ran a number of other atoms (H,C,Cl,Br) and the same thing occurs for
carbon i.e. it is missing the last projector function, for the other atoms
(H,Cl,Br) the files havethe correct numbers in them. I thought that it may
have something to do with multiplicity but if I explicitly set this, say
for carbon:
ELECTRON_CONFIGURATION (3) [He] 2s2 2p2
I'm still missing the last line.
This is the input for C that I used:
&GLOBAL
PROGRAM_NAME ATOM
PRINT_LEVEL DEBUG
&END GLOBAL
&ATOM
ATOMIC_NUMBER 6
ELEMENT C
RUN_TYPE PSEUDOPOTENTIAL_OPTIMIZATION
CORE [He]
ELECTRON_CONFIGURATION (3) [He] 2s2 2p2
MAX_ANGULAR_MOMENTUM 3
COULOMB_INTEGRALS ANALYTIC
EXCHANGE_INTEGRALS ANALYTIC
&METHOD
METHOD_TYPE KOHN-SHAM
RELATIVISTIC DKH(2)
&XC
&XC_FUNCTIONAL BLYP
&END XC_FUNCTIONAL
&END XC
&END METHOD
&OPTIMIZATION
EPS_SCF 1.e-10
&END
&AE_BASIS
BASIS_TYPE GEOMETRICAL_GTO
&END AE_BASIS
&PP_BASIS
BASIS_TYPE GEOMETRICAL_GTO
&END PP_BASIS
&POTENTIAL
POTENTIAL_FILE_NAME C_blyp_gth_q1
POTENTIAL_NAME C
PSEUDO_TYPE GTH
CONFINEMENT 0.0 4.0 2.0
>H_POTENTIAL
2 2
0.33806609 2 -9.13626871 1.42925956
2
0.30232223 1 9.66551228
0.28637912 0
&END
&END POTENTIAL
&POWELL
ACCURACY 1.e-6
STEP_SIZE 0.005
MAX_FUN 100
&END
&END ATOM
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20141121/0d0c0be2/attachment.htm>
More information about the CP2K-user
mailing list