【求助】关于siesta中赝势价电子轨道的问题
新学siesta，正在看网上lyon的教程，有问题想请教一下。
在/Pseudo/atom/Tutorial/Si中的;Si.tm2.inp中，为何要将价电子轨道定义为4个？一般来说，Si是3s2 3p2 ，顶多再加上一个3d轨道，为什么要将4f轨道设置上呢？价电子能越过那么多层到4f上吗？
谢谢。
--------------
$ cat Si.tm2.inp
#
# Pseudopotential generation for Silicon
# pg: simple generation
#
pg Silicon
tm2 3.0 # PS flavor, logder R
n=Si c=car # Symbol, XC flavor,{ |r|s}
0.0 0.0 0.0 0.0 0.0 0.0
3 4 # norbs_core, norbs_valence
3 0 2.00 0.00 # 3s2
3 1 2.00 0.00 # 3p2
3 2 0.00 0.00 # 3d0
4 3 0.00 0.00 # 4f0
1.90 1.90 1.90 1.90 0.00 0.00
#
# Last line (above):
# rc(s) rc(p) rc(d) rc(f) rcore_flag rcore
#
#23456789012345678901234567890123456789012345678901234567890
[ Last edited by abinitio on 2009-3-14 at 16:29 ]
返回小木虫查看更多
京公网安备 11010802022153号
友情帮顶，本人没用过siesta，请高手指点
我以前也疑惑这个事，后来没管他。我想那个4f0是没有什么用的。我们可以自己把4f0去掉，看看有没有4f0，到底有什么区别。我估计没什么区别。明天我先试试，再跟你汇报。
非常感谢楼上的两位回复
谢谢老虎大王 ，我试了一下，应该是没有区别。而且说明书上说的也比较诡异，
Even though they are nominally empty in the ground state,we include the 3d and 4f states in order to generate the corresponding pseudopotentials.
呵呵，稀里糊涂，不知自己明白没明白
--------------
$ grep '&v' OUT
ATM 3.2.2 14-MAR-09 Silicon &v&d
3s 0.5 2.0000 -0.79937161 0.00000000 -17.74263363 &v
3p -0.5 0.6667 -0.30807129 0.00000000 -13.66178958 &v
3p 0.5 1.3333 -0.30567134 0.00000000 -13.60785822 &v
3d -0.5 0.0000 0.00000000 0.00000000 -0.27407047 &v
3d 0.5 0.0000 0.00000000 0.00000000 -0.27407047 &v
---------------------------- &v
3s 0.5 2.0000 -0.79936061 0.50555315 -3.74113059 &v
3p -0.5 0.6667 -0.30804995 0.77243805 -3.26356669 &v
3p 0.5 1.3333 -0.30565760 0.76702460 -3.25197500 &v
3d -0.5 0.0000 0.00000000 0.00140505 -0.07847269 &v
3d 0.5 0.0000 0.00000000 0.00140505 -0.07847269 &v
---------------------------- &v
$ grep '&v' OUT
ATM 3.2.2 11-MAR-09 Silicon &v&d
3s 0.5 2.0000 -0.79937161 0.00000000 -17.74263363 &v
3p -0.5 0.6667 -0.30807129 0.00000000 -13.66178958 &v
3p 0.5 1.3333 -0.30567134 0.00000000 -13.60785822 &v
3d -0.5 0.0000 0.00000000 0.00000000 -0.27407047 &v
3d 0.5 0.0000 0.00000000 0.00000000 -0.27407047 &v
4f -0.5 0.0000 0.00000000 0.00000000 -0.26482365 &v
4f 0.5 0.0000 0.00000000 0.00000000 -0.26482365 &v
---------------------------- &v
3s 0.5 2.0000 -0.79936061 0.50555315 -3.74113059 &v
3p -0.5 0.6667 -0.30804995 0.77243805 -3.26356669 &v
3p 0.5 1.3333 -0.30565760 0.76702460 -3.25197500 &v
3d -0.5 0.0000 0.00000000 0.00140505 -0.07847269 &v
3d 0.5 0.0000 0.00000000 0.00140505 -0.07847269 &v
4f -0.5 0.0000 0.00000000 0.00243411 -0.07586534 &v
4f 0.5 0.0000 0.00000000 0.00243411 -0.07586534 &v
---------------------------- &v，
哎呀，今天周末，我忘了做这件事了。
应该是用pg把赝势生成之后，看看要不要4f0，所生成的赝势有多大区别。也就是两个psf文件之间的区别。如果它们有一小点区别，再用这两个赝势分别搞一个计算，看看计算结果有多大区别。
你现是分别比较ae与ps的eigenvalues，当然不会有变化。
好的，我自己再比较一下。
再请教个问题：
生成的pseudo.gps中，对波函数所求的对数导数（logarithmic derivatives）和傅立中变换，有何用处？有什么物理意义吗？
那些文件都是让你做比较用的吧。不管是什么东西，全电子AE计算和PS计算出来的东西总是重合度越高越好吧，但是我看什么都重合也不大可能。到底重合到什么程度算好，那大概是凭经验的。一般来说，做好了一个势，这些东西都是要“have a look"一下。但是我不是在这方面有很多经验的人，我也只是随便装模作样地have a look。呵呵。自己做出来的势好不好，归根到底还是要经过实际计算的检验。