24小时热门版块排行榜    

查看: 2560  |  回复: 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的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫


【答案】应助回帖

感谢参与,应助指数 +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的回帖

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的回帖

fdd096030079

新虫 (小有名气)

【答案】应助回帖

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

yuyangzeng

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by jpchou at 2015-04-09 23:23:06
應該是因為結構對稱性很高造成的
一個簡單的做法可嘗試

試著將某一顆原子微小的移動一下
破壞其對稱性
比方把第一顆原子位置       0.000000000         0.000000000         0.000000000
                  ...

好的我试试,这样子得到的CONTCAR就是稳定的结构了吧?
6楼2015-04-10 18:28:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyangzeng

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by fdd096030079 at 2015-04-10 09:34:33
ISIF=2,晶胞是不会变的,但应该原子坐标会有所变化。楼主检查一下,vasp有没有进行完整的计算流程。还有可能就是在第一个离子步就达到了收敛条件,然后退出,致使离子没有弛豫。

请问如何检查有没有进行完整流程呢?
7楼2015-04-10 18:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuyangzeng 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:11 by DgNGHc3h5tPl
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:35 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:59 by zNcEhTcH7ihq
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:15 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见