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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】溶剂化’有质子转移的过渡态寻找出错
51/1返回列表
查看: 1042  |  回复: 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-11 19:20:
在输入文件最后加上以下二行：
scfvac
radii=uahf

另：为什么一定要在溶剂条件下优化结构呢，这样很容易出错，一般在真空条件下优化后，只需计算加溶剂条件的单点能就可以了

[ Last edited by tutu2000  ...

再次请教：如果只计算真空条件下的优化，再在溶剂化下算单点能，那单点能零点能校正怎么办呢，用真空下算出来的那个单点能？
      那这样算出来的结果可靠吗，我想毕竟不全是在溶剂下得出的结果。在这方面很多困惑、、、谢谢啦

[ Last edited by ccc13 on 2008-9-12 at 09:38 ]
赞一下 回复此楼
9楼2008-09-12 09:36:47
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭