|
本帖最后由 飞翔的猪 于 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: 48)
|