查看: 1086  |  回复: 1

yayanoba

新虫 (小有名气)

[求助] WIEN2K LDA+U出错?是不是设置有问题?

第一次用wien2k做LDA+U计算,然后一计算就有错误。望高人指点,我计算的系统的是BaZnIrO,其中Ba和O分别有two inequivalent atoms,Zn和Ir只有一个,一共有30个原子。我现在希望在Ir的d轨道加u的作用。首先我设置的是spin-polar的计算,但是整个系统其实是非磁的,因此运用instgen lapw -ask,我都是选择的non-magnetic。然后再设置了.indm和.inorb文件。以下为文件信息:(原子编号分别为1,2:Ba, 3:Zn, 4:Ir, 5,6:O)。考虑Zn和Ir的d轨道,但是只给Ir加u作用。
.indm
-9.                    
2
3  1  2                     
4  1  2      
0 0
疑问:1,index指的是.struct里面对应的原子index吗?
           2,这个里面需要列上所有的原子吗?还是只用列上我要考虑密度矩阵的原子(主要是d或f电子)
.inorb
  1  1  0                    
PRATT  1.0                  
  4 1 2                        
  1                           
   0.22 0.00   
dayfile 文件:
start         (Thu Oct 17 06:18:35 EDT 2013) with lapw0 (40/99 to go)

    cycle 1         (Thu Oct 17 06:18:35 EDT 2013)         (40/99 to go)

>   lapw0         (06:18:35) 22.4u 0.3s 0:22.82 99.9% 0+0k 0+6360io 0pf+0w
>   lapw1  -up   -orb         (06:18:58) 0.0u 0.0s 0:00.00 0.0% 0+0k 0+16io 0pf+0w
error: command   /home/noba/wien2k/lapw1 uplapw1.def   failed

>   stop error
STDOUT文件
LAPW0 END
forrtl: severe (24): end-of-file during read, unit 7, file /home/noba/WIEN2k/BZI-U/BZI-U.vorbup
Image              PC                Routine            Line        Source            
lapw1              000000000051AEAE  Unknown               Unknown  Unknown
lapw1              0000000000519946  Unknown               Unknown  Unknown
lapw1              00000000004C9602  Unknown               Unknown  Unknown
lapw1              000000000047F86C  Unknown               Unknown  Unknown
lapw1              000000000047ED8C  Unknown               Unknown  Unknown
lapw1              00000000004A9640  Unknown               Unknown  Unknown
lapw1              000000000043DB63  inilpw_                   263  inilpw.f
lapw1              000000000043E74D  MAIN__                     42  lapw1_tmp_.F
lapw1              00000000004048CC  Unknown               Unknown  Unknown
libc.so.6          00002B1C6263BC8D  Unknown               Unknown  Unknown
lapw1              00000000004047C9  Unknown               Unknown  Unknown

>   stop error
望高人指点!
我初步是觉得我的.indm 或者.inorb文件不对.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wilsun7

银虫 (初入文坛)

你选 -dm 了么?
2楼2013-11-19 07:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yayanoba 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 8/400 2026-10-12 20:09 by caixiaoye1
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见