| 查看: 1269 | 回复: 0 | ||
jk射手座新虫 (著名写手)
|
[求助]
求教优化大晶格结构，优化出的能量总是非常大怎么办
|
|
如题，最近优化PKT 1:5~1:7的结构，晶胞比较大，调试过好几次参数经常计算出错，而且算出来的能量经常是这样的形式： free energy TOTEN = -16745.72548078 eV free energy TOTEN = -354385.714304 eV free energy TOTEN = -3343.64830070 eV free energy TOTEN = -271897.42515522 eV free energy TOTEN = -72412.16945720 eV free energy TOTEN = -105104.91766061 eV free energy TOTEN = -223893.13581697 eV 按理说，稍大一些的结构也就偶尔会出现这样的情况，之前做过1:1~1:4的结构，都没这样过，但是现在一下都是这样，只有个别会优化出正常的结构，比如下边： free energy TOTEN = -267.661869 eV free energy TOTEN = -112091.39839325 eV free energy TOTEN = -583495.94751343 eV free energy TOTEN = -316312.26362194 eV 所以正常的能量范围应该在-267电子伏附近。。。 我想请教各位前辈这是什么原因造成的（目前通过在网上查找各种原因，我个人猜测可能服务器的vasp5.0安装的时候就有问题，但是那是两年前安装的，当时我根本不在这边，所以都不知道如果真是那样要怎么解决这个问题） 我把我的文件贴出来一下： INCAR SYSTEM = PKT ENCUT=400 PREC=N ISMEAR =0 SIGMA=0.05 EDIFF=1E-3 ISTART=0 ICHARG=2 IBRION=2 ISIF=2 NSW=100 LREAL=Auto（这个是为了解决计算中的一个错误提示：BRMIX: very serious problems :the old and the new charge density differ） KPOINTS Monkhorst Pack 0 Monkhorst Pack 8 8 4 0 0 0 poscar我写的是脚本，变换晶格常数c： #!/bin/bash #for ps in 0.5 10 20 30 50 100 150 200 250 300 350 400 450 500 550 600; do for c in 23.80 23.81 23.82 23.83 23.84 23.85 23.86 23.87 23.88 23.89 23.90 ; do cat > POSCAR <<! data_PKT 1.0 3.90 0.00 0.00 0.00 3.90 0.00 0.00 0.00 $c 1 1 5 5 18 Selective dynamics Direct 0.0 0.0 0.000 F F T 0.5 0.5 0.083 F F T 0.0 0.0 0.167 F F T 0.0 0.0 0.333 F F T 0.0 0.0 0.500 F F T 0.0 0.0 0.667 F F T 0.0 0.0 0.833 F F T 0.5 0.5 0.250 F F T 0.5 0.5 0.417 F F T 0.5 0.5 0.583 F F T 0.5 0.5 0.750 F F T 0.5 0.5 0.917 F F T 0.5 0.0 0.083 F F T 0.0 0.5 0.083 F F T 0.5 0.5 0.000 F F T 0.5 0.0 0.250 F F T 0.0 0.5 0.250 F F T 0.5 0.5 0.167 F F T 0.5 0.0 0.417 F F T 0.0 0.5 0.417 F F T 0.5 0.5 0.333 F F T 0.5 0.0 0.583 F F T 0.0 0.5 0.583 F F T 0.5 0.5 0.500 F F T 0.5 0.0 0.750 F F T 0.0 0.5 0.750 F F T 0.5 0.5 0.677 F F T 0.5 0.0 0.917 F F T 0.0 0.5 0.917 F F T 0.5 0.5 0.833 F F T ! /public/jzh/vasp5 cp CONTCAR CONTCAR$c cp OUTCAR OUTCAR$c grep "TOTEN" OUTCAR | tail -1 >>E.dat done POTCAR用的是GGA型的，之前PAW 和PBE 的（都是vasp自带的赝势包）比较容易出错，GGA算的东西比较稳定。 跪求大神们提提意见，怎么才能得到结果差不多可以接受的优化结构参数，不然瞎猜一个a0 c0 也不是那回事儿呀！！！小弟先谢谢了！ |
回复此楼» 猜你喜欢
27博士申请
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
厉娜
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有15人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
找到一些相关的精华帖子，希望有用哦~
加自旋几何优化晶格常数偏差较大
已经有6人回复
- 优化后晶格常数变化很大原因
已经有10人回复
- VASP做结构优化时用ISIF=3只想对晶格常数中的一个维度放开经行优化
已经有7人回复
- VASP优化几何时，ISIF=3时，如何只改变晶格，而不改变位矢方向.
已经有11人回复
- 无法优化出文献中的晶格常数
已经有7人回复
- 求教各位大侠，优化石墨烯的晶格常数时，IBRION与ISIF是如何设置的呢？
已经有4人回复
- 优化完成，发现晶格常数跟实验的数据相差太大
已经有8人回复
- Castep优化结果的晶格常数最终结果怎么看？
已经有13人回复
- castep晶格优化问题
已经有15人回复
- 求教Mn的晶格结构，MS的metals文件夹里没有
已经有11人回复
- 一次型NiAl结构优化，晶格矢量变化非常大，求助~
已经有5人回复
- 用LDA赝势优化出来的晶格常数一定比实验值小吗？谢谢！
已经有10人回复
- VASP一次性优化锐钛矿TiO2晶格常数时出错！！
已经有6人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助成功】晶格优化的结果很诡异
已经有13人回复
- 【求助成功】晶格优化中，murn.f的使用方法
已经有7人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【求助】晶体结构优化，计算出问题，请大家帮忙！
已经有7人回复
- 【原创】晶格优化中ISIF取值的探究（六方MgO为例）
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
