24小时热门版块排行榜    

查看: 1516  |  回复: 5

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[求助] slab表面小分子反应,如何验算过渡态的虚频?

比如一个6层原子层的slab上表面吸附小分子。slab下部3层原子是固定的。验算虚频的时候,是不是还是固定下面三层原子,然后按照一般频率计算方法来算虚频?这样的话,可以移动的原子数在20~40左右,考虑三个自由度,及其组合,就有很多很多可能了。请问该怎么设置这样的过渡态虚频计算呢?

计算slab吸附小分子的振动时,我是固定除小分子外所有原子的,如果小分子是三原子的,就有3*3*3=27个离子步要跑了,时间上很长的。

求slab表面小分子反应的过渡态虚频计算的INCAR的和POSCAR设置方法。

[ Last edited by cenwanglai on 2011-8-16 at 00:50 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

贺仪

铁杆木虫 (著名写手)

帮顶!没算过虚频……
嘟啦啦嘟啦啦嘟
2楼2011-08-16 08:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

【答案】应助回帖

cenwanglai(金币+4): 谢谢回复! 2011-08-16 16:13:13
引用回帖:
1楼: Originally posted by cenwanglai at 2011-08-16 00:49:22:
比如一个6层原子层的slab上表面吸附小分子。slab下部3层原子是固定的。验算虚频的时候,是不是还是固定下面三层原子,然后按照一般频率计算方法来算虚频?这样的话,可以移动的原子数在20~40左右,考虑三个自由度 ...

就是用IBRION=5算就行了,坐标部分跟优化时一样
3楼2011-08-16 09:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

引用回帖:
3楼: Originally posted by lgxyz at 2011-08-16 09:24:34:
就是用IBRION=5算就行了,坐标部分跟优化时一样

我先试试看。呵呵
4楼2011-08-16 16:21:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

引用回帖:
3楼: Originally posted by lgxyz at 2011-08-16 09:24:34:
就是用IBRION=5算就行了,坐标部分跟优化时一样

确实是如你如说的。

我用一个CeO2晶胞(12个原子)做计算测试,固定所有Ce原子,EDIFF = 1E-5, POTIM = 0.01的时候,有24个频率项,全部都是实频,在3~17THz之间。
5楼2011-08-21 17:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-08-23 14:48:38
引用回帖:
5楼: Originally posted by cenwanglai at 2011-08-21 17:22:24:
确实是如你如说的。

我用一个CeO2晶胞(12个原子)做计算测试,固定所有Ce原子,EDIFF = 1E-5, POTIM = 0.01的时候,有24个频率项,全部都是实频,在3~17THz之间。

有时在算过渡态的虚频也不只是一个,一个大的对应于反应坐标,有时也会有若干小虚频,应该是可以忽略的。本身来讲,vasp的频率分析并不是很准确,毕竟只算一个gamma点。
6楼2011-08-23 10:46:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 01:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:52 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 00:28 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见