24小时热门版块排行榜    

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

wax54

铜虫 (小有名气)

[求助] vasp优化后的晶格结构做静态计算的时候力非常大,不知道是哪出错了。 已有4人参与

我在Li4Fe2Si2S8基础上去掉一个锂离子 模拟 Li1.5FeSiS4,晶格弛豫的标准是EDIFF=10^-6,EDIFFG=-0.01,弛豫过程中不保存WAVECAR和CHGCAR,并逐步提高精度,过程很顺利。弛豫后将CONTCAR复制成POSCAR,然后去掉INCAR当中的EDIFFG,NSW,IBRION和ISIF,做静态计算,计算后的能量竟然比弛豫后的能量高了足足400 meV,并且力都高得像没有弛豫过一样。
如果保留那四个弛豫有关的参数(在没有wavecar和chgcar的情况下),只进行一个离子步的优化(此时满足力的标准自动跳出循环),按说也是静态计算,这时候的能量就和弛豫时的能量一致。
静态计算时候我试验了ALGO=F和N,结果都一样。
各位有没有遇到过类似情况,能否给予什么建议?
回复此楼

» 猜你喜欢

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

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

wax54

铜虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by awdawdawdd at 2015-07-15 04:35:39
我猜也有可能是lz结构没有relax完全,尤其是如果lz做了体积形状都改变的relaxation(ISIF=3),这样的话pulay stress会造成影响的。貌似manual里面说的是当体积形状改变之后,原本的plane wave basis set的complete ...

谢谢回复
10楼2015-07-15 07:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

wmy8802217

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-13 07:55:26
是不是你弛豫的时候没有开自旋极化?

» 本帖已获得的红花(最新10朵)

凝聚态理论与计算物理
2楼2015-07-11 09:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wax54

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by wmy8802217 at 2015-07-11 09:06:43
是不是你弛豫的时候没有开自旋极化?

ISPIN和LDAU都开了,算了好几年了,头一次发现这样的情况。
3楼2015-07-11 09:11:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

512771485

金虫 (小有名气)

博士

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-13 07:55:35
你这个是静态计算,能量比优化后高出0.4eV,作用在原子上的力又从很小变大了,我没理解错吧?我最近在做强度计算,用力做收敛全面优化后,再静态计算,力确实会重新变大是正常的,这是因为计算时只考虑了能量的收敛,静态计算时原子位置不变并不是代表电子分布也不变。而且力作收敛的精度远远高于能量作收敛,因此我觉得能量升高某种程度上来说也是正常的。纯属个人见解。

» 本帖已获得的红花(最新10朵)

努力做研究,早些出成果
4楼2015-07-11 18:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见