[交流] 【求助】Siesta赝势文件的生成和测试判断标准 已有3人参与

# make

# ../pg.sh Cu.3dtm2.inp

#
# Cu
#
   pg      Cu with 3d states in valence (as in TM paper)
        tm2     2.60
n=Cu c=car     
       0.0       0.0       0.0       0.0       0.0       0.0
    5    3
    4    0      1.00      0.00
    4    1      0.00      0.00
    3    2     10.00      0.00
      2.10      2.30      2.10     

#23456789012345678901234567890123456789012345678901234567890

# ../pt.sh Cu.test.inp Cu.3dtm2.vps
#grep '&d' OUT

ATM 3.2.2  9-DEC-09     Cu GS 3d10 4s1                                  &v&d
ATM 3.2.2  9-DEC-09     Cu 3d9 4s1 4p1                                  &v&d
ATM 3.2.2  9-DEC-09     Cu 3d9 4s1                                      &v&d
ATM 3.2.2  9-DEC-09     Cu 3d9 4s2  *** sd transfer                     &v&d
&d total energy differences in series
&d          1        2        3        4
&d  1    0.0000
&d  2    0.5108   0.0000
&d  3    0.8759   0.3651   0.0000
&d  4    0.1495  -0.3613  -0.7264   0.0000
*----- End of series ----* spdfg &d&v
ATM 3.2.2  9-DEC-09     Cu GS 3d10 4s1                                  &v&d
ATM 3.2.2  9-DEC-09     Cu 3d9 4s1 4p1                                  &v&d
ATM 3.2.2  9-DEC-09     Cu 3d9 4s1                                      &v&d
ATM 3.2.2  9-DEC-09     Cu 3d9 4s2  *** sd transfer                     &v&d
&d total energy differences in series
&d          1        2        3        4
&d  1    0.0000
&d  2    0.5193   0.0000
&d  3    0.8844   0.3651   0.0000
&d  4    0.1548  -0.3645  -0.7296   0.0000
*----- End of series ----* spdfg &d&v

Originally posted by maomao1210 at 2009-12-14 18:01:

支持该贴，呵呵，赝势好坏最终是由计算的结果好与坏决定的，呵呵，
当然了，如果你要是依据赝势与全电子势之间的比较，或者全电子波函数与赝势波函数之间的比较的话，那仅仅只能说明该赝势可以重现单个原子的价 ...

Originally posted by fanchen021 at 2009-12-10 15:10:
用生成的赝势就算一个 dimer Cu的键长，跟实验值比一比，误差在20%以内就可以接受了。