当前位置：
&PWscf计算的声子谱的问题求助
PWscf计算的声子谱的问题求助
作者 刘仕晨
我们组是用的vasp，然后计算声子是用的vasp+phonopy，前段时间导师让学习PWscf，然后最近用PWscf计算了下CaTiO3的cubic相声子谱，和vasp计算的进行了对比，发现了一些问题，我也不是很明白，希望大家能帮忙指导下，谢谢！第一张图是用vasp+phonopy计算的，第二张图是用PWscf计算的。
1，CaTiO3用PWscf计算出来的声子谱Gama点居然没有虚频，这个不应该啊，CaTiO3是铁电材料啊，Gama点肯定有虚频的啊，而且vasp的计算结果也是有的，包括我看的一篇文献也是有的啊，J. Chem. Phys., 114, 5, 2001
2，用PWscf计算的声子色散曲线出来的数据需要用plotband.x进行处理，plotband.in文件里要输入费米能级，我觉得声子是不需要像band哪里减去费米能级的吧，而且CTO是绝缘体，占据方法我选的是fixed，自洽结果中也没有这个费米能级，所以我取的是0，不知道这样处理数据对不对，给出in文件：
CTO444.freq
CTOphon.dat
CTOphon.eps
3，最后是画图问题，有点纠结了，就是plotband.x处理出来的数据dat文件，导入xmgrace画图的时候不能像band那样给出高对称点的位置，所以我只能在大致的位置给出了高对称点，我就想知道，PWscf这样处理的数据，怎么才能导入xmgrace，就能在图中给出高对称点位置的那条平行Y轴的线，谢谢！
[ Last edited by 刘仕晨 on
at 11:44 ]
请教楼主一个问题：为什么铁电材料在Gamma点就必须要有虚频呢？
引用回帖:: Originally posted by future_wl at
请教楼主一个问题：为什么铁电材料在Gamma点就必须要有虚频呢？ 我也刚进组刚开始学习铁电材料,也不是很懂,但是铁电材料在非稳定相下计算声子谱,Gamma点是会出现虚频的，这是探讨其稳定相的一个方法，就是通过观察分析声子谱。当然如果有理解错的地方，麻烦指导，谢谢。
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。
引用回帖:: Originally posted by gemucai at
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。 好的,谢谢,我试试,我本来以为是势选的不好,换了个势正在算,算我我就按您的方法试试.那我处理数据的方法有问题么,就是处理数据的时候把费米能级写0这样处理,谢谢您的指导!
引用回帖:: Originally posted by gemucai at
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。 看了一下LZ说的VASP计算的结果（如果在Gamma处有虚频的是的话），三个声学支还是强制为零的，虚频其实是光学支～也就是加上ASR也问题不大才对～
引用回帖:: Originally posted by souledge at
看了一下LZ说的VASP计算的结果（如果在Gamma处有虚频的是的话），三个声学支还是强制为零的，虚频其实是光学支～也就是加上ASR也问题不大才对～... 那请问,我的问题是在哪里啊,谢谢!
引用回帖:: Originally posted by gemucai at
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。 我按照您的方式尝试了下，还是没有出现虚频，声子谱相差不大，不过Gamma点最低的能量离开了0，附上第一张图。同时我之前还换了个势重新计算了声子谱，没有做什么改变，就是重新做了test的ecut，换了ecut和收敛精度，结果计算下来虚频有了，但是LO-TO spliting却没有处理，这个就奇怪了，不是自动处理的么，而且我看ph.out文件里也已经计算出来了波恩有效电荷和介电张量啊。附上第二张图。
pz-hgh.PNG
24小时热帖
下载小木虫APP
与700万科研达人随时交流当前位置：
&【求助】计算声子谱时遇到的问题
【求助】计算声子谱时遇到的问题
作者 imunbg
1. 用Pwscf计算声子谱时出现了如下提示，不知何意：& &
& & warning: symmetry operation #&&7 not allowed.& &fractional translation:
& && &-0...0000000&&in crystal coordinates
2. 计算一系列特殊点的声子谱时，其中某个点没能算，出现的提示是
& &from phq_readin : error #& && &&&59
& & reading xq
&&不知什么原因并不知如何解决。
&&以上2个问题向高手们赐教了！
这个问题你可以去pwscf的官方论坛上找一找.
主要是你选的这个q点不满足对称性导致的.
第二个问题,是你给的q点不对.
声子普 能与材料的局域结构对应吗？ 怎么找出来啊
从没有遇到过这种问题啊
1，到pw_forum上搜索下
2，贴出详细的输入文件给大家看看
第二个问题是不是你的临时文件太大,磁盘占满了.
好的，按照二位的指点去看一下，并谢谢二位了！
回答楼上：磁盘没满。用连续提交的脚本提交的，某个点没算直接跳到下个点了。
引用回帖:Originally posted by limaojlu at
声子普 能与材料的局域结构对应吗？ 怎么找出来啊 声子谱主要还是反映离子之间的振动关系，根结构有一定关系，但是不能反推出结构信，
24小时热帖
下载小木虫APP
与700万科研达人随时交流24小时热门版块排行榜&&&&
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 identation 的 10 个金币
，回帖就立即获得 2 个金币，每人有 1 次机会
(著名写手)
散金: 2603
在线: 496.9小时
虫号: 842155
注册: 性别: GG专业: 凝聚态物性 II ：电子结构
【求助】关于PWSCF介电常数和声子谱的计算
PWSCF的手册里有一句话&In order to have the complete dynamical matrix at q = 0 including the non-analytic terms, you need to calculate effective charges by specifying
option epsil=.true. to ph.x. This is however not possible (because not physical!) for metals
(i.e. any system subject to a broadening).&
也就是说对于要设定自旋的体系,由于SCF计算采用了非零的高斯展宽系数和大于nelec/2的nbnd从而无法在声子谱计算中获得介电常数，那您能不能推荐一个可以计算磁性体系介电常数的从头算免费软件呢？
另外，在声子谱的计算中，pwscf的ph.x计算需要采用小的网格nq1,nq2,nq3，那如果对于原子数较多的体系，为了减少计算时间把nq1,nq2,nq3设定为nq1=1,nq2=1,nq3=1，这样作可以么？计算结果可靠么？
& 猜你喜欢
已经有10人回复
已经有20人回复
已经有3人回复
已经有3人回复
已经有151人回复
已经有185人回复
已经有34人回复
已经有31人回复
已经有162人回复
已经有117人回复
& 本主题相关价值贴推荐，对您同样有帮助:
已经有11人回复
已经有15人回复
已经有4人回复
已经有8人回复
已经有7人回复
已经有10人回复
已经有28人回复
已经有9人回复
已经有9人回复
已经有8人回复
已经有10人回复
已经有10人回复
已经有50人回复
已经有13人回复
(著名写手)
散金: 2603
在线: 496.9小时
虫号: 842155
注册: 性别: GG专业: 凝聚态物性 II ：电子结构
引用回帖:Originally posted by goldenfisher at
nq1=1,nq2=1,nq3=1
应该是不行的，除非你计算的是某一点的声子频率。这个跟能带计算是一个道理。 如果只关心布里渊区中心GAMMA点的特征振动频率呢?
只计算gamma点的频率有意义么?
相关版块跳转
第一性原理
我要订阅楼主
的主题更新
小木虫,学术科研互动社区,为中国学术科研免费提供动力
违规贴举报删除请联系客服电话： 邮箱：（全天候） 或者 QQ:
广告投放与宣传请联系 李想 QQ：
QQ：&&邮箱：
Copyright &
MuChong.com, All Rights Reserved. 小木虫 版权所有当前位置：
&PWscf计算的声子谱的问题求助
PWscf计算的声子谱的问题求助
作者 刘仕晨
我们组是用的vasp，然后计算声子是用的vasp+phonopy，前段时间导师让学习PWscf，然后最近用PWscf计算了下CaTiO3的cubic相声子谱，和vasp计算的进行了对比，发现了一些问题，我也不是很明白，希望大家能帮忙指导下，谢谢！第一张图是用vasp+phonopy计算的，第二张图是用PWscf计算的。
1，CaTiO3用PWscf计算出来的声子谱Gama点居然没有虚频，这个不应该啊，CaTiO3是铁电材料啊，Gama点肯定有虚频的啊，而且vasp的计算结果也是有的，包括我看的一篇文献也是有的啊，J. Chem. Phys., 114, 5, 2001
2，用PWscf计算的声子色散曲线出来的数据需要用plotband.x进行处理，plotband.in文件里要输入费米能级，我觉得声子是不需要像band哪里减去费米能级的吧，而且CTO是绝缘体，占据方法我选的是fixed，自洽结果中也没有这个费米能级，所以我取的是0，不知道这样处理数据对不对，给出in文件：
CTO444.freq
CTOphon.dat
CTOphon.eps
3，最后是画图问题，有点纠结了，就是plotband.x处理出来的数据dat文件，导入xmgrace画图的时候不能像band那样给出高对称点的位置，所以我只能在大致的位置给出了高对称点，我就想知道，PWscf这样处理的数据，怎么才能导入xmgrace，就能在图中给出高对称点位置的那条平行Y轴的线，谢谢！
[ Last edited by 刘仕晨 on
at 11:44 ]
请教楼主一个问题：为什么铁电材料在Gamma点就必须要有虚频呢？
引用回帖:: Originally posted by future_wl at
请教楼主一个问题：为什么铁电材料在Gamma点就必须要有虚频呢？ 我也刚进组刚开始学习铁电材料,也不是很懂,但是铁电材料在非稳定相下计算声子谱,Gamma点是会出现虚频的，这是探讨其稳定相的一个方法，就是通过观察分析声子谱。当然如果有理解错的地方，麻烦指导，谢谢。
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。
引用回帖:: Originally posted by gemucai at
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。 好的,谢谢,我试试,我本来以为是势选的不好,换了个势正在算,算我我就按您的方法试试.那我处理数据的方法有问题么,就是处理数据的时候把费米能级写0这样处理,谢谢您的指导!
引用回帖:: Originally posted by gemucai at
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。 看了一下LZ说的VASP计算的结果（如果在Gamma处有虚频的是的话），三个声学支还是强制为零的，虚频其实是光学支～也就是加上ASR也问题不大才对～
引用回帖:: Originally posted by souledge at
看了一下LZ说的VASP计算的结果（如果在Gamma处有虚频的是的话），三个声学支还是强制为零的，虚频其实是光学支～也就是加上ASR也问题不大才对～... 那请问,我的问题是在哪里啊,谢谢!
引用回帖:: Originally posted by gemucai at
暂时不确定是什么原因，你可以将所有输入文件里的asr和zasr都设置为'no'，再算一遍。我怀疑是你加入了声学支求和规则，它是强制让gamma点声学支频率为0的。 我按照您的方式尝试了下，还是没有出现虚频，声子谱相差不大，不过Gamma点最低的能量离开了0，附上第一张图。同时我之前还换了个势重新计算了声子谱，没有做什么改变，就是重新做了test的ecut，换了ecut和收敛精度，结果计算下来虚频有了，但是LO-TO spliting却没有处理，这个就奇怪了，不是自动处理的么，而且我看ph.out文件里也已经计算出来了波恩有效电荷和介电张量啊。附上第二张图。
pz-hgh.PNG
24小时热帖
下载小木虫APP
与700万科研达人随时交流 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
pwscf计算晶体的声子色散曲线和态密度
下载积分：400
内容提示：pwscf计算晶体的声子色散曲线和态密度
文档格式：DOCX|
浏览次数：105|
上传日期： 11:06:27|
文档星级：
全文阅读已结束，如果下载本文需要使用
 400 积分
下载此文档
该用户还上传了这些文档
pwscf计算晶体的声子色散曲线和态密度
关注微信公众号