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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】溶剂化’有质子转移的过渡态寻找出错
52/1返回列表
查看: 1040  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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就是转移的质子。
请问能不能给指导一下，该怎样修改才能往下算。心切啊。谢谢啦！！！
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2008-09-02 09:09:32
已阅   回复此楼   关注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的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭