24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 wangle_xq 的 5 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wangle_xq

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 317.9
红花: 3
帖子: 78
在线: 25小时
虫号: 559160
注册: 2008-05-16
性别: GG
专业: 凝聚态物性I:结构、力学和

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

最近一段时间正在研究Siesta。编译成功以后就是赝势的生成了。

我的问题就是如何判断Siesta的赝势是否合适。
按照手册上的说法，如果AE(全电子)-PT(测试)小于等于1mRy就是好的势。

从结果上看要如何判断？
我贴出了Cu的势文件生成测试过程，和测试结果。各位帮忙看看这样的结果算不算一个合适的势。

先来介绍我的赝势是如何生成的。
在~/atom/文件夹下运行

CODE:

# make

得到atm程序
进入文件夹~/atom/Tutorial/Cu运行命令

CODE:

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

Cu.3dtm2.inp文件内容

CODE:

#

# 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

建立文件Cu.test.inp(文件内容见附件)运行命令

CODE:

# ../pt.sh Cu.test.inp Cu.3dtm2.vps

#grep '&d' OUT

得到以下输出

CODE:

 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