24小时热门版块排行榜    

查看: 2937  |  回复: 13

burningdays

铜虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

charleslian

木虫 (小有名气)

【答案】应助回帖

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

发自小木虫IOS客户端
3楼2017-02-08 08:06:24
已阅   回复此楼   关注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的回帖
普通回帖

burningdays

铜虫 (小有名气)

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

burningdays

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by dxcharlary at 2017-02-08 10:08:28
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. Th ...

感谢您的回复!
我在手册里也看到了这个,但是由于我是刚接触phonopy,对它的功能还不太清楚。
想请教前辈,如果我加上了这个Mass的选项,我最后计算出的频率在哪里提取?我一直没有找到。。
6楼2017-02-08 12:59:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by dxcharlary at 2017-02-08 10:08:28
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. Th ...

原来是dxcharlary 前辈,我之前消虚频就是按着您的建议消掉的太感谢了!
这里我只是需要改变质量后的振动频率,并不需要声子谱,但是我在phonopy里缺找不到该在哪里提取频率,恳请前辈指点,万分感谢!
7楼2017-02-08 13:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by burningdays at 2017-02-08 12:51:58
感谢您的回复!在赝势中修改质量的方法我试过了,在计算中出现这样的提示:BRMIX: very serious problems the old and the new charge density differ
old charge density:    16.00000 new   14.00000(在之前未 ...

这个应该要么是改的什么东西触发了某些关于电子数的bug要么是你没找对地方。
至于下面说phonopy的那个声子谱肯定是可以算的。我觉得你可以尝试一下算一个单胞,然后用phonopy处理,反正你关心的应该只是gamma点处的声子谱。
8楼2017-02-08 15:21:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

burningdays

铜虫 (小有名气)

引用回帖:
8楼: Originally posted by dxcharlary at 2017-02-08 14:21:06
这个应该要么是改的什么东西触发了某些关于电子数的bug要么是你没找对地方。
至于下面说phonopy的那个声子谱肯定是可以算的。我觉得你可以尝试一下算一个单胞,然后用phonopy处理,反正你关心的应该只是gamma点处 ...

感谢前辈的解答!
修改赝势文件这个我就是在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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

引用回帖:
9楼: Originally posted by burningdays at 2017-02-08 16:06:44
感谢前辈的解答!
修改赝势文件这个我就是在POTCAR中修改了POMASS,别的都没动:
    TITEL  = PAW_PBE Mg 13Apr2007
   LULTRA =        F    use ultrasoft PP ?
   IUNSCR =        1    unscreen: 0-lin 1- ...

phonopy可以指定q点,你就指定gamma点就行了

» 本帖已获得的红花(最新10朵)

10楼2017-02-08 19:28:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 burningdays 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 322求调剂 +4 X123- 2026-03-14 4/200 2026-03-21 20:05 by 无际的草原
[考研] 311求调剂 +12 冬十三 2026-03-15 13/650 2026-03-21 19:07 by ColorlessPI
[考研] 本人考085602 化学工程 专硕 +20 不知道叫什么! 2026-03-15 22/1100 2026-03-21 19:03 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +3 石瑞0426 2026-03-19 3/150 2026-03-21 18:59 by RRRKKK
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 7/350 2026-03-21 16:09 by lemonzzn
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见