is there any basis sets available for the In and Sb correspond to their GTH poseudopotentials?
Matt W
MattWa... at gmail.com
Wed Jul 25 14:12:59 UTC 2007
Dear Jiangping,
I'm afraid you'll have to make the basis sets yourself. There has
been some discussion on how to do this on the Berlioz site. Search
for "Basis". I'm certainly no expert on this, so the following is
only some uninformed suggestions, but might get you started.
Basically this repeats Mattias Krack's help from the other forum.
The basis set generator is available on the cp2k homepage. The
following input files gave me basis sets that at least ran, but
there'll be a lot more work to create production quality basis sets.
Hope this gets you going, and maybe one of the experts will add some
extra information,
Matt
p.s. you need to remove the comments for the In input to run.
----------------------------------------------
Input In
&ATOM
NAME In
RPRB 2.43 # confinement radius
PPOT # use pseudo potential
OPTI 1.0D-9 # SCF convergence on gradient
ZVAL 49 # Z nuclear
ZEFF 3 # Z effective with pseudo
LMAX 2 # lmax gaussian function to be optimised
OCCU # occupations
1 2 # 1 s shell with 2 e-
1 1 # 1 p shell with 1 e-
1 0 # 1 d shell with - e-
MIXI 0.2 # mixing parameter
CONV 1.0D-14 # numerical precision, I think
ITER 200 # max iterations
IPPN 100 # something else to do with numerical precision
XCFN PADE # exchange correlation function
&END
# Initial guess at basis functions, 1 for each l, with the same
exponents you get
# family basis functions. In this case 5 primitive gaussians for each
shell.
&BASIS
GAUS
5
2.4648638319 1.3110027756 0.6076850366
0.2532242096
0.1028326656
1 2 3 4 5
5
2.4648638319 1.3110027756 0.6076850366
0.2532242096
0.1028326656
1 2 3 4 5
5
2.4648638319 1.3110027756 0.6076850366
0.2532242096
0.1028326656
1 2 3 4 5
&END
# Pseudo potential in CPMD format
&POTENTIAL
GOEDECKER
3 LMAX
.610000000 RC
1 2.865776757 #C C1
.770601968 3 1.256193598 .153856196 -.067905338 -.
397254989 .175330829 -.278329062 H(s) 11 12 13 22 23 33
.858131595 2 .494458640 .160912630 -.380788783 H(p) 11
12 22
1.088690904 1 .129207782 H(d) 11
&END
-----------------------------------------------------------------------
In output
In DZVP-GTH-PADE-CONFINED
1
5 0 2 5 2 2 2
1.80983616 -0.02361623 -0.05100125 0.00154821
-0.03461126 0.00167386 -0.02856158
0.77135763 0.43048454 0.00739332 0.11866656
0.06787267 0.00186930 0.13649253
0.34424447 -0.10922219 -1.51160799 -0.13272537
-0.54312914 -0.06189631 -0.40435526
0.17775024 -0.68338967 5.59806414 -0.46524215
3.14220095 -0.13985358 2.52924071
0.10195031 -0.46088589 -4.22065408 -0.50517650
-2.74337882 -0.83799060 -2.13026017
--------------------------------------------------------------------------
Sb input
&ATOM
NAME Sb
RPRB 2.17
PPOT
OPTI 1.0D-9
ZVAL 51
ZEFF 5
LMAX 2
OCCU
1 2
1 3
1 0
MIXI 0.2
CONV 1.0D-14
ITER 200
IPPN 100
XCFN PADE
&END
&BASIS
GAUS
5
2.4648638319 1.3110027756 0.6076850366
0.2532242096
0.1028326656
1 2 3 4 5
5
2.4648638319 1.3110027756 0.6076850366
0.2532242096
0.1028326656
1 2 3 4 5
5
2.4648638319 1.3110027756 0.6076850366
0.2532242096
0.1028326656
1 2 3 4 5
&END
&POTENTIAL
GOEDECKER
3 LMAX
.590000000 RC
1 6.680227830 #C C1
.597683560 3 1.951476722 -.014538121 -.191918210 .
037537267 .495530687 -.786630579 H(s) 11 12 13 22 23 33
.672122045 2 .970312849 .197229982 -.466731323 H(p) 11
12 22
.856557288 1 .300102616 H(d) 11
&END
-------------------------------------------------------------------------------
Sb output
Sb DZVP-GTH-PADE-CONFINED
1
5 0 2 5 2 2 2
1.50050365 -0.93613600 -0.12677315 -0.15435553
-0.26627111 -0.01938142 -0.32017878
1.29300062 1.36401128 0.08289019 0.30691973
0.26354506 0.04158988 0.40615025
0.45755778 -0.04181078 -1.85746994 -0.12557015
-0.66496501 -0.05462655 -0.52273226
0.24426580 -0.74611674 5.73522614 -0.49245600
3.12477450 -0.16797157 2.53741909
0.13342565 -0.46718080 -4.03860863 -0.49513308
-2.61030563 -0.82413089 -2.02590232
More information about the CP2K-user
mailing list