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 ]

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有21人回复
导师想让我从独立一作变成了共一第一已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
申请2026年博士已经有6人回复

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

第二次修改意见求助！已经有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的回帖

hyj0752

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 132.5
帖子: 10
在线: 15.7小时
虫号: 815890

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

小弟现在也遇到了同样的问题，我计算了NiAl
NiAl.inp的内容如下：
NiAl-B2
NiAl.out
2
85.69
-42.02795262
0
0 0.25
17 0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0 55.0 60.0 65.0 70.0 75.0 80.0
17 0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0 450.0 500.0 550.0 600.0 650.0 700.0 750.0 800.0
17
127.13877389  -41.99380467
131.32274628  -42.00269369
135.57420210  -42.01000686
139.96062476  -42.01588269
144.41453085  -42.02045468
149.00340380  -42.02384201
153.65976016  -42.02615536
158.38359996  -42.02749451
163.24240660  -42.02795159
168.16869668  -42.02761134
173.22995360  -42.02655258
178.35869395  -42.02484748
183.62240115  -42.02256243
188.95359177  -42.01975833
194.41974925  -42.01649132
200.02087357  -42.01281037
205.62199790  -42.00876436
运行gibbs 出现minbrack: monotonic function, there's no minimum
Static 65.000000000000000             1          9       16
  0.31503262247819380    0.31538715327285544    0.31746662360608613    0.32128161297226093    0.32654953187700148    0.33330028985887611    0.34127411733514057    0.35037123437205953    0.36064860044192321    0.37187238607245132    0.38411286073592693    0.39714876493797269    0.41106279819505526    0.42564497096861337    0.44098823281930893    0.45704361370295465
STOP 1
不知道怎么办呀，求高手解答呀。