24小时热门版块排行榜    

查看: 2516  |  回复: 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一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cHIDTWCTtLlG 2026-08-19 3/150 2026-08-19 15:28 by rHxJfXUaiknj
[基金申请] filecode,4个jtjc了 +8 ziyangfang 2026-08-19 9/450 2026-08-19 13:08 by blacksandy
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +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
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 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
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +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
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见