24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1177  |  回复: 11
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ccc13

[交流] 【求助】溶剂化’有质子转移的过渡态寻找出错

溶剂化效应计算时 多次出现类似的情况,请前辈们帮忙看看是什么毛病。谢谢!

输入命令(溶剂化效应下找过渡态):
%chk=402ts.chk
%rwf=x1,200MW,x2,200MW,x3,200MW,x4,200MW,x5,200MW,x6,200MW,x7,-1
%mem=500mb
%nproc=1
# Opt(QST3,Noeigen,maxcycle=1000) B3LYP/6-31+g(d,p) pop=nbo freq scrf=(IEFPCM,read,solvent=water)
……(下面是产物、反映物、过渡态的分子信息)(用过这组命令曾经计算出过过渡态)

输出文件最后部分:
United Atom Topological Model (UA0  parameters set).
UA0: Hydrogen   13 is unbound. Keep it explicit at all point on the
UA0: potential energy surface to get meaningful results.
Error termination via Lnk1e in e:\g03\l301.exe at Wed Jan 01 12:43:01 2003.
Job cpu time:  0 days  0 hours  0 minutes 10.0 seconds.
File lengths (MBytes):  RWF=   1338 Int=      0 D2E=      0 Chk=      1 Scr=      1


ps:真空状态下能算下去。

我比较了一下用同一行命令,那个曾经算出来的和这些算不出来的区别是:算不出来的共同点是在溶剂化效应下找过渡态都有质子氢的转移,出来结果的没有质子氢转移。而且out文件中也说“UA0: Hydrogen   13 is unbound. Keep it explicit at all point on the
UA0: potential energy surface to get meaningful results.”其中H13就是转移的质子。
请问能不能给指导一下,该怎样修改才能往下算。心切啊。谢谢啦!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ccc13

引用回帖:
Originally posted by tutu2000 at 2008-9-13 17:32:
结果可靠不可靠,要依据实验数据,理论计算本就用很多基组与方法的选择问题,这方面要多看文献才能了解,一般来说可以这样进行能量的比较:
△Gsol(溶液条件下的相关能)=△E(真空条件下总能量)+△E0(零点能 ...

好不容易找到个明白人,真诚的再次请教:
     我在溶剂化条件光opt就算不下去,因此想找个别的办法。看到好些文献也是像您说的那样,一般在真空条件下优化后,只需计算加溶剂条件的单点能这样做得。
   
     但不明白的就是如果在真空条件下优化后,只计算加溶剂条件的单点能,那这样△G‘sol(溶剂化能)如何得到呢?【因为要用到△Gsol(溶液条件下的相关能)=△E(真空条件下总能量)+△E0(零点能)+△G‘sol(溶剂化能)】

     那溶剂化下的单点能的零点能校正,那个零点能又该在哪个文件里寻找呢?
11楼2008-09-16 15:13:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

ccc13

似乎解决了一小点问题,贴出来请大家在此基础上给指导一下

  对某氢原子不同其他原子成键的结构,例如质子转移反应的过渡态,UA0模型将无法构建空穴。此时可以使用Gaussian 03的SPHEREONH=N关键字在该氢原子中心处添加球体,N为输入文件中该原子在原子列表中的序数。这个关键字应该是文件末尾一个空行之后另起一行输入。
   感谢大话西游中的“fado"大侠


    我这样做了,优化能算一圈了,出了四个“NO"以后,再进行下一轮计算又出现了楼上一样的错误。也就是只起了一次作用。哎,愁死了

[ Last edited by ccc13 on 2008-9-11 at 19:09 ]
2楼2008-09-05 08:58:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alick1102

铜虫 (小有名气)

引用回帖:
Originally posted by ccc13 at 2008-9-2 09:09:
溶剂化效应计算时 多次出现类似的情况,请前辈们帮忙看看是什么毛病。谢谢!

输入命令(溶剂化效应下找过渡态):
%chk=402ts.chk
%rwf=x1,200MW,x2,200MW,x3,200MW,x4,200MW,x5,200MW,x6,200MW,x7,-1
%mem ...

为什么用闭壳层方法优化过渡态? 一般都要用开窍层方法的。
到了博士的时候才知道,科学也不过是为了吃饭的。
3楼2008-09-05 09:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我不是马甲!

木虫 (小有名气)

请问一下楼上的,计算TS,都要用开壳层吗?
对这个问题一直很困惑,搜了很多这方面的帖子,有的用闭壳层,有的用开壳层。真不知按哪个做好。
4楼2008-09-05 23:50:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +8 花田半亩~白 2026-04-21 8/400 2026-04-22 14:59 by jiaoxg
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[考研] 0854求调剂 +24 门路摸摸 2026-04-15 28/1400 2026-04-22 09:54 by Sy199704!
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
信息提示
请填处理意见