24小时热门版块排行榜    

查看: 3141  |  回复: 6

gjh123

金虫 (正式写手)

[求助] vasp 扫描势能面

请问怎么用vasp 扫描出如下图所示的势能面啊?
回复此楼

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

VASP

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

stractor

金虫 (著名写手)


zzy870720z(金币+1): 谢谢提示 2011-11-06 10:13:15
这个不是VASP扫描,而需要处理VASP输出文件中的 LOCPOT 文件,一个三维矩阵。我没做个类似的图,不过Matlab和Oringin很好处理。
2楼2011-11-06 02:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjh123

金虫 (正式写手)

引用回帖:
2楼: Originally posted by stractor at 2011-11-06 02:53:05:
这个不是VASP扫描,而需要处理VASP输出文件中的 LOCPOT 文件,一个三维矩阵。我没做个类似的图,不过Matlab和Oringin很好处理。

谢谢您的回复!在文献中,作者是想通过这个势能面来计算氢原子在碳纳米管表面的迁移势垒。对于这个势能面的做法,他是这样描述的:
          The determination of the hydrogen migration barrier on CNT was, however, quite a non-trivial task, because the standard techniques of migration barrier determination did not work well due to the inconvenient geometry of the problem.It turned out methodologically easier to scan the
potential surface of hydrogen atom.For this purpose a surface location of the probe hydrogen atom has been fixed and only H distance from the nanotube surface could relax. In order to exclude the global CNT rotation, a pair of carbon atoms far from the hydrogen atom was fixed as well.
         一直没搞明白他是怎么做的。难道他就是用LOCPOT做的吗?
3楼2011-11-06 08:37:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjh123

金虫 (正式写手)

引用回帖:
2楼: Originally posted by stractor at 2011-11-06 02:53:05:
这个不是VASP扫描,而需要处理VASP输出文件中的 LOCPOT 文件,一个三维矩阵。我没做个类似的图,不过Matlab和Oringin很好处理。

没法送你金币啊?
4楼2011-11-06 08:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)

★ ★
zzy870720z(金币+2): 谢谢指点 2011-11-06 10:13:34
引用回帖:
3楼: Originally posted by gjh123 at 2011-11-06 08:37:50:
谢谢您的回复!在文献中,作者是想通过这个势能面来计算氢原子在碳纳米管表面的迁移势垒。对于这个势能面的做法,他是这样描述的:
          The determination of the hydrogen migration barrier on CNT was ...

这个不是LOCPOT做的,上面的话应该说得很清楚了:固定离H较远的两个C,将H放在管子附近,只优化H离管子的距离(其他方向保持不变),如此要做很多点(!),然后把这些点进行拟合就可以了。
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
5楼2011-11-06 09:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjh123

金虫 (正式写手)

引用回帖:
5楼: Originally posted by cavediger at 2011-11-06 09:02:48:
这个不是LOCPOT做的,上面的话应该说得很清楚了:固定离H较远的两个C,将H放在管子附近,只优化H离管子的距离(其他方向保持不变),如此要做很多点(!),然后把这些点进行拟合就可以了。

谢谢您的指点。我再研究一下。
6楼2011-11-06 14:36:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gjh123

金虫 (正式写手)

今天发邮件请教了一下侯柱锋老师,侯老师更详细的给我讲解了计算扩散势垒的方法。
        在此帖出侯老师的讲解,与大家共同学习一下:要构造整个势能面的计算量非常大,最直接的一个方法就是,你所模拟的多孔材料表面按画棋盘的方式定出一个二维网格点,然后计算H2分子吸附在这些二维网格点的顶位时的吸附能(此时固定H2的其中的一个H原子与格点的xy坐标,进行原子位置的优化)。这样的一组 x  y     吸附能的数据,然后画出等高线图。
        另外一种方式:考虑H2在多孔材料表面的不同位置(根据表面结构的特征,考虑一些不等价的位置以及对称性高的位置)上的吸附,得到最稳定的吸附位以及能量较高(对应着不稳定的吸附位)的吸附位置,然后对这两个吸附位当作扩散的末态和初态,做nudged elastic band (NEB)的模拟来得到扩散势垒。
         
        在此,再一次感谢侯老师的讲解!
7楼2011-11-06 14:39:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gjh123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 321求调剂 +3 大米饭! 2026-03-15 3/150 2026-03-15 17:48 by 哈哈哈哈嘿嘿嘿
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 085600材料与化工 326 求调剂 +5 热爱生活ing 2026-03-09 5/250 2026-03-14 02:39 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 7/350 2026-03-13 17:28 by xujiaoszu
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
信息提示
请填处理意见