24小时热门版块排行榜    

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

fangyongxinxi

新虫 (正式写手)

[求助] VASP优化结构遇到的问题。

在用vasp 优化一个体系,INCAR 文件如下。优化完后,出现CONTCAR 文件,即优化成功后的文件。重新建一个文件夹,把前面得到的CONTCAR文件改为POSCAR文件,用同样的参数和方法再进行一次优化,得到CONTCAR_2文件。
问题:
1. 理论上这两个CONTCAR文件中的值应该相差很小,但晶体的a b c 的值却相差不较大,而原子的分数坐标没有差别,不知道为什么,或是哪里参数设置出了问题。
2. 自己这样验证的目的是想看看VASP对结构优化的有效性,所以还有个问题,即如何正确的优化一个体系。


system
ISTART = 0
ICHARG = 2
ISPIN =2
VOSKOWN = 1
ENCUT = 302.00 eV

ISMEAR = 0
SIGMA = 0.1
EDIFF = 0.1E-05

NSW = 100
IBRION = 2
ISIF = 3

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

理想路

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by ufo991 at 2011-07-19 09:27:54
楼主你的问题解决了没有呀?我也遇到类似的问题,我用实验数据导入VASP用ISIF=3优化,发现优化出的各参数严重偏离了实验参数呢。用状态方程拟合却不会,但是ISIF=3优化出的能量要低很多很多哦

“我也遇到类似的问题,我用实验数据导入VASP用ISIF=3优化,发现优化出的各参数严重偏离了实验参数呢。”    请问最后您是怎么在VASP里解决这里遇到的优化问题的?谢谢!
10楼2015-09-21 11:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

fangyongxinxi

新虫 (正式写手)

引用回帖:
Originally posted by fangyongxinxi at 2011-06-06 10:39:57:
在用vasp 优化一个体系,INCAR 文件如下。优化完后,出现CONTCAR 文件,即优化成功后的文件。重新建一个文件夹,把前面得到的CONTCAR文件改为POSCAR文件,用同样的参数和方法再进行一次优化,得到CONTCAR_2文件。 ...

这是待优化的结构:
TiNi B2
   1.00000000000000
     5.7516         0.0000000000000000    0.0000000000000000
     0.0000000000000000     4.06762   0.0000000000000000
     0.0000000000000000    0.0000000000000000    4.06762
   Ni   Ti
   4   4
Direct
  0.0000000000000000  0.7500000000000000  0.2500000000000000
  0.5                   0.75                0.25
  0.0000000000000000  0.2500000000000000  0.7500000000000000
  0.5                   0.25                0.75
  0.25000000000000000  0.2500000000000000  0.2500000000000000
  0.75                0.25                0.25
  0.25000000000000000  0.7500000000000000  0.7500000000000000
  0.75                 0.75                0.75

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

这是用上面INCAR文件参数后,得到的CONTCAR
NiTi B2
   1.00000000000000
     5.6274121763945830    0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.9757576568911093    0.0000000000000000
     0.0000000000000000    0.0000000000000000    3.9757576568911093
   Ni   Ti
   4   4
Direct
  0.0000000000000000  0.7500000000000000  0.2500000000000000
  0.5000000000000000  0.7500000000000000  0.2500000000000000
  0.0000000000000000  0.2500000000000000  0.7500000000000000
  0.5000000000000000  0.2500000000000000  0.7500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

将上面第一次优化的CONTCAR 改名为 poscar ,用同样的INCAR, 再次优化得到contcar_2
NiTi B2
   1.00000000000000
     5.5636005577964589    0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.9265497075565161    0.0000000000000000
     0.0000000000000000    0.0000000000000000    3.9265497075565161
   Ni   Ti
   4   4
Direct
  0.0000000000000000  0.7500000000000000  0.2500000000000000
  0.5000000000000000  0.7500000000000000  0.2500000000000000
  0.0000000000000000  0.2500000000000000  0.7500000000000000
  0.5000000000000000  0.2500000000000000  0.7500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

显然两个CONTCAR的结果差别很大:
5.6274121763945830     vs   5.5636005577964589   
悲剧,不知道错在哪里。
2楼2011-06-06 10:45:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
fangyongxinxi(金币+6): 2011-06-06 18:59:44
franch(金币+2): 谢谢回帖交流 2011-06-06 20:05:57
体积比是88.9vs85.8,是有点大,看看两者优化的总能的大小,取小的一个再进行优化,有关优化可以看http://muchong.com/bbs/viewthread.php?tid=3234667&fpage=1
不要直接用isif=3,这样会很不准确
嘟啦啦嘟啦啦嘟
3楼2011-06-06 14:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)

经过再一次验证,悲剧再次出现
将第一次优化的结果,用作第二次优化的结构输入参数,
发现这两次优化的能力和结构参数差别很大。(如上面所说)

but
如果只优化初始结构,无论选什么参数,得到的能力和结构参数都差不多。

是不是可以这样说:
vasp的重复性好。
但如何解释我遇到的问题?
苍天,帮帮我。
4楼2011-06-08 09:40:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
信息提示
请填处理意见