应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Wien2k&FLAPW&ELK » 【求助】利用gibbs包计算热容的问题
北京石油化工学院2026年研究生招生接收调剂公告
141/1返回列表
查看: 3206  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shiy

木虫 (小有名气)

[交流] 【求助】利用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 ]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-01-06 10:24:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

souledge

专家顾问 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 多谢提示 2011-07-16 12:17:36
shiy(金币+1): 谢谢交流 2011-07-16 15:29:56
这个程序不怎么靠谱。第一个，你在点多的一边去几个点，可能就能得到结果；第二，你在点少的一边加几个点，可能就能得到结果。第三，280后面那个点有点跳，把那个点去掉，可能也能得到结果。
以上三个处理方法得到的结果会有一定的差距……
所以，一般认为这个Gibbs程序使用的德拜法，就是个不严谨的凑数法，用来发个低端论文糊弄一下不懂行的人还是可以，真拿来做第一性，还是声子法（德拜法不属于第一性）比较靠谱……
赞一下(4人) 回复此楼
6楼2011-07-15 22:55:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ufo991

铁虫 (正式写手)
★
zhang668(金币+1):多谢交流 2011-01-09 14:30:24
shiy(金币+1): 2011-01-09 14:53:05
你的输出文件什么都没有吗？换几个点试试，可能是数值处理时某一点没通过，然后给终止了吧。优化应该没错、
赞一下(1人) 回复此楼
2楼2011-01-09 13:34:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shiy

木虫 (小有名气)
引用回帖:
Originally posted by ufo991 at 2011-01-09 13:34:41:
你的输出文件什么都没有吗？换几个点试试，可能是数值处理时某一点没通过，然后给终止了吧。优化应该没错、

输出文件啥也没有，就一个title
赞一下 回复此楼
3楼2011-01-09 14:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kuner

铁虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 说说具体原因？谢谢！ 2011-07-16 12:18:15
程序给的不对，明显
赞一下(2人) 回复此楼
4楼2011-05-09 23:52:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomagh

铁杆木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
youzhizhe(金币+1): 谢谢提示。 2011-05-12 18:28:02
shiy(金币+1): 谢谢回复，有道理 2011-05-14 20:13:36
uuv2010(金币+1): 多谢提示 2011-07-16 12:17:57
这个问题我想是这样的，你在压缩的这一边再算2-3个点，加进去再用GIbbs，估计就可以了。
赞一下(3人) 回复此楼
5楼2011-05-12 17:19:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiangxue

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
uuv2010(金币+1): 多谢提示 2011-07-16 12:17:47
shiy(金币+1): 谢谢交流 2011-07-16 15:29:18
应该去几个点，点太密了
赞一下(2人) 回复此楼
7楼2011-07-16 00:04:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyj0752

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主的问题解决了吗，我刚学gibbs这个程序，也遇到了这样的问题，按照上面的方法做了几遍，也没解决，不知道怎么办才好，郁闷呀，求高手解答呀

[ Last edited by hyj0752 on 2012-9-15 at 22:52 ]
赞一下 回复此楼
8楼2012-09-15 22:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyj0752

新虫 (初入文坛)
★
小木虫: 金币+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
不知道怎么办呀，求高手解答呀。
赞一下 回复此楼
9楼2012-09-15 23:05:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shiy

木虫 (小有名气)
引用回帖:
1619950楼: Originally posted by hyj0752 at 2012-09-15 23:05:25
小弟现在也遇到了同样的问题，我计算了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  ...

减小体积应变量
回复此楼
10楼2012-09-17 16:54:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴空微岚

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
YCu.inp中，第五行，那个能量你怎么不转换成Hartree 能量 ？ 是不是这里出了问题
赞一下 回复此楼
11楼2013-04-28 10:59:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

晴空微岚

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你在Gibbs里，整个输入文件里的能量都需要转化为以Hartree 为单位的能量，然后，第五行填写的能量是转化后的Hartree单位的能量里面最小的能量（要记得分正负号）
赞一下 回复此楼
12楼2013-04-28 11:05:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ip0371

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主：

minbrack: monotonic function, there's no minimum

的问题找到原因了吗？

我也遇到相同的问题了。

如果你知道的话，还请你回复一下
赞一下 回复此楼
13楼2014-08-11 21:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pokm9000

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by souledge at 2011-07-15 22:55:02
这个程序不怎么靠谱。第一个，你在点多的一边去几个点，可能就能得到结果；第二，你在点少的一边加几个点，可能就能得到结果。第三，280后面那个点有点跳，把那个点去掉，可能也能得到结果。
以上三个处理方法得到 ...

大神，您所说的声子法是什么方法，能向您请教一下吗
赞一下 回复此楼
14楼2017-04-26 17:45:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shiy 的主题更新
141/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 080200学硕，机械工程专业277分，求带走！ +7 瓶子PZ 2026-03-31 7/350 2026-04-05 17:49 by liucky
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 学硕288调剂!!! +3 小王xw123 2026-04-03 3/150 2026-04-03 21:20 by 啵啵啵0119
[考研] 一志愿生物与医药，296分，求调剂 +8 66鹿 2026-04-03 9/450 2026-04-03 14:22 by 化学化工硕士招�
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 土木水利328分求调剂 +6 疾风知劲草666 2026-04-02 6/300 2026-04-03 11:38 by znian
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 一志愿山东大学，085600，344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 09:15 by olim
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[考研] 085601 329分调剂 +6 yzsa12 2026-03-31 6/300 2026-03-31 15:23 by yanflower7133
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
信息提示
关闭
请填处理意见
关闭