24小时热门版块排行榜    

查看: 2496  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了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 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +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
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +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
信息提示
请填处理意见