24小时热门版块排行榜

返回列表

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

伐木工

木虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1989.9
帖子: 33
在线: 56.6小时
虫号: 803034

[交流] 【求助】SIESTA 中ATOM生成赝势测试后如何看赝势的好坏

在atom生成赝势的说明上有测试的相关例子，测试后查看OUT出现一些信息，ATOM说明中的
$ grep '&v' OUT | grep s
ATM3 12-JUL-02 Si Test -- GS 3s2 3p2
3s 0.0 2.0000 -0.79662742 3.23745215 -17.68692611
.... ....
ATM3 12-JUL-02 Si Test -- GS 3s2 3p2
1s 0.0 2.0000 -0.79938037 0.50556261 -3.74114712
.... ....
请问上面的数据的具体含义，为什么第一组数据前面是3S 而第二组数据前面为1S呢？

回复此楼

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-08 22:25:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

★
zzy870720z(金币+1): 谢谢补充 2011-03-09 11:49:17
伐木工(金币+5): 2011-03-09 14:51:25

有点答非所问，主要是比较第三列，是轨道的能量，这个是最重要的，第四列跟第五列的意义，你应该能从输出里面看一下吧，我忘记了，嘻嘻不好意思。

赞一下(1人)

回复此楼

3楼2011-03-09 09:37:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 5 个回答

minmin_0082003

金虫 (正式写手)

★ ★
zzy870720z(金币+2): 谢谢指教 2011-03-09 11:49:10
伐木工(金币+5): 2011-03-09 14:51:20

All-electron 里面的3s轨道对应的是pseudopotential 的1s 轨道（你已经把1s2s2p当做是内核处理了），所以这个没有错。
判断赝式的好坏，就看他的能量差距了。我记得你测试transferbitlity 的时候可以grep “&d” 会得到一个
&d total energy differences in series
&d       1       2       3       4       5
&d  1 0.0000
&d  2 0.4308 0.0000
&d  3 0.4961 0.0653 0.0000
&d  4 0.9613 0.5305 0.4652 0.0000
&d  5 1.4997 1.0689 1.0036 0.5384 0.0000
*----- End of series ----* spdfg &d&v
&d total energy differences in series
&d       1       2       3       4       5
&d  1 0.0000
&d  2 0.4299 0.0000
&d  3 0.4993 0.0694 0.0000
&d  4 0.9635 0.5336 0.4642 0.0000
&d  5 1.5044 1.0745 1.0051 0.5409 0.0000
这是各个configuration 之间的差值，当赝式与全电子的差距不在10^-3 eV左右的时候。可以认为你的赝式很好。其实这个是其次的。
最好的赝式要在你测试一些小分子时得到你想要的结果才是最重要的

赞一下(1人)

回复此楼

2楼2011-03-09 09:29:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖