24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1694  |  回复: 11

schalke

至尊木虫 (著名写手)

[求助] 各位师兄师姐请进,为什么找完过渡态后IRC跑出来的峰顶部有问题?

模拟过渡态后,只有一个虚频,振动方向也对,但是IRC跑出的峰顶是平的,很奇怪,是什么原因?怎么改才是正确的?

输入文件如下:
%mem=8GB
%nprocshared=6
%chk=20120612_Dimer_acid_ts_freq_simplized.chk
#p irc=(maxpoints=10,rcfc) b3lyp/6-31+g(d,p) maxdisk=40GB

20120612_Dimer_acid_ts_irc_simplized

0 1
C                  0.40233500    3.21378000    1.96150600
H                  1.27882900    3.86279900    2.07424400
H                 -0.32329700    3.76141800    1.32885100
C                  2.08893900   -0.31068100   -0.77122100
H                  2.16849500   -0.60881300   -1.82569300
C                  3.23509200   -0.97147600    0.01411800
H                  3.17320900   -0.66520900    1.06593000
H                  3.09715400   -2.06200100    0.00465800
C                  4.62235200   -0.62845200   -0.54533100
H                  4.76005700    0.46287800   -0.53014900
H                  4.67482500   -0.92696400   -1.60287300
C                  5.77270100   -1.29211400    0.22351000
H                  5.71759400   -0.99927900    1.28216000
H                  5.63982200   -2.38379200    0.20308900
C                  7.16299200   -0.94348000   -0.32422400
H                  7.29699800    0.14706400   -0.30113500
H                  7.21798800   -1.23434100   -1.38263800
C                  8.30514200   -1.61338100    0.44852300
H                  8.29797200   -1.31457700    1.50359000
H                  9.28234900   -1.34373400    0.03316900
H                  8.21913400   -2.70607000    0.41444000
C                  0.73921800    1.91133000    1.32206100
H                  2.22621400    0.77479000   -0.75792400
C                  0.25690300    0.69658700    1.77251800
C                  0.70238900   -0.71931100   -0.26267300
C                  0.27353600   -0.52056000    1.07644100
H                 -0.06597500    3.07807400    2.94245300
C                 -0.54391900    1.73160400   -1.10104200
C                 -0.49407100    0.35047800   -1.44024500
C                 -1.78674100   -0.46014200   -1.63761200
H                 -1.54374700   -1.32679100   -2.26625500
H                 -2.46943600    0.15655000   -2.23969700
C                 -2.56147300   -0.97643100   -0.41167900
H                 -2.65381600   -0.19110600    0.34823600
H                 -1.99929300   -1.78831300    0.06378900
C                 -3.95763400   -1.49547500   -0.78372400
H                 -3.86376700   -2.25497000   -1.57423400
H                 -4.54650800   -0.67563900   -1.22031200
C                 -4.72514200   -2.09487100    0.40204300
H                 -4.81499000   -1.33914400    1.19636200
H                 -4.13957400   -2.92012100    0.83325200
C                 -6.12371800   -2.60882400    0.03527700
H                 -6.03428700   -3.36110300   -0.76087100
H                 -6.71169700   -1.78380700   -0.39045800
C                 -6.87941300   -3.21312200    1.22443600
H                 -7.87188700   -3.57056800    0.92886200
H                 -7.01619200   -2.47484700    2.02368300
H                 -6.33280400   -4.06310300    1.65015800
C                 -1.70533700    2.44913300   -0.65386200
H                 -2.56681600    1.86380600   -0.33504700
C                 -1.82616900    3.80196000   -0.62268100
H                 -0.98652000    4.40960000   -0.96381000
C                 -3.07050300    4.54248200   -0.22484700
H                 -3.47401200    5.12732800   -1.06329900
H                 -2.87823800    5.25724000    0.58635200
H                 -3.85424200    3.85553800    0.11160600
H                  1.47452100    1.94070000    0.53045900
H                 -0.32104500    0.71418300    2.69794200
H                 -0.31814600   -1.31824700    1.51877200
H                  0.46032600   -1.72375600   -0.62145300
H                  0.19027900    0.18268200   -2.27161900
H                  0.28564700    2.34742700   -1.44234600
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

人工置顶,求解啊
2楼2012-06-15 09:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

第一次 见 开眼了...我想知道 top的2个结构 是一样的么...把ts再优化一下?
3楼2012-06-15 09:25:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
3楼: Originally posted by meteoric30 at 2012-06-15 09:25:23
第一次 见 开眼了...我想知道 top的2个结构 是一样的么...把ts再优化一下?

过渡态已经优化完毕,而且过渡态结构对应的是top2个点的中间位置,但并没有在IRC中显示出过渡态的点,但计算结果中含有过渡态结构。。。我是先算出过渡态的频率,然后读取其chk文件做irc,是否和读取chk文件有关?
4楼2012-06-15 09:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

引用回帖:
4楼: Originally posted by schalke at 2012-06-15 09:39:25
过渡态已经优化完毕,而且过渡态结构对应的是top2个点的中间位置,但并没有在IRC中显示出过渡态的点,但计算结果中含有过渡态结构。。。我是先算出过渡态的频率,然后读取其chk文件做irc,是否和读取chk文件有关?...

2点的能量和过渡态的能量也不一样吗?我意思 还是再优化下ts 再freq
5楼2012-06-15 17:14:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
5楼: Originally posted by meteoric30 at 2012-06-15 17:14:48
2点的能量和过渡态的能量也不一样吗?我意思 还是再优化下ts 再freq...

不一样,过渡态的结构能量最高,应该在两点中间。我是做完ts和freq后,用结果做的irc。我现在改用rcfc再计算一遍看看。另外,我算过渡态的时候用了nosymm、noeigentest、maxcyc的指令,但算irc的时候都没有用,会不会影响结果?
6楼2012-06-15 19:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

引用回帖:
6楼: Originally posted by schalke at 2012-06-15 19:05:19
不一样,过渡态的结构能量最高,应该在两点中间。我是做完ts和freq后,用结果做的irc。我现在改用rcfc再计算一遍看看。另外,我算过渡态的时候用了nosymm、noeigentest、maxcyc的指令,但算irc的时候都没有用,会不 ...

应该不会~~你有没有重新 算一下freq
7楼2012-06-16 09:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
7楼: Originally posted by meteoric30 at 2012-06-16 09:18:29
应该不会~~你有没有重新 算一下freq...

我先算了初始构型的频率,然后利用chk文件优化过渡态,再计算过渡态的频率,然后利用过渡态chk文件计算irc
8楼2012-06-16 09:20:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

引用回帖:
8楼: Originally posted by schalke at 2012-06-16 09:20:52
我先算了初始构型的频率,然后利用chk文件优化过渡态,再计算过渡态的频率,然后利用过渡态chk文件计算irc...

没问题~~我意思 再优化下过渡态~算freq 再看看irc应该没什么问题了
9楼2012-06-16 09:32:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

schalke

至尊木虫 (著名写手)

引用回帖:
9楼: Originally posted by meteoric30 at 2012-06-16 09:32:57
没问题~~我意思 再优化下过渡态~算freq 再看看irc应该没什么问题了...

我在linux下面算了一次,windows下面算了一次,都是这结果,太诡异。。。
10楼2012-06-16 09:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 schalke 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿A区211,22408 321求调剂 +7 随心所欲☆ 2026-04-15 8/400 2026-04-21 08:22 by Equinoxhua
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
信息提示
请填处理意见