| 查看: 1873 | 回复: 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，结果都一样。 各位有没有遇到过类似情况，能否给予什么建议？ |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS优化后的晶格参数变大了
已经有8人回复
- 优化后晶格常数变化很大原因
已经有10人回复
- VASP做结构优化时用ISIF=3只想对晶格常数中的一个维度放开经行优化
已经有7人回复
- vasp 静态计算出错
已经有20人回复
- VASP优化几何时，ISIF=3时，如何只改变晶格，而不改变位矢方向.
已经有11人回复
- vasp计算三角晶系转换成primitive cell优化后，晶格常数和实验值怎么比较？急求
已经有4人回复
- 求教各位大侠，优化石墨烯的晶格常数时，IBRION与ISIF是如何设置的呢？
已经有4人回复
- 用vasp优化氮化镧的晶格常数，为什么优化出来的值与实验值偏差很大
已经有18人回复
- 请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例
已经有7人回复
- VASP结构优化ISIF=3，为什么晶格常数不变
已经有12人回复
- VASP静态自洽计算，能量升高7eV，然后局域磁矩消失，请问是什么问题？
已经有8人回复
- 一次型NiAl结构优化，晶格矢量变化非常大，求助~
已经有5人回复
- 请教我的用vasp优化Ni的晶格常数的时候怎么不收敛啊。内含INCAR
已经有10人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- 【求助】关于VASP计算材料的晶格热导率
已经有8人回复
- vasp计算在进行静态计算时，出错。
已经有3人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- VASP一次性优化锐钛矿TiO2晶格常数时出错！！
已经有6人回复
- 【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？
已经有16人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】晶格常数大了以后 VASP的计算值准确吗
已经有5人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
512771485
金虫 (小有名气)
博士
- 应助: 5 (幼儿园)
- 金币: 1020.2
- 红花: 2
- 帖子: 271
- 在线: 66.7小时
- 虫号: 1895409
- 注册: 2012-07-16
- 性别: GG
- 专业: 金属材料的合金相、相变及
wax54
4楼2015-07-11 18:53:18
wmy8802217
木虫 (正式写手)
- 应助: 89 (初中生)
- 金币: 2877
- 红花: 12
- 帖子: 350
- 在线: 461.9小时
- 虫号: 2935682
- 注册: 2014-01-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2015-07-11 09:06:43
wax54
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 151.2
- 帖子: 86
- 在线: 28.5小时
- 虫号: 507923
- 注册: 2008-02-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2015-07-11 09:11:04
awdawdawdd
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 679
- 散金: 200
- 红花: 6
- 帖子: 111
- 在线: 67小时
- 虫号: 2237233
- 注册: 2013-01-12
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
|
我猜也有可能是lz结构没有relax完全，尤其是如果lz做了体积形状都改变的relaxation（ISIF=3），这样的话pulay stress会造成影响的。貌似manual里面说的是当体积形状改变之后，原本的plane wave basis set的completeness就有问题了，所以会导致计算的stress有误差，从而得到错误的结构，这种误差叫做pulay stress。消除pulay stress的办法有加大Encut （1.3*Enmax）或者作E-V fitting with equation of state。 我觉得如果lz用ISIF=3的话可以尝试重复relax，直到最后relaxation只有一步就完成 ref：http://cms.mpi.univie.ac.at/vasp ... s_Pulay_Stress.html |
» 本帖已获得的红花（最新10朵）
5楼2015-07-15 04:35:39
wax54
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 151.2
- 帖子: 86
- 在线: 28.5小时
- 虫号: 507923
- 注册: 2008-02-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-07-15 07:02:51
wax54
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 151.2
- 帖子: 86
- 在线: 28.5小时
- 虫号: 507923
- 注册: 2008-02-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2015-07-15 07:07:32
wax54
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 151.2
- 帖子: 86
- 在线: 28.5小时
- 虫号: 507923
- 注册: 2008-02-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2015-07-15 07:08:23
wax54
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 151.2
- 帖子: 86
- 在线: 28.5小时
- 虫号: 507923
- 注册: 2008-02-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-07-15 07:08:39
wax54
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 151.2
- 帖子: 86
- 在线: 28.5小时
- 虫号: 507923
- 注册: 2008-02-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2015-07-15 07:08:49
30