应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP-DFPT &phonopy 计算，处理数据时报错
51/1返回列表
查看: 3071  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者15702449008将赠送您 30 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

15702449008

铜虫 (小有名气)

[求助] VASP-DFPT &phonopy 计算，处理数据时报错

本人刚学vasp，纯新手，想用VASP-DFPT & phonopy 方法跑个声子谱，先用石墨烯做个测试，扩了4*4*1的胞，通过vasp计算后得到FORCE_CONSTANTS，
然后运行命令：phonopy --dim="4 4 1" -c POSCAR-unitcell band.conf 就会出现报错，请大佬们帮忙瞅一眼。
Creating primitive cell is failed.
PRIMITIVE_AXIS may be incorrectly specified.
Traceback (most recent call last):
  File "/home/yxmin/software/anaconda3/bin/phonopy", line 4, in <module>
    __import__('pkg_resources').run_script('phonopy==1.13.2.107', 'phonopy')
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/pkg_resources/__init__.py", line 658, in run_script
    self.require(requires)[0].run_script(script_name, ns)
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/pkg_resources/__init__.py", line 1438, in run_script
    exec(code, namespace, namespace)
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/phonopy-1.13.2.107-py3.6-linux-x86_64.egg/EGG-INFO/scripts/phonopy", line 469, in <module>
    log_level=log_level)
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/phonopy-1.13.2.107-py3.6-linux-x86_64.egg/phonopy/api_phonopy.py", line 109, in __init__
    self._build_primitive_cell()
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/phonopy-1.13.2.107-py3.6-linux-x86_64.egg/phonopy/api_phonopy.py", line 1487, in _build_primitive_cell
    self._supercell, trans_mat, self._symprec)
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/phonopy-1.13.2.107-py3.6-linux-x86_64.egg/phonopy/structure/cells.py", line 48, in get_primitive
    return Primitive(supercell, primitive_frame, symprec=symprec)
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/phonopy-1.13.2.107-py3.6-linux-x86_64.egg/phonopy/structure/cells.py", line 301, in __init__
    self._primitive_cell(supercell)
  File "/home/yxmin/software/anaconda3/lib/python3.6/site-packages/phonopy-1.13.2.107-py3.6-linux-x86_64.egg/phonopy/structure/cells.py", line 340, in _primitive_cell
    raise ValueError
ValueError

band.conf文件设置如下：
ATOM_NAME = C
DIM = 4 4 1
PRIMITIVE_AXIS= 0.25 0.0 0.0 0.0 0.25 0.0 0.0 0.0 1.0
BAND = 0.0 0.0 0.0 0.08333333 0.08333333 0.0 0.125 0.0 0.0 0.0 0.0 0.0
FORCE_CONSTANTS = READ

报错信息是由于PRIMITIVE_AXIS的参数设置引起的，如果删除这个参数的话，能出结果，但是画出来的声子谱和文献中的不一样。
VASP-DFPT &amp;phonopy 计算，处理数据时报错
图丑，不忍直视

POSCAR-unitcell：
graphene
   1.00
     2.4673063800000000    0.0000000000000000    0.0000000000000000
    -1.2336531900000000    2.1367500000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000   20.0000000000000000
   C
   2
Direct
0.33333333 0.66666667 0.50000000
0.66666667 0.33333333 0.50000000
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼2018-09-04 15:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天飘雪

版主 (知名作家)

认真做事，踏实做人
引用回帖:
4楼: Originally posted by 15702449008 at 2018-09-05 09:23:46
额我改了路径后就能做出石墨烯的声子谱了（没加PRIMITIVE_AXIS），但是加了这个参数PRIMITIVE_AXIS，仍然会报错
BAND = 0.0 0.0 0.0 0.333 0.333 0.0 0.5 0.0 0.0 0.0 0.0 0.0
...

不要加
回复此楼

» 本帖已获得的红花（最新10朵）

15702449008
http://blog.sciencenet.cn/u/pfliu89
5楼2018-09-05 11:01:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

漫天飘雪

版主 (知名作家)

认真做事，踏实做人
BAND = 0.0 0.0 0.0 0.08333333 0.08333333 0.0 0.125 0.0 0.0 0.0 0.0 0.0
路径都错了
赞一下(1人) 回复此楼
http://blog.sciencenet.cn/u/pfliu89
2楼2018-09-04 15:45:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15702449008

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 漫天飘雪 at 2018-09-04 15:45:15
BAND = 0.0 0.0 0.0 0.08333333 0.08333333 0.0 0.125 0.0 0.0 0.0 0.0 0.0
路径都错了

额确实和MS提供的路径不一样，我是按照以前的帖子来做的，原贴：http://muchong.com/t-8503127-1
赞一下 回复此楼
3楼2018-09-05 08:46:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15702449008

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 漫天飘雪 at 2018-09-04 15:45:15
BAND = 0.0 0.0 0.0 0.08333333 0.08333333 0.0 0.125 0.0 0.0 0.0 0.0 0.0
路径都错了

额我改了路径后就能做出石墨烯的声子谱了（没加PRIMITIVE_AXIS），但是加了这个参数PRIMITIVE_AXIS，仍然会报错
BAND = 0.0 0.0 0.0 0.333 0.333 0.0 0.5 0.0 0.0 0.0 0.0 0.0
VASP-DFPT &amp;phonopy 计算，处理数据时报错-1
赞一下 回复此楼
4楼2018-09-05 09:23:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[电化学] 常用的国产电化学工作站有哪些？ +3 123明湘 2024-05-11 3/150 2024-05-11 18:18 by GFSHARON520
[基金申请] 科研之友老是给我发消息 +4 问君611 2024-05-11 5/250 2024-05-11 16:42 by 问君611
[基金申请] NSFC函评几月几日截止提交评审意见？ +7 瞬息宇宙 2024-05-11 7/350 2024-05-11 16:17 by 南海猪
[基金申请] 河南工业大学过年开学后，就不让用国基金报账，这合法合规吗？ 100+6 知己而已 2024-05-09 25/1250 2024-05-11 14:51 by Fluid小伟
[考博] 矿大本-中科大硕（推免）-24年科研助理求职-方向不限 +5 Kelaizhang 2024-05-08 6/300 2024-05-11 14:34 by 明刃
[硕博家园] 29 岁去读一个双非学校机械工程的博士 还有前途吗 +12 funaizheng 2024-05-09 14/700 2024-05-11 13:43 by 维多利亚123
[找工作] 985博士毕业，前几天通过了河南省科学院下面一个所的面试，值得去吗？ +11 wendao_2016 2024-05-07 14/700 2024-05-11 13:12 by 皮特上尉
[考研] 关于要不要三战 +18 煎饼果子007 2024-05-06 25/1250 2024-05-11 11:51 by 千山绝影
[硕博家园] 科研不是打打杀杀，科研就是人情世故 +10 杞天大圣 2024-05-08 10/500 2024-05-11 09:05 by 贪吃fish
[基金申请] E05送审了么？ +4 wly78528 2024-05-08 4/200 2024-05-10 16:51 by 路遥还有谁
[有机交流] 酰胺合成 +6 Spiralup 2024-05-06 11/550 2024-05-10 10:10 by Spiralup
[论文投稿] 论文投稿这种情况希望大吗？ +3 woshichuxuan 2024-05-08 6/300 2024-05-09 17:27 by topedit
[基金申请] 基金函评会不会跨学部的送啊 +5 cls512 2024-05-09 5/250 2024-05-09 17:08 by cauzcl
[教师之家] 请问有没有想要联合培养学生的老师？ +3 吴小懒 2024-05-08 4/200 2024-05-09 15:35 by ou0551
[教师之家] 校级优秀本科毕业论文，比例多少？看指导教师资历吗？ +9 河西夜郎 2024-05-04 13/650 2024-05-09 12:41 by huixiong0627
[有机交流] 胺醛缩合实验问题求助 20+3 坏蛋不叫鸟 2024-05-05 11/550 2024-05-09 12:14 by 坏蛋不叫鸟
[有机交流] 低沸点物质浓缩 100+3 桃桃PXS 2024-05-04 8/400 2024-05-08 19:41 by 谢大奇
[论文投稿] Optics Express一直卡在Final Decision +4 MissPro 2024-05-05 11/550 2024-05-08 08:21 by MissPro
[有机交流] 长烷基链氰基水解成羧酸 求助 5+4 Sheerio 2024-05-05 10/500 2024-05-07 12:06 by Sheerio
[考博] 求博导收留 +4 黑森林nono派 2024-05-06 4/200 2024-05-06 08:02 by 安塔瓦拉多
信息提示
关闭
请填处理意见
关闭