24小时热门版块排行榜    

查看: 10327  |  回复: 78
【奖励】 本帖被评价44次,作者lyt幻想泡沫增加金币 35

lyt幻想泡沫

金虫 (小有名气)


[资源] 自编vtstcode编译和CINEB计算教程

自编vtstcode编译和CINEB计算教程正式完成,其中CINEB计算中已将上次的*.xyz文件动画制作包含在内,以供算友们分享。不完善之处还望指正,谢谢!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : CINEB计算方法介绍.pdf
  • 2017-02-23 13:30:02, 1.89 M
  • 附件 2 : VTST编译方法.pdf
  • 2017-02-23 13:30:02, 192.74 K

» 收录本帖的淘帖专辑推荐

精华网帖收集 研究生资料集 VASP and MS 第一性原理相关文档
第一性原理计算-tg vasp 第一性原理 分子模拟 VASP使用

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

obaica

禁言 (著名写手)

本帖内容被屏蔽

3楼2017-02-23 17:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

抖抖豆豆

新虫 (小有名气)


★★★★★ 五星级,优秀推荐

赞一个, 先赞后读
4楼2017-02-23 18:17:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chally123

铁虫 (著名写手)


★★★★★ 五星级,优秀推荐

支持分享
8楼2017-02-24 11:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ak81

新虫 (初入文坛)


★★★★★ 五星级,优秀推荐

哈哈
17楼2017-03-07 12:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


关于gif动画图的制作也可以使用MS和Vedio Mach软件制作,效果完爆VMD

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

22楼2017-04-01 18:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)


送红花一朵
引用回帖:
22楼: Originally posted by lyt幻想泡沫 at 2017-04-01 18:40:43
关于gif动画图的制作也可以使用MS和Vedio Mach软件制作,效果完爆VMD

我按你的方法编译,然后重做你的例子,计算跑完200次迭代就结束了,结果发现01文件夹下得到的CONTCAR和POSCAR完全一样,只是下面多了一堆NaN,OUTCAR里也搜不到ICHAIN等关键词。
这是不是说明我的vtst根本就没有成功安装?
26楼2017-05-18 19:14:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)


送红花一朵
引用回帖:
22楼: Originally posted by lyt幻想泡沫 at 2017-04-01 18:40:43
关于gif动画图的制作也可以使用MS和Vedio Mach软件制作,效果完爆VMD

和vtst官网上的安装方法相比,你的编译方法多了(3)(4)(5)这三步,这个没问题吧?
27楼2017-05-18 19:30:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


引用回帖:
27楼: Originally posted by ra2ghgzh at 2017-05-18 19:30:31
和vtst官网上的安装方法相比,你的编译方法多了(3)(4)(5)这三步,这个没问题吧?...

我编译后计算过,没问题!

发自小木虫Android客户端
29楼2017-05-21 12:00:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)


引用回帖:
29楼: Originally posted by lyt幻想泡沫 at 2017-05-21 12:00:37
我编译后计算过,没问题!
...

“CINEB计算方法介绍.pdf”第十页,也就是附录一.4. CINEB 计算,里面的INCAR中 POTIM=0,这个没问题吗?
30楼2017-05-22 16:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


引用回帖:
30楼: Originally posted by ra2ghgzh at 2017-05-22 16:40:05
“CINEB计算方法介绍.pdf”第十页,也就是附录一.4. CINEB 计算,里面的INCAR中 POTIM=0,这个没问题吗?...

这个参数可能是做教程时弄错了,不为0。

发自小木虫Android客户端
31楼2017-05-22 17:26:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


引用回帖:
26楼: Originally posted by ra2ghgzh at 2017-05-18 19:14:01
我按你的方法编译,然后重做你的例子,计算跑完200次迭代就结束了,结果发现01文件夹下得到的CONTCAR和POSCAR完全一样,只是下面多了一堆NaN,OUTCAR里也搜不到ICHAIN等关键词。
这是不是说明我的vtst根本就 ...

有可能是的。但如果01,02,…,0N文件夹中都生成了OUTCAR文件,说明编译应该没问题。教程里NEB计算的INCAR文件中POTIM有误,应该不为0,看有没有设置错误。

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

32楼2017-05-23 00:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)


送红花一朵
引用回帖:
32楼: Originally posted by lyt幻想泡沫 at 2017-05-23 00:12:37
有可能是的。但如果01,02,…,0N文件夹中都生成了OUTCAR文件,说明编译应该没问题。教程里NEB计算的INCAR文件中POTIM有误,应该不为0,看有没有设置错误。
...

我把POTIM设成0.5又算了一次,这次的INCAR是:
SYSTEM = Pt O
ENCUT = 400
PREC = Normal
ISTART = 0
ICHARG = 2
EDIFF = 1E-5
EDIFFG = -0.05
ISMEAR = 0
SIGMA = 0.10
NELM=60
NSW = 200
ISIF = 2
LCHARG = .FALSE.
LWAVE = .FALSE.
LREAL = Auto
IBRION = 3
POTIM = 0.5
IOPT = 1
ICHAIN = 0
LCLIMB = .TRUE.
SPRING = -5
IMAGES = 1
NPAR = 4
你看这样没问题吧。
然后还是算了200次,01文件夹得到了一个结果,其中的CONTCAR基本就是00的POSCAR,说明过渡态优化到fcc上了。
另外,运行 nebresults.pl 的结果是:

Unziping the OUTCARs ... done
Do nebbarrier.pl ; nebspline.pl
sh: gnuplot: command not found
Do nebef.pl
Do nebmovie.pl
Do nebjmovie.pl
Do nebconverge.pl
mv: cannot stat `vaspout1.eps': No such file or directory
Zipping the OUTCARs again ... done

Forces and Energy:
   0         0.000000      -146.450864         0.000000
   1         0.000000      -146.382045         0.068819
   2         0.000000      -146.029768         0.421096

Extremum 1 found at image  0.007470 with energy: -0.000012
Extremum 2 found at image  1.000000 with energy:  0.068819

没有Extremum 3,这是不是说明hcp没优化好不是极值点?但是明明显示02的Force是0.000000了?
还有sh: gnuplot: command not found是不是说明有个程序没装好?
麻烦帮忙看一下,谢谢。
34楼2017-05-27 09:42:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


引用回帖:
34楼: Originally posted by ra2ghgzh at 2017-05-27 09:42:23
我把POTIM设成0.5又算了一次,这次的INCAR是:
SYSTEM = Pt O
ENCUT = 400
PREC = Normal
ISTART = 0
ICHARG = 2
EDIFF = 1E-5
EDIFFG = -0.05
ISMEAR = 0
SIGMA = 0.10
NELM=60
NSW = 200
ISIF = 2
...

你先将spline.dat导入origin软件作能垒图(教程有介绍),看一下曲线趋势对比数据就能看出是否合理了。那个gunplot.sh应该是没安装吧,但不影响得出结果文件spline.pl,建议你还是用单个命令(nebbarrier.pl和nebspline.pl等)生成结果文件。

发自小木虫Android客户端
36楼2017-05-27 10:14:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


引用回帖:
34楼: Originally posted by ra2ghgzh at 2017-05-27 09:42:23
我把POTIM设成0.5又算了一次,这次的INCAR是:
SYSTEM = Pt O
ENCUT = 400
PREC = Normal
ISTART = 0
ICHARG = 2
EDIFF = 1E-5
EDIFFG = -0.05
ISMEAR = 0
SIGMA = 0.10
NELM=60
NSW = 200
ISIF = 2
...

程序只找出了两个极值点,-0.000012eV的那个极值点00优化的结构能量差别那么小,接近0eV,结构差别可以忽略!

发自小木虫Android客户端
37楼2017-05-27 10:21:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)


引用回帖:
36楼: Originally posted by lyt幻想泡沫 at 2017-05-27 10:14:19
你先将spline.dat导入origin软件作能垒图(教程有介绍),看一下曲线趋势对比数据就能看出是否合理了。那个gunplot.sh应该是没安装吧,但不影响得出结果文件spline.pl,建议你还是用单个命令(nebbarrier.pl和nebspli ...

用origin看过了,不合理。是不是应该加大spring?
38楼2017-05-27 11:00:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)


引用回帖:
37楼: Originally posted by lyt幻想泡沫 at 2017-05-27 10:21:32
程序只找出了两个极值点,-0.000012eV的那个极值点00优化的结构能量差别那么小,接近0eV,结构差别可以忽略!
...

搞定了,把spring加到-10就得到了真正的过渡态。
谢谢您的教程,功德无量啊。
39楼2017-05-27 16:03:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyt幻想泡沫

金虫 (小有名气)


引用回帖:
39楼: Originally posted by ra2ghgzh at 2017-05-27 16:03:09
搞定了,把spring加到-10就得到了真正的过渡态。
谢谢您的教程,功德无量啊。...

共同进步!

发自小木虫Android客户端
40楼2017-05-29 10:13:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张家菜菜子

新虫 (初入文坛)


★★★★★ 五星级,优秀推荐

лл~
42楼2017-08-30 16:39:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

doudouyu198

铁杆木虫 (正式写手)


nebresults.pl

Unziping the OUTCARs ... done
Do nebbarrier.pl ; nebspline.pl
sh: gnuplot: command not found
Do nebef.pl
Do nebmovie.pl
Do nebjmovie.pl
Do nebconverge.pl
mv: cannot stat `vaspout1.eps': No such file or directory
Zipping the OUTCARs again ... done

这个错误是什么原因造成的mv: cannot stat `vaspout1.eps': No such file or directory
再有就是结果里面没有mep.eps这个文件
50楼2018-01-29 12:00:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2017-02-23 17:07   回复  
五星好评  顶一下,感谢分享!
2017-02-23 18:19   回复  
五星好评  顶一下,感谢分享!
mink6楼
2017-02-24 02:36   回复  
五星好评  顶一下,感谢分享!
ZZZ???7楼
2017-02-24 09:46   回复  
五星好评  顶一下,感谢分享!
chally1239楼
2017-02-24 11:37   回复  
顶一下,感谢分享!
chally12310楼
2017-02-24 11:38   回复  
顶一下,感谢分享!
zhangfrank11楼
2017-02-24 13:28   回复  
五星好评  顶一下,感谢分享!
chinaynjs12楼
2017-02-24 14:30   回复  
三星好评  顶一下,感谢分享!
bird00713楼
2017-02-24 14:44   回复  
五星好评  顶一下,感谢分享!
芳苓14楼
2017-02-27 09:43   回复  
gpwang52615楼
2017-02-27 13:42   回复  
五星好评  顶一下,感谢分享!
w47496426016楼
2017-02-28 14:37   回复  
五星好评  顶一下,感谢分享!
常国敬18楼
2017-03-15 16:14   回复  
五星好评  顶一下,感谢分享!
2017-03-18 14:12   回复  
五星好评  顶一下,感谢分享!
2017-03-18 14:18   回复  
顶一下,感谢分享!
2017-03-31 12:39   回复  
五星好评  顶一下,感谢分享!
xn王月23楼
2017-05-11 14:33   回复  
五星好评  顶一下,感谢分享!
ra2ghgzh24楼
2017-05-11 15:37   回复  
五星好评  顶一下,感谢分享!
冯亚娟25楼
2017-05-12 09:03   回复  
五星好评  顶一下,感谢分享!
lj29381528楼
2017-05-19 10:42   回复  
五星好评  顶一下,感谢分享!
fangchaoge33楼
2017-05-25 16:30   回复  
五星好评  顶一下,感谢分享!
2017-05-27 09:52   回复  
五星好评  顶一下,感谢分享!
xieruobing41楼
2017-06-01 16:01   回复  
五星好评  顶一下,感谢分享!
1234bubble43楼
2017-11-29 17:13   回复  
五星好评  顶一下,感谢分享!
2017-11-30 10:43   回复  
五星好评  顶一下,感谢分享!
2017-12-11 17:09   回复  
五星好评  顶一下,感谢分享!
2017-12-13 20:38   回复  
五星好评  顶一下,感谢分享!
cdg123147楼
2017-12-16 12:53   回复  
五星好评  顶一下,感谢分享!
2017-12-26 10:42   回复  
五星好评  顶一下,感谢分享!
2018-01-29 11:31   回复  
顶一下,感谢分享!
相关版块跳转 我要订阅楼主 lyt幻想泡沫 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见