24小时热门版块排行榜    

查看: 3234  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 04:45 by 4wMiSEwB6436
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见