| 查看: 836 | 回复: 11 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】热力学计算，temperature＝0K的情况
|
|||
|
诸位大侠，有没有计算过物质在0K时候的热力学？ 我在命令行中输入temperature＝0K的时候，发现输出的结果是298K的数值，并且输出行的开始也明确告诉你是temperature＝298K。 我更改了一下输入温度，比如设置temperature＝1K or 443K，最后输出的结果都分别是在1K和443K条件下的结果。 怎么temperature＝0K的时候和温度缺省条件下（temperature＝298K）的结果一样？ 今天遇到问题，有点糊涂了，请高手解答。 |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有153人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
不会，帮你顶一下 2楼2008-10-15 17:09:13
|
这个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
，你的电脑很好用啊，我用G3MP2这种方法的时候，我们的电脑算不动 4楼2008-10-17 09:09:40
5楼2008-10-17 10:25:42
|
还有就是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
7楼2008-10-17 16:05:23
8楼2008-10-17 16:24:20
9楼2008-10-17 16:33:01
10楼2008-10-17 16:34:07
