应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】溶剂化效应自由能的输出问题
51/1返回列表
查看: 2951  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hlq3175

木虫 (正式写手)

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

各位高手帮帮忙呀，我需要用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

木虫 (正式写手)
引用回帖:
楼: Originally posted by 雨晨田木 at 2012-03-14 16:47:01:
溶剂化自由能不就是溶解后的自由能减去溶解前的吗

我是这样理解的，但是我觉得用频率分析做出来的那个自由能不能包括所有的自由能，我按您说的已经试过了，输出文件SCF Done:  E(RB+HF-LYP) =  -399.126337746     A.U. after    7 cycles
             Convg  =    0.8999D-08             -V/T =  2.0039
             S**2   =   0.0000
KE= 3.975926270936D+02 PE=-1.499897705926D+03 EE= 4.215325170137D+02
--------------------------------------------------------------------
Variational PCM results
=======================
                     (a.u.) =    -399.104877
                     (a.u.) =    -399.126338
Total free energy in solution:
  with all non electrostatic terms            (a.u.) =    -399.114839
--------------------------------------------------------------------
(Polarized solute)-Solvent               (kcal/mol) =     -13.47
--------------------------------------------------------------------
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
--------------------------------------------------------------------
Partition over spheres:
Sphere  on Atom  Surface  Charge   GEl     GCav    GDR
    1      N1      14.71  -0.106   -3.05    2.51   -1.52
    2      C4      27.03  -0.023   -0.32    3.28   -2.30
    3      C7      24.01  -0.039   -0.52    3.47   -2.07
    4      O10     19.37   0.053   -6.10    2.53   -1.82
    5      C12     12.00  -0.109   -0.90    1.91   -1.10
    6      O13     18.00   0.100   -0.52    2.61   -1.56
    7      O14     15.46   0.084   -0.48    2.65   -1.38
Added spheres:    24.76   0.019   -1.57    0.00    0.00
--------------------------------------------------------------------
After PCM corrections, the SCF energy is  -399.126337746     a.u.
--------------------------------------------------------------------
Leave Link  502 at Wed Mar 14 11:20:35 2012, MaxMem=  838860800 cpu:      80.8
(Enter /home/apps/g03/l701.exe)这里面哪个是我要的自由能呢？难道是Total free energy in solution:
  with all non electrostatic terms            (a.u.) =    -399.114839减去RHF吗？我觉得这样做出来结果是不对的。
赞一下 回复此楼
5楼2012-03-14 20:24:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

雨晨田木

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +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=雨晨田木&say=不要太辛苦哦。[/img][/url]
2楼2012-03-14 10:33:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨晨田木

金虫 (小有名气)

【答案】应助回帖

对于???一个分??结构，分别计算有无溶剂效应时的???布斯自由能，其差值自然??就是溶剂化自由能????
赞一下 回复此楼
[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&say=不要太辛苦哦。[/img][/url]
3楼2012-03-14 16:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨晨田木

金虫 (小有名气)

【答案】应助回帖

★ ★
御剑江湖: 金币+2, 谢谢回帖交流 2012-03-17 18:58:41
溶剂化自由能不就是溶解后的自由能减去溶解前的吗
赞一下 回复此楼
[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&say=不要太辛苦哦。[/img][/url]
4楼2012-03-14 16:47:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 329求调剂 +5 想上学吖吖 2026-03-19 5/250 2026-03-19 23:56 by 23Postgrad
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 296求调剂 +3 www_q 2026-03-18 6/300 2026-03-19 22:28 by zhq0425
[考研] 085600材料与化工 求调剂 +14 enenenhui 2026-03-13 15/750 2026-03-19 21:53 by allen-yin
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
关闭
请填处理意见
关闭