24小时热门版块排行榜    

查看: 3073  |  回复: 8
【悬赏金币】回答本帖问题,作者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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15702449008

铜虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by 漫天飘雪 at 2018-09-05 17:53:40
看你是在原胞的基础上计算,还是在超胞的基础上计算。
官网给的例子是NaCl的Conventional Cell, 你用Primitive cell计算的时候不用加这个参数。

好的,谢谢大佬回复,我再研究研究
8楼2018-09-05 18:06:55
已阅   回复此楼   关注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
路径都错了
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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[催化] 齿球形催化剂的尺寸 +3 anndy1971 2024-05-08 5/250 2024-05-11 23:57 by 596699273
[基金申请] 河南工业大学过年开学后,就不让用国基金报账,这合法合规吗? 100+7 知己而已 2024-05-09 27/1350 2024-05-11 22:31 by nono2009
[基金申请] 基金委治打招呼顽疾越治越严重 +18 zzahkj 2024-05-10 26/1300 2024-05-11 21:44 by 星火12
[教师之家] 相比过去,现在高校导师水平越来越高(可见招聘条件),研究生越来越差 +7 苏东坡二世 2024-05-11 10/500 2024-05-11 21:12 by feng6531
[论文投稿] 投稿竟然倒在了Scientific Reports上。 +7 lizhengke06 2024-05-10 7/350 2024-05-11 20:53 by feng6531
[硕博家园] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +8 Kelaizhang 2024-05-08 16/800 2024-05-11 19:58 by Kelaizhang
[电化学] 常用的国产电化学工作站有哪些? +4 123明湘 2024-05-11 4/200 2024-05-11 19:01 by sakuraai
[基金申请] NSFC函评几月几日截止提交评审意见? +7 瞬息宇宙 2024-05-11 7/350 2024-05-11 16:17 by 南海猪
[考博] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +5 Kelaizhang 2024-05-08 6/300 2024-05-11 14:34 by 明刃
[考博] 准研三生态毒理学,爱化学生物新能源碳中和,电池等,25博士岗位求捞啦! +3 wangchensi 2024-05-10 11/550 2024-05-11 09:41 by 斟索研学
[硕博家园] 盲审 +13 暮色恋伊人 2024-05-07 15/750 2024-05-10 21:28 by 暮色恋伊人
[电化学] 硫化锂正极TEM测试如何制样? +3 发仔子呀 2024-05-07 3/150 2024-05-10 13:14 by 姜尚真。
[有机交流] 酰胺合成 +6 Spiralup 2024-05-06 11/550 2024-05-10 10:10 by Spiralup
[基金申请] 基金函评会不会跨学部的送啊 +5 cls512 2024-05-09 5/250 2024-05-09 17:08 by cauzcl
[电化学] 水系电池正极制作 +3 shangfujun 2024-05-09 3/150 2024-05-09 12:03 by jacky601
[考博] 24-25申博 +5 muscling 2024-05-06 10/500 2024-05-08 21:33 by lizhengke06
[基金申请] 2024国家自然科学基金资助额度 +14 i6277 2024-05-05 15/750 2024-05-08 16:37 by andywei1028
[考研] 要不要二战? +9 ajjznskk 2024-05-05 13/650 2024-05-08 15:00 by ckxbear
[基金申请] 想问下,F口信息的函评是不是还没开始? +3 wowhow 2024-05-07 5/250 2024-05-08 12:31 by cuicyge
[教师之家] 请教下各位,南京航空航天大学天目湖校区如何? +3 pp爱写论文 2024-05-05 5/250 2024-05-07 16:48 by 似水流年Youth
信息提示
请填处理意见