24小时热门版块排行榜    

查看: 2354  |  回复: 4

2019211142

捐助贵宾 (初入文坛)

[求助] 求教:vasp+phonopy计算声子谱 已有1人参与

初学声子谱,DFPT方法。计算声子谱,由于我的原胞很大,超过10A,就建了1*1*1的超胞,然后将SPOSCAR改为POSCAR开始计算,但是在计算完成后,我输入phonopy --fc vasprun.xml,出现了以下提示,求教各位,不知道什么原因。后面是我的incar
Traceback (most recent call last):
  File "/home-gg/soft/packagedir/anaconda3/install/bin/phonopy", line 4, in <module>
    __import__('pkg_resources').run_script('phonopy==2.2.0', 'phonopy')
  File "/home-gg/soft/packagedir/anaconda3/install/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-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/pkg_resources/__init__.py", line 1438, in run_script
    exec(code, namespace, namespace)
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/phonopy-2.2.0-py3.6-linux-x86_64.egg/EGG-INFO/scripts/phonopy", line 429, in <module>
    error_num = create_FORCE_CONSTANTS(filename, write_hdf5, log_level)
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/phonopy-2.2.0-py3.6-linux-x86_64.egg/phonopy/interface/vasp.py", line 110, in create_FORCE_CONSTANTS
    fc_and_atom_types = parse_force_constants(filename)
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/phonopy-2.2.0-py3.6-linux-x86_64.egg/phonopy/interface/vasp.py", line 151, in parse_force_constants
    return vasprun.read_force_constants()
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/phonopy-2.2.0-py3.6-linux-x86_64.egg/phonopy/interface/vasp.py", line 517, in read_force_constants
    return self._get_force_constants(vasprun)
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/phonopy-2.2.0-py3.6-linux-x86_64.egg/phonopy/interface/vasp.py", line 533, in _get_force_constants
    for event, element in vasprun_etree:
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/site-packages/phonopy-2.2.0-py3.6-linux-x86_64.egg/phonopy/interface/vasp.py", line 600, in _parse_by_etree
    for event, elem in etree.iterparse(fileptr):
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/xml/etree/ElementTree.py", line 1227, in iterator
    root = pullparser._close_and_return_root()
  File "/home-gg/soft/packagedir/anaconda3/install/lib/python3.6/xml/etree/ElementTree.py", line 1274, in _close_and_return_root
    root = self._parser.close()
xml.etree.ElementTree.ParseError: no element found: line 1017, column 0

INCAR
System = Co3Sn2S2
PREC = Accurate
LREAL = .F.
IBRION = 8
IALGO=38
ADDGRID=.T.
ISIF = 2
POTIM=0.2
NSW = 1
EDIFFG=-0.05
ISMEAR=0
SIGMA=0.1
NELM = 100
ENCUT = 400
GGA=PE
VOSKOWN=1
EDIFF=1E-7
LCHARG=.F.
LWAVE=.F.
LDAU=.TRUE.
LDAUTYPE = 2
LDAUL = -1 2 -1
LDAUU = 0 15 0
LDAUJ = 0 0 0
LMAXMIX = 6
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc

新虫 (知名作家)

没算好vasprun.xml

发自小木虫IOS客户端

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

2楼2020-04-29 13:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2019211142

捐助贵宾 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by 123dingyc at 2020-04-29 13:05:42
没算好vasprun.xml

谢谢您的回复,请问vasprun.xml算错了,是incar出问题了吗

发自小木虫IOS客户端
3楼2020-04-29 13:34:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

j_o_y

新虫 (小有名气)

你好,请问你的问题解决了吗?我也遇到同样的问题,想请教你解决办法
4楼2020-07-21 10:35:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1649972541

新虫 (初入文坛)

【答案】应助回帖

你是对Co元素加的U值吗?U值加到15有点大吧?而且Co有3d轨道,没有f轨道。设置LDA+U时,对于s和p元素, LM AXM IX应设置为2;对于d元素, LM AXM IX应设置为4;对于f元素, LM AXM IX应设置为6,所以你的LMAXMIX应该是4吧?如有错误请指出
5楼2020-09-27 10:44:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 2019211142 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +9 羊腰板 2026-08-21 9/450 2026-08-22 00:03 by anji999
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
请填处理意见