24小时热门版块排行榜    

查看: 11095  |  回复: 274
当前主题已经存档。

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

引用回帖:
Originally posted by yufeng1117 at 2009-9-3 19:03:
我搜索了mailing list ,为什么只出现 一下两行文字:You have searched for: argvl:subscript out of range
in the WIEN2k mailing list . Please follow the links below:
然后就没下文了
谢了

这个问题偶也没有碰到过。
不好意思
…………
261楼2009-09-03 19:37:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小胖妹yz

L2main - QTL-B Error


小木虫(金币+0.5):给个红包,谢谢回帖交流
楼主,我做的是个磁性多层膜结构,从下至上依次为NiO,Co,Cu,Co,结构优化已经结束,现在进行具体的计算,但这几天出现了ghost band,我自己看了半天,不知道从何下手,错误信息如下;
WARN:   QTL-B value eq.   4.50  in Band of energy  -0.65219   ATOM=    5   L=  1,You should change the E-parameter in case.in1 or use -in1new switch。
  case.in1文件如下:
       WFFIL        (WFPRI, SUPWF)
  7.00       10    4 (R-MT*K-MAX; MAX L IN WF, V-NMT
.37428   4   0      global e-param with N other choices, napw
0   -0.046     0.000 CONT 1
1   -0.003     0.000 CONT 1
1   -4.332     0.002 STOP 1
2    0.204     0.000 CONT 1
.37428   4   0      global e-param with N other choices, napw
0   -0.188     0.000 CONT 1
1   -0.153     0.000 CONT 1
1   -4.371     0.002 STOP 1
2    0.159     0.000 CONT 1
.37428   4   0      global e-param with N other choices, napw
0   -0.191     0.000 CONT 1
1   -0.149     0.000 CONT 1
1   -4.370     0.002 STOP 1
2    0.162     0.000 CONT 1
.37428   4   0      global e-param with N other choices, napw
0   -0.100     0.000 CONT 1
1   -0.010     0.000 CONT 1
1   -4.368     0.002 STOP 1
2    0.162     0.000 CONT 1
.37428   4   0      global e-param with N other choices, napw
0    0.100     0.000 CONT 1
0   -1.188     0.000 CONT 1
1   -0.130     0.000 CONT 1
2    0.142     0.000 CONT 1
.37428   3   0      global e-param with N other choices, napw
0    0.086     0.000 CONT 1
0   -1.180     0.000 CONT 1
1   -0.116     0.000 CONT 1
.37428   3   0      global e-param with N other choices, napw
0    0.098     0.000 CONT 1
0   -1.189     0.000 CONT 1
1   -0.127     0.000 CONT 1
     费米能为 0.57428
    Ni,o,Co.Cude RMT 分别为1.94, 1.72, 2.07, 2.35
麻烦您了,楼主^_^
262楼2009-09-08 15:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pxjiushun100

金虫 (小有名气)

Structure 问题


小木虫(金币+0.5):给个红包,谢谢回帖交流
我用附件中的structure file,initialize calculation的第一步X nn就出现错误:
Commandline: x nn
Program input is: "2 "

  please specify nn-bondlength factor: (usually=2)
DSTMAX:   20.0000000000000     
iix,iiy,iiz           3           3           2

    ATOM  1  Ba         ATOM  7  O         
RMT(  1)=2.50000 AND RMT(  7)=1.55000
SUMS TO 4.05000  LT.  NN-DIST= 5.23730

    ATOM  2  Y          ATOM  8  O         
RMT(  2)=2.34000 AND RMT(  8)=1.55000
SUMS TO 3.89000  LT.  NN-DIST= 4.52577

    ATOM  3  Cu         ATOM  7  O         
RMT(  3)=1.75000 AND RMT(  7)=1.55000
SUMS TO 3.30000  LT.  NN-DIST= 3.38602

    ATOM  4  Cu         ATOM  8  O         
RMT(  4)=1.75000 AND RMT(  8)=1.55000
SUMS TO 3.30000  LT.  NN-DIST= 3.66385

    ATOM  5  Cu         ATOM  9  O         
RMT(  5)=1.75000 AND RMT(  9)=1.55000
SUMS TO 3.30000  LT.  NN-DIST= 3.38602

    ATOM  6  Cu         ATOM  8  O         
RMT(  6)=1.75000 AND RMT(  8)=1.55000
SUMS TO 3.30000  LT.  NN-DIST= 3.66385

    ATOM  7  O          ATOM  3  Cu        
RMT(  7)=1.55000 AND RMT(  3)=1.75000
SUMS TO 3.30000  LT.  NN-DIST= 3.38602

    ATOM  8  O          ATOM  4  Cu        
RMT(  8)=1.55000 AND RMT(  4)=1.75000
SUMS TO 3.30000  LT.  NN-DIST= 3.66385

    ATOM  9  O          ATOM  5  Cu        
RMT(  9)=1.55000 AND RMT(  5)=1.75000
SUMS TO 3.30000  LT.  NN-DIST= 3.38602
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: Mult not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
WARNING: ityp not equal. PLEASE CHECK outputnn-file
NN created a new CASE.STRUCT_NN FILE
0.016u 0.009s 0:00.04 25.0%        0+0k 0+0io 0pf+0w


我再试着计算YBCO6 antiferromagnetic,所以把原来的cell double了一下
请问在structure文件中什么搞错了?
谢谢
263楼2009-09-25 09:38:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pxjiushun100

金虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuli8(金币+3,VIP+0):3ks 9-26 13:00
如果我去掉  ATOM  9  O,就变得正常了
264楼2009-09-25 10:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huoyang66

铜虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
用wien2k计算能带时,借助xcrysden在布里渊区中取K点,我用的是六角格子,hex,建的晶胞,但是用xcrysden取K点时,case.klist_band中K点坐标出现了如下的样子
    6354756     *********   *********  3  导致计算band时,出现错误

我觉得是K点坐标的有效数位超出了xcrysden的有效数位的限制,所以会出现*号,来代替坐标。这样怎么可以得到有效的klist,来计算能带呢?有没有什么计算k点的方法来借鉴一下呢?

谢谢
265楼2009-10-22 13:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fansw1129

铁杆木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuli8(金币+2,VIP+0):3块肉 very much 11-6 09:41
既然你计算的体系是六角格子问题(hcp)在计算能带动时候就用wien带动hcp就挺好的不要在用crystal产生k点了
引用回帖:
Originally posted by huoyang66 at 2009-10-22 13:33:
用wien2k计算能带时,借助xcrysden在布里渊区中取K点,我用的是六角格子,hex,建的晶胞,但是用xcrysden取K点时,case.klist_band中K点坐标出现了如下的样子
    6354756     *********   *********  3  导致计 ...

266楼2009-10-31 14:09:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fansw1129

铁杆木虫 (著名写手)

★ ★
wuli8(金币+2,VIP+0):3块肉 very much 11-6 09:42
出现这个问题往往是在构造结构的时候
对结构中的原子进行编号后就会出现
其实没有问题
可以继续下面的初始化
引用回帖:
Originally posted by pxjiushun100 at 2009-9-25 09:38:
我用附件中的structure file,initialize calculation的第一步X nn就出现错误:
Commandline: x nn
Program input is: "2 "

  please specify nn-bondlength factor: (usually=2)
DSTMAX:   20. ...

267楼2009-10-31 14:11:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fansw1129

铁杆木虫 (著名写手)

★ ★
wuli8(金币+2,VIP+0):3块肉 very much 11-6 09:42
增加下Ni的Rmt看怎么样?
引用回帖:
Originally posted by 小胖妹yz at 2009-9-8 15:36:
楼主,我做的是个磁性多层膜结构,从下至上依次为NiO,Co,Cu,Co,结构优化已经结束,现在进行具体的计算,但这几天出现了ghost band,我自己看了半天,不知道从何下手,错误信息如下;
WARN:   QTL-B value eq.    ...

268楼2009-10-31 14:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小胖妹yz


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by fansw1129 at 2009-10-31 14:13:
增加下Ni的Rmt看怎么样?


Ok, I'll try it,thasnk you!
269楼2009-11-02 08:59:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huoyang66

铜虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by fansw1129 at 2009-10-31 14:09:
既然你计算的体系是六角格子问题(hcp)在计算能带动时候就用wien带动hcp就挺好的不要在用crystal产生k点了






谢谢啊
计算能带时不是要在布里渊区中选取高对称性K点么.用xcrysden 产生K点很方便。要是用别的方法或是算的话,只是听别人说可以那么做,但是我不会。。
是六角格子( hexagon)不是(hcp)的,就是166号空间群,R-3m,这个格子是不是很特殊啊?
270楼2009-11-02 13:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuli8 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见