| 查看: 2940 | 回复: 10 | ||||
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说明书里没找到他们的解释呀，请各位高手给我指点指点。感激不尽呀! |
回复此楼» 收录本帖的淘帖专辑推荐
 计算和软件 |
» 猜你喜欢
26年储能、电池方向博士申请
已经有1人回复
-
求调剂，本科吉大，一志愿华工339
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有293人回复
-
意大利博洛尼亚大学招聘研究助理2名（2D材料/锂离子与钠离子电池方向）
已经有1人回复
-
关于碳量子点抗菌的研究
已经有15人回复
-
南大化学调剂
已经有11人回复
-
0703调剂
已经有2人回复
-
1
已经有20人回复
-
0703 调剂
已经有5人回复
-
西南某211招27届博士
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高斯溶剂化能求助
已经有3人回复
- 自定义溶剂的溶剂化能量计算如何输入
已经有3人回复
- 溶剂化能计算问题
已经有7人回复
- 溶剂化效应出现gaussian 9999错误
已经有7人回复
- 请教一下 关于表面应力 和 比表面自由能的问题 ！
已经有6人回复
- 用SMD算溶剂下单点能，自由能应该用哪个
已经有8人回复
- 请教高斯09算溶剂效应及能量校正问题，多送金币+小红花
已经有19人回复
- 求带有溶剂化效应计算的文章，不胜感激！
已经有6人回复
- 量化中的溶剂效应，什么情况加
已经有3人回复
- 为什么只有溴化氢有过氧化物效应
已经有6人回复
- 关于溶剂化效应的计算
已经有3人回复
- 怎么读取Siesta输出的Hamiltonian and overlap矩阵(*.HSX文件)？
已经有16人回复
- 溶剂化效应输入
已经有6人回复
- 急需帮助！！！！聚丙烯酸甲酯溶于苯乙烯吗？
已经有10人回复
- 能够溶于油的染料或者颜料，但不会给脲醛树脂染色
已经有4人回复
- 求助《钛合金金相学》中发现的问题
已经有6人回复
- 【求助】溶剂化模型中关于溶剂的问题？
已经有5人回复
- 【求助】DMOL 构建簇模型如何添加点电荷来消除边界效应
已经有15人回复
- 【求助】为什么在构象搜索等过程中总是寻求自由能最低？
已经有3人回复
- 【求助】Comsol多相流模块如何自定义自由能函数项？
已经有4人回复
- 【求助】吉普斯自由能在哪能查到啊？
已经有11人回复
- 【求助】吉普斯自由能在哪能查到啊？
已经有3人回复
- 【转帖】极化连续介质模型（PCM）
已经有12人回复
雨晨田木
金虫 (小有名气)
- 应助: 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
感谢参与，应助指数 +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 |
2楼2012-03-14 10:33:58
雨晨田木
金虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 1177.7
- 散金: 91
- 红花: 1
- 帖子: 219
- 在线: 49.6小时
- 虫号: 1012487
- 注册: 2010-05-06
- 性别: GG
- 专业: 理论和计算化学
3楼2012-03-14 16:44:41
雨晨田木
金虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 1177.7
- 散金: 91
- 红花: 1
- 帖子: 219
- 在线: 49.6小时
- 虫号: 1012487
- 注册: 2010-05-06
- 性别: GG
- 专业: 理论和计算化学
4楼2012-03-14 16:47:01
hlq3175
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 2378.8
- 散金: 107
- 红花: 4
- 帖子: 444
- 在线: 301.2小时
- 虫号: 870062
- 注册: 2009-10-12
- 性别: MM
- 专业: 高分子合成化学
|
我是这样理解的，但是我觉得用频率分析做出来的那个自由能不能包括所有的自由能，我按您说的已经试过了，输出文件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 ======================= 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
雨晨田木
金虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 1177.7
- 散金: 91
- 红花: 1
- 帖子: 219
- 在线: 49.6小时
- 虫号: 1012487
- 注册: 2010-05-06
- 性别: GG
- 专业: 理论和计算化学
6楼2012-03-15 09:29:41
hlq3175
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 2378.8
- 散金: 107
- 红花: 4
- 帖子: 444
- 在线: 301.2小时
- 虫号: 870062
- 注册: 2009-10-12
- 性别: MM
- 专业: 高分子合成化学
7楼2012-03-16 10:49:08
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
8楼2012-03-16 11:14:28
boylc789
木虫 (著名写手)
独孤
- 应助: 3 (幼儿园)
- 金币: 3144.7
- 散金: 82
- 红花: 5
- 帖子: 1775
- 在线: 539.1小时
- 虫号: 579648
- 注册: 2008-07-11
- 专业: 理论和计算化学
9楼2012-03-31 13:39:39
雨晨田木
金虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 1177.7
- 散金: 91
- 红花: 1
- 帖子: 219
- 在线: 49.6小时
- 虫号: 1012487
- 注册: 2010-05-06
- 性别: GG
- 专业: 理论和计算化学
10楼2012-03-31 15:42:11
50