应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】关于siesta中赝势价电子轨道的问题
111/2返回列表12下一页
查看: 2674  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

abinitio

版主 (著名写手)

[交流] 【求助】关于siesta中赝势价电子轨道的问题 已有2人参与

新学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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

siesta问题 SIESTA

» 本帖已获得的红花（最新10朵）

zyzoscar

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2009-03-11 20:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdhlover

荣誉版主 (职业作家)

小木虫之莲花佛
友情帮顶，本人没用过siesta，请高手指点
赞一下 回复此楼
2楼2009-03-13 20:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
wuli8(金币+3,VIP+0):(*^__^*) 3-14 10:35
我以前也疑惑这个事，后来没管他。我想那个4f0是没有什么用的。我们可以自己把4f0去掉，看看有没有4f0，到底有什么区别。我估计没什么区别。明天我先试试，再跟你汇报。
赞一下(10人) 回复此楼
3楼2009-03-13 20:45:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
非常感谢楼上的两位回复
赞一下 回复此楼
4楼2009-03-14 08:28:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
谢谢老虎大王  ，我试了一下，应该是没有区别。而且说明书上说的也比较诡异，
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
赞一下 回复此楼
5楼2009-03-14 20:43:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
qasd(金币+1,VIP+0):xiexie 3-14 21:28
zdhlover(金币+10,VIP+0):应楼主的要求转10个金币给你 3-20 22:18
哎呀，今天周末，我忘了做这件事了。

应该是用pg把赝势生成之后，看看要不要4f0，所生成的赝势有多大区别。也就是两个psf文件之间的区别。如果它们有一小点区别，再用这两个赝势分别搞一个计算，看看计算结果有多大区别。

你现是分别比较ae与ps的eigenvalues，当然不会有变化。
赞一下(16人) 回复此楼
6楼2009-03-14 20:55:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
引用回帖:
Originally posted by 老虎大王 at 2009-3-14 20:55:
哎呀，今天周末，我忘了做这件事了。

应该是用pg把赝势生成之后，看看要不要4f0，所生成的赝势有多大区别。也就是两个psf文件之间的区别。如果它们有一小点区别，再用这两个赝势分别搞一个计算，看看计算结果有 ...

好的，我自己再比较一下。

再请教个问题：
生成的pseudo.gps中，对波函数所求的对数导数（logarithmic derivatives）和傅立中变换，有何用处？有什么物理意义吗？
赞一下 回复此楼
7楼2009-03-14 21:06:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★
spur(金币+2,VIP+0): 3Q！！谢谢参与！欢迎常来！ 3-15 11:11
引用回帖:
Originally posted by abinitio at 2009-3-14 21:06:


好的，我自己再比较一下。

再请教个问题：
生成的pseudo.gps中，对波函数所求的对数导数（logarithmic derivatives）和傅立中变换，有何用处？有什么物理意义吗？

那些文件都是让你做比较用的吧。不管是什么东西，全电子AE计算和PS计算出来的东西总是重合度越高越好吧，但是我看什么都重合也不大可能。到底重合到什么程度算好，那大概是凭经验的。一般来说，做好了一个势，这些东西都是要“have a look"一下。但是我不是在这方面有很多经验的人，我也只是随便装模作样地have a look。呵呵。自己做出来的势好不好，归根到底还是要经过实际计算的检验。
赞一下(10人) 回复此楼
8楼2009-03-15 11:10:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
谢谢，大王
烦请版主转１０个金币给老虎大王，以示感谢。
赞一下 回复此楼
9楼2009-03-15 14:52:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by 老虎大王 at 2009-03-15 11:10:45:
那些文件都是让你做比较用的吧。不管是什么东西，全电子AE计算和PS计算出来的东西总是重合度越高越好吧，但是我看什么都重合也不大可能。到底重合到什么程度算好，那大概是凭经验的。一般来说，做好了一个势， ...

这两天也想做个赝势，就不知道调试怎么去做，难道就一个个修改截断半径？？
希望高手们对这个问题都发表一下看法。
赞一下(1人) 回复此楼
坠入无边的网~~
10楼2011-03-24 18:45:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 abinitio 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭