24小时热门版块排行榜    

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

burningdays

铜虫 (小有名气)

[求助] 求助:vasp如何计算获得同位素振动频率? 已有1人参与

请教各位前辈,vasp是否有像高斯中的freqchk程序这样功能,我用vasp求得振动频率之后,如何根据计算的力常数求得将某个原子替换为重同位素之后的振动频率?  我看phonopy之中有修改原子质量的参数选项,但是手册中又未提及,望高手指点,万分感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

burningdays

铜虫 (小有名气)

另外,有一种方法是通过在POTCAR中修改元素质量,这样需要再重做一遍频率计算,但是,这是修改了体系中所有一类原子的质量,如果只需要修改体系中某一个原子的质量则不能满足。所以,我想请教各位前辈,可不可以通过别的软件来处理获得重同位素替换之后的振动频率。
2楼2017-02-06 10:14:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charleslian

木虫 (小有名气)

【答案】应助回帖

蛋蛋小童鞋: 应助指数+1 2017-02-08 09:18:40
可以把想用x的同位素的原子换成另一个元素符号,比如x1,然后建立一个改了质量的x赝势给x1即可

发自小木虫IOS客户端
3楼2017-02-08 08:06:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by charleslian at 2017-02-08 07:06:24
可以把想用x的同位素的原子换成另一个元素符号,比如x1,然后建立一个改了质量的x赝势给x1即可

感谢您的回复!在赝势中修改质量的方法我试过了,在计算中出现这样的提示:BRMIX: very serious problems the old and the new charge density differ
old charge density:    16.00000 new   14.00000(在之前未修改赝势计算中,是没有这样的提示的),而且最后计算出的频率存在虚频。
按理说只是修改了质量,对别的性质都没有影响,不知您对这一错误有什么看法?
5楼2017-02-08 12:51:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +65 AI8RaGaPaSCR 2026-08-07 68/3400 2026-08-14 06:45 by QNSXhxdpnY4H
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 30/1500 2026-08-14 06:02 by zyjwdhcs
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见