24小时热门版块排行榜    

查看: 3589  |  回复: 6

gjh123

金虫 (正式写手)

[求助] 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] filecode,4个jtjc了 +4 ziyangfang 2026-08-19 4/200 2026-08-19 10:23 by JChen2024
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 今天放榜吗? +11 布布和一二 2026-08-19 11/550 2026-08-19 10:06 by xter9665
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +52 医学老男孩 2026-08-13 114/5700 2026-08-19 09:38 by 医学老男孩
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见