24小时热门版块排行榜    

查看: 6224  |  回复: 10

liuyunshenm

木虫 (正式写手)

[求助] phono3py计算晶格热导率提取数据phonon lifetime与声子频率的关系数据求助。 已有2人参与

大神,‘
    您好!
    我用的系统是Ubuntu 18.04; 安装了phono3py 1.16.0版本(最新版本),VASP+phono3py计算热导率的
     系统安装在WMware虚拟机中的。
      我在用phono3py-kdeplot提取数据画图时,总是报错,报错提示如下:
     (另外,Ubuntu 14.04+phono3py 1.12.3,  Ubuntu 16.04+phono3py 1.16.0 也基本是如下报错情况)
     请求大神能帮我从根本上解决问题。最好是能从文件中直接提取数据,我另外画图。


/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axes/_base.py:3215: MatplotlibDeprecationWarning:
The `xmin` argument was deprecated in Matplotlib 3.0 and will be removed in 3.2. Use `left` instead.
  alternative='`left`', obj_type='argument')
/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axes/_base.py:3221: MatplotlibDeprecationWarning:
The `xmax` argument was deprecated in Matplotlib 3.0 and will be removed in 3.2. Use `right` instead.
  alternative='`right`', obj_type='argument')
/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axes/_base.py:3604: MatplotlibDeprecationWarning:
The `ymin` argument was deprecated in Matplotlib 3.0 and will be removed in 3.2. Use `bottom` instead.
  alternative='`bottom`', obj_type='argument')
/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axes/_base.py:3610: MatplotlibDeprecationWarning:
The `ymax` argument was deprecated in Matplotlib 3.0 and will be removed in 3.2. Use `top` instead.
  alternative='`top`', obj_type='argument')
Traceback (most recent call last):
  File "/home/pkl/anaconda3/bin/phono3py-kdeplot", line 4, in <module>
    __import__('pkg_resources').run_script('phono3py==1.16.0.post1', 'phono3py-kdeplot')
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/pkg_resources/__init__.py", line 666, in run_script
    self.require(requires)[0].run_script(script_name, ns)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/pkg_resources/__init__.py", line 1446, in run_script
    exec(code, namespace, namespace)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/phono3py-1.16.0.post1-py3.6-linux-x86_64.egg-info/scripts/phono3py-kdeplot", line 294, in <module>
    main(get_options())
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/phono3py-1.16.0.post1-py3.6-linux-x86_64.egg-info/scripts/phono3py-kdeplot", line 290, in main
    fig.savefig("lifetime.png"
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/figure.py", line 2094, in savefig
    self.canvas.print_figure(fname, **kwargs)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/backend_bases.py", line 2075, in print_figure
    **kwargs)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/backends/backend_agg.py", line 510, in print_png
    FigureCanvasAgg.draw(self)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/backends/backend_agg.py", line 402, in draw
    self.figure.draw(self.renderer)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/artist.py", line 50, in draw_wrapper
    return draw(artist, renderer, *args, **kwargs)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/figure.py", line 1649, in draw
    renderer, self, artists, self.suppressComposite)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/image.py", line 138, in _draw_list_compositing_images
    a.draw(renderer)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/artist.py", line 50, in draw_wrapper
    return draw(artist, renderer, *args, **kwargs)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axes/_base.py", line 2628, in draw
    mimage._draw_list_compositing_images(renderer, self, artists)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/image.py", line 138, in _draw_list_compositing_images
    a.draw(renderer)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/artist.py", line 50, in draw_wrapper
    return draw(artist, renderer, *args, **kwargs)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axis.py", line 1187, in draw
    renderer)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/axis.py", line 1125, in _get_tick_bboxes
    extent = tick.label1.get_window_extent(renderer)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/text.py", line 930, in get_window_extent
    bbox, info, descent = self._get_layout(self._renderer)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/text.py", line 313, in _get_layout
    ismath=ismath)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/backends/backend_agg.py", line 209, in get_text_width_height_descent
    s, fontsize, renderer=self)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/texmanager.py", line 464, in get_text_width_height_descent
    dvifile = self.make_dvi(tex, fontsize)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/texmanager.py", line 328, in make_dvi
    texfile], tex)
  File "/home/pkl/anaconda3/lib/python3.6/site-packages/matplotlib/texmanager.py", line 298, in _run_checked_subprocess
    stderr=subprocess.STDOUT)
  File "/home/pkl/anaconda3/lib/python3.6/subprocess.py", line 356, in check_output
    **kwargs).stdout
  File "/home/pkl/anaconda3/lib/python3.6/subprocess.py", line 423, in run
    with Popen(*popenargs, **kwargs) as process:
  File "/home/pkl/anaconda3/lib/python3.6/subprocess.py", line 729, in __init__
    restore_signals, start_new_session)
  File "/home/pkl/anaconda3/lib/python3.6/subprocess.py", line 1364, in _execute_child
    raise child_exception_type(errno_num, err_msg, err_filename)
FileNotFoundError: [Errno 2] No such file or directory: 'latex': 'latex'
回复此楼
nopainnogian
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ftwoo

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
根据最后一句,不知是不是坐标的符号latex不能识别?如果有坐标的话,去掉看可否?

自己导出数据的话,你需要从一系列的h5py文件中提取数据的,需要自己写一个简单的程序,比如python。你可以从一个h5py文件中读取数据开始,可以从Togo的网站看他的例子。然后可以根据他的例子写一个简单的python程序读取所有的 kappa-mxxxxx.hdf5 文件中的所有mode的gamma
2楼2019-04-01 01:08:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ftwoo

银虫 (小有名气)

【答案】应助回帖

忘记贴Togo的例子了,地址在这里:

https://atztogo.github.io/phono3py/hdf5_howto.html
3楼2019-04-01 01:14:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyunshenm

木虫 (正式写手)

引用回帖:
2楼: Originally posted by ftwoo at 2019-04-01 01:08:15
根据最后一句,不知是不是坐标的符号latex不能识别?如果有坐标的话,去掉看可否?

自己导出数据的话,你需要从一系列的h5py文件中提取数据的,需要自己写一个简单的程序,比如python。你可以从一个h5py文件中读 ...

你好,我在Ubuntu 16.04+phono3py 1.16.0下,安装过latex,除了latex提示没有,其他错误提示依然存在,这个最明显的错误是在matplotlib上,没办法画图。另外,我不会pytho语言,更不会编写读取程序,所以望大神提供数据提取脚本。
nopainnogian
4楼2019-04-01 14:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyunshenm

木虫 (正式写手)

引用回帖:
3楼: Originally posted by ftwoo at 2019-04-01 01:14:17
忘记贴Togo的例子了,地址在这里:

https://atztogo.github.io/phono3py/hdf5_howto.html

python语言,我看不懂,还望指导
nopainnogian
5楼2019-04-01 14:26:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ftwoo

银虫 (小有名气)

引用回帖:
5楼: Originally posted by liuyunshenm at 2019-04-01 14:26:28
python语言,我看不懂,还望指导...

最好是使用python,Togo提供的例子就是python的,按照那个链接里的code基本可以解决。
如果你用其他编程语言,你可以照着那个链接里的程序一步步改成其他的语言。
6楼2019-04-02 03:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyunshenm

木虫 (正式写手)

引用回帖:
6楼: Originally posted by ftwoo at 2019-04-02 03:49:55
最好是使用python,Togo提供的例子就是python的,按照那个链接里的code基本可以解决。
如果你用其他编程语言,你可以照着那个链接里的程序一步步改成其他的语言。...

我试一试看看
nopainnogian
7楼2019-04-03 16:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc

新虫 (知名作家)

8楼2019-04-03 22:56:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyunshenm

木虫 (正式写手)

引用回帖:
6楼: Originally posted by ftwoo at 2019-04-02 03:49:55
最好是使用python,Togo提供的例子就是python的,按照那个链接里的code基本可以解决。
如果你用其他编程语言,你可以照着那个链接里的程序一步步改成其他的语言。...

已经试过多次了,还是提取不出声子寿命的数据。望你能提供其他解决办法
nopainnogian
9楼2019-04-04 19:17:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雷克萨斯哈哈

银虫 (小有名气)

楼主,你能告诉我怎么提取声子寿命吗
10楼2019-04-20 12:10:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuyunshenm 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题,想知道的进来 +7 且听虎啸 2026-08-07 10/500 2026-08-07 20:43 by 医学老男孩
[基金申请] 关于filecode +3 布布和一二 2026-08-07 6/300 2026-08-07 17:51 by 医学老男孩
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见