24小时热门版块排行榜    

查看: 1644  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 本人考085602 化学工程 专硕 +11 不知道叫什么! 2026-03-15 13/650 2026-03-16 15:10 by 了了了了。。
[考研] 环境工程调剂 +4 大可digkids 2026-03-16 4/200 2026-03-16 14:00 by laoshidan
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 一志愿安徽大学材料工程专硕313分,求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
信息提示
请填处理意见