24小时热门版块排行榜    

查看: 2160  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-15 7/350 2026-03-16 09:43 by 闲人终南山
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
信息提示
请填处理意见