24小时热门版块排行榜    

查看: 3605  |  回复: 16

Prettyswan

木虫 (小有名气)

[交流] 【求助成功】DOSCAR中的信息不全,没有site-projected DOS 已有9人参与

这是我以前算的一个结果,DOSCAR:
引用回帖:
10  10   1   0
  0.1786177E+02  0.3912680E-09  0.3912680E-09  0.1166745E-08  0.5000000E-15
  1.00000000000000005E-004
  CAR
CaFe2As2                                
     20.00000000    -20.00000000 3501      3.96332741      1.00000000
    -20.000  0.0000E+00  0.0000E+00 #总的态密度信息,下同
    -19.989  0.0000E+00  0.0000E+00
    -19.977  0.0000E+00  0.0000E+00
…… ……
…… ……
  19.977  0.0000E+00  0.6800E+02
     19.989  0.0000E+00  0.6800E+02
     20.000  0.0000E+00  0.6800E+02 #总的态密度信息到此为止
     20.00000000    -20.00000000 3501      3.96332741      1.00000000 #分的态密度信息,下同
    -20.000  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -19.989  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -19.977  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
…… ……
…… ……

这是这次的结果,就只有开始的那三列总的,下面的那些分的态密度就没有。这是INCAR:
引用回帖:
System = CaFe2As2 ThCr2Si2 DOS
PREC = Accurate
ISPIN = 1
ISMEAR = -5 #for DOS,tethraedron method with Blochl corrections
SGIMA = 0.02
ISTART = 0
EDIFF = 10E-5
ICHARG = 2 #SC-calculation
#ISIF = 2 #stress and relaxation
NSW = 0 #for DOS, the steps of ions relaxtions
IBRION = -1 #for DOS, ions no update, not moved

LREAL = Auto #for supercell or largecell
LORBIT = 11
EMIN = -10 #energy-range for DOS
EMAX = 15
NEDOS = 3501 #number of grid kpoints DOS
ENLMDL = -5 #no update charge for 3steps
NELMIN = 2 #min of two steps
NELM = 60

LPLANE = .TRUE. #parallelisation for linux cluster
NPAR = 4 #number of nodes
LSCALU = .FALSE.
NSIM = 4

请问大家支招。

[ Last edited by zzy870720z on 2010-6-5 at 13:11 ]
回复此楼

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

VASP 第一性

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by fanchen021 at 2010-06-03 23:04:34:
NPAR = 1 试试呢

这个NPAR是并行计算的控制参数,在不设置NPAR时,默认为1,我开始算的时候没有设置NAPR不过结果也是一样的,说明书上的建议是对于linux cluster来说是NPAR=number of nodes;而且以前有人这么说过,我也试过,不过确实不行。
6楼2010-06-04 10:05:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fanchen021

铁杆木虫 (正式写手)


cenwanglai(金币+1):呵呵 2010-06-03 23:39:11
Prettyswan(金币+3): 2010-06-05 09:48:57
NPAR = 1 试试呢
2楼2010-06-03 23:04:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher


cenwanglai(金币+1):RWIGS好象也不是非常必须的.因为对于含多种原子以及有真空层的结构,RWIGS的物理意义好象并不是很有意义. 2010-06-03 23:40:52
RWIGS
3楼2010-06-03 23:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu-jw09

铜虫 (初入文坛)

zzy870720z(金币+1):xiexie 2010-06-04 00:08:20
zzy870720z(金币-1):不好意思,昨天明明不是这个解释?为什么要修改呢?呵呵 2010-06-04 10:28:11
不懂,帮顶

[ Last edited by liu-jw09 on 2010-6-4 at 00:04 ]
4楼2010-06-04 00:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by liu-jw09 at 2010-06-04 00:01:28:
不懂,帮顶

[ Last edited by liu-jw09 on 2010-6-4 at 00:04 ]

在LORBIT=11时,RWIGS是不被读取的,所以设置也没有用。除了LORBIT=5时不往DOSCAR中写信息,其他时候都是写的。
引用回帖:
VASP.4.6 behaviour:
integer RWIGS line in INCAR  files written   
0         line required             DOSCAR and PROCAR file   
1         line required             DOSCAR and lm decomposed PROCAR file   
2         line required             DOSCAR and lm decomposed PROCAR file + phase factors   
5         line required             PROOUT file   
10       not read                  DOSCAR and PROCAR file   
11       not read                  DOSCAR and lm decomposed PROCAR file   
12       not read                  DOSCAR and lm decomposed PROCAR file + phase factors   

这是我贴的那个完整的DOSCAR的INCAR:
引用回帖:
System = CaFe2As2
PREC = accurate
ISPIN = 1
ISMEAR = -5
SGIMA = 0.2
ISTART = 0
EDIFF = 10E-4
ICHARG = 2
#ISIF = 2
#NSW = 0
#IBRION = 2
NELMDL = -6
LORBIT = 11
EMIN = -20
EMAX = 20
NEDOS = 3501

[ Last edited by Prettyswan on 2010-6-4 at 10:15 ]
5楼2010-06-04 10:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by 锐利的碎片 at 2010-06-03 23:35:43:
RWIGS

在LORBIT=11的时候,RWIGS是不被读取的,VASP会在POTCAR中寻找相关参数,只有当LORBIT为0、1、2时才被读取。4楼是我在第一次计算时的INCAR的信息,没有问题。
但是我在以前的计算中发现,如果当LORBIT为1和2时,使用doscar.f编译的程序读取分态密度时候会提示PROCAR中的数据错误;但是当11和12时,就不会有这个问题。

[ Last edited by Prettyswan on 2010-6-4 at 10:11 ]
7楼2010-06-04 10:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu-jw09

铜虫 (初入文坛)

昨天看错了,我以为是没有设置NEDOS的问题呢,可是发现楼主的INCAR里是设置了的,本想删帖的,可是不知道怎么删,就只有修改了
引用回帖:
Originally posted by liu-jw09 at 2010-06-04 00:01:28:
不懂,帮顶

[ Last edited by liu-jw09 on 2010-6-4 at 00:04 ]

8楼2010-06-04 12:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by liu-jw09 at 2010-06-04 12:44:59:
昨天看错了,我以为是没有设置NEDOS的问题呢,可是发现楼主的INCAR里是设置了的,本想删帖的,可是不知道怎么删,就只有修改了

如果没有设置的话,会使用默认值NEDOS=300的。
明察、慎思、笃行
9楼2010-06-04 13:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

Prettyswan(金币+7): 2010-06-05 09:49:19
引用回帖:
The parallel version has some restrictions: The site projected DOS is not evaluated in the parallel version in the following cases:
VASP.4.5, NPAR 6 =1 no site projected DOS
VASP.4.6, NPAR 6 =1, LORBIT=0-5 no site projected DOS

引用回帖:
The wave function character is calculated, either by projecting the wavefunctions onto spherical harmonics that are non zero within spheres of a radius RWIGS around each ion (LORIT=1, 2), or using a quick projection scheme relying that works only for the PAW method (LORBIT=10,11,12, see below).

引用回帖:
If the projector augmented wave method is used, LORBIT can also be set to 10, 11 or 12. This alternative setting selects a quick method for the determination of the spd- and site projected wave function character and does not require the specification of a Wigner-Seitz radius in the INCAR file (the RWIGS line is neglected in this case). The method works only for PAW POTCAR files and not for ultrasoft or norm conserving pseudopotentials.

必须要NPAR=1,LORBIT=11或者12,必须要使用PAW方法,采用PAW势。试着改动一下。

[ Last edited by gleerat on 2010-6-4 at 18:32 ]
明察、慎思、笃行
10楼2010-06-04 18:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 重要消息,中午系统在维护 +6 yuleib84 2026-08-18 7/350 2026-08-18 15:02 by sytp
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 5/250 2026-08-18 14:55 by LNP@mRNA
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见