24小时热门版块排行榜    

Znn3bq.jpeg
查看: 848  |  回复: 11
当前主题已经存档。

majun04

木虫 (著名写手)

majun

[交流] 【求助】热力学计算,temperature=0K的情况

诸位大侠,有没有计算过物质在0K时候的热力学?
我在命令行中输入temperature=0K的时候,发现输出的结果是298K的数值,并且输出行的开始也明确告诉你是temperature=298K。
我更改了一下输入温度,比如设置temperature=1K or 443K,最后输出的结果都分别是在1K和443K条件下的结果。
怎么temperature=0K的时候和温度缺省条件下(temperature=298K)的结果一样?
今天遇到问题,有点糊涂了,请高手解答。
回复此楼

» 猜你喜欢

乐观、执着、勤奋
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓儿165

铜虫 (正式写手)

不会,帮你顶一下
做个聪明优雅的女人
2楼2008-10-15 17:09:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by 皓儿165 at 2008-10-15 17:09:
不会,帮你顶一下

这个gaussian是不用单独计算0K的情况的,因为他的计算都是先计算0K,然后在这个基础上,再计算你指定的温度。所以你需要0K的数值可以直接在输出文件里面读取出来的。
例如:
Temperature=              298.150000 Pressure=                      1.000000
E(ZPE)=                     0.076807 E(Thermal)=                    0.081167
E(QCISD(T))=             -154.561808 E(Empiric)=                   -0.092790
DE(MP2)=                   -0.213161
G3MP2(0 K)=              -154.790952 G3MP2 Energy=               -154.786592
G3MP2 Enthalpy=          -154.785647 G3MP2 Free Energy=          -154.816364


求助了以后,仔细想了想,感觉是这样的,不知道理解的对不对。
乐观、执着、勤奋
3楼2008-10-15 17:30:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓儿165

铜虫 (正式写手)

引用回帖:
Originally posted by majun04 at 2008-10-15 17:30:


这个gaussian是不用单独计算0K的情况的,因为他的计算都是先计算0K,然后在这个基础上,再计算你指定的温度。所以你需要0K的数值可以直接在输出文件里面读取出来的。
例如:
Temperature=              298 ...

我觉得你理解的对了,我刚查了些资料,也是这么说的,你的电脑很好用啊,我用G3MP2这种方法的时候,我们的电脑算不动
做个聪明优雅的女人
4楼2008-10-17 09:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by 皓儿165 at 2008-10-17 09:09:

我觉得你理解的对了,我刚查了些资料,也是这么说的,你的电脑很好用啊,我用G3MP2这种方法的时候,我们的电脑算不动

我的电脑是单机,配置不是特别高,4月份买的,大约7000元。在计算过程中也有“算不动”的情况。感觉用G3MP2方法计算,算不动的情况是不一定就是电脑配置不行。我前一阵子用这个方法计算含有锌原子的络合物,结果也是非正常结束。后来才明白,是我选择的计算方法G3(G3MP2)里面设置的基组,不能用来计算锌原子,所以程序非正常结束。我觉得要看你电脑具体的配置和你调用的方法和基组来分析具体是什么原因计算失败。
乐观、执着、勤奋
5楼2008-10-17 10:25:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

还有就是G3MP2调用的方法和基组如下:

%  testrt
TestRt Revision:  
Please type in the route specification, terminated with a blank line:
g3MP2

-----
g3MP2
-----
1/18=20,38=1/1,3;
2/9=110,17=6,18=5,40=1/2;
3/5=1,6=6,7=1,11=9,16=1,25=1,30=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7//1,2,3,16;
1/18=20/3(3);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99//99;
2/9=110/2;
3/5=1,6=6,7=1,11=9,16=1,25=1,30=1/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/18=20/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;

%chk=temp
#P g3MP2

--Link1--
%chk=temp
#P Geom=AllCheck Guess=Read SCRF=Check HF/6-31G(d) Freq

--Link1--
%chk=temp
#P Geom=AllCheck Guess=Read SCRF=Check MP2(Full)/6-31G(d) Opt=RCFC

--Link1--
%chk=temp
#P Geom=AllCheck Guess=Read SCRF=Check QCISD(T)/6-31G(d)

--Link1--
%chk=temp
#P Geom=AllCheck Guess=Read SCRF=Check MP2/GTMP2Large

你看看是否适合你计算的分子。应用范围看手册P34。

[ Last edited by majun04 on 2008-10-17 at 10:37 ]
乐观、执着、勤奋
6楼2008-10-17 10:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓儿165

铜虫 (正式写手)

引用回帖:
Originally posted by majun04 at 2008-10-17 10:34:
还有就是G3MP2调用的方法和基组如下:

%  testrt
TestRt Revision:  
Please type in the route specification, terminated with a blank line:
g3MP2

-----
g3MP2
-----
1/18=20,38=1/1,3;
  ...

太感谢了哈
我没有找到类是g3mp2(0K)=......这几行,是不是因为我用的密度泛含,没有这几行输出,不过我觉得可以推出来:用 sum of electronic and zero-point energies减去zero-point  correction不就是总的电子能量吗,这样对不
做个聪明优雅的女人
7楼2008-10-17 16:05:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by 皓儿165 at 2008-10-17 16:05:

太感谢了哈
我没有找到类是g3mp2(0K)=......这几行,是不是因为我用的密度泛含,没有这几行输出,不过我觉得可以推出来:用 sum of electronic and zero-point energies减去zero-point  correction不就是总 ...

因为方法不一样,所以最后输出结果的信息也不一样。
输出结果各个能量代表什么,建议你看看这个书:
Exploring Chemistry with Electronic Structure Methods
这个书比较几本全面的介绍gaussian的应用,适合初学者。
乐观、执着、勤奋
8楼2008-10-17 16:24:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓儿165

铜虫 (正式写手)

引用回帖:
Originally posted by majun04 at 2008-10-17 16:24:


因为方法不一样,所以最后输出结果的信息也不一样。
输出结果各个能量代表什么,建议你看看这个书:
Exploring Chemistry with Electronic Structure Methods
这个书比较几本全面的介绍gaussian的应用,适 ...

好的,应该有电子版的吧,我找一下,但愿能找到
做个聪明优雅的女人
9楼2008-10-17 16:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

引用回帖:
Originally posted by 皓儿165 at 2008-10-17 16:33:

好的,应该有电子版的吧,我找一下,但愿能找到

信箱告诉我,发给你一个电子版的。
乐观、执着、勤奋
10楼2008-10-17 16:34:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 majun04 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 297,工科调剂?河南农业大学本科 +15 河南农业大学-能 2026-04-14 15/750 2026-04-20 12:39 by Equinoxhua
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
信息提示
请填处理意见