请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 184|回复: 7

[Gaussian/gview] 求助:走irc两边一样,想问还有什么建议。

[复制链接]

37

帖子

0

威望

532

eV
积分
569

Level 4 (黑子)

发表于 2019-7-10 14:49:20 | 显示全部楼层 |阅读模式
本帖最后由 哈哈哈 于 2019-7-11 10:08 编辑

优化出了过渡态,虚频的振动是我想要的。把过渡态拿去走irc,走出来两边一样。
把过渡态的输出结构弄成输入文件加calcall tight int=ultrafine 重新优化,结构偏离。
稍微修改了结构 仍加 calcall tight int=ultrafine优化,虚频的振动是我想要的。把这个结构走irc 步长为5 irc中用calcall int=ultrafine,两边仍是一样的。步长为20 reverse (看到之前irc 左边一步原子1 4距离缩短)calcall int=ultrafine走irc,感觉趋势仍然1 4在远离,像产物结构。想问还有什么其它方法能走正常的irc来找反应物和产物。 麻烦老师了

1万

帖子

25

威望

2万

eV
积分
43584

管理员

公社社长

发表于 2019-7-10 20:34:14 | 显示全部楼层
上传文件时一定要把压缩包里每个文件都对应什么说清楚。你的压缩包里还有ORCA的任务,令人糊涂

类似问题,该说的这里都说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400http://bbs.keinsci.com/thread-7736-1-1.html
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

51

帖子

0

威望

1538

eV
积分
1589

Level 5 (御坂)

发表于 2019-7-10 20:54:55 | 显示全部楼层
本帖最后由 mutron 于 2019-7-10 20:56 编辑

你的TS:
E(UB3LYP) =  -1796.27392904     A.U. after    1 cycles
<S**2>= 0.9945

你的IRC第一点:
E(UB3LYP) =  -1796.27324240     A.U. after   45 cycles
<S**2>= 0.7818

明显IRC第一点(对应TS)没有收敛到TS的波函数,IRC任务加上guess=read,读取你TS所得到的chk文件。顺便也需要考察一下TS波函数是否稳定



评分

参与人数 1eV +1 收起 理由
哈哈哈 + 1 谢谢

查看全部评分

37

帖子

0

威望

532

eV
积分
569

Level 4 (黑子)

 楼主| 发表于 2019-7-11 14:23:10 | 显示全部楼层
sobereva 发表于 2019-7-10 20:34
上传文件时一定要把压缩包里每个文件都对应什么说清楚。你的压缩包里还有ORCA的任务,令人糊涂

类似问题 ...

不好意思  之后会注意的

37

帖子

0

威望

532

eV
积分
569

Level 4 (黑子)

 楼主| 发表于 2019-7-11 14:37:11 | 显示全部楼层
mutron 发表于 2019-7-10 20:54
你的TS:
E(UB3LYP) =  -1796.27392904     A.U. after    1 cycles
= 0.9945

guess=read 走irc就对了,谢谢你哈。
我通常走irc都是提取ts的结构做输入来走,没有用guess=read。想问一下,走irc是不是最好以guess=read的过渡态来走?

1万

帖子

25

威望

2万

eV
积分
43584

管理员

公社社长

发表于 7 天前 | 显示全部楼层
哈哈哈 发表于 2019-7-11 14:37
guess=read 走irc就对了,谢谢你哈。
我通常走irc都是提取ts的结构做输入来走,没有用guess=read。想问 ...

一般没必要,但极个别情况下IRC任务的默认初猜波函数下SCF可能收敛到与过渡态优化最后不同的波函数,相当于走IRC时候和找过渡态时不在同一个势能面上,此时会出现IRC不正常
如果你想保险,每次跑IRC都加guess=read也可以,没坏处
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

37

帖子

0

威望

532

eV
积分
569

Level 4 (黑子)

 楼主| 发表于 7 天前 | 显示全部楼层
sobereva 发表于 2019-7-12 05:08
一般没必要,但极个别情况下IRC任务的默认初猜波函数下SCF可能收敛到与过渡态优化最后不同的波函数,相当 ...

感觉这个可以添到irc常见问题那,多一种尝试方式,嘿嘿。

1万

帖子

25

威望

2万

eV
积分
43584

管理员

公社社长

发表于 6 天前 | 显示全部楼层
哈哈哈 发表于 2019-7-12 09:33
感觉这个可以添到irc常见问题那,多一种尝试方式,嘿嘿。

目前就在文中
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-7-19 02:28 , Processed in 0.149550 second(s), 24 queries .

快速回复 返回顶部 返回列表