24小时热门版块排行榜

返回列表

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

hlq3175

木虫 (正式写手)

应助: 7 (幼儿园)
金币: 2378.8
散金: 107
红花: 4
帖子: 444
在线: 301.2小时
虫号: 870062
注册: 2009-10-12
性别: MM
专业: 高分子合成化学

[求助] 【求助】溶剂化效应自由能的输出问题

各位高手帮帮忙呀，我需要用g03算一下我的体系的溶剂化自由能 DeltaG (solv)，可是我按照常规的输入#p ccsd(t)/6-311++g(d,p) geom=modredundant scf=tight，它并没有那一项结果输出。根据小木虫里面找的的帖子，我的输入文件改成这样的#p ccsd(t)/6-311++g(d,p) scrf=(pcm,dovacuum,externaliteration) geom=modredundant scf=tight，可是总是出错类似QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.
t)/6-311++g(d,p) scrf=(pcm,dovacuum,exte
                        '
Last state="SCR2"
TCursr=38161 LCursr= 35
Error termination via Lnk1e in /home/apps/g03/l1.exe at Mon Mar 12 15:39:46 2012.
Job cpu time:  0 days  0 hours  0 minutes  1.2 seconds.
File lengths (MBytes):  RWF=    8 Int=    0 D2E=    0 Chk=    1 Scr=    1
这些输入文件里面的字母都是我一个一个敲进去的，没有任何空格，可是一提交就出现莫名其妙的空格，不知道是什么情况？另外能不能解释一下这两个dovacuum,externaliteration的意思，尤其是第一个，我在03说明书里没找到他们的解释呀，请各位高手给我指点指点。感激不尽呀!

回复此楼

1楼2012-03-12 22:47:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

雨晨田木

金虫 (小有名气)

应助: 17 (小学生)
金币: 1177.7
散金: 91
红花: 1
帖子: 219
在线: 49.6小时
虫号: 1012487
注册: 2010-05-06
性别: GG
专业: 理论和计算化学

引用回帖:

9楼: Originally posted by boylc789 at 2012-03-31 13:39:39:
能否解释下，这篇文章的内容，有点没看明白
作者的意思是说用方程4来计算溶剂化自由能吗

对，但是有个符号写错了，Esolv应该改成Esoln。就是说溶液下溶质的总自由能= Esoln（溶液下溶质的总能量）+Gnes（溶解增加的自由能）+deltaGcorr_gas（气相下溶质的热自由能校正）。Esoln+Gnes的值在PCM计算后在log文件中作为Total free energy in solution: with all non electrostatic terms 项自动给出，而deltaGcorr_gas是在气相下作freq计算后得到的值。注意，在这里溶液和气相下溶质的结构需要分别在液相和气相重作结构优化，很有可能是不一样的。

赞一下

回复此楼

[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&amp;amp;amp;say=不要太辛苦哦。[/img][/url]

10楼2012-03-31 15:42:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答

雨晨田木

金虫 (小有名气)

应助: 17 (小学生)
金币: 1177.7
散金: 91
红花: 1
帖子: 219
在线: 49.6小时
虫号: 1012487
注册: 2010-05-06
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gmy1990(金币+2): 2012-03-14 11:04:18
hlq3175(金币+1): ★有帮助, 谢谢您，不过我要的是溶剂化自由能，这样子算出来还是没有我要的自由能那一项。 2012-03-14 14:43:25

溶剂效应这样写就行了
%chk=H2O-S.chk
# opt rb3lyp/6-31+g(d,p) SCRF(CPCM,Read,Solvent=Acetonitrile)

H2O Structure in Solution-Phase

0 1
H                0.00000000 0.77125000 -0.46932000
H                0.00000000 -0.77125000 -0.46932000
O                0.00000000 0.00000000 0.11733000

RADII=UAKS

赞一下

回复此楼

[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&amp;amp;amp;say=不要太辛苦哦。[/img][/url]

2楼2012-03-14 10:33:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖