24小时热门版块排行榜    

查看: 2055  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

snowmami

新虫 (正式写手)


[交流] 【求助】用WIEN2k算DOS时,如何编辑*.int文件

各位高手:
       我用WIEN2k算2*2*2的ZnO:Ni的超胞,到了DOS这一步,在edit *.int文件时,我想得到TDOS、第一个Zn原子的s、p、d和第二个O的s、p以及位于第九个原子的Ni的tdos、s、p、d的态密度,下面这样设置对么?主要是不明白,DZ2,DX2Y2,DXY,DXZ,DYZ这几个的区别,究竟哪个对应于3d轨道呢?打问号处应该写几才对呢?我是新手,望详细指教!
ATOM  1  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3   
ATOM  2  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3   
ATOM  3  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3   
ATOM  4  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3   
ATOM  5  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3   
ATOM  6  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3
            ..........  
            ..........   
            ..........
ATOM 31  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3   
ATOM 32  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3

ZnNiO
-0.50 0.002 1.500 0.003   EMIN, DE, EMAX, Gauss-broadening(>;de)
    10                      NUMBER OF DOS-CASES specified below
    0    1   total         atom, case=column in qtl-header, label
    1    2   Zn s  
    1    3   Zn p
    1    ?   Zn d
    2    2   O s
    2    3   O p
    9    1   Ni total
    9    2   Ni s
    9    3   Ni p
    9    ?  Ni d
回复此楼

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

wien2k常见问题集锦

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+1):谢谢交流。 2010-12-03 01:25:46
运行configure_int_lapw -b total 1 s,p,d 2 s,p 9 tot,s,p,d end就可以自动生成case.int拉,祝好运!
5楼2010-12-02 21:24:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ellsaking(金币+1):感谢回答~ 2010-12-03 14:55:03
引用回帖:
Originally posted by snowmami at 2010-12-03 09:54:55:

感谢您的回答,有点不太相信平时都是手动设置的,有这么好的事可以自动生成我最头疼的文件?我打算试试再说,请原谅,我是新手,很多东西整不明白,呵呵!请问:上述命令-b后的轨道是如何设置呢?您写的是total ...

当然可以,total表示总态密度,1 s,p,d是说原子1的s,p,d轨道,以后一次类推!
7楼2010-12-03 14:39:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 snowmami 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见