24小时热门版块排行榜    

查看: 769  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见