24小时热门版块排行榜

返回列表

jk射手座

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 9882.6
散金: 1955
红花: 2
帖子: 2281
在线: 175.2小时
虫号: 2777009
注册: 2013-11-04
性别: GG
专业: 零件成形制造

[求助] 求教优化大晶格结构，优化出的能量总是非常大怎么办

如题，最近优化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 也不是那回事儿呀！！！小弟先谢谢了！

回复此楼