24小时热门版块排行榜    

查看: 652  |  回复: 2

byin

铜虫 (正式写手)

[求助] 请问如何从siesta的赝势文件里面看出价电子和半芯电子信息

比如Ti的LDA-CA赝势,用ATOM生成:
===========================
Ti ca rel nc  
ATM3       3-MAY-12 Troullier-Martins                       
3s 2.00r r= 1.30/3p 6.00r r= 1.30/3d 2.00r r= 1.30/4f 0.00r r= 1.98/  
   4  3 1110  0.112670553485E-03  0.125000000000E-01   12.0000000000   
Radial grid follows
  0.141722109707E-05  0.285226864102E-05  0.430536686097E-05  0.577674280646E-05
  0.726662638298E-05  0.877525038787E-05  0.103028505467E-04  0.118496655501E-04
  0.134159370911E-04  0.150019099027E-04  0.166078317965E-04  0.182339537010E-04
  0.198805297011E-04  0.215478170776E-04  0.232360763476E-04  0.249455713050E-04
……
===========================
在siesta的输出文件中,显示价电子为4s2,3d2,半芯电子为3s2,3p6,这条信息怎么从赝势文件中得出?
请高手指教!感激不尽
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

byin

铜虫 (正式写手)

2楼2012-05-05 22:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+3, 谢谢指教 2012-06-04 21:00:54
byin: 金币+10 2012-06-05 14:55:08
首先声明,下面只是我的臆测,仅当讨论而已。
Ti的价电子组态是3d24s2。
而在制作赝势的时候,是按照spdf这样的顺序写所要的电子组态的。所以没有芯态的赝势组态应该类似于这样:
4s 2.00r r= 1.30/4p 0.00r r= 1.30/3d 2.00r r= 1.30/4f 0.00r r= 1.98/
但是你现在要加入芯态,同时组态中不能有重复的角动量,那么就要这样写了
3s 2.00r r= 1.30/3p 6.00r r= 1.30/3d 2.00r r= 1.30/4f 0.00r r= 1.98/
这样写不代表你的4s2上的两个电子没有考虑,同样也考虑进来了。
进而言之,就是首先知道原子的基态电子组态,再去看赝势中给出的电子组态,如果基态的电子组态没有出现在里面的,就一起加上,这样就构成了赝势中考虑的电子。
弘德明志博学笃行
3楼2012-06-04 20:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 byin 的主题更新
信息提示
请填处理意见