24小时热门版块排行榜    

查看: 1886  |  回复: 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 ...

cut off energy满足1.3*Enmax,EDIFFG,NSW,IBRION,ISIF这四个参数 一旦注释掉或者删除就会出现力非常大的情况,如果不注释掉或者删除力就貌似正常了,且一步收敛。只是从没有用过状态方程来计算。
6楼2015-07-15 07:02:51
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 能否申诉? +6 echo8914667 2026-08-30 7/350 2026-08-31 08:41 by kingkocxr
[基金申请] 面上意见出来了 +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
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +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
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +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
信息提示
请填处理意见