【求助】關(guān)于siesta中贗勢價(jià)電子軌道的問題
新學(xué)siesta，正在看網(wǎng)上lyon的教程，有問題想請(qǐng)教一下。
在/Pseudo/atom/Tutorial/Si中的;Si.tm2.inp中，為何要將價(jià)電子軌道定義為4個(gè)？一般來說，Si是3s2 3p2 ，頂多再加上一個(gè)3d軌道，為什么要將4f軌道設(shè)置上呢？價(jià)電子能越過那么多層到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 ]
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
友情幫頂，本人沒用過siesta，請(qǐng)高手指點(diǎn)
我以前也疑惑這個(gè)事，后來沒管他。我想那個(gè)4f0是沒有什么用的。我們可以自己把4f0去掉，看看有沒有4f0，到底有什么區(qū)別。我估計(jì)沒什么區(qū)別。明天我先試試，再跟你匯報(bào)。
非常感謝樓上的兩位回復(fù)
謝謝老虎大王 ，我試了一下，應(yīng)該是沒有區(qū)別。而且說明書上說的也比較詭異，
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，
哎呀，今天周末，我忘了做這件事了。
應(yīng)該是用pg把贗勢生成之后，看看要不要4f0，所生成的贗勢有多大區(qū)別。也就是兩個(gè)psf文件之間的區(qū)別。如果它們有一小點(diǎn)區(qū)別，再用這兩個(gè)贗勢分別搞一個(gè)計(jì)算，看看計(jì)算結(jié)果有多大區(qū)別。
你現(xiàn)是分別比較ae與ps的eigenvalues，當(dāng)然不會(huì)有變化。
好的，我自己再比較一下。
再請(qǐng)教個(gè)問題：
生成的pseudo.gps中，對(duì)波函數(shù)所求的對(duì)數(shù)導(dǎo)數(shù)（logarithmic derivatives）和傅立中變換，有何用處？有什么物理意義嗎？
那些文件都是讓你做比較用的吧。不管是什么東西，全電子AE計(jì)算和PS計(jì)算出來的東西總是重合度越高越好吧，但是我看什么都重合也不大可能。到底重合到什么程度算好，那大概是憑經(jīng)驗(yàn)的。一般來說，做好了一個(gè)勢，這些東西都是要“have a look"一下。但是我不是在這方面有很多經(jīng)驗(yàn)的人，我也只是隨便裝模作樣地have a look。呵呵。自己做出來的勢好不好，歸根到底還是要經(jīng)過實(shí)際計(jì)算的檢驗(yàn)。