24小时热门版块排行榜    

查看: 2174  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yuyangzeng

新虫 (初入文坛)

[求助] VASP新手求助~不稳定的结构,结构优化后CONTCAR不变 已有3人参与

在数据库里找到材料的晶体结构,使用MS建模后,故意把Z轴方向压缩到了90%以及拉伸到了120%,按理说处于不稳定结构后,计算出来的CONTCAR会和原来的POSCAR有一定区别,但是计算后我发现CONTCAR与POSCAR一模一样,希望大家可以帮忙看看问题出在哪儿?谢谢。!!

INCAR
=============================
SYSTEM=Test
#Starting parameters
ISTART=0
#Electronic Relaxation
ENCUT=450
PREC=High
NELMIN=10
EDIFF=1.0E-5
LREAL=.FALSE.
#Ionic Relaxation
EDIFFG=-0.005
NSW=100
IBRION=2
ISIF=2
#Spin-polarized calculations
ISPIN=2
#DOS related values
ISMEAR=1
#Writing files
LWAVE=.FALSE.
LCHARG=.FALSE.
==============================

POSCAR
==============================
New structure
1.0
        7.1199998856         0.0000000000         0.0000000000
        0.0000000000         7.1199998856         0.0000000000
        0.0000000000         0.0000000000         6.4099998474(这个值就是已经压缩过z轴的,原来是7.12)
   Dy   Co
    8   16
Direct
     0.000000000         0.000000000         0.000000000
     0.000000000         0.500000000         0.500000000
     0.500000000         0.500000000         0.000000000
     0.500000000         0.000000000         0.500000000
     0.750000000         0.250000000         0.750000000
     0.250000000         0.250000000         0.250000000
     0.250000000         0.750000000         0.750000000
     0.750000000         0.750000000         0.250000000
     0.625000000         0.625000000         0.625000000
     0.375000000         0.875000000         0.125000000
     0.875000000         0.125000000         0.375000000
     0.125000000         0.375000000         0.875000000
     0.875000000         0.375000000         0.125000000
     0.375000000         0.125000000         0.875000000
     0.125000000         0.875000000         0.375000000
     0.625000000         0.125000000         0.125000000
     0.375000000         0.375000000         0.625000000
     0.875000000         0.625000000         0.875000000
     0.875000000         0.875000000         0.625000000
     0.375000000         0.625000000         0.375000000
     0.125000000         0.625000000         0.125000000
     0.625000000         0.375000000         0.375000000
     0.625000000         0.875000000         0.875000000
     0.125000000         0.125000000         0.625000000
==================================


然后得到的CONTCAR与POSCAR输出数据一样,根本没有一点改变。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
應該是因為結構對稱性很高造成的
一個簡單的做法可嘗試

試著將某一顆原子微小的移動一下
破壞其對稱性
比方把第一顆原子位置       0.000000000         0.000000000         0.000000000
                         改成      0.010000000         0.000000000         0.000000000
再去讓 VASP 做結構收斂即可
4楼2015-04-09 23:23:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫


【答案】应助回帖

感谢参与,应助指数 +1
去vasp官网查询INCAR中ISIF这个参数,其值等于2,优化过程晶胞基矢是不会变的,很正常。
2楼2015-04-09 21:14:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyangzeng

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 未岸亡灵 at 2015-04-09 21:14:27
去vasp官网查询INCAR中ISIF这个参数,其值等于2,优化过程晶胞基矢是不会变的,很正常。

你好,我是希望晶胞基矢不变的,但我认为下面的原子坐标应该是有明显改变的,但是没有一点变化。
3楼2015-04-09 21:32:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fdd096030079

新虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
ISIF=2,晶胞是不会变的,但应该原子坐标会有所变化。楼主检查一下,vasp有没有进行完整的计算流程。还有可能就是在第一个离子步就达到了收敛条件,然后退出,致使离子没有弛豫。
5楼2015-04-10 09:34:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见