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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求大神们帮忙算一个IRC
111/2返回列表12下一页
查看: 1987  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

计算化学新人

铜虫 (小有名气)

[求助] 求大神们帮忙算一个IRC 已有2人参与

小弟遇到一个IRC问题，前面的ts优化得到的振动很不错，频率-800左右，振动模式也是我想要的
但是在接下来做IRC后发现，IRC图像基本就是个对称的，产物和反应物一模一样
后来又分别用了forward和reverse去算，还是那样
所以求大神们帮个忙看一看到底是怎么回事
跪谢！

输入文件如下

%nprocshared=16
%chk=1.chk
# irc=(maxpoints=100,calcfc) ub3lyp/genecp geom=connectivity
pressure=1 temperature=673.15

Title Card Required

-1 2
C                  4.43756300   -0.96670700    0.94520300
H                  4.96127500   -1.31901600    0.04765900
H                  5.19942800   -0.59545500    1.64802500
H                  3.94498500   -1.82992000    1.40218600
C                  3.42655900    0.13453200    0.61628700
H                  2.88696000   -0.37553600   -0.61164700
H                  2.89376600    0.41782800    1.53024600
C                  4.07374500    1.38926100    0.02265800
H                  4.60025400    1.17255900   -0.91563800
H                  3.31524900    2.14724100   -0.18813500
H                  4.80779200    1.81367900    0.72460500
Pt                -1.17647800    0.17646700    1.43840500
Pt                 1.37489300   -0.50952800   -0.10824100
Pt                -1.12039500   -1.16475100   -0.74681000
Pt                -0.41439000    1.43660800   -0.75184400

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 7 1.0 8 1.0
6 13 1.0
7
8 9 1.0 10 1.0 11 1.0
9
10
11
12 14 1.0 15 1.0
13 14 1.0 15 1.0
14 15 1.0
15

H 0
6-31g**
****
C 0
6-31g**
****
Pt
lanl2dz
****

Pt 0
lanl2dz
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-12-23 08:20:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xm-xi

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
我之前遇到过这个起。概况，你把b3lyp改成mp2,试一下.或者设置step=3或8试试。忘记具体操作了。

发自小木虫Android客户端
赞一下 回复此楼
2楼2015-12-23 08:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

计算化学新人

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by xm-xi at 2015-12-23 08:31:29
我之前遇到过这个起。概况，你把b3lyp改成mp2,试一下.或者设置step=3或8试试。忘记具体操作了。

改方法吗？但是我这是一连串的体系，要改方法的话，那么所有的都要改了啊
step是指maxstep吗？
赞一下 回复此楼
3楼2015-12-23 09:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xm-xi

银虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 计算化学新人 at 2015-12-23 09:08:12
改方法吗？但是我这是一连串的体系，要改方法的话，那么所有的都要改了啊
step是指maxstep吗？...

应该是有印象。我毕业后改行了。好多年没用了。

发自小木虫Android客户端
赞一下 回复此楼
4楼2015-12-24 19:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

计算化学新人

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by xm-xi at 2015-12-24 19:23:48
应该是有印象。我毕业后改行了。好多年没用了。
...

这个东西我调过，没有用啊
赞一下 回复此楼
5楼2015-12-25 07:59:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不取俗名

木虫 (正式写手)
总觉得gjf哪里有点不太对劲

发自小木虫Android客户端
赞一下 回复此楼
6楼2015-12-25 08:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

计算化学新人

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by 不取俗名 at 2015-12-25 08:05:19
总觉得gjf哪里有点不太对劲

比如说？我这个关键词基本是我做IRC通用的，其他的都挺好，就这个体系做出来有问题
赞一下 回复此楼
7楼2015-12-25 15:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lankadesufei

木虫 (小有名气)

大嘻哈家

【答案】应助回帖

感谢参与，应助指数 +1
可能是过渡态定位精度不够，使用opt=(calcall,tight)，int=ultrafine关键词更精确的确定过渡态结构。也可以先尝试IRC计算中将步长设的大一些，比如stepsize=15。
赞一下 回复此楼
厚积薄发，厚德载物！
8楼2015-12-25 16:41:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

计算化学新人

铜虫 (小有名气)
引用回帖:
8楼: Originally posted by lankadesufei at 2015-12-25 16:41:35
可能是过渡态定位精度不够，使用opt=(calcall,tight)，int=ultrafine关键词更精确的确定过渡态结构。也可以先尝试IRC计算中将步长设的大一些，比如stepsize=15。

请问int=ultrafine这个关键词是什么意思，做什么的？谢谢~
赞一下 回复此楼
9楼2015-12-27 18:36:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lankadesufei

木虫 (小有名气)

大嘻哈家
引用回帖:
9楼: Originally posted by 计算化学新人 at 2015-12-27 18:36:01
请问int=ultrafine这个关键词是什么意思，做什么的？谢谢~...

DFT 计算在Hartree-Fock 计算的每一主要阶段上添加一个另外的步骤。这一步是泛函（或各种泛函的导数）的数值积分。因此除了来源于Hartree-Fock 计算的数值误差（积分的精度，SCF 收敛，CPHF 收敛）之外，DFT 计算的精度还依赖于数值积分中使用的网格点数。
Gaussian 09 默认为“精密的”积分网格（对应于Integral=FineGrid）。该网格以最少的额外耗时最大限度地提高计算精度。不推荐在正式的DFT 计算中使用更疏的网格。还要注意，在比较能量时（如计算能量差，生成热，等），所有计算需要使用相同的积分网格。需要的话可以使用更密的网格（如对某些分子体系进行严格的几何优化计算），int=ultrafine就是增加格点精度。
PS：这些在gaussian说明书里就有
赞一下 回复此楼
厚积薄发，厚德载物！
10楼2015-12-27 21:27:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 计算化学新人 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
关闭
请填处理意见
关闭