24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2185  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 11/550 2026-04-06 10:43 by 蓝云思雨
[考研] 085410人工智能 初试316分 求调剂 +7 残星拂曙 2026-03-31 7/350 2026-04-06 10:09 by 蓝云思雨
[考研] 298求调剂 +3 manman511 2026-04-05 3/150 2026-04-05 18:09 by kk112233
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 344材料与化工调剂 +9 调剂上岸玘 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 调剂0855-288 +5 x熊二a 2026-04-03 5/250 2026-04-04 00:19 by 猪会飞
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 0710生物学,325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[硕博家园] 考研调剂 +5 骆驼男人 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 0710生物学求调剂! +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
信息提示
请填处理意见