24小时热门版块排行榜    

查看: 2743  |  回复: 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的回帖

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的回帖

dujunhui

金虫 (正式写手)

最近也在看关于空位扩散的文章,感觉确实需要计算的东西很多,和楼主探讨一下。
1.我理解的一般计算振动频率,是需要计算声子谱的,而一般计算声子谱是需要用phonon的,但是好像是收费的。我看的文章上是用的ATAT。
2.三四百个数据一个一个的乘起来肯定是不现实的,应该可以通过自己编写程序实现。
3.计算声子的振动频率应该可以设置沿111方向。
5楼2014-10-28 15:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dujunhui

金虫 (正式写手)

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

lizixiaoxiao

木虫 (初入文坛)

你好!我最近也在计算振动熵,我的邮箱是lkfillz@sina.com,我也有一些问题想和你交流。期待你的回复!熵的公式文献里面有。
7楼2017-10-15 21:36:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuhaoli 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:34 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 06:50 by xSrQNSSjcgoX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 05:11 by xSrQNSSjcgoX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:38 by xSrQNSSjcgoX
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 uKOSCocCb6Fm 2026-11-12 5/250 2026-11-15 03:21 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 9/450 2026-11-15 02:04 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[考博] 27博士申请(市政/环境) +3 ladf 2026-11-12 3/150 2026-11-13 08:19 by Equinoxhua
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:36 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 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一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
信息提示
请填处理意见