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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】关于siesta中赝势价电子轨道的问题
111/2返回列表12下一页
查看: 3152  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
关闭
请填处理意见
关闭