24小时热门版块排行榜    

查看: 1687  |  回复: 10

计算化学新人

铜虫 (小有名气)

[求助] 求大神们帮忙算一个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
回复此楼
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 计算化学新人 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 一志愿河工大 081700 276求调剂 +3 地球绕着太阳转 2026-03-23 3/150 2026-03-25 19:10 by 雾散后相遇lc
[考研] 招08考数学 +8 laoshidan 2026-03-20 17/850 2026-03-25 17:52 by 一个红太阳
[考研] 材料学硕,求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 资源与环境 调剂申请(333分) +7 holy J 2026-03-21 7/350 2026-03-24 17:24 by xiaohai104
[考研] 求调剂 +5 林之夕 2026-03-24 5/250 2026-03-24 17:16 by dick_runner
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
信息提示
请填处理意见