| 查看: 3235 | 回复: 13 | ||||
burningdays铜虫 (小有名气)
|
[求助]
求助：vasp如何计算获得同位素振动频率？ 已有1人参与
|
|||
| 请教各位前辈，vasp是否有像高斯中的freqchk程序这样功能，我用vasp求得振动频率之后，如何根据计算的力常数求得将某个原子替换为重同位素之后的振动频率？ 我看phonopy之中有修改原子质量的参数选项，但是手册中又未提及，望高手指点，万分感谢！ |
回复此楼» 猜你喜欢
时间戳又变了8-15
已经有15人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！如何用Castep 计算材料的光学特性？
已经有3人回复
charleslian
木虫 (小有名气)
- 应助: 21 (小学生)
- 金币: 1725.5
- 散金: 100
- 红花: 7
- 帖子: 83
- 在线: 49.3小时
- 虫号: 1089067
- 注册: 2010-09-04
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
蛋蛋小童鞋: 应助指数+1 2017-02-08 09:18:40
|
可以把想用x的同位素的原子换成另一个元素符号，比如x1，然后建立一个改了质量的x赝势给x1即可 发自小木虫IOS客户端 |
3楼2017-02-08 08:06:24
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10506.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
burningdays: 金币+10, ★★★★★最佳答案 2017-02-10 19:26:36
感谢参与，应助指数 +1
burningdays: 金币+10, ★★★★★最佳答案 2017-02-10 19:26:36
|
MASS This tag is not necessary to use usually, because atomic masses are automatically set from the chemical symbols. Atomic masses of a primitive cell are overwritten by the values specified. The order of atoms in the primitive cell that is defined by PRIMITIVE_AXIS tag can be shown using -v option. It must be noted that this tag does not affect to the symmetry search. For example, when there are six atoms in a primitive cell, MASS is set as follows MASS = 28.085 28.085 16.000 16.000 16.000 16.000 phonopy的manual里面写的很清楚啊 https://atztogo.github.io/phonopy/setting-tags.html#pm |
4楼2017-02-08 11:08:28
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
2楼2017-02-06 10:14:17
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
5楼2017-02-08 12:51:58
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
6楼2017-02-08 12:59:57
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
太感谢了！7楼2017-02-08 13:23:48
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10506.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
8楼2017-02-08 15:21:06
burningdays
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 30.1
- 散金: 10
- 红花: 2
- 帖子: 142
- 在线: 32.6小时
- 虫号: 3437080
- 注册: 2014-09-23
- 专业: 材料物理化学
|
感谢前辈的解答！ 修改赝势文件这个我就是在POTCAR中修改了POMASS,别的都没动： TITEL = PAW_PBE Mg 13Apr2007 LULTRA = F use ultrasoft PP ? IUNSCR = 1 unscreen: 0-lin 1-nonlin 2-no RPACOR = 1.500 partial core radius POMASS = 26.982593; ZVAL = 2.000 mass and valenz RCORE = 2.000 outmost cutoff radius 前辈您看我这样修改是不是对的。 关于使用phonopy的这个，我没太理解前辈说的，我并不需要声子谱，我只要振动频率，就是vasp OUTCAR中输出的那3N个，不过我现在需要得到同位素替换后的振动频率，phonopy带有修改质量的选项，但是我不知道在如何操作来提取这个频率 |
9楼2017-02-08 16:06:44
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10506.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
burningdays10楼2017-02-08 19:28:11
