应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 在用Vasp计算时，遇到的问题
91/1返回列表
查看: 2436  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fy-dream

银虫 (初入文坛)

[求助] 在用Vasp计算时，遇到的问题

在用vasp计算DOS过程中，发现计算很慢，就用vasp命令看了一下，出现了下面的情况：Error reading item “images” from file INCAR
    FORTRAN STOP
不知道这是什么情况？
请指教，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-07-02 16:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)
★
fzx2008: 金币+1, 谢谢指教 2012-07-02 18:07:21
查一下：
1.四个输入文件是否齐全
2.文件名是否大写
3.POSCAR格式是否正确
4.INCAR参数是否合理
没看到你的文件，只能这样猜了
赞一下 回复此楼
我们都老得太快，却聪明得太迟。当我一路跌跌撞撞登上山巅，蓦然发现，本来就有一条笔直的路，直达山顶....
2楼2012-07-02 16:33:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 谢谢指教 2012-07-02 18:07:29
注意四个输入文件文件要在当前工作目录下...
基本上是vasp没有找到任何输入文件
赞一下 回复此楼
3楼2012-07-02 16:51:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fy-dream

银虫 (初入文坛)
引用回帖:
2楼: Originally posted by cavediger at 2012-07-02 16:33:13
查一下：
1.四个输入文件是否齐全
2.文件名是否大写
3.POSCAR格式是否正确
4.INCAR参数是否合理
没看到你的文件，只能这样猜了

四个输入文件齐全
文件名字是大写
POSCAR格式是正确的
INCAR参数设置如下：SYSTEM=TiO2
   ISTART=0
   ICHARG=11
   ENCUT=500
   ISMEAR=-5
   NSW=0
   LWAVE=.FALSE.
   IBRION=-1
   LORBIT=10
   LREAL=A
   PREC=Accurate
不知道是什么原因导致的？请指点。
赞一下 回复此楼
4楼2012-07-02 16:54:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fy-dream

银虫 (初入文坛)
引用回帖:
3楼: Originally posted by mousekingadv at 2012-07-02 16:51:08
注意四个输入文件文件要在当前工作目录下...
基本上是vasp没有找到任何输入文件

在当前工作目录下？这个不太懂   我已经把任务在当前文件夹下提交上去了  而且输出文件时有大小的  
嗯 这个问题该怎么解决了？谢谢
赞一下 回复此楼
5楼2012-07-02 16:57:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)
★
xueht987: 金币+1, 多谢指教！ 2012-07-03 11:53:25
引用回帖:
4楼: Originally posted by fy-dream at 2012-07-02 16:54:35
四个输入文件齐全
文件名字是大写
POSCAR格式是正确的
INCAR参数设置如下：SYSTEM=TiO2
   ISTART=0
   ICHARG=11
   ENCUT=500
   ISMEAR=-5
   NSW=0
   LWAVE=.FALSE.
   IBRION=-1
   LORBIT=1 ...

把ISTART改为1，把前面优化好的WAVECAR、CHGCAR拷贝到当前计算目录下，再试试看
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

fy-dream
我们都老得太快，却聪明得太迟。当我一路跌跌撞撞登上山巅，蓦然发现，本来就有一条笔直的路，直达山顶....
6楼2012-07-03 09:45:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fy-dream

银虫 (初入文坛)
送鲜花一朵
引用回帖:
6楼: Originally posted by cavediger at 2012-07-03 09:45:39
把ISTART改为1，把前面优化好的WAVECAR、CHGCAR拷贝到当前计算目录下，再试试看...

谢谢你耐心的帮助，我已经弄好了 。。。
赞一下 回复此楼
7楼2012-07-03 15:08:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fys131415

金虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
WDD880227: 金币+1, 多谢提示 2012-07-04 10:19:59
计算态密度的话，我给你个INCAR
SYSTEM =ureaf
    ISTART =     0    job   : 0-new  1-cont  2-samecut
    ISPIN  =      1    spin polarized calculation?
  NPAR   =  4
    RWIGS  =  1.323 1.217#这个是原子的波尔半径，从POTCAR里找
    LPLANE =   .TRUE.
    NBANDS = 120#不要取太多，用公式
    ICHARG = 2
#    LORBIT  =11#算态密度，不用
#    PREC=  Accurate #不用精确计算

    ENCUT = 300 eV
    IALGO = 48
    LREAL = .TRUE.
    NELMIN = 2
    NELM   = 100;
    EDIFF  =1E-05
DOS related values:
    ISMEAR = -5 ;   SIGMA  =    0.2
赞一下 回复此楼
8楼2012-07-04 09:49:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fys131415

金虫 (正式写手)

【答案】应助回帖

★ ★
WDD880227: 金币+2, 多谢耐心指导 2012-07-04 10:20:17
SYSTEM =ureaf

Startparameter for this Run:
#   NWRITE =      2;   LPETIM=F    write-flag & timer
    输出内容详细程度,做长时间动力学计算的话,最好选0或1(首末步/每步核运动输出),短时运算用2(缺省值),选3出错时给出说明。
    ISTART =      0    job   : 0-new  1-cont  2-samecut
    决定是否读取WAVECAR文件,0-根据INWAV初始化波函数,1、2-从WAVECAR读取波函数（重定义平面波集/平面波集不变）。
ISPIN  =      1    spin polarized calculation?
    是否进行自旋极化计算，1-no,2-yes,缺省值2.
#   ttISPIN  =      1    spin polarized calculation?
#  MAGMOM = 5 -5 46*0
   
NPAR   =   8
    控制VASP切换到带间并行模式，取值整数，缺省值：总结点数。
#   LAECHG = .TRUE.
RWIGS  =   0.6  1.3  1.3  1.2  1.2
   给出Wigner Seitz半径，dos计算用，缺省值：从POTCAR读取。（第一个为原子半径，第二个为波尔半径，单位 ）用第二个
LPLANE =   .TRUE.
#   LSCALU =   .FALSE.
#   SAXIS  =   0 0 1
    NBANDS = 39
    给出计算中实际的总能带数，取值整数，缺省值：NELECT/2+NIONS/2(非自旋)，0.6*NELECT+NMAG（自旋）。
#   LORBMOM = .TRUE.
#   LSORBIT  = .TRUE.
    ICHARG = 2
    决定如何建立初始电荷密度。0-由初始波函数计算电荷密度，1-从CHGCAR文件读取，2-使用原子电荷密度的叠加。
#   LNONCOLLINEAR = .TRUE.
#   NSIM   =    5
#   LWAVE  = .FALSE.
   决定是否把波函数（或电荷密度）写入外部文件WAVECAR(或CHGCAR和CHG)中，取值：.TRUE/.FALSE,缺省值：.TRUE。
  LORBIT  =11
    和适当的RWIGS一起决定是否创建PROCAR或PROOUT,取值：0/1/2/5/10/11/12,缺省值0（创建DOSCAR和PROCAR）
Electronic Relaxation 1
ENCUT = 500 eV
平面波基组的截断能量ev。取值实数，缺省值：受PREC设置影响，从POTCAR文件中找出相应的ENMAX/ＥＮＭＩＮ值来设置。
ＰＲＥＣ＝Ｌｏｗ／Ｍｅｄｉｕｍ／Ａｃｃｕｒａｔｅ／Ｎｏｒｍａｌ／
   ＥＮＣＵＴ＝ＥＮＭＩＮ/ＥＮＭＡＸ/ＥＮＭＡＸ/ＥＮＭＡＸ/
IALGO = 48
指定主算法（整数选择算法），8-共轭梯度算法,38-Davidson block iteration scheme,48-RMM-DIIS，缺省值38。
LREAL = .TRUE.
决定投射是在实空间还是倒异空间。取值：.TRUE.(实空间)/.FALSE.(倒异空间)，缺省值：.FALSE.  。
#   NELMIN = 2
    电子自洽循环最小次数，缺省值2。
#   BMIX  =  2.0
NELM   = 100;
电子自洽循环最大次数，缺省值60。
#   NELMDL = 5
    弛豫次数。
EDIFF  =1E-05
指定电子自洽循环的全局中断条件，用于控制收敛精度。取值实数，缺省值：1E-4。即使其等于0，NELM步电子自洽循环也会执行。
#   EDIFFG= -1E-2
    指定离子弛豫循环的中断条件，用于控制核运动的收敛精度，取值实数，缺省值：10*EDIFF。其>0,在两个离子步的总自由能之差小于EDIFFG时停止；其<0,在所有力都小于EDIFFG时停止；其=0，在NSW步弛豫后停止。
#  Ionic Relaxation
#   ISYM = 2
    决定是否使用对称性。-1、0不使用，1、2（效率高节省内存的电荷密度对称性）、3（仅考虑力和应力张量的对称性，而电荷密度是非对称的）使用。缺省值：如果使用US-PP，取1；如果使用PAW，取2。
#   NSW    =  300    number of steps for IOM
    最大离子步数，取值整数，缺省值0。
#   NBLOCK = 1;  
    在NBLOCK离子步后对成对相关函数和DOS进行计算，并把离子配置写入XDATCAR文件；缺省值1。
KBLOCK = 20    inner block; outer block
在KBLOCK *NBLOCK步主循环后平均的成对相关函数和DOS被写入PCDAT和DOSCAR。缺省值NSW。
#   IBRION =     2    ionic relax: 0-MD 1-quasi-New 2-CG
    决定离子怎样更新和运动。-1无更新，0-MD，1-RMM-DIIS,2-共轭梯度算法，3-Damped MD,5、6-有限差分，7、8-密度函数扰动理论，缺省值：if NSW=0/1,取-1,else 0。1-3是三种弛豫方法，根据ISIF决定是否固定离子位置、晶包大小和形状，在INCAR中必须设置POTIM；0是标准的ab-initio MD,不受ISIF影响，即不改变晶胞大小和形状。5-8支持Hessian Matrix和phonon frequency的计算以及部分固定的ＭＤ。





#   ISIF   =     2    stress and relaxation
决定是否计算应力张量以及弛豫中晶胞变化的自由度。取值０－６，缺省值：if IBRION=0,则0,else2。

#   SMASS  = 2.0
    控制从头MD中的速度，-3微正则系统（总自由能不变），-2保持初速度不变，-1每NBLOCK步调整速度，来保证动能的连续，>=0，Nose算法模拟正则系统。
#   POTIM = 1.00
    IBRION=0时，给出MD每步步长fs，IBRION=1-3时，给出最小化的度量常量，取值实数，缺省值：IBRION=0无缺省，必须指定。
#   TEBEG  573
#   ttIBRION =      2    ionic relax: 0-MD 1-quasi-New 2-CG
#   ttISIF   =      7    stress and relaxation

#   ttISYM   =      1
#   IWAVPR =      1    prediction:  0-non 1-charg 2-wave 3-comb
#   LCORR  =      T    Harris-correction to forces
#   POTIM  =   0.2    time-step for ion-motion
#   TEBEG  = 0.0 ;   TEEND= 0.0 initial temperature
DOS related values:
ISMEAR = 0 ;   


SIGMA  = 0.2   broadening in eV -4-tet -1-fermi 0-gaus
决定smearing的宽度eV 。取值实数，缺省值0.2 。
#   ttNEDOS = 500
# Electronic Relaxation 2
#   IALGO  = 48    algorithm
# PC-function
#   APACO = 10.0
这是详细过程，INCAR一定不能错，否则全白算，费时费力。
赞一下 回复此楼
9楼2012-07-04 09:51:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fy-dream 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭