24小时热门版块排行榜    

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

yuhaoli

铜虫 (初入文坛)

[求助] 求助,VASP如何计算分析空位型扩散的空位振动熵和有效频率

请问一下VASP怎么计算分析空位振动熵和有效频率。
要算的体系是包含N(N=128)个原子的超晶胞,完整晶格、含单空位的晶格以及鞍点(空位扩散的能垒处)已经找到,然后现在想得到自扩散的系数。然后需要计算一下振动熵以及有效频率。
问题一:在做频率计算的时候是通过弛豫好的构型作为POSCAR,设置INCAR为POTIM=0.015;IBRION=5或6;NFREE=2或4,然后可以得到体系的振动频率吗?
问题二:在计算得到振动频率之后,怎么得到振动熵以及有效频率。空位振动熵的公式中,对于完整晶格要对3N个声子的振动频率求积,对于含单空位的体系也需要对3N-3个声子的频率求积,这是要我把这三四百个数据一个一个的乘起来吗?还是说我对于公式的理解有误差?
问题三:在计算有效频率的时候,需要求初态(含单空位的晶格)和鞍点(空位扩散的能垒处)的振动频率的商。然后对于鞍点,在空位扩散的路径上,应该会有一个声子的振动频率是虚频,应该去掉。我想问的是,在实际计算中,因为我的空位扩散路径是沿111方向,但是我们计算声子的振动频率应该是分别沿x,y,z方向,那么对于鞍点应存在的虚频会有体现吗?或者说在x,y,z方向都有体现?那么应该如何分析呢?
求解答,谢谢。。。。
回复此楼
科学分两种,一种是物理,一种是集邮。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dujunhui

金虫 (正式写手)

另外,问题一中,关于POTIM,VASP中是设置ion-motion步长的,单位是fs,不明白你为什么设置的步长这么小?NFREE是系统的自由度,为什么是2或者4呢?
6楼2014-10-28 15:46:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

yuhaoli

铜虫 (初入文坛)

没有人么
科学分两种,一种是物理,一种是集邮。
2楼2014-09-23 17:21:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhaoli

铜虫 (初入文坛)

3楼2014-09-25 08:33:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhaoli

铜虫 (初入文坛)

4楼2014-09-26 18:21:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:18 by xSrQNSSjcgoX
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:07 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 uKOSCocCb6Fm 2026-11-12 10/500 2026-11-15 02:15 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-14 22:39 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-14 20:11 by KXypaGgBBlw7
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-14 16:09 by KXypaGgBBlw7
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1Nb0IVlmV6v0 2026-11-11 4/200 2026-11-13 11:18 by Equinoxhua
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:13 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:00 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 10:03 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:41 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
信息提示
请填处理意见