| 查看: 3705 | 回复: 13 | |||
[交流]
【求助】利用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 ] |
回复此楼» 猜你喜欢
娱乐
已经有3人回复
-
面上提前没消息，有中的吗
已经有17人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 第二次修改意见求助！
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
国自然发榜在即：此刻的你，是什么心情？
+2/948
-
助理研究员——吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/475
-
2027年博士生（硕博连读/直博）吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/472
-
山东征女友，坐标济南
+1/175
-
威格科技张德龙：二十年深耕，铸就国产惰性气氛设备新高度
+1/88
-
华东理工大学 生物工程学院 生物反应器工程国家重点室 博士后招聘
+1/76
-
★★★中山大学（深圳）院士团队招收27届推免硕士生（低维材料与器件方向）
+1/74
-
2026年University of South Carolina继续招收博士生
+1/62
-
中科院东莞材料所依托中科大招收2027届推免生
+1/33
-
有机合成高级研究员
+1/30
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘（长期有效）
+5/25
-
招收博士后（神经退行性疾病、认知障碍疾病）
+1/22
-
组委会诚挚滴邀请各界学者踊跃投稿参会
+1/18
-
国科金评审方式改革建议
+1/11
-
瑞典林雪平大学 高峰课题组 诚聘 氮化物基无机材料方向 博士后1名
+1/11
-
北京师范大学（珠海校区）赵亚然副教授招收2027年电催化推免硕士/博士研究生
+1/5
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/5
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
宁波诺丁汉大学新能源与环境材料技术研究中心招收2027年春季入学博士生
+1/1
-
致每一位等待国自然结果的科研人：我们想听听你的心声
+1/1
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 多谢提示 2011-07-16 12:17:36
shiy(金币+1): 谢谢交流 2011-07-16 15:29:56
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 多谢提示 2011-07-16 12:17:36
shiy(金币+1): 谢谢交流 2011-07-16 15:29:56
|
这个程序不怎么靠谱。第一个，你在点多的一边去几个点，可能就能得到结果；第二，你在点少的一边加几个点，可能就能得到结果。第三，280后面那个点有点跳，把那个点去掉，可能也能得到结果。 以上三个处理方法得到的结果会有一定的差距…… 所以，一般认为这个Gibbs程序使用的德拜法，就是个不严谨的凑数法，用来发个低端论文糊弄一下不懂行的人还是可以，真拿来做第一性，还是声子法（德拜法不属于第一性）比较靠谱…… |
6楼2011-07-15 22:55:02
2楼2011-01-09 13:34:41
3楼2011-01-09 14:54:06
4楼2011-05-09 23:52:37
5楼2011-05-12 17:19:39
7楼2011-07-16 00:04:48
8楼2012-09-15 22:51:20
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+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 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 不知道怎么办呀，求高手解答呀。 |
9楼2012-09-15 23:05:25
10楼2012-09-17 16:54:12
11楼2013-04-28 10:59:38
12楼2013-04-28 11:05:32
13楼2014-08-11 21:59:52
14楼2017-04-26 17:45:51
