24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4397  |  回复: 7

happy小海

新虫 (小有名气)

[求助] gibbs程序计算热力学性质的具体用法 已有1人参与

gibbs程序windows版在算热力学性质时,命令窗口应该输入一些什么指令?
按照英文用户说明书准备好的EV文件test01如下,输入readEVdata(test01)等没有反应,
test01.dat:
# volume ang^3
# energy ry
# nelectrons 2
# z 4
31.5176 -14.5664202
32.4533 -14.5804702
33.4073 -14.5938196
34.3799 -14.6065011
35.3711 -14.6185620
36.3812 -14.6300054
37.4103 -14.6408790
对程序运行、编译知识很匮乏,英文说明书看的听糊涂,我想问下我哪里出错了,该怎么做?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baowanrui110

铁杆木虫 (知名作家)

神虫浮云

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流 2013-12-02 22:17:33
happy小海: 金币+2, ★★★很有帮助 2013-12-06 15:47:45
首先,本软件基于Quasi-Harmonic Debye模型,与声子谱计算后的QHA模型起源于相同的公式。这是由于在统计热力学中,零点和振动自由能只能被定义为声子在FBZ内的积分。
其次,本软件使用其他第一性原理计算得到的E-V曲线和弹性泊松比来拟合EOS并且进一步拟合热力学性能。具体的是把声子与声速进行关联,分类成纵/横两种声速,然后纵/横声速又对应了弹性常数C11和C12(仅限各向同性材料,其只有两个相互独立的弹性常数C11和C12),从而完成了声子的积分。
第三,本软件计算效率非常高,无论多么大的体系,只要有E-V曲线,有泊松比,就可以在一秒内得到想要的热力学结果,不占空间不占时间。
第四,2004年左右发布的Gibbs程序使用Fortran 90写成,本此发行的使用bash脚本,不需要进行编译。具体的使用方法可参考Doc中的用户手册。
Gibbs第一版的发布文献:M.A.Blanco, E.Francisco, and V.Luana, Comput. Phys. Commu. 158, 57 (2004)
Gibbs第二版的发布文献:A.Otero-de-la-Roza and V. Luana, Comput. Phys. Commu. 182, 1708-1720 (2011)
           A.Otero-de-la-Roza and V. Luana, Comput. Phys. Commu. 182, 2232-2248 (2011)
文献和软件都附加在了附件中。
顺便提一下个人的观点,不见得正确,还希望有大牛来补充和说明:Debye模型提出的时间是非常早的,至少是在能够计算完整声子谱之前就已经存在的,不知道最早是怎么使用的,但是比较早的计算方面的应用,可以回溯到1988年,由纽约IBM工作的Moruzzi和Janak及维也纳科技大学的Schwarz发表在Phys Rev B, Volume 37, Number 2的论文Calculated thermal Properties of Metals。不知道是否是注意到了Debye模型中对弹性常数的各向同性近似,作者只在文章中计算了与各向同性最相近的立方结构的金属单质(立方体系中有三个相互独立的弹性常数C11,C12和C44)。后来在声子色散谱的计算大量普及之前,也有不少人结合LMTO、LCAO等计算方法,配合Debye模型计算热力学性能的。至于现在,基本主流都是使用声子色散谱+QHA的方法了,在Quantum-ESPRESSO中,甚至给出了三阶非谐波的计算(d3.x)。
另外,由于忽略了声子的计算,也就同时忽略了动力学稳定性的检验,比如虚频的问题。从原始公式可以知道,吉布斯函数中的零点能和振动自由能是需要对声子进行积分的,而虚频的存在则会让公式失去计算意义。但是Debye模型则可以“overcome”这个问题,在我个人看来,这就好像比较ln(-5)和ln(-2)的大小,使用ln(a)=0.5*2ln(a)=0.5*ln(a^2)然后得到ln(-5)>ln(-2)一样不靠谱~
至于E-V的取值,由于EOS对参数非常敏感,所以取值的密度和广度都对最终结果有着非常巨大的影响,如果是为了拟合已有的实验数据或者QHA计算数据,还是比较容易多快好省的做到的。作者使用自己的Gibbs2,在Phys. Rev. B, 84, 024109 (2011)上也发表了一篇文章。不过从我个人观察的来看,除非原创性的,否则比较少见使用该方法的高水平文章。
总之,这样的方法总的来说还是比较考验人的~当然,如果作为对实验的粗略预测,还是有不错的应用的。
最后一说的是,Debye模型一般不被计算领域的人承认为第一性原理方法的。
再补充个第一版的Gibbs,这样就从拟合到计算齐活了~
为照顾部分Windows用户,特此推出编译好的Windows版Gibbs v1.0程序,但是仅限于64-bit的Windows用户使用。
你必须非常努力,才能看起来还不费力。
2楼2013-12-01 21:59:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★ ★ ★ ★
franch: 金币+5, 捍卫原创。。。 加金币~ 2013-12-02 22:18:19
LS的抄得真爽……这不是我写的原文么……
思想重于技巧,内涵重于表象
3楼2013-12-01 22:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happy小海

新虫 (小有名气)

引用回帖:
3楼: Originally posted by souledge at 2013-12-01 22:57:54
LS的抄得真爽……这不是我写的原文么……

太赞了,你回我的帖子了,知道这个软件就是看的楼上贴出的你的那个帖子,附件我都下下来了,但源代码里有很多错误,自己改程序的能力约等于0,给老师看了一下老师说自己学一下试着慢慢改一改%>_<%,正在想办法呢,你是用源代码在matlab中运行出来的吗,还有一个包里有gibbs1.0exe程序,那是编译好的吗?可以直接用吗?用什么命令呢?
4楼2013-12-02 20:16:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happy小海

新虫 (小有名气)

引用回帖:
2楼: Originally posted by baowanrui110 at 2013-12-01 21:59:02
首先,本软件基于Quasi-Harmonic Debye模型,与声子谱计算后的QHA模型起源于相同的公式。这是由于在统计热力学中,零点和振动自由能只能被定义为声子在FBZ内的积分。
其次,本软件使用其他第一性原理计算得到的E-V ...

感谢回帖,这个帖子我看到过
5楼2013-12-02 20:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happy小海

新虫 (小有名气)

引用回帖:
4楼: Originally posted by happy小海 at 2013-12-02 20:16:37
太赞了,你回我的帖子了,知道这个软件就是看的楼上贴出的你的那个帖子,附件我都下下来了,但源代码里有很多错误,自己改程序的能力约等于0,给老师看了一下老师说自己学一下试着慢慢改一改%>_<%,正在想办 ...

6楼2013-12-11 18:00:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tonisam

银虫 (小有名气)

【答案】应助回帖

内容已删除
7楼2014-12-22 17:47:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hhq_nole

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by Tonisam at 2014-12-22 17:47:34
要用souledge提供的Gibbs,你得先学会cmd操作方法,这个在百度上搜一下就行很简单的。
你得通过cmd才能打开gibbs,并进行计算。
对于输入文件的编写,我觉得你好好钻研下软件包中的gibb.doc文件,这个就是输入文件 ...

你好,可以大致说说cmd运行gibbs的方法么,或者告知一下该搜些什么关键字,输入文件差不多明白了,但试了好久也不知道怎么打开gibbs,麻烦你了。

发自小木虫IOS客户端
8楼2017-11-30 19:26:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 happy小海 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿生物与医药,296分,求调剂 +9 66鹿 2026-04-03 10/500 2026-04-05 20:11 by lys0704
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +8 哇呼哼呼哼 2026-04-02 9/450 2026-04-05 17:14 by lbsjt
[考研] 求调剂 +4 wos666 2026-04-03 4/200 2026-04-05 11:48 by arrow8852
[考研] 270求调剂 +9 小杰pp 2026-03-31 11/550 2026-04-05 11:02 by 风雨无晴
[考研] 288求调剂,一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 0703求调剂 +6 zizimo 2026-03-31 6/300 2026-04-04 14:16 by 无际的草原
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦! 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 285求调剂 +7 AZMK 2026-04-02 9/450 2026-04-03 11:12 by wanwan00
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 13:03 by yulian1987
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[考研] 297 地理学070500 复试求调剂 +3 小圆圈圈ooo 2026-03-30 3/150 2026-03-30 21:05 by 余震yz
信息提示
请填处理意见