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的回帖

hlq3175

木虫 (正式写手)

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

引用回帖:

6楼: Originally posted by 雨晨田木 at 2012-03-15 09:29:41:
你看看这个吧Comment on the Correct Use of Continuum Solvent Models
dx.doi.org/10.1021/jp107136j

我要的自由能应该是Cavitation energy                      (kcal/mol) =    18.95
Dispersion energy                      (kcal/mol) =    -12.91
Repulsion energy                      (kcal/mol) =    1.17
Total non electrostatic                (kcal/mol) =    7.22这几项的和，不过我好像应该尝试一下在优化阶段就加上溶剂去试试，谢谢您了。

赞一下

回复此楼

7楼2012-03-16 10:49:08

已阅回复此楼关注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]