24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

shiy

木虫 (小有名气)

应助: 10 (幼儿园)
金币: 2825.4
帖子: 95
在线: 74.3小时
虫号: 693583

[交流] 【求助】利用gibbs包计算热容的问题

各位虫友，最近利用gibbs包算YCu的热容，测试例子通过，但是算YCu时没有结果，不知道问题出在哪儿？
YCu.inp如下：

YCu B2-CsCl phase
YCu.out
2
152.45185
-5040.755164
0
0 0.25
11 0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0
6 0.0 100.0 200.0 300.0 400.0 500.0
21
  311.7664          -5040.752158
  308.9322          -5040.752572
  306.0975          -5040.753117
  303.2640             -5040.753691
  300.4291          -5040.753866
  297.5946          -5040.754299
  294.7607          -5040.754544
  291.9264          -5040.754780
  289.0919          -5040.754972
  286.2576          -5040.755098
  283.4238          -5040.755124
  280.5894          -5040.755189
  277.7549          -5040.755125
  274.9217       -5040.754945
  272.0865          -5040.754740
  269.2532          -5040.754395
  266.4187          -5040.753873
  263.5843          -5040.753374
  260.7506          -5040.752824
  257.9153          -5040.752199
  255.0813          -5040.751396
数据是根据E-V曲线计算得到的，如图，图有拐点，结果没错，但是执行gibbs < YCu.inp，就会提示：
minbrack: monotonic function, there's no minimum
Static 10.0000000000000 1 9 15
9.046723932968691E-002 9.062748377090271E-002 9.096617006668994E-002
9.137931254146701E-002 9.184369293828577E-002 9.290311384214425E-002
9.366176615027956E-002 9.444473868047840E-002 9.534415306517782E-002
9.637253346193685E-002 9.845105633357000E-002 9.960647071853258E-002
0.100805817126108 0.102014027578541 0.103410441964146
1
大概意思是单调函数，没有最小值，不知道问题出在哪儿？请各位虫友帮忙分析一下。

[ Last edited by shiy on 2011-1-6 at 10:26 ]

回复此楼

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有9人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

第二次修改意见求助！已经有5人回复
德拜温度有什么用？已经有9人回复
Gaussian 热化学计算问题请教已经有8人回复
准谐德拜模型计算热力学性质已经有7人回复
关于等体热容Cv值得计算已经有5人回复
纯过氧化氢的热容是多少已经有5人回复
求助：计算出的热容单位换算已经有12人回复
求助材料专业同学：晶格热容、电子热容、电子热导率是可以测量的吗？已经有15人回复
锅炉热力性能计算软件已经有398人回复
四方晶系热膨胀系数的理论计算已经有5人回复
【讨论】沼液的开发利用问题已经有28人回复
【求助】linux下装adf源码包出现错误已经有3人回复
【求助】Gibbs反应器计算问题已经有6人回复
【求助】求2010年注册化工工程师基础考试的考试大纲已经有4人回复
【求助】函数最小值问题？已经有9人回复
【求助】热容流率概念已经有5人回复
【求助】热挤压锭包套已经有5人回复
【求助】如何计算ＴＤＢ中的Gibbs能已经有13人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-01-06 10:24:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 多谢提示 2011-07-16 12:17:36
shiy(金币+1): 谢谢交流 2011-07-16 15:29:56

这个程序不怎么靠谱。第一个，你在点多的一边去几个点，可能就能得到结果；第二，你在点少的一边加几个点，可能就能得到结果。第三，280后面那个点有点跳，把那个点去掉，可能也能得到结果。
以上三个处理方法得到的结果会有一定的差距……
所以，一般认为这个Gibbs程序使用的德拜法，就是个不严谨的凑数法，用来发个低端论文糊弄一下不懂行的人还是可以，真拿来做第一性，还是声子法（德拜法不属于第一性）比较靠谱……