计算化学公社

标题: 表面吸附物种熵的计算：脚本报错该如何解决 [打印本页]

作者
Author: 飞翔的猪 时间: 2019-8-18 22:02
标题: 表面吸附物种熵的计算：脚本报错该如何解决
本帖最后由飞翔的猪于 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。

作者
Author: ggdh 时间: 2019-8-19 00:02
装个anaconda3就行了

作者
Author: 卡开发发 时间: 2019-8-19 00:43

ImportError: No module named scipy

需要安装scipy。

作者
Author: 飞翔的猪 时间: 2019-8-19 08:57
本帖最后由飞翔的猪于 2019-8-19 10:22 编辑

ggdh 发表于 2019-8-19 00:02
装个anaconda3就行了

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

作者
Author: 飞翔的猪 时间: 2019-8-19 08:57
本帖最后由飞翔的猪于 2019-8-19 10:21 编辑

卡开发发发表于 2019-8-19 00:43
需要安装scipy。

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

欢迎光临计算化学公社 (http://ccc.keinsci.com/)