| 查看: 3208 | 回复: 13 | ||||
burningdays铜虫 (小有名气)
|
[求助]
求助：vasp如何计算获得同位素振动频率？ 已有1人参与
|
|||
| 请教各位前辈，vasp是否有像高斯中的freqchk程序这样功能，我用vasp求得振动频率之后，如何根据计算的力常数求得将某个原子替换为重同位素之后的振动频率？ 我看phonopy之中有修改原子质量的参数选项，但是手册中又未提及，望高手指点，万分感谢！ |
回复此楼» 猜你喜欢
重要来源：本周末出结果
已经有10人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
-
Filecode 又变了，巨变
已经有4人回复
-
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！如何用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 (大学生)
- 金币: 10505.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 (大学生)
- 金币: 10505.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 (大学生)
- 金币: 10505.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
burningdays10楼2017-02-08 19:28:11
