24小时热门版块排行榜    

查看: 2496  |  回复: 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的回帖

ufo991

铁虫 (正式写手)



mazuju028(金币+1):谢谢交流 2010-11-28 13:31:59
印象中手册上面有这一部分的介绍
2楼2010-11-28 11:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宝5403

铁虫 (初入文坛)



sunyang1988(金币+1):谢谢交流 2010-11-29 19:19:14
snowmami(金币+10): 2010-12-01 11:18:51
引用回帖:
Originally posted by snowmami at 2010-11-27 10:30:39:
各位高手:
       我用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
代表了d轨道的五个分态,,,可以认为是5个更细的轨道吧。。。。。


应该填7吧,,,从左往右数 1,2,3,4,。。。。
3楼2010-11-29 18:10:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snowmami

新虫 (正式写手)



youzhizhe(金币+1):鼓励交流。 2010-12-03 01:25:31
引用回帖:
Originally posted by 宝5403 at 2010-11-29 18:10:25:

DZ2,DX2Y2,DXY,DXZ,DYZ
代表了d轨道的五个分态,,,可以认为是5个更细的轨道吧。。。。。


应该填7吧,,,从左往右数 1,2,3,4,。。。。

我觉得你说的有道理,试了试,果然,如果输入“7”即对应于tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3 中的"2".算好的图中会自动显示为Zn d,Ni d,所以我觉得你说的是对的,谢谢!
但是有一点有些疑惑我算了一个Zn原子的TDOS,发现它和Znd的PDOS几乎重合,这是不是合理的?是不是可以认为Zn 3d对Zn的总态密度贡献最大?请指教,另外如果可以的话能否认识一下,我们都是用WIEN2k的,以后可以交流,QQ号我用短消息发给你了
4楼2010-12-01 11:25:34
已阅   回复此楼   关注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的回帖

snowmami

新虫 (正式写手)


引用回帖:
Originally posted by gzqdyouxia at 2010-12-02 21:24:38:
运行configure_int_lapw -b total 1 s,p,d 2 s,p 9 tot,s,p,d end就可以自动生成case.int拉,祝好运!

感谢您的回答,有点不太相信平时都是手动设置的,有这么好的事可以自动生成我最头疼的文件?我打算试试再说,请原谅,我是新手,很多东西整不明白,呵呵!请问:上述命令-b后的轨道是如何设置呢?您写的是total 1 s,p,d 2 s,p 9 tot,s,p,d end,看不太明白,请详细告诉我好么?
6楼2010-12-03 09:54:55
已阅   回复此楼   关注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的回帖
帖子真精彩!
已经收录到淘贴专辑《wien2k常见问题集锦
8楼2011-09-06 21:40:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 snowmami 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode,很负责任的告诉大家 +4 爱看书的可乐 2026-08-08 5/250 2026-08-08 19:06 by 超级海虾
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见