计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9171|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 表面吸附物种熵的计算:脚本报错该如何解决

[复制链接 Copy URL]

112

帖子

0

威望

851

eV
积分
963

Level 4 (黑子)

本帖最后由 飞翔的猪 于 2019-8-19 10:27 编辑

最近在重复大师兄网站的脚本解决表面吸附物种熵的计算https://www.bigbrosci.com/2018/11/07/ex69/时候,出现了如下报错:输入的脚本一样,不知道他这个提示是需要修改什么内容所以想请教大家
报错:
Traceback (most recent call last):
  File "/bin/entropy.py", line 7, in <module>
    from scipy import constants as con
ImportError: No module named scipy

自己电脑:#!/usr/bin/env python
# -*- coding: utf-8 -*-
#Statistical thermodynamics: the vibrational contribution,Page 642

import sys
import math
from scipy import constants as con

script, nu = sys.argv

nu = float(nu) * 100 #convert cm-1 to m-1

h_p = con.h # 6.62606975E-34 # J*s Plank constant
k_b = con.k # 1.38064852E-23 #m2*kg*s-2*k-1 Boltzman constant
R_gas = con.R # 8.3144598 # J*mol-1*K-1 Gas constant
l_s = con.c # 299792458 # light speed m * s-1
Tem = 298.15 # Temperature K
beta = 1/(k_b * Tem)

def get_pf(nu): #get partition function
    x_i = h_p * float(nu) * l_s * beta
    pf_l = x_i / (math.exp(x_i) - 1) # left part in the entropy equation
    pf_r = math.log(1 - math.exp(-x_i))
    pf  = pf_l - pf_r
    entropy = R_gas * pf
    return entropy

entropy = get_pf(nu) # J * K-1 * mol-1
ts      = entropy * Tem /1000 / 96.485  # in ev
print('%.4f\t %.4f' %(entropy , ts))

网站脚本:














该问题非常感谢ggdh和卡开发发(名字不分前后)的帮助下成功解决,把解决的方法放在这里分享,非常感谢大家的帮助

该报错的出现主要是缺少scipy库,所以解决方法:先安装python,再安装numpy,最后安装scipy,具体参考帖子链接:https://blog.csdn.net/shomy_liu/article/details/48543449




ex69-6.png (130.12 KB, 下载次数 Times of downloads: 49)

ex69-6.png

908

帖子

37

威望

5435

eV
积分
7083

Level 6 (一方通行)

2#
发表于 Post on 2019-8-19 00:02:20 来自手机 | 只看该作者 Only view this author
装个anaconda3就行了

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

3#
发表于 Post on 2019-8-19 00:43:43 | 只看该作者 Only view this author
ImportError: No module named scipy

需要安装scipy。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

112

帖子

0

威望

851

eV
积分
963

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2019-8-19 08:57:13 | 只看该作者 Only view this author
本帖最后由 飞翔的猪 于 2019-8-19 10:22 编辑
ggdh 发表于 2019-8-19 00:02
装个anaconda3就行了

谢谢,这个开源的包挺好的 里面包含很多的科学包,可能需要的内存有些大还是非常感谢您的帮助

112

帖子

0

威望

851

eV
积分
963

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-8-19 08:57:34 | 只看该作者 Only view this author
本帖最后由 飞翔的猪 于 2019-8-19 10:21 编辑

感谢您的提示,已经找到对应的资源成功安装,谢谢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 18:24 , Processed in 0.205290 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list