24小时热门版块排行榜

返回列表

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

pokm9000

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 59
帖子: 16
在线: 20.4小时
虫号: 3323013
注册: 2014-07-14
专业: 凝聚态物性I:结构、力学和

[求助] 想把MS得出的数据利用Gibbs算热力学性质，单位应该如何转化？

如下表所示是小弟用GGA-PBE优化了一个立方晶系的物质，设置了一系列晶格常数所得到的数据，左列为晶体的晶格常数，单位是埃；中间列是晶体体积，单位是埃的三次方；右列是得到的能量，单位是电子伏特Ev。
6.3 192.486 -12999.5469
6.35 197.106 -13000.17144
6.4 201.799 -13000.62863
6.43 204.65 -13000.82637
6.45 206.565 -13000.92935
6.48 209.461 -13001.04055
6.5 211.406 -13001.08584
6.515 212.873 -13001.10947
6.53 214.347 -13001.11608
6.55 216.323 -13001.11066
6.58 219.309 -13001.06313
6.6 221.315 -13001.01029
6.63 224.346 -13000.896
6.65 226.383 -13000.79594
6.68 229.46 -13000.61508
6.7 231.527 -13000.47305
而Gibbs软件所用的体积单位是波尔的三次方V(bohr^3)，能量单位是hartree。这样的话上述数据就需要转化。1bohr=0.529177A，1Ev=0.0366993672hartree，这样的话我的体积数据就需要乘以（1/0.529177）^3也就是6.74833，能量需要乘以0.036699，由此可以由Gibbs软件的形式表示如下：
ZnBO rutile phase. aiPI-uCHF calculation.
ZnBO.out
23
534.3782
-477.132856
2
0 0.280347
7 0.0 5.0 10.0 15.0 20.0 25.0 30.0
6 0.0 100.0 200.0 300.0 400.0 500
16
1298.959626 -477.0751453
1330.136924 -477.0980654
1361.806851 -477.1148439
1381.046348 -477.122101
1393.969406 -477.12588
1413.512579 -477.1299612
1426.638086 -477.1316234
1436.537891 -477.1324905
1446.484934 -477.1327329
1459.81964 -477.1325341
1479.970162 -477.1307897
1493.507318 -477.1288507
1513.961515 -477.1246564
1527.70787 -477.1209841
1548.47249 -477.1143468
1562.421294 -477.1091343
输出的结果是说我是单调函数，无最小值。当然得到的输出文件里面也是什么都没有。

大幅度.png

但是我抱着试一试的态度，不转化直接用MS得到的数据输入Gibbs，这时却是运行正常的，虽然得到的数据是错的。
所以求教大神，是不是之前我数据转化错了（我有这种想法的原因之一是感觉转化后的数值跟我看到过别人贴出来的相比太大了，别人都是几百，我是上千）。如果不是，求大神指点我的问题出在了哪里？

回复此楼