24小时热门版块排行榜    

查看: 2581  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:46 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:35 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 06:02 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +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
[基金申请] 国合可查了 +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 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见